MATLAB选修论文二.docx

MATLAB选修论文二.docx

《MATLAB选修论文二.docx》由会员分享，可在线阅读，更多相关《MATLAB选修论文二.docx（13页珍藏版）》请在冰豆网上搜索。

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进行幂率变换后的图像'）;

具体图形变换如图