基于Segformer網(wǎng)絡(luò)的地塊尺度作物種植結(jié)構(gòu)精細(xì)化識(shí)別與分類

摘要：防止耕地“非糧化”、穩(wěn)定糧食生產(chǎn)是中國(guó)糧食安全的基石。為實(shí)現(xiàn)地塊破碎化地區(qū)作物類型及種植結(jié)構(gòu)精細(xì)化識(shí)別和分類，本研究以江蘇省泰興市為研究區(qū)，基于高分辨率遙感影像和多尺度融合特征顯著的Segformer語(yǔ)義分割模型，實(shí)現(xiàn)地塊尺度的耕地信息精細(xì)化提取;同時(shí)結(jié)合多源遙感數(shù)據(jù)構(gòu)建主要植被類型歸一化植被指數(shù)（NDVI）時(shí)序曲線及植被生長(zhǎng)關(guān)鍵時(shí)間節(jié)點(diǎn)的光譜反射特征，開展地塊尺度的作物種植結(jié)構(gòu)分類。結(jié)果表明：基于Segformer模型的分割方法可有效識(shí)別耕地，F(xiàn)1系數(shù)達(dá)92.4%;基于主要植被類型多時(shí)相NDVI時(shí)序特征及植被生長(zhǎng)關(guān)鍵時(shí)間節(jié)點(diǎn)光譜反射特征的作物種植結(jié)構(gòu)分類方法能夠?qū)崿F(xiàn)地塊尺度的種植結(jié)構(gòu)分類，總體分類精度達(dá)82.38%。因此，本研究建立的方法可有效實(shí)現(xiàn)地塊尺度耕地信息的精細(xì)化提取及種植結(jié)構(gòu)識(shí)別和分類，為耕地保護(hù)提供技術(shù)支持。

關(guān)鍵詞：種植制度;地塊尺度；精細(xì)化識(shí)別和分類;遙感

中圖分類號(hào)：S127文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1000-4440（2024）02-0293-10

Refined identification and classification of crop planting structure at plot scale based on Segformer network

GU Yu-qing1，LI Xiao-wen2，3，CAO Wei1，WANG Ya-hua2，3，ZHOU Xin-xin3，4，ZHAO Bi1

（1.Nanjing Guotu Information Industry Co.， Ltd.， Nanjing 210036， China;2.School of Geography， Nanjing Normal University， Nanjing 210023， China;3.Key Laboratory of Virtual Geographic Environment， Nanjing Normal University， Ministry of Education， Nanjing 210023， China;4.School of Geographic and Biologic Information， Nanjing University of Posts and Telecommunications， Nanjing 210023， China）

Abstract：Preventing the “non-grain” of cultivated land and stabilizing food production are the cornerstones of China’s food security. In order to realize the fine identification and classification of crop types and planting structure in the area of land fragmentation， this study took Taixing City， Jiangsu province as the research area， and realized the fine extraction of cultivated land information at the plot scale based on the high-resolution remote sensing images and the Segformer semantic segmentation model with significant multi-spatial scale fusion features. At the same time， the normalized difference vegetation index （NDVI） time series curve of the main vegetation types and the spectral reflectance characteristics at the key time nodes of vegetation growth were constructed by combining multi-source remote sensing data， and the classification of crop planting structure at the plot scale was carried out. The results showed that the segmentation method based on Segformer model could effectively identify cultivated land， and the F1 was 92.4%. The classification method of crop planting structure based on multi-temporal NDVI time series characteristics of main vegetation types and spectral reflection characteristics at key time nodes of vegetation growth could realize the classification of planting structure at plot scale， and the overall classification accuracy was 82.38%. Therefore， the method established in this study could effectively realize the fine extraction of cultivated land information at the plot scale and the identification and classification of planting structure， and provide technical support for cultivated land protection.

Key words：planting system;plot scale;refined identification and classification;remote sensing

