《基于TM影像提取城区.doc》由会员分享,可在线阅读,更多相关《基于TM影像提取城区.doc(5页珍藏版)》请在金锄头文库上搜索。
1、基于TM影像提取城市区域提取说明Landsat5影像条带号124、036;Landsat5影像 TM3、TM4、TM5波段分别为红、中红、近红;数据与方法NDVIlandsat影像二值处理NDBI栅格计算建设用地数据采用landsat5的TM3、TM4、TM5波段。提取方法为指数提取。即在归一化裸露指数(NDBI)的基础上,用归一化植被指数(NDVI)去掉噪声点,通过对图像的二值处理、栅格计算。最终提取城市建设用地。研究表明提取的平均精度达到90%。提取的流程如下图所示:图1建设用地提取流程图Landsat5影像栅格计算NDBINDBINDBINDBI图2局部放大影像提取流程信息提取解释归一化
2、植被指数(NDVI)可以间接表达地表植被的覆盖。利用TM影像提取NDVI的公式为:归一化裸露植被指数(NBVI)可以很好的反映出地表面的裸露程度。利用TM影像提取NDBI的公式为:一般来说,NDVI0表示城市建设区域和低密度植被覆盖的裸地。将NDBI大于0的区域赋值为1,小于0的区域赋为0.然后再进行栅格计算,将二值化后的NDVI和NDBI相乘,值为1的区域即为城市建设用地。地物二值化处理表格如下:表1部分地物二值化处理表建设用地水体低植被覆盖区NDVI110NDBI101NDVI*NDBI100阈值改进考虑到影像成像时间处在春夏季,低密度植分类被容易混淆,因此可以通过设置分类的阈值来改进提取
3、区域的精度。根据精度需要,将NDVI的阈值设为0.1,即原值在00.1的范围内二值化后值为1,其余为0。设置NDBI的阈值为-0.1,即原值在-0.11的范围内二值化后值为1,其余为0。城市区域提取成果展示根据指数提取方法,2006年建设用地提取结果如下图:图3阈值为0,提取城市区域图4NDVI阈值为0.1,NDBI阈值为-0.1图5不同阈值提取结果对比图3、图4中,白色区域表示城市区域,图3表示阈值为0时,提取城市区域成果。图4表示NDVI阈值为0.1,NDBI阈值为-0.1时城市区域提取成果。图5为不同阈值提取成果叠加。改进与总结对比图3、4、5发现,设置阈值的不同对结果的重要影响。建议实际研究中修改阈值,做精度评价,客观选择阈值;图3中部分道路信息清晰可见。但提取结果与实际城区任有很大差别;图4中将部分水体裸地纳入城区,可以用归一化水体指数进一步修正。该方法克服了单一NDBI指数提取信息,排除了监督分类的主观干扰,效果客观可信。但提取阈值较难把握,需反复改进。提取成果影像仅能做参考,不能全面代表实际面积。