landsat8波段叠加（layer stacking）
2024-10-08 07:13:19
许久没更。最近一直在看IDL,忽略了gdal的学习。
今天做了landsat8的辐射定标，需要通过reflectance gains/bias来进行波段运算。由于landsat8 oli未提供一个完整的叠加数据，所以首先第一步就是要将所有波段进行一个叠加。这里分享一下代码。
解释一下：
首先获取tiff驱动，通过第一个波段读取prj\geotransform,然后创建tiff文件导入这两个信息。
随后遍历文件，每次都打开波段文件，读取数组，随后通过GetRasterBand方法来新增波段，用WriteArray写入数组。
此外，ds.RasterCount可以获取栅格影像的波段数（ds为数据集）。
此外，其实不借助gdal也可以进行波段合成。只要有read/write功能的函数就可以，每次读取一个band,就readasarray，然后存储为一个多维数组再创建tif文件即可。
下次会补充IDL进行波段叠加的方法，毕竟对于笔者这么一个不会写正则表达式的人，相比python，还是在IDL里面进行辐射定标方便一点。
最新文章
- snmp ubuntu/centos--
- asp.net(C#)利用QRCode生成二维码(续)-在二维码图片中心加Logo或图像
- Config文件
- oracle学习笔记系列------oracle操作例子的专用表
- Connection的使用
- [zz]论程序员
- 脚本乐园 Shell中命令行选项和参数的处理
- 【BZOJ 3172】 [Tjoi2013]单词
- P2032 「Poetize9」升降梯上
- CNTV视频回看下载地址
- Ubuntu 16.04 设置MySQL远程访问权限
- Linux指令--watch，at
- 初学python之路-day06
- PHPStorm从入门到精通
- linux中安装gcc
- 斯坦福大学公开课机器学习：advice for applying machine learning | diagnosing bias vs. variance（机器学习：诊断偏差和方差问题）
- exp导出数据时丢表
- Codeforces Round #553 (Div. 2) D. Stas and the Queue at the Buffet 贪心+公式转化
- loglog()函数
- iframe解决ajax主域和子域之间的跨域问题