function?newim?=?createimage（?imuFlag?）
%?CREATEIMAGE?returns?a?color?image?（?3?channel?uint8?image）
%?in?which?the?current?zero_level_set?is?in?red?overlaied
%?onto?the?original?image?（im.）.And?then?re_initializes?u.
[nnynnx]=size（im）;
curvIndex=zeros（5*nnx*nny2）;
curvImag=zeros（[nnynnx]）;
num=0;
for?i?=?2?:?nnx?-?1
?????????for?j?=?2?:?nny?-?1
????????????if?u（ij）<0??&?（u（i+1j）>=0?|?u（i-1j）>=0?|?u（ij+1）>=0?|?u（ij-1）>=0）
???????????????num=num+1;?curvIndex（num1）=i;curvIndex（num2）=j;
???????????????curvImag（ij）=255;
????????????end
?????????end
??????end
tempim?=?im;
tempim（curvImag>0）?=?255;
newim（?:?:?1?）?=?tempim;
tempim（curvImag>0）?=?0;
newim（?:?:?2?）?=?tempim;
newim（?:?:?3?）?=?tempim;
%?convert?to?uint8
newim?=?uint8（?newim?）;
%If?Flag==1re_initialize?u
if?Flag
???new_u?=?zeros（nnxnny）;
???dist=zeros（1num）;
???for?j=1:nny
?????for?i=1:nnx
???????for?k=1:num
??????????dist（k）=sqrt（（i-curvIndex（k1））.^2+（j-curvIndex（k2））.^2）;
???????end
???????new_u（ij）=min（dist）;
???????if?u（ij）<0
??????????new_u（ij）=-new_u（ij）;
???????end
?????end
???end
???u=new_u;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2012-05-18?15:39??GAC\
?????文件???????11274??2012-03-12?15:49??GAC\03.jpg
?????文件??????120054??2006-04-22?13:07??GAC\3.bmp
?????文件????????4063??2012-05-17?19:28??GAC\brain.jpg
?????文件????????1103??2007-11-20?12:07??GAC\createimage.m
?????文件????????7554??2012-05-17?19:29??GAC\flower.jpg
?????文件????????3514??2012-05-18?14:16??GAC\gacgaijinf.m
?????文件????????1073??2006-07-13?18:21??GAC\gauss.m
?????文件????????4282??2012-05-17?19:04??GAC\images?（1）.jpg
?????文件????????8134??2004-07-05?18:59??GAC\twoObj.bmp
?????文件???????16434??2004-07-23?20:36??GAC\twocells.bmp