我正在尝试将WGS84系统中的WGS84坐标重新投影到SIRGAS 2000 coord系统中的UTM坐标。使用GDAL,我首先将已知的utm坐标更改为相同坐标系(29 N)中的lat/lon对应的,以检查是否编写了正确的代码(这里省略了错误检查):monUtm.SetWellKnownGeogCS2000 UTM
// Rodovia dos Tamoios, B
我使用pyproj将一些坐标从WGS84转换为OSGB36。当映射出来时,转换的和未转换的不会映射到地图上的相同位置。OSGB36坐标被绘制在它们应该在的位置西北约100米处。下面是我在代码中所做的一个示例wgs84=pyproj.Proj("+init=EPSG:4326")row = pyproj.transform(wgs84, os