批量处理LAI指数并求栅格内的均值--arcgis

寒江雪（王训

# Import system modules
from arcpy import env
from arcpy.sa import *
import os
arcpy.env.overwriteOutput = True
env.workspace="E:/model_data/LAI"
arcpy.CheckOutExtension("Spatial")

# read list name for LAI
flai=open("E:/model_data/LAI/lai.txt","r")
dirs_lai=flai.read().splitlines()
flai.close()

#read list name for landuse
fland=open("E:/model_data/LAI/land.txt","r")
dirs_land=fland.read().splitlines()
fland.close()

for i in dirs_lai:
    if   os.path.splitext(i)[1]=='.tif':
         a=i
         # remove 0 values for LAI
         c=Con(Raster(str(a))>0,Raster(str(a)))
         for j in dirs_land:
                 if os.path.splitext(j)[1]=='.tif':
                         b=j
                         outname1=os.path.join("E:/model_data/LAI/result",a.split('.')[0]+"_L"+b)
                         print outname1
                         #raster calculator
                         outtimes=Times(c,b)
                         # Zonal analysis
                         inZoneData ="36km_lambert.shp"
                         zoneField="FID"
                         inValueRaster=outtimes
                         outTable=a.split('.')[0]+b.split('.')[0]+".dbf"
                         print outTable
                         outZSaT = ZonalStatisticsAsTable(inZoneData, zoneField,inValueRaster,outTable,"DATA","MEAN")