首页 > 其他 > 详细

landsat8波段叠加(layer stacking)

时间:2019-11-28 01:05:43      阅读:414      评论:0      收藏:0      [点我收藏+]

许久没更。最近一直在看IDL,忽略了gdal的学习。

今天做了landsat8的辐射定标,需要通过reflectance gains/bias来进行波段运算。由于landsat8 oli未提供一个完整的叠加数据,所以首先第一步就是要将所有波段进行一个叠加。这里分享一下代码。

技术分享图片

 

解释一下:

首先获取tiff驱动,通过第一个波段读取prj\geotransform,然后创建tiff文件导入这两个信息。

随后遍历文件,每次都打开波段文件,读取数组,随后通过GetRasterBand方法来新增波段,用WriteArray写入数组。

 

此外,其实不借助gdal也可以进行波段合成。只要有read/write功能的函数就可以,每次读取一个band,就readasarray,然后存储为一个多维数组再创建tif文件即可。

下次会补充IDL进行波段叠加的方法,毕竟对于笔者这么一个不会写正则表达式的人,相比python,还是在IDL里面进行辐射定标方便一点。

 

landsat8波段叠加(layer stacking)

原文:https://www.cnblogs.com/ljwgis/p/11946336.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!