基于Python的批量簡單柵格計(jì)算

需求

現(xiàn)有一些tif格式的MODIS EVI柵格數(shù)據(jù)，需要進(jìn)行單位轉(zhuǎn)換乘以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)讀取已有的柵格數(shù)據(jù)來建立柵格對象
   inRaster = arcpy.Raster(raster)
   outraster = inRaster * 0.0001
   outraster.save(out)
print("OK!")

小貼士

這里的save保存tif格式的柵格成功了，不明白setnull怎么就成功不了了。。。

最后編輯于 ：
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成，瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明：文章內(nèi)容（如有圖片或視頻亦包括在內(nèi)）由作者上傳并發(fā)布，文章內(nèi)容僅代表作者本人觀點(diǎn)，簡書系信息發(fā)布平臺，僅提供信息存儲服務(wù)。

友情鏈接更多精彩內(nèi)容