import sys, os import arcpy from arcpy import * import glob input_path = r"E:\phdpaper\city1.shp" # shp路径 tif_path = r"E:\phdpaper\谋远的夜间灯光\【立方数据学社】DMSP-like2020.tif" # tif图像的路径 output_path = r"E:\phdpaper\gzp" # 输出路径 # # # 创建文件夹 if os.path.exists(output_path) == False: os.mkdir(output_path) # # # 定义工作空间 arcpy.env.workspace = r"E:\phdpaper\city1.shp" # 不用改,每一年的shp文件所在文件夹 # # # # 得到所有shp格式的城市矢量图层 rasterlist = arcpy.ListFiles("*.shp") # 收集每个shp图层到rasterlist变量里 shps = glob.glob(os.path.join(input_path, "*.shp")) # 收集每个shp图层的路径到shp变量里 # # for shp in shps: name = os.path.basename(shp).split(".")[0] + ".tif" # print (path_input) name_output = os.path.join(output_path, name) print (name_output) arcpy.gp.ExtractByMask_sa(tif_path ,shp ,name_output)