Get imagery samples from imagery layer with raster function applied

258
0
09-17-2018 05:11 AM
AlinPlesoianu3
New Contributor

Hi everyone!

I encountered a problem while attempting to use the get_samples method of the Imagery Layer from the Python API.

I accessed the Landsat Imagery layer, applied the NDVI raster function, but the get_samples method does not return the correct values.

This is the code

import arcgis

from arcgis.raster.functions import *

from arcgis.geocoding import geocode

## skipped connection lines

# access Landsat Imagery

landsat_item = gis.content.search('Multispectral Landsat', 'Imagery Layer', outside_org=True)[0]

landsat = landsat_item.layers[0]

# Create geometry

area = geocode('Bucharest, RO', out_sr=landsat.properties.spatialReference)[0]

geom = arcgis.geometry.Geometry(area['extent'])

# Get samples

ndvi_lyr = ndvi(landsat, '5 4')
sam = ndvi_lyr.get_samples(geom, sample_count=1)
sam[0]['value']  # get only the pixel values

The output is:

'1405 1198 1012 927 1367 1146 857 9 2720'

I expected to get the values as NDVI values (-1 to +1 interval), not unprocessed data.

Am I doing something wrong?

Thanks for any help,
Alin

0 Kudos
0 Replies