區(qū)域尺度作物種植結(jié)構(gòu)的精細(xì)識(shí)別和分類是農(nóng)情監(jiān)測(cè)、作物生產(chǎn)力預(yù)測(cè)、農(nóng)業(yè)政策制定的信息基礎(chǔ)，對(duì)農(nóng)業(yè)精細(xì)化管理、農(nóng)業(yè)結(jié)構(gòu)調(diào)整、國(guó)家糧食安全保障具有重要意義[1-3]。2020年，國(guó)務(wù)院辦公廳出臺(tái)《關(guān)于防止耕地“非糧化”穩(wěn)定糧食生產(chǎn)的意見》，要求優(yōu)化區(qū)域布局和生產(chǎn)結(jié)構(gòu)，科學(xué)合理利用耕地資源，保障國(guó)家糧食安全。大尺度的作物種植類型和種植制度監(jiān)測(cè)雖然有一定的參考價(jià)值，但其精度一般不高。為此，實(shí)現(xiàn)地塊尺度上的作物種植類型及種植制度的精細(xì)化識(shí)別和分類勢(shì)在必行。然而，要實(shí)現(xiàn)地塊尺度上的作物種植類別和種植制度精細(xì)化識(shí)別和分類，不僅需要高時(shí)效、大尺度的耕地提取技術(shù)，還需要地塊尺度的高精度作物種植信息[4-5]。近年來，隨著遙感與深度學(xué)習(xí)技術(shù)的飛速發(fā)展，高時(shí)效高精度的耕地監(jiān)測(cè)、作物類別和種植制度精細(xì)識(shí)別和分類已成為農(nóng)業(yè)遙感領(lǐng)域的熱點(diǎn)問題[6-9]。

光譜、時(shí)相與空間特征是作物種植結(jié)構(gòu)遙感提取的理論基礎(chǔ)[1，9]。已有的作物遙感分類研究中[10-15]，主要聚焦在省級(jí)、國(guó)家級(jí)、全球尺度上，分類尺度過大，難以實(shí)現(xiàn)地塊尺度的作物提取分類。如以MODIS衛(wèi)星影像作為數(shù)據(jù)源，利用植被指數(shù)時(shí)間序列數(shù)據(jù)并結(jié)合植被生長(zhǎng)周期獲取作物信息的方法，可以實(shí)現(xiàn)大尺度種植區(qū)的作物種植結(jié)構(gòu)識(shí)別與分類，但這種方法只適用于規(guī)?；a(chǎn)地區(qū)[10-12]。對(duì)于地塊破碎化嚴(yán)重的長(zhǎng)江中下游地區(qū)，這種方法精度較低。高分辨率遙感影像可以提供地塊尺度作物類型的空間分布與紋理信息[16-22]，為地塊尺度的作物高精度提取與分類提供了可能。但高分辨率遙感衛(wèi)星的幅寬較窄與重返周期較長(zhǎng)，增加了作物識(shí)別與精細(xì)分類的難度。

鑒于遙感影像時(shí)空分辨率的矛盾性，本研究充分發(fā)揮中低分辨率遙感影像高時(shí)間分辨率與高分辨率影像的高空間分辨率優(yōu)勢(shì)，結(jié)合遙感與深度學(xué)習(xí)技術(shù)，提出一種基于深度學(xué)習(xí)語(yǔ)義分割的多源遙感作物類型及種植制度精細(xì)識(shí)別和分類方法。即通過Segformer深度學(xué)習(xí)網(wǎng)絡(luò)提取地塊尺度的耕地信息，以地塊尺度為觀測(cè)單元，結(jié)合主要植被類型周年多時(shí)相歸一化植被指數(shù)（NDVI）特征和關(guān)鍵時(shí)間節(jié)點(diǎn)的光譜反射特征，實(shí)現(xiàn)地塊尺度作物類型及種植制度的識(shí)別與精細(xì)分類和制圖。

1材料與方法

1.1研究區(qū)概況

選取江蘇省泰興市作為研究區(qū)。泰興市位于江蘇省中部、長(zhǎng)江下游北岸，介于119°54′E～120°22′E、31°58′N～32°23′N，面積1 172.27 km2，是典型的平原地貌區(qū)（圖1）。泰興市地處北亞熱帶海洋性季風(fēng)氣候區(qū)，夏季高溫多雨，年平均氣溫14.9 ℃，年平均降水量1 027 mm，年平均日照時(shí)數(shù)2 125 h，無霜期220 d。泰興市河網(wǎng)密集、地塊破碎度高、作物類別多、種植結(jié)構(gòu)復(fù)雜，主要作物包括小麥、水稻、油菜、玉米、大豆等。

1.2多源遙感數(shù)據(jù)

