许久没更.最近一直在看IDL,忽略了gdal的学习. 今天做了landsat8的辐射定标,需要通过reflectance gains/bias来进行波段运算.由于landsat8 oli未提供一个完整的叠加数据,所以首先第一步就是要将所有波段进行一个叠加.这里分享一下代码. 解释一下: 首先获取tiff驱动,通过第一个波段读取prj\geotransform,然后创建tiff文件导入这两个信息. 随后遍历文件,每次都打开波段文件,读取数组,随后通过GetRasterBand方法来新增波段,用Wr