opencv 笛卡尔坐标转换为极坐标 cvLogPolar的应用一

对于二维图像，Log-Polar转换表示从笛卡尔坐标到极坐标的变换。
OpenCV中用于这个变换的函数是cvLogPolar，其函数原型如下：
cvLogPolar：把图像映射到极指数空间

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(x^2+y2))； phi=atan(y/x)

此函数模仿人类视网膜中央凹视力，并且对于目标跟踪等可用于快速尺度和旋转变换不变模板匹配。
参考链接：https://blog.csdn.net/qiuqchen/article/details/26228507
参考链接：https://blog.csdn.net/gdut2015go/article/details/46397999
https://blog.csdn.net/zcg1942/article/details/80886999

代码如下：

#include   
#include   

int main(int argc, char** argv)
{
    IplImage* src = cvLoadImage("C:\\分割的虹膜图像.bmp", 1);
    IplImage* dst = cvCreateImage(cvSize(256,256), 8, 3);
    IplImage* src2 = cvCreateImage(cvGetSize(src), 8, 3);
    cvLogPolar(src, dst, cvPoint2D32f(src->width / 2, src->height / 2), 40, 
        CV_INTER_LINEAR + CV_WARP_FILL_OUTLIERS);  //把图像映射到极指数空间  
    //CV_WARP_FILL_OUTLIERS 填充输出图像所有像素，如果这些点有和外点对应的，则置零。  
    cvLogPolar(dst, src2, cvPoint2D32f(src->width/2 , src->height/2 ), 40, 
        CV_INTER_LINEAR + CV_WARP_FILL_OUTLIERS + CV_WARP_INVERSE_MAP);
    //CV_WARP_INVERSE_MAP表示矩阵由输出图像到输入图像的逆变换，并且因此可以直接用于像素插值。否则，函数从map_matrix中寻找逆变换。  
    cvNamedWindow("log-polar", 1);
    cvShowImage("log-polar", dst);
    cvNamedWindow("inverse log-polar", 1);
    cvShowImage("inverse log-polar", src2);
    cvWaitKey();
    return 0;
}

图像变形严重，总结了一下是代码里面那个点选取的有失误，应该用虹膜瞳孔的圆心作为坐标点，然后图像画质不清楚，应该是像素的插值填充没选好