#?-*-?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