Re: [gdal-dev] WMTS services error

2021-02-15 Thread jratike80
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:


  EPSG:4326
  urn:ogc:def:crs:EPSG::4326
  
EPSG:4326:0
2.795411320143589E8
90.0 -180.0

GeoServer at IGN:


 EPSG:4326
 EPSG:4326
 #commented out! --  crs="urn:ogc:def:crs:EPSG::4326"  -->
 
   EPSG:4326:0
   2.795411320143589E8
   90.0 -180.0

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:
> 
>   
> 
> http://www.ign.es/wmts/pnoa-ma?request=GetCapabilities&service=WMTS&version=1.0.0
> 
>   
> 
> OI.OrthoimageCoverage
> 
>   

>   
> 
> EPSG:4326
> 
>   
> 
> 
> 
> -180
> 
> 
> 
> 90
> 
> 
> 
> 180
> 
> 
> 
> -90
> 
>   
> 
>   
> 
> 3
> 
>   
> 
>   
> 
> true
> 
>   
> 
> 404
> 
>   
> 
> true
> 
> 
> -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@.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 cas

Re: [gdal-dev] WMTS services error

2021-02-12 Thread Rahkonen Jukka (MML)
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:


  
http://www.ign.es/wmts/pnoa-ma?request=GetCapabilities&service=WMTS&version=1.0.0
  OI.OrthoimageCoverage
  default
  EPSG:4326
  
-180
90
180
-90
  
  3
  
  true
  404
  true


-Jukka-

Lähettäjä: Travis Kirstine 
Lähetetty: perjantai 12. helmikuuta 2021 20.43
Vastaanottaja: Rahkonen Jukka (MML) 
Kopio: gdal dev 
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 
mailto:jukka.rahko...@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@.osgeo<mailto:gdal-dev@.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.
>>
>>
>>
>>
> 
>>
> 
> InspireCRS84Quad
> 
>>
> 
> http://www.opengis.net/def/crs/OGC/1.3/CRS84
>>
> 
>>
> 
>>
> 
> 0
> 
>>
> 
> 2.79541132014358E8
> 
>>
> 
> -180.0 90.0
> 
>>
> 
> 256
> 
>>
> 
> 256
> 
>>
> 
> 2
> 
>>
> 
> 1
> 
>>
> 
>>
>>
>>
>>
> 
>>
> 
> EPSG:4326
> 
>>
> 
> EPSG:4326
> 
>>

>>
> 
>>
> 
> EPSG:4326:0
> 
>>
> 
> 2.795411320143589E8
> 
>>
> 
> 90.0 -180.0
> 
>>
> 
> 256
> 
>>
> 
> 256
> 
>>
> 
> 2
> 
>>
> 
> 1
> 
>>
> 
>>
>>
>>
>>
>>
>> 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:
>> 

Re: [gdal-dev] WMTS services error

2021-02-12 Thread Travis Kirstine
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 
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/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@.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.
> >>
> >>
> >>
> >>
> > 
> >>
> > 
> > InspireCRS84Quad
> > 
> >>
> > 
> > http://www.opengis.net/def/crs/OGC/1.3/CRS84
> >>
> > 
> >>
> > 
> >>
> > 
> > 0
> > 
> >>
> > 
> > 2.79541132014358E8
> > 
> >>
> > 
> > -180.0 90.0
> > 
> >>
> > 
> > 256
> > 
> >>
> > 
> > 256
> > 
> >>
> > 
> > 2
> > 
> >>
> > 
> > 1
> > 
> >>
> > 
> >>
> >>
> >>
> >>
> > 
> >>
> > 
> > EPSG:4326
> > 
> >>
> > 
> > EPSG:4326
> > 
> >>
>
> >>
> > 
> >>
> > 
> > EPSG:4326:0
> > 
> >>
> > 
> > 2.795411320143589E8
> > 
> >>
> > 
> > 90.0 -180.0
> > 
> >>
> > 
> > 256
> > 
> >>
> > 
> > 256
> > 
> >>
> > 
> > 2
> > 
> >>
> > 
> > 1
> > 
> >>
> > 
> >>
> >>
> >>
> >>
> >>
> >> On Thu, 11 Feb 2021 a

Re: [gdal-dev] WMTS services error

2021-02-12 Thread jratike80
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/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@.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.
>>
>>
>>
>> 
> 
>> 
> 
> InspireCRS84Quad
> 
>> 
> 
> http://www.opengis.net/def/crs/OGC/1.3/CRS84
>> 
> 
>> 
> 
>> 
> 
> 0
> 
>> 
> 
> 2.79541132014358E8
> 
>> 
> 
> -180.0 90.0
> 
>> 
> 
> 256
> 
>> 
> 
> 256
> 
>> 
> 
> 2
> 
>> 
> 
> 1
> 
>> 
> 
>>
>>
>>
>> 
> 
>> 
> 
> EPSG:4326
> 
>> 
> 
> EPSG:4326
> 
>> 

>> 
> 
>> 
> 
> EPSG:4326:0
> 
>> 
> 
> 2.795411320143589E8
> 
>> 
> 
> 90.0 -180.0
> 
>> 
> 
> 256
> 
>> 
> 
> 256
> 
>> 
> 
> 2
> 
>> 
> 
> 1
> 
>> 
> 
>>
>>
>>
>>
>>
>> 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&

Re: [gdal-dev] WMTS services error

2021-02-12 Thread Travis Kirstine
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  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
> er...@aplitop.com
> *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 
> *Enviado el:* jueves, 11 de febrero de 2021 20:20
> *Para:* Elena Ruiz 
> *CC:* gdal-dev@lists.osgeo.org
> *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.
>
>
>
> 
> InspireCRS84Quad
> http://www.opengis.net/def/crs/OGC/1.3/CRS84
> 
> 
> 0
> 2.79541132014358E8
> -180.0 90.0
> 256
> 256
> 2
> 1
> 
>
>
>
> 
> EPSG:4326
> EPSG:4326
> 
> 
> EPSG:4326:0
> 2.795411320143589E8
> 90.0 -180.0
> 256
> 256
> 2
> 1
> 
>
>
>
>
>
> On Thu, 11 Feb 2021 at 11:51, Elena Ruiz  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
> er...@aplitop.com
> 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),

Re: [gdal-dev] WMTS services error

2021-02-12 Thread Elena Ruiz
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

[https://www.aplitop.com/descargas/logo.png]
Sofware Development & Technical Support
Tel.: +34 952 43 97 71
er...@aplitop.com
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 
Enviado el: jueves, 11 de febrero de 2021 20:20
Para: Elena Ruiz 
CC: gdal-dev@lists.osgeo.org
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.


InspireCRS84Quad
http://www.opengis.net/def/crs/OGC/1.3/CRS84

0
2.79541132014358E8
-180.0 90.0
256
256
2
1



EPSG:4326
EPSG:4326


EPSG:4326:0
2.795411320143589E8
90.0 -180.0
256
256
2
1



On Thu, 11 Feb 2021 at 11:51, Elena Ruiz 
mailto:er...@aplitop.com>> 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


[https://www.aplitop.com/descargas/firma1.jpg]

Elena Ruiz
Sofware Development & Technical Support
Tel. +34 952 43 97 71
er...@aplitop.com<mailto:er...@aplitop.com>
Sumatra, 9 - 29190 Málaga (Spain)
www.aplitop.com<http://www.aplitop.com>

[cid:image002.jpg@01D70148.0DA21C80]<https://www.facebook.com/AplitopSL/>   
[cid:image004.jpg@01D70148.0DA21C80] <https://twitter.com/aplitop>
[cid:image006.jpg@01D70148.0DA21C80] <https://www.linkedin.com/company/aplitop> 
   [cid:image008.png@01D70148.0DA21C80] <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@lists.osgeo.org<mailto:gdal-dev@lists.osgeo.org>
https://lists.osgeo.org/mailman/listinfo/gdal-dev
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev


Re: [gdal-dev] WMTS services error

2021-02-11 Thread Travis Kirstine
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.


InspireCRS84Quad
http://www.opengis.net/def/crs/OGC/1.3/CRS84


0
2.79541132014358E8
-180.0 90.0
256
256
2
1



EPSG:4326
EPSG:4326


EPSG:4326:0
2.795411320143589E8
90.0 -180.0
256
256
2
1



On Thu, 11 Feb 2021 at 11:51, Elena Ruiz  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
> er...@aplitop.com
> Sumatra, 9 - 29190 Málaga (Spain)
>
>
> *www.aplitop.com 
>    
> 
> *
>
> 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@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev
>
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev


[gdal-dev] WMTS services error

2021-02-11 Thread Elena Ruiz
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


[https://www.aplitop.com/descargas/firma1.jpg]


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

[cid:SocialLink_Facebook_32x32_a745adb8-c250-4372-8480-e77f503dad3d.png]
   [cid:SocialLink_Twitter_32x32_858247f6-9d4a-4740-a6b3-9b5ebfd98a97.png] 

[cid:SocialLink_Linkedin_32x32_154d63de-bf68-468f-9eb0-ffed102fac77.png] 

[cid:SocialLink_Youtube_32x32_0dea15dd-3774-4086-be9d-f62560b70c86.png] 



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@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev