medical.m

来自「提供的文件包含matlab的多个图象处理程序,其中有图象分割,以及图象边缘提取的」· M 代码 · 共 22 行

I = imread('rice.tif');
figure, imshow(I), title('原图');
BWs = edge(I, 'sobel', (graythresh(I) * .1));%边缘提取
figure, imshow(BWs), title('binary gradient mask');
se90 = strel('line', 3, 90);%建立结构函数
se0 = strel('line', 3, 0);%建立结构函数
BWsdil = imdilate(BWs, [se90 se0]);
figure, imshow(BWsdil), title('dilated gradient mask');
BWdfill = imfill(BWsdil, 'holes');
figure, imshow(BWdfill);title('binary image with filled holes');
BWnobord = imclearborder(BWdfill, 4);
figure, imshow(BWnobord), title('清除了边界的图像');
seD = strel('diamond',1);
BWfinal = imerode(BWnobord,seD);
BWfinal = imerode(BWfinal,seD);
figure, imshow(BWfinal), title('segmented image');
BWoutline = bwperim(BWfinal);
Segout = I; 
Segout(BWoutline) = 255; 
figure, imshow(Segout), title('outlined original image');