盧輝雄， 田言亮， 董雙發(fā)， 汪 冰， 李名松， 牛海威， 張 恩， 薛 慶， 章新益， 程思思

(1.核工業(yè)航測(cè)遙感中心，石家莊 050002；2.東華理工大學(xué) 地球科學(xué)學(xué)院，南昌 330013； 3.中國(guó)地質(zhì)科學(xué)院水文地質(zhì)環(huán)境地質(zhì)研究所，石家莊 050061)

基于Landsat 8 OLI數(shù)據(jù)的華北平原農(nóng)作物種植結(jié)構(gòu)遙感解譯研究

盧輝雄1,2， 田言亮3， 董雙發(fā)1， 汪 冰1， 李名松1， 牛海威1， 張 恩1,2， 薛 慶1,2， 章新益1， 程思思1

(1.核工業(yè)航測(cè)遙感中心，石家莊 050002；2.東華理工大學(xué) 地球科學(xué)學(xué)院，南昌 330013； 3.中國(guó)地質(zhì)科學(xué)院水文地質(zhì)環(huán)境地質(zhì)研究所，石家莊 050061)

利用華北平原2013年65景Landsat 8 OLI數(shù)據(jù)，通過(guò)ENVI軟件對(duì)影像進(jìn)行預(yù)處理，采用面向?qū)ο蟮男畔⑻崛『腿藱C(jī)交互式解譯方法，對(duì)華北平原小麥、玉米、水稻、大豆等主要作物分布及種植結(jié)構(gòu)進(jìn)行了系統(tǒng)的解譯，劃分了研究區(qū)小麥玉米兩熟區(qū)、小麥水稻兩熟區(qū)、水稻兩熟區(qū)、小麥大豆兩熟區(qū)、玉米一熟區(qū)、水稻一熟區(qū)分布范圍，并分析了研究區(qū)農(nóng)作物種植結(jié)構(gòu)影響因素。結(jié)果表明，基于Landsat 8 OLI數(shù)據(jù)對(duì)大區(qū)域主要作物種植面積及種植結(jié)構(gòu)遙感解譯的技術(shù)和方法是可行的，能為研究農(nóng)作物種植結(jié)構(gòu)提供有效方法。

種植結(jié)構(gòu)； 遙感解譯； Landsat8 OLI； 華北平原

0 前言

在當(dāng)前市場(chǎng)經(jīng)濟(jì)及技術(shù)條件下，利用遙感技術(shù)可以在短期內(nèi)連續(xù)獲取大范圍的地面信息，能夠?qū)崿F(xiàn)對(duì)糧食作物種植面積快速、客觀監(jiān)測(cè)。但由于不同作物物候歷不同，且不同作物具有交錯(cuò)、插花種植等現(xiàn)象，對(duì)于大面積不同作物種植結(jié)構(gòu)區(qū)，如何利用快速、經(jīng)濟(jì)可行的遙感解譯技術(shù)準(zhǔn)確獲取不同作物面積已成為熱點(diǎn)和難點(diǎn)[1-5]。

目前,隨著傳感器空間分辨率地提高，利用SPOT、TM等數(shù)據(jù)提取作物播種面積的研究開(kāi)展較多，但這些數(shù)據(jù)因價(jià)格、分辨率等因素，多用于小空間尺度地研究[6-8]；適用于大尺度范圍農(nóng)作物信息提取研究的 NOAA AVHRR、MODIS等具有低空間分辨率、高時(shí)間分辨率的遙感影像又很難保證提取結(jié)果的準(zhǔn)確性[9-15]。高分辨率的遙感影像可以精確地提取不同作物信息，但是其覆蓋面積小，遙感解譯的工作量大，費(fèi)時(shí)費(fèi)力，并且重訪周期長(zhǎng)，不能及時(shí)監(jiān)測(cè)；低分辨率遙感影像具有覆蓋面積大，重訪周期短的優(yōu)點(diǎn)，但是精度低，容易導(dǎo)致提取信息的不準(zhǔn)確。

