阿爾金山草地鼠害發生區及鼠荒地面積遙感估算

李培先 ,鄭江華,3,倪亦非 ,吳建國 ,吾買爾·吾守 ,艾合買江·阿吉 ,那松曹克圖

(1.新疆大學資源與環境科學學院，烏魯木齊 830046； 2.新疆大學智慧城市與環境建模普通高校重點實驗室，烏魯木齊 830046；3.綠洲生態教育部重點實驗室，烏魯木齊 830046； 4.新疆維吾爾自治區治蝗滅鼠指揮部辦公室，烏魯木齊 830001；5.新疆巴音郭楞蒙古自治州草原工作站，新疆庫爾勒 841000)

?

阿爾金山草地鼠害發生區及鼠荒地面積遙感估算

李培先1,2,鄭江華1,2,3,倪亦非4,吳建國4,吾買爾·吾守5,艾合買江·阿吉5,那松曹克圖5

(1.新疆大學資源與環境科學學院，烏魯木齊 830046； 2.新疆大學智慧城市與環境建模普通高校重點實驗室，烏魯木齊 830046；3.綠洲生態教育部重點實驗室，烏魯木齊 830046； 4.新疆維吾爾自治區治蝗滅鼠指揮部辦公室，烏魯木齊 830001；5.新疆巴音郭楞蒙古自治州草原工作站，新疆庫爾勒 841000)

【目的】運用“3S”技術對阿爾金山草地鼠害發生區和鼠荒地進行監測和解譯，獲得鼠害發生面積和分布，為研究區鼠害蔓延的防治及鼠荒地的治理提供參考。【方法】研究區Landsat8遙感影像，采用最大似然監督分類方法，初步分出鼠害發生區、草地、裸地、水體等地物，再添加海拔高度和植被覆蓋度的訓練規則，運用決策樹分類方法，進一步從鼠害發生區中提取鼠荒地、鼠害發生區，鼠害潛在發生區，最后生成專題圖。【結果】遙感估算得阿爾金山地區鼠害發生區面積為3 419.621 km2，其中鼠荒地面積1 376.476 km2，分別占到研究區的6.92%和2.78%，與2014年11月新疆巴州草原站實地調查估算結果相近，且有擴大趨勢。【結論】鼠荒地或鼠害潛在發生區有沿水源邊緣分布特征，鼠荒地與草地，裸地多呈鑲嵌分布。鼠荒地和裸地、草地共存易形成混合像元，造成解譯誤差；最大似然法和決策樹方法的結合使結果更符合實際調查。

“3S”技術; 阿爾金山；監督分類；決策樹；鼠荒地

0 引 言

