IDL實現植被覆蓋度計算

本實例實現植被覆蓋度計算，利用NDVI數據像元二分模型。

像元二分計算方法是將像元植被覆蓋結構分爲純像元和混合像元兩種。純像元被植被完全覆蓋，覆蓋度爲1；混合像元由植被和非植被部分組成，其NDVI值是植被覆蓋部分非NDVI值與非植被覆蓋部分的NDVI值線性加權，公式如下：

                                               

    上式中：NDVI爲像元的NDVI值；fv爲像元的植被覆蓋度；NDVIv和NDVI0分別爲植被覆蓋部分和非植被覆蓋部分的NDVI值，上式可推導爲：

                                                                

注意：在水體較多的區域爲了避免水體面積過大對NDVI取值的影影響，可將水體去除(判定條件NDVI<0)。

                                          

IDL核心代碼：

function cal_VegFraction,NDVI,NDVI0,NDVIv
   ;計算植被覆蓋度
   ;參數NDVI爲歸一化植被指數，NDVI0和NDVIv分別爲純裸地和純植被的NDVI值
   ;像元二分法
   result=(NDVI-NDVI0)/(NDVIv-NDVI0)
   ;將NDVI小於，NDVI0（即fv<0）或者水體區域(NDVI<0)的像元植被覆蓋度值改爲0
   w=where(NDVI le NDVI0 or NDVI le 0)
   result[w]=0
   w=where(NDVI ge NDVIv)
   result[w]=1
   return,result
end

數據來源：地理空間數據雲樣例數據(參考：http://www.rscloudmart.com/)

參考：1.徐永明老師的遙感二次開發語言IDL。