侯玉婷 李令軍 姜 磊 武鳳霞

(1.十堰市環(huán)境保護監(jiān)測站，湖北 十堰 442000；2.北京市環(huán)境保護監(jiān)測中心，北京 100048)

當前，灰霾污染已成為公眾關(guān)注的主要環(huán)境問題之一。冬小麥收獲季節(jié)的秸稈焚燒是導(dǎo)致每年6月前后主產(chǎn)區(qū)嚴重灰霾污染事件的重要原因[1-2]。秸稈焚燒排放大量的污染物，主要包括PM2.5、NOX、NH3、CH4、揮發(fā)性有機化合物(VOC)、CO、SO2，嚴重影響區(qū)域空氣質(zhì)量。獲取冬小麥秸稈焚燒污染物的排放量及空間分布，有助于客觀了解冬小麥秸稈焚燒空氣污染現(xiàn)象。

目前，計算區(qū)域尺度秸稈焚燒大氣污染物排放量的方法主要有兩種：一是依靠秸稈焚燒火點監(jiān)測的過火面積進行測算[3]1591，[4]，獲取過火面積、單位面積作物產(chǎn)量、谷草比、秸稈污染物排放因子等參數(shù)，得到最終的秸稈焚燒污染物排放量；二是利用統(tǒng)計數(shù)據(jù)估算[5]1826-1827，[6]3046-3048，獲取年鑒作物種植面積統(tǒng)計數(shù)據(jù)、單位面積作物產(chǎn)量、谷草比、秸稈焚燒比例、秸稈污染物排放因子等參數(shù)。然而，火點監(jiān)測受云遮蓋及衛(wèi)星回訪周期影響，僅能反映部分時段、部分區(qū)域的秸稈焚燒事件[3]1592-1593，得到的污染物排放量不完整；利用統(tǒng)計數(shù)據(jù)估算的秸稈焚燒污染物排放量，僅能從宏觀上了解秸稈焚燒污染排放總體情況，不能提供更精細化的網(wǎng)格污染物排放分布信息。

本研究基于MODIS衛(wèi)星250 m分辨率的16 d合成歸一化植被指數(shù)(NDVI)數(shù)據(jù)，提出了一種冬小麥種植區(qū)遙感快速提取方法，并建立了基于解譯的冬小麥種植區(qū)秸稈焚燒大氣污染排放的測算方案。以冬小麥種植第一大省河南省為例，遙感提取了2010年全省冬小麥的種植面積和空間分布，對比統(tǒng)計年鑒中的數(shù)據(jù)，對提取結(jié)果進行了驗證、分析；在此基礎(chǔ)上，結(jié)合調(diào)研獲取的大氣污染物排放因子等參數(shù)，估算了冬小麥秸稈焚燒的污染排放情況，進一步得到了秸稈焚燒污染物排放量及其空間分布。本研究設(shè)計的方法具有數(shù)據(jù)易獲取、過程簡便的特點，結(jié)果客觀、可靠，能為秸稈焚燒監(jiān)管工作以及其他區(qū)域的相關(guān)研究提供參考應(yīng)用和技術(shù)支持。

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

總體技術(shù)流程分為兩個部分(見圖1)：一是提取冬小麥種植范圍，設(shè)計了一種簡便、高效的提取方法，過程采用物候規(guī)則、非監(jiān)督分類和人工識別結(jié)合的方式，選用250 m分辨率的MODIS NDVI數(shù)據(jù)，并參考統(tǒng)計年鑒數(shù)據(jù)對提取結(jié)果進行驗證分析；二是估算冬小麥種植區(qū)秸稈焚燒大氣污染物排放量，在冬小麥種植區(qū)遙感提取結(jié)果的基礎(chǔ)上，結(jié)合調(diào)研獲取的秸稈產(chǎn)量、秸稈焚燒大氣污染物排放量等參數(shù)，進一步計算得到冬小麥種植區(qū)秸稈焚燒大氣污染物排放量及其空間分布。

