16-21.m

来自「matlab6.5图形图像处理一书中所有的源程序」· M 代码 · 共 19 行

I = imread('bonemarr.tif');
imshow(I),title('原图像')
I2 = imhmin(I,50);
% H-minima transform
figure,subplot(221)
imshow(I2),title('H-minima transform')
I2 = imhmax(I,50);
% H-maxima transform
subplot(222),imshow(I2),title('H-maxima transform')
BW = imextendedmax(I,40);
% 扩展极大值变换
subplot(223),imshow(BW),title('Extended-maxima transform')
bw = zeros(size(I));
bw(98:102,101:105) = 1;
% 选择区域
J = I;J(bw ~= 0) = 255;K = imimposemin(I,bw);
% 抑制极小值
subplot(224), imshow(K),title('Impose minima')