¿Son estas dos proyecciones OGC WKT iguales? Y si es así, ¿por qué son diferentes sus cadenas de referencia?

Question

Así que aquí está la .prj de texto de una ley federal de shapefile conjunto de datos de todos los aeropuertos públicos en su nativa OGC WKT formato:

GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]

Aquí está el OGC WKT de texto de un spatialreference.org partido (http://spatialreference.org/ref/epsg/4326/ogcwkt/) que creo que es, efectivamente, la misma proyección. La mayoría del tiempo las dos cadenas son una coincidencia exacta, así que no estoy completamente seguro de que:

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.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]

Como se puede ver, de un alto nivel, estos dos aspecto astucia similares. Pero a partir de un bajo nivel, parece posible que un desajuste en la cadena de ids de referencia podría ser una señal de que efectivamente son diferentes..

Answer 1

El oficial de la OGC "el conocido Texto de la Representación de Sistemas de Referencia Espacial" para EPSG 4326 (http://spatialreference.org/ref/epsg/4326/ogcwkt/) es (su segunda proyección):

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

y no

GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]

que es el "ESRI WKT" dialecto representación que usted puede encontrar en http://spatialreference.org/ref/epsg/4326/esriwkt/

Confirmación con gdalsrsinfo:

$ gdalsrsinfo  -o wkt_all  "EPSG:4326"

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]]

....

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

Hay otros dialectos de WKT. El más común dialecto es este "ESRI WKT" dialecto. Pero para las proyecciones de otros que WGS84, ésta no incluye la TOWGS84 y la AUTORIDAD de los nodos de parámetros:

$ gdalsrsinfo  -o wkt_all  "EPSG:4632"

 OGC WKT :
 GEOGCS["Combani 1950",
     DATUM["Combani_1950",
         SPHEROID["International 1924",6378388,297,
             AUTHORITY["EPSG","7022"]],
         TOWGS84[-382,-59,-262,0,0,0,0],
         AUTHORITY["EPSG","6632"]],
     PRIMEM["Greenwich",0,
         AUTHORITY["EPSG","8901"]],
     UNIT["degree",0.0174532925199433,
         AUTHORITY["EPSG","9122"]],
     AUTHORITY["EPSG","4632"]]

 OGC WKT (simple) :
 GEOGCS["Combani 1950",
     DATUM["Combani_1950",
         SPHEROID["International 1924",6378388,297],
         TOWGS84[-382,-59,-262,0,0,0,0]],
     PRIMEM["Greenwich",0],
     UNIT["degree",0.0174532925199433]]

 .....

 ESRI WKT :
 GEOGCS["Combani 1950",
     DATUM["D_Combani_1950",
         SPHEROID["International_1924",6378388,297]],
     PRIMEM["Greenwich",0],
     UNIT["Degree",0.017453292519943295]]

lo que crea problemas para las aplicaciones o bibliotecas que no tienen ESRI WKT en cuenta (como libLAS) (los basados en GDAL/OGR (QGIS, GRASS GIS, PostGIS ...), potencialmente, sí).