利用Landsat8 數(shù)據(jù)進(jìn)行遙感解譯研究，主要應(yīng)用在耕地信息提取或單一作物的大面積提取中[16]，針對(duì)大尺度多類(lèi)型的作物解譯應(yīng)用較少。筆者在前人的研究的基礎(chǔ)上，針對(duì)華北平原研究區(qū)，利用Landsat8 OLI數(shù)據(jù)，使用面向?qū)ο蟮男畔⑻崛『腿藱C(jī)交互式解譯相結(jié)合的方法，對(duì)華北平原不同類(lèi)型作物的面積及其種植結(jié)構(gòu)進(jìn)行了系統(tǒng)解譯，旨在探討Landsat8 OLI數(shù)據(jù)在大區(qū)域作物解譯中的應(yīng)用。

1 研究區(qū)概況

研究區(qū)位于華北平原平原區(qū)，包括北京、天津、河北、河南、江蘇、山東、山西、安徽8個(gè)省，79個(gè)地級(jí)市，位于113°00′E～121°30′E，32°00N′～40°30′N(xiāo)之間，總面積約32×104km2。

研究區(qū)三面環(huán)山呈簸箕狀傾向渤海、黃海，平原區(qū)地勢(shì)平坦，海拔大部分在 50 m以下，從南到北貫穿有淮海、黃河、海河三大水系；屬暖溫帶半濕潤(rùn)、濕潤(rùn)季風(fēng)氣候，年平均降水量400 mm～1 200 mm，具有從南到北逐漸減少趨勢(shì)。研究區(qū)農(nóng)作物主要以冬小麥、玉米為主，其次為大豆、水稻、棉花、油料等，是“北糧南調(diào)”的重要的農(nóng)產(chǎn)品基地，也是保障我國(guó)未來(lái)食物安全的主要區(qū)域。

研究區(qū)作物主要為一年兩熟、其次為一年一熟。糧食作物主要以冬小麥、玉米、水稻、大豆為主，各類(lèi)作物的生長(zhǎng)期及解譯時(shí)相見(jiàn)表1。

表1 研究區(qū)內(nèi)各類(lèi)作物的生長(zhǎng)期

2 數(shù)據(jù)與方法

2.1 數(shù)據(jù)源與數(shù)據(jù)預(yù)處理

根據(jù)研究區(qū)作物的生長(zhǎng)期及最佳解譯時(shí)相，同時(shí)考慮到研究區(qū)面積較大，本次研究選取了兩期Landsat8影像，共65景數(shù)據(jù)。其中第一期影像時(shí)相的選擇為夏季農(nóng)作物最佳解譯時(shí)相，成像時(shí)間主要集中在2013年3月-5月；第二期影像的選擇為秋季農(nóng)作物最佳解譯時(shí)相，成像時(shí)間主要集中在2013年6月-9月。通過(guò)兩期時(shí)相的綜合對(duì)比，能夠有效確定各類(lèi)作物的解譯特征。

數(shù)據(jù)預(yù)處理主要包括輻射校正、大氣校正、影像融合、彩色合成處理、信息增強(qiáng)處理等。影像融合經(jīng)過(guò)多次對(duì)比試驗(yàn)，確定Gram-Schmidt融合處理方法，通過(guò)對(duì)遙感全色數(shù)據(jù)與多光譜數(shù)據(jù)進(jìn)行融合，形成兼具高空間分辨率空間信息和多光譜彩色信息的融合影像。在此基礎(chǔ)上，利用7(短波紅外)、5(近紅外)、4(紅)波段進(jìn)行了彩色合成處理。同時(shí)，為了突出主要糧食作物(小麥、玉米、水稻等)信息，運(yùn)用了線性及非線性拉伸，對(duì)比度拉伸、飽和度拉伸等方法對(duì)影像進(jìn)行了增強(qiáng)處理，最后獲取了空間分辨率為15 m的影像。

2.2 解譯方法

2.2.1 面向?qū)ο蟮男畔⑻崛?/p>

