光学图像处理作业----图像增强

《光学图像处理作业----图像增强》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、光学图像处理作业学号：05S001127姓名：马泽峰matlab实现对一幅图像增强处理：1.灰度的对数变换；2.保持边界的滤波(有选择保持边缘平滑法)；3.LOG算子实现高通滤波处理；4.作付立叶变换与DCT变换；5.实现频域内的梯度高通滤波。程序代码：A=imread('gray.bmp');subplot(241)imshow(A);title('变换前');[m,n]=size(A);M=medfilt2(A,[55]);fori=1:nM(1,i)=A(1,i);M(m,i)=A(m,i);endforj=1:mM(j,1)=A(j,1);M(j,n

2、)=A(j,n);endsubplot(242)imshow(M);title('保持边界中值虑波');D=dct2(A);D(abs(DCTA)<10)=0;D1=idct2(D)/255;subplot(243),imshow(D1);title('DCT变换');M=im2double(A);H=[-2-4-4-4-2,-4080-4,-48248-4,-4080-4,-2-4-4-4-2];LG=conv2(M,H,'same');subplot(244)imshow(mat2gray(LG));title('LOG算子高通滤波');L=double

3、(A);c=255/log(256);fori=1:mforj=1:ng(i,j)=c*log(L(i,j)+1);endendsubplot(245)imshow(mat2gray(g));title('灰度对数变换');F=fftshift(fft2(double(A)));F(abs(F)<10)=0;F1=ifftshift(F);F1=ifft2(F1);F2=double(real(F1));subplot(246),imshow(mat2gray(F2));title('傅里叶变换');f=fftshift(fft2(double(A)));[

4、N1,N2]=size(f);d0=6;d1=3;n1=fix(N1/2);n1=fix(N2/2);fori=1:N1forj=1:N2d=sqrt((i-n1)^2+(i-n1)^2);ifd<=d1h=0;elseifd>d1&d<=d0h=(d-d1)/(d0-d1);elseh=1;endr(i,j)=h*f(i,j);endendr=ifftshift(r);r1=ifft2(r);r2=double(real(r1));subplot(247)imshow(mat2gray(r2));title('梯度高通滤波')变换前后图像：