arcgis裁剪多个shp文件_基于arcpy批量利用shapfile数据裁剪遥感数据并计算均值

我有3000多个美国县级矢量数据，每个县的矢量数据在单独文件夹中，表述如下：

其中文件夹与shp文件同名。

面临一个需求是，用shp裁剪遥感数据，大概有60张遥感影像，一个个裁剪显然不现实，对此查阅资料后设计了实现方式：

import arcpy
import os
import glob
from arcpy.sa import *

arcpy.CheckOutExtension("spatial")
arcpy.env.workspace = 'E:/data/'
inMaskfile = "E:/AmericaAllStateShapefile/"
rasterfile = "E:/COVID_19_data/TEMP_now/"
outputTXTfile = "E:/COVID_19_data/ceshi.txt"

for rastername in os.listdir(rasterfile):
    print rastername
    for filename in os.listdir(inMaskfile):
        # print filename
        for shpname in glob.glob(inMaskfile + filename + '/' + '*.shp'):
            outfile = arcpy.sa.ExtractByMask(rasterfile + rastername,shpname)

            mean = arcpy.GetRasterProperties_management(outfile,"MEAN")
            print filename + "_" + str(mean)
            f = open(outputTXTfile,'a+')
            f.write(filename + "_" + str(mean) + 'n')

这里用到了arcpy库，值得一提的是，该库不能在pip下载了，只能安装arcgis利用其自带的python2.7文件裁剪，具体为在pycharm setting中选择python2.7作为编译器即可。

再利用arcpy.GetRasterProperties_management获取裁剪后的均值，保存在txt即可。