圖像分割

閾值分割方法一般可以分爲全局閾值法和局部閾值法。全局閾值法指利用全

局信息(例如整幅圖像的灰度直方圖)對整幅圖像求出最優分割閾值,可以是單閾

值,也可以是多閾值:局部閾值法是把原始的整幅圖像分爲若干小的子圖像,再

對每個子圖像應用全局閾值法分別求出最優分割閾值。閾值法的最大特點是計算

簡單,在重視運算效率的場合得到了廣泛的應用。

imagesc是軟件Matlab中的一個函數。imagesc(A) 將矩陣A中的元素數值按大小轉化爲不同顏色,並在座標軸對應位置處以這種顏色染色imagesc(x,y,A) x,y決定座標範圍,x,y應是兩個二維向量,即x=[x1 x2],y=[y1 y2],matlab會在[x1,x2]*[y1,y2]的範圍內染色。 如果xy超過兩維,則座標範圍爲[x(1),x(end)]*[y(1),y(end)]

B = sum(A)返回數組中各維的元素和。如果A是一個向量(即一個n1列的矩陣),sum(A)返回這個向量中所有元素的和;如果A是一個矩陣,sum(A)A的各列看做一個向量,並返回一個行向量(即一個1n列的矩陣),這個行向量的第n個元素是A中第n列元素的和;如果A是一個多維數組,sum(A)僅僅計算A中第一個非奇異維,並把它看成一個向量,計算後返回一個行向量組。B = sum(A,dim) 只對A中第dim維的元素進行計算。如果dim1,計算各列的元素之和;如果dim2,則計算A中各行元素之和。

最速下降法

下降法亦稱極小化方法,是一類重要的迭代法。這類方法將方程組求解問題轉化爲求泛函極小問題。

設給出方程組 F(x)=0,其中 https://gss0.bdstatic.com/-4o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D122/sign=6e6c4ce5a3d3fd1f3209a638024c25ce/0824ab18972bd407d91d0c3670899e510eb30927.jpg ,令https://gss2.bdstatic.com/-fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D243/sign=c0168098fd246b607f0eb570d8f91a35/0df431adcbef76093529475425dda3cc7dd99edf.jpg F(x)=0 的充分必要條件是φ(x)=0. φ(x) 的極小點 x* 使φ(x*)=0,則 x* 也是方程組 F(x)=0 的解。只要構造迭代序列 {xk} 使https://gss0.bdstatic.com/94o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D125/sign=da2de00f10d5ad6eaef960e8b4ca39a3/42a98226cffc1e1782b457fa4190f603728de9f6.jpg且滿足https://gss2.bdstatic.com/-fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D107/sign=3dfb12f89945d688a702b6a493c37dab/aa64034f78f0f73651c597ca0155b319eac413c6.jpg就可求得方程組的足夠好的近似解。具體做法是選初始近似 x0,沿一個使φ(x) 下降的方向 p0,令 https://gss0.bdstatic.com/94o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D136/sign=0b1d5f85c18065387feaa010a1dca115/fcfaaf51f3deb48f145a7b15fb1f3a292cf578ce.jpg 然後選步長因子 https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D45/sign=44014fbfa6345982c18ae4970ef4dc88/c9fcc3cec3fdfc0341993873df3f8794a5c2261c.jpg ,使 https://gss3.bdstatic.com/7Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D154/sign=2256ee123201213fcb334ad960e636f8/b2de9c82d158ccbfc3c76ff112d8bc3eb0354186.jpg ,得 https://gss2.bdstatic.com/-fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D104/sign=8e5008e28d18367aa9897bdd1a728b68/b3fb43166d224f4a1ea5b6d802f790529822d139.jpg 一般情況是從 xk 出發,沿φ(x) 下降方向 pk 求出https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D214/sign=d7ee686fafc27d1ea1263cc52fd7adaf/0b55b319ebc4b74540786dddc4fc1e178b821508.jpg https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D124/sign=3d00ed6e017b020808c93be356d8f25f/eaf81a4c510fd9f9842c90a22e2dd42a2934a4c6.jpg 直到

 https://gss1.bdstatic.com/9vo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D71/sign=baf36c79da58ccbf1fbcb73b18d80154/1f178a82b9014a90ff2e8376a2773912b21beeed.jpg 爲止,xm 就作爲方程組解 x* 的近似,上述算法中也可選 uk 使https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D233/sign=4f956015fb1f3a295ec8d2cdaa27bce3/71cf3bc79f3df8dc12c14e5ac611728b46102879.jpg這是一個求一維的極小問題。以上算法即爲下降法。如果 pk 選擇不同,就可得到不同的下降法,特別地,若選 p φ 的負梯度,即 https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D105/sign=9c13b6528926cffc6d2abbb28c004a7d/10dfa9ec8a136327498602f89a8fa0ec09fac796.jpg 則得梯度算法https://gss2.bdstatic.com/9fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D159/sign=6a30e9cbe224b899da3c7d3d57071d59/0df3d7ca7bcb0a46b01748f26063f6246a60afc0.jpg其中https://gss1.bdstatic.com/9vo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D232/sign=c26b0d09b9119313c343f8b357390c10/fcfaaf51f3deb48f1d7c6215fb1f3a292cf578e0.jpg此算法也成最速下降法,此法的優點是計算量少,程序簡單,但收斂慢。在下降法中可去下降方向 pk 爲牛頓方向,即https://gss1.bdstatic.com/-vo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D155/sign=126b0525bf8f8c54e7d3c12a0f282dee/c8177f3e6709c93dcdf5cd3c943df8dcd00054f2.jpg特別地,當 https://gss0.bdstatic.com/-4o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D42/sign=a98fa300b9fb43161e1f7b7822a4938d/b58f8c5494eef01faecbdb0eebfe9925bd317d1a.jpg 時,就得到牛頓法,此外還可取沿座標方向下降的方法,實際上就是一步的 SOR 牛頓法。

