Re: [gdal-dev] gdalwarp across the antimeridian with a source extent >180E

2024-05-31 Thread Michael Sumner via gdal-dev
I forgot to mention that it needs earthdata credentials set up, basically your "Authorization: Bearer " in GDAL_HTTP_HEADERS or similar config. https://urs.earthdata.nasa.gov/documentation/for_users/user_token You can't download or stream these files without that set or logging in (the file used

[gdal-dev] gdalwarp across the antimeridian with a source extent >180E

2024-05-31 Thread Michael Sumner via gdal-dev
This data source has an odd georeferencing, it's a 36000x17999 raster in the extent -179.995 -89.995 180.005 89.995. vrt://NETCDF:/vsicurl/ https://archive.podaac.earthdata.nasa.gov/podaac-ops-cumulus-protected/MUR-JPL-L4-GLOB-v4.1/2024010109-JPL-L4_GHRSST-SSTfnd-MUR-GLOB-v02.0-fv04.1.nc:analy

Re: [gdal-dev] GDAL WMTS get tiles no documentation

2024-05-31 Thread Michał Kowalczuk via gdal-dev
Thank you again. Now I know which direction should I go.This is my goal to have a control over this process using the low-level RasterIO. I browse GDAL github repository a lot and it's very helpful, especially the tests. This time I was confused with WMTS and I needed a support at the begining. Ha

Re: [gdal-dev] GDAL WMTS get tiles no documentation

2024-05-31 Thread Rahkonen Jukka via gdal-dev
Hi, You do not need to use GDALTranslate even if may be convenient. In the RasterIO() function GDALRasterBand C++ API — GDAL documentation the region is selected by using XOff, YOff, XSize and YSize which map directly with -srcwin. But I was thinki

Re: [gdal-dev] GDAL WMTS get tiles no documentation

2024-05-31 Thread Rahkonen Jukka via gdal-dev
(resending, need to cut old discussion due to too many bytes in the body) Hi, gdalinfo "WMTS:http://maps.wien.gv.at/wmts/1.0.0/WMTSCapabilities.xml,layer=fmzk3dviewer"; Driver: WMTS/OGC Web Map Tile Service Files: none associated Size is 305719, 257379 The unit of size is pixels. You used -srcw

Re: [gdal-dev] GDAL WMTS get tiles no documentation

2024-05-31 Thread Michał Kowalczuk via gdal-dev
As you proposed I tested *gdal_translate *utility but for the following sample command I get an empty image (tested with jpg and png) gdal_translate "WMTS: http://maps.wien.gv.at/wmts/1.0.0/WMTSCapabilities.xml,layer=fmzk3dviewer"; -srcwin 1820157 6144777 512 512 tile.jpg [image: image.png] pt.,

Re: [gdal-dev] GDAL WMTS get tiles no documentation

2024-05-31 Thread Michał Kowalczuk via gdal-dev
Thank you Rukka, Since *-srcwin* and *-projwin* are options only for *gdal_translate* utility, your suggestion leads me to use *GDALTranslate()* function with these parameters. This means that I cannot use directly created WMTS dataset to retrieve data (using GDALDatasetRasterIO), but create a new

Re: [gdal-dev] GDAL WMTS get tiles no documentation

2024-05-31 Thread Rahkonen Jukka via gdal-dev
Hi, Before starting to program you can familiar yourself with the topic by using gdal_translate gdal_translate — GDAL documentation. Play with the -srcwin and -projwin / -projwin_srs options and gather some trust that is is possible to select and d

Re: [gdal-dev] GDAL WMTS get tiles no documentation

2024-05-31 Thread Michał Kowalczuk via gdal-dev
Hi, thank you for your replies! As I'm still a newbie especially with WMTS. I need to take a step back to the basics. Let's forget about multithreading for now. - I'm in the stage that I can open a dataset using WMTS service path from data provider - I can list subdatasets and open the se

Re: [gdal-dev] OGR SQL CAST geometry example

2024-05-31 Thread Rahkonen Jukka via gdal-dev
Hi, As far I know the the Spatialite functions https://www.gaia-gis.it/gaia-sins/spatialite-sql-4.2.0.html by design work row by row and one input row returns one output row. Therefore DissolvePoints creates one row with a multipoint geometry, not several rows like ST_Dump does. You can discus