Elena This appears to be a bug in wmts driver in gdal, I would suggest filing a issue
The below is from the OGC Two Dimensional Tile Matrix Set specs, section d2 http://docs.opengeospatial.org/is/17-083r2/17-083r2.html "NOTE4: Some implementers prefer to define this TileMatrixSet using the CRS http://www.opengis.net/def/crs/EPSG/0/4326. The definition is the same except that CRS coordinates are expressed in latitude, longitude order, affecting the TopLeftCorner and the BBox encoding only". On Fri, 12 Feb 2021 at 08:05, Elena Ruiz <[email protected]> wrote: > Hello, my goal is to obtain images of both WMS and WMTS web services using > GDAL, with the CRS that each service provides, in the case of the example, > how can I do then to obtain images in EPSG: 4326 or any other CRS that have > changed axes? Should a parameter be added to the gdal_translate call? What > would be the correct way to obtain this type of images using gdal > applications like gadlinfo and gdaltranslate? > > I have found very little information about this on the web, many thanks > and best regards > > > > > > ------------------------------ > > *Elena Ruiz* > > > Sofware Development & Technical Support > Tel.: +34 952 43 97 71 > [email protected] > *www.aplitop.com <http://www.aplitop.com>* > > > > <http://www.aplitop.com/Product/es/6/5/tcptunnel> > > <http://www.aplitop.com/New/es/281/mdt-75-a-la-venta> > > > > > > *De:* Travis Kirstine <[email protected]> > *Enviado el:* jueves, 11 de febrero de 2021 20:20 > *Para:* Elena Ruiz <[email protected]> > *CC:* [email protected] > *Asunto:* Re: [gdal-dev] WMTS services error > > > > Elena, > > > > I'm not exactly sure what is causing the issue but if you look at the > capabilities of the 4326 tilematrixset compared to the InspireCRS84Quad the > topleftcorner coordinates are reversed. My guess gdal is expecting a X Y > order. > > > > <TileMatrixSet> > <ows:Identifier>InspireCRS84Quad</ows:Identifier> > <ows:SupportedCRS>http://www.opengis.net/def/crs/OGC/1.3/CRS84 > </ows:SupportedCRS> > <TileMatrix> > <ows:Identifier>0</ows:Identifier> > <ScaleDenominator>2.79541132014358E8</ScaleDenominator> > <TopLeftCorner>-180.0 90.0</TopLeftCorner> > <TileWidth>256</TileWidth> > <TileHeight>256</TileHeight> > <MatrixWidth>2</MatrixWidth> > <MatrixHeight>1</MatrixHeight> > </TileMatrix> > > > > <TileMatrixSet> > <ows:Identifier>EPSG:4326</ows:Identifier> > <ows:SupportedCRS>EPSG:4326</ows:SupportedCRS> > <!-- crs="urn:ogc:def:crs:EPSG::4326" --> > <TileMatrix> > <ows:Identifier>EPSG:4326:0</ows:Identifier> > <ScaleDenominator>2.795411320143589E8</ScaleDenominator> > <TopLeftCorner>90.0 -180.0</TopLeftCorner> > <TileWidth>256</TileWidth> > <TileHeight>256</TileHeight> > <MatrixWidth>2</MatrixWidth> > <MatrixHeight>1</MatrixHeight> > </TileMatrix> > > > > > > On Thu, 11 Feb 2021 at 11:51, Elena Ruiz <[email protected]> wrote: > > Hello, I´m trying call web services (WMTS and WMS) using GDAL 3.0.2, but I > have a problem with this example: > > > > gdal_translate "WMTS: > http://www.ign.es/wmts/pnoa-ma?request=GetCapabilities&service=WMTS,layer=OI.OrthoimageCoverage,tilematrixset=EPSG:4326" > d://wmts_bm.xml -of WMTS > > ERROR 1: Invalid dataset dimensions : -135056383 x -167901439 > > > > gdal_translate "WMTS: > http://www.ign.es/wmts/pnoa-ma?request=GetCapabilities&service=WMTS,layer=OI.OrthoimageCoverage,tilematrixset=EPSG:25830" > d://wmts_bm.xml -of WMTS > > Input file size is 20616192, 26817505 > > > > In this server http://www.ign.es/wmts/pnoa-ma have available the next SRS > : > > SUBDATASET_1_NAME=WMTS: > http://www.ign.es/wmts/pnoa-ma?request=GetCapabilities&service=WMTS,layer=OI.OrthoimageCoverage,tilematrixset=InspireCRS84Quad > > SUBDATASET_2_NAME=WMTS: > http://www.ign.es/wmts/pnoa-ma?request=GetCapabilities&service=WMTS,layer=OI.OrthoimageCoverage,tilematrixset=EPSG:4326 > > SUBDATASET_3_NAME=WMTS: > http://www.ign.es/wmts/pnoa-ma?request=GetCapabilities&service=WMTS,layer=OI.OrthoimageCoverage,tilematrixset=EPSG:4258 > > SUBDATASET_4_NAME=WMTS: > http://www.ign.es/wmts/pnoa-ma?request=GetCapabilities&service=WMTS,layer=OI.OrthoimageCoverage,tilematrixset=EPSG:25830 > > SUBDATASET_5_NAME=WMTS: > http://www.ign.es/wmts/pnoa-ma?request=GetCapabilities&service=WMTS,layer=OI.OrthoimageCoverage,tilematrixset=GoogleMapsCompatible > > > > But when I want Know the information or download a image in certain > coordinates with EPSG:4326, it returns an error “ERROR 1: Invalid dataset > dimensions : -135056383 x -167901439”, with EPSG:25830 works fine. I don`t > understand this error, someone could help me? > > > > Thank you 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 <http://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 > >
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