共軛梯度法

另一類較重要的下降法爲共軛梯度法。共軛梯度法是最簡單的下降法,早在 1847 年就由法國數學家、力學家柯西 (Cauchy,A.-L.)提出,以後坦普爾 (Temple,G.)、柯里 (Curry,H.) 等人也進行過研究並證明了方法的收斂性。20世紀5060年代,又有不少學者對下降法做了很多研究,提出不少具體算法並建立了收斂性理論,使這類算法在解非線性方程組和最優化計算中得到廣泛應用。 [1] 

參考資料

  • 1.  《數學辭海》總編輯委員會.《數學辭海》第3卷.南京:東南大學出版社,2002.8

 

符號距離函數sign distance function),簡稱SDF,又可以稱爲定向距離函數(oriented distance function),在空間中的一個有限區域上確定一個點到區域邊界的距離並同時對距離的符號進行定義:點在區域邊界內部爲正,外部爲負,位於邊界上時爲0

在歐幾里得空間屬性

如果 https://gss0.bdstatic.com/94o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D14/sign=70d69ce2a718972ba73a04cee7cd4a3b/9922720e0cf3d7cadb0b9f80fe1fbe096a63a943.jpg 是具有分段平滑邊界的歐幾里德空間 https://gss1.bdstatic.com/9vo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D18/sign=686a2771dca20cf44290fad77409d9e7/d53f8794a4c27d1ee3f39b3217d5ad6edcc4386a.jpg 一個子集,那麼帶符號的距離函數幾乎在任何地方都是可微的,並且它的梯度滿足方程

https://gss0.bdstatic.com/-4o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D58/sign=95e16964c111728b342d8c2acafc64df/1c950a7b02087bf44e30b11cfed3572c10dfcf63.jpg