【研究意義】草原鼠害是有害嚙齒類動物對草原資源及草原畜牧業產生的危害性影響[1]。草原鼠荒地指主要因嚙齒類動物活動和超載過牧等原因引起草原嚴重退化的次生裸地，其植被覆蓋度低于20%，地上植物生物量低于原生草原的20%[2]。草原鼠害監測預報僅依靠人工的采樣調查雖然有較高的準確性，但費時耗力，存在時效性差等缺點，難以適應大范圍實時監測的需求。結合多種現代技術有效監測鼠害發生區及鼠荒地面積與分布對研究區鼠害蔓延的防治、鼠荒地的治理以及草原植被生態恢復等具有重要意義。【前人研究進展】一些學者已運用“3S”技術(Geographic Information System，GIS；Remote Sensing，RS；Global Positional System，GPS)進行鼠害監測：黃建文等[3]利用TM影像，分析了天然梭梭林大沙鼠鼠害防治前后的土壤調節植被指數和植被生長狀況，為鼠害防治管理和動態監測提供科學依據。何詠琪等[4]運用“3S”技術，確定了高程、坡度、坡向、草地類型、土壤類型、5～10月增強型植被指數6個鼠害監測模型因子，計算因子權重值和量化值，建立鼠害監測模型，提取了不同鼠害發生區和鼠害危害區的閾值，最終生成青海省草原鼠害模擬分布區，結果與實地調查的鼠害區域基本一致。徐正剛等[5]基于MODIS數據，提出通過分析洞庭湖區植被指數變動規律以確定評價時間的方法，選擇增強型植被指數EVI進行鼠害的敏感性分析，確定了災害閾值，最終采用植被指數差值法獲得了研究區作物危害分布圖與危害等級分布圖。王瑋等[6]指出滿足草地鼠害監測要求的有TM和SPOT、Quick Bird等衛星遙感數據, TM影像是進行草原鼠蟲害的研究主要數據。阿爾金山鼠害發生區屬于新疆巴音郭楞蒙古自治州四大鼠害發生區之一[7]。該發生地分布有5個鼠害區，害鼠有灰尾兔(Lepusoiostolus)，大耳鼠兔(Ochotonamacrotis)、高原鼠兔(O.curzoniae)、柯氏鼠兔(O.koslowi)、喜馬拉雅旱獺(Marmotahimalayana)、白尾松田鼠(Microtus(Phaiomys)leucurus Blyth(1863))等 17 種，均為寒旱型種類[8]。其中優勢種白尾松田鼠有效洞口數 187.5 個/hm2，高原鼠兔有效洞口數 260個/hm2。相關高原鼠兔的生態學研究涉及與草甸植物群落的關系[9-12]、食物選擇[13]、冬季死亡率[14]和行為模式[15]等，以及阿爾金山自然保護區高原鼠兔夏季微生境選擇的主導因子分析[16]與食性顯微組織分析[17]。賈婷婷等[18]得出：隨高原鼠兔有效洞穴密度增加，樣方內出現的植物種數逐漸增多，說明高原鼠兔活動增加了高寒草甸植物群落的物種豐富度，但增加物種主要是有毒的乳漿大戟、圓葉筋骨草、藍花棘豆，以及適口性差的矮火絨草和蘭石草，草甸質量變劣。馬波等[19]得出：NDVI與洞穴數有顯著正相關的線性回歸關系:在無家畜競爭的自然環境中, 高原鼠兔有選擇地利用植被較好的生境而避免使用植被過度退化的生境。【本研究切入點】對高原鼠兔生態學特性的研究可為更準確分析鼠害的空間分布提供參考依據，然而對阿爾金山鼠害發生區及鼠荒地的發生面積以及空間分布的研究卻少有涉及。研究運用“3S技術”對阿爾金山草地鼠害發生區和鼠荒地進行監測和解譯，獲得鼠害發生面積和分布。【擬解決的關鍵問題】研究結合實地采集的GPS數據等，參照鼠害發生海拔、植被覆蓋度的生境特點，利用研究區的Landsat8 遙感影像，采用最大似然監督分類方法和決策樹分類方法提取鼠荒地、鼠害發生區，鼠害潛在發生區并統計其面積，分析鼠荒地和鼠害發生區的空間分布特征。

1 材料與方法

1.1 研究區概況

阿爾金山鼠害發生區主要位于阿爾金山自然保護區境內低于海拔4 800 m的區域，行政區劃隸屬于巴音郭楞蒙古自治州若羌縣和且末縣境內。在保護區中是以高原有蹄類野生動物為主的高原荒漠生態系統[20]。年平均降雨量 100～300 mm，主要集中在 6～8 月，形成明顯的雨季。年平均氣溫0～1 ℃，最熱月平均氣溫8 ℃，最低溫度在-30 ℃以下；土壤類型主要有山地棕漠土、高寒荒漠土、高寒荒漠草原土、高寒草原土和高寒草甸土等[16]。該發生區海拔在3 500 m以上，植物異常稀疏。4 000 m 以上的高山帶為墊狀優若黎、藏艾菊構成的寒漠，4 000 m以下的溝谷有水柏枝和小蔚構成的稀疏灌叢[21]。山地荒漠帶的上限高達 3 100～3 300 m，以合頭草荒漠為主，蒿屬荒漠次之[7]。主要植被類型包括紫花針茅、墊狀駝絨藜、硬葉苔草，以及棘豆、蒿子、冰草、草地早熟禾、芨芨草、合頭草、黑枸杞、沙拐棗、琵琶柴、駝絨藜、鳶尾等。隨著海拔升高和植被類型變化草地類型依次為、溫性荒漠草原、荒漠草原、高山草原、高山草原和高寒草原的過渡帶、高寒草原。表1，圖1

表1 阿爾金山鼠害發生研究區坐標范圍

注：研究區確定依據：1)包含實地調查路線和采樣點，包含2015年新疆畜牧廳治蝗辦提供的阿爾金山鼠害發生區。2)由高原鼠兔生境多發生在3 500～4 000 m海拔高度，研究區需覆蓋阿爾金山該海拔區間

Note: the study range was determined by two reasons:(1) The study area need to contain the routes of field survey, sampling points, and the data of rodent damage rangeland in Altun Mountain grassland in 2015 provided by Locust and Rodent Control Headquarters of Xinjiang. (2) The study area need to cover the altitude interval from 3500m to 4000m in Altun Mountain grassland where the habitats of plateau pika(Ochotonacurzoniae) exist

