学习opencv第六章第五题（6_5）

#include "cv.h"
#include "highgui.h"
#include <IOSTREAM.H>

#pragma comment(lib,"cv.lib")
#pragma comment(lib,"highgui.lib")
#pragma comment(lib,"cxcore.lib")

int main(int argc, char** argv)
{
	const CvSize size = cvSize(512,512);
	const int aperture[] = {3,5,7};
	IplImage* src = cvCreateImage(size,8,1);

	cvZero(src);

	cvLine(src,cvPoint(0,0),cvPoint(size.width -1, size.height-1),CV_RGB(255,255,255),3);

	cvNamedWindow("src",CV_WINDOW_AUTOSIZE);
	cvShowImage("src",src);

	IplImage*  deriv_x = cvCreateImage(size,IPL_DEPTH_32F,1); //x方向导数形成的图像
	IplImage*  deriv_y = cvCreateImage(size,IPL_DEPTH_32F,1); //y方向导数形成的图像

	IplImage* angle = cvCreateImage(size, IPL_DEPTH_32F,1);  //进行极坐标变换后形成的角度图像
	IplImage* magnitude = cvCreateImage(size,IPL_DEPTH_32F,1);//进行极坐标变换后形成的幅度图像
	IplImage* mask = cvCreateImage(size,IPL_DEPTH_8U,1);//进行比较过后（选取幅度超过阈值的点）形成的模板图像
    
	cvNamedWindow("src",CV_WINDOW_AUTOSIZE);
	cvNamedWindow("mask",CV_WINDOW_AUTOSIZE);

	for (int i = 0; i<sizeof(aperture)/sizeof(aperture[0]); i++)
	{
		cvSobel(src,deriv_x,1,0,aperture[i]);//x求导
		cvSobel(src,deriv_y,0,1,aperture[i]);//y求导

		cvCartToPolar(deriv_x,deriv_y,magnitude,angle,1);//极坐标变换

		cvSave("x.xml",deriv_x);
		cvSave("y.xml",deriv_y);
		cvSave("magnitude.xml",magnitude);
		cvSave("angle.xml",angle);

		double maxmag, minmag, maxangle, minangle;
		cvMinMaxLoc(magnitude,&minmag,&maxmag);
		cvMinMaxLoc(angle,&minangle,&maxangle);

		cout<<"magnitude: max="<<maxmag<<" ,min="<<minmag<<endl;
		cout<<"angle: max="<<maxmag<<" ,min="<<minangle<<endl;

		//生成模板
		cvCmpS(magnitude,maxmag*3/4,mask,CV_CMP_GT);//寻找最大幅值附近的点
		CvScalar scalar = cvAvg(magnitude,mask);
		cout<<"aperture="<<aperture[i]<<", line angle ="<<scalar.val[0]<<endl;

		cvShowImage("src",src);
		cvShowImage("mask",mask);

		if (27 == cvWaitKey(0))
		{
			continue;
		}

	}

	cvReleaseImage(&src);
	cvReleaseImage(&deriv_x);
	cvReleaseImage(&deriv_y);
	cvReleaseImage(&magnitude);
	cvReleaseImage(&angle);
	cvReleaseImage(&mask);
	cvDestroyAllWindows();
	return 1;
}