朗伯定律的寬觀測帶SAR海冰圖像分割

趙慶平

(1.淮北師范大學物理與電子信息學院，淮北 235000； 2.淮北師范大學信息學院，淮北 235000)

朗伯定律的寬觀測帶SAR海冰圖像分割

趙慶平1,2

(1.淮北師范大學物理與電子信息學院，淮北 235000； 2.淮北師范大學信息學院，淮北 235000)

入射角效應是寬觀測帶SAR海冰圖像分割的主要障礙之一。基于寬觀測帶SAR海冰圖像數據，提出了一種集成余弦朗伯定律的分割算法。為了提高分割算法對SAR海冰圖像的適應性，充分考慮了斑點噪聲和入射角效應因素，并在區域K-means聚類、朗伯定律校正之后，進行區域合并。分別針對合成SAR海冰圖像和星載SAR海冰圖像的實驗結果表明，該算法可有效提高分割的準確性。

入射角效應； 余弦朗伯定律； 寬觀測帶； 海冰； 圖像分割

0 引言

海冰是各種冰類的復雜組合，具有較強的季節性變化，影響著海上交通和海上作業，也是全球氣候變化的一個預測因子，為了更好地了解地球、改善和保護高緯度地區的航線安全，開展海冰的精確化研究和業務化監測十分重要。目前的星載合成孔徑雷達(synthetic aperture Radar，SAR)已經具有了業務化海冰監測的能力，但對于特定的星載SAR成像系統，海冰成像仍存在著不確定性，主要是因為相干斑噪聲和距離向入射角等的影響，前者通常被作為影響SAR海冰圖像分割的主要研究對象[1-3]； 而后者的影響往往被忽略，但是針對于寬觀測帶SAR海冰圖像，該影響不可忽略[4]。由于入射角逐漸變大，會極大地造成寬觀測帶SAR海冰圖像的衰退。基于星載寬觀測帶模式SAR數據，本文提出了一種集成朗伯定律的分割方法。

對于寬觀測帶SAR海冰圖像而言，近距(小入射角)要比遠距(大入射角)亮得多。例如，一幅由ENVISAT ASAR獲取的SAR海冰圖像，其入射角變化范圍達30°，隨著入射角的增加，海冰后向散射每度衰減0.36 dB[5]； 而對于航天飛機獲取的SAR海冰數據而言，入射角的跨度高達40°，入射角的影響是顯而易見的。Ulaby等[6-7]的研究表明，入射角效應不僅取決于入射角，而且取決于海冰的含鹽量、粗糙度和海冰的類型。

本文選擇寬觀測帶SAR海冰圖像為研究對象，綜合考慮了斑點噪聲和入射角效應等不確定因素，從各類海冰固有的特性出發，將冰況圖所包含的信息轉換為海冰的類內特征信息，并利用基于馬爾可夫隨機場(Markov random field，MRF)對二者進行建模，將海冰解譯轉化為目標函數的優化問題。通過區域聚類、余弦朗伯定律校正以及區域合并等步驟實現寬觀測帶SAR海冰圖像的分割。

1 分割算法

若把圖像分成n類。設S是柵格上的一組位置，s是該組中的某個位置，s∈S，xs表示賦給位置s的(冰)類。圖像分割是把原始圖像Is轉換成符號表示[2]

К：{Is|s∈S}→{xs|s∈S}，

(1)

通過上述方式把圖像空間劃分成n個分割段Ω1，…，Ωn，得到

Ωi={s|xs=i，s∈S}，

(2)

(3)

?i≠j，Ωi∩Ωj=φ。

(4)

具體算法流程如圖1所示，其中實線表示本文算法的主要步驟，虛線代表基于區域MRF分割的算法步驟[8]。

圖1 算法流程

1.1 區域化表示

首先,采用分水嶺算法對圖像進行分割，把圖像分割成許多小的區域，每個區域具有相對一致的后向散射值。每個區域Ω由一組位置SΩ構成，這些位置屬于該區域。通過把這些位置組合到對應區域來降低斑點噪聲效應。其次，確定局部強度最小值，定義分水嶺的底，該最小值所對應像素的強度要小于8個近鄰像素。然后，計算圖像梯度。由于在一幅圖像中可能有多個通道，使用矢量場梯度(vector field gradient，VFG)方法利用所有通道來計算聯合圖像梯度。并通過除以場景中的最大梯度進行歸一化。最后，針對每一個殘留像素，利用梯度信息找出圖像中的某個強度最小值。把該最小值的標示符賦給輸出圖像中對應的像素位置。基于平均強度、大小和梯度，把每一對鄰近分水嶺合并在一起。

