arcpy对影像批量求和

最新推荐文章于 2024-07-25 14:17:35 发布

贾村拓哉

最新推荐文章于 2024-07-25 14:17:35 发布

阅读量125

点赞数

文章标签： arcgis

本文链接：https://blog.csdn.net/Jschao/article/details/129258067

版权

arcpy计算LST白天和黑夜的温度和

import os
import re
import arcpy
from arcpy.sa import *
import os

arcpy.env.workspace = "D:/LST_8day_eartdata"
for i in range(408):
    files_day = os.listdir("D:/LST_8day_eartdata/Day")
    file_day = files_day[i]
    raster_day = Raster("D:/LST_8day_eartdata/Day/"+file_day)
    print(file_day)
    file_out = re.split('.tif$', file_day)[0]+ "_A.tif"

    files_night = os.listdir("D:/LST_8day_eartdata/Night")
    file_night = files_night[i]
    raster_night = Raster("D:/LST_8day_eartdata/Night/"+file_night)
    print(file_night)
    # Check out the ArcGIS Spatial Analyst extension license
    arcpy.CheckOutExtension("Spatial")
    ave_day_night = ((raster_day+raster_night)*0.01-273.15)*8
    #save
    out_raster = os.path.join("D:/LST_8day_eartdata/tif/", file_out)
    print(out_raster)
    ave_day_night.save(out_raster)