Standardising rasters in arcpy

DanielFuller

Hello

 

I have to to standardise some rasters using a loop and following equation

X1= (X - X.min())/(X.max()-X.min())

where X is the original dataset, and X1 the new standardized dataset, which has a minimum value of 0 and a maximum value of 1.

I was wondering how I would go about doing this? the current code I have is 

for calc in projected_rast:
    print calc
    minLSTresult = arcpy.GetRasterProperties_management(calc, "MINIMUM")
    maxLSTresult = arcpy.GetRasterProperties_management(calc, "MAXIMUM")
    standardised = ((calc - minLSTresult)/(maxLSTresult-minLSTresult))
    standardised.save()

 

were projected_rast is a list of strings to my rasters

Luke_Pinner

a Raster object has maximum and minimum properties (if statistics have been calculated).

for calc in projected_rast:
    print(calc)
    raster = arcpy.Raster(calc)
    minLSTresult = raster.minimum
    maxLSTresult = raster.maximum
    standardised = ((raster - minLSTresult)/(maxLSTresult-minLSTresult))
    standardised.save()