Application: gvSIG desktop

It seems that also other prj files are not recognized.

Update (build 2414)

1. its geometries are not correctly reprojected on the fly (see #2739);
2. EPSG:3857 (but also other CRSs) is not identified in the Add layer window, near the "(reprojected on the fly" text, because the following error occurs
   

   ERROR AWT-EventQueue-0 org.gvsig.crs.Crs - Can't autoidentifyEPSG.
   java.lang.RuntimeException: OGR Error: Unsupported SRS
       at org.gdal.osr.osrJNI.SpatialReference_AutoIdentifyEPSG(Native Method)
       at org.gdal.osr.SpatialReference.AutoIdentifyEPSG(SpatialReference.java:247)
       at org.gvsig.crs.Crs.initFromSpatialReference(Crs.java:113)
       at org.gvsig.crs.Crs.<init>(Crs.java:100)
       [...]
   

Hi all,
gdalsrsinfo could be our friend here... A great advantage of using GDAL/OGR libraries is that we can do all it does. So if a projection file (usually ESRI WKT, sometimes OGC WKT) is parsed by gdalsrsinfo correctly, it should be readable by gvSIG. But the AUTHORITY tag is available only in the full OGC WKT.

For instance:

gdalsrsinfo -o wkt_all non_working_3857.prj

OGC WKT :
PROJCS["WGS 84 / Pseudo-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["Geodetic longitude",EAST],
        AXIS["Geodetic latitude",NORTH],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Popular Visualisation Pseudo Mercator",
        AUTHORITY["EPSG","1024"]],
    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["Easting",EAST],
    AXIS["Northing",NORTH],
    AUTHORITY["EPSG","3857"]]

OGC WKT (simple) :
PROJCS["WGS 84 / Pseudo-Mercator",
    GEOGCS["WGS 84",
        DATUM["World Geodetic System 1984",
            SPHEROID["WGS 84",6378137.0,298.257223563]],
        PRIMEM["Greenwich",0.0],
        UNIT["degree",0.017453292519943295]],
    PROJECTION["Popular Visualisation Pseudo Mercator"],
    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]]

OGC WKT (no CT) :
PROJCS["WGS 84 / Pseudo-Mercator",
    GEOGCS["WGS 84",
        DATUM["World Geodetic System 1984",
            SPHEROID["WGS 84",6378137.0,298.257223563]],
        PRIMEM["Greenwich",0.0],
        UNIT["degree",0.017453292519943295]],
    PROJECTION["Popular Visualisation Pseudo Mercator"],
    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]]

ESRI WKT :
PROJCS["WGS_84_Pseudo_Mercator",
    GEOGCS["GCS_WGS_1984",
        DATUM["D_World Geodetic System 1984",
            SPHEROID["WGS_1984",6378137.0,298.257223563]],
        PRIMEM["Greenwich",0.0],
        UNIT["Degree",0.017453292519943295]],
    PROJECTION["Popular Visualisation Pseudo Mercator"],
    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]]

Francisco Diaz wrote:

Gracias Antonio,
pero el problema no está cuando viene un prj cuyo contenido es OGC-WKT,
el problema está cuando en el prj está el ESRI-WKT.

Hi Francisco,
I'm agree with you... OGC WKT is easy to identify thanks to the AUTHORITY tag.

El archivo que adjuntó Andrea era del primer tipo, OGC-WKT, pero era otro el problema que hacía que no se cargase correctamente (el parseo de un número creo recordar).

Ok, the prj attached by Andrea (OGC WKT) is well recognized as EPSG:3857 in build 2417, but not reprojected on the fly because gvSIG isn't able to get the proj4 string from EPSG:3857 actually, so it's still not working.

Adjunto un archivo de este tipo, ESRI_3857.prj con el que podemos comprobar que gdalsrsinfo no es capaz de descubrir la AUTHORITY.

This is true for EPSG:3857 (ESRI WKT), because EPSG:3857 and EPSG:3395 (WGS 84 / World Mercator) have the same proj4 definition and their ESRI WKTs differ only for the name of the PROJCS. So it's not possible to identify the CRS univocally. Instead, EPSG:23030 (ESRI WKT) should be identified by gdalsrsinfo without problems.

Es con este tipo de prj con los que tenemos problemas para identificar la AUTHORITY y, por tanto, hemos decidido que si no podemos hacerlo, mejor no hacer caso al archivo para asignar la proyección.

It's a first approximation that solves the issue only in part. The problem is to convert EPSG:3857 (and also the other CRSs) in proj4 strings now and gdalsrsinfo is able to do it, without reinventing the wheel. ;-)

Antonio Falciano wrote:

This is true for EPSG:3857 (ESRI WKT), because EPSG:3857 and EPSG:3395 (WGS 84 / World Mercator) have the same proj4 definition and their ESRI WKTs differ only for the name of the PROJCS. So it's not possible to identify the CRS univocally. Instead, EPSG:23030 (ESRI WKT) should be identified by gdalsrsinfo without problems.

EPSG:3857 and EPSG:3395 have different proj4 definition, sorry for the mistakes

Related to #4288 and #4290.