matlab对图像进行频率域图像增强，计算并画出此图像的中心化频率谱，用拉普拉斯算子对此图像锐化

clear; close all; clc A=0; a=imread('E:\DIP\PremierZhou.png');%读入图像 b=reshape(a,591,2760);%转换为二维矩阵 bb=double(b); figure,imshow(a); %画出原始图像 title('原始图像'); %原图像频谱 ffta=fft2(bb);%对图像进行二维FFT变换 fftac=fftshift(ffta);%频谱中心化 if(A==0) BB=abs(ffta); BB=BB/max(max(BB)); BB=256*log2(BB+1); figure,imshow(BB); title('原图的频谱'); end if(A==0) BBB=abs(fftac); BBB=BBB/max(max(BBB)); BBB=256*log2(BBB+1); figure,imshow(BBB); title('中心化的频谱'); end Laplace=[0 1 0 1 -4 1 0 1 0 ]; im1=double(a(:,:,1)); laplaceimage=conv2(im1,Laplace); figure,imshow(uint8(laplaceimage)); title('Laplace图像'); figure,imshow(uint8(im1)); title('laplace 锐化后的图像'); h=-fspecial('laplacian');%h为拉普拉斯滤波器 g=imfilter(b,h); n=imadd(b,immultiply(g,1));%原图像与拉普拉斯图像叠加 N=reshape(n,591,920,3); figure,imshow(N); title('锐化增强后的图像');