基于Python的批量简单栅格计算

需求

现有一些tif格式的MODIS EVI栅格数据,需要进行单位转换乘以0.0001。

代码

import arcpy
arcpy.env.workspace = "I:\\MODIS\\EVI\\test"
rasterlist = arcpy.ListRasters("*","tif")
output_path = "I:\\MODIS\\EVI\\test2\\"
for raster in rasterlist:
   print str(raster)
   out = output_path + raster
   #Raster(raster)读取已有的栅格数据来建立栅格对象
   inRaster = arcpy.Raster(raster)
   outraster = inRaster * 0.0001
   outraster.save(out)
print("OK!")

小贴士

这里的save保存tif格式的栅格成功了,不明白setnull怎么就成功不了了。。。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。