本研究使用覆蓋泰興市全域的高分辨率遙感影像與多時(shí)相Sentinel-2 A/B衛(wèi)星遙感影像作為數(shù)據(jù)源。高分辨率遙感影像為Google Earth的18級(jí)產(chǎn)品，影像分辨率為1.2 m。利用2020年6月泰興市全市遙感影像作為樣本來源，進(jìn)行耕地樣本的勾繪工作，用于耕地模型的訓(xùn)練。利用2022年8月泰興市1.2 m影像作為推理影像，運(yùn)用模型對(duì)耕地進(jìn)行推理，實(shí)現(xiàn)泰興市耕地的精細(xì)化提取。

多源數(shù)據(jù)結(jié)合使用具有更高的時(shí)間分辨率，有利于作物類型和種植制度的精確分類。本研究采集2022年1月至2023年4月Sentinel-2 A/B衛(wèi)星中分辨率遙感數(shù)據(jù)作為泰興市農(nóng)田作物類型和種植結(jié)構(gòu)識(shí)別和分類數(shù)據(jù)源。 Sentinel-2數(shù)據(jù)共13個(gè)光譜波段，影像分辨率包含10 m、20 m和60 m，選取其中10 m分辨率的可見光（B2、B3、B4）和近紅外（B8）波段圖像用于計(jì)算歸一化植被指數(shù)（NDVI）。Sentinel-2 A/B數(shù)據(jù)下載于歐洲空間局官網(wǎng)，級(jí)別為L(zhǎng)evel-2A級(jí)，原始地表反射率數(shù)據(jù)經(jīng)過輻射定標(biāo)、幾何校正和大氣校正后使用與分析。

光學(xué)遙感影像通常易受云霧遮擋的影響，這對(duì)耕地作物類型和種植制度的識(shí)別分類產(chǎn)生了很大干擾。因此研究中首先根據(jù)泰興市行政邊界文件裁剪得到泰興市2022年1月-2023年4月所有Sentinel-2遙感影像，然后選擇云霧覆蓋率≤5%的影像作為分析數(shù)據(jù)。最終共選取31景Sentinel-2影像，影像詳細(xì)信息如表 1所示。

1.3實(shí)地采樣數(shù)據(jù)

為掌握泰興市種植的主要作物類型以及空間分布，2023年6月1日-6月10日對(duì)泰興市耕地地塊進(jìn)行了隨機(jī)采樣。為避免因土地流轉(zhuǎn)以及采樣數(shù)據(jù)與遙感影像數(shù)據(jù)時(shí)間滯后等因素造成的樣本地塊屬性變化，選取的采樣地塊均位于永久基本農(nóng)田內(nèi)，且盡量滿足均勻分布。調(diào)查過程中，準(zhǔn)確記錄采樣地塊的經(jīng)緯度坐標(biāo)，并向當(dāng)?shù)卣胺N植戶咨詢地塊全年耕作情況以及種植作物類型等。研究中共調(diào)查982個(gè)地塊樣本，其中種植小麥-玉米的樣本183個(gè)，小麥-水稻的樣本234個(gè)、油菜-水稻的樣本163個(gè)、鮮食大豆-水稻的樣本52個(gè)、油菜-休耕樣本61個(gè)，水稻-休耕樣本117個(gè)，其他（蔬菜等）地塊169個(gè)，撂荒地塊3個(gè)。

1.4研究方法

