вопросы о самоопределенном гауссовом размытии с использованием свертки - PullRequest
Новая
       31

вопросы о самоопределенном гауссовом размытии с использованием свертки

4 голосов
/ 12 октября 2011

Я пишу функцию Conv2ByFFT (), чтобы сделать размытие по Гауссу, которое похоже на GaussianBlur () в API openCV.Но, сравнивая эффекты между функциями и API GaussianBlur (), я обнаружил, что первый не такой «размытый», как второй, и я не знаю, почему.

это «правильный»this is the

это результат использования моего 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;
}

1 Ответ

5 голосов
/ 12 октября 2011

getGaussianKernel возвращает вектор коэффициентов, а не двумерное ядро.

Поскольку двумерное гауссовское ядро ​​отделимо, в методе свертки этот вектор применяется в обоих направлениях, что имеет тот же эффект, что и полное ядро ​​сразу.Функция просто сворачивает вектор с изображением.Я думаю, если вы посмотрите внимательно, размытие только в одном направлении.

Вам нужно создать полное двумерное гауссовское ядро ​​и применить его.В качестве альтернативы, я думаю, вы можете использовать отделимость и применить вектор дважды.

...