python求投影距离_python地图投影转换

一、投影包osr与proj4的使用

1.osr投影转换示例

from osgeo import osr,ogr

#定义投影

#wgs84

source=osr.SpatialReference()

source.ImportFromEPSG(4326)

#google

target=osr.SpatialReference()

target.ImportFromEPSG(3857)

#简单投影转换

coordTrans=osr.CoordinateTransformation(source,target)

#点转换

coordTrans.TransformPoint(117,40)

#点数组转换

coordTrans.TransformPoints([(117,40),(117.5,39.5)])

#SF几何对象转换

g=ogr.CreateGeometryFromWkt("POINT(117 40)")

#转换前wgs84

print(g.ExportToWkt())

print(g.GetX(),g.GetY())

#转换后google

g.Transform(coordTrans)

print(g.ExportToWkt())

print(g.GetX(),g.GetY())

2.投影转换示例

from pyproj import Proj,Geod,transform

# projection 1: UTM zone 15, grs80 ellipse, NAD83 datum

# (defined by epsg code 26915)