Extract a single band from a multi-band image

Posted on In
Symbols count in article: 0 Reading time ≈ 1 mins. 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
# -*- coding: utf-8 -*-
import arcpy
import os

def extractSingleBand(multiBandImgPath,singleBandSavePath):
    f_list = os.listdir(multiBandImgPath)
    for file in f_list:
        if os.path.splitext(file)[1] == ".tif":
            tif = str(multiBandImgPath) + "\\" + str(file)
            arcpy.env.workspace = tif
            for raster in arcpy.ListRasters():
                #The name of band to extract
                if str(raster) == "bandName":
                    bandrasteroutputname = singleBandSavePath+"bandName" + file
                    arcpy.CopyRaster_management(raster, bandrasteroutputname, "", -999, -999)
    print  "done"