DRLSE 水平集演算法總結

2021-06-22 07:53:03 字數 3373 閱讀 4381

背景:

level set方法是美國數學家osher(加州大學洛杉磯分校)和sethian(加州大學伯克利分校)合作提出的。後者因為對level set的貢獻獲得了去年美國數學會與工業應用數學會聯合頒發的維納獎。遺憾的是這兩位level set的開創者現在正為爭奪level set的名譽而對簿公堂。

level set方法是他們在98年文章「front propagation with curvature depedent speed: algorithms based on hamilton-jacobi formulation"中第提出的。這個方法提出以後被成功地應用於流體力學,計算機圖形學,材料科學等領域。應用於影象處理和計算機視覺始於93年caselles等人和95年malladi等人的兩篇著名文章。他們用level set來表示snakes,想法很簡單:乙個平面上的曲線可以表示成乙個二元函式z=f(x,y)的零點集合(zero level set),既這個二元函式z=f(x,y)所表示的三維曲面與xy平面的交線。更一般地,任何n維曲面都可以表示為乙個n+1維曲面與乙個n維超平面的交集,或稱為n+1維曲面在乙個n維超平面上的投影。相對最早的snake(用引數化的曲線,所以也叫parametric active contour),用level set表示的活動曲線叫著geometric active contours。

用二維曲面與二維平面的交線表示曲線,這在微積分甚至中學數學裡都是很平凡的。但是,當我們要描述曲線運動的時候,用level set表示曲線就有很明顯的優點。比如說,幾條曲線在運動中merge成一條曲線,或一條曲線**成幾條曲線,這樣的拓撲變化是不可能表示成一條連續的引數化曲線的運動。原因很簡單,一條連續的引數化曲線是用乙個一元連續函式來卞表示的,它顯然不能表示幾條分開的曲線(這與連續性矛盾)。

然而,以上所說的曲線的拓撲變化卻可以簡單地表示成乙個連續變化的的曲面與乙個固定的平面的交線。這個曲面本身可以不發生拓撲變化,它可以始終是乙個連續的二元函式z=f(x,y)的圖象。這樣,複雜的曲線運動就可以簡單地表示成乙個更高一維的函式的演化,這可以用乙個發展方程(evolution equation)來描述,數學裡已經有很多任務具可以用了。

再接著說說level set吧。我這次參加cvpr,遇到很多人都說level set很難。我甚至還聽說有個做snake很有名的教授對level set很反感。說實話,我在去年開始實現別人的level set方法時,也對level set越來越討厭。原因是:現有的level set方法實現跟理論不一致(這在gomes和faugeras的文章裡已經指出),需要攙雜不少remedies,比如最令人討厭的re-initialization。還有一些步驟,比如velocity extension,都不是讓人賞心悅目的東西。這些步驟就象長在人身上的贅肉或瘤,也許是良性的,但總是讓人看著極不舒服,甚至怕它惡化。

但是,還是有些人能把這樣複雜而且不是那麼優美的方法實現出來,而且用得不錯,具有其它很多方法不具備的優點,比如讓曲線自然地split和merge。象任何一種理論和方法一樣,儘管有缺點(比如re-initialization),但畢竟大家還是看到了它的巨大潛力。所以這還是成為乙個很熱的研究方向。 

影象分割問題:

偏微分方程影象分割:

近年來,基於偏微分方程的影象分割方法稱為研究熱點,其中kass等在2023年發表的**,snakes:active contour models 在影象分割領域影響深遠。該模型將影象分割問題轉化為乙個能量泛函極小值問題。通過變分法,將泛函極值問題轉化為對偏微分方程的求解。然後把偏微分方程的極小解作為影象分割的結果。目前,偏微分方程的分割方法主要有引數活動輪廓模型和幾何活動輪廓模型。

levelset方法:

level set方法的基本思想是將平面閉合曲線隱含地表達為二維曲面函式的水平集,即具有相同函式值的點集,通過level set函式曲面的進化隱含地求解曲線的運動.儘管這種轉化使得問題在形式上變得複雜,但在問題的求解上帶來很多優點,其最大的優點在於曲線的拓撲變化能夠得到很自然的處理,而且可以獲得唯一的滿足熵條件的弱解.

結合:幾何活動輪廓模型與水平集的結合的曲線演化方法。

該方法利用高維函式曲面,來表達低維的演化曲線或曲面,即,將演化的曲面或曲線嵌入到高維函式表示的曲面中。將演化方程轉化為高維水平集函式的偏微分方程。從未避免引數化過程。故,水平及方法將幾何活動輪廓模型的演化轉化為,水平集函式的偏微分方程的表示式的,數值解的過程。

基本理論:just plain easy to understand: there is a su***ce, it intersects a plane, that gives us a contour and that's it. with image segmentation, the su***ce is updated with forces derived from the image.

符號距離函式的重新初始化:

finet=sgn(fine0)(1-|deltafine|);

fine(x,0)=fine0;

目的:數字糾錯

重新初始化:迭代fine符號距離函式

速度函式:水平集方法進行邊界輪廓提取的關鍵是根據實際問題的需要選取合適的速度函式f,f一般為與影象相關的項(梯度資訊)以及與輪廓曲線的幾何形狀有關的項(曲率)的函式

其中,f表示曲線上各點的演化速度,方向沿著曲線的法線方向,通常與影象梯度和曲線曲率有關。我們想要分析和計算的是在速度f的作用下,曲面隨後的演化情況。從式(2)可以看出,只要速度f變化平滑,則妒(菇,y,t)始終保持是乙個光滑函式,零水平集始終與運動曲面相對應,曲面的拓撲變化可以很容易地描述。對於不同的分割模型,速度f表示式也不相同,從而出現了多種基於水平集的影象分割方法。另外,可以看出式(2)對任意維數的曲面演化都適用。曲線的幾何特性也可以很容易地從水平集函式妒(髫,,r,t)得出,譬如,曲線c上各點的

限差分法,梯度:求導

ms-mumford-shah模型:

包括對區域和邊界的描述

h是hausdorff測度:h0表示點數,h1:表示長度,h2:表示面積

其中第一項稱為資料項或忠誠項,保證近似影象u 保留觀察影象i 的主要資訊;第二項稱為近似影象的正則項,保證影象的光滑性,在觀察影象有雜訊時,正則項有去除雜訊的作用;c 表示目標輪廓的長度;α 、β 為非負常數,平衡各項在能量中的權重。由於mumford-shah 能量泛函難以求解,許多學者給出了近似模型,最簡單的是chan 和vese 提出的兩相分片常數能量泛函:

不需要重新初始化的

DRLSE 水平集演算法總結

背景 level set方法是美國數學家osher 加州大學洛杉磯分校 和sethian 加州大學伯克利分校 合作提出的。後者因為對level set的貢獻獲得了去年美國數學會與工業應用數學會聯合頒發的維納獎。遺憾的是這兩位level set的開創者現在正為爭奪level set的名譽而對簿公堂。l...

演算法總結 並查集

搞了兩年acm,終於準備寫一些部落格了,總得留下點什麼。不了解並查集的話,推薦這個部落格學習,講的很有意思,當時笑出聲。下面寫一下自己的理解吧 並查集最簡單的用法,就是把不同的元素歸併到乙個集合裡。在集合裡選出乙個老大作為代表,集合裡的其他元素都是他的小弟。這樣,就可以判斷任意兩個元素是不是屬於同乙...

並查集演算法總結

並查集可以稱之為不相交集合,在處理查詢幾個元素是否在同一集合時使用並查集可以達到非常快的處理速度,並查集的思想有一點很重要,就是利用樹的思想,表示不同的不相交的集合狀態,利用根節點作為代表元素,來對整個集合的資料進行處理。並查集是一種非常精巧而實用的資料結構,主要用於處理一些不相交集合的合併問題,一...