这是提出NDBI的那个人的文章,有兴趣的可以下载下来看 这是他的技术路线图,他有做了验证,我们这里只是实现,不做验证。 这是NDBI方法的公式:我们这里用L8 OLI数据做实验,也就是对应的MIR是第6波段,NIR是第5波段。 后面实验操作就是按照他的步骤! 原始数据需要经过辐射定标和大气校正,如上图所示,我使用整景的L8数据做实验,上面是辐射定标结果。 为了节省时间,我使用快速大气校正,但是提醒大家:最好使用FLAASH那个完整模块进行校正,完整大气校正结果是地表发射率数据。 这是大气校正结果,数据预处理结果。 但是简单FLAASH的结果并不是地表反射率,需要通过转化:b1*0.0001. ![在这里插入图片描述](https://img-blog.csdnimg.cn/20210331104710429.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3l1YW5ZVzc1NTY=,size_16,color_FFFFFF,t_70)
使用BAND MATH 工具,如上图所示。 结果如图,可以看出已经转化为地表反射率数据。 ![在这里插入图片描述](https://img-blog.csdnimg.cn/20210331104501467.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3l1YW5ZVzc1NTY=,size_16,color_FFFFFF,t_70)
band math 工具位置在上图所示的地方。 使用band math 计算NBDI,千万不要选错数据。 ![在这里插入图片描述](https://img-blog.csdnimg.cn/20210331104922651.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3l1YW5ZVzc1NTY=,size_16,color_FFFFFF,t_70)
这是计算的归一化建筑指数(NDBI)的结果。 ![在这里插入图片描述](https://img-blog.csdnimg.cn/202103311050070.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3l1YW5ZVzc1NTY=,size_16,color_FFFFFF,t_70)
还是使用Band math 结合IDL的运算符进行后续处理。 上图所示:白色是提取的建筑物,黑色是植被、水体等其他的东西。 使用上图所示工具查看,提取的是否准确。 个人感觉:提取结果有点粗糙,因为好像把裸地和路面也给提取出来了。肯定还需要后续操作,才能保证提取精度。
|