面向?qū)ο蟮倪b感影像分類(lèi)技術(shù)是以像元為基本單元進(jìn)行遙感信息提取，采用決策支持的模糊分類(lèi)方法，建立不同尺度的分類(lèi)層次，在每一層次上分別定義對(duì)象的光譜特征、形狀特征、紋理特征、上下文關(guān)系特征和相鄰關(guān)系特征，通過(guò)對(duì)影像對(duì)象定義多種特征并指定模糊化函數(shù)，給出每個(gè)對(duì)象隸屬于某一類(lèi)的概率，建立分類(lèi)標(biāo)準(zhǔn)，最終按照最大概率產(chǎn)生確定分類(lèi)結(jié)果。采用面向?qū)ο蠓椒ㄌ崛≈饕魑锏牧鞒倘鐖D1所示。

圖1 面向?qū)ο蛐畔⑻崛×鞒虉DFig.1 Oriented to the flow chart of information extraction

1)影像分割。影像分割分兩個(gè)步驟：①對(duì)整個(gè)影像進(jìn)行尺度空間的構(gòu)建，即根據(jù)所預(yù)設(shè)的尺度分割參數(shù)，將影像中的像元?jiǎng)澐譃椴煌炼戎档膮^(qū)域；②根據(jù)其他的分割參數(shù)提取出不同的對(duì)象多邊形。本次采用多尺度分割方法，根據(jù)研究區(qū)影像紋理特征與農(nóng)作物實(shí)際生長(zhǎng)特征，采用不同尺度進(jìn)行分割，根據(jù)分割結(jié)果對(duì)比分析( 圖2)，最終確立的研究區(qū)的最優(yōu)分割尺度(Scale:60；Color:0.9)。

圖2 多尺度分割效果對(duì)比圖Fig.2 Multiscale segmentation comparison chart(a)Scale:60;Color:0.6；(b)Scale:60;Color:0.9； (c)Scale:100;Color:0.6；(d)Scale:100;Color:0.9

2)信息提取。信息提取，即基于分割結(jié)果，對(duì)某一類(lèi)的典型特征或條件進(jìn)行描述，影像對(duì)象按照是否滿足這些屬性進(jìn)行分類(lèi)，本次采用最鄰近分類(lèi)方法進(jìn)行信息提取。

最鄰近分類(lèi)，主要包括定義特征空間、定義分類(lèi)體系、選擇典型樣本、執(zhí)行分類(lèi)四個(gè)步驟。①通過(guò)對(duì)研究區(qū)目標(biāo)地物進(jìn)行特征分析，選擇多維的典型特征作為特征空間；②接著通過(guò)目視判讀，結(jié)合其他相關(guān)資料，定義分類(lèi)體系，選擇典型的樣本對(duì)象；③在特征空間中，按照公式尋找最近的樣本對(duì)象，如果某一對(duì)象與某類(lèi)地物距離最近，則將其分為該類(lèi)。通過(guò)該方法，對(duì)研究區(qū)多尺度分割的影像進(jìn)行了信息提取，獲得了各類(lèi)作物分類(lèi)結(jié)果，分類(lèi)效果見(jiàn)圖3～圖4。

2.2.2 人機(jī)交互式解譯

利用研究區(qū)信息提取結(jié)果，根據(jù)建立的解譯標(biāo)志，進(jìn)行人機(jī)交互式解譯，對(duì)提取的每一類(lèi)圖斑賦屬性，對(duì)信息提取中局部存在錯(cuò)分、漏分的圖斑進(jìn)行完善，最終形成解譯成果。

圖3 多尺度分割效果圖Fig.3 Signal multi-scale segmentation effect

圖4 信息提取效果示意圖Fig.4 Information extraction effect diagram

3 結(jié)果及分析

3.1 研究區(qū)作物空間分布特征

研究區(qū)解譯作物主要以冬小麥、玉米為主，其次為水稻、大豆、果園、蔬菜及其他類(lèi)作物(棉花、花生等)，其空間分布見(jiàn)圖5～圖6。

3.2 研究區(qū)作物種植結(jié)構(gòu)特征

通過(guò)兩期解譯，研究區(qū)解譯作物總播面積為37.86×104km2，其中冬小麥種植面積為14.41×104km2，占解譯作物總播面積的38.06%；玉米種植面積為12.15×104km2，占解譯作物總播面積的32.10%；水稻種植面積為3.51×104km2，占解譯作物總播面積的9.27%；大豆種植面積為1.14×104km2，占解譯作物總播面積的3.01%；果園種植面積為0.44×104km2，占解譯作物總播面積的1.17%；蔬菜種植面積為3.67×104km2萬(wàn)畝，占解譯作物總播面積的9.70%；其他作物種植面積為2.53×104km2，占解譯作物總播面積的6.69%，解譯的其他類(lèi)作物主要以棉花為主。解譯的各類(lèi)作物結(jié)構(gòu)見(jiàn)圖7。

通過(guò)對(duì)研究區(qū)解譯的主要作物進(jìn)行路線效驗(yàn)、實(shí)驗(yàn)區(qū)調(diào)查工作，小麥、玉米等大面積種植類(lèi)作物解譯精度在90%以上，水稻、棉花等種植面積較小、具有插花種植特征的作物解譯進(jìn)度在70%～80%之間。研究區(qū)各種作物總體解譯精度為86.01%(表2)，對(duì)大區(qū)域作物解譯而言，精度滿足需要。

圖5 研究區(qū)夏季農(nóng)作物布局結(jié)構(gòu)解譯分布圖Fig.5 Summer crops in the study area layout structure interpretation map

圖6 研究區(qū)秋季農(nóng)作物布局結(jié)構(gòu)分布解譯圖Fig.6 Autumn crops in the study area layout structure interpretation map

序號(hào) 效驗(yàn)?zāi)康膱D斑效驗(yàn)圖斑數(shù)正確圖斑數(shù)準(zhǔn)確率/%1冬小麥12211090．22玉米13512290．33水稻4030754棉花171270．65裸地蔬菜423071．46設(shè)施蔬菜(大棚)343191．27果園393487．2合計(jì)42936986．01

圖7 研究區(qū)解譯作物結(jié)構(gòu)圖Fig.7 The structure interpretation of crops

3.3 研究區(qū)作物種植結(jié)構(gòu)劃分

根據(jù)研究區(qū)解譯結(jié)果，對(duì)研究區(qū)內(nèi)解譯的小麥、玉米、水稻、大豆進(jìn)行種植結(jié)構(gòu)劃分。按作物耕作制度研究區(qū)內(nèi)作物可劃分為小麥玉米兩熟區(qū)、小麥水稻兩熟區(qū)、水稻兩熟區(qū)、小麥大豆兩熟區(qū)、玉米一熟區(qū)、水稻一熟區(qū)(圖8)。

1)小麥玉米兩熟區(qū)，主要分布地有：①研究區(qū)河北省清苑縣至邯鄲縣、大名縣北部至吳橋；②山東省樂(lè)棱市—禹城縣—東明縣一帶；③河南省南樂(lè)縣至衛(wèi)輝市、開(kāi)封縣南部至淮陽(yáng)縣、寧棱縣至夏邑縣一帶；④安徽省利辛縣東部至宿縣南部一帶。

2)小麥水稻兩熟區(qū)，主要集中分布在3個(gè)區(qū)域：①黃河沿線，主要分布在延津縣南部至東明縣東部一帶；②淮海沿線，主要位于阜南縣南部至壽縣一帶；③主要位于江蘇省境內(nèi)，主要位于淮安市北部至東?？h南部一帶。

3)水稻兩熟區(qū)，主要分布在黃河與淮河之間，集中分布在3片區(qū)域：①淮河沿線固始縣及壽縣北東部一帶；②江蘇省境內(nèi)的宿遷市至泗洪縣一帶；③江蘇省濱?？h一帶。

圖8 研究區(qū)種植結(jié)構(gòu)區(qū)劃圖Fig.8 Planting structure zoning map in the study area

4)小麥大豆兩熟區(qū)，主要集中分布在安徽省界首市至渦陽(yáng)縣一帶。

5)玉米一熟區(qū)，主要集中在黃河以北區(qū)域，在河北省獻(xiàn)縣至文安縣一帶、河北省灤縣至遼寧省綏中縣一帶、山東省利津縣一帶分布較為廣泛。

6)水稻一熟區(qū)，主要分布在河北省唐?？h至天津市靜?？h一帶。

3.4 研究區(qū)作物種植結(jié)構(gòu)影響因素分析

