mirror of
https://github.com/opencv/opencv.git
synced 2024-12-13 16:09:23 +08:00
97681bdfce
imgproc: add optimized warpAffine kernels for 8U/16U/32F + C1/C3/C4 inputs #25984 Merge wtih https://github.com/opencv/opencv_extra/pull/1198. Merge with https://github.com/opencv/opencv_contrib/pull/3787. ### Pull Request Readiness Checklist See details at https://github.com/opencv/opencv/wiki/How_to_contribute#making-a-good-pull-request - [x] I agree to contribute to the project under Apache 2 License. - [x] To the best of my knowledge, the proposed patch is not based on a code under GPL or another license that is incompatible with OpenCV - [x] The PR is proposed to the proper branch - [x] There is a reference to the original bug report and related work - [x] There is accuracy test, performance test and test data in opencv_extra repository, if applicable Patch to opencv_extra has the same branch name. - [x] The feature is well documented and sample code can be built with the project CMake
260 lines
8.8 KiB
C++
260 lines
8.8 KiB
C++
// This file is part of OpenCV project.
|
|
// It is subject to the license terms in the LICENSE file found in the top-level directory
|
|
// of this distribution and at http://opencv.org/license.html.
|
|
#include "perf_precomp.hpp"
|
|
|
|
namespace opencv_test {
|
|
|
|
enum{HALF_SIZE=0, UPSIDE_DOWN, REFLECTION_X, REFLECTION_BOTH};
|
|
|
|
CV_ENUM(BorderMode, BORDER_CONSTANT, BORDER_REPLICATE)
|
|
CV_ENUM(InterType, INTER_NEAREST, INTER_LINEAR)
|
|
CV_ENUM(InterTypeExtended, INTER_NEAREST, INTER_LINEAR, WARP_RELATIVE_MAP)
|
|
CV_ENUM(RemapMode, HALF_SIZE, UPSIDE_DOWN, REFLECTION_X, REFLECTION_BOTH)
|
|
|
|
typedef TestBaseWithParam< tuple<Size, InterType, BorderMode, MatType> > TestWarpAffine;
|
|
typedef TestBaseWithParam< tuple<Size, InterType, BorderMode, int> > TestWarpPerspective;
|
|
typedef TestBaseWithParam< tuple<Size, InterType, BorderMode, MatType> > TestWarpPerspectiveNear_t;
|
|
typedef TestBaseWithParam< tuple<MatType, Size, InterTypeExtended, BorderMode, RemapMode> > TestRemap;
|
|
|
|
void update_map(const Mat& src, Mat& map_x, Mat& map_y, const int remapMode, bool relative = false );
|
|
|
|
PERF_TEST_P( TestWarpAffine, WarpAffine,
|
|
Combine(
|
|
Values( szVGA, sz720p, sz1080p ),
|
|
InterType::all(),
|
|
BorderMode::all(),
|
|
Values(CV_8UC3, CV_16UC3, CV_32FC3, CV_8UC1, CV_16UC1, CV_32FC1, CV_8UC4, CV_16UC4, CV_32FC4)
|
|
)
|
|
)
|
|
{
|
|
Size sz, szSrc(512, 512);
|
|
int type, borderMode, interType;
|
|
sz = get<0>(GetParam());
|
|
interType = get<1>(GetParam());
|
|
borderMode = get<2>(GetParam());
|
|
type = get<3>(GetParam());
|
|
Scalar borderColor = Scalar::all(150);
|
|
|
|
Mat src(szSrc,type), dst(sz, type);
|
|
switch (src.depth()) {
|
|
case CV_8U: {
|
|
cvtest::fillGradient<uint8_t>(src);
|
|
if(borderMode == BORDER_CONSTANT) cvtest::smoothBorder<uint8_t>(src, borderColor, 1);
|
|
break;
|
|
}
|
|
case CV_16U: {
|
|
cvtest::fillGradient<uint16_t>(src);
|
|
if(borderMode == BORDER_CONSTANT) cvtest::smoothBorder<uint16_t>(src, borderColor, 1);
|
|
break;
|
|
}
|
|
case CV_32F: {
|
|
cvtest::fillGradient<float>(src);
|
|
if(borderMode == BORDER_CONSTANT) cvtest::smoothBorder<float>(src, borderColor, 1);
|
|
break;
|
|
}
|
|
}
|
|
Mat warpMat = getRotationMatrix2D(Point2f(src.cols/2.f, src.rows/2.f), 30., 2.2);
|
|
declare.in(src).out(dst);
|
|
|
|
TEST_CYCLE() warpAffine( src, dst, warpMat, sz, interType, borderMode, borderColor );
|
|
|
|
SANITY_CHECK(dst, 1);
|
|
}
|
|
|
|
PERF_TEST_P( TestWarpPerspective, WarpPerspective,
|
|
Combine(
|
|
Values( szVGA, sz720p, sz1080p ),
|
|
InterType::all(),
|
|
BorderMode::all(),
|
|
Values(1, 3, 4)
|
|
)
|
|
)
|
|
{
|
|
Size sz, szSrc(512, 512);
|
|
int borderMode, interType, channels;
|
|
sz = get<0>(GetParam());
|
|
interType = get<1>(GetParam());
|
|
borderMode = get<2>(GetParam());
|
|
channels = get<3>(GetParam());
|
|
|
|
Scalar borderColor = Scalar::all(150);
|
|
|
|
Mat src(szSrc, CV_8UC(channels)), dst(sz, CV_8UC(channels));
|
|
cvtest::fillGradient<uint8_t>(src);
|
|
if(borderMode == BORDER_CONSTANT) cvtest::smoothBorder<uint8_t>(src, borderColor, 1);
|
|
Mat rotMat = getRotationMatrix2D(Point2f(src.cols/2.f, src.rows/2.f), 30., 2.2);
|
|
Mat warpMat(3, 3, CV_64FC1);
|
|
for(int r=0; r<2; r++)
|
|
for(int c=0; c<3; c++)
|
|
warpMat.at<double>(r, c) = rotMat.at<double>(r, c);
|
|
warpMat.at<double>(2, 0) = .3/sz.width;
|
|
warpMat.at<double>(2, 1) = .3/sz.height;
|
|
warpMat.at<double>(2, 2) = 1;
|
|
|
|
declare.in(src).out(dst);
|
|
|
|
TEST_CYCLE() warpPerspective( src, dst, warpMat, sz, interType, borderMode, borderColor );
|
|
|
|
SANITY_CHECK(dst, 1);
|
|
}
|
|
|
|
PERF_TEST_P( TestWarpPerspectiveNear_t, WarpPerspectiveNear,
|
|
Combine(
|
|
Values( Size(640,480), Size(1920,1080), Size(2592,1944) ),
|
|
InterType::all(),
|
|
BorderMode::all(),
|
|
Values( CV_8UC1, CV_8UC4 )
|
|
)
|
|
)
|
|
{
|
|
Size size;
|
|
int borderMode, interType, type;
|
|
size = get<0>(GetParam());
|
|
interType = get<1>(GetParam());
|
|
borderMode = get<2>(GetParam());
|
|
type = get<3>(GetParam());
|
|
Scalar borderColor = Scalar::all(150);
|
|
|
|
Mat src(size, type), dst(size, type);
|
|
cvtest::fillGradient<uint8_t>(src);
|
|
if(borderMode == BORDER_CONSTANT) cvtest::smoothBorder<uint8_t>(src, borderColor, 1);
|
|
int shift = static_cast<int>(src.cols*0.04);
|
|
Mat srcVertices = (Mat_<Vec2f>(1, 4) << Vec2f(0, 0),
|
|
Vec2f(static_cast<float>(size.width-1), 0),
|
|
Vec2f(static_cast<float>(size.width-1), static_cast<float>(size.height-1)),
|
|
Vec2f(0, static_cast<float>(size.height-1)));
|
|
Mat dstVertices = (Mat_<Vec2f>(1, 4) << Vec2f(0, static_cast<float>(shift)),
|
|
Vec2f(static_cast<float>(size.width-shift/2), 0),
|
|
Vec2f(static_cast<float>(size.width-shift), static_cast<float>(size.height-shift)),
|
|
Vec2f(static_cast<float>(shift/2), static_cast<float>(size.height-1)));
|
|
Mat warpMat = getPerspectiveTransform(srcVertices, dstVertices);
|
|
|
|
declare.in(src).out(dst);
|
|
declare.time(100);
|
|
|
|
TEST_CYCLE()
|
|
{
|
|
warpPerspective( src, dst, warpMat, size, interType, borderMode, borderColor );
|
|
}
|
|
|
|
SANITY_CHECK(dst, 1);
|
|
}
|
|
|
|
PERF_TEST_P( TestRemap, remap,
|
|
Combine(
|
|
Values( CV_8UC1, CV_8UC3, CV_8UC4, CV_32FC1 ),
|
|
Values( szVGA, sz1080p ),
|
|
InterTypeExtended::all(),
|
|
BorderMode::all(),
|
|
RemapMode::all()
|
|
)
|
|
)
|
|
{
|
|
int type = get<0>(GetParam());
|
|
Size size = get<1>(GetParam());
|
|
int interpolationType = get<2>(GetParam());
|
|
int borderMode = get<3>(GetParam());
|
|
int remapMode = get<4>(GetParam());
|
|
unsigned int height = size.height;
|
|
unsigned int width = size.width;
|
|
Mat source(height, width, type);
|
|
Mat destination;
|
|
Mat map_x(height, width, CV_32F);
|
|
Mat map_y(height, width, CV_32F);
|
|
|
|
declare.in(source, WARMUP_RNG);
|
|
|
|
update_map(source, map_x, map_y, remapMode, ((interpolationType & WARP_RELATIVE_MAP) != 0));
|
|
|
|
TEST_CYCLE()
|
|
{
|
|
remap(source, destination, map_x, map_y, interpolationType, borderMode);
|
|
}
|
|
|
|
SANITY_CHECK_NOTHING();
|
|
}
|
|
|
|
void update_map(const Mat& src, Mat& map_x, Mat& map_y, const int remapMode, bool relative )
|
|
{
|
|
for( int j = 0; j < src.rows; j++ )
|
|
{
|
|
for( int i = 0; i < src.cols; i++ )
|
|
{
|
|
switch( remapMode )
|
|
{
|
|
case HALF_SIZE:
|
|
if( i > src.cols*0.25 && i < src.cols*0.75 && j > src.rows*0.25 && j < src.rows*0.75 )
|
|
{
|
|
map_x.at<float>(j,i) = 2*( i - src.cols*0.25f ) + 0.5f ;
|
|
map_y.at<float>(j,i) = 2*( j - src.rows*0.25f ) + 0.5f ;
|
|
}
|
|
else
|
|
{
|
|
map_x.at<float>(j,i) = 0 ;
|
|
map_y.at<float>(j,i) = 0 ;
|
|
}
|
|
break;
|
|
case UPSIDE_DOWN:
|
|
map_x.at<float>(j,i) = static_cast<float>(i) ;
|
|
map_y.at<float>(j,i) = static_cast<float>(src.rows - j) ;
|
|
break;
|
|
case REFLECTION_X:
|
|
map_x.at<float>(j,i) = static_cast<float>(src.cols - i) ;
|
|
map_y.at<float>(j,i) = static_cast<float>(j) ;
|
|
break;
|
|
case REFLECTION_BOTH:
|
|
map_x.at<float>(j,i) = static_cast<float>(src.cols - i) ;
|
|
map_y.at<float>(j,i) = static_cast<float>(src.rows - j) ;
|
|
break;
|
|
} // end of switch
|
|
|
|
if( relative )
|
|
{
|
|
map_x.at<float>(j,i) -= static_cast<float>(i);
|
|
map_y.at<float>(j,i) -= static_cast<float>(j);
|
|
}
|
|
}
|
|
}
|
|
}
|
|
|
|
PERF_TEST(Transform, getPerspectiveTransform_1000)
|
|
{
|
|
unsigned int size = 8;
|
|
Mat source(1, size/2, CV_32FC2);
|
|
Mat destination(1, size/2, CV_32FC2);
|
|
Mat transformCoefficient;
|
|
|
|
declare.in(source, destination, WARMUP_RNG);
|
|
|
|
PERF_SAMPLE_BEGIN()
|
|
for (int i = 0; i < 1000; i++)
|
|
{
|
|
transformCoefficient = getPerspectiveTransform(source, destination);
|
|
}
|
|
PERF_SAMPLE_END()
|
|
|
|
SANITY_CHECK_NOTHING();
|
|
}
|
|
|
|
PERF_TEST(Transform, getPerspectiveTransform_QR_1000)
|
|
{
|
|
unsigned int size = 8;
|
|
Mat source(1, size/2, CV_32FC2);
|
|
Mat destination(1, size/2, CV_32FC2);
|
|
Mat transformCoefficient;
|
|
|
|
declare.in(source, destination, WARMUP_RNG);
|
|
|
|
PERF_SAMPLE_BEGIN()
|
|
for (int i = 0; i < 1000; i++)
|
|
{
|
|
transformCoefficient = getPerspectiveTransform(source, destination, DECOMP_QR);
|
|
}
|
|
PERF_SAMPLE_END()
|
|
|
|
SANITY_CHECK_NOTHING();
|
|
}
|
|
|
|
} // namespace
|