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(x2+y2)); phi=atan(y/x)
此函数模仿人类视网膜中央凹视力,并且对于目标跟踪等可用于快速尺度和旋转变换不变模板匹配。
参考链接:https://blog.****.net/qiuqchen/article/details/26228507
参考链接:https://blog.****.net/gdut2015go/article/details/46397999
https://blog.****.net/zcg1942/article/details/80886999
代码如下:
#include <highgui.h>
#include <cv.h>
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;
}
图像变形严重,总结了一下是代码里面那个点选取的有失误,应该用虹膜瞳孔的圆心作为坐标点,然后图像画质不清楚,应该是像素的插值填充没选好