影像做除法时,要先进行float()运算,在进行Divide运算。以下将指定文件夹下的.tif文件首先转换为float并保存在指定文件,再进行divide运算。
import arcpy,os,glob arcpy.CheckOutExtension("spatial") datafiles=glob.glob("F:\\Modis_16\\1Moasic\\*.tif") for datafile in datafiles: outRaster_float ="F:\\Modis_16\\2float\\"+"flo_"+datafile[20:28] arcpy.gp.Float_sa(datafile, outRaster_float) outRaster_divide="F:\\Modis_16\\3divide\\"+"div_"+outRaster_float[23:32] arcpy.gp.Divide_sa(outRaster_float,"10000",outRaster_divide) print datafile+" has done!" print 'done....'
说明:此处没有用到Raster Calculator,对于栅格计算器的使用,还有待于学习。
时间: 2024-10-22 16:56:35