影像做除法时,要先进行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....'