顧玲嘉，趙 凱，任瑞治，孫 健，

1.吉林大學(xué)電子科學(xué)與工程學(xué)院，長(zhǎng)春 130012

2.中國(guó)科學(xué)院東北地理與農(nóng)業(yè)生態(tài)研究所，長(zhǎng)春 130012

0 引言

被動(dòng)微波遙感技術(shù)在全球遙感觀測(cè)中有著廣泛的應(yīng)用，它能夠穿透云層獲得地表一定深度的信息[1]。由于星載微波輻射計(jì)無(wú)法裝載大型天線，被動(dòng)微波數(shù)據(jù)的空間分辨率都較低[2]。通常使用的星載被動(dòng)微波輻射計(jì)，如SSM／I、AMSR-E和FY3BMWRI的空間分辨率都是幾十千米[3-6]。對(duì)于觀測(cè)地區(qū)的數(shù)據(jù)，低空間分辨率導(dǎo)致獲得的微波數(shù)據(jù)是其像元內(nèi)多類地物的混合結(jié)果。

經(jīng)前期調(diào)研[7]，針對(duì)星載被動(dòng)微波數(shù)據(jù)空間分辨率低導(dǎo)致的混合像元分問(wèn)題，國(guó)內(nèi)外研究者基于圖像處理技術(shù)解決的方案主要有2類：一種是圖像超分辨率增強(qiáng)方法，另一種是混合像元分解方法。2類方法的實(shí)現(xiàn)都可以提高被動(dòng)微波地物反演參數(shù)，但是處理的方法和最終獲得的數(shù)據(jù)有所不同。

被動(dòng)微波數(shù)據(jù)超分辨率增強(qiáng)方法可以直接提高亮溫的分辨率。其主要的理論依據(jù)是衛(wèi)星的觀測(cè)數(shù)據(jù)腳印之間有很多部分的重疊，利用這些重疊觀測(cè)可以獲得更高分辨率的亮溫?cái)?shù)據(jù)，這種方法也是目前國(guó)內(nèi)外研究者多數(shù)采用的方法。與之相比，混合像元分解方法對(duì)原有被動(dòng)微波數(shù)據(jù)處理后，獲得的是原混合像元中的各組分亮溫，其組分分類根據(jù)具體的應(yīng)用有所不同，組分的獲得需要借助于高于被動(dòng)微波空間分辨率的地物分類數(shù)據(jù)。將組分亮溫?cái)?shù)據(jù)代入到地物參數(shù)反演中，可以提高直接采用混合數(shù)據(jù)進(jìn)行反演的結(jié)果。例如若直接使用水陸邊界處的被動(dòng)微波混合數(shù)據(jù)進(jìn)行海水反演，由于空間分辨率低，該混合像元數(shù)據(jù)是陸地和海洋的綜合結(jié)果，相對(duì)于海洋表面，陸地表面呈現(xiàn)高發(fā)射率，導(dǎo)致反演精度減低。如果通過(guò)微波混合像元分解，將該混合像元內(nèi)的海水組分亮溫提取出來(lái)并進(jìn)行海水參數(shù)的反演，其反演數(shù)據(jù)結(jié)果的精度肯定是要高于前者。

基于地物分類的被動(dòng)微波混合像元分解方法取決于微波像元內(nèi)地物的分類。目前國(guó)內(nèi)外采用的被動(dòng)微波混合方法主要針對(duì)沿海地區(qū)水陸分界比較明顯的區(qū)域進(jìn)行分解，像元內(nèi)的地物分類主要分為海洋和陸地，代表方法是Bellerby等[8]結(jié)合海岸線數(shù)據(jù)利用相鄰像元信息的混合像元分解方法。Bennartz[9]采用周?chē)耆顷懙仡愋偷牧翜財(cái)?shù)據(jù)平均值來(lái)代替陸地組分亮溫[10-13]，然后代入建立的陸地水體混合像元分解模型，求得天線足印中海洋的組分亮溫。MaaB等[14]研究并提出了一種類似于Bennartz的混合像元分解方法，兩者的主要差異在于應(yīng)用領(lǐng)域。MaaB應(yīng)用該方法來(lái)提高海冰參數(shù)，而B(niǎo)ennartz利用該方法實(shí)現(xiàn)大氣參數(shù)的反演[15]。在Bellerby提出的微波混合像元分解算法的基礎(chǔ)上，筆者[16]曾提出基于水陸分類信息的被動(dòng)微波混合像元分解方法，并對(duì)2010年中國(guó)吉林省永吉縣洪澇災(zāi)害進(jìn)行了有效評(píng)估。

