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

[Even Rouault]

Hi,

That sounds like a cool idea. I've had a crack at it in 
https://github.com/OSGeo/gdal/pull/6877

Even

Le 08/12/2022 à 22:10, Meyer, Jesse R. (GSFC-618.0)[SCIENCE SYSTEMS AND 
APPLICATIONS INC] via gdal-dev a écrit :
>
> 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
>
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev

-- 
http://www.spatialys.com
My software is free, but my time generally not.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20221208/b21f0755/attachment.htm>