Python读取栅格并计算

# -*- coding: utf-8 -*-
“““
Created on Wed Jan 21 19:29:08 2015

@author: Alvin
“““

import sys
import numpy
#output=open（r‘output.txt‘‘w‘）
#sys.stdout=output
from osgeo import gdal
from osr import SpatialReference
from functools import reduce
from operator import add
import itertools

def GetRasterInfo（rasterFile）:
    ds = gdal.Open（rasterFile）
    band = ds.GetRasterBand（1）
    nRows = band.YSize
    nCols = band.XSize

    projWkt = ds.GetProjection（）
    srs = SpatialReference（）
    srs.ImportFromWkt（projWkt）
    proj4Str = srs.ExportToProj4（）

    geotrans = ds.GetGeoTransform（）
    dx = geotrans[1]
    dy = geotrans[5]
    xmin = geotrans[0]
    ymax = geotrans[3]
    extent = （xmin ymax + nRows*dy xmin + nCols*dx ymax）

    data = band.ReadAsArray（）
    noDataValue = band.GetNoDataValue（）

    valRange = （band.GetMinimum（） band.GetMaximum（））

    return nRows nCols srs geotrans data noDataValue

def WriteGTiffFile（filename nRows nCols data geotransform srs noDataValue