上述被動(dòng)微波混合像元分解方法都是利用已知地物分類的先驗(yàn)知識(shí)，建立被動(dòng)微波混合像元分解模型，通過(guò)求解方程組獲得被動(dòng)微波混合像元中水體和陸地的組分亮溫，其有效性的驗(yàn)證都是根據(jù)具體的研究地區(qū)展開(kāi)，在各自研究的地區(qū)內(nèi)可以保證其組分亮溫求解的精度。其中，Bennartz提出的基于陸地組分亮溫的混合像元分解模型相對(duì)簡(jiǎn)單，其陸地組分亮溫是由周?chē)耆顷懙仡愋偷牧翜財(cái)?shù)據(jù)平均值計(jì)算得到；因此，這種模型的適用范圍較小。筆者提出基于地物分類信息的被動(dòng)微波混合像元分解方法，并將其與Bellerby提出的分解方法進(jìn)行對(duì)比研究?？紤]到實(shí)際測(cè)試場(chǎng)地中較難獲得不同地物的被動(dòng)微波組分亮溫?cái)?shù)據(jù)，筆者通過(guò)模擬具有不同組分亮溫分布的觀測(cè)地區(qū)，設(shè)計(jì)具有不同組分亮溫分布的觀測(cè)地區(qū)，對(duì)2種分解方法進(jìn)行評(píng)價(jià)，進(jìn)一步分析2種分解方法的適用范圍和影響其求解精度的主要因素，為建立具有更高求解精度和適用范圍的被動(dòng)微波混合像元分解方法提供理論基礎(chǔ)。

1 基于地物分類信息的被動(dòng)微波混合像元分解方法

Bellerby假設(shè)相鄰9個(gè)像元區(qū)域中陸地和海洋組分的亮溫值相同，這樣會(huì)使計(jì)算結(jié)果產(chǎn)生一定誤差[16]；筆者在 Bellerby算法[17]的基礎(chǔ)上，針對(duì)其算法的局限性進(jìn)行改進(jìn)。已知被動(dòng)微波混合像元亮溫值TB，混合像元中水體分布比例矩陣PW和陸地分布比例矩陣PL，考慮到鄰近區(qū)域內(nèi)水體被動(dòng)微波組分亮溫變化不大的特征，可對(duì)m×n窗口尺寸內(nèi)的被動(dòng)微波混合數(shù)據(jù)進(jìn)行分解，通過(guò)構(gòu)建方程組求解水體組分亮溫TW和陸地組分亮溫TL：

令L（i，j）為水體陸地分類結(jié)果。當(dāng)點(diǎn)（i，j）是陸地上的點(diǎn)時(shí)，L（i，j）為1；當(dāng)點(diǎn)（i，j）是水體的點(diǎn)時(shí)，L（i，j）為0。水體陸地分類結(jié)果來(lái)自于觀測(cè)地區(qū)土地利用分類數(shù)據(jù)。則P（i，j）為天線增益函數(shù)ˉG與水體陸地分類結(jié)果的卷積：

此外，還需設(shè)計(jì)一個(gè)目標(biāo)函數(shù)作為判斷可能解優(yōu)劣的標(biāo)準(zhǔn)。選取F作為目標(biāo)函數(shù)，ξ為事先給定的閾值，其取值取決于對(duì)解的精度要求。目標(biāo)函數(shù)F定義為

該混合像元分解模型中，m×n方程組中有m×n＋1個(gè)未知量，即TL中的m×n個(gè)量和TW。這種方程數(shù)小于未知量數(shù)的方程組屬于欠定方程組，有無(wú)窮多組解，求得最優(yōu)解是該算法需要解決的問(wèn)題。

筆者提出在整個(gè)解空間中搜索最優(yōu)解的計(jì)算方法，對(duì)不同類型的組分亮溫值設(shè)定其對(duì)應(yīng)的取值區(qū)間。具體采用k-means聚類算法對(duì)所有的被動(dòng)微波亮溫?cái)?shù)據(jù)進(jìn)行統(tǒng)計(jì)分類，聚類后的中心值代表不同組分的綜合初值，即CL和CW。根據(jù)陸地被動(dòng)微波混合像元的統(tǒng)計(jì)結(jié)果如均值、最大值、最小值、方差，以及混合像元中不同組分的比例，計(jì)算陸地亮溫的變化閾值RL，確定TL（i，j）初值選取范圍是[CLRL，CL＋RL]；同理，計(jì)算水體亮溫的變化閾值RW，確定TW初值選取范圍是[CW-RW，CW＋RW]。由于水體亮溫應(yīng)小于陸地亮溫，需滿足CW＋RW＜CL-RL，對(duì)初值選取范圍進(jìn)行重新設(shè)定。

