对于二维图形,Log-polar转换表示从笛卡尔坐标到极坐标的变化,广泛应用在计算机视觉中。此函数模仿人类视网膜中央凹视力,并且对于目标跟踪等可用于快速尺度和旋转变换不变模板匹配。
本例程实现极坐标变化,并反转。
----------------------------------------------------------------------------------------------
LogPolar
把图像映射到极指数空间
void cvLogPolar( const CvArr* src, CvArr* dst,
CvPoint2D32f center, double M,
int flags=CV_INTER_LINEAR+CV_WARP_FILL_OUTLIERS );
src
输入图像。
dst
输出图像。
center
变换的中心,输出图像在这里最精确。
M
幅度的尺度参数,见下面公式。
flags
插值方法和以下选择标志的结合
· CV_WARP_FILL_OUTLIERS-填充输出图像所有像素,如果这些点有和外点对应的,则置零。
· CV_WARP_INVERSE_MAP- 表示矩阵由输出图像到输入图像的逆变换,并且因此可以直接用于像素插值。否则,函数从map_matrix中寻找逆变换。
fillval
用于填充外点的值。
函数cvLogPolar用以下变换变换输入图像:
正变换 (CV_WARP_INVERSE_MAP 未置位):
dst(phi,rho)<-src(x,y)
逆变换 (CV_WARP_INVERSE_MAP 置位):
dst(x,y)<-src(phi,rho),
这里,
rho=M*log(sqrt(x2+y2))
phi=atan(y/x)
----------------------------------------------------------------------------------------------
/*code*/
#include <opencv/cv.h>
#include <opencv/highgui.h>
#include <opencv2/opencv.hpp>
using namespace cv;
int main(int argc,char** argv)
{
CvCapture *capture=0;
capture=cvCreateFileCapture("../test.avi");
if(!capture)
{
return -1;
}
IplImage *bgr_frame=cvQueryFrame(capture);
double fps=cvGetCaptureProperty(capture,CV_CAP_PROP_FPS);
CvSize size=cvSize((int)cvGetCaptureProperty(capture,CV_CAP_PROP_FRAME_WIDTH),(int)cvGetCaptureProperty(capture,CV_CAP_PROP_FRAME_HEIGHT));
CvVideoWriter *writer=cvCreateVideoWriter("../heibai.avi",CV_FOURCC('M','J','P','G'),fps,size);
IplImage *logpolar_frame=cvCreateImage(size,IPL_DEPTH_8U,3);
IplImage *logpolar_frame1=cvCreateImage(size,IPL_DEPTH_8U,3);
cvNamedWindow("Video test",CV_WINDOW_AUTOSIZE);
while((bgr_frame=cvQueryFrame(capture))!=NULL)
{
cvLogPolar(bgr_frame,logpolar_frame,cvPoint2D32f(bgr_frame->width/2,bgr_frame->height/2),50,CV_INTER_LINEAR+CV_WARP_FILL_OUTLIERS);//把图像映射到极指数空间
//CV_WARP_FILL_OUTLIERS -填充输出图像所有像素,如果这些点有和外点对应的,则置零。
//cvLogPolar(bgr_frame,logpolar_frame1,cvPoint2D32f(bgr_frame->width/2,bgr_frame->height/2),10,CV_INTER_LINEAR + CV_WARP_INVERSE_MAP);
//CV_WARP_INVERSE_MAP - 表示矩阵由输出图像到输入图像的逆变换,并且因此可以直接用于像素插值。否则,函数从map_matrix中寻找逆变换。
cvShowImage("Video test",logpolar_frame);
cvWaitKey((int)fps);
cvWriteFrame(writer,logpolar_frame);
}
cvReleaseVideoWriter(&writer);
cvReleaseImage(&logpolar_frame);
cvReleaseImage(&logpolar_frame1);
cvReleaseImage(&bgr_frame);
cvReleaseCapture(&capture);
cvDestroyWindow("Video test");
return 0;
}
/*result*/
依旧水果图
log-polar:
对于这一转换在实际应用中有何用途不清楚!