二十二、Arcpy批量波段组合——结合Landat数据城市建成区提取

一、前言

其实波段组合和GIS中栅格计算有点类似，实质上就是对每个像素点对应的DN值进行数学计算，也就是可以进行运算表达式是三个或多个变量相加、相减……每一个变量对应于一个图像数据，对这三个或多个图像数据求值并输出结果图像。

 二、具体操作

1、实验具体目标

将2021年Landast8数据波段进行组合

其实就是Band=band1+band2+.....band9

（2）筛选数据

只需要对文件夹中的band1-band9波段进行计算，所以你现在只需要把三景文件夹影像的对应波段筛选出来，分别放置一个文件夹里面即可

特别提醒：在这里我只是以三景影像进行波段合成计算，如果数据很多，你完全可以放置很多个波段影像至文件夹 

（3）创建Python文件，新建好程序文件之后就是写入代码，创建Python具体操作就不演示。

# -*- coding: utf-8 -*-
import os, arcpy, datetime
# 读取和存放文件夹的位置
layers