根據(jù)目標(biāo)函數(shù)調(diào)整搜索窗口尺寸m×n，采用帶約束條件的非負(fù)最小二乘法迭代運(yùn)算求得最優(yōu)解。由于初值的選取范圍考慮到不同地物的亮溫變化空間，且存在多個(gè)條件限制方程組的求解范圍，因此采用筆者提出的方法可以得到較好的最優(yōu)解。

2 仿真實(shí)驗(yàn)與分析

仿真實(shí)驗(yàn)中，首先設(shè)計(jì)星載被動(dòng)微波天線增益函數(shù)，將其與具有不同模擬組分亮溫值分布的觀測(cè)地區(qū)進(jìn)行卷積，獲得該地區(qū)的微波混合像元亮溫；然后，分別采用Bellerby和筆者提出的分解方法對(duì)混合像元進(jìn)行分解；最后，將其各自分解后得到的組分亮溫與原始組分亮溫進(jìn)行比較，進(jìn)而驗(yàn)證2種被動(dòng)微波混合像元分解方法的有效性。

2.1 微波天線增益函數(shù)

通常，星載微波天線增益函數(shù)近似滿足二維高斯分布。實(shí)際的天線增益方向圖是一個(gè)包括旁瓣的復(fù)雜三維結(jié)構(gòu)，Drusch等[18]研究表明天線函數(shù)可以近似為高斯分布函數(shù)。因此，天線增益函數(shù)的公式定義為

其中：g0和g1是描述高斯分布的常數(shù)；x和y是當(dāng)前位置到測(cè)量中心點(diǎn)的水平和垂直距離；a和b表示-3dB橢圓的半軸，可以根據(jù)被動(dòng)微波數(shù)據(jù)中不同頻率對(duì)應(yīng)的足印來(lái)確定；橢圓半徑r定義為

-3dB點(diǎn)定義為從中心點(diǎn)到天線增益降低到最大值一半的位置，Drusch等[18]將式（6）代入式（5），推導(dǎo)得到

本研究使用-3dB點(diǎn)作為橢圓半徑，則微波天線增益函數(shù)為

2.2 試驗(yàn)樣區(qū)設(shè)計(jì)

假設(shè)觀測(cè)區(qū)域?yàn)?5km×75km的范圍，組分亮溫網(wǎng)格的空間分辨率為5km×5km，則共包含有15×15個(gè)網(wǎng)格。令被動(dòng)微波數(shù)據(jù)的采樣分辨率為25km×25km，考慮到相對(duì)于其他地表類型，水體類型的被動(dòng)微波遙感數(shù)據(jù)亮溫值變化的動(dòng)態(tài)范圍較小，根據(jù)被動(dòng)微波遙感數(shù)據(jù)AMSR-E和MWRI的主要下墊面輻射亮溫統(tǒng)計(jì)特性，分別采用如下方案進(jìn)行樣地組分亮溫設(shè)計(jì)。

方案一：假設(shè)觀測(cè)地區(qū)內(nèi)所有陸地組分亮溫相同，為260K，所有水體組分亮溫相同，為120K。水體組分亮溫在區(qū)域呈十字分布，令水體比例逐漸增加，構(gòu)成被動(dòng)微波組分亮溫。具體組分亮溫分布比例如圖1a、b所示，其中藍(lán)色的圖像位置代表水體，紅色的圖像位置代表陸地，都呈十字分布。由于定義的天線增益函數(shù)中橢圓的長(zhǎng)短半軸均為25km，所以用每5×5網(wǎng)格的組分亮溫與天線增益函數(shù)作卷積計(jì)算，計(jì)算后得到的模擬被動(dòng)微波混合亮溫?cái)?shù)據(jù)為3像元×3像元，具體值如圖1c、d所示。由圖1c、d可見(jiàn)，水體組分亮溫比例的增加使混合像元亮溫?cái)?shù)據(jù)的十字中心位置更接近水體組分的亮溫值，即表現(xiàn)為藍(lán)色圖像。

