上篇文章【Arcpy】相邻时相栅格数据批量作差的代码构建针对的是同一文件夹下的栅格文件,在env.workspace中设置工作空间,使用ListRasters获取栅格列表然后利用Python切片功能构造相邻文件作差的效果
本篇将解决不同文件夹下栅格文件的提取,以及代码构建过程和常见bug处理
需求
两个文件夹分别有地温和气温,求其对应日期的差,同时要注意单位换算和栅格大小
代码
# coding=utf-8
# 不同文件夹数据囧计算
import arcpy
from arcpy.sa import *
import glob
import os
output_path = r'H:\240329\LSTmiTa201703\tset'
input_path1 = r'H:\240329\Clip_LST201703\*.tif'
input_path2 = r'H:\240329\Tavg201703\*.img' # 这里最好不要把反斜杠放最后,文章会解释
rasterlist1