clear;
close?all;
clc;

%Img=imread（‘three.bmp‘）;??
Img=imread（‘twocells.bmp‘）;??


U=Img（::1）;


Img=double（U）;

[nrowncol]=size（Img）;?

%根據paper設定演化方程中的一些權重參數
epsilon=1.5;
timestep=5;
mu=0.2/timestep;
lambda=5;
nu=1.5;

???
initial=4*ones（nrowncol）;
initial（6:end-56:end-5）=-4;
figure;
imagesc（Img）;colormap（gray）;
hold?on;
[~~]=contour（initial[0?0]‘r‘）;title（‘初始輪廓‘）;


Img_new=conv2（Imgfspecial（‘gaussian‘151.5）‘same‘）;
[xy]=gradient（Img_new）;
g=1./（1+x.^2+y.^2）;
[vxvy]=gradient（g）;

u=initial;
for?n=1:200
????u=evolution（uvxvyglambdamunuepsilontimestep）;?????
????if?mod（n2）==0
????????pause（0.07）;%程序暫停0.07秒用于觀察測試現象，然后繼續執行
????????imagesc（Img）;colormap（gray）;
????????hold?on;
????????[~~]=contour（u[0?0]‘r‘）;?
????????m

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

?????文件?????515899??2018-11-02?13:11??上傳到網上的數字圖像大作業\CV_code_Sample.zip

?????文件???????1761??2019-01-18?17:10??上傳到網上的數字圖像大作業\LevelSet.m

?????文件??????11766??2004-06-09?12:07??上傳到網上的數字圖像大作業\three.bmp

?????文件??????16434??2004-07-23?20:36??上傳到網上的數字圖像大作業\twocells.bmp

?????文件??????77549??2019-01-18?17:06??上傳到網上的數字圖像大作業\水平集變分推導.docx

?????目錄??????????0??2019-01-18?17:07??上傳到網上的數字圖像大作業

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

???????????????623409????????????????????6