分別采用Bellerby和筆者提出的方法對(duì)微波混合像元進(jìn)行分解，筆者提出的模型采用的窗口尺寸與Bellerby方法一致，為3像元×3像元，其各組分亮溫的初值選取由地物分類矩陣中各地物的比例動(dòng)態(tài)估計(jì)。將2種模型分解后得到的組分亮溫與原始組分亮溫進(jìn)行比較。實(shí)驗(yàn)結(jié)果表明：在方案一條件下，2種方法求得的組分亮溫與原組分亮溫的誤差為0K，即滿足該條件下，2種方法都可以完全獲得混合像元的組分亮溫。

方案二：假設(shè)觀測(cè)地區(qū)內(nèi)所有陸地組分亮溫相同，為260K；所有水體組分亮溫相同，為120K。水體組分亮溫位置在區(qū)域呈隨機(jī)分布。令水體比例逐漸增加，構(gòu)成被動(dòng)微波組分亮溫。其具體分布如圖2a、b所示。計(jì)算后得到的被動(dòng)微波混合亮溫?cái)?shù)據(jù)如圖2c、d所示。從圖2c、d可見(jiàn)，水體組分亮溫的隨機(jī)分布使得混合像元亮溫?cái)?shù)據(jù)的分布也較隨機(jī)，與方案一中混合像元亮溫呈十字分布有所不同。

實(shí)驗(yàn)結(jié)果表明：在方案二條件下，采用Bellerby方法和筆者提出的方法求得的組分亮溫與原組分亮溫的誤差也為0K。

方案一和方案二的結(jié)論證明：對(duì)于地物類型僅為水體和陸地2類的觀測(cè)地區(qū)，若模型求解區(qū)域?qū)?yīng)的觀測(cè)地區(qū)中陸地亮溫組分相同、水體亮溫組分相同，只要已知地物分類信息，無(wú)論水體的位置和比例如何變化，采用這2種混合像元分解模型均可以從被動(dòng)混合像元中準(zhǔn)確地獲得對(duì)應(yīng)的水體和陸地組分亮溫。

方案三：假設(shè)當(dāng)觀測(cè)地區(qū)沒(méi)有水體時(shí)，陸地組分亮溫每列相同，且從第1列到第15列由246K到260K遞增變化，變化精度為1K（圖3a）。當(dāng)觀測(cè)地區(qū)含有水體時(shí)，根據(jù)相鄰微波亮溫?cái)?shù)據(jù)中水體組分亮溫變化不大的特點(diǎn)，仍假設(shè)觀測(cè)區(qū)域內(nèi)水體組分亮溫相同，為120K，且水體組分亮溫位置與方案二中的分布位置和比例變化相同，即呈隨機(jī)分布并逐漸增加，構(gòu)成如圖3b、c所示的被動(dòng)微波組分亮溫。將亮溫組分?jǐn)?shù)據(jù)和天線增益方程進(jìn)行卷積計(jì)算，得到其對(duì)應(yīng)的被動(dòng)微波混合亮溫如圖3d、e、f所示。圖3d中，由于水體組分亮溫的比例為0，因此構(gòu)成的混合像元亮溫圖仍是陸地分布。隨著水體組分亮溫比例增加，其對(duì)應(yīng)的混合像元亮溫圖3e、f由黃色向藍(lán)色變化，即水體比例的增加使得構(gòu)成的混合像元亮溫?cái)?shù)據(jù)值降低。

觀測(cè)地區(qū)陸地和水體組分亮溫的數(shù)據(jù)統(tǒng)計(jì)結(jié)果和采用2種方法計(jì)算得到的組分亮溫與實(shí)際組分亮溫差值的統(tǒng)計(jì)結(jié)果見(jiàn)表1和表2。

