HOW TO
The instructions provided describe how to add individual raster bands from a multiband raster to the Extract Multi Values to Points geoprocessing tool using Python.
Note:
This article can be used as a workaround for the bug #NIM085128 - The Extract Multi Values to Points tool reports only the first band and repeats the value for subsequent bands at 10.1, when the input is a multiband raster.
Code:
import arcpy, os
from arcpy.sa import *
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
multibandRaster = r'C:\arcgis\ArcTutor\Raster\Data\Amberg_tif\090160.tif'
desc = arcpy.Describe(multibandRaster)
bands = desc.bandCount
inPointFeatures = r"C:\arcgis\ArcTutor\Raster\Data\points.shp"
in_rasters = []
for band in desc.children:
bandName = band.name
#append each band to the in_rasters list
in_rasters.append(os.path.join(multibandRaster, bandName))
# Execute ExtractValuesToPoints
ExtractMultiValuesToPoints(inPointFeatures, in_rasters)
Get help from ArcGIS experts
Download the Esri Support App