【ArcGIS遇上Python】栅格影像批量除以10000

最新推荐文章于 2023-12-20 15:29:43 发布

刘一哥GIS

最新推荐文章于 2023-12-20 15:29:43 发布

阅读量5.7k

点赞数 5

分类专栏：《ArcGIS遇上Python》文章标签： ArcGIS Python 影像 divide

本文为博主「刘一哥GIS 」原创文章，未经博主允许不得转载或爬取，侵权必究！

本文链接：https://blog.csdn.net/lucky51222/article/details/38931211

版权

《ArcGIS遇上Python》专栏收录该内容

该专栏为热销专栏榜第18名

39 篇文章 772 订阅 ¥89.90 ¥99.00

订阅专栏

本文介绍如何利用Python脚本，对ArcGIS中的tif影像文件进行浮点转换，然后批量除以10000的操作。脚本将分别完成影像的浮点转换保存和除法运算。

摘要由CSDN通过智能技术生成

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