本發(fā)明屬于圖像分割技術(shù)領(lǐng)域，具體涉及一種基于馬爾可夫模型的合成孔徑雷達圖像快速分割方法。

背景技術(shù)：

合成孔徑雷達(syntheticapertureradar，以下簡稱sar)能夠在全氣候條件下全天提供高分辨率圖像。當前sar圖像分割的主要方法是基于模型驅(qū)動的，通過對分割問題構(gòu)建數(shù)學模型實現(xiàn)sar圖像分割，往往能得到較為理想的分割結(jié)果。但隨著sar成像技術(shù)的發(fā)展，海量的、高分辨率的雷達圖像給分割解譯工作的效率提出了嚴苛的要求，所以如何提高sar圖像分割速度是當前sar圖像分割方法研究的重點。

基于模型驅(qū)動的圖像分割算法，以圖像的先驗知識為基礎，對分割問題進行建模，通常將分割問題轉(zhuǎn)化為優(yōu)化問題。那么提高分割速度的方法就主要為優(yōu)化迭代算法和減少待處理目標數(shù)量兩方面了。優(yōu)化迭代算法是從迭代運算過程入手，對算法本身或者整個流程進行簡化，以達到加速的目的，一般而言難度較大，且適用性不高，容易出現(xiàn)重大問題。減少待處理目標數(shù)量，是算法加速的最直接的方法，其以不改變算法核心部分，實現(xiàn)簡單，成為算法加速中最常用的方法。當確定以減少待處理目標對象數(shù)量來實現(xiàn)算法加速后，如何有效的、合理的實現(xiàn)對象減少就成為了主要問題。

圖像的鄰域系統(tǒng)，也就是圖像中每一像素鄰域關(guān)系的集合，一般而言表現(xiàn)為網(wǎng)格化。網(wǎng)格化的鄰域系統(tǒng)也有利于大量圖像處理算法的實現(xiàn)，尤其是基于模型驅(qū)動的圖像分割算法。對圖像像素進行操作減少算法待處理對象，一般而言有提取代表性對象、歸納圖像先驗信息等方法，也就是說，如何實現(xiàn)利用更少的信息描述原有圖像的信息是加速的關(guān)鍵。而該過程，往往伴隨著原有網(wǎng)格化系統(tǒng)的打破，進而影響到算法的實現(xiàn)。合理的鄰域系統(tǒng)的建立，使得算法加速的同時，最大限度的保有原算法的分割效果。反之，未能充分描述圖像信息則使得分割效果出現(xiàn)嚴重的誤分等情況。對于算法本身而言，如何建立新的鄰域系統(tǒng)，即描述每一對象的鄰域關(guān)系，是算法實現(xiàn)部分難度最大的問題。

技術(shù)實現(xiàn)要素：

本發(fā)明的目的是為了克服上述針對sar圖像分割的不足之處，以達到sar圖像分割速度更快，提供一種基于馬爾可夫模型的sar圖像快速分割方法。馬爾可夫模型，也稱為馬爾可夫隨機場模型(markovrandomfieldmodel，以下簡稱mrf模型)，是指滿足正定性和馬爾可夫性的隨機場模型，它結(jié)合bayes理論，提供了不確定性描述與先驗知識聯(lián)系的紐帶，并利用觀測圖像，根據(jù)統(tǒng)計決策和估計理論中的最優(yōu)準則確定分割問題的目標函數(shù)，求解滿足這些條件的最大可能分布，從而將分割問題轉(zhuǎn)化為最優(yōu)化問題。mrf模型可以將元素的空間關(guān)系緊密結(jié)合起來，能夠充分反應圖像的潛在結(jié)構(gòu)以及圖像的隨機性；通過嚴密的數(shù)學推導和實現(xiàn)，可以得到效果良好的分割結(jié)果。對于基于mrf模型的sar圖像分割算法的加速，以基于聚類思想的超像素算法，作為減少待處理對象的方法，并針對超像素預處理結(jié)果生成新的鄰域系統(tǒng)并結(jié)合mrf模型實現(xiàn)最后的分割。

第一步，利用簡單線性迭代聚類(simplelineariterativeclustering，以下簡稱slic)超像素算法對sar圖像進行預處理，選擇合適的分割步長，保證預處理結(jié)果不出現(xiàn)欠分割和嚴重的過分割現(xiàn)象，該方法為傳統(tǒng)技術(shù)，本發(fā)明直接使用其傳統(tǒng)方法，在此就不再贅述，通過該處理方法獲得分割步長。

本文使用的slic算法，利用顏色相似度和空間距離關(guān)系，以局部迭代聚類實現(xiàn)快速的圖像分割。

第二步，對slic預處理結(jié)果生成偽網(wǎng)格化鄰域系統(tǒng)。

偽網(wǎng)格化(pseudo-meshing)，就是假設slic預處理結(jié)果的所有超像素塊的節(jié)點信息，是按照網(wǎng)格化排列的。首先，根據(jù)超像素分割結(jié)果中圖像像素點信息，即每個像素點的橫縱坐標rk和ck，以及像素點的灰度值gk，對超像素塊的節(jié)點信息進行描述：

d(i,j)＝max(|r(i)-r(j)|,|c(i)-c(j)|)

