MATLAB数字图像处理(二)图像增强
1 图像增强
1.1 直方图均衡化
对于灰度图像,可以使用直方图均衡化的方法使得原图像的灰度直方图修正为均匀的直方图。
代码如下:
I2=histeq(I1);
figure,imshow(I2);
figure,imhist(I2);
原图像为lena的图片,经过直方图均衡化后的图片为Figure 1。
Figure 1
原图灰度直方图请参见:MATLAB数字图像处理(一)基础操作和傅立叶变换,经过修正后的灰度直方图为Figure 2:
Figure 2
可以看出,与原图对比,灰度分布更加均匀。
1.2 灰度变换
在MATLAB中,可以使用imadjust函数对图像进行线性变换。常用的为线性变换。使用的函数为imadjust(I, [low_in, high_in], [low_out high_out],gamma),其中,in为输入的图片变化的灰度范围,out为输出的图片变化的灰度范围,均为[0 1]之间。也就是将low_in和high_in之间的值映射到low_out和high_out之间。gamma为线性值,若为1则为线性变换。
代码和效果如下:
I3=imadjust(I1,[0.3 0.7],[0.1 0.9],1);
figure,imshow(I3);
figure,imhist(I3);
Figure 3
Figure 4
其中,突出了原图0.3-0.7之间的灰度元素。
1.3 空间域平滑
图像平滑是为了抑制图像噪声,改善图像质量而进行的图片处理。在这里,笔者将首先介绍MATLAB中自带的均值滤波和中值滤波实现,之后设计自己的模版算子, 再次分别实现了均值滤波和中值滤波。
MATLAB中自带的均值滤波可以采用fspecial(‘params’,n)函数实现,其中,params代表的是滤波的种类,常用的有average(均值滤波),还有可以在空间域锐化中用到的高斯、拉普拉斯、prewitt、sobel算子等。n代表的是模版的在这里仅仅介绍均值滤波。
MATLAB中,均值滤波代码如下:
avgModel=fspecial('average',3);
Iavg=filter2(avgModel,I2)/255;
figure,imshow(Iavg);
由于均值滤波对于高斯噪声效果较好,故对高斯噪声的图(参见MATLAB数字图像处理(一)基础操作和傅立叶变换)进行了图像平滑。效果如下(Figure 5):
Figure 5
而笔者根据均值滤波原理,自己写了一段代码,如下:
n=3;
model(1:n,1:n)=1;
Iavg=I2;
for i=2:length(Iavg)-1
for j=2:length(Iavg(2,:))-1
x=I2(i-(n-1)/2:i+(n-1)/2,j-(n-1)/2:j+(n-1)/2);
xuint=uint8(model).*uint8(x);
x_avg=mean(xuint(:));
Iavg(i,j)=x_avg;
end
end
figure,imshow(Iavg);
效果(Figure 6):
Figure 6
可以发现,除了运行速度稍慢,效果与MATLAB自带的函数相同。
MATLAB中自带的中值滤波可以采用medfilt2(I,[n n])函数实现,其中,I代表的是需要处理的图片,n代表的是模版的大小。
MATLAB中,中值滤波代码如下:
Imid=medfilt2(I3,[3,3]);
figure,imshow(Imid);
由于中值滤波对于椒盐噪声效果较好,故对椒盐噪声的图(参见MATLAB数字图像处理(一)基础操作和傅立叶变换)进行了图像平滑。效果如下(Figure 7):
Figure 7
而笔者根据中值滤波原理,自己写了一段代码,如下:
n=3;
Imid=I3;
for i=2:length(I1)-1
for j=2:length(I1(2,:))-1
x=I3(i-(n-1)/2:i+(n-1)/2,j-(n-1)/2:j+(n-1)/2);%3*3
x_mid=median(x(:));
Imid(i,j)=x_mid;
end
end
figure,imshow(Imid);
效果(Figure 8):
Figure 8
可以发现,除了运行速度稍慢,效果与MATLAB自带的函数相同。
1.4 空间域锐化
图像锐化有多种方法,最常用的锐化算子包括梯度算子、Roberts算子、Prewitt算子、Sobel算子以计算梯度的梯度锐化法,Laplacian算子等。接下来,笔者将一一介绍相关的算子。
笔者根据梯度算子的计算原理,编写了如下代码:
Igrad=I1;
for i=1:length(Igrad)-1
for j=1:length(Igrad(2,:))-1
x=I1(i,j+1)-I1(i,j);
y=I1(i+1,j)-I1(i,j);
grad=max(abs(x),abs(y));
Igrad(i,j)=grad;%¸³Öµ
end
end
figure,imshow(Igrad);
效果(Figure 9):
Figure 9
MATLAB中自带了Prewitt算子和Sobel算子的函数,代码及效果如下:
model=fspecial('prewitt');
Iprewitt=filter2(model,I1);
figure,imshow(Iprewitt);
Figure 10
model=fspecial('sobel');
Isobel=filter2(model,I1);
figure,imshow(Isobel);
Figure 11
事实上,可以发现这几种算子之间不同的地方只是在于模版,因此笔者编写了一个可以自定义各种模版的代码,如下:
n=2
modelx=[-1 0;0 1];
modely=[0 -1;1 0];
Iend=I1;
Idouble=double(I1);
for i=1:length(Idouble)-5+n
for j=1:length(Idouble(2,:))-5+n
area=Idouble(i:i+n-1,j:j+n-1);
x=area.*modelx;
y=area.*modely;
grad=max(abs(sum(x(:))),abs((sum(y(:)))));
Iend(i,j)=grad;%¸³Öµ
end
end
figure,imshow(Iend);
可以显示出Roberts算子的效果:
Figure 12
与Prewitt算子、Sobel算子经过对比后发现效果相同,说明有效。
对于Laplacian算子,同样可以使用MATLAB中自带的函数,也可以进行自定义。自带函数使用方法如下:
model=fspecial('laplacian');
Ilaplacian=filter2(model,I1);
figure,imshow(Ilaplacian,[]);
效果:
Figure 13
套用前面自定义模版后的代码如下:
n=3;
model=[0 -1 0;-1 5 -1;0 -1 0];
Iend=I1;
Idouble=double(I1);
for i=1:length(Idouble)-5+n
for j=1:length(Idouble(2,:))-5+n
area=Idouble(i:i+n-1,j:j+n-1);
x=area.*model;
grad=sum(x(:));
Iend(i,j)=grad;
end
end
figure,imshow(Iend);
效果:
Figure 14
与教科书上的效果相同。
1.5 频率域图像增强
频率域增强包括高通滤波器和低通滤波器,主要是使用傅立叶变换和逆变换对图像进行处理。在MATLAB中没有直接的滤波器函数,但是可以结合原理与MATLAB中的数学公式,构造各种滤波器。在本文中,笔者分别构造了理想、Butterworth、指数的高通和低通滤波器,并一一进行了实验。
首先需要提前设定截止频率D0和阶数n。之后就可以套用具体的公式进行计算了。部分代码如下:
n=6;%滤波器的阶数
D0=80;%滤波器的截止频率
[u,v]=meshgrid(-M/2:(M/2-1),-N/2:(N/2-1));%产生离散数据
D=sqrt(u.^2+v.^2);
%各种不同的滤波器
H1=double(D<=D0);%理想低通滤波器
H2=1./(1+(D/D0).^(2*n));%Butterworth低通滤波器
H3=exp(-(D/D0).^n);%指数低通滤波器
H4=double(D>=D0);%理想高通滤波器
H5=1./(1+(D0./D).^(2*n));%Butterworth高通滤波器
H6=exp(-(D0./D).^n);%指数高通滤波器
之后,可以使用各种不同的滤波器对图片进行处理。效果如下:
Figure 15
Figure 16
Figure 17
Figure 18
Figure 19
Figure 20
通过以上一系列的图片可以看出,低通滤波器可以抑制图像噪声,改善图像质量,高通滤波器可以突出图像的边界。
具体代码请参考:https://github.com/kkyyhh96/DigitalImageProcessing/tree/master/code中kyh_Three.m和kyh_Four.m的内容。
MATLAB数字图像处理(二)图像增强的更多相关文章
- matlab数字图像处理 入门基础
本代码基于张铮主编的<数字图像处理与机器视觉>一书. 源图片:lena A = imread ('lena.jpg'); %读入图像lena.jpg,赋给变量A %imwrite(A,'l ...
- MATLAB数字图像处理(一)基础操作和傅立叶变换
数字图像处理是一门集计算机科学.光学.数学.物理学等多学科的综合科学.随着计算机科学的发展,数字图像处理技术取得了巨大的进展,呈现出强大的生命力,已经在多种领域取得了大量的应用,推动了社会的发展.其中 ...
- 第一章 MATLAB数字图像处理编程基础
1 为什么用MATLAB MATLAB的图像处理工具箱(Image Processing Toolbox,IPT)封装了一系列不同图像处理需求的标准算法,它们都是通过直接或间接调用MATLAB中矩阵运 ...
- MATLAB数字图像处理学习笔记
我们都知道一幅图片就相当于一个二维数组,可以用一个矩阵来表示,而MATLAB可以说就是为矩阵运算而生的,所以学习图像处理,学习MATLAB势在必行! 一. MATLAB基础知识 1. 读取图像 %im ...
- MATLAB数字图像处理基础
图像的输入.输出和显示 1.图像的输入 imread('filename'), 实际中写的是 >> f = imread('sky.jpg'); 2.图像的显示 imshow ...
- 数字图像处理实验(5):PROJECT 04-01 [Multiple Uses],Two-Dimensional Fast Fourier Transform 标签: 图像处理MATLAB数字图像处理
实验要求: Objective: To further understand the well-known algorithm Fast Fourier Transform (FFT) and ver ...
- MatLab数字图像处理实战(赵小川)-sift原理
- matlab数字图像处理-冈萨雷斯-数据类和图像类之间的转换
亮度图像 二值图像 属于注释 数据类间的转换 图像类和类型间的转化 把一个double类的任意数组转换成[0,1]的归一化double类数组----->mat2gray 图像类和类型间的转化例题 ...
- matlab数字图像处理-冈萨雷斯-读取,显示,保存图像
图像读取:imread(filename) 显示图像‘ 显示多幅图像 保存图像 b 计算压缩比
随机推荐
- 如何卸载Python通过setup.py安装的模块
0,pip uninstall xxxx 1.找到egg sudo easy_install -m BitVector .... Using /usr/local/lib/python2./ ...
- Restoring Road Network(Floyd算法的推广)
个人心得:看懂题目花费了不少时间,后面实现确实时间有点仓促了,只是简单的做出了判断是否为真假的情况, 后面看了题解发现其实在判断时候其实能够一起解决的,算了,基础比较差还是慢慢的来吧. 题意概述: 就 ...
- Milking Time(DP)
个人心得:一开始自己找状态,是这么理解的,只要前面一个满足就等于此时的值加上d(n-1),否则就是不挖此时的比较d(n-1)和 d(n-2)+cost,不过仔细一想忽略了很多问题,你无法确定n-2和此 ...
- 5、Selenium+Python自动登录163邮箱发送邮件
1.Selenium实现自动化,需要定位元素,以下查看163邮箱的登录元素 (1)登录(定位到登录框,登录框是一个iframe,如果没有定位到iframe,是无法定位到账号框与密码框) 定位到邮箱框( ...
- 安装成功的nginx如何添加未编译模块?
在重启nginx后发生了错误,错误如下: nginx: [emerg] the "ssl" parameter requires ngx_http_ssl_module in /u ...
- Java基础--枚举Enum
Java中的枚举是一种特殊的类,可以将一组固定常量的集合组成一种类型,使用方便且类型安全.使用enum关键字定义. enum类型父类为Enum,通过Enum.class可见Enum为抽象类,实现了Co ...
- python风格指南
还在让别人一眼就看出你是一只野生程序猿嘛,快来看看谷歌的python风格指南提升逼格吧! http://zh-google-styleguide.readthedocs.io/en/latest/go ...
- Mongodb 副本集的节点详细操作
副本集操作 官方文档:https://docs.mongodb.com/v3.2/reference/method/js-replication/ 1 rs.add(){ _id: <int&g ...
- 清除SUN服务器部件的suspect状态
对于suspect状态的部件,可以用setchs命令清除其状态.如果ScApp的版本在5.20.15之前,需要进入service模式后才能使用setchs命令.如果ScApp版本 升级到5.20.15 ...
- 11-23网页基础--JavaScript基础知识
第一课 JavaScript简介 一.定义:JavaScript是脚本语言,需要宿主文件,它的宿主文件是html文件. JavaScript 是一种轻量级的编程语言. JavaScript 是可插入 ...