研究區(qū)地形地貌、氣候、地下水資源存在區(qū)域不均衡性，對(duì)作物種植具有較大的影響。本次在解譯的基礎(chǔ)上，綜合分析了研究區(qū)作物種植結(jié)構(gòu)的影響因素。

研究區(qū)山前洪積—沖積扇形平原主要以冬小麥—玉米為主，具有一年兩熟特征；沖積平原主要為小麥玉米一年兩熟、小麥水稻兩熟，其次為玉米一年一熟、水稻一年兩熟、小麥大豆一年兩熟。海積平原具有棉花和玉米一年一熟的種植結(jié)構(gòu)(圖9)。

研究區(qū)面積較大，氣溫、降水等氣候因素在該區(qū)不同地域具有較大差異，其中降水量因子對(duì)農(nóng)作物的種植結(jié)構(gòu)影響較大?？v向上，因降水量因素影響，南北區(qū)域農(nóng)作物種植類(lèi)型不同。研究區(qū)黃河以北區(qū)域，年均降水量為400 mm～600 mm，主要為小麥玉米一年兩熟，局部存在棉花一年一熟的種植結(jié)構(gòu)(圖10)。橫向上，濱海平原區(qū)至山前沖洪積平原區(qū)，因降水量不同，作物種植存在差異。濱海平原區(qū)年降水在650 mm以上，除河北省唐?？h至天津市靜?？h以種植水稻為主，其他區(qū)域主要種植棉花為主；在研究區(qū)西部山前沖洪積平原區(qū)，降水量較少，降水量年均約為500 mm，該區(qū)域以種植小麥玉米為主。

研究區(qū)不同的地下水資源類(lèi)型區(qū)，作物種植結(jié)構(gòu)存在差異。研究區(qū)微咸水分布區(qū)域，以小麥、玉米種植為主，具有小麥玉米兩熟特征；半咸水分布區(qū)域，土壤鹽堿化較強(qiáng)烈，以種植棉花為主；南部淮河流域松散巖類(lèi)孔隙水補(bǔ)給較高區(qū)域，主要以水稻種植為主，其次為小麥，水稻具有一年兩熟特征。

4 結(jié)論

1)本次利用Landsat 8 OLI數(shù)據(jù)對(duì)華北平原小麥、玉米、水稻、大豆等主要作物進(jìn)行了解譯，通過(guò)效驗(yàn)，總體解譯精度達(dá)到86.01%，對(duì)大區(qū)域作物解譯而言，精度基本滿足需要。

圖9 研究區(qū)不同地形地貌區(qū)作物種植結(jié)構(gòu)圖Fig.9 The structure of crops planting structure in different areas of the study area

圖10 研究區(qū)不同降水量區(qū)作物種植結(jié)構(gòu)圖Fig.10 The crop planting structure in different rainfall areas in the study area

2)利用Landsat 8 OLI數(shù)據(jù)在大區(qū)域開(kāi)展遙感解譯，相比于高空間分辨率遙感數(shù)據(jù)，更經(jīng)濟(jì)、快速。

3)基于解譯結(jié)果，能夠有效劃分研究區(qū)作物種植結(jié)構(gòu)區(qū)劃圖，獲得大區(qū)域作物種植結(jié)構(gòu)特征，并通過(guò)分析影響農(nóng)作物種植結(jié)構(gòu)的主要因素，進(jìn)一步闡述解譯結(jié)果的吻合性。

4)研究表明，利用Landsat 8 OLI數(shù)據(jù)，采用面向?qū)ο笮畔⑻崛『腿藱C(jī)交互式解譯方法，能高效快速獲取大區(qū)域主要作物的種植結(jié)構(gòu)特征，該方法具有可行性。

[1] 屈寶香，張華，李剛. 中國(guó)糧食生產(chǎn)布局與結(jié)構(gòu)區(qū)域演變分析[J].中國(guó)農(nóng)業(yè)資源與區(qū)劃，2011，32(1) : 1-6. QU B X, ZHANG H, LI G. Analysis on regional development of China's grain production distribution and structure [J]. China Agricultural Resources and Regional Planning, 2011,32 (1): 1-6.(In Chinese)

