Elena,

the difference is likely due to the PROJ version used. I'd assume in the first case to be PROJ < 4.9.3 and in the second PROJ >= 4.9.3

You actually get a better result in the first case. If you add -oo EXTENT_METHOD=MOST_PRECISE_TILE_MATRIX, you should get the best result in any case

I've filed an issue https://github.com/OSGeo/gdal/issues/4461 about that

Even

Le 09/09/2021 à 18:03, Elena Ruiz a écrit :

Hello, I have two versions of GDAL 2.4.4 with the dll gdal204.dll dated 2020/01/08, but one is 17.4Mb and the other 20.3 Mb, with Gdal.VersionInfo I get the version of both files and is the same, but they give different results with the following command:

DLL gdal204.dll (20,3 MB)

>*gdalinfo "WMTS:http://www.ign.es/wmts/pnoa-ma?request=GetCapabilities&service=WMTS&username=&password=,layer=OI.OrthoimageCoverage,tilematrixset=EPSG:25830"*

ERROR 1: latitude or longitude exceeded limits

Driver: WMTS/OGC Web Map Tile Service

Files: none associated

Size is 20616192, 28663296

Coordinate System is:

PROJCS["ETRS89 / UTM zone 30N",

GEOGCS["ETRS89",

DATUM["European_Terrestrial_Reference_System_1989",

SPHEROID["GRS 1980",6378137,298.257222101,

AUTHORITY["EPSG","7019"]],

TOWGS84[0,0,0,0,0,0,0],

AUTHORITY["EPSG","6258"]],

PRIMEM["Greenwich",0,

AUTHORITY["EPSG","8901"]],

UNIT["degree",0.0174532925199433,

AUTHORITY["EPSG","9122"]],

AUTHORITY["EPSG","4258"]],

PROJECTION["Transverse_Mercator"],

PARAMETER["latitude_of_origin",0],

PARAMETER["central_meridian",-3],

PARAMETER["scale_factor",0.9996],

PARAMETER["false_easting",500000],

PARAMETER["false_northing",0],

UNIT["metre",1,

AUTHORITY["EPSG","9001"]],

AXIS["Easting",EAST],

AXIS["Northing",NORTH],

AUTHORITY["EPSG","25830"]]

Origin = (-1968157.094999999972060,7271532.000000000000000)

Pixel Size = (0.149291070899963,-0.149291070899963)

Metadata:

ABSTRACT=Imagen de sat├®lite Sentinel2 a escalas menores de 1:70.000 y las ortofotograf├¡as PNOA de m├íxima actualidad para escalas mayores, para toda Espa├▒a. Las teselas se pregeneran en formato JPEG y hasta el nivel 19 (correspondiente con una escala aproximada 1:1.000) en el Sistema de Referencia por Coordenadas WGS84 Web Mercator (EPSG:3857).

  TITLE=Im├ígenes de sat├®lite Sentinel y ortofotos PNOA

Image Structure Metadata:

INTERLEAVE=PIXEL

Corner Coordinates:

Upper Left  (-1968157.095, 7271532.000) ( 45d 9'51.15"W, 57d59'16.66"N)

Lower Left  (-1968157.095, 2992357.845) ( 27d 0' 0.04"W, 24d59'58.98"N)

Upper Right ( 1109656.287, 7271532.000) ( 10d 0' 2.21"E, 65d 0' 0.77"N)

Lower Right ( 1109656.287, 2992357.845) (  3d 8' 0.50"E, 26d55'11.11"N)

Center      ( -429250.404, 5131944.922) ( 14d56'30.30"W, 45d42'45.53"N)

DLL gdal204.dll (17,4 MB)

>*gdalinfo "WMTS:http://www.ign.es/wmts/pnoa-ma?request=GetCapabilities&service=WMTS&username=&password=,layer=OI.OrthoimageCoverage,tilematrixset=EPSG:25830"*

Driver: WMTS/OGC Web Map Tile Service

Files: none associated

Size is 432948, 28663296

Coordinate System is:

PROJCS["ETRS89 / UTM zone 30N",

GEOGCS["ETRS89",

DATUM["European_Terrestrial_Reference_System_1989",

SPHEROID["GRS 1980",6378137,298.257222101,

AUTHORITY["EPSG","7019"]],

TOWGS84[0,0,0,0,0,0,0],

AUTHORITY["EPSG","6258"]],

PRIMEM["Greenwich",0,

AUTHORITY["EPSG","8901"]],

UNIT["degree",0.0174532925199433,

     AUTHORITY["EPSG","9122"]],

AUTHORITY["EPSG","4258"]],

PROJECTION["Transverse_Mercator"],

PARAMETER["latitude_of_origin",0],

PARAMETER["central_meridian",-3],

PARAMETER["scale_factor",0.9996],

PARAMETER["false_easting",500000],

PARAMETER["false_northing",0],

UNIT["metre",1,

AUTHORITY["EPSG","9001"]],

AXIS["Easting",EAST],

AXIS["Northing",NORTH],

AUTHORITY["EPSG","25830"]]

Origin = (404318.828097780002281,7271532.000000000000000)

Pixel Size = (0.149291070899963,-0.149291070899963)

Metadata:

ABSTRACT=Imagen de sat├®lite Sentinel2 a escalas menores de 1:70.000 y las ortofotograf├¡as PNOA de m├íxima actualidad para escalas mayores, para toda Espa├▒a. Las teselas se pregeneran en formato JPEG y hasta el nivel 19 (correspondiente con una escala aproximada 1:1.000) en el Sistema de Referencia por Coordenadas WGS84 Web Mercator (EPSG:3857).

  TITLE=Im├ígenes de sat├®lite Sentinel y ortofotos PNOA

Image Structure Metadata:

INTERLEAVE=PIXEL

Corner Coordinates:

Upper Left  (  404318.828, 7271532.000) (  5d 4'19.71"W, 65d33' 6.57"N)

Lower Left  (  404318.828, 2992357.845) ( 3d57'53.21"W, 27d 3' 0.59"N)

Upper Right (  468954.099, 7271532.000) ( 3d40'21.33"W, 65d33'52.09"N)

Lower Right (  468954.099, 2992357.845) ( 3d18'47.02"W, 27d 3'11.25"N)

Center      (  436636.463, 5131944.922) ( 3d49'24.02"W, 46d20'17.22"N)

As you can see, the command is exactly the same but in the first case an ERROR1 appears and also the coordinates are not the same as in the second case. What could the error be due to? Thanks and regards

------------------------------------------------------------------------

        
*Elena Ruiz *
Sofware Development & Technical Support
Tel. +34 952 43 97 71
[email protected]
Sumatra, 9 - 29190 Málaga (Spain)
*www.aplitop.com

<https://www.facebook.com/AplitopSL/> <https://twitter.com/aplitop> <https://www.linkedin.com/company/aplitop> <https://www.youtube.com/user/aplitop>*

In accordance with the provisions of the European Regulation of Data Protection 2016/679 (Reglamento Europeo de Protección de Datos 2016/679), we inform you that the data and the information you provide us through this medium will be used by APLITOP, S.L., with C.I.F. B-92543396 and with address at C / Sumatra, 9, Malaga, 29190, in order to answer your questions and inform you about our products. The data provided will be kept as long as it does not request its cessation and will not be transferred to third parties except in cases where there is a legal obligation. You have the right to access your personal data, correct inaccurate data or request its deletion when the data is no longer necessary for the purposes that were collected, as well as any rights recognized in the RGPD 2016/679




_______________________________________________
gdal-dev mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/gdal-dev

--
http://www.spatialys.com
My software is free, but my time generally not.

_______________________________________________
gdal-dev mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to