gpt4 book ai didi

image-processing - 关于使用卷积自定义高斯模糊的问题

转载 作者:太空宇宙 更新时间:2023-11-03 20:57:41 24 4
gpt4 key购买 nike

我正在编写一个函数 Conv2ByFFT() 来执行高斯模糊,这类似于 openCV api 中的 GaussianBlur()。但是当我比较 by function 和 GaussianBlur() api 之间的效果时,我发现前者没有后者那么“模糊”,我不知道为什么。

这是“正确的” this is the "correct" one

这是使用我的 Conv2ByFFT() 的结果 this is the result using my Conv2ByFFT()

这是一些代码

void Conv2ByFFT(const Mat& f,const Mat& g,Mat& result)
{
result.create(abs(f.rows-g.rows)+1,abs(f.cols-g.cols)+1,f.type());

//pad the images and get optimal FFT size
Size dftSize;
dftSize.width = getOptimalDFTSize(f.cols + g.cols - 1);
dftSize.height = getOptimalDFTSize(f.rows + g.cols - 1);

Mat tmpF(dftSize,f.type(),Scalar::all(0));
Mat tmpG(dftSize,g.type(),Scalar::all(0));

Mat roiF(tmpF, Rect(0,0,f.cols,f.rows));
f.copyTo(roiF);
Mat roiG(tmpG, Rect(0,0,g.cols,g.rows));
g.copyTo(roiG);

//perform Fourier Transform
dft(tmpF,tmpF,0,f.rows);
dft(tmpG,tmpG,0,g.rows);

//perform per-element multiplication of two Fourier spectrums
mulSpectrums(tmpF,tmpG,tmpF,0);

//perform inverse Fourier Transform
dft(tmpF,tmpF,DFT_INVERSE+DFT_SCALE,result.rows);

tmpF(Rect(0,0,result.cols,result.rows)).copyTo(result);
}

int main()
{
//read image
const char* imagename = "c:\\lena.bmp";
Mat img = imread(imagename);

//check image
if(img.empty())
{ fprintf(stderr, "Can not load image %s\n", imagename);
return -1;
}
if( !img.data )
return -1;

Mat src;

//convert the rgbimage into grayimage
cvtColor(img,src,CV_BGR2GRAY);

//save the grayimage
imwrite("lenagray.bmp",src);

//convert the image into float type
src.convertTo(src,CV_64FC1);


//******************************************************************************
// use GaussianBlur() in openCV
//******************************************************************************

//use Gaussian filter to blur the image
Mat dst = src.clone();
GaussianBlur(src,dst,Size(11,11),2);

//show and save the result
dst.convertTo(dst,CV_8U);
imshow("image",dst);
imwrite("lenablur.bmp",dst);

//******************************************************************************
// use GaussianBlur() in openCV
//******************************************************************************




//******************************************************************************
// use self-defining Conv2ByFFT()
//******************************************************************************

Mat result;

Mat gaussianFilter = getGaussianKernel(11,2,CV_64FC1);

//do the convolution to blur the image
Conv2ByFFT(src,gaussianFilter,result);

//show and save the result
result.convertTo(result,CV_8U);
//imshow("image1",result);
imwrite("lenablur1.bmp",result);

//******************************************************************************
// use self-defining Conv2ByFFT()
//******************************************************************************


cvWaitKey();
return 0;
}

最佳答案

getGaussianKernel 返回系数向量,而不是二维内核。

由于二维高斯核是可分离的,在卷积方法中,这个向量在两个方向上应用,与一次应用完整核具有相同的效果。

您的 FFT 函数只是将矢量与图像进行卷积。我想如果你仔细看,模糊只是在一个方向上。

您需要创建一个完整的二维高斯核,然后应用它。或者,我认为您可以利用可分离性并将向量应用两次。

关于image-processing - 关于使用卷积自定义高斯模糊的问题,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/7735719/

24 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com