¿Cómo encontrar el nombre real de proyección (código) en un dato de texto PROJCS?

Question

Dado un PROJCS de datos

PROJCS["unnamed",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223560493,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",-25.91105555555576],
    PARAMETER["standard_parallel_2",-32.72172222222201],
    PARAMETER["latitude_of_origin",-29.51611111111111],
    PARAMETER["central_meridian",-52.09266666666708],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (-725301.460735953995027,572835.914701772970147)
Pixel Size = (169.335263023642995,-169.335263023642995)

¿Cómo puedo encontrar la proyección actual nombre? Puedo ver EPSG:7030, EPSG:6326, EPSG:4326 y EPSG:9001. ¿Cuál es la real EPSG:xxxx de esta proyección?

En comparación con el conocido EPSG:900913 puedo concluir la última AUTORIDAD es la nomenclatura de la Proyección:

PROJCS["WGS84 / Google Mercator", 
  GEOGCS["WGS 84", 
    DATUM["World Geodetic System 1984", 
      SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]], 
      AUTHORITY["EPSG","6326"]], 
    PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], 
    UNIT["degree", 0.017453292519943295], 
    AXIS["Longitude", EAST], 
    AXIS["Latitude", NORTH], 
    AUTHORITY["EPSG","4326"]], 
  PROJECTION["Mercator_1SP"], 
  PARAMETER["semi_minor", 6378137.0], 
  PARAMETER["latitude_of_origin", 0.0], 
  PARAMETER["central_meridian", 0.0], 
  PARAMETER["scale_factor", 1.0], 
  PARAMETER["false_easting", 0.0], 
  PARAMETER["false_northing", 0.0], 
  UNIT["m", 1.0], 
  AXIS["x", EAST], 
  AXIS["y", NORTH], 
  AUTHORITY["EPSG","900913"]]

Es esto correcto? Así que mi ejemplo será EPSG:9001?

Estoy en troble para crear una capa en mi geoserver utilizando un GeoTIFF con esta proyección. Mientras yo puede obtener una vista previa de la capa fina, no puedo ponerlo a trabajar en un mapa real o en agrupación con otras capas. Siempre ir a 0,0

Answer 1

Como se puede ver a partir de los soportes, la EPSG:9001 sólo es válida para la UNIDAD de la variable.

Por lo general, un distribuidor de datos puede configurar CRS parámetros como él quiere. Algunos de ellos están en la base de datos EPSG (principalmente aquellos que se utilizan para estudiar), y tienen un código EPSG número. Pero tienes que aceptar que un CRS no puede ser representado por un código EPSG.

Si usted trabaja con PROJ.4/GDAL basado en software, puede que desee que el proj.4 definición de la WKT. Usted puede conseguir que si pones el WKT en un archivo de texto, y ejecutar gdalsrsinfo en eso.

Usted recibirá:

 PROJ.4 : '+proj=lcc +lat_1=-25.91105555555576 +lat_2=-32.72172222222201 +lat_0=-29.51611111111111 +lon_0=-52.09266666666708 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs '

Usted puede crear una costumbre de CRS con esa cadena, y trabajar con eso. http://docs.geoserver.org/2.7.1/user/advanced/crshandling/customcrs.html da las instrucciones para que en Geoserver.