圖1 總體技術(shù)流程

1.1 數(shù)據(jù)資料

主要采用MODIS NDVI數(shù)據(jù)，其他資料還包括高分衛(wèi)星影像、統(tǒng)計年鑒中的數(shù)據(jù)。MODIS NDVI數(shù)據(jù)來自美國國家航空航天局(NASA)MODIS數(shù)據(jù)分發(fā)中心，空間分辨率為250 m，本研究采用16 d合成NDVI數(shù)據(jù)，該數(shù)據(jù)采用降低云覆蓋影響的合成方法。本研究共獲取2010年覆蓋河南全境的23期16 d合成NDVI數(shù)據(jù)。

統(tǒng)計數(shù)據(jù)主要取自《河南省2011年統(tǒng)計年鑒》，用于冬小麥種植區(qū)遙感提取面積驗證。數(shù)據(jù)分省、市(18個)、縣(108個)3個尺度。其中，省、市和10個縣的冬小麥種植區(qū)面積為統(tǒng)計年鑒結(jié)果；由于年鑒中未提供其他98個縣的冬小麥種植區(qū)面積，故采用年鑒中提供的相關(guān)數(shù)據(jù)進行折算。

1.2 冬小麥種植范圍的提取

1.2.1 基于物候規(guī)則提取冬小麥種植區(qū)

冬小麥種植面積的估算是遙感和農(nóng)業(yè)領(lǐng)域的熱點課題[6]3046，[7]。主流方法之一即是以MODIS NDVI數(shù)據(jù)為基礎(chǔ)進行分類提取。與面向?qū)ο蠓诸?、混合像元分解[8]等方法相比，基于MODIS NDVI數(shù)據(jù)的方法雖然在分辨率和精度上有所不足，但具有簡便易行、數(shù)據(jù)資源豐富等優(yōu)點，因而被廣泛采用。具體方法有：采用冬小麥關(guān)鍵生長期的幾期NDVI圖像做簡單運算，或設(shè)定適當閾值的提取方法[9-12]；建立物候規(guī)則模型計算的方法[13-14]；NDVI閾值和物候規(guī)則模型兩者相結(jié)合的方法[15-20]；引入其他來源的土地利用分類結(jié)果，劃定農(nóng)業(yè)種植區(qū)范圍，再在范圍內(nèi)利用NDVI閾值或物候規(guī)則提取冬小麥種植面積分布信息[21-22]的方法。

但現(xiàn)有方法存在以下問題：(1)選用的MODIS NDVI數(shù)據(jù)時相少，信息量不豐富；(2)設(shè)置NDVI閾值需要豐富的野外調(diào)查工作經(jīng)驗，而且受時間差異性和地域局限性影響，不能推廣應(yīng)用[23]；(3)現(xiàn)有土地利用分類結(jié)果內(nèi)再分類的方法實用性差，因為高質(zhì)量的土地利用分類結(jié)果不易獲取，而且冬小麥種植區(qū)的提取結(jié)果精度和時效性也直接依賴外源土地利用分類數(shù)據(jù)。

本研究設(shè)計提出的冬小麥種植區(qū)面積提取方法，保留了數(shù)據(jù)易獲取、方法簡單高效的特點，且提高了分類準確性，降低了人工經(jīng)驗差別對分類結(jié)果的影響，具體步驟為：

(1) 根據(jù)冬小麥生長典型物候特征，結(jié)合全年MODIS NDVI的數(shù)據(jù)情況，制定物候規(guī)則，據(jù)此提取冬小麥種植區(qū)掩膜。

(2) 利用全年的MODIS NDVI數(shù)據(jù)做非監(jiān)督分類，參考高分衛(wèi)星影像，人工識別篩選出含有冬小麥種植區(qū)的有效類別。

(3) 疊加利用物候規(guī)則提取的冬小麥種植區(qū)掩膜及非監(jiān)督分類并人工識別得到的冬小麥種植區(qū)有效類別，最終確定冬小麥種植區(qū)范圍和面積。

