KUAN滤波和Lee滤波的C++结合opencv实现

KUAN滤波和Lee滤波的opencv实现。

#include "stdafx.h"  
#include "cv.h"
#include "highgui.h"
#include <stdio.h>  
#include <string.h>  
#include <math.h>  
#include <opencv2/objdetect/objdetect.hpp>  
#include "opencv2/imgproc/imgproc.hpp"    
#include <opencv2/opencv.hpp>  
using namespace cv;
using namespace std;
double **shuju(Mat src)//将图像转化为二维数组
{
	int row = src.rows;
	int col = src.cols;
	double **a;
	a = new double*[row];
	for (int i = 0;i < row;i++)
	{
		a[i] = new double[col];
	}

	for (int i = 0;i < src.rows;i++)
		for (int j = 0;j < src.cols;j++)
		{
			a[i][j] = (double)src.at < uchar >(i, j);
		}
	return a;
}
void KUAN(Mat &src, Mat &dst)
{
	double **i1, **i2;
	i1=shuju(src);
	double ybar;
	double ystad;
	double ENL;
	double sx2;
	double xcap;
	for (int i = 1; i < src.rows - 1; ++i)
		for (int j = 1; j < src.cols - 1; ++j)
		{
			//ybar = (i1[i][j] + i1(i - 1, j - 1) + i1(i - 1, j) + i1(i, j - 1)
				//+ i1(i + 1, j) + i1(i, j + 1) + i1(i + 1, j + 1) + i1(i - 1, j + 1) + i1(i + 1, j - 1)) / 9;
			ybar = (i1[i - 1][j - 1] + i1[i - 1][j] + i1[i - 1][j + 1] + i1[i][j - 1] + i1[i][j] + i1[i][j + 1] + i1[i + 1][j - 1] + i1[i + 1][j] + i1[i + 1][j + 1])/9;
			/*ystad = sqrt(((src.at<double>(i, j) - ybar)*(src.at<double>(i, j) - ybar)
				+ (src.at<double>(i - 1, j - 1) - ybar)*(src.at<double>(i - 1, j - 1) - ybar)
				+ (src.at<double>(i - 1, j) - ybar)*(src.at<double>(i - 1, j) - ybar)

				+ (src.at<double>(i, j - 1) - ybar)*(src.at<double>(i, j - 1) - ybar)
				+ (src.at<double>(i + 1, j) - ybar)*(src.at<double>(i + 1, j) - ybar)
				+ (src.at<double>(i, j + 1) - ybar)*(src.at<double>