最近要做个遥感相关的小系统,需要波段组合功能,网上找了可以使用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
1