电子科大图像处理第二次作业

第二次作业对一幅灰度图像：(1)计算并画出此图像的中心化频率谱。(2)分别用高斯低通和高斯高通滤波器对图像进行频域处理。(3)用频域拉普拉斯算子对此图像进行锐化处理。代码：clear;clc;Data=imread('C:\Pdog.JPG');DataGray=rgb2gray(Data);figure(1),imshow(Data);title('原始图像');%**************计算并画出此图像的中心化频率谱*************Data1=double(DataGray);FFT2=fft2(Data1);FFTcenter=fftshift(FFT2);%频谱中心化FFT2abs=abs(FFT2);FFTresult=256*log2(FFT2abs/max(max(FFT2abs))+1);figure(2),subplot(1,2,1);imshow(FFTresult),title('原图频谱');FFTc_abs=abs(FFTcenter);FFTc_result=256*log2(FFTc_abs/max(max(FFTc_abs))+1);subplot(1,2,2);imshow(FFTc_result),title('中心化频谱');%*******分别用低通滤波和高通滤波对此图像进行频域处理****[m,n]=size(FFTcenter);x_center=round(m/2);y_center=round(n/2);d=10;%半径取10LF=FFTcenter;HF=FFTcenter;%**************低通滤波器*************fori=1:m;forj=1:ndistance=sqrt((i-x_center)^2+(j-y_center)^2);ifdistance=dflag=1;elseflag=0;endLF(i,j)=flag*FFTcenter(i,j);endend%逆变换转换成对应图像LF=uint8(real(ifft2(ifftshift(LF))));figure(3),subplot(1,2,1);imshow(LF);title('低通滤波后图像');%**************高通滤波器*************fori=1:m;forj=1:ndistance=sqrt((i-x_center)^2+(j-y_center)^2);ifdistancedflag=1;elseflag=0;endHF(i,j)=flag*FFTcenter(i,j);endendHF=uint8(real(ifft2(ifftshift(HF))));subplot(1,2,2),imshow(HF);title('高通滤波后图像');%*****************用拉普拉斯算子对此图像锐化*************Laplace=[0-10;-14-1;0-10];LaplaceImage=conv2(Data1,Laplace,'same');figure(4),subplot(1,2,1);imshow(uint8(LaplaceImage));title('Laplace图像');DataLap=imadd(Data1,immultiply(LaplaceImage,1));%原图像与拉普拉斯图像叠加subplot(1,2,2),imshow(uint8(DataLap));title('锐化增强后的图像');截图：