function?fd=blanket（A）
M=double（A）;
m=size（M）;
e=101;
U=zeros（m（1）m（2）e）;
B=zeros（m（1）m（2）e）;
%初始化雙毯為圖像灰度表面
for?i=1:1:m（1）
????for?j=1:1:m（2）
????????U（ij1）=M（ij）;
????????B（ij1）=M（ij）;
????end
end
%雙毯生長
for?p=2:1:e-1
????for?i=1:1:m（1）
????????for?j=1:1:m（2）
????????????if?i==1
????????????????a=i;
????????????else?a=i-1;
????????????end
????????????if?j==1
????????????????b=j;
????????????else?b=j-1;
????????????end
????????????if?i==m（1）
????????????????c=i;
????????????else?c=i+1;
????????????end
????????????if?j==m（2）
????????????????d=j;
????????????else?d=j+1;
????????????end
??U（ijp）=max（U（ijp-1）+1max（U（ijp-1）max（max（U（ajp-1）U（ibp-1））max（U（cjp-1）U（idp-1）））））;
??B（ijp）=min（B（ijp-1）-1min（B（ijp-1）min（min（B（ajp-1）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1672??2011-10-06?17:04??blanket.m

-----------?---------??----------?-----??----

?????????????????1672????????????????????1