注：研究區范圍(綠色邊框)；鼠害發生區估算范圍，由巴州草原站和治蝗辦提供(淡藍色面)；野外調查路線和采樣點(紅線和黃色標記)；高原鼠兔和鼠害發生區照片(紅色邊框)

Note: the scope of the study area (green border).The estimation range of rodent damage rangeland was provided by Grassland Workstation of Bayinguoleng Mongolian Autonomous Prefecture (light blue area). Routes of field investigation and sampling points were shown with red lines and yellow marks separately. The photos of Plateau pika (Ochotonacurzoniae) and rodent damage rangeland were presented with red border

圖1 研究區位置及采樣點分布

Fig.1 The location and sampling point of study area

1.2 材 料

1.2.1 遙感數據及預處理

遙感數據為Landsat8 OLI數據，DEM數字高程數據為ASTER GDEM 30 m數據。表2

表2 研究區影像數據信息

遙感影像空間分辨率為30 m，時間分辨率為16 d，地圖投影坐標系為 UTM-WGS84。影像選擇考慮到其適用性(研究區云量<10%，無陰影條帶),與調查時間2015年9月13～15日最接近的影像。選擇日期在9月14日和9月5日的影像覆蓋研究區的主體部分，與調查時間很接近。而8月22日與8月27日影像僅覆蓋在研究區的兩個邊沿，與調查時間間隔小于1個月，對研究區鼠荒地解譯的影響可以接受。

歷史野外數據：新疆畜牧廳治蝗辦提供的2014年阿爾金山鼠害發生區和鼠荒地GPS數據。

圖像處理軟件為ENVI5.1，ArcGIS10.0和Google Earth Free 7.1.4。

遙感數據預處理：對已進行系統輻射校正和地面控制點幾何校正的Landsat8 OLI的影像進行輻射定標，圖像鑲嵌與裁剪，FLAASH大氣校正；對14景DEM數據進行鑲嵌與裁剪。

1.2.2 實地調查數據

經實地調查，鼠害發生區的優勢鼠種為高原鼠兔，處于高寒草原和高山草原的過渡帶，海拔基本在3 500 m以上，鼠害發生區主要植被類型：針茅、硬葉苔草、蒿子、冰草、草地早熟禾。徐瑤[22]得出在草場發生輕度退化的樣方中，鼠害的發生率為 100%。草場退化加劇，鼠害發生的概率則隨之減少。與此次調查的情況相同：在此次調查中發現的鼠害發生區，植被覆蓋度都較高，部分達到70%，對鼠害發生區不同草種光譜數據的采集與分析，發現有較明顯退化跡象。采用GPS采集研究區鼠害發生地、草地、鼠荒地、裸地等地物樣本，作為訓練樣本輔助后期遙感影像解譯。

1.3 方 法

1.3.1 最大似然監督分類(MLC)

基于參數化密度分布函數判別的最大似然方法MLC(Maximum Like-hood Classification)是遙感影像監督分類最常用的統計方法之一[23-24]。最大似然法與其他非參數方法比較具有易于與先驗知識融合和算法簡單而便于實施等優點[25],且具有較高的分類精度。

根據此次調查的實際情況和影像特點，將研究區地物分為：鼠荒地(粗分)、草地、裸地、水體、雪地、云和陰影。通過GPS數據、結合Google Earth進行訓練樣本的選擇。具體為將部分實地采集的草地、裸地的GPS面作為驗證樣本，用于精度評價；將實地采集的鼠害發生區樣本，結合巴州治蝗辦提供的鼠荒地坐標點在Google Earth繪制60 m×60 m鼠荒地樣本， 隨機均勻選擇一半為訓練樣本，另一半為驗證樣本，運用最大似然監督分類初步進行阿爾金山鼠荒地(粗分)解譯和精度驗證。

1.3.2 決策樹分類

決策樹方法主要是決策樹學習和決策樹分類兩個過程。決策樹學習過程是通過對訓練樣本進行歸納學習(Inductive learning) ,生成以決策樹形式表示的分類規則的機器學習(Machine learning)過程[26]。通過對訓練樣本進行決策樹學習生成決策樹,決策樹可以根據屬性的取值對一個未知樣本集進行分類,就是決策樹分類[27-28]。

巴州復雜的地理自然環境，為各種鼠類提供了有利的生存條件，隨著海拔高度和草地類型的垂直變化，鼠類的分布呈現出明顯的垂直分異[7]。李葉等[16]指出對阿爾金山高寒荒漠草原高原鼠兔生境選擇的主成分分析表明影響高原鼠兔生境選擇的主要因子依次是海拔、蓋度、植物種數和距道路距離。說明海拔高度是提取阿爾金山鼠害發生區的最重要因子。分析GPS數據對海拔高度的訓練規則如下。

(1)鼠害發生區(3 500～4 000 m)：由此次調查，以及新疆畜牧廳治蝗辦、巴州治蝗辦提供的鼠害發生區GPS坐標數據確定。(2)高海拔潛在發生區 (4 000～4 400 m):指發現鼠害有局部發生的高海拔地區。由此次實際調查有鼠害局部發生(4 026～4 069 m),結合巴州草原站與治蝗辦提供的鼠害局部發生GPS坐標(海拔4 409 m)確定。以及李葉等[16]于2012年7～8月在阿爾金山海拔高度為(4 122.25±139.48)m處做高原鼠兔穴居地樣本實驗。鼠害潛在發生區說明鼠害發生有局部擴大趨勢。(3)低海拔潛在發生區(3 250～3 500 m): 指發現鼠害有局部發生的低海拔地區。由巴州草原站與治蝗辦提供的鼠害局部發生GPS坐標確定(最低海拔3 238 m)。(4)未確定發生區(>4 250 m or <3 250 m)：高原鼠兔生境選擇在植被高度、植被密度、坡度、海拔高度、植被類型和土壤硬度等生態因子方面差異都有統計學意義(P<0.01)[16]。所以僅依靠海拔高度無法確定是否發生鼠害，定為不確定發生區。在該海拔高度內有適合高原鼠兔等鼠種適宜的生境，也可能變成鼠害潛在發生區。

由草原鼠荒地“植被覆蓋度低于20%”的定義[2]作為鼠荒地的提取規則。

2 結果與分析

2.1 最大似然監督分類方法粗分鼠害發生區

鼠荒地(粗分)包括鼠害發生區，鼠害潛在(局部)發生區，鼠荒地等，因為訓練樣本中包含鼠害發生區，鼠害潛在(局部)發生區，鼠荒地信息；另外分類中可能造成的鼠荒地與裸地混分，鼠害發生區與草地之間的混分等。作為初步分類。

從混淆矩陣中得出總體分類精度為86.67%，Kappa系數為0.81。鼠荒地的生產者精度為85.37%；用戶精度低，為31.25%，錯分誤差高。混合像元造成了鼠荒地與裸地、草地的錯分現象。從影像解譯圖中鼠荒地空間分布上看，也說明鼠荒地和草地，鼠荒地與裸地都有交集，鑲嵌分布，符合實際的調查情況。圖2

2.2 決策樹分類方法細分鼠害發生區和提取鼠荒地

對鼠荒地(粗分)結果進行基于海拔訓練規則的決策樹分類，分離具體鼠害發生區和未確定發生區，減小錯分影響。進一步對結果運用決策樹分類，添加植被覆蓋度小于20%的訓練規則，減小草地對鼠荒地造成的錯分誤差，從鼠害發生區中提取鼠荒地。圖3

圖2 2015年新疆阿爾金山地區鼠荒地潛在分布衛星影像解譯

圖3 2015新疆阿爾金山地區鼠害發生區分布衛星影像解譯

首先運用基于NDVI的像元二分法選擇不同的閾值計算植被覆蓋度，用ArcGIS計算出從0%～10%到90%～100% 10個等級；再結合實地調查目測法，照片比對，進行實際調查結果和不同閾值下計算結果的誤差分析。選擇的NDVIsoil和NDVIveg的閾值分別為累計貢獻率的5.2%和99.5%時各誤差率相對最低。表3

然后繪出植被覆蓋度等級圖。圖4

表3 不同閾值下基于NDVI像元二分法計算的草地覆蓋度及其誤差

圖4 2015年新疆阿爾金山地區植被覆蓋等級

最后，按決策樹方法進一步得出鼠荒地解譯圖并進行面積統計。得出2015年新疆阿爾金山地區鼠荒地空間分布圖，阿爾金山地區鼠害發生區面積為3 419.621 km2(512.943 17萬畝)，其中鼠荒地面積1 376.476 km2(206.471 43萬畝)，以及潛在鼠害發生區面積。將圖2、圖5與Google Earth影像進行比對發現鼠荒地及鼠害潛在發生區有沿阿牙克庫木湖與阿其克庫勒湖周邊分布以及沿山地河流邊緣分布的特征。 圖5，表4

2014年巴州草原站由實地調查GPS點勾畫初步計算的阿爾金山草原鼠害面積已達3 466.667 km2(520萬畝)，嚴重危害面積(鼠荒地)約1 333.333 km2(200萬畝)，且有繼續擴大的趨勢。這與2015年遙感估算得到的面積相近，但鼠害發生區及鼠荒地在空間分布上并非成片的面狀分布(圖1)，而是由生境等多因素決定的點狀或斑塊狀分布，呈現出與草地，裸地互有交集，鑲嵌分布特點。圖5

