河北工业大学《计算机图像处理》实验指导书Word下载.docx
- 文档编号:20642469
- 上传时间:2023-01-24
- 格式:DOCX
- 页数:20
- 大小:939.95KB
河北工业大学《计算机图像处理》实验指导书Word下载.docx
《河北工业大学《计算机图像处理》实验指导书Word下载.docx》由会员分享,可在线阅读,更多相关《河北工业大学《计算机图像处理》实验指导书Word下载.docx(20页珍藏版)》请在冰豆网上搜索。
clear%执行本命令将会清除内存中的全部变量
figure
(1);
%生成一个图像窗口1
I=imread('
e:
'
);
%将硬盘e:
根目录上的图像文件的数据读入矩阵变量I中
—
imshow(I);
%在当前的图像窗口中显示图像矩阵I
title('
原始图像'
%在当前的图像窗口中加上标题
但为了能够对程序进行调试和重复应用,我们要求用M文件的方式完成实验中各个程序的编写。
2、数字图像的表示和类别
根据图像数据矩阵解释方法的不同,MATLAB把其处理为四类:
(1)亮度图像(Intensityimages)
一幅亮度图像是一个数据矩阵,其归一化的取值表示亮度。
若亮度图像的像素都是uint8类或uint16类,则它们的整数值范围分别是[0,255]和[0,65536]。
若图像是double类,则像素取值就是浮点数。
规定双精度型归一化亮度图像的取值范围是[0,1]。
(2)二值图像(Binaryimages)
一幅二值图像是一个取值只有0和1的逻辑数组。
而一幅取值只包含0和1的uint8类数组,在MATLAB中并不认为是二值图像。
使用logical函数可以把数值数组转化为逻辑数组,其语法为B=logical(A);
其中A是由0和1构成的数值数组。
要测试一个数组是否为逻辑数组,可以使用函数islogical(C);
若C是逻辑数组,则该函数返回1;
否则,返回0。
(3)索引图像(Indexedimages)
索引颜色通常也称为映射颜色,在这种模式下,颜色都是预先定义的,并且可供选用的一组颜色也很有限,索引颜色的图像最多只能显示256种颜色。
一幅索引颜色图像在图像文件里定义,当打开该文件时,构成该图像具体颜色的索引值就被读入程序里,然后根据索引值找到最终的颜色。
、
(4)RGB图像(RGBimages)
一幅RGB图像就是彩色像素的一个M×
N×
3数组,其中每一个彩色像素点都是在特定空间位置上相对应的红、绿、蓝三个分量。
按照惯例,形成一幅RGB彩色图像的三个矩阵常称为红、绿或蓝分量图像。
令fR,fG和fB分别代表三种RGB分量图像。
一幅RGB图像就利用cat(级联)操作将这些分量图像组合成彩色图像:
rgb_image=cat(3,fR,fG,fB)在操作中,图像按顺序放置。
3、数据类和图像类型间的转化
表1-2中列出了MATLAB和ImageProcessingToolbox为表示图像所支持的各种数据类。
表中的前8项称为数值数据类,第9项称为字符类,最后一项称为逻辑数据类。
表1-3列出了工具箱中提供的一些必要函数,以便在图像类和数据类之间进行转化。
表1-2MATLAB和IPT支持数据类型
名称
double
双精度浮点数,范围为
uint8
无符号8比特整数,范围为[0255]
"
uint16
无符号16比特整数,范围为[065536]
uint32
无符号32比特整数,范围为[095]
int8
有符号8比特整数,范围为[-128127]
int16
有符号16比特整数,范围为[-3276832767]
int32
有符号32比特整数,范围为[-4847]
single
单精度浮点数,范围为
char
字符
logical
值为0或1
表1-3格式转换函数
将输入转化为
有效的输入图像数据类
im2uint8
logical,uint8,uint16和doulbe
im2uint16
mat2gray
double,范围为[01]
im2double
im2bw
uint8,uint16和double
4、常用函数和命令
(1)读写图像文件
a)imread:
用于读入各种图像文件,返回一个矩阵,如:
a=imread('
d:
\'
)
b)whos:
用于读取图像的基本信息,如:
whosa
c)imwrite:
用于写入图像文件,如:
imwrite(a,'
’tif’)
d)imfinfo:
用于读取图像文件的有关信息,如:
imfinfo('
(2)图像的显示
a)image:
image函数是MATLAB提供的最原始的图像显示函数,如:
a=[1,2,3,4;
4,5,6,7;
8,9,10,11];
image(a);
b)
【
c)imshow:
用于图像文件的显示,如:
i=imread('
imshow(i);
d)colorbar:
用于显示图像的颜色条,如:
colorbar;
e)figure:
用于设定图像显示窗口,如:
figure
(1);
figure
(2);
%n为图形窗口号数
f)subplot:
把图形窗口分成多个矩形部分,每个部分可以分别用来进行显示,如:
subplot(m,n,p)
分成m×
n个小窗口,在第p个窗口中创建坐标轴为当前坐标轴,用于显示图形
g)plot:
绘制二维图形,如:
plot(y);
plot(x,y);
%x、y可以是向量、矩阵
(3)图像类型转换
a)rgb2gray:
把真彩图像转换为灰度图像,如:
i=rgb2gray(j)
b)im2bw:
通过阈值化方法把图像转换为二值图像,如:
I=im2bw(j,level)
Level表示灰度阈值,取值范围0~1,即表示阈值取自原图像灰度范围的n%
%
(4)其它运算函数
a)zeros:
生成全0数组或矩阵,如:
B=zeros(m,n)或B=zeros([mn]),返回一个m×
n的全0矩阵
b)取整函数:
round四舍五入取整函数;
floor最小取整函数;
ceil最大取整函数,
如:
a=[+];
I=round(a)
输出I=[-20362+4i]
三、实验内容
1、读入一幅RGB图像(),分别转换为灰度图像和二值图像,并在同一个窗口内分成三个子窗口来分别显示这三幅图像,注上文字标题,即“原图像”、“灰度图像”、“二值图像”。
subplot(1,3,1);
I=imread('
imshow(I);
!
title('
原图像'
subplot(1,3,2);
J=rgb2gray(I);
imshow(J);
灰度图像'
subplot(1,3,3);
Q=im2bw(J,;
imshow(Q);
二值图像'
¥
2、说明以下程序的显示结果为何是一幅几乎全黑的图像。
clear;
closeall;
myi=zeros(20,20);
myi(2:
2:
18,2:
18)=1;
myi=uint8(myi);
figure,imshow(myi,'
notruesize'
3、阅读分析下列程序,对每条语句给出注释,运行并解释最终执行结果;
若将for循环改为方框中的程序,有何不同
clearall;
closeall;
[height,width]=size(I);
figure;
;
subplot(2,3,1);
L=1;
form=1:
5
L=2*L;
quartimage=I(1:
L:
height,1:
width);
subplot(2,3,m+1);
imshow(uint8(quartimage));
end
?
quartimage=zeros(ceil(height/L),ceil(width/L));
k=1;
n=1;
fori=1:
height
forj=1:
width
quartimage(k,n)=I(i,j);
n=n+1;
end
k=k+1;
end
*
4、编写程序将一幅灰度图像()每两行之间增加一行,即将图像的高度拉伸一倍、宽度不变,然后将其存储为图像文件。
要求拉伸后图像的奇数行对应原图像的各行,偶数行是原图像相邻两行的灰度平均值(即偶数行由插值得到),并将结果与MATLAB图像处理工具箱中提供的改变图像大小的函数imresize(imresize(I,[mn])将图像I大小调整为m行n列)的处理结果进行比较。
T1=maketform('
affine'
[100;
020;
001]);
c:
I1=imtransform(I,T1);
imshow(I),figure,imshow(I1)
5、熟悉数字图像处理常用函数的使用,调出帮助文档查看其各种不同用法。
MatLab具有完善的帮助系统,包括命令行帮助、联机帮助和演示帮助等,应学会充分利用帮助系统来解决问题。
获得帮助的途径有多种,如:
(1)在MatLab界面中单击工具条上的问号,或单击Help菜单中的MATLABHelp选项;
或按F1键打开Help窗口。
(2)选择函数(函数所在区变暗),点右键弹出菜单,选择HelponSelection;
或在command窗口键入help函数名(或工具箱名)。
实验二图像直方图及灰度变换
1.掌握图像灰度直方图的概念及其计算方法,编写灰度直方图统计程序。
2.通过对图像直方图的分析,学习应用直方图法解决诸如图像二值化等具体问题。
!
3.熟悉直方图均衡化的计算过程及其应用。
4.掌握图像灰度变换技术,通过调整图像的对比度和亮度等参数,改善视觉效果。
直方图是图像最基本的统计特征,是图像亮度分布的概率密度函数,反映了图像灰度值的分布情况。
直方图是多种空间域处理技术的基础。
直方图操作能有效地用于图像增强,如通过直方图均衡化处理,可使图像在整个灰度级范围内的分布均匀化,即在每个灰度级上都具有相同的像素点数,从而获得较好的视觉效果。
另外,直方图固有的信息也可用在图像分割等其它图像处理的应用中。
灰度级变换技术可用g(x,y)=T[f(x,y)]的形式表示,其中f(x,y)为输入图像,g(x,y)为输出图像,T是对图像f进行某种处理的操作。
由于(x,y)处的g值仅由f在该点处的亮度决定,T也称为一个亮度或灰度级变换函数,它与位置坐标(x,y)无关,所以通常写成如下的简化形式s=T(r),其中r和s分别表示图像f和g在相应点(x,y)的亮度。
利用灰度变换可调整一幅图像的明暗、对比度等。
1、编写一个图像灰度直方图统计函数my_imhist,选择一幅图像利用my_imhist显示其直方图,将结果与MATLAB图像处理工具箱中提供的灰度直方图函数imhist的处理结果进行比较,并在同一窗口中显示出来。
C:
J=a;
L=256;
my=zeros(L,1);
[row,col]=size(a);
n=row*col;
for
i
=
1:
row
j
col
num
double(a(i,j))+1;
my(num)
my(num)+1;
%计算直方图概率估计
i=1:
L
subplot(1,2,2);
imhist(a);
imhist处理'
subplot(1,2,1);
bar(my);
自编函数处理'
2、利用以上编写的函数my_imhist或imhist,估算图像中瞳孔的半径(以像素为单位)。
100~200之间的单位个数为8,则瞳孔面积大概为80000像素,利用求半径公式S=^2得到半径R=√ ̄(S/=√ ̄(80000/≈160像素
img=imread('
imshow(img);
[M,N]=size(img);
img=double(img);
h=zeros(256,256);
M,
j=1:
N,
f=img(i,j);
h(f+1)=h(f+1)+1;
end
bar(h);
grid;
I
imread('
D:
imhist(I)
3、利用MATLAB提供的图像直方图修正函数histeq,选择一幅直方图不均匀的图像如,对其进行直方图均衡化处理,显示处理前后的图像以及它们的灰度直方图,体会直方图均衡化算法的特点。
H=histeq(I);
subplot(2,2,1),imshow(I);
源图像'
subplot(2,2,2),imshow(H);
均衡化后的图像'
subplot(2,2,3),imhist(I);
源图像的直方图'
subplot(2,2,4),imhist(H);
均衡化后的图像直方图'
4、按照教材68页上的公式(4.1.6)
(4.1.6)
以及公式
,分别编程实现图像的分段线性灰度变换和非线性灰度变换。
利用自己编写的代码对图像进行灰度调整,使拍摄主体瓶子的花纹更加清晰地表现出来,以改善视觉效果。
将同一个图形窗口分成四个子窗口来分别显示灰度变换前、后的图像以及它们的灰度直方图,注上文字标题。
并考虑利用MATLAB提供的imadjust函数是否能完成同样的工作。
~
F=imread('
[row,col]=size(F);
max=double(F(1,1))+1;
fori=1:
row
k=double(F(i,j))+1;
ifmax<
k
max=k;
disp('
灰度最大值:
max
--------------首先进行线性变换--------------'
fo=input('
请输入源图像的灰度范围:
a=fo
(1);
b=fo
(2);
ifb>
disp('
输入源图像的范围不能大于最大值!
re=input('
请输入变换后的灰度范围:
c=re
(1);
d=re
(2);
ifd>
256%判断d是否超过256
输入变换后的范围不能超过255!
G=zeros(row,col);
ifk<
a
G(i,j)=(c/a)*F(i,j);
elseifk<
b
G(i,j)=((d-c)/(b-a))*(F(i,j)-a)+c;
else
G(i,j)=((255-d)/(max-b))*(F(i,j)-b)+d;
out=uint8(round(G-1));
figure,imshow(F);
figure,imshow(out);
线性变换后'
figure,imhist(F);
figure,imhist(out);
线性变换后的直方图'
--------下面进行非线性变换,其变换公式为g=c*f^y--------'
y=input('
请输入指数y:
c=input('
请输入乘数c:
g=zeros(row,col);
k=double(F(i,j))+1;
g(i,j)=c*(k^y);
max=g(1,1);
g(i,j)
max=g(i,j);
scale=max/256;
g(i,j)=round(g(i,j)/scale);
out2=uint8(round(g-1));
figure,imshow(out2);
非线性变换后'
figure,imhist(out2);
非线性变换后的直方图'
实验三图像增强与边缘检测
1.加深对图像增强及边缘检测技术的感性认识,应用MATLAB工具箱自带的处理函数或自己编程完成相关的工作,分析处理结果,巩固所学理论知识。
2.熟练掌握空域滤波中常用的平滑和锐化滤波器,针对不同类型和强度的噪声,进行滤波处理,体会并正确评价滤波效果,了解不同滤波方式的使用场合,能够从理论上作出合理的解释。
3.掌握基于边缘检测的图像分割方法,通过实验体会主要的边缘检测算子在处理结果上的异同,探讨不同因素对轮廓提取效果的影响,培养处理实际图像的能力。
图像增强是指按特定的需要突出一幅图像中的某些有用信息,同时消弱或去除某些不需要的信息的处理方法,其主要目的是使处理后的图像对某些特定的应用比原来的图像更加有效。
图像平滑与锐化处理是图像增强的主要研究内容。
图像分割是由图像处理过渡到图像分析的关键步骤,一般是按照一定的规则,把图像分成互不重叠的若干区域或子集,并提取出感兴趣目标的技术和过程。
本实验主要研究利用边缘提取算子的分割方法。
和本实验有关的几个常用Matlab函数:
·
(1)imnoise:
用于对图像生成模拟噪声,如:
j=imnoise(i,'
gaussian'
0,%在图像i上叠加均值为0、方差为的高斯噪声,得到含噪图像j
salt&
pepper'
%在图像i上叠加密度为的椒盐噪声,得到含噪图像j
(2)fspecial:
用于产生预定义滤波器,如:
h=fspecial('
average'
3);
%产生3×
3模板的均值滤波器
sobel'
%产生sobel水平边缘增强的滤波器
可选项还有:
高斯低通滤波器、'
laplacian'
拉普拉斯滤波器、'
log'
高斯拉普拉斯滤波器等
(3)imfilter、filter2、conv2:
均是基于卷积的图像滤波函数,都可用于图像滤波,用法类似,如:
prewitt'
%产生prewitt算子的水平方向模板
j1=imfilter(i,h);
%或者j2=filter2(h,i);
或者j3=conv2(i,h);
。
(4)medfilt2:
用于图像的中值滤波,如:
j=medfilt2(i,[MN]);
%对矩阵i进行二维中值滤波,邻域为M×
N,缺省值为3×
3
(5)edge:
检测灰度或二值图像的边缘,返回一个二值图像,1像素是检测到的边缘,0像素是非边缘
用法:
BW=edge(I,'
thresh,direction);
%I为检测对象;
边缘检测算子可用sobel、roberts、prewitt、zerocross、log、canny;
thresh指定阈值,检测时忽略所有小于阈值的边缘,默认自动选择阈值;
direction指定方向,可选项有horizontal、vertical或both,在指定的方向上用算子进行边缘检测
1、图像平滑(去噪):
选择一幅图像,例如,分别叠加零均值高斯噪声、椒盐噪声,然后分别采用3×
3、5×
5不同大小的窗口,再分别利用邻域平均法、中值滤波法对噪声图像进行滤波,最后采用便于观察比较的形式,显示出滤波前后的各个图像,分析不同滤波器对不同噪声的处理效果及其优缺点。
思考如何寻找合适的滤波器,使之能较好地滤除噪声,同时又能保持细节和边缘清晰。
J0=imnoise(I,'
J1=imnoise(I,'
J0J3=imfilter(J0,fspecial('
));
J0J5=imfilter(J0,fspecial('
[55]
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算机图像处理 河北 工业大学 计算机 图像 处理 实验 指导书