SRAD算法C++实现

SRAD算法平滑去噪的效果非常理想，能在去除噪声的同时保持或增强目标边缘的清晰度。本程序使用OpenCV1.0 在VC6.0下运行成功。这个程序是我读研时候的一个学弟写的，这方面的资料比较少，所以拿出来和大家分享，希望对您有帮助！

#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）;// 创建窗口
    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）;
    // 将图像像素变为32位的单精度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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       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