GDAL에는 PROJ4 라이브러리에 대한 편리한 인터페이스가 있습니다.
GDAL Python 바인딩을 사용하여 Python에 확신이있는 경우 osr 클래스를 가져 오면 프로젝션 표현을 읽고 PROJ4, WKT, Esri .PRJ와 같은 다양한 형식으로 내보내는 매우 편리한 방법이 있습니다.
예를 들어이 스크립트는 shapefile의 .PRJ 파일을 WKT 및 PROJ4로 변환합니다 (마지막은 PostGIS에서 사용됨).
#! /usr/bin/env python
import sys
from osgeo import osr
def esriprj2standards(shapeprj_path):
prj_file = open(shapeprj_path, 'r')
prj_txt = prj_file.read()
srs = osr.SpatialReference()
srs.ImportFromESRI([prj_txt])
print 'Shape prj is: %s' % prj_txt
print 'WKT is: %s' % srs.ExportToWkt()
print 'Proj4 is: %s' % srs.ExportToProj4()
srs.AutoIdentifyEPSG()
print 'EPSG is: %s' % srs.GetAuthorityCode(None)
esriprj2standards(sys.argv[1])
명령 행에서 이것을 실행하십시오 :
$ python esriprj2standards.py /home/pcorti/data/shapefile/country.prj
Shape prj is: GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]
WKT is: GEOGCS["GCS_WGS_1984",DATUM["WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]
Proj4 is: +proj=longlat +datum=WGS84 +no_defs
EPSG is: 4326