rst实验

前两天用到这块径向对称变换，就用vs2010+opencv2.4.2实验了下，效果和原文差不多，可以

图：

                      (radius={1,3,5,7,9}的结果}

       

 (radius=3,下面的是相应colormap)

下面这段代码是我对3种存取元素操作的对比： 

分别是ptr指针， at方法， 直接指针读取。测试图片就是上图， 运行半径是{1,3,5,7,9},belta分别是0和0.2

	for(int i=0;i<On.rows;i++)
	{
		//* 52ms-g.jpg belta=0 || 32ms-g.jpg belta=0.2
		float *magVal	= mag.ptr<float>(i);
		float *gradxVal	= unitGradX.ptr<float>(i);
		float *gradyVal	= unitGradY.ptr<float>(i);
		/**/
		for(int j=0;j<On.cols;j++)
		{
			/* 56ms-g.jpg belta=0|| 34ms-g.jpg belta=0.2

			float magVal=mag.at<float>(i,j);
			float gradxVal=unitGradX.at<float>(i,j);
			float gradyVal=unitGradY.at<float>(i,j);
			/**/
			/* 39ms-g.jpg belta=0 || 33ms-g.jpg bleta=0.2
			float magVal   = ((float*)(mag.data))[i*mag.cols+j];
			float gradxVal = ((float*)(unitGradX.data))[i*unitGradX.cols+j];
			float gradyVal = ((float*)(unitGradY.data))[i*unitGradY.cols+j];
			/**/
		}
	}

可见，用指针直接操作还是很有诱惑力的···

后来我把程序中无关的显示输出都去掉，速度可以达到论文中说的程度