Dog(Difference of Gaussian: 高斯差分)

///////////////////////////////////////////////////////////////////// //Author��Michael.Chen //Description: Iplement the DoG(Difference of Gaussian). // It's just a demo to show the algorithm. // You can change the code as you like. // // I use OpenCV2.1 // You have to config the environment firs. // // // I would to discuss with you. You can contact me by QQ or Email. // Thanks for your reading. // //Date: 2011-09-25 //QQ: 1241159842 //Email: 1241159842@qq.com ///////////////////////////////////////////////////////////////////// #include "cv.h" #include "cxcore.h" #include "highgui.h" int main(void) { IplImage* pSrc = cvLoadImage("e:\\lena.jpg"); IplImage* pGray = cvCreateImage(cvGetSize(pSrc), IPL_DEPTH_8U, 1); cvCvtColor(pSrc, pGray, CV_RGB2GRAY); IplImage* pGray32 = cvCreateImage(cvGetSize(pSrc), IPL_DEPTH_32F, 1); cvCvtScale(pGray, pGray32, 1.00 / 255.00); IplImage* pGaussian1 = cvCreateImage(cvGetSize(pSrc), IPL_DEPTH_32F, 1); IplImage* pGaussian2 = cvCreateImage(cvGetSize(pSrc), IPL_DEPTH_32F, 1); double sigma = 1.6; double k = 1.414; cvSmooth(pGray32, pGaussian1, CV_GAUSSIAN, 0, 0, sigma, sigma); cvSmooth(pGaussian1, pGaussian2, CV_GAUSSIAN, 0, 0, sigma * k, sigma * k); IplImage* pDog1 = cvCreateImage(cvGetSize(pSrc), IPL_DEPTH_32F, 1); cvSub(pGaussian2, pGaussian1, pDog1); IplImage* pShow = cvCreateImage(cvGetSize(pSrc), IPL_DEPTH_8U, 1); cvNormalize(pDog1, pShow, 0, 255, CV_MINMAX); cvNamedWindow("Dog"); cvShowImage("Dog", pShow); cvWaitKey(); return 0; }