Ch-9.morphological-image-processing.rar_478ch .com_Morphological

%% 第九章 形态学图像处理 %% 例9.1 imdilate 膨胀 clc clear A = imread('..\Pictures\images_ch09\Fig0906(a)(broken-text).tif'); info = imfinfo('..\Pictures\images_ch09\Fig0906(a)(broken-text).tif') B = [0 1 0 1 1 1 0 1 0]; A2 = imdilate(A, B); A3 = imdilate(A2, B); % 二次膨胀 A4 = imdilate(A3, B); % 三次膨胀 % A2 = imdilate(B,A); % 哈哈... imshowMy(A) title('原始图像') imshowMy(A2) title('使用结构元素[B]一次膨胀后的图像') imshowMy(A3) title('使用结构元素[B]二次膨胀后的图像') imshowMy(A4) title('使用结构元素[B]三次膨胀后的图像') imshowMy(A2-A) % 显示增加的部分 title('使用结构元素[B]一次膨胀后和原图像相比较增加的部分') %% 例9.3 imerode clc clear A = imread('..\Pictures\images_ch09\Fig0908(a)(wirebond-mask).tif'); se = strel('disk', 10); imshowMy(A) title('原始图像') A2 = imerode(A, se); imshowMy(A2) title('使用结构元素[disk（10）]腐蚀后的图像') se = strel('disk', 5); A3 = imerode(A, se); imshowMy(A3) title('使用结构元素[disk（5）]腐蚀后的图像') A4 = imerode(A, strel('disk', 20)); imshowMy(A4) title('使用结构元素[disk（20）]腐蚀后的图像') %% 例9.4.1 imopen imclose 注意当结构元素分别使用 square 和 disk 的巨大区别！！！ clc clear f = imread('..\Pictures\images_ch09\Fig0910(a)(shapes).tif'); % se = strel('square', 5); % 结构元素 方型 se = strel('disk', 5); % 结构元素 圆盘形 imshowMy(f) % 原始图 图1 title('原始图像') fo = imopen(f, se); % 开 图2 imshowMy(fo) title('使用结构元素[disk（5）]开操作后的图像') fc = imclose(f, se); % 闭 图3 imshowMy(fc) title('使用结构元素[disk（5）]闭操作后的图像') foc = imclose(fo, se); % 先开再闭 图4 imshowMy(foc) title('使用结构元素[disk（5）]先开操作再闭操作后的图像') fco = imopen(fc, se); % 先闭再开 图5 imshowMy(fco) title('使用结构元素[disk（5）]先闭操作再开操作后的图像') % 先膨胀再腐蚀 图6 fse = imdilate(f, se); figure,set(gcf,'outerposition',get(0,'screensize')) subplot(211),imshow(fse) title('使用结构元素[disk（5）]先膨胀后的图像') fes = imerode(fse, se); subplot(212),imshow(fes) title('使用结构元素[disk（5）]先膨胀再腐蚀后的图像') % 先腐蚀再膨胀 图7 fse = imerode(f, se); figure, set(gcf,'outerposition',get(0,'screensize')) subplot(211),imshow(fse) title('使用结构元素[disk（5）]先腐蚀后的图像') fes = imdilate(fse, se); subplot(212),imshow(fes) title('使用结构元素[disk（5）]先腐蚀再膨胀后的图像') %% 例9.4.2 imopen imclose 指纹 clc clear f = imread('..\Pictures\images_ch09\Fig0911(a)(noisy-fingerprint).tif'); se = strel('square', 3); % 结构元素 % se = strel('disk', 2); % 结构元素 圆盘形 imshowMy(f) % 原始图 title('原始图像') A = imerode(f, se); % 腐蚀 imshowMy(A) title('使用结构元素[square（3）]腐蚀后的图像') fo = imopen(f, se); imshowMy(fo) % 开 title('使用结构元素[square（3）]开操作后的图像') fc = imclose(f, se); % 闭 imshowMy(fc) title('使用结构元素[square（3）]闭操作后的图像') foc = imclose(fo, se); % 先开再闭 imshowMy(foc) title('使用结构元素[square（3）]先开操作再闭操作后的图像') fco = imopen(fc, se); % 先闭再开 imshowMy(fco) title('使用结构元素[square（3）]先闭操作再开操作后的图像') %% 例9.5 bwhitmiss 击中或击不中变换 clc clear f = imread('..\Pictures\images_ch09\Fig0913(a)(small-squares).tif'); imshowMy(f) imfinfoMy(f) title('原始图像') B1 = strel([0 0 0; 0 1 1; 0 1 0]); % 击中：要求“击中”所有的“1”位置，不需要考虑“0”位置 B2 = strel([1 1 1; 1 0 0; 1 0 0]); % 击不中：要求“击不中”所有的“1”位置，不需要考虑“0”位置 B3 = strel([0 1 0; 1 1 1; 0 1 0]); B4 = strel([1 0 1; 0 0 0; 0 0 0]); B5 = strel([0 0 0; 0 1 0; 0 0 0]); % 击中 B6 = strel([1 1 1; 1 0 0; 1 0 0]); % 击不中 g = bwhitmiss(f, B1, B2); imshowMy(g) title('使用结构元素组[1]击中击不中变换后的图像') g2 = bwhitmiss(f, B3, B4); imshowMy(g2) title('使用结构元素组[2]击中击不中变换后的图像') g3 = bwhitmiss(f, B5, B6); imshowMy(g3) title('使用结构元素组[3]击中击不中变换后的图像') %% 利用定义来实现“击中击不中”运算 clc clear f = imread('..\Pictures\images_ch09\Fig0913(a)(small-squares).tif'); imshowMy(f) imfinfoMy(f) title('原始图像') B1 = strel([0 0 0; 0 1 1; 0 1 0]); % 击中 B2 = strel([1 1 1; 1 0 0; 1 0 0]); % 击不中 B3 = strel([0 1 0; 1 1 1; 0 1 0]); B4 = strel([1 0 1; 0 0 0; 0 0 0]); g = imerode(f,B1) & imerode(~f,B2); imshowMy(g) title('使用结构元素组[1]击中击不中变换后的图像') %% makelut clc clear f = inline('sum(x(:)) >= 3'); lut2 = makelut(f,2) lut3 = makelut(f,3) %% 例9.6 makelut 大笑猫之生存游戏 P267 clc clear lut = makelut(@conwaylaws, 3); bw1 = [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 0 0 0 0 1 1 1 1 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 1 1 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ]; imshowMy(bw1, 'notruesize'), title('Generation 1') % 'n' == 'notruesize' % imshowMy(bw1), title('Generation 1') bw2 = applylut(bw1, lut); imshowMy(bw2, 'n'), title('Generation 2') bw3 = applylut(bw2, lut); imshowMy(bw3, 'n'), title('Generation 3') bw4 = applylut(bw3, lut); imshowMy(bw4, 'n'), title('Generation 4') bw5 = applylut(bw4, lut); imshowMy(bw5, 'n'), title('Generation 5') bw6 = applylut(bw5, lut); imshowMy(bw6, 'n'), title('Generation 6') bw7 = applylut(bw6, lut); imshowMy(bw7, 'n'), title('Generation 7') bw8 = applylut(bw7, lut); imshowMy(bw8, 'n'), title('Generation 8') temp = bw1; for i = 2:100; bw100 = applylut(temp, lut); temp = bw100; end imshowMy(bw100, 'n'), title('Generation 100') %% getsequence clc clear se = strel('diamond', 5) decomp = getsequence(se) decomp(1) decomp(2) decomp(3) decomp(4) % decomp(1)*decomp(2)*decomp(3)*decomp(4) 是错误的 % rse = imdilate(imdilate(imdilate(decomp(1), decomp(2)), decomp(3)),... % decomp(4)) 是错误的 %% endpoints clc clear f = imread('..\Pictures\images_ch09\Fig0914(a)(bone-skel).tif'); imshowMy(f) title('原始形态骨骼的图像') g = endpoints(f); imshowMy(g) title('使用函数[endpoints]后得到的端点图像') f = imread('..\Pictures\images_ch09\Fig0916(a)(bone).tif'); imshowMy(f) title('原始骨头图像') g = endpoints(f); imshowMy(g) title('使用函数[endpoints]后得到的端点图像（什么也没有）') %% bwmorph 'remove' 'skel' 'shrink' ‘thin’ clc clear BW = imread('circles.png'); imshowMy(BW); BW2 = bwmorph(BW,'remove'); imshowMy(BW2) BW3 = bwmorph(BW,'skel',Inf); imshowMy(BW3) BW4 = bwmorph(BW,'shrink',Inf); % 去掉毛刺的“骨骼化” imshowMy(BW4) ginf = bwmorph(BW, 'thin', Inf); imshowMy(ginf) title('使用函数[bwmorph]细化到稳定状态后的图像') %% P433 边界提取 bwmorph 'remove' 'skel' 'shrink' ‘thin’ clc clear BW = imread('..\DIP_SourceImage\images_chapter_09\Fig9.14(a).jpg'); imshowMy(BW) title('原始图像') BW2 = bwmorph(BW,'remove'); % 有时达不到理想的提取边界效果 imshowMy(BW2) BW3 = bwmorph(BW,'skel',Inf); imshowMy(BW3) BW4 = bwmorph(BW,'shrink',Inf); % 去掉毛刺的“骨骼化” imshowMy(BW4) ginf = bwmorph(BW, 'thin', Inf); imshowMy(ginf) title('使用函数[bwmorph]细化到稳定状态后的图像') %% bwmorph 'bridge' 'clean' 'hbreak' clc clear f = imread('..\Pictures\images_ch09\Fig0911(a)(noisy-fingerprint).tif'); imshowMy(f) % 原始图 title('