#include?“cv.h“
#include?“highgui.h“
#include?“cxcore.h“
#include?
#include?
#include?
using?namespace?std;

#pragma?comment（lib?“cv.lib“）
#pragma?comment（lib?“cxcore.lib“）
#pragma?comment（lib?“highgui.lib“）

const?float?eps?=?0.00000000000000001f;
void?CvMat_Sqrt（const?CvArr?*src?const?CvArr?*dst）;

IplImage*?Srad（IplImage*?src?float?delta_t?float?q0?float?rho?int?option?int?ITER）;

int?main（int?argc?char**?argv）
{
????IplImage*?pImg?*pDst;
????if?（（pImg?=?cvLoadImage（“P6.jpg“?0））?==?0）//?如果載入圖片失敗，退出
	{
?????????return?-1;
	}
	pDst?=?Srad（pImg?1.0/7?1.0f?1.0f?2?30）;

	cvNamedWindow（“SrcImage“?1）;//?創(chuàng)建窗口
????cvShowImage（“SrcImage“?pImg）;//?顯示圖片
	cvNamedWindow（“DstImage“?1）;
????cvShowImage（“DstImage“?pDst）;
????cvWaitKey（0）;?
????cvDestroyWindow（“SrcImage“）;//?銷毀窗口
	cvDestroyWindow（“DstImage“）;//?銷毀圖片指針
????cvReleaseImage（&pImg）;
	cvReleaseImage（&pDst）;

????return?0;		
}

void?InitConstMat（CvMat?*mat?float?num）//?初始化元素全部相同的矩陣
{
	int?h?w;
	for?（h=0;?hrows;?h++）
		for?（w=0;?wcols;?w++）
		{
			（mat->data.fl+h*mat->cols）[w]?=?num;
		}
}
void?CvMat_Sqrt（const?CvMat?*src?const?CvMat?*dst）
{
	int?i?j?r?c;
	r?=?src->rows;
	c?=?src->cols;
	dst?=?cvCreateMat（r?c?CV_32FC1）;
	float?*psrc;
	float?*pdst;
	for?（i=0;?i	{
		psrc?=?src->data.fl?+?i*c;
		pdst?=?dst->data.fl?+?i*c;
		for?（j=0;?j		{
			pdst[j]?=?（float）sqrt（psrc[j]）;
		}
	}
}

IplImage*?Srad（IplImage*?src?float?delta_t?float?q0?float?rho?int?option?int?ITER）
{
????CvMat*?Im?=?cvCreateMat（src->height?src->width?CV_32FC1）;
????//?cvConvert（srcIm）;
????//?將圖像像素變?yōu)?2位的單精度float型
	uchar?*pchar;
	int?i?j;
	for?（i=0;?iheight;?i++）
	{
		pchar?=?（uchar*）src->imageData?+?i*（src->widthStep）;?
		for?（j=0;?jwidth;?j++）
		{
			（Im->data.fl+i*src->width）[j]?=?（float）pchar[j]?/?255;
		}
	?}
	//?初始化
????CvMat*?gradRx?*gradRy?*gradLx?*gradLy;?
	gradRx?=?cvCreateMat（Im->height+2?Im->width+2?CV_32FC1）;
	gradRy?=?cvCreateMat（Im->height+2?Im->width+2?CV_32FC1?）;
	gradLx?=?cvCreateMat（Im->height+2?Im->width+2?CV_32FC1）;
	gradLy?=?cvCreateMat（Im->height+2?Im->width+2?CV_32FC1?）;
????cvZero（gradRx）;?cvZero（gradRy）;	cvZero（gradLx）;	cvZero（gradLy）;
	float?k1[9]?=?{0?0?0?0-1?1?0?0?0};?
	float?k2[9]?=?{0?1?0?0?-1?0?0?0?0};?
	float?k3[9]?=?{0?0?0?1?-1?0?0?0?0};?
	float?k4[9]?=?{0?0?0?0?-1?0?0?1?0};
	float?mm1[9]?=?{0?0?0?0?0?1?0?0?0};
	float?mm3[9]?=?{0?1?0?0?0?0?0?0?0};
	CvMat?kRx?=?cvMat（3?3?CV_32FC1?k1）;
	CvMat?kRy?=?cvMat（3?3?CV_32FC1?k2）;
	CvMat?kLx?=?cvMat（3?3?CV_32FC1?k3）;
	CvMat?kLy?=?cvMat（3?3?CV_32FC1?k4）;
	CvMat?km1?=?cvMat（3?3?CV_32FC1?mm1）;
	CvMat?km3?=?cvMat（3?3?CV_32FC1?mm3）;

	CvMat?*epsMat;
	epsMat?=?cvCreateMat（Im->height+2?Im->width+2?CV_32FC1）;
	InitConstMat（epsMat?eps）;
	CvMat?*oneMat;
	oneMat?=?cvCreateMat（Im->height+2?I

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????6805??2011-12-25?21:50??SRAD+VC6.0\Debug\P6.jpg

?????文件?????303190??2013-06-10?16:12??SRAD+VC6.0\Debug\Srad.exe

?????文件?????428920??2013-06-10?16:12??SRAD+VC6.0\Debug\Srad.ilk

?????文件??????64090??2013-06-10?16:12??SRAD+VC6.0\Debug\Srad.obj

?????文件????5875464??2013-06-10?16:12??SRAD+VC6.0\Debug\Srad.pch

?????文件?????623616??2013-06-10?16:12??SRAD+VC6.0\Debug\Srad.pdb

?????文件?????214016??2013-06-10?16:12??SRAD+VC6.0\Debug\vc60.idb

?????文件?????151552??2013-06-10?16:12??SRAD+VC6.0\Debug\vc60.pdb

?????文件???????8095??2014-05-22?09:32??SRAD+VC6.0\Srad.cpp

?????文件???????3462??2013-06-10?16:13??SRAD+VC6.0\Srad.dsp

?????文件????????533??2012-06-15?12:31??SRAD+VC6.0\Srad.dsw

?????文件??????41984??2013-06-10?16:13??SRAD+VC6.0\Srad.ncb

?????文件??????48640??2013-06-10?16:13??SRAD+VC6.0\Srad.opt

?????文件???????1342??2013-06-10?16:12??SRAD+VC6.0\Srad.plg

?????目錄??????????0??2014-05-22?09:53??SRAD+VC6.0\Debug

?????目錄??????????0??2014-05-22?09:53??SRAD+VC6.0

-----------?---------??----------?-----??----

??????????????7771709????????????????????16