由表1和表2可見(jiàn)，當(dāng)水體比例由0逐漸增加時(shí)：對(duì)于陸地組分亮溫，Bellerby方法得到的組分亮溫絕對(duì)誤差的平均值約為3.1K，而筆者提出方法獲得的組分亮溫絕對(duì)誤差的平均值約為0.1K；對(duì)于水體組分亮溫，Bellerby方法得到的組分亮溫絕對(duì)誤差的平均值約為5.9K，筆者提出方法獲得的組分亮溫?fù)?jù)對(duì)誤差的平均值約為0.2K。當(dāng)觀測(cè)地區(qū)陸地組分亮溫變化時(shí)，由于Bellerby采用的混合像元分解方法假設(shè)3×3混合像元內(nèi)所有陸地組分亮溫相同，所有水體組分亮溫也相同，因此采用最小二乘法求得的陸地組分亮溫精度與3×3混合像元區(qū)域即75km×75km范圍內(nèi)的組分亮溫變化有關(guān)。與之相比，筆者采用的混合像元分解模型假設(shè)3×3混合像元中每個(gè)混合像元內(nèi)的陸地組分亮溫相同，即每25km×25km范圍內(nèi)的陸地組分亮溫一致，同時(shí)假設(shè)所有混合像元中的水體組分亮溫相同，因此采用該方法求得的陸地組分亮溫的誤差主要來(lái)自于25km×25km范圍（即5×5個(gè)組分亮溫網(wǎng)格）內(nèi)陸地組分亮溫的變化。因此，筆者提出的方法能得到更加準(zhǔn)確的分解結(jié)果。

圖1 方案一組分亮溫分布和混合像元亮溫Fig.1 The distribution of component brightness temperature and mixed data in the first project

圖2 方案二中組分亮溫分布和混合像元亮溫Fig.2 The distribution of component brightness temperature and mixed data in the second project

圖3 方案三中組分亮溫分布和混合像元亮溫Fig.3 The distribution of component brightness temperature and mixed data in the third project

表1 方案三陸地組分亮溫?cái)?shù)據(jù)統(tǒng)計(jì)Table1 Statistics results of land component brightness temperature in the third project

表2 方案三水體組分亮溫?cái)?shù)據(jù)統(tǒng)計(jì)Table2 Statistics results of water component brightness temperature in the third project

3 結(jié)論與建議

1）由于陸地、水體的亮度溫度差異很大，且水體亮度溫度比較穩(wěn)定，因此對(duì)于地物類型僅為水體和陸地2類的觀測(cè)地區(qū)，若模型求解區(qū)域?qū)?yīng)的觀測(cè)地區(qū)中陸地、水體亮溫組分相同，只要已知地物分類信息，無(wú)論水體的位置和比例如何變化，采用Bellerby和筆者提出的混合像元分解模型均可以從被動(dòng)混合像元中準(zhǔn)確地獲得實(shí)際的水體組分和陸地組分亮溫。

2）當(dāng)觀測(cè)地區(qū)陸地組分亮溫變化時(shí)，筆者提出的方法進(jìn)行混合像元分解后獲得的組分亮溫的精度要高于Bellerby提出的方法。筆者提出的微波混合像元分解方法中需要對(duì)于欠定方程組進(jìn)行求解，其中各組分亮溫的初值變化范圍和條件對(duì)于求解的精度尤為重要。欠定方程組有無(wú)窮多組解，實(shí)驗(yàn)結(jié)果表明，如果定義好求解條件，可以從多個(gè)解中找到最優(yōu)解。

3）與水體組分亮溫相比，陸地上不同類型地表的亮度溫度差異較大，而且觀測(cè)區(qū)域內(nèi)陸地類型的微波輻射亮度溫度也不完全一致，這是導(dǎo)致目前提出的微波混合像元分解方法求解精度降低的主要原因。通過(guò)對(duì)實(shí)驗(yàn)結(jié)果進(jìn)行分析，如果對(duì)微波混合像元中陸地類型進(jìn)行再次分類，即將陸地根據(jù)地表分類數(shù)據(jù)劃分為不同地表類型，進(jìn)一步求得微波混合像元中多類地物的比例，建立基于多類地物分類的被動(dòng)微波混合像元分解模型，通過(guò)合理定義每種類型地物組分亮溫范圍和設(shè)定求解約束條件，將有可能進(jìn)一步提高被動(dòng)微波混合像元求解的精度，這也是未來(lái)需要研究的工作。

長(zhǎng)春凈月潭遙感實(shí)驗(yàn)站工作人員對(duì)本文仿真實(shí)驗(yàn)方案設(shè)計(jì)提供了幫助和支持，謹(jǐn)致謝忱！

（References）：

[1]Thomas J J.Passive Microwave Remote Sensing for Land Applications[J].Earth and Environmental Science，2008（1）：9-18.

