整理实验四形态学Word格式文档下载.docx
《整理实验四形态学Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《整理实验四形态学Word格式文档下载.docx(16页珍藏版)》请在冰豆网上搜索。
BW2=imdilate(BW,se);
figure,imshow(BW2,'
B)改变上述结构元素类型(如:
line,diamond,disk等),重新进行膨胀操作。
Line:
BW=zeros(9,10);
line'
3,3);
C)对图像‘text.png’进行上述操作,观察不同结构元素膨胀的效果。
BW3=imread('
text.png'
);
imshow(BW3)
se2=strel('
11,90);
%线型结构元素
BW4=imdilate(BW3,se2);
figure,imshow(BW4)
2.图像腐蚀
A)对图像‘circbw.tif’进行腐蚀操作。
BW1=imread('
circbw.tif'
arbitrary'
eye(5));
BW2=imerode(BW1,se);
imshow(BW1)
figure,imshow(BW2)
B)对图像‘text.tif’进行腐蚀操作。
BW=imread('
text.tif'
BW2=imerode(BW3,se);
imshow(BW)
3.膨胀与腐蚀的综合使用
A)从原始图像‘circbw.tif’中删除电流线,仅保留芯片对象。
方法一:
先腐蚀(imerode),再膨胀(imdilate);
rectangle'
[4030]);
%选择适当大小的矩形结构元素
%先腐蚀,删除较细的直线
BW3=imdilate(BW2,se);
%再膨胀,恢复矩形的大小
figure,imshow(BW3)
方法二:
使用形态开启函数(imopen)。
[30,20]);
BW2=imopen(BW1,se);
%开启操作
B)改变结构元素的大小,重新进行开启操作,观察处理结果。
se=strel(‘rectangle’,[2010]);
se=strel(‘rectangle’,[5040]);
C)置结构元素大小为[43],同时观察形态开启(imopen)与闭合(imclose)的效果,总结形态开启与闭合在图像处理中的作用。
I=imread('
imshow(I)
[43]);
I1=imopen(I,se);
I2=imclose(I,se);
%闭合操作
figure,imshow(I1)
figure,imshow(I2)
4.*高帽与低帽变换
A)读入图像‘pearlite.tif’,分别显示其高帽变换与低帽变换结果,并与原图像比较。
(设se=strel('
disk'
5);
)。
pearlite.tif'
subplot(221),imshow(I)
J=imtophat(I,se);
subplot(222),imshow(J)
K=imbothat(I,se);
subplot(223),imshow(K)
L=imsubtract(imadd(J,I),K);
subplot(224),imshow(L)
B)要求显示在一个窗口中。
5.图像极值的处理方法
A)对于下图所示的图像矩阵A,利用函数imregionalmax寻找其局部极大值
A=[10101010101010101010;
10131313101011101110;
10131313101010111010;
10101010101010101010;
10111010101818181010;
10101011101818181010;
10101110101818181010;
10111011101010101010;
10101010101011101010];
B=imregionalmax(A)
B=
0000000000
0111001010
0111000100
0100011100
0001011100
0010011100
0101000000
0000001000
B)利用函数imextendedmax寻找像素值大于其邻域像素值2个单位的局部极大值。
C=imextendedmax(A,2)
C=
0111000000
0000011100
6.*创建一幅图像,求其距离矩阵。
bw=zeros(5,5);
bw(2,2)=1;
bw(4,4)=1;
D=bwdist(bw)
center1=-10;
center2=-center1;
dist=sqrt(2*(2*center1)^2);
radius=dist/2*1.4;
lims=[floor(center1-1.2*radius)ceil(center2+1.2*radius)];
[x,y]=meshgrid(lims
(1):
lims
(2));
bw1=sqrt((x-center1).^2+(y-center2).^2)<
=radius;
bw2=sqrt((x-center2).^2+(y-center2).^2)<
bw1=sqrt((x-center1).^2+(y-center1).^2)<
bw=bw1|bw2;
imshow(bw)
D=bwdist(bw);
figure,imshow(D,[])
D1=bwdist(~bw);
figure,imshow(D1,[])
D=
1.41421.00001.41422.23613.1623
1.000001.00002.00002.2361
1.41421.00001.41421.00001.4142
2.23612.00001.000001.0000
3.16232.23611.41421.00001.4142
7.*使用伪彩色显示标记矩阵。
BW=[00000000;
01100111;
01100011;
01100000;
00011000;
00000000]
X=bwlabel(BW,4)
RGB=label2rgb(X,@jet,'
k'
imshow(RGB,'
8.利用选择控制点实现图像匹配。
Matlab图像匹配的步骤:
将标准图像和待匹配图像读入Matlab;
指定图像中的控制点对并保存;
使用互相关性进一步协调控制点对(可选);
制定所需变换类型并根据控制点对推断变换参数;
变换未匹配的图像。
%读入待匹配图像和标准图像
(五)建设项目环境影响评价文件的审批unregistered=imread('
westconcordaerial.png'
imshow(unregistered)
1.依法评价原则;
orthophoto=imread('
westconcordorthophoto.png'
报告内容有:
建设项目基本情况、建设项目所在地自然环境社会环境简况、环境质量状况、主要环境保护目标、评价适用标准、工程内容及规模、与本项目有关的原有污染情况及主要环境问题、建设项目工程分析、项目主要污染物产生及预计排放情况、环境影响分析、建设项目拟采取的防治措施及预期治理效果、结论与建议等。
figure,imshow(orthophoto)
在可行性研究时应进行安全预评价的建设项目有:
%选择图像中对应的控制点
cpselect(unregistered(:
:
1),orthophoto)
每名环境影响评价工程师申请登记的类别不得超过2个。
%保存控制点对
(6)列出选定的评价方法,并作简单介绍。
input_points
(二)规划环境影响评价的技术依据和基本内容base_points
%使用相关性进一步协调控制点对
input_points_corr=cpcorr(input_points,base_points,unregistered(:
1),orthophoto);
input_points_corr
(2)是否符合国家产业政策和清洁生产标准或要求。
%根据控制点对推断空间变换参数
6.提出安全对策措施建议mytform=cp2tform(input_points,base_points,'
projective'
%变换未匹配的图像
registered=imtransform(unregistered,mytform);
(2)防护支出法figure,imshow(registered)
要求做前5个,6-8为选择实验。