1.2.2 非監(jiān)督分類提取冬小麥有效類別

基于2010年共23期的MODIS NDVI數(shù)據(jù)，利用ArcGIS軟件做非監(jiān)督分類，將所有柵格劃分為20類。將分類結(jié)果與河南省境內(nèi)2010年冬小麥生長期內(nèi)的高分衛(wèi)星影像疊加，人工識別剔除完全不屬于冬小麥種植區(qū)的類別，保留其余屬于冬小麥種植區(qū)的類別或部分含有冬小麥種植區(qū)的混合類別。

1.2.3 冬小麥種植區(qū)范圍的確定

通過對比高分衛(wèi)星影像發(fā)現(xiàn)，采用物候規(guī)則得到的冬小麥種植區(qū)掩膜(見圖2)將部分水體、水體周邊地塊、城市建筑用地、林地等誤判為冬小麥種植區(qū)，是由于這些地塊沒有明顯的NDVI物候變化特征，或是受到其他隨機因素、混合像元等的影響。同時，采用非監(jiān)督分類的方法可以有效剔除部分水體、林地和城市建筑用地，但因分類總數(shù)較少，且未加入先驗樣本，部分冬小麥像元和其他植被覆蓋像元不能完全分開，而被判為同一個混合類別，因而非監(jiān)督分類得到的冬小麥種植區(qū)范圍偏大。

圖2 基于物候規(guī)則提取的冬小麥種植區(qū)掩膜

將基于物候規(guī)則提取的冬小麥種植區(qū)掩膜(見圖2)和非監(jiān)督分類提取的冬小麥種植區(qū)有效類別(見圖3)疊加，以取其交集的方式，可以同時剔除被物候規(guī)則誤判的屬于水體、林地、城市建筑用地的柵格，以及非監(jiān)督分類混合類別中非冬小麥種植區(qū)的柵格。

1.3 冬小麥秸稈焚燒大氣污染物排放估算

圖3 非監(jiān)督分類提取的冬小麥種植區(qū)有效類別

在提取冬小麥種植范圍的基礎(chǔ)上，設(shè)計了一套測算秸稈焚燒大氣污染物排放清單的模型。首先基于調(diào)研獲取的單位面積冬小麥產(chǎn)量、谷草比、秸稈焚燒比例和排放因子參數(shù)，測算單位面積冬小麥秸稈焚燒大氣污染物排放量，得到研究區(qū)域內(nèi)冬小麥秸稈焚燒大氣污染物排放量。計算方法見式(1)至(2)。

Pi=k×Gi×A

(1)

(2)

式中：Pi為冬小麥秸稈焚燒大氣污染物i排放量，g；k為面積校正系數(shù)；Gi為單位面積冬小麥秸稈焚燒大氣污染物i排放量，g/hm2；A為冬小麥種植面積，hm2；Y為單位面積冬小麥產(chǎn)量，kg/hm2；η為秸稈焚燒比例，%；Ri為污染物i排放因子，g/kg；Fgs為谷草比。

在250 m的中分辨率尺度上提取冬小麥種植區(qū)，因為田間地壟、鄉(xiāng)村道路甚至袖珍村落無法徹底剔除以及混合像元的影響，使遙感結(jié)果仍會略大于實際的冬小麥種植面積。假設(shè)面積的高估程度在研究區(qū)內(nèi)相對均勻，則可根據(jù)統(tǒng)計數(shù)據(jù)對遙感提取結(jié)果進行校正。根據(jù)遙感提取的冬小麥種植面積與統(tǒng)計年鑒中種植面積，計算高估的數(shù)量。相關(guān)計算見式(3)。

(3)

式中：S0為統(tǒng)計年鑒中的冬小麥種植面積，hm2；A0為遙感提取得到的冬小麥種植面積，hm2。

2 結(jié)果與討論

