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;
% 選擇區(qū)域
J = I;J(bw ~= 0) = 255;K = imimposemin(I,bw);
% 抑制極小值
subplot(224), imshow(K),title('Impose minima')