计算MPP的M函数

%计算MPP的M函数
B=imread('mapleleaf.tif');
subplot(5,2,1),imshow(B),title('原图')

b = boundaries(B, 4, 'cw')
b = b{1};
[M, N] = size(B);
xmin = min(b(:, 1));
ymin = min(b(:, 2));
bim = bound2im(b, M, N, xmin, ymin);
subplot(5,2,2),imshow(bim),title('4连接边界')

[x, y] =  minperpoly(B, 2);
b2 = connectpoly(x, y);