map - 如何从 shapefile 的 .prj 文件中获取 Proj4 详细信息？

Question

我正在为我们的 gis 应用程序使用 mapdotnet 服务来加载 shapefile，这个 mapdotnet 服务需要 proj4 的详细信息。我从spatialreference.org 获取它们，但对于this projection proj4 详细信息是 blank .如何从 .prj 文件或 shapefile 中获取 proj4 详细信息？

下面是 shapefile 的 .prj:

PROJCS["NAD_1983_HARN_WISCRS_EauClaire_County_Feet",GEOGCS["GCS_North_American_1983_HARN",DATUM["D_North_American_1983_HARN",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",394000.0],PARAMETER["False_Northing",300812.797],PARAMETER["Central_Meridian",-91.28888888888889],PARAMETER["Standard_Parallel_1",45.87228112638889],PARAMETER["Scale_Factor",1.000035079],PARAMETER["Latitude_Of_Origin",45.87228112638889],UNIT["Foot_US",0.3048006096012192]]

Answer 1

您还可以使用此 Python 脚本(在 Internet 上的其他地方可见):

#!/usr/bin/env python

import osr
import sys

def main(prj_file):
    prj_text = open(prj_file, 'r').read()
    srs = osr.SpatialReference()
    if srs.ImportFromWkt(prj_text):
        raise ValueError("Error importing PRJ information from: %s" % prj_file)
    print srs.ExportToProj4()
    #print srs.ExportToWkt()

if __name__=="__main__":
    main(sys.argv[1])