本研究主要包括耕地地塊信息精細(xì)化提取和地塊尺度作物種植結(jié)構(gòu)識(shí)別與分類兩個(gè)內(nèi)容，技術(shù)路線如圖2所示。首先基于高分辨率遙感影像與泰興當(dāng)?shù)厥噶空{(diào)查數(shù)據(jù)，構(gòu)建泰興市耕地與非耕地光譜樣本庫(kù);基于Segformer深度學(xué)習(xí)網(wǎng)絡(luò)模型訓(xùn)練耕地樣本，建立耕地智能解譯模型，實(shí)現(xiàn)耕地地塊的精細(xì)化識(shí)別。在上述得到的耕地地塊信息的基礎(chǔ)上，基于多時(shí)序中分辨率遙感影像，構(gòu)建地塊多時(shí)相光譜物特征及不同作物生長(zhǎng)關(guān)鍵時(shí)間節(jié)點(diǎn)光譜曲線，套合語(yǔ)義分割的耕地圖斑，確定耕地地塊種植作物，實(shí)現(xiàn)地塊尺度上的作物種植結(jié)構(gòu)精細(xì)分類。

1.4.1耕地地塊信息精細(xì)化提取為提高耕地地塊識(shí)別的準(zhǔn)確率，基于泰興市2020年高分辨率遙感影像，繪制耕地與背景兩類語(yǔ)義分割樣本（背景指除耕地外的其他用地），構(gòu)建針對(duì)耕地對(duì)象的譜系化樣本庫(kù)，共計(jì)完成人工標(biāo)注耕地樣本14 496張，樣本尺寸均為512×512像素，并將其按6∶2∶2的比例隨機(jī)分為訓(xùn)練集（8 698張）、驗(yàn)證集（2 899張）與測(cè)試集（2 899張），樣本制作完成后按公共數(shù)據(jù)集PASCAL VOC2012格式存儲(chǔ)。

考慮到泰興市耕地存在破碎度高、空間分布與幾何形狀復(fù)雜、影像分辨率不一致等特征，本研究選取Segformer深度學(xué)習(xí)模型提取耕地地塊信息。該模型結(jié)合了Transformer結(jié)構(gòu)與多層感知器（MLP）解碼器。在編碼器階段，每一層Transformer結(jié)構(gòu)輸出的特征尺度逐層遞減，確保能夠獲取不同尺度的特征信息，可以有效減少幾何形狀、紋理、影像分辨率等因素對(duì)耕地信息解譯結(jié)果的影響。在獲取不同尺度耕地地塊特征信息的同時(shí)，Segformer網(wǎng)絡(luò)在位置編碼部分，擯棄了ViT網(wǎng)絡(luò)中使用的固定分辨率，彌補(bǔ)卷積神經(jīng)網(wǎng)絡(luò)（CNN）在不同分辨率影像上特征信息提取疏漏的缺陷，其網(wǎng)絡(luò)結(jié)構(gòu)圖如圖3所示。高分辨率遙感圖像耕地地塊語(yǔ)義分割的具體過程如下：

（1） 輸入耕地地塊樣本圖片，通過數(shù)據(jù)處理將耕地地塊圖片分成4×4的圖像塊（patches）;

（2） 將4×4的patches輸入到編碼器（Encoder）模塊中，生成多層級(jí)特征圖，然后基于圖像塊疊加方式合并特征圖，降低特征圖的分辨率，進(jìn)一步生成分辨率為原先1/2、1/4、1/8、1/16的特征圖，并在解碼器（Decoder）模塊中融合不同分辨率特征圖，形成耕地地塊的多層級(jí)特征。

（3）將耕地地塊多層級(jí)特征輸入MLP層，統(tǒng)一通道維度，同時(shí)在不同分辨率的特征圖上采集原圖的1/2，連接上采樣特征圖并利用MLP聚合，輸出耕地地塊目標(biāo)。

Segformer模型訓(xùn)練環(huán)境及參數(shù)配置如下：訓(xùn)練環(huán)境為Pytorch1.10，CUDA版本11.3，顯卡為NVIDIA-A40-48Q，顯存48 GB，內(nèi)存128 G;使用AdamW優(yōu)化器，初始學(xué)習(xí)率設(shè)置為1×10-6，采用指數(shù)變換策略動(dòng)態(tài)調(diào)整學(xué)習(xí)率，共訓(xùn)練300代，批大小設(shè)置為16，其他參數(shù)采用默認(rèn)值。

基于測(cè)試數(shù)據(jù)集，選取F1系數(shù)（F1）來評(píng)價(jià)耕地地塊信息提取的精度;采用總體精度（Accuracy）衡量作物分類識(shí)別精度，相關(guān)計(jì)算公式如下：

F1=2×Precision×RecallPrecision+Recall（1）

Precision=TPTP+FP（2）

Recall=TPTP+FN（3）

Accuracy=TP+TNTP+FP+TN+FN（4）

式中，TP、TN、FN、FP分別為混淆矩陣中真正類、真負(fù)類、假負(fù)類和假正類。

1.4.2作物種植模式識(shí)別與分類首先利用Sentinel-2遙感數(shù)據(jù)構(gòu)建泰興市全域歸一化植被指數(shù)（NDVI）時(shí)序庫(kù)。NDVI計(jì)算如下：

NDVI=RNIR-RREDRNIR+RRED（5）

其中，RNIR與RRED分別表示近紅外波段、紅波段的地表反射率。

為消除NDVI異常值對(duì)種植制度識(shí)別和分類的影響，本研究設(shè)計(jì)了一套NDVI異常值的剔除方案：以任意NDVI觀測(cè)時(shí)間節(jié)點(diǎn)為中心建立一個(gè)滑動(dòng)窗口，對(duì)其前時(shí)相和后時(shí)相的NDVI作線性內(nèi)插，并設(shè)置一個(gè)差異閾值（β，β取值為0.2）;若內(nèi)插值與該時(shí)間節(jié)點(diǎn)觀測(cè)值相差超過差異閾值（β），則判定該時(shí)間節(jié)點(diǎn)的觀測(cè)值為異常值，并用內(nèi)插值代替觀測(cè)值，重建NDVI時(shí)序曲線，否則保留原觀測(cè)值。

然后根據(jù)調(diào)查樣地?cái)?shù)據(jù)，建立不同種植制度下植被周年NDVI多時(shí)相特征曲線及關(guān)鍵發(fā)育期不同波段光譜特征曲線。通過NDVI時(shí)序特征曲線，確定作物種植類型，對(duì)于生長(zhǎng)周期相似的作物，以關(guān)鍵生育期的光譜曲線作為輔助信息進(jìn)一步判斷，確定泰興市的種植結(jié)構(gòu)類型。最終以上述特征曲線為依據(jù)，套合耕地地塊精細(xì)化圖斑，利用決策樹分類法對(duì)泰興市所有像素點(diǎn)的種植結(jié)構(gòu)進(jìn)行識(shí)別和分類，并通過實(shí)地調(diào)研數(shù)據(jù)對(duì)種植結(jié)構(gòu)分類識(shí)別結(jié)果進(jìn)行精度評(píng)價(jià)，進(jìn)而實(shí)現(xiàn)泰興市種植制度的精細(xì)化識(shí)別與分類。

