【正文】 黑白轉(zhuǎn)換figure, % 新建窗口imshow(BW), % 顯示處理后的圖像title(39。Binary Image39。)。 % 設(shè)置圖像標(biāo)題%step 4 將黑白圖像反色BW = ~ BW。 % 圖像反色figure, % 新建窗口imshow(BW), % 顯示處理后的圖像title(39。Inverted Binary Image39。)。 % 設(shè)置圖像標(biāo)題%step 5 尋找邊界線，參數(shù)‘noholes’將加快圖像的處理速度，它可以防止bwboundaries搜索內(nèi)部邊界 [B,L] = bwboundaries(BW, 39。noholes39。)。 % 尋找邊界%step 6 測定形狀屬性STATS = regionprops(L, 39。all39。)。 % 其實只需要參數(shù)39。BoundingBox39。和39。Extent39。，這樣寫程序簡單點。%step 7 根據(jù)形狀屬性對形狀進(jìn)行分類識別figure, % 新建窗口imshow(RGB), % 顯示原始圖像title(39。Results39。)。 % 設(shè)置圖像標(biāo)題hold on % 繼續(xù)繪圖for i = 1 : length(STATS) W(i) = uint8(abs(STATS(i).BoundingBox(3)STATS(i).BoundingBox(4)) )。 W(i) = W(i) + 2 * uint8((STATS(i).Extent 1) == 0 )。 centroid = STATS(i).Centroid。 switch W(i) case 1 plot(centroid(1),centroid(2),39。wO39。)。 case 2 plot(centroid(1),centroid(2),39。wX39。)。 case 3 plot(centroid(1),centroid(2),39。wS39。)。 endendreturn實驗15：Surface Evolver在材料科學(xué)研究中的應(yīng)用在材料科學(xué)領(lǐng)域, 許多現(xiàn)象取決于材料的表面或分界面的形狀 。 材料的許多特性由其表面的形狀 、 穩(wěn)定性以及表面勢能所決定, 而確定材料表面的形狀通常并非容易 。 為確定材料表面形狀, 有時可采取減少維數(shù)或其它的一些簡化方法來獲得可用的表面信息 。然而這些簡化方法是不可靠的 。 在大多數(shù)情況下, 當(dāng)材料表面的能量已減小到最小時, 其表面將具有確定的形狀 。 因此, 可以通過建立表面能量與表面形狀間的函數(shù)方程來確定表面形狀 。 但建立這樣的函數(shù)方程是非常困難的 。 而且在建立函數(shù)方程的過程中表面張力 、 表面勢能及穩(wěn)定性也應(yīng)當(dāng)考慮到方程中去, 這就進(jìn)一步增加了確定表面形狀的難度 。 對于一些具有簡單幾何形狀的表面, 目前己可以通過建立表面能量減至最小時的函數(shù)方程這一方法來確定表面形狀 。 但對于求解一般形狀問題該方法并不適用 。一種通用的解決具有最小能量的表面問題的方法是使用 Su rface E vo l ve r 軟件( 該軟件也適用于解其它最小值問題) 。Su rface E vo l ve r 軟件是一種通用軟件, 在材料學(xué)中解決上述問題具有靈活性和通用性 。Su rface E vo l ve r 軟件是由美國M i n n e s t o ta 大學(xué)國家科學(xué)和幾何結(jié)構(gòu)計算與可視化技術(shù)研究中心Ken n e th , A. B rak k e 教授研究開發(fā)的 。該軟件是受美國國家科學(xué)基金資助的幾何超級計算工程項目中的能量部分軟件 。 主要應(yīng)用于計算在表面勢能和其他形式的能量作用下的液態(tài)靜平衡形態(tài) 。 目前已經(jīng)得到了廣泛的應(yīng)用, 如: 用于計算火箭燃料箱中低重力作用下的燃料形態(tài) 、 計算毛細(xì)作用表面等 。Su rface E vo l ve r 軟件基于最小能量原理和有限元數(shù)值分析方法, 針對表面成形演變過程分析的一般問題用 C 語言編寫而成, 是一種通用型軟件工具, 軟件本身不具備處理特定問題的能力, 但它提供了靈活的函數(shù)定義方法和能量及約束條件的描述方法 。 使用者可針對自己的特定問題, 按照不同的數(shù)學(xué)模型, 采用基本的幾何元素點 、 邊 、 面 、體，在輸入數(shù)據(jù)文件中對研究對象進(jìn)行合理的定義和描述, 確定研究對象的初始液態(tài)表面以及相關(guān)初始條件 。在應(yīng)用時, 定義初始的幾何條件 、 能量形式和相關(guān)的約束條件, 該軟件計算系統(tǒng)能量最小時的靜平衡形態(tài)和相關(guān)參數(shù) 。 初始幾何條件不僅可以是一維 、 二維 、 三維表面, 而且可以是任意維數(shù)的空間表面。 能量形式包括表面勢能 、 重力勢能 、 用戶定義的各類面積分 、 節(jié)點能等多種形式。 可處理體積約束 、 邊界約束 、方 程約束 、 預(yù)設(shè)平均曲率等約束條件 。 通過運行Su rface E vo l ve r 軟件, 可以用有限元方法求解包含表面勢能 、 重力勢能以及外力作用勢能等作用下的系統(tǒng)在能量最小時的靜平衡狀態(tài), 得出表面形態(tài)的估計結(jié)果, 并能在運行中有效地跟蹤和控制迭代過程 。（或自行copy）。在surface evolver的命令提示符下運行，觀察圖形變化。并試驗運行其它程序例子，在surface evolver ，再自行編寫橢球水滴收縮成圓球形的程序并運行。// // Evolver data for drop of prescribed volume sitting on plane with gravity.// Contact angle with plane can be varied.PARAMETER angle = 90 // interior angle between plane and surface, degreesgravity_constant 0 // start with gravity offdefine WALLT (cos(angle*pi/180)) // virtual tension of facet on plane constraint 1 /* the table top */formula: x3 = 0energy: // for contact anglee1: (WALLT*y)e2: 0e3: 0 vertices1 constraint 1 /* 4 vertices on plane */2 constraint 13 constraint 14 constraint 15 6 7 8 9 fixed /* for table top */10 fixed11 fixed12 fixededges /* given by endpoints and attribute */1 1 2 constraint 1 /* 4 edges on plane */2 2 3 constraint 13 3 4 constraint 14 4 1 constraint 15 5 66 6 7 7 7 8 8 8 59 1 5 10 2 6 11 3 7 12 4 813 9 10 no_refine fixed /* for table top */14 10 11 no_refine fixed15 11 12 no_refine fixed16 12 9 no_refine fixedfaces /* given by oriented edge loop */1 1 10 5 92 2 11 6 103 3 12 7 114 4 9 8 125 5 6 7 87 13 14 15 16 no_refine density 0 fixed /* table top for display */bodies /* one body, defined by its oriented faces */1 1 2 3 4 5 volume 1 density 1readre := {refine edges where on_constraint 1 }// Typical evolutiongogo := { re。 g 5。 r。 g 5。 r。 g 5。 hessian。 hessian。 }// Evolution with 45 degree contact anglegogo2 := { angle := 45。 re。 g 5。 V。V。 r。 g 5。 V。V。 r。 g 5。 hessian。 hessian。 }// Evolution with 90 contact and high gravitygogo3 := { angle := 90。 G 5。 re。 g 5。 r。 g 5。 r。 g 5。 hessian。 hessian。 }// Evolution with 90 contact and negative gravity, . pendant dropgogo4 := { angle := 90。 G 2。 re。 g 5。 r。 g 5。 r。 g 5。 hessian。 hessian。 }// Pendant drop falling off ceilinggogo5 := { angle := 90。 G 5。 re。 g 10。 t .1。 unset vertex constraint 1。 g 3。 }