形态学处理代码matlab,（8）MATLAB形态学处理

形态学处理代码matlab,（8）MATLAB形态学处理

腐蚀，膨胀，开运算：先腐蚀后膨胀 ；

腐蚀

I = imread('../erode_'); %读⼊8位灰度图像

%⼆值形态学处理中将灰度图像中所有⾮0值都看作是1，即前景物体

figure, imshow(I);

= strel('square', 3) %3*3的正⽅形结构元素

Ib= imerode(I, ); %腐蚀

figure, imshow(Ib);

I = imread('../'); %读⼊图像

figure, imsho世界第二高峰是什么 w(I);

Ie1 = imerode(I, [1 1 1; 1 1 1; 1 1 1]) ; %3*3正⽅形结构元素的腐蚀

figure, imshow(Ie1); %得到图10.10(b)

Ie2 = imerode(Ie1, [0 1 0; 1 1 1; 0 1 0]); %3*3⼗字形结构元素的腐蚀

figure, imshow(Ie2);

获取纹理⽅法

I = imread('../head_'); %读⼊原图像

figure, imshow(I); %得到图10.18(a)中的图像

= strel('square', 3); %3*3的正⽅形结构元素

Ie = imerode(I, ); %腐蚀得到内部点

Iout = I - Ie; %减去内部点留下边界点

figure, imshow(Iout); %得到图10.18(b)中的图像

% locateMouth.m

I = imread('../'); %读⼊图像

Id = im2double(I);

figure, imshow(Id) % 得到10.24(a)

Ibw = im2bw(Id, 0.38); %以0.38为阈值⼆值化

Ibw = 1 - Ibw; %为在Matlab中进⾏处理，将图像反⾊

figure, imshow(Ibw) % 得到10.24(b)

hold on

[L, num] = bwlabel(Ibw, 8); % 标注连通分量

disp(['图中共有' num2str(num) '个连通分量'])

% 找出最⼤的连通分量(嘴)

max = 0; % 当前最⼤连通分量的⼤⼩

indMax = 0; % 当前最⼤连通分量的索引

for k = 1:num

[y x] = find(L == k); % 找出编号为k的连通区的⾏索引集合y和列索引集合x

nSize = length(y); % 计算该连通区的像素数⽬

if(nSize > max)

max 涯组词 = nSize;

indMax = k;

end

end

if indMax == 0

disp('没有找到连通分量')

return

end

% 计算并显⽰最⼤连通分量(嘴)的中⼼

[y x] = find(L == indMax);

yMean = mean(y);

xMe樱桃造句 an = mean(x);

plot(xMean, yMean, 'Marker', 'o', 'MarkerSize', 14, 'MarkerEdgeColor', 'w', 'MarkerFaceColor', 'w');

plot(xMean, yMean, 'Marker', '*', 'MarkerSize', 12, 'Mark压力线 erEdgeColor', 'k'); % 得到10.24(c)

f = [0 1 2 3 4 5 4 3 2 1 0];

figure, h_f = plot(f);

Flat = strel([1 1 1]) % 构造平坦(⾼度为0)的结构元素

fd1 = imdilate(f, Flat); % 使⽤平坦的结构元素灰度膨胀

hold on, h_fd1 = plot(fd1, '-ro');

axis([1 11 0 8])

Height = strel([1 1 1], [1 1 1]) % 注意此处strel的⽤灯笼的谜语 法，第⼀个参数的元素为0或1，表⽰结构元素的区域范围(形状)，第⼆个参数表

⽰结构元素中各个元素的⾼度

fd2 = imdilate(f, Height); %使⽤具有⾼度的结构元素的灰度膨胀

hold on, h_fd2 = plot(fd2, '-g*');

legend('原灰度1维函数f', '使⽤平坦结构元素膨胀后', '使⽤⾼度为1的结构元素膨胀后');

I = imread('');

subplot(2, 4, 1), imshow(I, []);%得到图10.39(a)

thresh = graythresh(I) %⾃适应确定阈值

Ibw = im2bw(I, thresh);

subplot女朋友过生日 (2, 4, 2), imshow(Ibw, []);%得到图10.39(b)

subplot(2, 4, 3), surf(double(I(1:8:end,1:8:end))),zlim([0 255]),colormap gray;%显⽰I的3维可视化效果，(c)图

bg = imopen(I,strel('disk',15));%半径为15的圆形结构元素进⾏灰度开微电影大全 运算提取背景曲⾯

subplot(2, 4, 4), surf(double(bg(1:8:end,1:8:end))),zlim([0 255]), colormap gray; %显⽰背搞笑歇后语大全 景曲⾯的三维可视化效果，图10.39(d)

Itophat = imsu青菜疙瘩汤 btract(I, bg); %顶帽变换

subplot(2, 4, 5), imshow(Itophat); %得到图10.39(e)

subplot(2, 4, 6), surf(double(Itophat(1:8:end,1:8:end))),zlim([0 255]); %显⽰顶帽变换图像的三维可视化效果

I2 = imadjust(Itophat);%对⽐度拉伸

subplot(2, 4, 7), imshow(I2); %得到图10.39(f)

thresh2 = graythresh(I2) %⾃适应确定阈值

Ibw2 = im2bw(I2, thresh2); %得到图10.39(g)

subplot(2, 4, 8), 清蒸石斑鱼 imshow(Ibw2); %得到图10.39(h)