2結(jié)果與分析

2.1耕地地塊信息提取結(jié)果

Segformer模型經(jīng)過170代訓(xùn)練后，模型達(dá)到最優(yōu)水平。利用最優(yōu)模型對(duì)測(cè)試集進(jìn)行驗(yàn)證，結(jié)果如表2所示。最優(yōu)模型的F1達(dá)到92.4%。

耕地地塊信息的預(yù)測(cè)結(jié)果如圖4所示。從圖中可以看出，模型能有效區(qū)分耕地地塊與背景地塊（水系、林地、建設(shè)用地等其他地物地塊），兩者邊界清晰，整體提取效果較好。

2.2不同植被NDVI時(shí)序特征及關(guān)鍵生長(zhǎng)時(shí)間節(jié)點(diǎn)光譜分析

982個(gè)實(shí)地調(diào)研地塊，共篩選得到8種典型的種植制度：小麥-玉米、油菜-休耕、小麥-水稻、水稻-休耕、油菜-水稻、鮮食大豆-水稻、其他、撂荒。不同種植制度下，耕地周年NDVI時(shí)序特征曲線如圖5所示。從圖中可以看出，不同種植制度下由于作物種植類型及季節(jié)的差異，耕地周年NDVI時(shí)序特征存在明顯的差異。小麥與油菜是主要的越冬作物，一般在10-11月播種或移栽，此時(shí)NDVI較低，此后隨著小麥、油菜的生長(zhǎng)，NDVI逐漸增加，至4月達(dá)峰值。此后隨著小麥、油菜的成熟和收獲，NDVI呈現(xiàn)下降特征（圖5A、圖5B、圖5C、圖5E）。夏季種植的作物主要為水稻，其NDVI曲線一般從6月開始逐漸增加，9月前后達(dá)峰值，此后隨著水稻的成熟，NDVI逐漸減少（圖5C～圖5F）。玉米的生長(zhǎng)季與水稻類似，但NDVI峰值出現(xiàn)的時(shí)間遲、持續(xù)時(shí)間短（圖5A）。作為水稻的前茬，鮮食大豆NDVI從5月中下旬開始逐漸增加，6月中下旬達(dá)峰值，此后隨著鮮食大豆的收獲，NDVI迅速下降（圖5F）。其他植被（如蔬菜、果園等）其NDVI曲線波動(dòng)無明顯的周期性規(guī)律（圖5G）。而休耕季節(jié)及撂荒耕地，NDVI一般始終保持較低狀態(tài)（圖5B、圖5D、圖5H）。不同植被生長(zhǎng)關(guān)鍵時(shí)間節(jié)點(diǎn)，Sentinel-2多光譜數(shù)據(jù)不同波段的反射率曲線如圖6所示。從圖6中可以看出，同一生長(zhǎng)季不同作物生長(zhǎng)關(guān)鍵時(shí)間節(jié)點(diǎn)不同波段反射率曲線存在一定的差異，可以作為種植制度識(shí)別的輔助信息。2022年4月8日，小麥、油菜生長(zhǎng)旺盛，植被特征明顯。小麥在近紅外波段（B8）的光譜反射峰明顯高于油菜，而其他波段的反射率則是油菜高于小麥。2022年6月17日，此時(shí)鮮食大豆生長(zhǎng)旺盛、植被特征明顯，而休耕地?zé)o明顯的植被特征。2022年10月10日，玉米、水稻植被特征明顯，水稻光譜在近紅外波段（B8）的反射率明顯高于玉米，在紅波段（B4）的吸收率低于玉米。除上述常見種植制度外，還有部分蔬菜、果園、草地等（本文統(tǒng)稱為其他類），其輪作特征不明顯。

根據(jù)不同種植制度下的耕地周年NDVI時(shí)序特征曲線和作物生長(zhǎng)關(guān)鍵時(shí)間節(jié)點(diǎn)Sentinel-2多光譜數(shù)據(jù)不同波段的反射率曲線對(duì)982個(gè)樣地的種植制度進(jìn)行檢驗(yàn)，結(jié)果如表3所示。從種植制度分類的混淆矩陣可以得到，泰興市種植結(jié)構(gòu)分類的總體精度為82.38%，這說明泰興市種植制度的識(shí)別和分類較為準(zhǔn)確。

2.3泰興市地塊尺度耕地種植結(jié)構(gòu)分類

泰興市地塊尺度的耕地種植結(jié)構(gòu)如圖7所示。可以看出，泰興市主要的種植結(jié)構(gòu)為小麥-水稻輪作，面積約為371 km2，約占總耕地面積的32.8%，其他種植結(jié)構(gòu)分布較為零散。部分地塊，存在休耕，休耕總面積約181 km2，占總耕地面積的16.0%。撂荒地塊面積占比不到1.0%。根據(jù)泰興市作物種植結(jié)構(gòu)分類結(jié)構(gòu)圖，可以看出，泰興市存在少量耕地“非糧化”傾向，呈零散式分布，未出現(xiàn)大規(guī)?！胺羌Z化”行為，種植結(jié)構(gòu)方式較為穩(wěn)定。

3結(jié)論

利用單一遙感數(shù)據(jù)源進(jìn)行區(qū)域種植結(jié)構(gòu)識(shí)別難度大、精度低，本研究提出了一種基于深度學(xué)習(xí)和多源遙感數(shù)據(jù)的地塊破碎化嚴(yán)重地區(qū)的種植結(jié)構(gòu)分類方法，得出以下結(jié)論：

（1）利用Segformer深度學(xué)習(xí)語(yǔ)義分割模型和高空間分辨率遙感數(shù)據(jù)，構(gòu)建泰興市耕地樣本，可以有效識(shí)別耕地圖斑，實(shí)現(xiàn)地塊尺度耕地信息的精細(xì)提取，耕地與非耕地邊界清晰，驗(yàn)證集分類的F1達(dá)92.4%。

（2）在耕地信息提取的基礎(chǔ)上，進(jìn)一步利用中分辨率多源遙感數(shù)據(jù)，構(gòu)建耕地不同種植制度周年NDVI時(shí)序特征曲線及典型植被生長(zhǎng)關(guān)鍵時(shí)間節(jié)點(diǎn)的光譜特征曲線，實(shí)現(xiàn)了泰興市種植制度的精細(xì)化分類，調(diào)查樣地的分類總體精度達(dá)82.38%。泰興市主要的種植制度為小麥-水稻輪作，約占耕地總面積的32.8%。油菜是主要油料作物，泰興市并未出現(xiàn)大規(guī)?！胺羌Z化”行為。

參考文獻(xiàn)：

[1]胡瓊，吳文斌，宋茜，等. 農(nóng)作物種植結(jié)構(gòu)遙感提取研究進(jìn)展[J]. 中國(guó)農(nóng)業(yè)科學(xué)，2015，48（10）：1900-1914.

[2]SALMON J M， FRIEDL M A， FROLKING S， et al. Global rain-fed，irrigated， and paddy croplands：a new high resolution mapderived from remote sensing， crop inventories and climate data[J]. International Journal of Applied Earth Observation amp; Geoinformation，2015，38：321-334.

