本文使用Python的GDAL库读取prj标识，构建空间参考对象。再将空间参考对象导出成WKT、Proj4和EPSG标识。
请参考如下代码：
#! /usr/bin/env python
import sys
from osgeo import osr
def prj2standards():
with open(prj_filepath, 'r') as f:
prj_txt = f.read()
srs = osr.SpatialReference()
srs.ImportFromESRI([prj_txt])
print('Shape prj is: ' + prj_txt)
print('WKT is: ' + str(srs.ExportToWkt()))
print('Proj4 is : ' + str(srs.ExportToProj4()))
srs.AutoIdentifyEPSG()
print('EPSG is: ' + str(srs.GetAuthorityCode(None)))
prj2standards(sys.argv[1])