python实现栅格数据最大值合成法

今天需要用python实现栅格数据的最大值合成法，辗转颇久，偶然搜到了CellStatistics。哈哈，先用arcpy实现吧，有余力了再自己写（估计猴年马月了）。
根据多个栅格数据计算每个像元的统计数据。可用的统计数据有：众数、最大值、均值、中位数、最小值、少数、范围、标准差、总和及变异度。
在本示例中，将针对多个输入 Grid 栅格中的每个像元计算标准差，并将结果输出为 Grid 栅格。
CellStatistics （独立脚本）

# Name: CellStatistics_Ex_02.py
# Description: Calculates a per-cell statistic from multiple rasters
# Requirements: Spatial Analyst Extension

# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *