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.rahko...@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.me...@nasa.gov>, "gdal-dev@lists.osgeo.org" <gdal-dev@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-boun...@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@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@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev
