Re: [gdal-dev] VRT file: Recursion detected

2023-04-13 Thread Ari Jolma
Even Rouault kirjoitti 13.4.2023 klo 12.44: Le 13/04/2023 à 11:02, Laurențiu Nicola via gdal-dev a écrit : Hi, Just a guess, but maybe your VRT is including itself? That's indeed the purpose of this check (or detecting a VRT A including a VRT B including a VRT A, etc.). But the check is a

[gdal-dev] VRT file: Recursion detected

2023-04-13 Thread Ari Jolma
What can cause "Recursion detected" error when opening a VRT-file? Is there something I can do in gdal.BuildVRT (I'm using gdal python bindings) to avoid that? Ari ___ gdal-dev mailing list gdal-dev@lists.osgeo.org

Re: [gdal-dev] AWS Lambda layer with GDAL python bindings

2023-02-08 Thread Ari Jolma
in 2019 with a few hacks: https://github.com/developmentseed/geolambda *From:*gdal-dev on behalf of Ari Jolma *Date:*Wednesday, February 8, 2023 at 1:42 AM *To:*gdal-dev@lists.osgeo.org *Subject:*[gdal-dev] AWS Lambda layer with GDAL python bindings Hi, I need an AWS Lambda layer with GDAL

Re: [gdal-dev] AWS Lambda layer with GDAL python bindings

2023-02-08 Thread Ari Jolma
access to one AWS account. How do I get those layers? Ari Joe Lee kirjoitti 8.2.2023 klo 16.35: This worked quite well in 2019 with a few hacks: https://github.com/developmentseed/geolambda *From: *gdal-dev on behalf of Ari Jolma *Date: *Wednesday, February 8, 2023 at 1:42 AM *To: *gdal

Re: [gdal-dev] AWS Lambda layer with GDAL python bindings

2023-02-08 Thread Ari Jolma
I'm giving up for now and will move the task into an EC2, where it is really easy to set up the required python environment. Ari Howard Butler kirjoitti 8.2.2023 klo 16.27: On Feb 8, 2023, at 1:42 AM, Ari Jolma wrote: Hi, I need an AWS Lambda layer with GDAL python bindings, for "

[gdal-dev] AWS Lambda layer with GDAL python bindings

2023-02-07 Thread Ari Jolma
Hi, I need an AWS Lambda layer with GDAL python bindings, for "from osgeo import gdal", but it seems difficult to make such a layer. In fact all I need is gdal.Open and gdal.BuildVRT - I have not dug deeper into rasterio, can it do those? Anyway, I installed GDAL on a machine with python

Re: [gdal-dev] gdallocationinfo + WCS question...

2021-10-21 Thread Ari Jolma
I don't understand the problem :) Maybe it's because I ran gdalinfo on the WCS first. For the gdallocationinfo I get gdallocationinfo "WCS:https://elevation.nationalmap.gov:443/arcgis/services/3DEPElevation/ImageServer/WCSServer?version=2.0.1=DEP3Elevation; -wgs84 -122.086 42.948 Report:  

Re: [gdal-dev] OpenFileGDB: Error occurred in filegdbtable.cpp at line 696

2021-10-12 Thread Ari Jolma
The problem was insufficient policy for the IAM of the EC2, it had no rights for the objects in the bucket. Ari Jolma kirjoitti 12.10.2021 klo 12.40: This works: ari@xps:~/lahde$ ogrinfo -ro /vsis3/lahde-test-test/kuviot01102021.gdb INFO: Open of `/vsis3/lahde-test-test/kuviot01102021.gdb

[gdal-dev] OpenFileGDB: Error occurred in filegdbtable.cpp at line 696

2021-10-12 Thread Ari Jolma
This works: ari@xps:~/lahde$ ogrinfo -ro /vsis3/lahde-test-test/kuviot01102021.gdb INFO: Open of `/vsis3/lahde-test-test/kuviot01102021.gdb'   using driver `OpenFileGDB' successful. 1: kuviot (Multi Polygon) 2: invpsto_osite (None) ari@xps:~/lahde$ ogrinfo --version GDAL 3.3.2, released

Re: [gdal-dev] gdal-dev Digest, Vol 203, Issue 38

2021-04-21 Thread Ari Jolma
Elena, It seems that the 2 x 2 image which GDAL first fetches to get information about the data (BBOX=101082,4289059,101102,4289079) triggers an error on the server. I don't know at this point if that is an error on the server side (wrong extent) or on the client side (gdal). Best regards,

Re: [gdal-dev] libgdal depends on two libprojs, what's going on?

2021-03-16 Thread Ari Jolma
Ari Jolma kirjoitti 16.3.2021 klo 16.53: Even Rouault kirjoitti 16.3.2021 klo 16.49: Ari, Maybe it's because libgdal seems to use two libprojs: ldd /usr/lib/libgdal.so.28.0.1 | grep libproj     libproj.so.19 => /usr/lib/x86_64-linux-gnu/libproj.so.19 (0x7f70d3343000)     libproj.so

Re: [gdal-dev] libgdal depends on two libprojs, what's going on?

2021-03-16 Thread Ari Jolma
Even Rouault kirjoitti 16.3.2021 klo 16.49: Ari, Maybe it's because libgdal seems to use two libprojs: ldd /usr/lib/libgdal.so.28.0.1 | grep libproj     libproj.so.19 => /usr/lib/x86_64-linux-gnu/libproj.so.19 (0x7f70d3343000)     libproj.so.15 => /usr/lib/x86_64-linux-gnu/libproj.so.15

[gdal-dev] libgdal depends on two libprojs, what's going on?

2021-03-16 Thread Ari Jolma
I needed later libgdal so I added ubuntugi-unstable and installed libgdal-dev. Now I have both libgdal26 and libgdal28. But I assume something is wrong because I get from trying to use the libgdal28 (gdb) bt #0  __GI_raise (sig=sig@entry=6) at ../sysdeps/unix/sysv/linux/raise.c:50 #1 

Re: [gdal-dev] Wiki page with current maintainers per sub-system

2021-03-06 Thread Ari Jolma
Even Rouault kirjoitti 5.3.2021 klo 20.31: Hi, following recent discussions about maintenance, I've initiated https://github.com/OSGeo/gdal/wiki/Maintainers-per-sub-system with "just" drivers and binding languages You're invited to check (I've put a few names that I had in mind in front of

Re: [gdal-dev] Driver maintenance - long-term solution ?

2021-01-13 Thread Ari Jolma
Thank you Howard, Daniel, Jukka, and others. When I google for qgis I get its page and six sub pages in the way Google does it. One of those is "Donations". When I do the same for GDAL, there's no such page. I guess the first and practical thing to do would be to make regular donations

Re: [gdal-dev] gdal-dev Digest, Vol 198, Issue 14

2020-11-13 Thread Ari Jolma
. But this is a sidepath because I believe that the WCS driver wants to read a small area from WCS as a sample, not just one pixel. -Jukka- Ari Jolma-2 wrote 2) gdalinfo "WCS:https://demo.geo-solutions.it/geoserver/wcs?version=2.0.1=nurc__Img_Sample; .. takes again some time as it makes

Re: [gdal-dev] gdal-dev Digest, Vol 198, Issue 14

2020-11-11 Thread Ari Jolma
developer of the WCS plugin probably could tell if it is at all supported to use the XML file as input like you tried. I wonder if Ari Jolma happens to follow the list. And what QGIS does is all mystery for me. -Jukka Rahkonen- *Lähettäjä:* 1520 gis *Lähetetty:* keskiviikko 11. marraskuuta

Re: [gdal-dev] Retrieving 3DEP data via GDAL's WCS driver...

2020-11-02 Thread Ari Jolma
The WCS driver stores data in folder $HOME/.gdal/wcs_cache. The folder contains a text file db which contains keys and URLs and responses from the server // key.xml = service file // key.xml.aux.xml = metadata file // key.xml = Capabilities response // key.aux.xml = Global metadata //

Re: [gdal-dev] Vector Tile Service

2020-05-12 Thread Ari Jolma
Sorry, no. I'm using PostGIS via tilestrata-postgismvt. The code that combines the single-theme tiles it makes is mine, though. The tilestrata-postgismvt works under tilestrata, which works under express. Ari rmaddu kirjoitti 12.5.2020 klo 14.33: Hi, Are you using GDAL binaries for vector

Re: [gdal-dev] Vector Tile Service

2020-05-12 Thread Ari Jolma
I don't know about C#. In JavaScript I have a use case where I combine several single theme vector tiles into one multi theme tile each in its own thread async and then combine them before sending. I use the async.series method for that. Ari rmaddu kirjoitti 12.5.2020 klo 13.28: Hi, We

Re: [gdal-dev] Vector Tile Service

2020-05-12 Thread Ari Jolma
rmaddu, VRT is raster format and MVT is vector format, so they are incompatible without raster to vector transformation. So, I think there's something missing from your problem statement. Best, Ari rmaddu . kirjoitti 12.5.2020 klo 8.14: I need some help on usage of GDAL C# bindings to

Re: [gdal-dev] Algorithm to clip the inner loop of a polygon

2020-04-15 Thread Ari Jolma
Andreas Oxenstierna kirjoitti 15.4.2020 klo 9.03: Hi Have you tested PostGIS? ST_ExteriorRing((ST_Dump(ST_LineMerge(ST_Union(ST_Force2D(geometry>).geom) should extract the points you want to keep by splitting at intersections - but it may fail in very complex cases. Then I would start

Re: [gdal-dev] Algorithm to clip the inner loop of a polygon

2020-04-14 Thread Ari Jolma
Hi, I wrote a method to a Polygon class (in Python) that solves that problem because the ones in GEOS via Shapely did not seem good enough. My method was based on no published algorithm but just my own reasoning, it is probably not optimal speed-wise. I did not perform formal testing on it

Re: [gdal-dev] GPKG to PostgreSQL, version 2.2.3 vs version 3.0.2

2020-04-10 Thread Ari Jolma
Ari Jolma kirjoitti 10.4.2020 klo 9.54: This seems to be the case. v 3.0.2 sends SELECT setval(pg_get_serial_sequence('"stand"', 'id'), MAX("id") FROM "stand"; and similar for other tables but v 2.2.3 does not. The fix was committed 9 Sep 2017 as a fix to

Re: [gdal-dev] GPKG to PostgreSQL, version 2.2.3 vs version 3.0.2

2020-04-10 Thread Ari Jolma
Ari Jolma kirjoitti 9.4.2020 klo 17.23: Same outcome. There's something very strage going on as the max of the id is 305926 (after creating the table with the first gpkg) but when I try to manually insert a row to the table it says "key id=39 already exists". If I then manua

Re: [gdal-dev] GPKG to PostgreSQL, version 2.2.3 vs version 3.0.2

2020-04-09 Thread Ari Jolma
the table. So somehow the sequence is lagging behind. Ari jratike80 kirjoitti 9.4.2020 klo 16.14: Hi, I would have a try with -unsetfid. That should be the default when in append mode but the test is easy. -Jukka Rahkonen- Ari Jolma-2 wrote I have this command ogr2ogr -update -append -f

[gdal-dev] GPKG to PostgreSQL, version 2.2.3 vs version 3.0.2

2020-04-09 Thread Ari Jolma
I have this command ogr2ogr -update -append -f PostgreSQL PG:... merged.gpkg which works ok in version 3.0.2 (self-compiled) but not in 2.2.3 (from ubuntu apt). The error is ERROR: duplicate key value violates unique constraint "stand_pkey" DETAIL: Key (id)=(1) already exists. CONTEXT: COPY

Re: [gdal-dev] WCS: disk cache?

2019-03-31 Thread Ari Jolma
Tobias Wendorff kirjoitti 29.3.2019 klo 0.09: Hi there, right now, when using WCS, the image cache goes to /vsimem/wcs/, can I also write the files to disk? GDALDriver objecs have a method CreateCopy just for that purpose. Best regards, Ari Sure, I could play with a local HTTPS proxy,

Re: [gdal-dev] Segfault in gdalrasterize

2019-03-27 Thread Ari Jolma
with 2.4.1 the segfault comes from gvBurnPoint, so that may need similar fix than gvBurnScanline Ari Even Rouault kirjoitti 26.3.2019 klo 20.39: On mardi 26 mars 2019 17:50:19 CET Ari Jolma wrote: We're seeing transient segfaults in gvBurnScanline with a huge in-memory raster created

Re: [gdal-dev] Segfault in gdalrasterize

2019-03-27 Thread Ari Jolma
Thanks Even, Ari Even Rouault kirjoitti 26.3.2019 klo 20.39: On mardi 26 mars 2019 17:50:19 CET Ari Jolma wrote: We're seeing transient segfaults in gvBurnScanline with a huge in-memory raster created by rasterio. The raster size is (67264, 51776). The segfaults happen in cloud servers

[gdal-dev] Segfault in gdalrasterize

2019-03-26 Thread Ari Jolma
We're seeing transient segfaults in gvBurnScanline with a huge in-memory raster created by rasterio. The raster size is (67264, 51776). The segfaults happen in cloud servers, not in laptops or local virtual machines. Debugging is hard because rasterio installs its own libgdal but with gdb

[gdal-dev] gdalinfo of a WMS layer

2019-02-26 Thread Ari Jolma
I'm wondering why gdalinfo "WMS:http://tiles.kartat.kapsi.fi/taustakartta?layers=taustakartta; and gdalinfo "WMS:http://tiles.kartat.kapsi.fi/taustakartta?layers=taustakartta=epsg:3857; both report corner coordinates in degrees although the second one reports EPSG 3857? What am I

Re: [gdal-dev] Problem with Perl module Geo:GDAL: undefined symbol: GTIFNewEx

2019-02-11 Thread Ari Jolma
Hernán De Angelis kirjoitti 11.2.2019 klo 19.01: Dear all (but particularly users and developer of Perl module Geo::GDAL), I have resumed work in a project that has been in standstill for about two years. The workflow involves some Perl scripts that read and process properties of satellite

[gdal-dev] libgdal runpath

2018-11-26 Thread Ari Jolma
I'm building compartmentalized systems on servers. There may be several systems with different setups, thus library management is important. I'm wondering why libgdal looks like this readelf -d parts/gdal/lib/libgdal.so  0x0001 (NEEDED) Shared library: [libgeos_c.so.1]

Re: [gdal-dev] problem with loading libraries

2018-11-22 Thread Ari Jolma
I figured out a solution to the problem. The library is loaded by the postgres process, so it needs to have its library path appended. The simplest way is to start it with command LD_LIBRARY_PATH= pg_ctl -D database start Ari Ari Jolma kirjoitti 22.11.2018 klo 17.06: Sorry

[gdal-dev] problem with loading libraries

2018-11-22 Thread Ari Jolma
Sorry, this is not strictly a GDAL question, but it is related. I have built postgis with raster support, which required gdal. Now, when I CREATE EXTENSION postgis; I get ERROR:  could not load library "/home/ajolma/ods/server/parts/postgresql/lib/postgis-2.5.so": libgeos_c.so.1: cannot

Re: [gdal-dev] WCS: Bug in WCSDataset110::GetCoverageRequest()

2018-10-17 Thread Ari Jolma
t;WCS:http://openwms.statkart.no/skwms1/wcs.dtm?VERSION=1.1.1=land_utm33_10m; is the correct approach and it does not fail. Ari Odd-Ragnar< -Original Message----- From: Ari Jolma Sent: onsdag 10. oktober 2018 15:02 To: gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] WCS: Bug in

Re: [gdal-dev] WCS: Bug in WCSDataset110::GetCoverageRequest()

2018-10-10 Thread Ari Jolma
Hi Odd Ragnar, Without going very deep into this, I would do (and tried) gdalinfo "WCS:https://wms.geonorge.no/skwms1/wcs.hoyde-dtm1_32; This tells me that the preferred version of the service is 2.0.1 and it serves coverage Coverage1. However, if I want to use service version 1.1.1, I do

Re: [gdal-dev] C++ API for WCS: GetCapabilities and the new cache

2018-09-07 Thread Ari Jolma
Odd Ragnar Lydersen kirjoitti 07.09.2018 klo 12:57: Hi, I'm struggling to understand the new setup for WCS. In the documentation it says: - If the URL does not contain a coverage name, the driver attempts to fetch the capabilities document from the server, parse it, and show the resulting

Re: [gdal-dev] GDAL dependencies

2018-05-22 Thread Ari Jolma
Even Rouault kirjoitti 22.05.2018 klo 11:05: On mardi 22 mai 2018 07:59:08 CEST Lars wrote: > Hello, > > Does GDAL have an updated list of required and optional dependencies? Lars, Required: none currently, except a C++ compiler With C++ 11 support. I just run into this as I had a

[gdal-dev] installed header files

2018-05-17 Thread Ari Jolma
In 2.3 the installed header files is limited to a fixed list. That is a good thing. Even better would be to install the header files into a GDAL specific subdir. However, that change broke one thing: the bindings cannot now be built against the public header files. At least cpl_vsi_error.h is

Re: [gdal-dev] [OSGeo/gdal] Fix the datadir in gdal.pc. (#609)

2018-05-12 Thread Ari Jolma
Even Rouault kirjoitti 12.05.2018 klo 11:18: > > Now I did it by creating a local clone of OSGeo/gdal and > editing/committing/pushing its release/2.3 branch. You can use a single clone for all branches and uses git checkout to switch branches. I know that but the issue here was that I

Re: [gdal-dev] [OSGeo/gdal] Fix the datadir in gdal.pc. (#609)

2018-05-12 Thread Ari Jolma
Ari Jolma kirjoitti 12.05.2018 klo 10:18: Even Rouault kirjoitti 11.05.2018 klo 23:19: cherry-pick in release/2.3 How do I merge a PR into another branch than master? I made a try to -- PR #614 -- which was a mistake. Now I did it by creating a local clone of OSGeo/gdal and editing

Re: [gdal-dev] PHP and Ruby bindings gone [was: Anyone interested in resurrecting PHP (7) bindings ?]

2018-05-02 Thread Ari Jolma
Updated. I also added the new FFI based Perl API that I'm working on. Ari Ben Elliston kirjoitti 02.05.2018 klo 12:58: On 02/05/18 19:37, Even Rouault wrote: just to inform that I've acceepted 2 pull requests in master to remove the PHP and Ruby bindings, which haven't seen any activity in

Re: [gdal-dev] #include tiffio.h in gtiff/libgeotiff/

2018-04-12 Thread Ari Jolma
On 12.04.2018 19:25, Even Rouault wrote: On jeudi 12 avril 2018 09:56:35 CEST Ari Jolma wrote: The file frmts/gtiff/libgeotiff/xtiffio.h #include's "tiffio.h", which is in frmts/gtiff/libtiff. However, when compiling in this directory frmts/gtiff/libtiff is not included.

Re: [gdal-dev] #include tiffio.h in gtiff/libgeotiff/

2018-04-12 Thread Ari Jolma
ets of .lo and .o: one in current directory and one in ../../frmts/o? Ari On 12.04.2018 09:56, Ari Jolma wrote: The file frmts/gtiff/libgeotiff/xtiffio.h #include's "tiffio.h", which is in frmts/gtiff/libtiff. However, when compiling in this directory frmts/gtiff/libtiff is not incl

[gdal-dev] #include tiffio.h in gtiff/libgeotiff/

2018-04-12 Thread Ari Jolma
The file frmts/gtiff/libgeotiff/xtiffio.h #include's "tiffio.h", which is in frmts/gtiff/libtiff. However, when compiling in this directory frmts/gtiff/libtiff is not included. In my Linux Mint the compiler picks up tiffio.h from /usr/include/x86_64-linux-gnu/ but when I compile in freebsd I

[gdal-dev] Adding ignored fields to layers

2018-04-03 Thread Ari Jolma
When fields are added to a layer, the field definition that is stored in the layer definition is created as a clone of the provided field definition. Neither field definitions nor geometry field definitions set the ignored flag of the clone. Is this a feature and not a bug? I was already

Re: [gdal-dev] Wondering gdal_grid moving average interpolation

2018-03-13 Thread Ari Jolma
On 13.03.2018 09:54, Rutger wrote: It looks alright to me. The moving average algorithm searches (for each pixel) for points within the specified radius, and then averages the values of all those points. "gdal_grid" also supports inverse distance or linear interpolation if you want a more

[gdal-dev] Wondering gdal_grid moving average interpolation

2018-03-10 Thread Ari Jolma
I'm translating things on QGIS and now looking at GDAL tool rasterizing point data set with moving average interpolation. That's making a call to gdal_grid with -a average. I've got the attached result from a point set (circles with red border and inside inverse grayscale with interpolated

Re: [gdal-dev] layer union returns a geometry collection

2018-02-12 Thread Ari Jolma
Matthew, The documentation(1) says "The result layer contains features whose geometries represent areas that are either in the input layer, in the method layer, or in both.". I assume you a polygon in the input layer B which has a common singular point and common area with a polygon in layer

Re: [gdal-dev] Map Algebra revisited

2018-01-23 Thread Ari Jolma
an older discussion. At several points in the last year there has been a discussion of raster algebra / map algebra functionality on this list. There still is an associated RFC by Ari Jolma (https://trac.osgeo.org/gdal/wiki/rfc62_raster_algebra). The reason to re-open the discussion is that I have

Re: [gdal-dev] FFI Interface (to Perl)

2017-12-20 Thread Ari Jolma
o:schw...@gmail.com>> wrote: Cool!  Should be useful for other languages to see what is possible with libffi On Dec 20, 2017 3:23 PM, "Ari Jolma" <ari.jo...@gmail.com <mailto:ari.jo...@gmail.com>> wrote: Folks, I've been experimenting

[gdal-dev] FFI Interface (to Perl)

2017-12-20 Thread Ari Jolma
Folks, I've been experimenting with the Perl module FFI::Platypus, to create a Perl interface to GDAL with its help(1) and it's progressing well. FFI::Platypus itself uses libffi, which is a portable interface library. That, and Alien::gdal, which is a Perl module to use GDAL in a system

Re: [gdal-dev] port/md5 cvs_MD5*

2017-12-20 Thread Ari Jolma
uch but issues such as unique filenames and simultaneous access from threads/processes add complication. Ari Best regards, Dmitry 21.12.2017 0:32, Ari Jolma пишет: My comments on this: 1) The MD5 function seems to be needed only(?) for the WMS cache and probably also WFS. To

Re: [gdal-dev] port/md5 cvs_MD5*

2017-12-20 Thread Ari Jolma
My comments on this: 1) The MD5 function seems to be needed only(?) for the WMS cache and probably also WFS. To me that highlights the need for some integration work regarding caching in GDAL. I also implemented yet another caching code for the WCS (without MD5). Looking at my $HOME/.gdal I

[gdal-dev] gdal.pc

2017-12-14 Thread Ari Jolma
When gdal is built in the Strawberry Perl environment (which includes MinGW C/C++ compiler gcc) the subsequent use of libgdal in building other tools (the Perl bindings) is based on using gdal.pc and pkgconfig. Also, in that environment the linking of gdal with the dependent tools requires

Re: [gdal-dev] Shapefiles and z dimensions

2017-12-09 Thread Ari Jolma
Nyall Dawson kirjoitti 09.12.2017 klo 18:04: Hi all, I'm trying to track down some odd behaviour in QGIS, and determine whether this is a qgis bug, gdal bug, or shapefile limitation. If I create a new shapefile layer using GDALDatasetCreateLayer, with a geometry type of wkbPoint25D (actually

Re: [gdal-dev] pass gdal dataset as an argument from Python to C

2017-11-28 Thread Ari Jolma
Sean Gillies kirjoitti 27.11.2017 klo 20:46: I believe it's possible, yes, but it will be challenging. I'm not aware of anybody else doing this. I'm doing it but in the Perl side. It's not impossible but Swig does wrap the C pointers quite deep into the Perl object (it must be similar in

[gdal-dev] Build with MinGW in Linux

2017-11-27 Thread Ari Jolma
My PR https://github.com/OSGeo/gdal/pull/265 currently fails on MinGW cross compilation and on Windows (AppVeyor) while running a test involving gdal.Translate (python code). The error is "ValueError: Received a NULL pointer." which is from the Swig bindings. How could I study the issue on my

Re: [gdal-dev] Standardization of Bounding Box structure

2017-11-24 Thread Ari Jolma
Roberto Ribeiro kirjoitti 23.11.2017 klo 20:43: In OGR, the GetGeometry() method from the Geometry class returns a bounding box in the (xmin, xmax, ymin, ymax) structure. However, the SetSpatialFilterRect() method from the Layer class expects a bounding box in the (xmin, ymin, xmax, ymax)

Re: [gdal-dev] GDAL WCS driver: testing range and domain subsetting

2017-11-23 Thread Ari Jolma
Ari Jolma kirjoitti 23.11.2017 klo 15:40: I guess this goes to the Rasdaman people (My application to the Rasdaman email list is pending). Why does this query http://ows.rasdaman.org/rasdaman/ows?SERVICE=WCS=GetCoverage=2.0.1=test_irr_cube_2=image%2Ftiff=E%2875042.7273594,95042.7273594%29=N

[gdal-dev] GDAL WCS driver: testing range and domain subsetting

2017-11-23 Thread Ari Jolma
I guess this goes to the Rasdaman people (My application to the Rasdaman email list is pending). Why does this query http://ows.rasdaman.org/rasdaman/ows?SERVICE=WCS=GetCoverage=2.0.1=test_irr_cube_2=image%2Ftiff=E%2875042.7273594,95042.7273594%29=N%285434865.557937,5454865.557937%29=b2

Re: [gdal-dev] Multidimensional raster support in GDAL

2017-11-22 Thread Ari Jolma
Even Rouault kirjoitti 22.11.2017 klo 15:11: On mercredi 22 novembre 2017 14:38:55 CET Ari Jolma wrote: > Even Rouault kirjoitti 21.11.2017 klo 16:06: > > gdal_translate (actually GDALDatasetCopyWholeRaster()) uses the > > INTERLEAVE=PIXEL/BAND metadata item of the IMAGE_S

Re: [gdal-dev] Multidimensional raster support in GDAL

2017-11-22 Thread Ari Jolma
Even Rouault kirjoitti 21.11.2017 klo 16:06: gdal_translate (actually GDALDatasetCopyWholeRaster()) uses the INTERLEAVE=PIXEL/BAND metadata item of the IMAGE_STRUCTURE domain as a hint whether to read band by band or all bands together So if you add the following, it should request all

Re: [gdal-dev] Multidimensional raster support in GDAL

2017-11-21 Thread Ari Jolma
Even Rouault kirjoitti 21.11.2017 klo 18:51: > Somehow the user must be able to write > > dataset_1 = GDAL.Open("WCS:URL?coverage=coverage_a") > dataset_2 = GDAL.Open("WCS:URL?coverage=coverage_a") > > So that for example dataset_1 is band c from coverage_a and dataset_2 is > band d from

Re: [gdal-dev] Multidimensional raster support in GDAL

2017-11-21 Thread Ari Jolma
Even Rouault kirjoitti 21.11.2017 klo 16:06: > and different time steps as bands. My reading of the code and doc of the current driver is that time steps are treated as subdatasets, and not bands. Yes. My mistake. > I'm not sure if there are many people > using it since I believe the

[gdal-dev] Multidimensional raster support in GDAL

2017-11-21 Thread Ari Jolma
I'd like to get back to this topic a bit since I learned something from the standards that I believe was not discussed in the earlier discussion https://lists.osgeo.org/pipermail/gdal-dev/2017-October/047464.html In WCS 1 the range of the data - remember the main conceptual division to the

Re: [gdal-dev] The new WCS driver

2017-11-17 Thread Ari Jolma
Even Rouault kirjoitti 16.11.2017 klo 22:07: We have already various "linting" tools that run: cppcheck, clang static analyzer, coverity scan and a few ad-hoc scripts in scripts/ . But indeed we don't enforce things like maximum size of a function, etc. clang looks interesting. I

[gdal-dev] The new WCS driver

2017-11-16 Thread Ari Jolma
I'm happy to suggest that https://github.com/ajolma/gdal/tree/trunk/gdal/frmts/wcs would be accepted into the GDAL trunk, thus to be included into the coming version 2.3. I'm not making a PL since the GDAL at github is not a primary source. There is a python test program in the above

Re: [gdal-dev] Logics of CPLGetXMLValue and CPLFetchBool

2017-11-15 Thread Ari Jolma
Ari Jolma kirjoitti 15.11.2017 klo 15:51: Even Rouault kirjoitti 15.11.2017 klo 13:40: On mercredi 15 novembre 2017 10:58:23 CET Ari Jolma wrote: > > Thus the return value of CPLGetXMLBoolean can't use CPLGetXMLValue if it > follows the logic of CPLFetchBool. That is true is

Re: [gdal-dev] Logics of CPLGetXMLValue and CPLFetchBool

2017-11-15 Thread Ari Jolma
Even Rouault kirjoitti 15.11.2017 klo 13:40: On mercredi 15 novembre 2017 10:58:23 CET Ari Jolma wrote: > > Thus the return value of CPLGetXMLBoolean can't use CPLGetXMLValue if it > follows the logic of CPLFetchBool. That is true is existence and no denial. > Why not CPLG

[gdal-dev] Logics of CPLGetXMLValue and CPLFetchBool

2017-11-15 Thread Ari Jolma
I thought it would be perhaps good to have CPLGetXMLBoolean since we have CPLFetchBool and many of the option values needed by the new WCS driver are boolean flags and the option values end up in the service XML. It would twist the brain a bit less if one could simply use CPLGetXMLBoolean.

Re: [gdal-dev] WCS driver work

2017-11-14 Thread Ari Jolma
Ari Jolma kirjoitti 14.11.2017 klo 17:41: I've got very close to being able to get valid data (visual comparison in QGIS) from various server types and WCS versions. I've documented the results here: https://trac.osgeo.org/gdal/wiki/WCS%2Binteroperability It seems that I need eight option

Re: [gdal-dev] Wiki page for WCS discussion and best practices created

2017-11-14 Thread Ari Jolma
jratike80 kirjoitti 13.11.2017 klo 15:39: Ari Jolma-2 wrote I made a bold attempt to begin describing the fundamentals of coordinate systems in WCS and how those relate to the fundamentals of coordinate systems in GDAL. Ari And I have - Materialized my 6 pixel test image into a GeoTIFF

[gdal-dev] WCS driver work

2017-11-14 Thread Ari Jolma
I've got very close to being able to get valid data (visual comparison in QGIS) from various server types and WCS versions. I've documented the results here: https://trac.osgeo.org/gdal/wiki/WCS%2Binteroperability It seems that I need eight option flags to cover all cases. Mostly the

Re: [gdal-dev] Wiki page for WCS discussion and best practices created

2017-11-13 Thread Ari Jolma
jratike80 kirjoitti 13.11.2017 klo 15:39: Ari Jolma-2 wrote I made a bold attempt to begin describing the fundamentals of coordinate systems in WCS and how those relate to the fundamentals of coordinate systems in GDAL. Ari And I have - Materialized my 6 pixel test image into a GeoTIFF

Re: [gdal-dev] Numeric accuracy

2017-11-11 Thread Ari Jolma
Ari Jolma kirjoitti 11.11.2017 klo 12:05: I have been bitten twice, once with ArcGIS and now with Rasdaman WCS, with numeric accuracy. I'm making a data request on the corner of the bounding box, let's say it has minimum X of 75042.7273594. I'm setting my minX to that value and I'm enforcing

[gdal-dev] Numeric accuracy

2017-11-11 Thread Ari Jolma
I have been bitten twice, once with ArcGIS and now with Rasdaman WCS, with numeric accuracy. I'm making a data request on the corner of the bounding box, let's say it has minimum X of 75042.7273594. I'm setting my minX to that value and I'm enforcing it to that value with MAX (this was

Re: [gdal-dev] Wiki page for WCS discussion and best practices created

2017-11-11 Thread Ari Jolma
Rahkonen Jukka (MML) kirjoitti 10.11.2017 klo 10:44: Hi, As an initial step I created a new page about WCS into the developers section of the GDAL wiki https://trac.osgeo.org/gdal/wiki/WCS%2Binteroperability I made a bold attempt to begin describing the fundamentals of coordinate systems in

Re: [gdal-dev] WCS GetCoverage with AxisOrder swap

2017-11-09 Thread Ari Jolma
Ari Jolma kirjoitti 09.11.2017 klo 15:52: jratike80 kirjoitti 09.11.2017 klo 15:18: Piero Campa wrote And for WCS subsetting, when you write something like SUBSET=AXIS_NAME(min,max), ​ ​ where AXIS_NAME should come from ? From the RectifiedGrid.axisLabels I guess ? ​ ​NO, subsetting

Re: [gdal-dev] WCS GetCoverage with AxisOrder swap

2017-11-09 Thread Ari Jolma
jratike80 kirjoitti 09.11.2017 klo 15:18: Piero Campa wrote And for WCS subsetting, when you write something like SUBSET=AXIS_NAME(min,max), ​ ​ where AXIS_NAME should come from ? From the RectifiedGrid.axisLabels I guess ? ​ ​NO, subsetting/subspacing​ is meant to work aligned with the

Re: [gdal-dev] WCS GetCoverage with AxisOrder swap

2017-11-08 Thread Ari Jolma
Even Rouault kirjoitti 08.11.2017 klo 17:31: On mercredi 8 novembre 2017 15:50:48 CET Piero Campalani wrote: > This order is what can be referred to in the GridFunction: +1 --> first > GRID axis, +2 --> second GRID axis, etc. Ok, this tells me that the axisOrder attribute in sequenceRule

Re: [gdal-dev] WCS GetCoverage with AxisOrder swap

2017-11-08 Thread Ari Jolma
Ari Jolma kirjoitti 08.11.2017 klo 13:48:  for MapServer and GeoServer I need to have at least two hack options (NoOffsetSwap and NoGridEnvelopeSwap, maybe there could be only one). I have not yet checked if the rasters returned from those two are ok. The rasters that those servers return

[gdal-dev] WCS GetCoverage with AxisOrder swap

2017-11-08 Thread Ari Jolma
I've got to the point in the new WCS driver development, where I'm trying to issue GetCoverage request to a GeoServer, which is serving data in EPSG that has axis order swapped. This is related to the earlier discussion https://lists.osgeo.org/pipermail/gdal-dev/2017-April/046366.html The

Re: [gdal-dev] WCS 2.0, getting (multidimensional) data

2017-11-07 Thread Ari Jolma
Peter Baumann kirjoitti 07.11.2017 klo 16:29: Ari- On 11/07/2017 09:56 AM, Ari Jolma wrote: To get multidimensional data into GDAL 2D domain, one may need to slice the data at some point in the non x/y dimensions. I'm testing this with the Rasdaman server. My code currently creates a KVP

Re: [gdal-dev] WCS 2.0, getting (multidimensional) data

2017-11-07 Thread Ari Jolma
ntralMed-MCol=image%2Ftiff=long(9.83125,9.83958334)=lat(46.18124999,46.18958333)=1 works. Ari Ari Jolma kirjoitti 03.11.2017 klo 09:57: I have two things that I need to ask/discuss regarding WCS 2.0 1) My understanding is that there are potentially three CRS in operation: the CR

Re: [gdal-dev] WCS 2.0, getting (multidimensional) data

2017-11-03 Thread Ari Jolma
Even Rouault kirjoitti 03.11.2017 klo 11:07: > 1) My understanding is that there are potentially three CRS in > operation: the CRS of the bbox of the coverage (the native CRS of the > server); the CRS of the grid; and the CRS where the user wants the data. > I can imagine that it is

Re: [gdal-dev] WCS 2.0, getting (multidimensional) data

2017-11-03 Thread Ari Jolma
Peter Baumann kirjoitti 03.11.2017 klo 10:54: On 11/03/2017 08:57 AM, Ari Jolma wrote: I have two things that I need to ask/discuss regarding WCS 2.0 1) My understanding is that there are potentially three CRS in operation: the CRS of the bbox of the coverage (the native CRS of the server

Re: [gdal-dev] WCS 2.0, getting (multidimensional) data

2017-11-03 Thread Ari Jolma
one detail forgotten Ari Jolma kirjoitti 03.11.2017 klo 09:57: I have two things that I need to ask/discuss regarding WCS 2.0 1) My understanding is that there are potentially three CRS in operation: the CRS of the bbox of the coverage (the native CRS of the server); the CRS of the grid

[gdal-dev] WCS 2.0, getting (multidimensional) data

2017-11-03 Thread Ari Jolma
I have two things that I need to ask/discuss regarding WCS 2.0 1) My understanding is that there are potentially three CRS in operation: the CRS of the bbox of the coverage (the native CRS of the server); the CRS of the grid; and the CRS where the user wants the data. I can imagine that it is

Re: [gdal-dev] Reading sentinel 2 L2A datasets

2017-11-03 Thread Ari Jolma
Just a comment on this post. It was initially blocked since it is > 50 kB. I've now increased the size limit to 100 kB. Ari Róbert Špir kirjoitti 03.11.2017 klo 08:55: Hi all, I'm using GDAL to process sentinel-2 data (in c++). I have no problem processing S2B and S2A L1C datasets, I can

Re: [gdal-dev] WCS 1.1 in GDAL

2017-11-01 Thread Ari Jolma
Even Rouault kirjoitti 01.11.2017 klo 13:20: On mercredi 1 novembre 2017 13:05:09 CET Ari Jolma wrote: > Anyway, the question is, should gdalgmlcoverage.cpp be in the gcore? Feel free to move it to frmts/wcs as it only used by it. You'd probably want to remove the declaration of it f

Re: [gdal-dev] Proposed changed to solve inconsistant spatial filtering for null or empty geometries

2017-11-01 Thread Ari Jolma
Even Rouault kirjoitti 01.11.2017 klo 13:02: Hi, Nyall Dawson pointed to me in https://github.com/qgis/QGIS/pull/5501#issuecomment-340934947 about a weird behaviour of OGR spatial filtering regarding null and empty geometries. ... Any thought ? Yes. Deep ones. If something is

Re: [gdal-dev] WCS 1.1 in GDAL

2017-11-01 Thread Ari Jolma
Ari Jolma kirjoitti 01.11.2017 klo 11:04: It seems to me that WCS 1.1 has been dormant and perhaps nobody has really used it ever(?). I just found the gdalgmlcoverage.cpp file in gcore, which is only used by the WCS driver and only in its 1.1 part Actually, in 1.0, since 1.1 does not seem

[gdal-dev] WCS 1.1 in GDAL

2017-11-01 Thread Ari Jolma
It seems to me that WCS 1.1 has been dormant and perhaps nobody has really used it ever(?). I just found the gdalgmlcoverage.cpp file in gcore, which is only used by the WCS driver and only in its 1.1 part and it hasn't really been changed since Frank wrote it ten years ago. I basically

[gdal-dev] Dangers or letting a domain name expire

2017-10-31 Thread Ari Jolma
I had the domain name ajolma.net at some point and I kept there the documentation of the Perl bindings. Now they are at arijolma.org. At some point I let the name expire but I have not been very good at changing it everywhere. It is still for example in

[gdal-dev] SRS definition

2017-10-31 Thread Ari Jolma
Rasdaman writes SRS definitions such as "http://ows.rasdaman.org/def/crs/EPSG/0/32633; in WCS coverage descriptions. The URL returns a gml:ProjectedCRS XML. GDAL fetches the XML but then gives up in SetFromUserInput since it only expects "..that this method will be extended in the future to

Re: [gdal-dev] WCS driver

2017-10-31 Thread Ari Jolma
Even Rouault kirjoitti 30.10.2017 klo 22:41: On lundi 30 octobre 2017 10:54:56 CET Ari Jolma wrote: > > XML namespaces can be declared in any element and then they are in > effect in that element and in all of its children unless they override > the namespace declarations.

Re: [gdal-dev] OGC announces a new standard that improves the way information is referenced to the earth

2017-10-31 Thread Ari Jolma
That also caught my eye. The text sounds a bit like marketing talk but maybe there is something. From a quick look my understanding is that the idea is to create a grid that divides the whole earth into cells of similar shape in a sequence of increasing cell size. And that sounds to me like a

  1   2   3   4   5   6   >