1.2 基于區域的K-means聚類初始化

作為初始化過程的一部分，區域的K-means聚類是一種推動該算法找出局部可接受方案的技術。利用該聚類算法把標記賦給每個所得區域v。對于每個區域v，K-means聚類算法都能找出最佳的標記，尋找最佳標記的公式為

(5)

1.3 基于朗伯定律校正模型

基于余弦朗伯定律[6，10-11]提出按類校正模型。如果散射目標為完美的朗伯碰撞截面，輻射波全反射回天線，則后向散射系數服從一階余弦定律，即

(6)

對于表面粗糙的散射物，輻射波不能完全反射至天線，則σ0服從余弦定律的平方，即

(7)

如果散射物是分布式和獨立的散射目標，那么

(8)

然而，單一的物理模型不能適用于所有的SAR海冰圖像，綜合式(6)—(8)，定義了一種通用的描述模型，即

(9)

式中m的估計值可以通過規格化曲線得到，該值使得曲線更加平緩。

如果目標物的類型為平地，也就是θa=0，那么式(9)可以簡化為

(10)

為了解決入射角對目標物解譯的影響，參數m可以通過線性回歸模型得到。

s=mt+k，

(11)

式中k為常數。

因此，基于朗伯定律校正模型為

(12)

1.4 區域合并

為了產生最優分割xr*，最小化的成本函數，即

(13)

式中：Xr為所有可能的xr集合；Ef是特征模型能量；Es是空間語境模型能量。

Ef定義為

(14)

式中：c是圖像中的通道數；n表示海冰分割區域的數量;Σi是區域Ωi的特征協方差矩陣;Ef是與假設有關的能量。每一類分水嶺區域的圖像值遵守多變量高斯分布，該分布能給出合理的結果[12-13]。

Es定義為

(15)

式中：β為正權重值，在成本函數中決定了空間語境模型的強度； ?Ωi和?Ωj分別為一組沿著i類和j類邊緣的位置；s為SAR海冰圖像在位置s處的

歸一化梯度；g(s)為邊界懲罰函數，定義為

(16)

式中變量K為SAR圖像的邊界強度[1，14]。

2 實驗分析及評價

2.1 評價標準

整體準確率和Kappa系數[15]是2種常用的分割精度評價參數。整體準確率是被正確分割像素的百分比。Kappa系數定義為

(17)

式中：P(A)為觀測一致率；P(E)為期望一致率。Kappa代表精度評價度量，范圍為[-1，1]。假設分割結果隨機分配。當Kappa= 0時，分割結果與隨機分配一樣好。當Kappa=1時，分割結果是完美的。Kappa為負值表明分割結果相對于正確分割有偏差。

2.2 實驗數據

為了驗證分割的有效性，針對合成SAR海冰圖像和真實星載SAR海冰圖像分別采用本文算法和基于區域MRF分割算法進行對比實驗。合成SAR海冰圖像入射角范圍約為15°，大小為600像素×595像素，并將其劃分為3個相等區段，區段大小為200像素×595像素。該圖含有4類海冰： 尼羅冰、密集冰、稀疏冰和固定冰。真實星載SAR海冰圖像為2000年3月6日獲取的波斯尼亞灣影像，入射角在35°～46°之間，其中海冰入射角的變化范圍為9°，圖像大小為615像素×938像素，同樣劃分為3個相等區段，區段大小為205像素×938像素，波斯尼亞灣一般為短期冰(夏季融化)，參考芬蘭海洋研究所提供的冰況圖和文獻[16]，同樣將其海冰分為尼羅冰、密集冰、稀疏冰和固定冰4類。

2.3 實驗結果分析

針對整體準確率和Kappa系數2個參數，對比利用不同算法2組SAR海冰圖像的分割結果。合成SAR海冰圖像、不同算法分割結果及海面實況如圖2所示。

(a) 合成SAR海冰圖像(b) 區域MRF分割結果(c) 本文算法分割結果(d) 海面實況

圖2 合成SAR海冰圖像分割結果

從圖2可以發現，針對合成SAR海冰圖像，基于區域MRF的分割結果顯示出較多的歧義性，特別是在左、中區段； 而本文算法在該區域兼顧了邊緣定位和區域一致性。該圖像各區段的分割精度對比如表1所示。

表1 合成SAR海冰圖像分割精度對比

從表1可以看出，區域MRF分割算法右區段準確率最高為91.97%； 本文算法的中區段準確率比區域MRF算法最高提升了2.23%，Kappa系數改進了0.029 2。對于整幅SAR海冰圖像來說，由于入射角的跨度增大和其他不確定因素(形狀)的影響，區域MRF分割算法的整幅準確率比該算法各區段最高準確率降低了1.63%； 但是本文分割算法的整幅準確率為91.88%，與各區段準確率基本持平，Kappa系數也提高到了0.914 2，與區域MRF分割算法相比增幅為0.014 6。結果表明，采用本文算法整幅合成海冰圖像的分割精度得到明顯提高。

整幅波斯尼亞灣SAR海冰圖像、不同算法分割結果及海面實況如圖3所示。

(a) 原始SAR海冰圖像(b) 區域MRF分割結果 (c) 本文算法分割結果(d) 海面實況

圖3 波斯尼亞灣SAR海冰圖像分割結果

從圖3可以發現，針對波斯尼亞灣SAR海冰圖像，與基于區域MRF分割結果相比，本文算法分割結果具有更好的區域連通性，圖像中孤立區域明顯減少，更接近海面實際情況，特別是在左、中區段。該圖像各區段的分割精度如表2所示。

表2 波斯尼亞灣SAR圖像海冰分割精度對比

由表2可以得出與合成SAR海冰圖像分割類似的結果。但2幅圖像有所不同的是，對于波斯尼亞灣各區段圖像來說，本文分割算法在中區段準確度上最高僅提升了1.84%，相應的Kappa系數只提升了0.012 1，明顯低于合成SAR圖像的改進精度。然而對于整幅波斯尼亞灣SAR海冰圖像而言，本文算法在整體準確率上提升了3.23%，Kappa系數也提升了0.030 5。波斯尼亞灣海冰分布對應的入射角范圍小于合成SAR海冰圖像可能是該圖像整體分割精度提升低于合成SAR海冰圖像的主要原因。

綜上所示，本文提出的分割算法通過集成余弦朗伯定律，能夠較好地平滑同質海冰的像素值，有效提高SAR海冰圖像的分割精度，但是伴隨著入射角尺度的提高，分割精度提高的效果會逐漸降低。

3 結論

本文針對寬觀測帶SAR海冰圖像提出了一種集成余弦朗伯定律的分割算法，通過對合成SAR海冰圖像和真實星載SAR海冰圖像的對比實驗，證明了本文算法對寬觀測帶SAR圖像分割是有效的，為實現海冰的準確分割提供了一個新途徑。雖然大多數的研究忽視了入射角效應對SAR圖像分割的影響或忽視不同海冰類型對入射角效應的補償，但是筆者認為對于寬觀測帶SAR海冰圖像應該考慮更完美的校正方案，如入射角效應及不同海冰類型的分布。在未來的研究中，將會把本文提出的分割算法應用于其他SAR平臺(如ALOS PALSA和ASAR等采集的)，以增強其理論的可靠性。

[1] Clausi D A,Qin K,Chowdhury M S,et al.MAGIC:Map-guided ice classification system[J].Canadian Journal of Remote Sensing,2010,36(S1):S13-S25.

[2] Karvonen A J.Baltic Sea ice SAR segmentation and classification using modified pulse-coupled neural networks[J].IEEE Transactions on Geoscience and Remote Sensing,2004,42(7):1566-1574.

[3] Soh L K,Tsatsoulis C,Gineris D,et al.ARKTOS:An intelligent system for SAR sea ice image classification[J].IEEE Transactions on Geoscience and Remote Sensing,2004,42(1):229-248.

[4] Leclerc G,Beaulieu N,Bonn F A.A simple method to account for topography in the radiometric correction of radar imagery[J].International Journal of Remote Sensing,2001,22(17):3553-3570.