2015年鼠荒地遙感估算面積比2014年多出3.24%，為43.143 km2(6.471 43萬畝)，尤其鼠害潛在(局部)發生或活動區的大面積出現，說明鼠害有一定擴大蔓延的趨勢。

圖5 2015年新疆阿爾金山地區鼠荒地及鼠害發生區分布衛星影像解譯

鼠害地類型Thetypeofrodentdamagerangelands面積(km2)Area面積(萬畝)Area(104mu)占研究區比(%)Ratio鼠荒地137647620647143278鼠害發生區(3500～4000m)204314530647174413高海拔潛在發生區(4000～4400m)6334533950179951281低海拔潛在發生區(3250～3500m)87092313063842176未確定發生區(>4400mor<3250m)471203170680465953其它341078365116175466898

注：研究區總面積為49 444.944 km2( 7 416.741 65萬畝)。Landsat8 OLI影像單個像元的實際面積為900 m2≈0.001 km2故面積(km2)保留三位小數。實際工作中多使用(萬畝)單位，故進行換算；0.001 km2=0.000 15萬畝，面積(萬畝)保留五位小數

Note: the whole study area is 49,444.944 km2(7,416.741,65×104mu).The single pixel area of Landsat8 OLI image is 900 m2≈0.001 km2, so the area/ km2keep three decimal. The unit of mu is used more often in practical work because of 0.001 km2equal to 0.000,15×104mu. So the area/ 104mu keep five decimal

3 討 論

此次監測完成了研究區鼠害發生區和鼠荒地分布及面積的提取。但從影像方面看，30 m空間分辨率下只能通過適當波段組合后不同地物的顏色和深淺來區別，造成訓練樣本選擇的目視誤判。區分草地、裸地、鼠荒地的四景鑲嵌影像時間不同也對影像的自動解譯造成誤差。由于阿爾金山鼠害研究區屬無人區，海拔高，空氣較稀薄，山路崎嶇或無法通行，給樣本的采集帶來很大困難，造成嚴重鼠荒地訓練樣本少，對鼠荒地的解譯精度造成一定影響。可考慮引進無人機航拍技術對鼠害區進行監測，提高解譯精度。

4 結 論

4.1 從影像解譯效果和方法上:30 m中空間分辨率的landsat8 OLI數據針對鼠荒地和裸地的用戶精讀依然較低，最主要原因是鼠害發生區與草地、鼠荒地與裸地存在共存現象。30 m中空間分辨率的影像對阿爾金山草地和裸地區分基本可行，但共存現象和混合像元，造成較多錯分。針對鼠荒地高原鼠兔為主要鼠害，進行生境因子分析，將海拔高度以及植被覆蓋度作為訓練規則，運用最大似然監督分類和決策樹方法相結合提取鼠害發生區、鼠害潛在發生區以及鼠荒地。使結果更符合實際調查。

4.2 從鼠害發生區及鼠荒地的空間分布來看，鼠荒地或鼠害潛在發生區有沿水源邊緣分布的特征，且鼠荒地與草地，裸地互相有交集。與黃倩等[1]得出的高原鼠兔和大耳鼠兔主要沿河谷階地，陽坡或避風的洼地分布，鼠害草原多呈鑲嵌分布的結論相近。鼠害潛在分布區(4 000～4 400 m)、(3 250～3 500 m)的局部鼠害發生，其危害程度和面積估算等需要進一步實地驗證工作，為有效預防鼠害蔓延擴散提供依據。

4.3 從鼠害防治上看，阿爾金山草地鼠害發生區和鼠荒地面積大并有一定擴大趨勢，對當地草原生態環境勢必造成負面影響，需進行及時有效的鼠害監測和相應防治工作，同時鼠荒地宜進行合理的草地植被生態修復，以實現阿爾金山草原畜牧業的可持續發展。

致謝：感謝新疆維吾爾自治區產學研聯合培養研究生基地—新疆維吾爾自治區畜牧廳治蝗滅鼠指揮部辦公室和新疆巴州草原站與治蝗辦的大力支持與協助！

References)

[1]黃倩,花立民,曹慧,等.甘肅草原鼠害區劃研究[J]. 草業科學,2009,26(2):91-99.

HUANG Qian, HUA Li-ming, CAO Hui, et al. (2009). Study on grassland rodent pest division in Gansu Province [J].PrataculturalScience, 26(2):91-99. (in Chinese)

[2]唐川江,周俗,張新躍,等.四川省草原鼠荒地調查報告[J].草業科學,2007,24(7):58-61.

TANG Chuan-jiang, ZHOU Su, ZHANG Xin-yue, et al.(2007). Survey On rodent damaged rangeland in Sichuan Province [J].PrataculturalScience, 24(7):58-61. (in Chinese)

[3]黃建文,鞠洪波,特木欽,等. 阿拉善左旗天然梭梭林鼠害防治的遙感監測[J].林業科學, 2004,40(3):107-110.

HUANG Jian-wen, JU Hong-bo, TE Mu-qin ,et al.(2004).Monitoring on Controlling Effect of Rhombomys opimus in Haloxylon Ammodendron Bunge Stands Using Remote TM Imagery in Alashan [J].ScientiaSilvaeSinicae, 40(3):107-110. (in Chinese)

[4]何詠琪, 黃曉東, 侯秀敏,等. 基于3S技術的草原鼠害監測方法研究[J]. 草業學報, 2013,22(3):33-40.

HE Yong-qi，HUANG Xiao-dong，HOU Xiu-min, et al. (2013). Monitoring grassland rodents with 3S technologies [J].ActaPrataculturaeSinica, 22(3):33-40. (in Chinese)

[5]徐正剛, 趙運林, 李波,等. 基于MODIS植被指數評估洞庭湖區東方田鼠大暴發的危害[J]. 生態學報, 2014,34(23):7 101-7 109.

XU Zheng-gang, ZHAO Yun-lin, LI Bo, et al.(2014).Assessing Yangtze vole damage in Dongting Lake region of outbreak year based on MODIS imagery [J].AetaEeologieaSinica, 34(23):7,101-7,109. (in Chinese)

[6]王瑋, 馮琦勝, 于惠,等. "3S"技術在草地鼠蟲害監測與預測中的應用[J]. 草業科學, 2010,27(3):31-39.

WANG Wei, FENG Qi-sheng, YU Hui, et al.(2010). Application of "3S" technology on monitoring and prediction of the grassland rodents and insects [J].PrataculturalScience, 27(3):31-39.(in Chinese)

[7]李文利. 新疆巴音郭楞蒙古自治州草地優勢種害鼠發生危害區域劃分研究[J]. 草食家畜, 2013,(1):58-61.

LI Wen-li.(2013).Distribution of Dominant Rodents on Grassland in the Bayinguoleng Mongulia Autonomous Prefecture of Xinjiang [J].Grass-FeedingLivestock, (1):58-61. (in Chinese)

[8]鄭智民，姜志寬，陳安國.嚙齒動物學[M].上海：上海交通大學出版社，2008:83-114.

ZHENG Zhi-min, JIANG Zhi-kuan, CHEN An-guo. (2008).RodentsZoology[M]. Shanghai: Shanghai Jiaotong University Press:83-114. (in Chinese)

[9]孫飛達, 龍瑞軍, 路承香. 高原鼠兔不同洞穴密度對高寒草地植物群落組成及多樣性的影響[J]. 干旱區資源與環境, 2010,24(7):181-186.

SUN Fei-da, LONG Rui-jun, LU Chen-xiang. (2010). Effects of plateau pikas (Ochotonacurzoniae) burrow densities on plant community composition and population diversity in alpine meadow [J].JournalofAridLandResourcesandEnvironment, 24(7):181-186.(in Chinese)

[10] 劉偉, 張毓, 王溪,等. 高原鼠兔刈割行為與棲息地植物群落的關系[J]. 獸類學報, 2009,29(1):40-49.

LIU Wei, ZHANG Yu, WANG Xi, et al.(2009).The relationship of the harvesting behavior of plateau pikas with the plant community[J].ActaTheriologicaSinica, 29(1):40-49.(in Chinese)

[11] 劉偉, 李里, 嚴紅宇,等. 高原鼠兔挖掘活動對植物種的豐富度和地上生物量的影響[J]. 獸類學報, 2012,32(3):216-220.

LIU Wei, LI Li, YAN Hong-yu, et al.(2012).Effects of burrowing activity of plateau pika (Ochotonacurzoniae) on plant species richness and aboveground biomass[J].ActaTheriologicaSinica,32(3):216-220.(in Chinese)

[12] 周雪榮, 郭正剛, 郭興華. 高原鼠兔和高原鼢鼠在高寒草甸中的作用[J]. 草業科學, 2010,27(5):38-44.

ZHOU Xue-rong, GUO Zheng-gang, GUO Xing-hua, et al.(2010).The role of plateau pika and plateall zokor in alpine meadow [J].PrataculturalScience, 27(5):38-44.(in Chinese)

