图像直接滤波 opencv下的问题显示为黑色

2013-04-30 02:36:09 327 阅读 AI 写同款 · GEO 优化AI 选题数据监测一键分发 ›

【OpenCV篇】OpenCV——03day.图像预处理（2）本文介绍了图像处理中的常用技术，包括多种插值方法（最近邻、双线性、双三次等）、边缘填充方式（边界复制、反射等）、图像矫正、掩膜操作、ROI切割、水印添加以及噪声消除滤波方法（均值、高斯、中值滤波等）。这些技术涵盖了图像处理的核心操作，为图像分析、增强和修复提供了基础工具和方法论支持。文章通过系统分类和详细说明，帮助读者快速掌握各类图像处理技术的原理与应用场景。阅读详情

#include "cv.h"
#include "highgui.h"

using namespace cv;
using namespace std;

Mat scr_img ;
Mat gray_img ,blur_img;

Mat InverseFliter();

Mat InverseFliter()
{

Mat padded;
int dft_h = getOptimalDFTSize(blur_img.rows);
int dft_w = getOptimalDFTSize(blur_img.cols);
copyMakeBorder(blur_img,padded,0,dft_h-blur_img.rows,0,dft_w-blur_img.cols,BORDER_CONSTANT,Scalar::all(0));
Mat planes[] = {Mat_<float>(padded) , Mat::zeros(padded.size(),CV_32F)};
Mat complexI;
merge(planes,2,complexI);
dft(complexI,complexI);
double a,b,c,d;
int height = complexI.rows;
int width =complexI.cols;
int ch = complexI.channels();
Mat H_Function = (Mat_<float>(dft_h, dft_w ));
Mat h_fuction;
for(int j=0;j <dft_h;j++)
{
uchar *H_ptr = H_Function.ptr<uchar>(j);
for(int i=0; i<dft_w ;i++)
{
//for(int k=0; k<ch;k++)
//{
if(i < 5 && j< 5)
{
H_ptr[i ] =0.04;
}
else
{
H_ptr[i] =0.0;
}
//}
}
}
Mat planes1[] = { H_Function , Mat::zeros(H_Function.size(),CV_32F)};
merge(planes1,2,h_fuction);

dft(h_fuction,h_fuction);
Mat dst_img;
dst_img.create(dft_h, dft_w, CV_64FC2);
for(int row=0; row <height;row++)
{
uchar *Scr_ptr =complexI.ptr<uchar>(row);
uchar *H_ptr1 = h_fuction.ptr<uchar>(row);
double *Dst_ptr = dst_img.ptr<double>(row);
for (int col = 0; col< width;col++)
{

a = Scr_ptr[col*ch];
b = Scr_ptr[col*ch +1];
c = H_ptr1[col*ch];
d = H_ptr1[col*ch+1];
if(c*c +d*d >1e-3)
{
Dst_ptr[col *ch] = (a*c +b*d)/(c*c + d*d);
Dst_ptr[col *ch +1] = (b*c-a*d)/(c*c+d*d);
}

}
}
idft(dst_img, dst_img); //或者 dft(dst_img,dst_img,CV_DXT_INVERSE_SCALE);

split(dst_img,planes1);
Mat idft_img2 = planes1[0];
return idft_img2;

int main()
{
scr_img = imread("lena.jpg",1);
if(scr_img.empty ())
{
cout<<"Miss the Image File"<<endl;
return -1;
}
cvtColor(scr_img,gray_img,CV_BGR2GRAY);
blur(gray_img,blur_img,Size(4,4));
Mat scr = InverseFliter();
//s_img.convertTo(blur_img,CV_8UC3,1,0);
imshow("window2",scr);
waitKey(0);
}

...全文