[5] M?kynen M P,Manninen A T,Simil? M H,et al.Incidence angle dependence of the statistical properties of C-Band HH-Polarization backscattering signatures of the Baltic sea ice[J].IEEE Transactions on Geoscience and Remote Sensing,2002,40(12):2593-2605.

[6] Ulaby F T,Batlivala P P,Dobson M C.Microwave backscatter dependence on surface roughness,soil moisture,and soil texture:Part I-bare soil[J].IEEE Transactions on Geoscience Electronics,1978,16(4):286-295.

[7] Ulaby F T,Bush T F,Batlivala P P.Radar response to vegetation II:8-18 GHz band[J].IEEE Transactions on Antennas and Propagation,1975,23(5):608-618.

[8] Yu Q Y,Clausi D A.SAR sea-ice image analysis based on iterative region growing using semantics[J].IEEE Transactions on Geoscience and Remote Sensing,2007,45(12):3919-3931.

[9] Vincent L,Soille P.Watersheds in digital spaces:An efficient algorithm based on immersion simulations[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,1991,13(6):583-598.

[10]Menges C H,Van Zyl J J,Hill G J E,et al.A procedure for the correction of the effect of variation in incidence angle on AIRSAR data[J].International Journal of Remote Sensing,2001,22(5):829-841.

[11]Menges C H,Hill G J E,Ahmad W,et al.Incidence angle correction of AirSAR data to facilitate land-cover classification[J].Photogrammetric Engineering and Remote Sensing,2001,67(4):479-489.

[12]Yu Q,Clausi D A.IRGS:Image segmentation using edge penalties and region growing[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,2008,30(12):2126-2139.

[13]李旭超,朱善安.圖像分割中的馬爾可夫隨機場方法綜述[J].中國圖象圖形學報,2007,12(5):789-798. Li X C,Zhu S A.A survey of the Markov random field method for image segmentation[J].Journal of Image and Graphics,2007,12(5):789-798.

[14]Qin A K,Clausi D A.Multivariate image segmentation using semantic region growing with adaptive edge penalty[J].IEEE Transactions on Image Processing,2010,19(8):2157-2170.

[15]Richards J A,Jia X P.Remote Sensing Digital Image Analysis[M].Berlin Heidelberg:Springer-Verlag,2006:304-305.

[16]Peter Y.Segmentation of RADARSAT-2 Dual-Polarization Sea Ice Imagery[D].Waterloo:University of Waterloo,2009.

(責任編輯： 陳理)

Wide-swath SAR ice images segmentation based on Lambert’s law

ZHAO Qingping1，2

(1.SchoolofPhysicsandElectronicInformation,HuaibeiNormalUniversity,Huaibei235000,China; 2.InformationCollege,HuaibeiNormalUniversity,Huaibei235000,China)

Incidence angle effect of the SAR images is a major obstacle to the automatic interpretation of SAR sea ice image. Based on wide-swath SAR ice data, this paper proposes a new segmentation algorithm which integrates Lambert’s law correction step. The segmentation algorithm considers the effects of speckle noise and the angle of incidence of factors. The Lambert’s law correction and region merging will be combined. The efficiency of the proposed method has been demonstrated on the segmentation of synthetic SAR sea ice image and gulf of Bothnia SAR sea ice image, where the segmentation accuracy has been substantially improved in contrast to area-based Markov random field(MRF) algorithm.

incident angle effect； Lambert’s cosine law； wide-swath ； sea ice； image segmentation

10.6046/gtzyyg.2017.02.10

趙慶平.朗伯定律的寬觀測帶SAR海冰圖像分割[J].國土資源遙感,2017,29(2)：67-71(Zhao Q P.Wide-swath SAR ice images segmentation based on Lambert’s law[J].Remote Sensing for Land and Resources,2017,29(2)：67-71.)

2016-09-07；

2016-11-04

安徽省高校自然科學研究重點項目“距離相校正的SAR海冰圖像自動分割研究” (編號： KJ2016A650)資助。

趙慶平(1972- ),男,碩士,副教授,主要研究方向為合成孔徑遙感圖像處理與分析、計算機網絡通信等。Email： zhaoqingping1215@163.com。

TP 751.1

A

1001-070X(2017)02-0067-05