Arcpy对栅格数据的处理

       个人数据处理基本流程：转tif→乘因子系数→掩膜提取→统一投影→栅格转点→插值→掩膜提取→统计分析

1.hdf数据转tif

import os
import arcpy
from arcpy import env

inpath = "inpath"
outpath = "outpath\\"
arcpy.CheckOutExtension("Spatial")
env.workspace = inpath
arcpy.env.scratchWorkspace = inpath

hdfList = arcpy.ListRasters('*','HDF')
for hdf in hdfList:
    eviName = hdf[:-4] + ".tif"
    arcpy.ExtractSubDataset_management(hdf, outpath+ eviName)
    print hdf + "_ok"	
print 'All Done'

2.栅格计算

import arcpy
arcpy.env.workspace = "inpath"
rasterlist = arcpy.ListRasters("*","tif")
outpath = "outpath\\"

for raster in rasterlist:
  out = outpath + raster
  inRaster = arcpy.Raster(raster)
  outRaster = inRaster * 6
  outRaster.save(out)
  print out + "_