import arcpy
from arcpy import env
from arcpy.sa import *
arcpy.CheckOutExtension("Spatial")
env.workspace = "H:\\Nth_MasterPaper\\nth_yfx\\nth"
rasterList = arcpy.ListRasters("*")
output_path = "H:\\Nth_MasterPaper\\nth_yfx\\gyhresult"
for raster in rasterList:
print raster
inRaster = raster
maxValueDS = arcpy.GetRasterProperties_management(inRaster, "MAXIMUM")
maxValue = maxValueDS.getOutput(0)
print "MAXIMUM" + str(maxValue)
minValueDS = arcpy.GetRasterProperties_management(inRaster, "MINIMUM")
minValue = minValueDS.getOutput(0)
print "MINIMUM" + str(minValue)
NormalizationRaster = (Raster(inRaster) - float(minValue)) / (float(maxValue)- float(minValue))
out = output_path + inRaster + '.tif'
NormalizationRaster.save(out)
print 'finish: ' + out
arcpy栅格归一化
最新推荐文章于 2024-04-14 18:02:29 发布