[gdal-dev] WMTS services error
jratike80
jukka.rahkonen at maanmittauslaitos.fi
Mon Feb 15 04:45:53 PST 2021
Some more information into this issue.
The server at IGN is obviously a GeoServer. GDAL usually does not have any
trouble with GeoServers. Try for example
gdal_translate
"http://demo.geo-solutions.it/geoserver/gwc/service/wmts?REQUEST=GetCapabilities"
-of WMTS geoserver.xml
No errors with that URL. Conclusion: there must be some difference. By
comparing the GetCapabilities responses there seems to be one here:
GeoServer at Geosolutions:
<TileMatrixSet>
<ows:Identifier>EPSG:4326</ows:Identifier>
<ows:SupportedCRS>urn:ogc:def:crs:EPSG::4326</ows:SupportedCRS>
<TileMatrix>
<ows:Identifier>EPSG:4326:0</ows:Identifier>
<ScaleDenominator>2.795411320143589E8</ScaleDenominator>
<TopLeftCorner>90.0 -180.0</TopLeftCorner>
GeoServer at IGN:
<TileMatrixSet>
<ows:Identifier>EPSG:4326</ows:Identifier>
<ows:SupportedCRS>EPSG:4326</ows:SupportedCRS>
#commented out! -- 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>
IGN has changed the SupportedCRS from urn:ogc:def:crs:EPSG::4326 (that works
with GDAL) into EPSG:4326 (that does not work with GDAL).
Travis did find this note from the OGC standard:
"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"
Ok, GeoServer seems to one of those implementers. Now there is a question
what was the alternative. In the standard the CRS used in "World CRS84 Quad
TileMatrixSet definition" is
CRS: http://www.opengis.net/def/crs/OGC/1.3/CRS84, CRS84
Final conclusion:
OGC has defined http://www.opengis.net/def/crs/OGC/1.3/CRS84 for the
longitude-latitude order and http://www.opengis.net/def/crs/EPSG/0/4326 for
the latitude-longitude order. For some reason IGN wants to use a third
alternative EPSG:4326, also for latitude-longitude order. This is something
that is undefined in the standard. I would say that the bug is on the IGN
side and I would like to know what problem they have had with
http://www.opengis.net/def/crs/EPSG/0/4326.
What GDAL might do better is to have an option to flip the axes with the
WMTS XML definition file.
-Jukka Rahkonen-
jratike80 wrote
> 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>
>
>
> <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@
> >
> Lähetetty: perjantai 12. helmikuuta 2021 20.43
> Vastaanottaja: Rahkonen Jukka (MML) <
> jukka.rahkonen@
> >
> Kopio: gdal dev <
> gdal-dev at .osgeo
> >
> 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@
> <mailto:
> jukka.rahkonen@
> >> 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 .osgeo
> <mailto:
> gdal-dev at .osgeo
> >
> https://lists.osgeo.org/mailman/listinfo/gdal-dev
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at .osgeo
> https://lists.osgeo.org/mailman/listinfo/gdal-dev
--
Sent from: http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html
More information about the gdal-dev
mailing list