何澤，李世華，2

1.電子科技大學(xué) 資源與環(huán)境學(xué)院,成都 611731;

2.電子科技大學(xué) 長三角研究院（湖州）,湖州 313001

1 引言

水稻生產(chǎn)關(guān)系國計民生，是經(jīng)濟社會發(fā)展和糧食安全的基礎(chǔ)保障，及時準(zhǔn)確的監(jiān)測水稻耕種狀況是農(nóng)業(yè)生產(chǎn)管理的重中之重。傳統(tǒng)的種植范圍劃定和長勢觀測依賴實地調(diào)查，時效性和準(zhǔn)確性較低。衛(wèi)星遙感技術(shù)具有觀測頻次高和覆蓋范圍廣的特點，通過數(shù)理方法和計算機技術(shù)對遙感獲取的水稻波譜響應(yīng)進行分析，可以有效地識別水稻耕作和植株生長狀態(tài)（姚飛和何隆華，2016）。陸表遙感監(jiān)測研究中應(yīng)用最為廣泛的光學(xué)衛(wèi)星影像具有數(shù)據(jù)源多樣、波段信息豐富、成像噪聲小和解譯難度低等優(yōu)點。然而水稻生長在氣候濕熱的季節(jié)，光學(xué)遙感在這一時期往往受云霧氣象干擾，尤其在中低緯度地區(qū)難以獲取高質(zhì)量的光學(xué)影像。以成都平原為例，檢索該地區(qū)2021 年4 月中旬至9月水稻耕作季常用的Landsat 8和Sentinel-2光學(xué)衛(wèi)星數(shù)據(jù)，區(qū)域云量高于30%的影像占比接近90%，并且在6 月至8 月這一水稻快速生長期幾乎沒有可用的光學(xué)數(shù)據(jù)。星載合成孔徑雷達SAR（Synthetic Aperture Radar）通過主動發(fā)射電磁波，記錄地表散射回波特征實現(xiàn)目標(biāo)觀測，可以不受天氣狀況影響持續(xù)對地成像，因此成為水稻遙感監(jiān)測的重要數(shù)據(jù)源。在過去的三十年（1991 年—2021年），星載SAR 傳感器持續(xù)升級迭代，成像性能和數(shù)據(jù)獲取手段不斷優(yōu)化，基于SAR 數(shù)據(jù)的水稻遙感監(jiān)測研究取得了長足的進步（唐鵬欽 等，2009）。

波段、極化和時空分辨率是SAR 成像的關(guān)鍵指標(biāo)，很大程度上決定了SAR 數(shù)據(jù)的應(yīng)用潛力（Zhao 等，2021）。波段表征雷達發(fā)射的電磁波頻率，目前常見的有L波段（如PALSAR-2）、C波段（如Sentinel-1、Radarsat-2 和高分三號）和X 波段（如TerraSAR-X 和Cosmo-Skymed）。不同波段的電磁波擁有不同的波長，因而對探測目標(biāo)的尺寸敏感性不同。C 波段和L 波段電磁波的波長與水稻植株器官（莖桿、葉片和稻穗等）的尺寸接近，對植株生長過程的散射響應(yīng)最為明顯，因此被廣泛用于水稻遙感監(jiān)測（Suga 和Konishi，2008）。極化表征電磁波的電場取向，分為水平極化H（Horizontal）和垂直極化V（Vertical），根據(jù)收發(fā)信號的配置組合產(chǎn)生同極化（HH 和VV）和交叉極化（HV 和VH）通道。星載SAR 傳感器經(jīng)歷了從單極化到全極化成像的發(fā)展歷程。不同極化信號對地物有不同的電磁作用機制，極化方式的增加使SAR 獲取的水稻散射信息更豐富（Yang 等，2012）。時空分辨率表征衛(wèi)星成像的重訪周期和空間精細度，水稻生長迅速且地塊狹小，故而需要高時空分辨率SAR 對地觀測。受相干成像機理影響，SAR 數(shù)據(jù)不可避免地存在斑點噪聲，時空分辨率的提高可以為噪聲抑制提供更多的散射信息參考。

SAR 影像特征參數(shù)主要有后向散射系數(shù)和各種極化參數(shù)，后向散射系數(shù)是最基本的SAR 參數(shù)，表征地物后向散射強度，與水稻介電和幾何特性密切相關(guān)。極化參數(shù)主要來源于緊縮極化模擬、極化分解（Freeman 和Durden，1998；Cloude 和Pottier，1997），以及其他極化運算，如極化相干和相位差（Lopez-Sanchez 等，2012，2014）。通過相干或協(xié)方差矩陣運算可以綜合SAR 數(shù)據(jù)的幅度和相位信息，得到數(shù)理意義明確的極化參數(shù)，更為清晰地描述水稻散射的能量組成。SAR 干涉測量技術(shù)InSAR（Interferometric SAR）通過對不同軌道的SAR 數(shù)據(jù)進行干涉處理獲取地表高程的微小變化，多被用于估算水稻植株高度（Erten 等，2015，2016；Rossi 和Erten，2015；Lopez-Sanchez等，2017）。受數(shù)據(jù)質(zhì)量和反演精度的制約，基于SAR 干涉技術(shù)的水稻監(jiān)測應(yīng)用仍處于實驗探索階段。

目前水稻雷達監(jiān)測研究主要聚焦種植面積提取、生理參數(shù)反演、物候與熟制識別等。種植面積提取研究關(guān)注水稻耕作的空間分布，基于水稻與其他地表覆蓋類型的后向散射特征時空差異進行SAR 影像分類，經(jīng)歷了從先驗閾值決策判別到時空特征機器學(xué)習(xí)的發(fā)展過程。生理參數(shù)反演研究關(guān)注水稻生育相關(guān)參量（如株高、葉面積指數(shù)和生物量），基于植株生理參數(shù)與后向散射特征的時空關(guān)聯(lián)性進行參量映射，方法主要有經(jīng)驗回歸、機器學(xué)習(xí)、機理建模和數(shù)據(jù)同化。物候和熟制識別研究關(guān)注水稻的生育節(jié)律和復(fù)種情況，基于水稻后向散射特征的時間變化規(guī)律進行生長情況判讀，方法主要有時序檢測法和機器學(xué)習(xí)法。

雖然目前水稻遙感監(jiān)測精度和廣度不斷提升，但各研究案例中水稻種植情況、SAR 數(shù)據(jù)屬性、算法模型配置、訓(xùn)練和驗證樣本質(zhì)量均存在較大差異，導(dǎo)致不同技術(shù)思路的實驗結(jié)果難以橫向比較。隨著多源數(shù)據(jù)獲取能力的提高和水稻監(jiān)測對象的多元化，協(xié)同多種技術(shù)方法的復(fù)合監(jiān)測策略被廣泛采用，這種應(yīng)用趨勢雖然可以有效地整合已有研究成果，發(fā)揮不同技術(shù)手段的優(yōu)勢，但也造成算法模型的發(fā)展脈絡(luò)越來越模糊。為了明晰水稻雷達遙感監(jiān)測的研究歷程、現(xiàn)狀與趨勢，本文系統(tǒng)梳理了各研究方向的理論基礎(chǔ)和技術(shù)瓶頸，探討了SAR 數(shù)據(jù)用于水稻遙感監(jiān)測的優(yōu)缺點和存在的科學(xué)問題，對未來的研究重點進行了討論與總結(jié)，以期為相關(guān)方向的科研人員提供參考。

2 水稻種植面積提取