2.1 冬小麥種植區(qū)提取結(jié)果

根據(jù)上述方法，將基于物候規(guī)則提取的冬小麥種植區(qū)掩膜和非監(jiān)督分類得到的冬小麥種植區(qū)有效類別求交集，即可得到本研究冬小麥種植面積的最終監(jiān)測結(jié)果。圖4為2010年河南省冬小麥種植區(qū)遙感最終結(jié)果。結(jié)果顯示，2010年河南省冬小麥種植區(qū)主要集中在豫東平原和南陽盆地，在豫西山地、豫南山地以及鄭州等分布較少。此外，大范圍種植區(qū)內(nèi)散落著面積大小不等的城市及村落，形成冬小麥種植的空白區(qū)。

圖4 2010年河南省冬小麥種植區(qū)遙感最終結(jié)果

對比河南全省、各市、各縣的冬小麥種植面積統(tǒng)計年鑒結(jié)果(《河南省2011年統(tǒng)計年鑒》顯示為52 879 km2)及本研究提取的冬小麥種植區(qū)面積統(tǒng)計結(jié)果，分析冬小麥種植區(qū)提取結(jié)果的準確性，見表1，其中A1、A2、A3分別代表僅采用物候規(guī)則的提取結(jié)果、僅采用非監(jiān)督分類并篩選出有效類別的結(jié)果、A1與A2求交集得到的冬小麥種植區(qū)遙感最終結(jié)果；縣域尺度-1包含統(tǒng)計年鑒中直接提供冬小麥種植面積的10個縣，縣域尺度-2包含河南省全部108個縣。

表1 冬小麥種植區(qū)提取結(jié)果評價

結(jié)果顯示，在市、縣尺度上，A3與統(tǒng)計數(shù)據(jù)的相關(guān)系數(shù)都在0.9以上、平均偏差均在13%以內(nèi)，說明提取結(jié)果總體上精度較好。對比A3與A1、A2統(tǒng)計的偏差面積、相關(guān)系數(shù)，可見通過疊加物候規(guī)則和非監(jiān)督分類結(jié)果的方式，可以大幅減少面積偏差。

參考中國生態(tài)系統(tǒng)與生態(tài)功能區(qū)劃數(shù)據(jù)庫中的河南地貌矢量數(shù)據(jù)，按面積最大法將河南省內(nèi)的108個縣分為4類：平原區(qū)、丘陵區(qū)、臺地區(qū)、山地區(qū)，具體劃分見圖5。統(tǒng)計了不同地貌類型上冬小麥種植區(qū)面積提取結(jié)果與統(tǒng)計年鑒值之間的關(guān)系，見表2。

圖5 河南省地貌分類

地貌類型項目A1A2A3平原區(qū)相關(guān)系數(shù)0．960．930．91平均偏差/%614732丘陵區(qū)相關(guān)系數(shù)0．930．950．94平均偏差/%524825臺地區(qū)相關(guān)系數(shù)0．900．940．93平均偏差/%52661山地區(qū)相關(guān)系數(shù)0．810．830．76平均偏差/%309741

結(jié)果表明，平原區(qū)A3提取結(jié)果與統(tǒng)計年鑒值的相關(guān)系數(shù)達到0.91，丘陵區(qū)達到0.94，臺地區(qū)達到0.93，山地區(qū)較低，但也能達到0.76，證明本研究冬小麥種植區(qū)的遙感提取方法在不同地貌類型上都是有效的。對比A3與A1、A2統(tǒng)計的平均偏差、相關(guān)系數(shù)，說明在平原區(qū)、丘陵區(qū)、臺地區(qū)均可以通過疊加物候規(guī)則和非監(jiān)督分類結(jié)果的方式，達到在不影響相關(guān)性的情況下有效減少面積偏差的效果。在山地區(qū)，本研究方法對面積影響較大，致使面積偏小41%，原因是在MODIS 250 m的分辨率下，山地區(qū)坡耕地呈碎片狀分布，耕地信息受混合像元的影響過重而不能有效反演所致。因為我國冬小麥種植區(qū)主要分布在江、淮、海河流域的平原區(qū)，本研究方法在冬小麥的主產(chǎn)平原區(qū)提取結(jié)果與統(tǒng)計年鑒值的相關(guān)系數(shù)達到0.91，表明該方法對于我國冬小麥整體種植面積的估計有非常好的借鑒意義。