[3]陳仲新，任建強(qiáng)，唐華俊，等. 農(nóng)業(yè)遙感研究應(yīng)用進(jìn)展與展望[J]. 遙感學(xué)報(bào)，2016，20（5）：748-767.

[4]劉煥軍，邱政超，孟令華，等. 黑土區(qū)田塊尺度遙感精準(zhǔn)管理分區(qū)[J]. 遙感學(xué)報(bào)，2017，21（3）：470-478.

[5]劉哲，劉帝佑，朱德海，等. 作物遙感精細(xì)識(shí)別與自動(dòng)制圖研究進(jìn)展與展望[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2018，49（12）：1-12.

[6]WU M， HUANG W， NIU Z， et al. Fine crop mapping by combining high spectral and high spatial resolution remote sensing data in complex heterogeneous areas[J]. Computers and Electronics in Agriculture，2017，139：1-9.

[7]YU B， SHANG S. Multi-year mapping of maize and sunflower in Hetao irrigation district of China with high spatial and temporal resolution vegetation index series[J]. Remote Sensing，2017，9（8）：855.

[8]楊穎頻，吳志峰，駱劍承，等. 時(shí)空協(xié)同的地塊尺度作物分布遙感提取[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2021，37（7）：166-174.

[9]張鵬，胡守庚. 地塊尺度的復(fù)雜種植區(qū)作物遙感精細(xì)分類[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2019，35（20）：125-134.

[10]ZHENG Y， WU B， ZHANG M， et al. Crop phenology detection using high spatio-temporal resolution data fused from SPOT5 and MODIS products[J]. Sensors，2016，16（12）：2099.

[11]GAO F， ANDERSON M C， ZHANG X， et al. Toward mapping crop progress at field scales through fusion of Landsat and MODIS imagery[J]. Remote Sensing of Environment，2017，188：9-25.

[12]HAO P， ZHAN Y， WANG L， et al. Feature selection of time series MODIS data for early crop classification using random forest：a case study in Kansas， USA[J]. Remote Sensing，2015，7（5）：5347-5369.

[13]王娜，李強(qiáng)子，杜鑫，等. 單變量特征選擇的蘇北地區(qū)主要農(nóng)作物遙感識(shí)別[J]. 遙感學(xué)報(bào)，2017，21（4）：519-530.

[14]SONOBE R， TANI H， WANG X， et al. Parameter tuning in the support vector machine and random forest and their performances in cross-and same-year crop classification using TerraSAR-X[J]. International Journal of Remote Sensing，2014，35（23）：7898-7909.

[15]PELLETIER C， VALERO S， INGLADA J， et al. Assessing the robustness of Random Forests to map land cover with high resolution satellite image time series over large areas[J]. Remote Sensing of Environment，2016，187：156-168.

[16]DONG J， XIAO X， MENARGUEZ M A， et al. Mapping paddy rice planting area in northeastern Asia with Landsat 8 images， phenology-based algorithm and Google Earth Engine[J]. Remote Sensing of Environment，2016，185：142-154.

[17]鄧劉洋，沈占鋒，柯映明，等. 基于地塊尺度多時(shí)相遙感影像的冬小麥種植面積提取[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2018，34（21）：157-164.

[18]田海峰，鄔明權(quán)，牛錚，等. 基于Radarsat-2 影像的復(fù)雜種植結(jié)構(gòu)下旱地作物識(shí)別[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2015，31（23）：154-159.

[19]黃啟廳，覃澤林，曾志康. 多星數(shù)據(jù)協(xié)同的地塊尺度作物分類與面積估算方法研究[J]. 地球信息科學(xué)學(xué)報(bào)，2016，18（5）：708-717.

[20]韓衍欣，蒙繼華. 面向地塊的農(nóng)作物遙感分類研究進(jìn)展[J]. 國(guó)土資源遙感，2019，31（2）：1-9.

[21]肖艷，姜琦剛，王斌，等. 基于ReliefF 和PSO 混合特征選擇的面向?qū)ο笸恋乩梅诸怺J]. 農(nóng)業(yè)工程學(xué)報(bào)，2016，32（4）：211-216.

[22]XIE E， WANG W， YU Z， et al. SegFormer：simple and efficient design for semantic segmentation with transformers[J]. Advances in Neural Information Processing Systems，2021，34：12077-12090.

（責(zé)任編輯：石春林）