[gdal-dev] WMTS services error

Rahkonen Jukka (MML) jukka.rahkonen at maanmittauslaitos.fi
Fri Feb 12 10:48:46 PST 2021


Hi,

That’s not the main issue because hand written XML gives the same error. I guess that GDAL is reading the EPSG:4326 tile matrix set from GetCapabilities and does not like it. This is what I used for testing:

<GDAL_WMTS>
  <GetCapabilitiesUrl>http://www.ign.es/wmts/pnoa-ma?request=GetCapabilities&service=WMTS&version=1.0.0</GetCapabilitiesUrl>
  <Layer>OI.OrthoimageCoverage</Layer>
  <Style>default</Style>
  <TileMatrixSet>EPSG:4326</TileMatrixSet>
  <DataWindow>
    <UpperLeftX>-180</UpperLeftX>
    <UpperLeftY>90</UpperLeftY>
    <LowerRightX>180</LowerRightX>
    <LowerRightY>-90</LowerRightY>
  </DataWindow>
  <BandsCount>3</BandsCount>
  <Cache />
  <UnsafeSSL>true</UnsafeSSL>
  <ZeroBlockHttpCodes>404</ZeroBlockHttpCodes>
  <ZeroBlockOnServerException>true</ZeroBlockOnServerException>
</GDAL_WMTS>

-Jukka-

Lähettäjä: Travis Kirstine <traviskirstine at gmail.com>
Lähetetty: perjantai 12. helmikuuta 2021 20.43
Vastaanottaja: Rahkonen Jukka (MML) <jukka.rahkonen at maanmittauslaitos.fi>
Kopio: gdal dev <gdal-dev at lists.osgeo.org>
Aihe: Re: [gdal-dev] WMTS services error

Jukka,

I think the main issue is that the gdal_translate command throws an error when generating the WMTS XML config file, so there is nothing to edit.

gdal_translate "WMTS:http://www.ign.es/wmts/pnoa-ma?request=GetCapabilities&service=WMTS,layer=OI.OrthoimageCoverage,tilematrixset=EPSG:4326" wmts_bm.xml -of WMTS
ERROR 1: Invalid dataset dimensions : -135056383 x -167901439

On Fri, 12 Feb 2021 at 10:32, jratike80 <jukka.rahkonen at maanmittauslaitos.fi<mailto:jukka.rahkonen at maanmittauslaitos.fi>> wrote:
Hi,

By reading the documentation https://gdal.org/drivers/raster/wmts.html it
could be possible to edit the automatically created XML service definition
file or write it by hand from a scratch.

-Jukka Rahkonen-


Travis Kirstine wrote
> 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 <

> eruiz@

> > 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
>>

> eruiz@

>> *www.aplitop.com<http://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 <

> traviskirstine@

> >
>> *Enviado el:* jueves, 11 de febrero de 2021 20:20
>> *Para:* Elena Ruiz <

> eruiz@

> >
>> *CC:*

> gdal-dev at .osgeo<mailto:gdal-dev at .osgeo>

>> *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>
>>

>>
> <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 <

> eruiz@

> > 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
>>

> eruiz@

>> Sumatra, 9 - 29190 Málaga (Spain)
>> *www.aplitop.com<http://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
>>

> gdal-dev at .osgeo<mailto:gdal-dev at .osgeo>

>> https://lists.osgeo.org/mailman/listinfo/gdal-dev
>>
>>
>
> _______________________________________________
> gdal-dev mailing list

> gdal-dev at .osgeo<mailto:gdal-dev at .osgeo>

> https://lists.osgeo.org/mailman/listinfo/gdal-dev
>
>
> image002.jpg (984 bytes)
> <http://osgeo-org.1560.x6.nabble.com/attachment/5457213/0/image002.jpg>
> image004.jpg (980 bytes)
> <http://osgeo-org.1560.x6.nabble.com/attachment/5457213/1/image004.jpg>
> image006.jpg (1008 bytes)
> <http://osgeo-org.1560.x6.nabble.com/attachment/5457213/2/image006.jpg>
> image008.png (918 bytes)
> <http://osgeo-org.1560.x6.nabble.com/attachment/5457213/3/image008.png>





--
Sent from: http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html
_______________________________________________
gdal-dev mailing list
gdal-dev at lists.osgeo.org<mailto:gdal-dev at lists.osgeo.org>
https://lists.osgeo.org/mailman/listinfo/gdal-dev
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20210212/0e173e37/attachment-0001.html>


More information about the gdal-dev mailing list