表3 秸稈焚燒污染物排放因子

2.2 污染物排放估測結(jié)果

根據(jù)式(3)，統(tǒng)計年鑒中2010年冬小麥種植面積為52 879 km2，提取的冬小麥種植面積為66 759 km2，故面積校正系數(shù)取0.79。根據(jù)文獻調(diào)研，單位面積冬小麥產(chǎn)量取5 838 kg/hm2(來自河南省2010年的統(tǒng)計年鑒)；冬小麥谷草比取1.366[24]。秸稈焚燒比例取值參考農(nóng)業(yè)部科技教育司2010年發(fā)布的《全國農(nóng)作物秸稈資源調(diào)查與評價報告》，取18.72%[25]。秸稈污染物排放因子由文獻[26、27]調(diào)研獲得，選取國內(nèi)常用的取值，詳見表3。

根據(jù)上述參數(shù)，采用式(1)、式(2)計算得到250 m分辨率的2010年河南省冬小麥秸稈焚燒污染物排放量空間分布。結(jié)果表明，250 m分辨率下冬小麥種植區(qū)的秸稈焚燒單位柵格污染物年排放量：PM2.5為154.1 kg，NOX為9.9 kg，NH3為5.1 kg，CH4為10.7 kg，VOC為62.0 kg，CO為363.4 kg，SO2為1.6 kg。

基于2010年河南省秸稈焚燒污染物排放量空間分布結(jié)果，計算2010年河南省冬小麥秸稈焚燒排放總量，分別如下：PM2.5為16.482×104t，NOX為1.057×104t，NH3為0.549×104t，CH4為1.141×104t，VOC為6.635×104t，CO為38.880×104t，SO2為0.169×104t。

表4對比了鄭有飛等[3]1592計算的2007年夏季河南省秸稈焚燒大氣污染物排放量與本研究得到的河南省2010年秸稈焚燒大氣污染物排放量、曹國良等[5]1827-1829計算的2000年河南省秸稈焚燒大氣污染物排放量。結(jié)果顯示，在數(shù)量級上，本研究的計算結(jié)果與文獻[5]的結(jié)果接近，但比之更小，因為后者結(jié)果包含了冬小麥以外其他作物的秸稈焚燒大氣污染物排放。文獻[3]的結(jié)果與本研究結(jié)果存在數(shù)量級的差異，分析原因是火點遙感監(jiān)測手段僅能反映部分時段、部分區(qū)域的秸稈焚燒事件，導(dǎo)致估計的秸稈焚燒大氣污染物排放量測算嚴重偏小。

表4 河南省秸稈焚燒大氣污染物排放情況

2.3 討 論

本研究提出一種新的秸稈焚燒大氣污染物排放量的遙感測算方法，基于MODIS衛(wèi)星250 m分辨率的16 d合成的NDVI數(shù)據(jù)，采用物候規(guī)則、非監(jiān)督分類和人工識別結(jié)合的方式，解譯提取冬小麥的種植面積與分布范圍。單純物候規(guī)則提取的冬小麥種植區(qū)結(jié)果與統(tǒng)計值的相關(guān)性較好，但面積總量偏大；單純采用非監(jiān)督分類結(jié)合識別有效類別的方法，能區(qū)分冬小麥、其他地物類別以及混合像元，但結(jié)果受數(shù)據(jù)質(zhì)量、聚類中心的位置、聚類數(shù)目的影響很大，平均偏差很高，面積也有所高估。采用疊加兩者并取交集的方式，使最終結(jié)果的平均偏差為3組中最小，與統(tǒng)計結(jié)果較接近。