獲取水稻種植面積是生理參數(shù)反演、物候熟制識別和產(chǎn)量預(yù)估工作的基礎(chǔ)，早期水稻雷達監(jiān)測應(yīng)用也主要著力于水稻空間分布信息提取。種植面積提取的本質(zhì)是SAR 影像中水稻像元識別，根據(jù)特定時期水稻與其他地表覆蓋類型的后向散射差異進行影像分類。

波譜信息較少和斑點噪聲較大是制約SAR 影像分類精度的主要問題，波譜信息較少導(dǎo)致水稻與其他植被的影像差異不顯著，斑點噪聲較大導(dǎo)致分類結(jié)果椒鹽噪聲明顯。從數(shù)據(jù)處理的角度，可以通過參數(shù)合成（如波段運算和極化分解）增加波譜信息，利用濾波技術(shù)（空間和時間濾波）平滑噪聲。而從技術(shù)策略的角度，則可以通過多時相和多源觀測來彌補波譜信息不足，利用空間聚類來降低像元異質(zhì)性。

早期的水稻監(jiān)測研究一般選擇在生長季獲取2—5 景影像，通過變化閾值檢測法識別水稻，SAR 觀測頻次增加后，該方法被時序特征識別法取代。隨著機器學(xué)習(xí)技術(shù)的發(fā)展，監(jiān)督分類和非監(jiān)督分類算法被引入SAR 影像分類，其中監(jiān)督分類算法因其較高的靈活性和精確度成為水稻種植面積提取的主流算法，此后學(xué)者分別從算法優(yōu)化（深度學(xué)習(xí)）和特征輔助（多源協(xié)同）的角度不斷提升模型分類精度。非監(jiān)督分類算法精度較低，但基于空間聚類思想發(fā)展出的面向?qū)ο蠓诸惙椒?，是抑制分類噪聲的有效手段（圖1）。

圖1 水稻種植面積提取方法發(fā)展歷程Fig.1 Development of rice planting area identification methods

2.1 時序檢測法

作為人工培育的作物，水稻有著特殊的耕作生長過程，在SAR 影像中表現(xiàn)出區(qū)別于其他地物的散射變化規(guī)律，隨氣候和耕種管理的差異表現(xiàn)出地域性的時序特征。很多研究利用水稻生長的物候先驗信息，針對性地獲取相應(yīng)時相的SAR 影像，通過對時序變化特征的統(tǒng)計分析，識別與水稻生長規(guī)律相符的像元，實現(xiàn)水稻種植面積的提取。

2.1.1 變化閾值檢測

變化閾值檢測是最經(jīng)典的水稻識別方法，通過分析研究區(qū)特定時段內(nèi)水稻像元和其他像元的散射變化差異，確定時變閾值進行決策分類。水稻生育周期可分為移栽期、營養(yǎng)生長期、生殖生長期和成熟期4 個階段（圖2），其中移栽期到營養(yǎng)生長期水稻植株快速生長，是變化閾值監(jiān)測法最關(guān)注的時段。

圖2 水稻生長關(guān)鍵物候期（He等，2018）Fig.2 Photographs for the key rice phenological periods（He et al.，2018）

移栽是水稻特有的栽種方式，水稻先在育秧田中完成從種子到秧苗的發(fā)育，期間對稻作田進行澆灌蓄水，隨后將水稻秧苗從育秧田移栽入稻作田，即為水稻的移栽期。隨后水稻將經(jīng)歷快速的植株生長，地面觀測特征主要有莖稈數(shù)目增加（分蘗）、葉片增多、莖稈直徑和長度增加（拔節(jié)）、植株孕穗，同時田間積水逐漸減少，這一階段稱為水稻的營養(yǎng)生長期。

從移栽期到營養(yǎng)生長期田間水分和植被覆蓋的獨特變化趨勢造就了這一時期水稻遙感信號的顯著特點，在光學(xué)遙感信號中表現(xiàn)為水稻像元綠度加深，植被指數(shù)快速增大，土壤水分指數(shù)逐步降低。對雷達觀測而言，移栽期間的稻田主要表現(xiàn)出田間積水的散射特征，水對微波能量的吸收和側(cè)視成像造成的水體鏡面散射，導(dǎo)致此時雷達接收的回波信號較為微弱。而營養(yǎng)生長期的水稻散射則由植株和土壤的散射能量組成，植株冠層的體散射、粗糙地表的表面散射、植株—土壤之間的多次散射，共同貢獻出較強的回波能量（Yonezawa 等，2012；Corcione 等，2016）。很多學(xué)者通過獲取水稻生長初始階段和高峰階段的雷達數(shù)據(jù)，識別水稻回波信號由弱到強的快速增長特征（圖3），實現(xiàn)了總體精度在80%左右的水稻種植面積提取（Bouvet和Le Toan，2011）。

圖3 各物候階段的水稻SAR散射特征變化（He等，2019）Fig.3 Phenological variation of the SAR polarimetric parameters（He et al.，2019）

變化閾值檢測法簡單高效，適用于地表覆蓋類型簡單、地形平坦的稻作區(qū)，主要的缺點有：（1）依賴研究區(qū)域的水稻物候先驗知識，需要確定適合閾值檢測的最優(yōu)前后時相；（2）需要足夠數(shù)量的水稻及非水稻樣本，用于散射變化差異的統(tǒng)計分析；（3）識別精度受成像噪聲和田間管理（如移栽日期、灌溉方式和施肥選種等）差異的顯著影響。因此，變化閾值檢測法多用于足夠先驗信息支撐下的小樣區(qū)水稻面積快速提取。

2.1.2 時序特征識別

水稻耕作物候歷呈現(xiàn)區(qū)域差異和年際差異，農(nóng)村土地分散經(jīng)營還導(dǎo)致同一地區(qū)不同田塊的水稻生長節(jié)律并不完全同步。高時間分辨率SAR 數(shù)據(jù)可以實現(xiàn)水稻復(fù)雜物候進程的精細追蹤（圖4），通過分析水稻散射的時間變化規(guī)律，提取散射時序曲線特征，如極值、曲率和拐點日期間隔（Phan等，2018），進行閾值規(guī)則判別或相似性度量，從而實現(xiàn)水稻屬性識別。其中相似性度量通過測試未知類型曲線與水稻標(biāo)準(zhǔn)曲線的形態(tài)匹配度，進行水稻散射模式判別，對經(jīng)驗閾值的依賴較小，是更為自動化的水稻識別方法，典型的算法有波譜相似性算法（Yang 等，2018）和時間動態(tài)規(guī)劃算法（Wang等，2020）。

圖4 水稻后向散射密集時序觀測（Nguyen等，2016）Fig.4 Time series VV and VH polarized backscatter over rice fields（Nguyen et al.，2016）

時序特征識別法是高時間分辨率SAR 數(shù)據(jù)應(yīng)用的有效方法，借助水稻生育周期的高頻次雷達觀測，避免了分類時相的人工選擇，水稻識別精度略高于變化閾值檢測法。變化閾值檢測法需要水稻物候先驗知識，分類效果易受閾值選取的影響。而時序相似性度量方法需要采集大量代表性樣本，在構(gòu)建時序標(biāo)準(zhǔn)特征曲線時需考慮復(fù)雜的耕種實際情況，此外逐像元的長時序特征匹配也導(dǎo)致運算量比較高，因此該方法的區(qū)域推廣還需要開展大量的樣本標(biāo)注、分析和驗證工作，同時需要高性能存儲與運算平臺的支撐。

2.2 機器學(xué)習(xí)

水稻種植面積提取的本質(zhì)是圖像分類，因此數(shù)字圖像處理領(lǐng)域的各種機器學(xué)習(xí)算法也被廣泛引入水稻識別研究，并成為目前最常用的面積提取方法，主要包括不需要人工干預(yù)的非監(jiān)督分類方法、人工選擇樣本的監(jiān)督分類方法和自動提取數(shù)據(jù)特征的深度學(xué)習(xí)方法。

2.2.1 非監(jiān)督分類

非監(jiān)督分類不需要地類先驗樣本，在對SAR影像特征進行聚集統(tǒng)計分析的基礎(chǔ)上，自動確定特征參數(shù)和決策規(guī)則，從而實現(xiàn)水稻像元集群的聚類。目前應(yīng)用于水稻種植面積提取的非監(jiān)督分類算法包括K 均值聚類（Gandharum 等，2021）、層次聚類分析（Rudiyanto 等，2019）和迭代自組織數(shù)據(jù)分析（ISODATA）算法（Karila等，2014）。

非監(jiān)督分類方法的優(yōu)點是無須人為干預(yù)、參數(shù)設(shè)置簡單和類內(nèi)同質(zhì)性高，缺點是水稻識別精度不高（一般在80%左右），且分類結(jié)果與地表覆蓋類型的對應(yīng)關(guān)系不明確，需對分類結(jié)果進行人工分析判讀才能確定水稻聚類。SAR 影像“同物異譜—同譜異物”現(xiàn)象和斑點噪聲顯著，水稻極易與其他地類混淆，進一步增加了類別匹配的難度，因此非監(jiān)督方法常用于水稻識別前的地物類別粗提取，作為非水稻（如水體和城鎮(zhèn)）像元剔除和異質(zhì)水稻像元歸并的屬性參考（Rudiyanto 等，2019）。

2.2.2 監(jiān)督分類

監(jiān)督分類是在已知樣本提供類別先驗信息的基礎(chǔ)上，篩選SAR 影像特征參數(shù)進行分離度測試，通過線性或非線性判別算法實現(xiàn)水稻和非水稻像元的分離，通常有以下幾個步驟：（1）分類樣本選取。通過實地調(diào)查或高分辨率光學(xué)影像目視解譯獲得水稻和非水稻樣本點；（2）影像獲取與特征選擇。根據(jù)實驗需求獲取各類成像模式（如多頻段、多極化和多時相）的SAR 數(shù)據(jù)集，提取并篩選出對水稻特征敏感的SAR 參數(shù)（如后向散射系數(shù)和極化分解參數(shù)）；（3）分類器選擇與算法參數(shù)調(diào)整。根據(jù)研究區(qū)和樣本情況選擇合適的分類算法，包括閾值決策（Choudhury 和Chakraborty，2006；Choudhury等，2012）、最大似然（Panigrahy等，1997）、人工神經(jīng)網(wǎng)絡(luò)（Chen 和McNairn，2006）、隨機森 林（Hütt 等，2016；Bazzi 等，2019）和支持向量機（Tan 等，2007；Zhang 等，2009）等；（4）分類后處理與精度驗證。對分類結(jié)果中可能存在的錯分和漏分等制圖誤差進行合理修正，并基于驗證樣本測試分類精度。通過對研究區(qū)的大量代表性地物樣本選擇和SAR 影像特征優(yōu)選，監(jiān)督分類器可以較為準(zhǔn)確的識別出雷達水稻像元。

監(jiān)督分類算法是最主流的水稻種植面積提取方法，其優(yōu)點是樣本特征可人工調(diào)整優(yōu)化、分類過程可控、類別屬性明確且精度較高（一般在90%以上），缺點是分類前需要進行專業(yè)細致的地面調(diào)查或光學(xué)影像判讀，訓(xùn)練樣本的采集和特征參數(shù)的評估工作量較大，樣本的完備性和可靠性易受人為因素影響。此外，SAR 影像波段信息較少，不同地類的散射特征易混淆，往往需要補充多時相數(shù)據(jù)才能獲得較高的分類精度，而最優(yōu)成像日期的選擇需要經(jīng)驗知識和反復(fù)試驗。盡管如此，在遙感大樣本和大數(shù)據(jù)蓬勃發(fā)展的背景下，監(jiān)督分類仍是水稻種植面積提取的首選方法，未來應(yīng)重點完善水稻散射特征的時空譜多維描述，同時提高監(jiān)督分類模型對小樣本的泛化能力。

2.2.3 深度學(xué)習(xí)

作為目前機器學(xué)習(xí)領(lǐng)域的熱點分支，深度學(xué)習(xí)算法在多種信號處理和模式識別任務(wù)中表現(xiàn)卓越，因此也被引入SAR 影像分類和水稻種植面積提取研究。傳統(tǒng)的機器學(xué)習(xí)算法需要人工提取篩選對水稻生長狀態(tài)敏感的SAR 分類特征，而深度學(xué)習(xí)模型是含有多個隱藏層的神經(jīng)網(wǎng)絡(luò)，能夠通過分層處理自動學(xué)習(xí)水稻特征。通過對多極化多時相影像的樣本學(xué)習(xí)，以及對模型的結(jié)構(gòu)參數(shù)優(yōu)化，深度學(xué)習(xí)算法往往能取得極高的水稻識別精度（Ndikumana 等，2018a；Wei 等，2021）。此外還可以補充光學(xué)數(shù)據(jù)，降低斑點噪聲對SAR 影像分類結(jié)果的影響，Zhou 等（2019）使用光學(xué)影像和多尺度圖像分割算法進行地塊識別，在此基礎(chǔ)上將SAR 散射時序特征輸入長短時記憶網(wǎng)絡(luò)LSTM（Long Short-Term Memory network），實現(xiàn)田塊級的水稻種植面積提?。傮w精度=84.2%，Kappa系數(shù)=0.81），優(yōu)于支持向量機（總體精度=75.4%，Kappa系數(shù)=0.72）和隨機森林（總體精度=77.4%，Kappa系數(shù)=0.73）算法。

深度學(xué)習(xí)算法在樣本和模型優(yōu)化的基礎(chǔ)上提升了水稻識別效能，在提取水稻種植面積時能保留田塊的邊界細節(jié)，辨別散射特征混淆的地物，抑制分類結(jié)果的椒鹽噪聲（de Castro 等，2020）。深度學(xué)習(xí)算法的缺點主要是模型設(shè)計復(fù)雜、計算量龐大，對訓(xùn)練數(shù)據(jù)的數(shù)量和質(zhì)量要求高，因而目前多是樣區(qū)尺度的研究。從訓(xùn)練數(shù)據(jù)角度來看，深度學(xué)習(xí)的理論基礎(chǔ)仍然是概率統(tǒng)計，在樣區(qū)研究中取得不錯效果的模型，很難保證基準(zhǔn)數(shù)據(jù)集之外的分類精度，復(fù)雜多樣的水稻生長場景需要使用龐大的標(biāo)注特征來描述，數(shù)據(jù)量過小會造成模型的過擬合。從模型設(shè)計角度來看，深度學(xué)習(xí)超參數(shù)眾多，參數(shù)調(diào)整在很大程度上依賴實驗經(jīng)驗，設(shè)置不同的學(xué)習(xí)率、梯度下降方法、參數(shù)初始化方法、樣本組織方法和歸一化方法等都會影響模型效果，驗證過程需要很高的計算成本（GPU 等硬件算力）和時間成本（數(shù)天至數(shù)月）。此外，深度學(xué)習(xí)的模型結(jié)構(gòu)復(fù)雜，分類結(jié)果的可解釋性不強，模型參數(shù)的含義及其對分類結(jié)果的影響很難明確，因此深度學(xué)習(xí)應(yīng)用于水稻遙感研究時，學(xué)者更為重視雷達數(shù)據(jù)使用策略的改進。傳統(tǒng)機器學(xué)習(xí)的模型結(jié)構(gòu)直觀明晰，還有特征提取和篩選的工作流程，反而帶來了較好的可解釋性。

已有一些方法可以減少深度學(xué)習(xí)算法對訓(xùn)練數(shù)據(jù)的依賴，比如遷移學(xué)習(xí)、少樣本學(xué)習(xí)、無監(jiān)督學(xué)習(xí)和弱監(jiān)督學(xué)習(xí)，雖然其分類性能尚弱于監(jiān)督學(xué)習(xí)，但對于干擾因素眾多、多樣化耕作場景難以窮舉的水稻遙感監(jiān)測而言，仍是具有重要意義的研究方向。隨著衛(wèi)星遙感數(shù)據(jù)獲取和處理能力的快速提升，各種高時空分辨率和覆蓋范圍的土地利用分類產(chǎn)品（如FROM-GLC10 和ESAWorldCover10m）陸續(xù)推出，結(jié)合眾源地理數(shù)據(jù)和高分辨率光學(xué)影像，水稻標(biāo)簽數(shù)據(jù)集的制作與驗證工作將得到有力支撐。模型的云端訓(xùn)練與調(diào)用也是很有前景的應(yīng)用方向，利用遙感云平臺獲取大規(guī)模水稻標(biāo)簽數(shù)據(jù)集，導(dǎo)入深度學(xué)習(xí)云平臺進行模型訓(xùn)練與測試，最后返回遙感云平臺上進行模型調(diào)用與大區(qū)域制圖，可以在一定程度上緩解數(shù)據(jù)集獲取、存儲和計算壓力，提升遙感監(jiān)測時效性。

2.3 面向?qū)ο蟮姆诸愃惴?/h3>

高空間分辨率SAR 影像可以提供更為精細的水稻散射信息，但由于波段信息較少且存在相干斑噪，SAR 影像空間異質(zhì)性也一定程度增加。使用傳統(tǒng)基于像元的影像分類方法進行水稻制圖，難以完整提取地塊邊界輪廓，且識別結(jié)果存在顯著噪聲。面向?qū)ο蟮姆诸愃惴ㄊ墙档头诸愒肼?、提高目?biāo)地物識別完整性的有效手段，目前在水稻識別研究中越來越受到關(guān)注。該方法首先根據(jù)散射信息同質(zhì)性對SAR 像元進行合并和分割，形成由同質(zhì)像元組成的目標(biāo)對象，然后在對象尺度進行散射特征提取與篩選，最后選用適當(dāng)?shù)挠跋穹诸愃惴ǎㄈ绺鞣N機器學(xué)習(xí)方法）實現(xiàn)水稻田塊的整體識別（翟鵬飛 等，2021）。

為保證圖像分割效果，往往需要使用多時相SAR 數(shù)據(jù)，提取多種極化特征（汪小欽 等，2008），或者輔助使用光學(xué)數(shù)據(jù)（Erasmi 和Twele，2009）和地理信息數(shù)據(jù)（杜燁 等，2014）。面向?qū)ο蟮姆诸愃惴梢暂^好地實現(xiàn)SAR 影像中的水稻像元的聚合和判別，提供田塊尺度的水稻制圖結(jié)果（圖5），總體精度略高于基于像素的監(jiān)督分類方法（Pierdicca等，2017）。

圖5 面向?qū)ο蟮乃痉N植面積提取效果（Zhang等，2018）Fig.5 The performance of object-oriented classification for rice planting area identification（Zhang et al.，2018）

面向?qū)ο蟮膱D像分類方法可以充分利用地物的波譜、紋理和拓?fù)湫畔?，刻畫地物的散射、幾何和鄰接屬性，在識別破碎田塊和稻田邊界方面顯示出優(yōu)勢（Gstaiger等，2012；Clauss等，2018）。該方法的缺點是存在最優(yōu)分割尺度選擇的問題，對象分割的質(zhì)量對分類效果有較大的影響。未來應(yīng)關(guān)注不同種植模式下的影像多尺度分割方法，優(yōu)化尺度參數(shù)調(diào)整的效率。

2.4 光學(xué)和SAR數(shù)據(jù)協(xié)同

SAR 數(shù)據(jù)對水稻的介電特性和幾何結(jié)構(gòu)特性敏感，光學(xué)數(shù)據(jù)對水稻的光合組分特性敏感，在成像條件支持的情況下獲取高質(zhì)量光學(xué)衛(wèi)星數(shù)據(jù)，利用光學(xué)精細成像和SAR全天候觀測的優(yōu)勢，可以有效提高水稻識別精度（Villa 等，2015；Mansaray等，2017）。

數(shù)據(jù)協(xié)同策略分為特征級協(xié)同和決策級協(xié)同，其中典型的特征級協(xié)同是將光學(xué)和SAR 影像特征參數(shù)聯(lián)合輸入分類決策器（Torbick 等，2011；Onojeghuo等，2018），這是最常見的多源數(shù)據(jù)使用方案。而決策級協(xié)同一般是使用光學(xué)遙感信息對SAR 分類過程進行指導(dǎo)或修正，比如先使用光學(xué)影像進行地表覆蓋粗分類（Tian 等，2018；Lee 和Lee，2003）或?qū)ο蠓指睿–ai 等，2019），在此基礎(chǔ)上使用SAR 影像特征識別水稻。決策級協(xié)同在技術(shù)路徑上考慮了兩種傳感器的成像信息差異，是更具機理性的方法。此外，非遙感數(shù)據(jù)（如DEM 數(shù)據(jù)和水文數(shù)據(jù)）的補充也對識別精度的提升有積極貢獻（Li 等，2003；Park 等，2018）。光學(xué)與SAR 成像優(yōu)勢互補，二者協(xié)同使用可以彌補單一類型傳感器的觀測局限。由于光學(xué)成像易受云霧干擾，連續(xù)高質(zhì)量的光學(xué)數(shù)據(jù)獲取難以保證，因此未來應(yīng)重點研究如何使用SAR 數(shù)據(jù)對光學(xué)信息時空缺失進行補償，以及如何在決策級協(xié)同框架下使用有限的光學(xué)數(shù)據(jù)最大程度提高SAR 數(shù)據(jù)的水稻制圖效能。

3 水稻生理參數(shù)反演

水稻長勢監(jiān)測是掌握水稻時空分布后的又一項重要工作，水稻生理參數(shù)反演是實現(xiàn)水稻長勢監(jiān)測的具體手段，關(guān)注對象包括葉面積指數(shù)LAI（Leaf Area Index）、光合有效輻射吸收率FPAR（Fraction of Absorbed Photosynthetically Active Radiation）、株高、生物量和含水量等，通過經(jīng)驗擬合和輻射傳輸建模等方法，建立SAR 特征參數(shù)和水稻生理參數(shù)之間的映射關(guān)系和預(yù)測模型（Inoue 和Sakaiya，2013；Inoue 等，2014）。水稻生理參數(shù)和長勢信息不僅能適時指導(dǎo)田間管理（如灌溉、施肥和收割），還能及時對農(nóng)業(yè)病蟲、干旱、洪澇和倒伏災(zāi)害進行預(yù)警和評估，通過對關(guān)鍵指標(biāo)的動態(tài)監(jiān)測，還能實現(xiàn)水稻產(chǎn)量預(yù)測，為糧食政策的制定提供技術(shù)參考。

早期雷達反演研究主要基于概率統(tǒng)計理論，使用經(jīng)驗?zāi)Ｐ停ɑ貧w擬合和機器學(xué)習(xí)）反演水稻生理參數(shù)。為提高反演過程機理性，各種基于輻射傳輸理論的物理模型和半經(jīng)驗?zāi)Ｐ捅惶岢?，揭示了水稻生理參?shù)和微波散射機制的內(nèi)在聯(lián)系。此后又出現(xiàn)了耦合水稻散射模型和生長模型的數(shù)據(jù)同化方法，實現(xiàn)了植株生長機理驅(qū)動下的參數(shù)動態(tài)反演?？紤]到SAR 數(shù)據(jù)無法反映水稻光合組分變化，也有研究在反演過程中加入光學(xué)遙感特征，彌補了雷達的觀測局限性（圖6）。

圖6 水稻生理參數(shù)反演方法發(fā)展歷程Fig.6 Development of rice biophysical parameter retrieval methods

3.1 經(jīng)驗?zāi)Ｐ?/h3>

經(jīng)驗建模是反演水稻生理參數(shù)的經(jīng)典方法，通過分析SAR 影像參數(shù)和水稻生理參數(shù)的統(tǒng)計特征和數(shù)值關(guān)聯(lián)性，建立映射模型實現(xiàn)參數(shù)估算（Chakraborty 等，2005；Chen 等，2009），統(tǒng)計回歸分析和機器學(xué)習(xí)是最常用于確定映射關(guān)系的方法。He 等（2019）使用多時相全極化Radarsat-2數(shù)據(jù)獲取后向散射系數(shù)和極化分解參數(shù)，分析各物候階段與全生育期的水稻LAI 與SAR 參數(shù)的相關(guān)性（圖7），發(fā)現(xiàn)基于Freeman-Durden 分解的雷達植被指數(shù)與水稻LAI的相關(guān)性最佳，可以使用指數(shù)擬合方程描述二者的映射關(guān)系（R2=0.84，RMSE=0.81）。Ndikumana等（2018b）獲取雙極化Sentinel-1數(shù)據(jù)，使用回歸模型和機器學(xué)習(xí)算法反演水稻株高和干生物量，發(fā)現(xiàn)支持向量回歸（株高：R2=0.89，RMSE=9.2 cm；生物量：R2=0.86，RMSE=193 g/m2）和隨機森林（株高：R2=0.91，RMSE=8.2 cm；生物量：R2=0.89，RMSE=174 g/m2）算法的參數(shù)反演效果優(yōu)于多元線性回歸方法（株高：R2=0.79，RMSE=12.4 cm；生物量：R2=0.83，RMSE=213 g/m2），且生物量與VH極化的相關(guān)性很高。

圖7 基于RVIFD（Freeman-Durden分解雷達植被指數(shù)）和HV/VV的水稻LAI經(jīng)驗反演模型（He等，2019）Fig.7 The relationship between rice LAI and RVIFD（Radar Vegetation Index based on Freeman-Durden decomposition）or HV/VV（He et al.，2019）

經(jīng)驗反演方法簡單易行，在訓(xùn)練樣本質(zhì)量足夠高的情況下，能獲得不錯的參數(shù)反演精度，適用于種植模式單一的小樣區(qū)，其弱點是缺乏對散射過程的物理認(rèn)知和機理解釋，普適性較差。當(dāng)研究區(qū)域、成像時期和數(shù)據(jù)來源發(fā)生變化時，需要重新進行參數(shù)訓(xùn)練，因而無法適應(yīng)大范圍復(fù)雜環(huán)境的水稻遙感監(jiān)測，但從經(jīng)驗?zāi)Ｐ瞳@得的參數(shù)響應(yīng)規(guī)律仍可以為水稻散射的機制理解和物理建模提供參考。

3.2 物理模型

物理模型通過前向運算模擬不同空間結(jié)構(gòu)和理化狀態(tài)下的水稻植株和土壤雷達回波，根據(jù)水稻冠層場景建模和輻射傳輸機理建立水稻生理參數(shù)和雷達散射的映射關(guān)系，結(jié)合真實SAR 影像和最優(yōu)化算法，反向運算獲得生理參數(shù)估計值。物理模型能夠較為精確地模擬目標(biāo)散射過程，在解釋水稻后向散射機制方面有重要意義，其對水稻植株結(jié)構(gòu)及散射機制進行不同程度的簡化，將水稻地表器官（如莖稈、葉片和稻穗）描述為具有一定尺寸、形狀、數(shù)量和空間分布規(guī)律的規(guī)則幾何體，根據(jù)生長階段的不同，通常莖稈和稻穗被描述為尺寸不同的圓柱體，而葉片被描述為橢圓盤或針狀體（Koay等，2007；Jia等，2013）。在此基礎(chǔ)上進行輻射傳輸方程解算（Wang 等，2009；Liu 等，2016b）或蒙特卡洛模 擬（Le Toan 等，1997；Yuzugullu 等，2017），從而獲得關(guān)鍵散射類型成分，包括冠層體散射、下墊面散射、冠層與下墊面之間的雙次散射等。

水稻微波散射模型最初多是針對森林場景構(gòu)建，如Karam 模型（Karam 等，1995）、Sun 模型（Sun 等，1991）和MIMICS（Michigan Microwave Canopy Scattering Model）（McDonald和Ulaby，1993），在水稻場景建模中將森林散射模型中的樹干層描述去除，同時調(diào)整冠層組分的介電和結(jié)構(gòu)參數(shù)，即可用于水稻等低矮作物的散射模擬。物理模型描述了電磁傳播過程和能量輻射機理，數(shù)理基礎(chǔ)完備，使水稻生理參數(shù)反演不再盲目，其缺點是算法復(fù)雜，需要獲取植株和土壤的介電參數(shù)，同時對冠層散射體的大小、形狀、方向和空間分布等幾何結(jié)構(gòu)參數(shù)進行測量和假設(shè)，模型的使用難度較大，為適應(yīng)大區(qū)域監(jiān)測中復(fù)雜的水稻耕作情況，物理模型的結(jié)構(gòu)應(yīng)進行適當(dāng)簡化。

3.3 半經(jīng)驗?zāi)Ｐ?/h3>

半經(jīng)驗?zāi)Ｐ蛯λ旧⑸溥^程做更為簡化的描述，既保留對散射機制的物理表達，又兼有易于計算的優(yōu)點，是介于經(jīng)驗?zāi)Ｐ秃臀锢砟Ｐ椭g的折中建模方案。水云模型（Water Cloud Model）是最常用的半經(jīng)驗散射模型（Attema和Ulaby，1978），其將冠層描述為空間中均勻分布的水分粒子云團，將雷達回波分解為植被體散射和植被雙層衰減后的土壤散射兩部分。水云模型中冠層屬性可以借由長勢相關(guān)參數(shù)（如LAI和株高）來表達，下墊面散射貢獻可以使用土壤含水量經(jīng)驗擬合獲得。張曉倩等（2014）考慮水稻生長過程中田間持續(xù)積水的情況，將水云模型中土壤散射成分視為常量，將LAI 作為水稻冠層的量化參數(shù)，使用全生育期LAI 和HH/VV 散射比進行水云模型參數(shù)訓(xùn)練，使用最小二乘法成功反演水稻LAI（R2=0.73，RMSE=0.63）。