其中用塊內(nèi)均值表示超像素塊節(jié)點信息，包括橫縱坐標r(i)和c(i)；對于節(jié)點間關(guān)系，考慮到處理速度原因，使用最大距離而不是歐式距離表示節(jié)點間距離d(i,j)。在得到每個超像素塊節(jié)點信息后，按照網(wǎng)格化排列方式，計算每個超像素塊節(jié)點在網(wǎng)格中的位置所在，得到待處理圖像。由于該待處理圖像相比于原圖像要小得多，故簡稱為“小圖”。待所有節(jié)點位置計算完畢，對小圖中出現(xiàn)的空白點進行處理，即將其2階鄰域所有點灰度值均值賦予該空白點。另，對節(jié)點位置計算過程中的重復落點情況進行判斷：若重復落點的兩個節(jié)點的灰度差值在一個識別區(qū)間，則認為這兩個節(jié)點屬同一分割類別，否則相反。識別區(qū)間以圖像灰度最大差值和分割總像素塊數(shù)而決定。

第三步，對slic預處理結(jié)果生成下采樣鄰域系統(tǒng)。

與偽網(wǎng)格化鄰域系統(tǒng)類似，下采樣鄰域系統(tǒng)同樣假設slic預處理結(jié)果中的超像素塊節(jié)點是按照網(wǎng)格化進行排列的。以超像素塊平均面積大小和圖像整體大小為依據(jù)，選擇合適的采樣步長，對slic預處理結(jié)果進行下采樣，直接得到無空白點的待處理小圖。對slic預處理結(jié)果中遺漏的節(jié)點信息進行判斷處理，使遺漏點與其最近中心點進行對比判斷：

其中像素塊大小為ss，識別區(qū)間步長設為j，識別區(qū)間灰度設為g，最近的中心點的距離為lm，最近的中心點的灰度值為gm，該“遺漏”點的灰度值為gt。根據(jù)遺漏節(jié)點的大小和灰度值，以及其與最近節(jié)點的距離和灰度差值，再加上整個分割的識別區(qū)間和步長，做出分割類別的判斷。

第四步，將兩種鄰域系統(tǒng)分別帶入基于mrf模型的分割算法中進行分割操作。

本文采用的mrf模型采用了有限正態(tài)混合模型(finitegaussmaturedmodel，以下簡稱fgmm)描述的特征場，potts模型描述的標號場，迭代條件模式(iterativeconditionmodel，以下簡稱icm)算法作為分割算法。

fgmm模型，是對圖像的每個具有相似特征區(qū)域的局部，建立混合高斯模型。這種方法可以看作是基于局部區(qū)域和基于高斯函數(shù)的聯(lián)合方法，在適用性、復雜性等方面做到了均衡。potts模型的取值范圍是多個值。但和ising模型一樣的是，該模型也規(guī)定每個位置僅有兩種取值。假設圖像中某一點像素為xk，則potts模型中的條件概率可以表示為：

其中，是i像素其鄰域像素的集合，α和β是對應的勢參數(shù)。potts模型的局部概率為：

icm算法主要是對圖像像素逐點更新達到圖像分割的目的。假設圖像y＝{y1,y2,…,yn}的每個像素yi在給定初始分割結(jié)果x的條件下是相互獨立的，且yi關(guān)于x的條件分布只依賴于該像素的標號xi，即

f(yi|x)＝f(yi|xi)

因此，y關(guān)于x的條件分布可以表示為

于是最后的分割結(jié)果就可以表示為

第五步，對分割結(jié)果進行后處理。將小圖分割結(jié)果按照之前鄰域系統(tǒng)生成時的對應關(guān)系，對照回slic預處理結(jié)果，得到兩種鄰域情況下的分割結(jié)果。

第六步，評估兩種鄰域系統(tǒng)不同的分割情況并進行取舍。

對分割結(jié)果進行監(jiān)督的評估，采用以diceratio和sa為主要評估參數(shù)。監(jiān)督評估是將分割結(jié)果和人工標注的真值圖進行對比評估。其中，diceratio是由分割的正確性precision和敏感性recall綜合評估而來。為評判該次分割結(jié)果，引入下列幾個參數(shù)：真陽性(truepositive，簡記為tp)即陽性區(qū)域中的真目標數(shù)量，假陽性(falsepositive，簡記為fp)即陽性區(qū)域中的假目標數(shù)量，真陰性(truenegative，簡記為tn)即陰性區(qū)域中的真目標數(shù)量，假陰性(falsenegative，簡記為fn)即陰性區(qū)域中的假目標數(shù)量。由此，精度precision和敏感性recall就可以定義為：

分割精度sa為

其中，tk表示區(qū)域rk的真值。

根據(jù)上述四個評估參數(shù)的評估情況選取最后的分割結(jié)果。

本發(fā)明的有益效果為：相比于現(xiàn)有技術(shù)，本發(fā)明的方法大幅度提高了圖像分割速度且最大程度地保留了分割質(zhì)量。

附圖說明

圖1為mrf模型中網(wǎng)格化鄰域系統(tǒng)結(jié)構(gòu)；

圖2為slic算法原理圖；

圖3為diceratio監(jiān)督評估原理圖。

具體實施方式

發(fā)明內(nèi)容部分已經(jīng)對本發(fā)明的技術(shù)方案做了詳細描述，在此不再贅述。