Python遥感图像处理应用篇(十四)：GDAL 读取多光谱数据为二维数组并存入csv文件

之前写过一遍文章是采用Arcgis读取多光谱数据为csv文件的，相对来说有一些麻烦，本文有个更简单一点的方法，使用python直接读取。

本文使用Landsat05数据，采用GEE下载的多光谱数据，该数据下载下来包含19个波段。这些波段信息可以查看数据该数据详细说明。1-7是其对应的多光谱波段，第六波段是地表温度数据（ST),其他的是(SR) 地表反射率产品数据。数据显示如下（321波段）:

from osgeo import gdal
import pandas as pd

gdal库的安装可以参考：python中如何导入gdal包？_空中旋转篮球的博客-CSDN博客_gdal包

gdal.UseExceptions()
inputPathFile='O:/GLAS_GLH01_2003-2009/Shenlongjia/shenlongjia_landsat5_20030328.tif'
ds = gdal.Open(inputPathFile)
band1