[13] 劉偉, 張毓, 王溪,等. 高原鼠兔冬季的食物選擇[J]. 獸類學報, 2009,29(1):12-19.

LIU Wei, ZHANG Yu, WANG Xi, et al. (2009). Food selection pattern for plateau pika in winter [J].ActaTheriologicaSinica, 29(1):12-19.(in Chinese)

[14] 楊虎虎, 曲家鵬, 劉力華,等. 食物及洞道密度對高原鼠兔冬季死亡率的影響[J].獸類學報,2011,31(3):257-264.

YANG Hu-hu, QU Jia-peng，LIU Li-hua, et al.(2011).Effects of food supply and burrow density on mortality of the plateau pika (Ochotonacurzoniae) in winter [J].ActaTheriologicaSinica, 31(3):257-264.(in Chinese)

[15] 張衛國, 劉蓉, 江小雷. 風險性聲訊信號對高原鼠兔行為模式的影響[J]. 草地學報, 2010,18(1):115-120.

ZHANG Wei-guo，LIU Rong，JIANG Xiao-lei. (2010). Influence of risk sound signal on behavior pattern of Pika [J].ActaAgrestiaSinica, 18(1):115-120.(in Chinese)

[16] 李葉, 王振宇, 張翔,等. 阿爾金山自然保護區高原鼠兔夏季微生境選擇的主導因子分析[J]. 中國媒介生物學及控制雜志, 2014,25(1):28-31.

Li Ye, WANG Zhen-yu, Zhang Xiang, et al. (2014). Analysis of dominant factors affecting microhabitat selection of plateau pika (Ochotonacurzoniae) during summer in Altun Mountain National Nature Reserve, Xinjiang Uygur Autonomous Region, China [J].ChineseJournalofVectorBiologyandControl, 25(1):28-31. (in Chinese)

[17] 熱娜古麗·艾合麥提, 李葉, 張翔,等. 阿爾金山國家級自然保護區高原鼠兔的食性顯微組織分析[J]. 中國媒介生物學及控制雜志,2015,(2):164-167.

Renagül Exmet, LI Ye, ZHANG Xiang, et al. (2015). Microscopic histological analysis of diet ofOchotonacurzoniaein the Altun Mountain National Nature Reserve [J].ChineseJournalofVectorBiologyandControl, 26(2):164-167.(in Chinese)

[18] 賈婷婷, 毛亮, 郭正剛. 高原鼠兔有效洞穴密度對青藏高原高寒草甸群落植物生態位的影響[J]. 生態學報, 2014, 34(4):869-877.

JIA Ting-ting, MAO Liang, GUO Zheng-gang.(2014).Effect of available burrow densities of plateau pika (Ochotonacurzoniae) on plant niche of alpine meadow communities in the Qinghai-Tibet Plateau [J].ActaEcologicaSinica, 34(4):869-877.(in Chinese)

[19]馬波, 王小明, 劉曉慶,等. 高原鼠兔洞穴數量與其棲息地植被分布格局的GIS分析[J]. 生物多樣性, 2011,19(1):71-78.

MA Bo, WANG Xiao-ming, LIU Xiao-qin, et al.(2011).GIS analysis of the spatial relationship between plateau pika burrow distribution and vegetation distributional patterns [J].BiodiversityScience, 19(1):71-78.(in Chinese)

[20] 劉世梁, 趙海迪, 董世魁,等. 基于SPOT NDVI的阿爾金山自然保護區植被動態變化研究[J]. 干旱區研究, 2014,31(5):832-837.

LIU Shi-liang, ZHAO Hai-di, DONG Shi-kui, et al.(2014).Dynamic of vegetation in the Altun Mountain Nature Reserve based on SPOT NDVI [J].AridZoneResearch, 31(5):832-837.(in Chinese)

[21]許鵬. 新疆草地資源及其利用[M].烏魯木齊：新疆科技衛生出版社，1993:349-357.

XU Peng.(1993).GrasslandResourcesandUtilizationinXinjiang[M].Urumqi: Health Science and Technology Publishing House in Xinjiang:349-357.(in Chinese)

[22]徐瑤. 藏北草地退化遙感監測與生態安全評價[D]. 成都:成都理工大學博士論文,2014.

XU Yao.(2014).RemotesensingdynamicmonitoringandEco-securityAssessmentofGrasslandDegradationinNorthernTibet[D]. PhD Thesis. Chengdu University of Technology，Chengdu. (in Chinese)

[23] Richards, John A. (1986). Remote Sensing Digital Image Analysis: An Introduction. Remote sensing digital image analysis:Springer-Verlag.

