图像边缘锐化和标准图像压缩过程

《图像边缘锐化和标准图像压缩过程》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、1.实现两个高频加强的图像边缘锐化的方法，原始图像是灰度图像，给出高频加强公式、程序和实验效果。(注意程序不能是直接一个函数的调用，要把变换过程和公式体现出来。)(1)理想高通滤波器functionout=imideal(I,freq)%I参数输入的灰度图像%freq参数高通滤波器的截止频率%返回值：out-指定的理想高通滤波器f=fft2(I);%傅立叶变换f=fftshift(f);%转换数据矩阵[M，N]=size(I);out=ones(M,N);fori=l:Mforj=l:Nif(sqrt(((i-M/2)A2+(j-N/2)A2))

2、,j)=0;endendendout=fout;%对应元素相乘实现频域滤波(2)巴特沃斯高通滤波器巴特沃斯高频加强公式H(u,v)=a+b/(1+D0/D(u,v))2nfunctionout=imibtws(I,dO)%I参数输入的灰度图像%d0参数截止频率距远点距离%返冋值：out-指定的巴特沃斯高通滤波器f=fft2(I);%傅立叶变换f=fftshift(f)；%转换数据矩阵[M,N]=size(I);out=zeros(M,N);nn=2;fori=l:Mforj=l:Nd=sqrt((i-M/2)A2+(j-N/2)A2);if(d)out(i,j)=l/

3、(l+0.414*(d0/d)A(2*nn));1.实现两个高频加强的图像边缘锐化的方法，原始图像是灰度图像，给出高频加强公式、程序和实验效果。(注意程序不能是直接一个函数的调用，要把变换过程和公式体现出来。)(1)理想高通滤波器functionout=imideal(I,freq)%I参数输入的灰度图像%freq参数高通滤波器的截止频率%返回值：out-指定的理想高通滤波器f=fft2(I);%傅立叶变换f=fftshift(f);%转换数据矩阵[M，N]=size(I);out=ones(M,N);fori=l:Mforj=l:Nif(sqrt(((i-M/2)A

4、2+(j-N/2)A2))

5、N/2)A2);if(d)out(i,j)=l/(l+0.414*(d0/d)A(2*nn));endendendout=f.*out;%对应元素相乘实现频域滤波(1)归一化函数functionout=guiyi(I)%guiyi图像频域显示前的预处理%I参数图像频域%返回值：out-图像频域显示I=ifftshift(I);I=ifft2(I);I=uint8(real(I));out=I;(2)主函数closeallclearclcI=imread(,LENA.bmp,);figure(1);subplot(1,2,1);imshow(I);title(’原图像’

6、)；f=fft2(I);%做fft变换f=fftshift(f);%将零点移到中心t=log(l+abs(f));%对幅值做对数变换，以压缩动态范围figure(1);subplot(l，2,2);imshow(t，[]);titlef原图像FFT);%t是doublearray,是浮点数，需要[].%理想高通过滤h=imideal(I，60);out=f+0.9*h;%归一化处理out=guiyi(out);figure(2);subplot(1，2,1);imshow(out);title(’理想高频加强锐化’)；h=guiyi(h);figure(2);subp

7、lot(l，2,2);imshow(h);titlef理想高频加强锐化’);%巴特沃斯高通过滤器h=imibtws(I,15);out=f+0.9*h;%高频加强锐化%归一化处理out=guiyi(out);figure(3);subplot(l，2，l);imshow(out);title('巴特沃斯高频加强锐化•);h=guiyi(h);figure(3);subplot(l,2,2);imshow(h);title(’巴特沃斯高频加强滤波');(5)实验结果原图像原图像FFT理想高频加強锐化理想高频加強锐化巴特沃斯高频加強锐化巴特沃斯高频加強滤