水云模型是對植被散射場景的巨大簡化，未考慮到水稻等莖稈作物的冠層結(jié)構(gòu)在水平和垂直方向具有異質(zhì)性，也忽略了水稻植株與下墊面之間的多次散射，因此參數(shù)反演精度受到一定影響，Yang 等（2016）針對這一缺陷對水云模型進行改進，細化了不同物候階段水稻冠層的非均勻空間分布描述，同時考慮了莖稈和稻穗的散射貢獻（圖8），將改進的水云模型和極化分解參數(shù)耦合，基于多時相Radarsat-2數(shù)據(jù)和遺傳算法反演了水稻LAI（R2=0.82，RMSE=0.48）、株高（R2=0.89，RMSE=10.37 cm）、植株體積含水量（R2=0.86，RMSE=0.57 kg/m3）和穗生物量（R2=0.85，RMSE=0.22 kg/m3）等參數(shù)，精度優(yōu)于原始的水云模型，Guo等（2018）驗證了這種改進模型與緊縮極化參數(shù)耦合用于參數(shù)反演的有效性。

圖8 考慮水稻冠層垂直和水平結(jié)構(gòu)異質(zhì)性的改進水云模型散射場景描述（Yang等，2016）Fig.8 The modified water cloud model considering the heterogeneity of the rice canopy in both horizontal and vertical direction（Yang et al.，2016）

半經(jīng)驗?zāi)Ｐ秃啙嵉孛枋隽怂旧⑸錂C制，植株和下墊面理化特性的量化指標(biāo)較少，參數(shù)測量和模型訓(xùn)練都更易實施，在算法復(fù)雜度和機理性方面取得了平衡，可以靈活地與光學(xué)遙感指數(shù)、土壤散射模型以及水稻生長模型耦合使用，在大區(qū)域監(jiān)測中有很大的應(yīng)用潛力。

3.4 數(shù)據(jù)同化

水稻生理參數(shù)反演需要開展與SAR 成像同步的地面觀測實驗，對植株生長狀態(tài)進行動態(tài)采樣，獲得生理參數(shù)變化規(guī)律的統(tǒng)計表達，往往需要等到生長季結(jié)束才能完成全生育周期的參數(shù)采集和模型訓(xùn)練，因此無法及時預(yù)測衛(wèi)星成像下一時刻的水稻生長。

水稻生長模型基于作物生長原理，機理性地描述水稻生長過程與環(huán)境條件的關(guān)系，先于實際觀測給出下一時刻的水稻長勢預(yù)估，可以在單點尺度較好地模擬水稻生理參數(shù)變化。為獲得區(qū)域尺度的生長模型預(yù)測，需要引入數(shù)據(jù)同化技術(shù)，將多時相SAR 數(shù)據(jù)、水稻生長模型和散射模型耦合起來，利用遙感觀測信息不斷進行生長模型的參數(shù)校正或替代，動態(tài)縮小預(yù)測誤差（Setiyono等，2018；宋麗娟 等，2020）。Shen等（2009）以水稻LAI 為媒介耦合水云模型和ORYZA2000 水稻生長模型（Bouman 和van Laar，2006），基于敏感性分析和SCE-UA 全局優(yōu)化方法對ORYZA2000 模型進行參數(shù)調(diào)整和初始化，實現(xiàn)了水稻產(chǎn)量預(yù)估（總體精度=87%）。Liu 等（2016a）使用水云模型和Radarsat-2 數(shù)據(jù)反演水稻生物量，將其同化到WOFOST 生長模型（Savin等，1997），模擬水稻根系干重并計算根系質(zhì)量占比，實現(xiàn)水稻植株重金屬脅迫的監(jiān)測。

數(shù)據(jù)同化方法將SAR 時序觀測與水稻生長模型有效地結(jié)合起來，利用動態(tài)迭代不斷降低生長模型的預(yù)測偏差，使模擬結(jié)果接近真實過程，還能通過生長模型獲得眾多作物狀態(tài)參量。目前遙感數(shù)據(jù)同化研究大多基于光學(xué)數(shù)據(jù)，使用SAR 數(shù)據(jù)與水稻生長模型進行同化反演的研究成果較為零散。這一方面是因為SAR 數(shù)據(jù)的波譜信息較少，與水稻長勢因子的相關(guān)性弱于光學(xué)指數(shù)，另一方面也是受SAR 數(shù)據(jù)較長的重訪周期和較高的獲取成本限制。隨著SAR 成像獲取能力的提高、水稻生長機制和散射機制的研究深入，數(shù)據(jù)同化應(yīng)用將迎來新的發(fā)展契機。

3.5 光學(xué)和SAR數(shù)據(jù)協(xié)同

光學(xué)影像有較高的波譜分辨率和較小的成像噪聲，對水稻植株光譜反射率變化的響應(yīng)較為敏感，而SAR 數(shù)據(jù)對水稻及其下墊面的介電性質(zhì)和幾何結(jié)構(gòu)變化敏感，兩種成像方式優(yōu)勢互補，可以從不同角度反映水稻信息。目前光學(xué)和SAR 數(shù)據(jù)協(xié)同反演水稻生理參數(shù)的研究尚不深入，以特征參數(shù)的簡單組合為主。Alebele 等（2020）基于Sentinel-1/2數(shù)據(jù)構(gòu)造了SAR 和光學(xué)組合植被指數(shù)，作為高斯過程回歸的輸入，用以估計水稻植株器官的生物量，結(jié)果表明使用組合植被指數(shù)的反演精度（R2=0.87，RMSE=0.16 kg/m2）優(yōu)于只使用光學(xué)指數(shù)（R2=0.45，RMSE=1.08 kg/m2）或SAR 指數(shù)（R2=0.60，RMSE=0.24 kg/m2）。

協(xié)同光學(xué)和SAR 數(shù)據(jù)能夠充分利用水稻在不同觀測手段下表現(xiàn)出的遙感特性，更為全面地反映植株及下墊面的理化性質(zhì)。目前數(shù)據(jù)協(xié)同策略多是特征組合，光學(xué)和SAR 特征被視為獨立變量，聯(lián)合輸入回歸方程或機器學(xué)習(xí)模型，這類協(xié)同方式缺乏對異構(gòu)數(shù)據(jù)的關(guān)聯(lián)性挖掘，未考慮兩種遙感特征的機理互補性。同時，特征組合需要在光學(xué)衛(wèi)星可以高質(zhì)量成像的前提下實現(xiàn)，因而不適用于氣象干擾頻繁的地區(qū)。如何綜合利用多源數(shù)據(jù)的時空譜變化規(guī)律，考慮光學(xué)和SAR 特征在時間趨勢和波譜響應(yīng)上的相關(guān)性，耦合光學(xué)和微波輻射傳輸機理，在有限光學(xué)信息輔助下優(yōu)化SAR參數(shù)反演效果值得進一步研究。

4 水稻物候與熟制識別

物候是水稻生長的時間節(jié)律，表征植株發(fā)育和耕作條件的規(guī)律性變化，包括育苗、移栽、分蘗、拔節(jié)、孕穗、抽穗、揚花、灌漿、成熟和收割等具體階段。熟制是稻作區(qū)受氣候和社會發(fā)展水平的影響，耕地田塊的水稻年度種植次數(shù)，主要分為單季稻（如川渝貴地區(qū)）、雙季稻（如粵桂地區(qū)）和三季稻（海南部分地區(qū)）。通過物候信息可以判斷水稻的發(fā)育狀況，指導(dǎo)田間種植管理，而熟制信息可以反映水稻種植強度，為區(qū)域估產(chǎn)提供依據(jù)。物候和熟制信息相互補充，共同描述了一年中各時刻的水稻長勢。

水稻物候和熟制的遙感提取需要借助多時相對地觀測，通過分析水稻散射的時間變化規(guī)律，篩選能夠指示水稻生長狀態(tài)的SAR 特征因子。早期的研究通常將不同物候或熟制的田塊視作不同的地物類型，采用影像分類的思路進行物候和熟制判別。隨著高時間分辨率SAR 數(shù)據(jù)的出現(xiàn)，通過識別季內(nèi)和年內(nèi)的散射時序波動特征提取水稻物候和熟制成為可能（圖9）。

圖9 水稻物候與熟制識別方法發(fā)展歷程Fig.9 Development of rice phenology and cropping intensity identification methods

物候識別方法主要有閾值分割法（Lopez-Sanchez 等，2014）、監(jiān)督分類法（Yuzugullu 等，2015）和時序特征法（Nguyen 等，2016）。閾值分割法分析水稻在多個關(guān)鍵生長期的SAR 特征參數(shù)聚集分布情況，在參數(shù)空間設(shè)置階段閾值實現(xiàn)不同物候期的分割，Yang 等（2014）分析緊縮極化SAR 參數(shù)對水稻物候變化的時域響應(yīng)，獲取物候分割最優(yōu)參數(shù)及經(jīng)驗閾值，采用決策樹方法實現(xiàn)水稻物候識別（總體精度=86%）。監(jiān)督分類法將不同物候期的水稻視為不同地類，篩選對水稻物候變化敏感的SAR 特征，通過對監(jiān)督分類模型進行樣本訓(xùn)練實現(xiàn)物候識別。Kucuk等（2016）從多時相TerraSAR-X 數(shù)據(jù)中提取極化和紋理特征，將水稻生育周期分階段進行類別標(biāo)注，使用支持向量機（總體精度=81.1%，Kappa 系數(shù)=0.76）、k近鄰（總體精度=81.2%，Kappa 系數(shù)=0.76）和決策樹（總體精度=73.2%，Kappa 系數(shù)=0.66）等方法實現(xiàn)了水稻物候期分類識別。時序特征法通過密集SAR 觀測獲取水稻散射時序，使用平滑濾波和多階求導(dǎo)等方法提取時序曲線的峰值、谷值和曲率極值等作為物候指示節(jié)點，Nguyen 等（2016）分析了Sentinel-1 VH 和VV 極化后向散射在水稻生長周期的時間變化規(guī)律，發(fā)現(xiàn)VH 極化受田間積水變化的影響更小，對VH 極化時序進行平滑求導(dǎo)，提取了水稻生長的移栽期（時序谷值）、抽穗期（時序峰值）和耕作季長度（峰谷間隔）等物候特征。

在對大量觀測樣本進行統(tǒng)計分析和模型訓(xùn)練的基礎(chǔ)上，閾值分割法和監(jiān)督分類法能取得不錯的物候識別效果，但其普適性易受各種經(jīng)驗因素的制約。時序特征法依賴高頻次衛(wèi)星觀測，通過對散射時序曲線的構(gòu)建與分析，能夠較好地追蹤水稻的物候發(fā)育過程，然而SAR 時序特征與水稻精細物候的對應(yīng)關(guān)系，目前尚缺乏統(tǒng)一認(rèn)識。

熟制識別方法主要分為多時相分類法和時序特征法，多時相分類法在一年中適宜栽植水稻的時間段內(nèi)獲取多期SAR 影像，通過季節(jié)分段依次分類（Bouvet 等，2009）和多季影像疊合分類（Singha 等，2019）兩種思路識別水稻多熟種植情況。Lasko 等（2018）通過地面調(diào)查獲取水稻樣本熟制信息，對Sentinel-1 數(shù)據(jù)進行多季疊合，使用隨機森林算法實現(xiàn)了水稻熟制提?。傮w精度=93.5%）。時序特征法分析單、雙和三季稻的散射時序變化規(guī)律（圖10），提取時序曲線的物候節(jié)點特征（如峰谷值）以識別單個生長周期，統(tǒng)計物候節(jié)點重復(fù)次數(shù)實現(xiàn)耕作輪次的判斷，精度一般在85%左右（Nguyen 等，2015）。多時相分類法對SAR 觀測頻次的要求不高，且影像分類算法較為成熟，因而實際應(yīng)用更為普遍。時序特征法適用于成像時間密集的SAR 數(shù)據(jù)，在峰形識別中需要解決成像噪聲和輪種、間種作物造成的偽波峰問題，而多時相分類法的技術(shù)關(guān)鍵是分類時相和分類樣本的合理選擇。

圖10 多季稻SAR散射時序規(guī)律（Nguyen等，2015）Fig.10 Time series SAR backscatter behavior in relation to rice growing season（Nguyen et al.，2015）

水稻散射與植株生育狀態(tài)息息相關(guān)，高頻次對地觀測可以提供較為連續(xù)的水稻散射時序信息，時序曲線的季相波動可以表征水稻物候和熟制，因此未來應(yīng)重點開展時序特征提取和模式識別工作，分析水稻生長機理和散射機制的內(nèi)在聯(lián)系，完善復(fù)雜種植條件下水稻散射時序的機理性表達，逐步降低對經(jīng)驗樣本的依賴。

5 討論

基于SAR 數(shù)據(jù)的水稻遙感監(jiān)測技術(shù)在過去三十年（1991年—2021年）取得了長足的發(fā)展，這主要歸功于SAR 成像能力增強和監(jiān)測方法進步?？臻g分辨率的提高使監(jiān)測尺度更加精細，時間分辨率的提高使生長動態(tài)捕獲更加密集，極化模式的增加使水稻散射機制理解更加全面。監(jiān)測方法的進步不僅使建模決策從經(jīng)驗分析發(fā)展到機理推導(dǎo)，更將數(shù)學(xué)和計算機領(lǐng)域的前沿算法引入了水稻遙感監(jiān)測。

5.1 困難與挑戰(zhàn)

雖然SAR 數(shù)據(jù)已成功應(yīng)用于水稻生長監(jiān)測的諸多方面，但要準(zhǔn)確描述水稻散射的時空譜特性仍需解決多項難題：

（1）破碎地塊和復(fù)雜地形。平原地區(qū)規(guī)?；N植是水稻遙感監(jiān)測的理想場景，然而以中國南方丘陵地區(qū)為代表的稻作區(qū)地塊破碎且地形復(fù)雜，這類區(qū)域同樣有大量稻田分布。這些稻田形狀大小各異，面積最小處寬度僅數(shù)米，毗鄰道橋、屋舍、溝渠、園林和坡埂等非農(nóng)耕地類，同時耕地分散經(jīng)營模式導(dǎo)致種植類型多樣、稻田分布零散。目前常用的SAR 衛(wèi)星影像空間分辨率多為10 m 量級，稻田像元極易受到非水稻地物的散射干擾，且SAR 斑點噪聲會進一步劣化破碎田塊的水稻散射特性表達。

丘陵地區(qū)稻作耕地受地形起伏影響，鄰近田塊存在明顯的高程差異，山嶺地區(qū)為提高土地利用率還存在依托山體開墾梯田的情況，鄰近高差不僅使相鄰田塊的雷達局部入射角存在差異，還導(dǎo)致田塊邊界出現(xiàn)高度不一的坡埂，這些土質(zhì)立面與水稻下墊面形成朝向各異的二面角結(jié)構(gòu)，對SAR 側(cè)視成像信號造成了嚴(yán)重干擾，此外SAR 影像會在地勢起伏區(qū)域出現(xiàn)疊掩和陰影等異常散射，進一步干擾水稻散射特征獲取。

