bu_te_wo_si_gao_tong.rar_图像增强_布特沃斯高通滤波器

%布特沃斯高通滤波器 I=imread('树.jpg'); I=im2double(I); A=rgb2gray(I); M=2*size(A,1);N=2*size(A,2); u=-M/2:(M/2-1);v=-N/2:(N/2-1); [U,V]=meshgrid(u, v); D=sqrt(U.^2+V.^2); D0=30;n=6; H=1./(1+(D0./D).^(2*n)); J=fftshift(fft2(A, size(H, 1), size(H, 2))); K=J.*H; T=ifft2(ifftshift(K)); T=T(1:size(I,1), 1:size(I, 2)); figure; subplot(221);imshow(A);title('原始图像'); subplot(222);imshow(T);title('布特沃斯高通滤波器滤波后图像'); subplot(223);imshow(A+T);title('高频增强图像'); [g w]=size(T); L=imresize(F,[floor(g/2) floor(w/2)]); L=imresize(T,[g w]); M=im2uint8(F); N=im2uint8(T); B=8; MAX6=2^B-1; MES6=sum(sum((M-N).^2))/(g*w); PSNR6=20*log10(MAX6/sqrt(MES6)) [g w]=size(A+T); L=imresize(F,[floor(g/2) floor(w/2)]); L=imresize(A+T,[g w]); M=im2uint8(F); N=im2uint8(A+T); B=8; MAX=2^B-1; MES=sum(sum((M-N).^2))/(g*w); PSNR=20*log10(MAX/sqrt(MES))