无忧文档

MODIS影像处理流程ok

同时期影像拼接、重投影、转格式(MRT)

黄土高原地区

Projected coordinate system name: Beijing_1954_3_Degree_GK_CM_108E

Geographic coordinate system name: GCS_Beijing_1954

Map Projection Name: Transverse Mercator

Scale Factor at Central Meridian: 1.000000

Longitude of Central Meridian: 108.000000

Latitude of Projection Origin: 0.000000

False Easting: 500000.000000

False Northing: 0.000000

北京54参考的椭圆:

椭球坐标参数:长半轴a=6378245m;短半轴=6356863.0188m

2、在Erdas中裁剪,得到研究区的原值NDVI -2000~9663

方法:aoi格式文件裁剪:在ERDAS图标面板工具条中单击Data Prep图标,Subset,打开Subset对话框。在Subset对话框中需要设置下列参数:

①输入文件名(input file)

②输出文件名(output file)

③单击aoi选中裁剪文件名:在choose aoi选项中选中文件名

④其余的我选择默认,然后ok等结果

3、得到范围在-1~1的NDVI

方法:ENVI —Basic Tools —Band Math(波段运算,把负值去掉)

在弹出的对话框Band Math中,Enter an expression:(b1 lt 0)*0+(b1 ge 0)*(b1*0.0001)。

这个公式意思就是:要是值小于0 就乘以0,使其变为0;同时,值大于等于0的话就乘以0.0001这个系数。这样就OK了。关于B1 代表所选的那个影像。

(对于陆地表面覆盖来说,云、雨、雪在可见光比近红外波段有较高的反射作用,所以NDVI 为负值;岩石、裸土的NDVI一般为0;有植被覆盖的地方一般大于0。)(b1*0.0001)

然后重新输出为tif格式,gis中classfied下查看。

在Erdas中波计算:

4、如何通过NDVI来计算植被覆盖度?

方法:一般都是用的像元二值模型。利用NDVI值得到植被覆盖度。

公式:fc = (NDVI - NDVIsoil)/( NDVIveg - NDVIsoil)

取NDVIsoil = NDVImin, NDVIveg = NDVImax。

(1)得到NDVI数据

(2)查出NDVI影像的最大最小值:在scroll 窗口右键—Quick Statistics

对于最小值,由于是计算植被覆盖度,植被的NDVI应该是大于0 ,因此应该选择大于0的最小值。

对于最大值,最大值不是1(红光波段反射率不等于0),因此应该选择数量大于0值小于1的最大值。

相关文档
热门文档
评论