[2] 李正國(guó)，唐華俊，揚(yáng)鵬. 植被物候特征的遙感提取與農(nóng)業(yè)應(yīng)用綜述[J].中國(guó)農(nóng)業(yè)資源與區(qū)劃，2012，33(5) : 20-28. LI Z G, TANG H J, YANG P. Remote sensing extraction and agricultural application of vegetation phenology [J]. China Agricultural Resources and Regional Planning, 2012, 33 (5): 20-28. (In Chinese)

[3] 李強(qiáng)子，吳炳方.作物種植成數(shù)的遙感監(jiān)測(cè)精度評(píng)價(jià)[J].遙感學(xué)報(bào)，2004，8(6) : 581-587. LI Q Z, WU B F. Precision evaluation of remote sensing monitoring of crop planting [J]. Journal of Remote Sensing, 2004, 8 (6): 581-587. (In Chinese)

[4] 吳炳方，范錦龍，田亦陳,等.全國(guó)作物種植結(jié)構(gòu)快速調(diào)查技術(shù)與應(yīng)用[J].遙感學(xué)報(bào)，2004，8(6) : 618-627. WU B F, FAN J L, TIAN Y C,et al. Rapid survey technology and application of crop planting structure in China [J]. Journal of RemoteSensing, 2004, 8 (6): 618-627. (In Chinese)

[5] 曾瀾，張勇，張雪，等.基于遙感估產(chǎn)和農(nóng)業(yè)統(tǒng)計(jì)的糧食供需平衡模型研究[J].遙感學(xué)報(bào)，2004，8(6) : 645-654. ZENG L, ZHANG Y, ZHANG X, et al. Based on yield estimation by remote sensing and the agricultural statistics of grain supply and demand equilibrium model [J]. Journal of remote sensing, 2004, 8 (6): 645-654. (In Chinese)

[6] 鄭長(zhǎng)春,王秀珍,黃敬峰.基于特征波段的SPOT-5衛(wèi)星影像水稻面積信息自動(dòng)提取的方法研究[J].遙感技術(shù)與應(yīng)用,2008,23(3):294-299. ZHENG C C, WANG X Z, HUANG J F. Research on automatic rice area information extraction based on characteristic band SPOT-5 [J]. Remote Sensing Technology and Application, 2008,23 (3): 294-299. (In Chinese)

[7] 鄧勁松，王坷.基于特征波段的spot5衛(wèi)星影像耕地信息自動(dòng)提取的方法研究[J].農(nóng)業(yè)工程學(xué)報(bào),2004，20(6):145-147. DENG J S, WANG K. Research on the automatic extraction of cultivated land information based on the SPOT5 satellite image in the feature band [J]. Journal of Agricultural Engineering,2004, 20 (6): 145-147. (In Chinese)

[8] 祝民強(qiáng)，周萬(wàn)蓬，吳仁貴，等.贛中丘陵區(qū)TM影像的耕地信息提取模型研究[J].遙感技術(shù)與應(yīng)用，2004，19(3):173-176. ZHU M Q, Z W P, WU R G,et al.TM images of the hill region in central Jiangxi province cultivated land information extraction model [J]. Remote Sensing Technology and Application, 2004, 19 (3): 173-176. (In Chinese)

[9] 林忠輝，莫興國(guó).NDVI時(shí)間序列諧波分析與地表物候信息獲取[J].農(nóng)業(yè)工程學(xué)報(bào)，2006，22(12):138-144. LIN Z H, MO X G. The harmonic analysis of the.NDVI time series and the land surface phenology information acquisition [J]. Journal of Agricultural Engineering, 2006, 22 (12): 138-144. (In Chinese)

[10]陳水森，柳欽火，陳良富，等.糧食作物播種面積遙感監(jiān)測(cè)研究進(jìn)展[J].農(nóng)業(yè)工程學(xué)報(bào)，2005，21(6):166-171. CHEN S S, LIU Q H,CHEN L F, et al. Research progress on remote sensing monitoring of sown area of grain crops [J]. Journal of Agricultural Engineering, 2005, 21 (6): 166-171. (In Chinese)

[11]張莉，吳文斌，左麗君，等.基于 EOS/MODIS 數(shù)據(jù)的南方水稻面積提取技術(shù)[J].中國(guó)農(nóng)業(yè)資源與區(qū)劃，2011，32(4):39-44. ZHANG L, WU W B, ZUO L J, et al. Rice area extraction technology based on EOS/MODIS data [J]. China Agricultural Resources and Regional Planning, 2011, 32 (4): 39-44. (In Chinese)

