python将tif投影转换至basemap

平面转经纬度：osr.CoordinateTransformation

ct = osr.CoordinateTransformation(inproj, outproj)
ct.TransformPoint(x,y)

经纬度或平面转平面时，（用GCS wgs84转merc，纬度坐标可正常转，但经度需加2089780m左右，才与basemap底图重叠。

http://blog.sina.com.cn/s/blog_6e51df7f0100ui7n.html

http://pcjericks.github.io/py-gdalogr-cookbook/projection.html

经纬度转平面：basemapobject(x,y)

m = Basemap(projection=proje,，，，)
x,y=m(lon,lat)