[gdal-dev] [EXTERNAL] Re: SQL query to set column to GPKG raster pixel value of point intersection result

[Meyer, Jesse R. (GSFC-618.0)[SCIENCE SYSTEMS AND APPLICATIONS INC]]

Hi,

I want to perform this on tens of thousands of GPKG databases, each possessing upwards of millions of points, so automation (and performance) is a high value.  The hope was that given both data sources leverage the same driver, there was some hope of this working with the proper incantation, especially considering the library routines available by Spatialite in SQL queries.  But without something like burn available in a SQL query then this can’t work at the query level.  This is easy to write in C/C++ but wanted to see if the database engine could handle it directly.

Jesse

From: Rahkonen Jukka <jukka.rahkonen at maanmittauslaitos.fi>
Date: Thursday, December 8, 2022 at 3:45 PM
To: "Meyer, Jesse R. (GSFC-618.0)[SCIENCE SYSTEMS AND APPLICATIONS INC]" <jesse.r.meyer at nasa.gov>, "gdal-dev at lists.osgeo.org" <gdal-dev at lists.osgeo.org>
Subject: [EXTERNAL] Re: SQL query to set column to GPKG raster pixel value of point intersection result

Hi,

It feels like quite special tool. Have you considered to use the Sample Raster Values tool in QGIS? The description of the tool is:

Sample raster values
This algorithm creates a new vector layer with the same attributes of the input layer and the raster values corresponding on the point location.
If the raster layer has more than one band, all the band values are sampled.

-Jukka Rahkonen-

Lähettäjä: gdal-dev <gdal-dev-bounces at lists.osgeo.org> Puolesta Meyer, Jesse R. (GSFC-618.0)[SCIENCE SYSTEMS AND APPLICATIONS INC] via gdal-dev
Lähetetty: torstai 8. joulukuuta 2022 20.07
Vastaanottaja: gdal-dev at lists.osgeo.org
Aihe: [gdal-dev] SQL query to set column to GPKG raster pixel value of point intersection result

Hi,

I have GPKG database with features represented by point geometries.  I’d like to update these features, in a single SQL statement, to add a new column whose value is the pixel value of an external GPKG raster for which the point geometry intersects with.  I was hoping to find an OGR SQL function similar to burn for this, but I’ve come up empty handed.

Is this functionality available? Ideally, without resorting to manually updating each feature individually in a user application.  Otherwise, perhaps a feature request is in order.

Thanks,
Jesse
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20221208/a1e3e74e/attachment-0001.htm>