本研究在提取冬小麥種植區(qū)的基礎(chǔ)上，結(jié)合調(diào)研獲取的秸稈焚燒污染物排放因子等參數(shù)，測算了網(wǎng)格化的冬小麥秸稈焚燒大氣污染物排放量及其空間分布，得到的數(shù)據(jù)結(jié)果為250 m分辨率下空間分布特征。但結(jié)果受限于柵格間污染物排放量均勻的假設(shè)，不能體現(xiàn)污染物排放強度的空間差異性。未來可以通過計算單位柵格的冬小麥秸稈生物量等，細化秸稈焚燒參數(shù)的空間差異，從而進一步提高秸稈焚燒大氣污染物排放量的測算精度。

3 結(jié) 論

(1) 本研究以MODIS衛(wèi)星250 m分辨率的16 d合成的NDVI數(shù)據(jù)為基礎(chǔ)，設(shè)計了一種冬小麥種植區(qū)提取方法，并以河南為例進行了全省域的2010年種植面積與分布范圍的遙感提取。結(jié)果表明，遙感提取的冬小麥面積分布與統(tǒng)計年鑒中的數(shù)據(jù)相關(guān)系數(shù)都在0.9以上、平均偏差均在13%以內(nèi)，且在各種地貌類型上都有較好的效果。

(2) 250 m分辨率下2010年河南省冬小麥種植區(qū)的秸稈焚燒單位柵格污染物年排放量：PM2.5為154.1 kg，NOX為9.9 kg，NH3為5.1 kg，CH4為10.7 kg，VOC為62.0 kg，CO為363.4 kg，SO2為1.6 kg。

(3) 2010年河南省冬小麥秸稈焚燒排放的污染物總量：PM2.5為16.482×104t，NOX為1.057×104t，NH3為0.549×104t，CH4為1.141×104t，VOC為6.635×104t，CO為38.880×104t，SO2為0.169×104t。

[1] 嚴文蓮，劉端陽，孫燕，等.秸稈焚燒導(dǎo)致的江蘇持續(xù)霧霾天氣過程分析[J].氣候與環(huán)境研究，2014，19(2)：237-247.

[2] 楊占婷，白永清，陳英英，等.鄂中東部兩次秸稈焚燒引發(fā)的嚴重灰霾天氣分析[J].氣象科技，2013，41(6)：1147-1152.

