% Click mouse to specify initial contour/region

figure('name','Manually initialization');

imshow(img);

BW = roipoly;  % get a region R inside a polygon, BW is a binary image with 1 and 0 inside or

              % outside the polygon;

说明：运行到 BW = roipoly;这句时，figure窗口处于当前状态，用鼠标在窗口中点击一些点，这些点自动连接成一系列直线段，形成一个多边形。最后一个点双击一下，就结束了。得到的结果存储在二值矩阵 BW 中。