（2）稻田栽培條件多樣性。稻田栽培條件包括田間管理和植株性狀兩個方面，田間管理多樣性主要指植秧技術(shù)（育苗、耕犁和移栽）、栽作技術(shù)（株簇密度和壟行朝向）、管水技術(shù)（淹水、灌溉和排水）、施肥技術(shù)（肥料類型和用量）、除草技術(shù)（除草方式和頻次）、曬田技術(shù)（曬田日期和時長）、病害管控（蟲害和倒伏防治）水平的田塊間差異。植株性狀多樣性主要是由水稻品種和水肥差異導(dǎo)致的植株地表器官（莖葉穗）數(shù)量、尺寸、含水量和質(zhì)地的差異，以及發(fā)育速度和生育周期的區(qū)別。

稻田散射特征受下墊面和植株的理化特性共同影響，田間管理和植株性狀的多樣化影響因子會顯著增加水稻散射的描述難度。雖然物理散射模型和作物生長模型的同化技術(shù)可以模擬SAR 信號對稻田栽培條件的響應(yīng)，但受限于模型參數(shù)獲取和算法復(fù)雜度，其應(yīng)用于生產(chǎn)實際的難度還很大。

（3）季候異步和輪作間作。受地形、水文、氣候、社會發(fā)展水平以及耕作分散經(jīng)營等因素的影響，不同地區(qū)的水稻熟制（單、雙、三季稻和再生稻等）、輪作（同一田塊輪流種植不同種類的作物）、間作（田壟邊界種植不同種類的作物，或田間混合種植不同品種的水稻）方式存在差異，耕作節(jié)律并不同步，導(dǎo)致水稻散射的季候規(guī)律極為復(fù)雜，比如單季中稻區(qū)（如川渝貴地區(qū)）通常在5 月移栽、7 月抽穗、9 月收割，輪作油菜，間作玉米，而雙季稻區(qū)（如粵桂地區(qū)）的早稻一般在4月移栽、6月抽穗、7月收割，晚稻一般在8月移栽、9 月抽穗、11 月收割，輪作間作各類蔬菜。不同于光學(xué)遙感植被指數(shù)典型的時序曲線峰型特征，SAR 數(shù)據(jù)對水稻生育過程的響應(yīng)往往局限于移栽期至抽穗期，且極易受田間管理因素和成像噪聲干擾，難以形成穩(wěn)定的物候動態(tài)響應(yīng)。

對于大區(qū)域水稻種植遙感監(jiān)測而言，先驗信息的獲取是極為復(fù)雜繁重的工作，水稻多樣化的耕作實踐也導(dǎo)致先驗信息的代表性和可靠性難以保證，因此現(xiàn)有研究成果的適用范圍一般局限于樣本訓(xùn)練區(qū)，水稻監(jiān)測的精細程度不夠理想。

5.2 前景與展望

數(shù)據(jù)與方法的進步，為水稻種植面積提取、生理參數(shù)反演、物候與熟制識別研究提供了具有潛力的探索方向：

（1）依賴更少先驗信息的種植面積精細識別。隨著SAR 數(shù)據(jù)空間分辨率的提高，面向?qū)ο蠓诸愃惴ㄔ谝种瞥上裨肼暫头诸愒肼暦矫鎯?yōu)勢凸顯，高空間分辨率SAR 影像和面向?qū)ο蠓诸愃惴ǖ慕Y(jié)合，是破碎稻田精細提取的有效方案。在分析水稻生長機理和散射變化規(guī)律的基礎(chǔ)上，生成多種典型時序曲線，進行包容度可調(diào)節(jié)的模式識別，有望降低分類過程對復(fù)雜類型稻田樣本的依賴。

（2）顧及建模效率與精度的生理參數(shù)動態(tài)反演。兼顧機理描述和計算簡化的半經(jīng)驗?zāi)Ｐ?，目前看來是平衡反演精度和效率的不錯選擇，此外模型能夠根據(jù)實時觀測進行參數(shù)自適應(yīng)調(diào)整也很重要。數(shù)據(jù)同化方法是耦合水稻生長模型和散射模型的動態(tài)反演方法，可以做到模型實時修正和長勢預(yù)測，隨著SAR 數(shù)據(jù)時間分辨率不斷提高，該方法具有極大的應(yīng)用潛力，未來對水稻生長模型的參數(shù)簡化將為數(shù)據(jù)同化方法的區(qū)域應(yīng)用提供便利。

（3）結(jié)合生長機理和時序觀測的物候熟制自動識別。隨著SAR 歷史數(shù)據(jù)的積累，耕地種植狀態(tài)的年際監(jiān)測成為可能，未來的研究應(yīng)結(jié)合水稻的生長機理，考慮多種典型耕作環(huán)境，進一步細化散射時序與水稻物候的對應(yīng)關(guān)系，同時優(yōu)化雷達時序降噪算法，關(guān)注復(fù)雜種植結(jié)構(gòu)下水稻散射的季候變化規(guī)律，通過時序模式識別實現(xiàn)水稻物候和熟制的自動判讀，降低因地而異的閾值依賴，提高識別算法的普適性和靈活性。

縱觀水稻雷達遙感監(jiān)測研究的發(fā)展進程，成像技術(shù)的進步一直是水稻監(jiān)測技術(shù)發(fā)展的重要驅(qū)動因子。從最初的單極化到現(xiàn)在的全極化，SAR數(shù)據(jù)實現(xiàn)了波譜信息量的提升；從早期的空間分辨率數(shù)十米到目前的數(shù)米，SAR 數(shù)據(jù)實現(xiàn)了空間信息量的提升；從早期的時間分辨率數(shù)十天到目前的數(shù)天，SAR 數(shù)據(jù)實現(xiàn)了時間信息量的提升。新一代的星載SAR 衛(wèi)星的空間分辨率能達到米級，重訪頻率也大大提高，Sentinel-1 雙星重訪周期縮短至6 d，高分三號單側(cè)視重訪周期小于3 d，有效地提高了觀測頻率和動態(tài)監(jiān)測能力，這些衛(wèi)星向用戶免費便捷地分發(fā)數(shù)據(jù)，將通過增加時間信息量有力彌補空間和波譜信息短板。

6 結(jié)論

糧食安全是經(jīng)濟發(fā)展和社會穩(wěn)定的基石，水稻種植是中國糧食生產(chǎn)的關(guān)鍵組成部分，在多云霧稻作區(qū)使用雷達數(shù)據(jù)進行水稻遙感監(jiān)測，對于及時把握水稻耕作生長信息和指導(dǎo)農(nóng)業(yè)管理決策有重要意義。水稻遙感監(jiān)測主要分為種植面積提取、生理參數(shù)反演、物候和熟制識別3個部分，其中種植面積提取關(guān)注水稻的空間分布信息，生理參數(shù)反演關(guān)注水稻的植株長勢信息，物候識別關(guān)注水稻的生育演化信息，熟制識別關(guān)注水稻的耕作輪次信息，這些信息共同構(gòu)成水稻耕種狀態(tài)的綜合描述，可為水稻的耕作管理和產(chǎn)量預(yù)估提供參考依據(jù)。

伴隨著電子技術(shù)和計算機技術(shù)的發(fā)展，SAR成像能力不斷進步，各領(lǐng)域的前沿技術(shù)方法也在與遙感技術(shù)深度融合，水稻雷達遙感監(jiān)測朝著更大的覆蓋面積、更長的時間跨度和更精細的時空感知能力推進。這一切都需要對水稻生長機制和散射機理的深入理解、對模型復(fù)雜度和靈活性的改進、對SAR 數(shù)據(jù)獲取和處理能力的提升。高質(zhì)量SAR 數(shù)據(jù)的開放獲取和云計算平臺的推廣普及，將助力水稻種植的長時序大區(qū)域精準(zhǔn)監(jiān)測。