[2]尹紅剛，張德海.使用BG反演算法細(xì)分重建微波輻射計(jì)圖像[J].遙感技術(shù)與應(yīng)用，2006，21（2）：120-124.Yin Honggang，Zhang Dehai.Subdivided Reconstruction of Microwave Radiometer Images with BG Inversion Method[J].Remote Sensing Technology and Application，2006，21（2）：120-124.

[3]Derksen C，Walker A，Ledrew E，et al.Combining SMMR and SSM／I Data for Time Series Analysis of Central North American Snow Water Equivalent[J].Journal of Hydrometeorology，2003，4（2）：304-316.

[4]Paloscia S，Macelloni G，Pampaloni P，et al.Global Scale Monitoring of Soil and Vegetation by Using AMSR-E Multi-Temporal Data [C ]／／International Geoscience and Remote Sensing Symposium（IGARSS）.Denver：IEEE，2006：1609-1612.

[5]楊虎，施建成.FY-3微波成像儀地表參數(shù)反演研究[J].遙感技術(shù)與應(yīng)用，2005，20（1）：194-200.Yang Hu，Shi Jiancheng.On the Estimation of Land Surface Parameters by Using FY-3AMicrowave Radiometer Imager（MWRI）[J].Remote Sensing Technology and Application，2005，20（1）：194-200.

[6]Gu Lingjia，Zhao Kai，Zhang Shuwen，et al.Comparative Analysis of Microwave Brightness Temperature Data in Northeast China Using AMSR-E and MWRI Products[J].Chinese Geographical Science，2011，21（1）：84-93.

[7]顧玲嘉，趙凱，孫健，等.被動(dòng)微波遙感數(shù)據(jù)超分辨率增強(qiáng)與混合像元分解研究綜述[J].遙感技術(shù)與應(yīng)用，2012，27（1）：25-31.Gu Lingjia，Zhao Kai，Sun Jian，et al.Review on the Super Resolution Enhancement and Unmixing Method for Passive Microwave Remote Sensing Data[J].Remote Sensing Technology and Application.2012，27（1）：25-31.

[8]Bellerby T，Taberner M，Wilmshurst A.Passive Microwave Retrieval of Ocean Surface Windspeeds for British Coastal Waters[J].IEEE Transactions on Geoscience and Remote Sensing，1999，37（5）：2578-2584.

[9]Bennartz R.On the Use of SSM／I Measurements in Coastal Regions[J].Journal of Atmospheric and Oceanic Technology，1999，16（4）：417-431.

[10]Bennartz R，Paape K，Gutman G G.Global Land Surface Characterization Using Combined Passive Microwave and Visible／Near-Infrared Satellite Data[J].Eos，Transactions，American Geophysical U-nion，1997，78（1）：59-60.

[11]Grody N C.Surface Identification Using Satellite Microwave Radiometers[J].IEEE Transactions on Geoscience and Remote Sensing，1988，26（6）：850-859.

[12]Jones A S，Vonder Haar T H.Retrieval of Microwave Surface Emittance Over Land Using Coincident Microwave and Infrared Satellite Measurements[J].Journal of Geophysical Research Atmospheres，1997，102（D12）：13609-13626.

[13]Prigent C，Rossow B W，Matthews E.Microwave Land Surface Emissivities Estimated from SSM／I Observations[J].Journal of Geophysical Research，1997，102（D18）：21867-21890.

[14]MaaB N，Kaleschke L.Improving Passive Microwave Sea Ice Concentration Algorithms for Coastal Areas：Applications to the Baltic Sea[J].Tellus Series A：Dynamic Meteorology and Oceanography，2010，62（4）：393-410.

[15]Wessel P，Smith W H F.A Global，Self-Consistent，Hierarchical，High-Resolution Shoreline Database[J].J Geophys Res：B Solid Earth，1996，101（B4）：8741-8743.

[16]Gu LingJia，Zhao Kai，Zhang Shuang，et al.An AMSR-E Data Unmixing Method for Monitoring Flood and Waterlogging Disaster[J].Chinese Geographical Science，2011，21（6）：666-675.

[17]Bellerby T，Taberner M，Wilmshurst A.Retrieval of Land and Sea Brightness Temperatures from Mixed Coastal Pixels in Passive Microwave Data[J].IEEE Transactions on Geoscience and Remote Sensing，1998，36（6）：1844-1851.

[18]Drusch M，Wood E F，Lindau R.The Impact of the SSM／I Antenna Gain Function on Land Surface Parameter Retrieval[J].Geophysical Research Letter，1999，26（23）：3481-3484.