또 다른 옵션은 ogr2ogr에서 재 투영 옵션을 사용하는 것입니다. JEQL, Fiona 또는 GeoScript 접근 방식보다 해커 접근 방식이지만, 그럼에도 불구하고 효과적입니다. s_srs와 t_srs에 사용 된 투영 사이에 잘못된 이스트와 노스 링이 변경되는 한, 시작 및 투영은 실제로 원본 모양 파일의 실제 투영 일 필요는 없습니다. 이 예에서는 Google Mercator를 사용하고 있습니다. 기본으로 사용할 훨씬 간단한 좌표계가 있다고 확신하지만 복사 / 붙여 넣기가 바로 앞에 있습니다.
ogr2ogr -s_srs EPSG:900913 -t_srs 'PROJCS["Google Mercator",GEOGCS["WGS 84",DATUM["World Geodetic System 1984",SPHEROID["WGS 84",6378137.0,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0.0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.017453292519943295],AXIS["Geodetic latitude",NORTH],AXIS["Geodetic longitude",EAST],AUTHORITY["EPSG","4326"]],PROJECTION["Mercator_1SP"],PARAMETER["semi_minor",6378137.0],PARAMETER["latitude_of_origin",0.0],PARAMETER["central_meridian",0.0],PARAMETER["scale_factor",1.0],PARAMETER["false_easting",1000.0],PARAMETER["false_northing",1000.0],UNIT["m",1.0],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","900913"]]' -f "ESRI Shapefile" shift.shp original.shp
또는 입력 / 붙여 넣기를 저장하려면 다음을 저장하십시오 projcs.txt
(위와 동일하지만 작은 따옴표를 제거).
-s_srs EPSG:900913
-t_srs PROJCS["Google Mercator",GEOGCS["WGS 84",DATUM["World Geodetic System 1984",SPHEROID["WGS 84",6378137.0,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0.0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.017453292519943295],AXIS["Geodetic latitude",NORTH],AXIS["Geodetic longitude",EAST],AUTHORITY["EPSG","4326"]],PROJECTION["Mercator_1SP"],PARAMETER["semi_minor",6378137.0],PARAMETER["latitude_of_origin",0.0],PARAMETER["central_meridian",0.0],PARAMETER["scale_factor",1.0],PARAMETER["false_easting",1000.0],PARAMETER["false_northing",1000.0],UNIT["m",1.0],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","900913"]]
그런 다음 다음을 실행하십시오.
ogr2ogr --optfile projcs.txt shifted.shp input.shp