如何从shapefile .prj文件中获取proj4字符串或EPSG代码？ [重复]

#1 楼

shapefile应该具有定义投影的.prj文件。您可以将其与以下3个选项之一一起使用来获取proj4字符串，WKT定义或EPSG代码。

要获取proj4定义：

如果您有gdal安装在系统上后，可以使用gdalsrsinfo命令行应用程序将proj4定义作为OGC WKT定义获取：

$ gdalsrsinfo <shapefile>.prj

示例：

$ gdalsrsinfo 7490.prj 

PROJ.4 : '+proj=aea +lat_1=-18 +lat_2=-32 +lat_0=0 +lon_0=24 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs '

OGC WKT :
PROJCS["Albers_Equal_Area_Conic_South_Africa",
    GEOGCS["GCS_WGS_1984",
        DATUM["WGS_1984",
            SPHEROID["WGS_84",6378137,298.257223563]],
        PRIMEM["Greenwich",0],
        UNIT["Degree",0.017453292519943295]],
    PROJECTION["Albers_Conic_Equal_Area"],
    PARAMETER["False_Easting",0],
    PARAMETER["False_Northing",0],
    PARAMETER["longitude_of_center",24],
    PARAMETER["Standard_Parallel_1",-18],
    PARAMETER["Standard_Parallel_2",-32],
    PARAMETER["latitude_of_center",0],
    UNIT["Meter",1]]

要获取EPSG编号：

您可以使用prj2epsg.org服务上载.prj文件并获取匹配的EPSG代码列表。警告：prj2epsg.org返回最匹配的代码列表。

使用python：

如果您熟悉Python，则可以在Shapefile PRJ上使用非常有用的脚本进行PostGIS SRID查找表？以获得WKT，proj4字符串和EPSG代码。

另请参见：

您可以在spatialreference.org上找到大多数标准投影的定义。

#2 楼

gdal应用程序gdalsrsinfo生成srs的许多不同表示形式：

kyle@kyle-hp-desktop:~$ gdalsrsinfo epsg:4326 -o all

PROJ.4 : '+proj=longlat +datum=WGS84 +no_defs '

OGC WKT :
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AUTHORITY["EPSG","4326"]]

OGC WKT (simple) :
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433]]

OGC WKT (no CT) :
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433]]

ESRI WKT :
GEOGCS["GCS_WGS_1984",
    DATUM["D_WGS_1984",
        SPHEROID["WGS_1984",6378137,298.257223563]],
    PRIMEM["Greenwich",0],
    UNIT["Degree",0.017453292519943295]]

...

您可以输入proj.4字符串，wkt字符串或epsg代码。从gdal 1.9开始（我认为）。

#3 楼

我一直使用grep来获取此信息：

ogrinfo -so <shapefile>.shp | grep -e "ESPG"