MATLAB选修论文二.docx
《MATLAB选修论文二.docx》由会员分享,可在线阅读,更多相关《MATLAB选修论文二.docx(13页珍藏版)》请在冰豆网上搜索。
![MATLAB选修论文二.docx](https://file1.bdocx.com/fileroot1/2023-6/15/e3607dc6-a386-4b1d-83ad-2b3f9a12539b/e3607dc6-a386-4b1d-83ad-2b3f9a12539b1.gif)
MATLAB选修论文二
2014111079刘青青测控二班
二、利用混合空间法增强锐化图像仿真
将原始图像
增强,最终增强为如下类似图像
。
可参考如下过程
首先用拉普拉斯法突出图像中的小细节,然后用梯度法突出其边。
平滑过的梯度图像将用于掩蔽拉普拉斯图像。
最后,用灰度变换来增强图像的灰度动态范围
clearall;clc;
JMa=imread('C:
\Users\qingqing\Pictures\tuxing.tif')
JMa=im2double(JMa);
%JMa=rgb2gray(JMa);
%a.原始图像
subplot(2,4,1);
imshow(JMa);
title('a:
原始图像');
%b.使用模板为[-1,-1,-1;-1,8,-1;-1,-1,-1]的滤波器对原图像进行拉普拉斯操作
h=[-1,-1,-1;-1,8,-1;-1,-1,-1];
JMb=imfilter(JMa,h);
subplot(2,4,2);imshow(JMb);
title('b:
拉普拉斯操作后图像');
%c.原图a和图b相加
subplot(2,4,3);
JMc=JMa+JMb;
imshow(JMc)
%d.对原图像试用Sobel梯度操作
subplot(2,4,4);
hx=[-1,-2,-1;0,0,0;1,2,1];%生产sobel垂直梯度模板
hy=[-1,0,1;-2,0,2;-1,0,1];%生产sobel水平梯度模板
gradx=filter2(hx,JMa,'same');
gradx=abs(gradx);%计算图像的sobel垂直梯度
grady=filter2(hy,JMa,'same');
grady=abs(grady);%计算图像的sobel水平梯度
JMd=gradx+grady;%得到图像的sobel梯度
imshow(JMd,[]);
title('d:
图sobel梯度处理后图像');
%e.使用大小为5*5的一个均值滤波器得到平滑后的Sobel梯度图像。
subplot(2,4,5);
h1=fspecial('average',5);
JMe=imfilter(JMd,h1);
imshow(JMe);
title('e:
使用5*5均值滤波器平滑后的sobel图像');
%f.将拉普拉斯图像(c)与平滑后的梯度图像(e)进行点乘。
subplot(2,4,6);
JMf=(JMc).*JMe;
imshow(JMf);
title('f:
图c和图e相乘相乘的掩蔽图像');
%g,将乘积图像(即图f)与原图像相加。
subplot(2,4,7);
JMg=JMa+JMf;
imshow(JMg);
title('乘积图像(即图f)与原图像相加得到的锐化图像');
%h.对图7进行幂率变换处理,r=0.5,c=1
subplot(2,4,8);
gamma=0.5;
c=1;
JMh=c.*JMg.^gamma;
imshow(JMh);
title('h:
图g进行幂率变换后的图像');
具体图形变换如图