python將tif投影轉換至basemap

yanlisuo發表於2018-03-21

平面轉經緯度:osr.CoordinateTransformation

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

經緯度或平面轉平面時,(用GCS wgs84轉merc,緯度座標可正常轉,但經度需加2089780m左右,才與basemap底圖重疊。

經緯度轉平面:basemapobject(x,y)

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

相關文章