[24] Islam Abou EL-Magd Corresponding author, & T. W. Tanton. (2003). Improvements in land use mapping for irrigated agriculture from satellite sensor data using a multi-stage maximum likelihood classification.InternationalJournalofRemoteSensing, 24(21):4,197-4,206.

[25]陳亮, 劉希, 張元. 結合光譜角的最大似然法遙感影像分類[J]. 測繪工程, 2007,16(3):40-42.

CHEN Liang, LIU Xi, ZHANG Yuan. (2007). A study of image classification based on MLC combined with spectral angle [J].EngineeringofSurveyingandMapping, 16(3):40-42. (in Chinese)

[26]李德仁, 王樹良, 李德毅,等.論空間數據挖掘和知識發現的理論與方法[J].武漢大學學報(信息科學版), 2002,27(3):221-233.

LI De-ren, WANG Shu-liang, LI De-yi, et al.(2002).Theories and Technologies of Spatial Data Mining and Knowledge Discovery [J].GeomaticsandInformationScienceofWuhanUniversity, 27(3):221-233. (in Chinese)

[27] Hansen, M., Dubayah, R., & Defries, R. (2007). Classification trees: an alternative to traditional land cover classifiers.InternationalJournalofRemoteSensing, 17(5):1,075-1,081.

Fund project:Supported by the project commissioned by Locust Rodent Headquarters Office of Xinjiang Uygur Autonomous Region "Remote sensing monitoring of grassland biological disasters in Xinjiang" (2014-2015)and High level talent plan program of Xinjiang Uygur Autonomous Region (104-40002)

Estimating Area of Grassland Rodent Damage Rangeland and Rat Wastelands Based on Remote Sensing in Altun Mountain, Xinjiang, China

LI Pei-xian1,2, ZHENG Jiang-hua1,2,3, NI Yi-fei4, WU Jian-guo4, Wumaier Wushou5,Aihemaijiang Aji5,Nasongcaoketu5

(1. College of Resources & Environmental Sciences，Xinjiang University, Urumqi 830046, China;2.KeyLaboratoryofCityIntellectualizingandEnvironmentModelling,XinjiangUniversity,Urumqi830046,China; 3.KeyLaboratoryofOasisEcology,Urumqi830046,China; 4.LocustandRodentControlHeadquartersofXinjiangUygurAutonomousRegion,Urumqi830004,China；5.GrasslandWorkstationofBayingolinMongolAutonomousPrefecture,KorlaXinjiang841000,China)

【Objective】 Monitoring and interpreting grassland rodent damage rangeland and rat wastelands in Altun Mountain with 3S technologies can acquire area and distribution of grassland rodent in order to provide a reference in grassland rodents control and management.【Method】Landsat 8 remote sensing image of study area were classified preliminarily into pest occurrence area, grassland, bare land, body of water and others by supervising classification of maximum likelihood classification. Then rat wastelands, rodent damage rangeland and rodent damage potential (local) occurrence area were extracted from preliminary pest occurrence area by decision tree classifier based on rules about altitude and vegetation coverage further. Finally, the thematic map of grassland rodent damage rangeland and rat wastelands was made and occurrence area was calculated.【Result】The area of grassland rodent damage rangeland and rat wastelands in Altun Mountain were 3,419.621 km2and 1,376.476 km2, accounting for 6.92% and 2.78% of the study area, which is similar with the estimated result of Bazhou Prairie Station in November, 2014. And the area of grassland rodent damage rangeland has the trend of expansion.【Conclusion】It is a characteristic that rat wastelands, rodent damage rangeland and rodent damage potential(local) occurrence area are along with the edge of the water, and the rat wasteland, the bare land and the grassland are mostly of mosaic distribution. The coexistence of the wasteland and bare land, grassland is easy to lead to a mixed pixel interpretation. The combination of the maximum likelihood supervised classification and the decision tree method to extract rat wastelands, rodent damage rangeland and rodent damage potential (local) occurrence area makes the results more consistent with the actual survey

"3S" technologies; Altun Mountain; supervised classification; decision tree; rat wasteland

10.6048/j.issn.1001-4330.2016.07.023

2016-03-20

新疆維吾爾自治區治蝗滅鼠指揮部辦公室委托項目“新疆草原生物災害遙感監測”(2014～2015)；新疆維吾爾自治區高層次人才計劃(104-40002)。

李培先(1991-)，男，新疆伊寧人，碩士研究生，研究方向為草原生物災害遙感監測，(E-mail)15894623722@163.com

鄭江華(1973-)，男，浙江江山人，教授，研究方向為遙感與地理信息系統應用，(E-mail)zheng_jianghua@126.com

S443

A

1001-4330(2016)07-1346-10