[12]胡宗辰.基于MODIS的中國(guó)主要糧食作物種植時(shí)空分布信息提取方法研究[D].成都:電子科技大學(xué)，2013:1-10. HU Z C. Based on the spatial temporal distribution information extraction method of main grain crops in China [D]. Chengdu:Electronic Science and Technology University, 2013:1-10. (In Chinese)

[13]呂婷婷,劉闖.基于MODIS數(shù)據(jù)的泰國(guó)耕地信息提取[J].農(nóng)業(yè)工程學(xué)報(bào),2010,26(2):244-250. LV T T, LIU C. Information extraction of cultivated land in Thailand based on MODIS data [J]. Journal of agricultural engineering, 2010,26 (2): 244-250. (In Chinese)

[14]左麗君,董婷婷,汪瀟,等.基于MODIS/EVI的中國(guó)北方耕地復(fù)種指數(shù)提取[J].農(nóng)業(yè)工程學(xué)報(bào),2009,25(8):141-146. ZUO L J,DONG T T,WANG X, et al. Transactions of the CSAE evi in northern China for multiple cropping index of arable land extraction based on MODIS/EVI [J]. Journal of Agricultural Engineering, 2009,25 (8): 141-146. (In Chinese)

[15]TOMITA A, INOUE Y, OGAWA S, et al. Vegetation patterns in the ChaoPhraya Delta, 1997 dry season using satellite image data[A].Proceedings of the international conference: The Chao PhrayaDelta: Historical development, dynamics and challenges of Thailandrice boel[C]. Thailand,2000:35-38.

[16]牛魯燕，張曉艷，鄭繼業(yè)，等.基于Landsat8 OLI數(shù)據(jù)的山東省耕地信息提取研究[J].中國(guó)農(nóng)學(xué)通報(bào),2014,30(34):264-269. NIU LUYAN,ZHANG XIAOYAN, ZHENG JIYE et al. Extraction of cultivated land lnformation in Shandong province based on Landsat8 OLI Data[J]. Chinese Agricultural Science Bulletin,2014,30(34):264-269. (In Chinese)

Interpretation study of crop planting structure in north China plain based on landsat 8 OLI data

LU Huixiong1,2, TIAN Yanliang3, DONG Shuangfa1,WANG Bin1,LI Mingsong1, NIU Haiwei1, ZHANG En1,2, XUE Qin1,2， ZHANG Xinyi1, CHENG Sisi1

(1.Airborne Survey and RemoteSensing Center of Nuclear Industry，Shijiazhuang 050002， China； 2.Faculty of Earth Sciences，East China Institute of Technology，Nanchang 330013， China; 3.The Institute of Hydrogeology and Environmental Geology ，Shijiazhuang 050061， China)

The use of the north China plain in 2013 65 king Landsat 8 OLI data, by ENVI software image preprocessing, sampling object-oriented man-machine interactive information extraction and interpretation methods, the north China plain (including wheat, corn, rice, soybeans distribution) and planting of crops and other major structures were interpreted system, the division of the study area of wheat and corn cropping area, wheat rice cropping area of rice cropping area, wheat soy cropping area, a ripe area corn, rice and a cooked zone distribution, and analysis of the impact of crop planting structure factors in the study area. The results show that the techniques and methods based on Landsat 8 OLI data for major crop acreage and planting large areas of remote sensing interpretation structure is feasible. It can provide an effective method for the study of crop planting structure.

planting structure; remote sensing interpretation; Landsat8 OLI; north China plain

2016-07-01 改回日期：2016-07-21

中國(guó)地質(zhì)調(diào)查局資助項(xiàng)目(1212011220941)

盧輝雄(1988-)，男，碩士，工程師，主要從事礦產(chǎn)資源調(diào)查評(píng)價(jià)、遙感技術(shù)應(yīng)用研究等工作，E-mail：1551310706@qq.com。

1001-1749(2017)03-0416-09

S 127

A

10.3969/j.issn.1001-1749.2017.03.19