最近要做个遥感相关的小系统,需要波段组合功能,网上找了可以使用ArcGIS安装时自带的arcpy包,但是Python3.7不能使用现有ArcGIS10.2版本,也不想再装其他版本,所以只能自己想了个办法解决。不过有点笨啊。 思路是: 1.读取需要组合遥感影像波段(此处用OLI) 2.创建数组,把读取的波段按序放进去 3.写入文件,写成tif多波段数据 上代码: from osgeo import gdal import os import numpy as np class GRID: #读图像文件 def read_img(self,filename): datas