如果 https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D14/sign=ac1206297ac6a7efbd26ac22fcfa3ef8/c995d143ad4bd11338f314c456afa40f4bfb050d.jpg 的邊界對於 https://gss0.bdstatic.com/-4o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D36/sign=4245ee4f92510fb37c197191d8334401/7acb0a46f21fbe09cfeb63e267600c338644ad8b.jpg  https://gss3.bdstatic.com/7Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D17/sign=f62e05949422720e7fcee6fd79cbebd6/5fdf8db1cb134954f6fa48925a4e9258d0094a64.jpg ,那麼在充分接近 https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D14/sign=ac1206297ac6a7efbd26ac22fcfa3ef8/c995d143ad4bd11338f314c456afa40f4bfb050d.jpg 邊界的點上, https://gss2.bdstatic.com/9fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D13/sign=fd706518b88f8c54e7d3c12c3b29379b/a50f4bfbfbedab64516eee36fb36afc379311e23.jpg  https://gss3.bdstatic.com/7Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D17/sign=f62e05949422720e7fcee6fd79cbebd6/5fdf8db1cb134954f6fa48925a4e9258d0094a64.jpg 。尤其在邊界 https://gss3.bdstatic.com/7Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D16/sign=764c2571dca20cf44290fad97409d9cb/d53f8794a4c27d1efdd5993217d5ad6edcc43850.jpg 滿足https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D98/sign=181d51c756afa40f38c6c2d5a96459eb/f7246b600c3387448015c9895d0fd9f9d62aa077.jpg其中 https://gss2.bdstatic.com/9fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D14/sign=3623ead8de1b0ef468e89c5adfc476da/279759ee3d6d55fb201ad4d061224f4a21a4dd6c.jpg 是內部法向矢量場。有符號距離函數因此是法向量場的可微分延伸。特別地, https://gss0.bdstatic.com/-4o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D14/sign=8f98f516fb03918fd3d139ce503d9c3b/a6efce1b9d16fdfa43ee6418b88f8c5495ee7b42.jpg 邊界上的帶符號距離函數的Hessian給出了Weingarten映射。進一步,如果是一個足夠接近Ω邊界的區域,那麼 https://gss3.bdstatic.com/7Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D16/sign=fff895323dadcbef05347a00adafd757/a5c27d1ed21b0ef4348f34a8d1c451da80cb3ee4.jpg 是兩次連續可微分的,那麼就存在一個明確的公式,該公式涉及Weingarten映射 https://gss0.bdstatic.com/-4o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D21/sign=91b21cf3e9dde711e3d244f7a5ef21c6/0eb30f2442a7d9338c454885a14bd11372f00152.jpg 用於根據有符號距離函數改變變量的雅可比矩陣和最近的邊界點。具體地,如果 https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D12/sign=180c9a8974f0f736dcfe48030b55535d/023b5bb5c9ea15ce85bb4cffba003af33b87b2eb.jpg是距離內的點的集合 https://gss3.bdstatic.com/7Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D8/sign=eaec1cf3e9dde711e3d24ec7a778a4/d50735fae6cd7b89a1df9675032442a7d8330e90.jpg 的邊界的 https://gss3.bdstatic.com/7Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D14/sign=f5d33dc947fbfbedd859327b79f05777/242dd42a2834349bac47c273c5ea15ce37d3be87.jpg (即管狀鄰域半徑的 https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D8/sign=40ad5a3bf0039245a5b5ec3e872451/63d9f2d3572c11df2a8f13d56f2762d0f603c2ee.jpg ),並 https://gss1.bdstatic.com/9vo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D12/sign=72616234be119313c343fbb26438c59b/d000baa1cd11728b663220d4c4fcc3cec3fd2c2c.jpg 是上一個絕對積函數Γ,然後

https://gss0.bdstatic.com/-4o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D432/sign=2058bac311950a7b71354fc738d0625c/0d338744ebf81a4c3d5f0d3cdb2a6059242da6d2.jpg

其中 https://gss2.bdstatic.com/9fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D23/sign=e80afb1ffed3572c62e29bdf881383d0/9a504fc2d562853507067ed79cef76c6a6ef6366.jpg 表示行列式 https://gss2.bdstatic.com/-fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D25/sign=090e3d04711ed21b7dc929e0ac6e7710/09fa513d269759eeb3f7cc3dbefb43166c22dfa5.jpg 表示表面積分 [2]  

應用

例如,在計算機視覺中使用簽名的距離函數。如使用了一種方法(由Valve Corporation開發),使用GPU加速來渲染大尺寸(或者高DPI)的平滑字體。Valve的方法計算了光柵空間中的有符號距離場,以避免解決(連續)向量空間中問題的計算複雜性。最近已經提出了分段逼近解決方案(例如,近似具有弧形樣條的Bézier),但即使如此,計算對於實時渲染而言可能太慢,並且必須通過基於網格的離散化技術來近似(並從計算中剔除)到距離太遠的點的距離。

參考資料

  • 1.  [1]董吉文,楊海英,張冰. 水平集方法中符號距離函數的快速生成[J]. 計算機工程與應用,2009,45(34):180-182.
  • 2.  [2] Stanley J. Osher and Ronald P. Fedkiw (2003). Level Set Methods and Dynamic Implicit Surfaces. Springer.

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章