Problemas de definición del registro CRS (EPSG:3035)

Question

Estoy realizando un análisis sobre uno de los mapas ráster de riesgo de inundación que se encuentran en https://data.jrc.ec.europa.eu/dataset/1d128b6c-a4ee-4858-9e34-6210707f3c81#dataaccess .

Cuando intento recortar los datos en Austria a través de "recortar para ver", obtengo la siguiente advertencia:

Advertencia 1: La definición del CRS proyectado EPSG:3035 obtenida de las claves GeoTIFF no es la misma que la del registro EPSG, lo que puede causar problemas durante las operaciones de reproyección. Establezca la opción de configuración GTIFF_SRS_SOURCE a EPSG para utilizar los parámetros oficiales (anulando los de las claves GeoTIFF), o a GEOKEYS para utilizar los valores personalizados de las claves GeoTIFF y abandonar el código EPSG.

Como soy nuevo en QGIS, no entiendo muy bien qué pasa y dónde se origina este problema. Una búsqueda en Internet no ha servido de nada. ¿Debo tomar en serio la advertencia? ¿O puedo esperar que todo funcione en la línea descendente si ignoro la advertencia? Si debo solucionarlo, ¿cómo puedo hacerlo?

Answer 1

Mirando floodMap_RP200 de su enlace, el problema está en los datos.

gdalinfo floodmap_EFAS_RP200_C.tif

Reports that the data is in EPSG:3035 like:

Driver: GTiff/GeoTIFF
Files: floodmap_EFAS_RP200_C.tif
       floodmap_EFAS_RP200_C.tif.ovr
       floodmap_EFAS_RP200_C.tif.aux.xml
Size is 63976, 45242
Coordinate System is:
PROJCRS["ETRS_1989_LAEA",
    BASEGEOGCRS["ETRS89",
        DATUM["European Terrestrial Reference System 1989",
            ELLIPSOID["GRS 1980",6378137,298.257222101004,
                LENGTHUNIT["metre",1]]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["degree",0.0174532925199433]],
        ID["EPSG",4258]],
    CONVERSION["Lambert Azimuthal Equal Area",
        METHOD["Lambert Azimuthal Equal Area",
            ID["EPSG",9820]],
        PARAMETER["Latitude of natural origin",52,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",10,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["False easting",4321000,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",3210000,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["(E)",east,
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["(N)",north,
            ORDER[2],
            LENGTHUNIT["metre",1]],
    ID["EPSG",3035]]
Data axis to CRS axis mapping: 1,2

pero el registro EPSG nos dice que CRS 3035 tiene ejes de Northing, y luego Easting, por lo que si los datos están en orden de easting/northing (como se informa) no están de hecho en EPSG:3035.

Por lo tanto, debe

Establezca GEOKEYS para utilizar los valores personalizados de las claves GeoTIFF y elimine el código EPSG.

Answer 2

Lo que se almacena en las etiquetas GeoTIFF se puede listar con gdalinfo:

gdalinfo floodmap_EFAS_RP010_C.tif
...
Coordinate System is:
PROJCRS["ETRS_1989_LAEA",
    BASEGEOGCRS["ETRS89",
        DATUM["European Terrestrial Reference System 1989",
            ELLIPSOID["GRS 1980",6378137,298.257222101004,
                LENGTHUNIT["metre",1]]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["degree",0.0174532925199433]],
        ID["EPSG",4258]],
    CONVERSION["Lambert Azimuthal Equal Area",
        METHOD["Lambert Azimuthal Equal Area",
            ID["EPSG",9820]],
        PARAMETER["Latitude of natural origin",52,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",10,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["False easting",4321000,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",3210000,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["(E)",east,
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["(N)",north,
            ORDER[2],
            LENGTHUNIT["metre",1]],
    ID["EPSG",3035]]

La definición actual de EPSG:3035 puede consultarse en el registro de EPSG https://epsg.org/crs/wkt/id/3035 y también se puede imprimir con gdalsrsinfo.

gdalsrsinfo epsg:3035

PROJCRS["ETRS89-extended / LAEA Europe",
    BASEGEOGCRS["ETRS89",
        DATUM["European Terrestrial Reference System 1989",
            ELLIPSOID["GRS 1980",6378137,298.257222101,
                LENGTHUNIT["metre",1]]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["degree",0.0174532925199433]],
        ID["EPSG",4258]],
    CONVERSION["Europe Equal Area 2001",
        METHOD["Lambert Azimuthal Equal Area",
            ID["EPSG",9820]],
        PARAMETER["Latitude of natural origin",52,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",10,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["False easting",4321000,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",3210000,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["northing (Y)",north,
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["easting (X)",east,
            ORDER[2],
            LENGTHUNIT["metre",1]],
    ID["EPSG",3035]]

Puedes ver que hay algunas diferencias:

PROJCRS["ETRS_1989_LAEA" vs. PROJCRS["ETRS89-extended / LAEA Europe", CONVERSION["Lambert Azimuthal Equal Area" vs. CONVERSION["Europe Equal Area 2001"

y

AXIS["(E)",east,
             ORDER[1],
             LENGTHUNIT["metre",1]],
         AXIS["(N)",north,
             ORDER[2],

contra.

AXIS["northing (Y)",north,
             ORDER[1],
             LENGTHUNIT["metre",1]],
         AXIS["easting (X)",east,
             ORDER[2],

Si quiere estar seguro de que la imagen recortada coincide correctamente con la original ahora y en el futuro, utilice la opción "GEOKEYS para utilizar valores personalizados de claves GeoTIFF". O si encuentra que la diferencia en la georreferenciación es pequeña, puede considerar actualizar el CRS de la imagen original con gdal_edit https://gdal.org/programs/gdal_edit.html -a_srs epsg:3035 .