专题七谢俊鸿.docVIP

基于遥感的泰安市生态环境评价 一、数据源 ETM+影像（含第8波段），1:5万数字地形图，1:5万行政界线数据。 二、数据预处理 影像裁切，数据配准和影像融合。 影像裁切 erdas→Dataprep→subset image→input file→out file→选择AOI文件→OK 数据配准 对遥感数据进行几何精校正 影像融合 erdas→interpreter→utilities→layer stack→input file添加landsat1、2、3、4、5、7、8波段→output file→OK 三、生态因子生成 1.植被因子：植被覆盖度（由NDVI转换） （1）计算NDVI erdas→interpreter→spectral enhancement→indices→input file(融合后的文件)→output file→select function（NDVI）→OK （2）计算植被覆盖度 计算植被覆盖度Fv采用的是混合像元分解法，将整景影像的地类大致分为水体、植被和建筑，具体的计算公式如下： FV = (NDVI- NDVIS)／(NDVIV - NDVIS) 其中，NDVI为归一化差异植被指数，当NDVI小于等于0时，fv=0;当0NDVI=1时fv=图像原值。 2.土壤因子： 土壤指数（GRABS=VI-0.09178*BI-5.58959，VI、BI指缨帽变换的绿度和亮度） （1）缨帽变换 erdas→interpreter→spectral enhancement→tasseled cap→input file(融合后的文件)→output file→OK 计算土壤指数 建立模型计算土壤指数： 算法语言为：$n1_cap(2)-0.09178*$n1_cap(1)-5.58959 3.坡度生成 （1）由等高线生成TIN 用arcmap打开等高线文件添加地理坐标和投影信息： arcview→file→extensions→3d analyst(选中) 添加等高线文件→surface→create TIN from features→height source(选major7)→OK 由TIN生成DEM、坡度 surface→Derive slope→output grid extent(选same as display)→output grid cell size(输入15)→OK 四、生态环境评价 1.生态因子的归一化处理 （1）植被因子，据植被覆盖度分为10个等级。 Erdas建模分级 算法语言为： 结果如图所示 （2）土壤因子，土壤指数依据大小划分10等级。（下表进攻参考，不作标准） 数值 -28157 ~ -20000 -20000 ~ -15000 -15000 ~ -10000 -10000 ~ -5000 -5000 ~ -3500 -3500 ~ -2500 -2500 ~ -1500 -1500~ -500 -500 ~ 300 300 ~ 772.43 编码 1 2 3 4 5 6 7 8 9 10 2.生态环境指数计算 采用指数加权叠加方法，计算环境指数，其中植被因子的权重为0.7，土壤因子的权重为0.2，地形因子的权重为0.1。 Erdas建模 算法语言为： $n1_fvfenjicut * 0.7+$n2_turangzhishufenji*0.2+$n3_slopefenjicut*0.1 3.结果输出 等价结果分级，研究区掩膜。 成果展示，如图所示