[3] 鄭有飛，田宏偉，陳懷亮，等.河南省夏季秸稈焚燒污染物排放量的估算與分析[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報，2010，29(8).

[4] DENNIS A，F(xiàn)RASER M，ANDERSON S，et al.Air pollutant emissions associated with forest，grassland，and agricultural burning in texas[J].Atmospheric Environment，2002，36(23)：3779-3792.

[5] 曹國良，張小曳，王亞強，等.中國區(qū)域農(nóng)田秸稈露天焚燒排放量的估算[J].科學(xué)通報，2007，52(15).

[6] 朱佳雷，王體健，鄧君俊，等.長三角地區(qū)秸稈焚燒污染物排放清單及其在重霾污染天氣模擬中的應(yīng)用[J].環(huán)境科學(xué)學(xué)報， 2012，32(12).

[7] 權(quán)文婷，王釗.冬小麥種植面積遙感提取方法研究[J].國土資源遙感，2013，25(4)：8-15.

[8] 潘耀忠，李樂，張錦水，等.基于典型物候特征的MODIS-EVI時間序列數(shù)據(jù)農(nóng)作物種植面積提取方法——小區(qū)域冬小麥實驗研究[J].遙感學(xué)報，2011，15(3)：578-594.

[9] 鄒金秋，陳佑啟，SATOSHI U，等.利用Terra/MODIS數(shù)據(jù)提取冬小麥面積及精度分析[J].農(nóng)業(yè)工程學(xué)報，2007，23(11)：195-200.

[10] 王云峰，沈潤平，田國珍.基于MODIS數(shù)據(jù)光譜突變法提取冬小麥種植面積研究[J].內(nèi)蒙古氣象，2009，32(6)：18-21.

[11] 競霞，劉良云，張超，等.利用多時相NDVI監(jiān)測京郊冬小麥種植信息[J].遙感技術(shù)與應(yīng)用，2005，20(2)：238-242.

[12] 趙文亮，賀振，賀俊平，等.基于MODIS-NDVI的河南省冬小麥產(chǎn)量遙感估測[J].地理研究，2012，31(12)：2310-2320.

[13] 王茂新，裴志遠，吳全，等.用NOAA圖像監(jiān)測冬小麥面積的方法研究[J].農(nóng)業(yè)工程學(xué)報，1998，14(3)：89-93.

[14] 李穎，陳秀萬，段紅偉，等.多源多時相遙感數(shù)據(jù)在冬小麥識別中的應(yīng)用研究[J].地理與地理信息科學(xué)，2010，26(4)：47-49.

[15] 黃青，吳文斌，鄧輝，等.2009年江蘇省冬小麥和水稻種植面積信息遙感提取及長勢監(jiān)測[J].江蘇農(nóng)業(yè)科學(xué)，2010，28(6)：508-511.

[16] 楊小喚，張香平，江東.基于MODIS時序NDVI特征值提取多作物播種面積的方法[J].資源科學(xué)，2004，26(6)：17-22.

[17] 喬紅波，張慧，程登發(fā).不同時序EOS/MODIS-NDVI監(jiān)測河南省冬小麥面積[J].安徽農(nóng)業(yè)科學(xué)，2008，36(27)：11940-11941.

[18] 康凌艷，雷玉平，鄭力，等.在GIS支持下利用MODIS數(shù)據(jù)監(jiān)測多種作物和果樹種植面積[J].遙感技術(shù)與應(yīng)用，2007，22(3)：361-366.

[19] 李紅梅，張樹譽，王釗.MODIS衛(wèi)星NDVI時間序列變化在冬小麥面積估算中的應(yīng)用分析[J].氣象與環(huán)境科學(xué)，2011，34(3)：46-49.

[20] 隋學(xué)艷，朱振林，李少昆，等.基于種植制度利用MODIS數(shù)據(jù)提取冬小麥種植面積[J].農(nóng)業(yè)工程學(xué)報，2010，26(增刊1)：225-229.

[21] 李寄，黃進良，許文波.湖北省冬小麥種植面積遙感估算方法研究[J].世界科技研究與發(fā)展，2008，30(5)：597-599.

[22] 馮美臣，楊武德，張東彥，等.基于TM和MODIS數(shù)據(jù)的水旱地冬小麥面積提取和長勢監(jiān)測[J].農(nóng)業(yè)工程學(xué)報，2009，25(3)：103-109.

[23] YOU Xingzhi，MENG Jihua，ZHANG Miao，et al.Remote sensing based detection of crop phenology for agricultural zones in China using a new threshold method[J].Remote Sensing，2013，5(7)：3190-3211.

[24] 曹國良，張小曳，鄭方成，等.中國大陸秸稈露天焚燒的量的估算[J].資源科學(xué)，2006，28(1)：9-13.

[25] 農(nóng)業(yè)部科技教育司.全國農(nóng)作物秸稈資源調(diào)查與評價報告[R].北京：農(nóng)業(yè)部科技教育司，2010.

[26] STREETS D G，BOND T C，CARMICHAEL G R，et al.An inventory of gaseous and primary aerosol emissions in Asia in the year 2000[J].Journal of Geophysical Research：Atmospheres，2003，108(D21)：8809-8831.

[27] 祝斌，朱先磊，張元勛，等.農(nóng)作物秸稈燃燒PM2.5排放因子的研究[J].環(huán)境科學(xué)研究，2005，18(2)：29-33.