学习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;}