gdal vector pipeline read --layer: make sure ExecuteSQL() forwards to the source dataset

Author: rouault

apps/gdalalg_vector_read.cpp

@@ -61,6 +61,10 @@ class GDALVectorPipelineReadOutputDataset final : public GDALDataset
 
     void ResetReading() override;
 
+    OGRLayer *ExecuteSQL(const char *pszStatement, OGRGeometry *poSpatialFilter,
+                         const char *pszDialect) override;
+    void ReleaseResultSet(OGRLayer *poResultsSet) override;
+
     OGRFeature *GetNextFeature(OGRLayer **ppoBelongingLayer,
                                double *pdfProgressPct,
                                GDALProgressFunc pfnProgress,
@@ -137,6 +141,28 @@ void GDALVectorPipelineReadOutputDataset::ResetReading()
     m_srcDS.ResetReading();
 }
 
+/************************************************************************/
+/*          GDALVectorPipelineReadOutputDataset::ExecuteSQL()           */
+/************************************************************************/
+
+OGRLayer *
+GDALVectorPipelineReadOutputDataset::ExecuteSQL(const char *pszStatement,
+                                                OGRGeometry *poSpatialFilter,
+                                                const char *pszDialect)
+{
+    return m_srcDS.ExecuteSQL(pszStatement, poSpatialFilter, pszDialect);
+}
+
+/************************************************************************/
+/*       GDALVectorPipelineReadOutputDataset::ReleaseResultSet()        */
+/************************************************************************/
+
+void GDALVectorPipelineReadOutputDataset::ReleaseResultSet(
+    OGRLayer *poResultsSet)
+{
+    m_srcDS.ReleaseResultSet(poResultsSet);
+}
+
 /************************************************************************/
 /*        GDALVectorPipelineReadOutputDataset::GetNextFeature()         */
 /************************************************************************/

autotest/utilities/test_gdalalg_vector_pipeline.py

@@ -1091,3 +1091,18 @@ def test_gdalalg_vector_pipeline_read_ds_take_ref(tmp_vsimem):
     ds = ogr.Open(f"{tmp_vsimem}/out.shp")
     lyr = ds.GetLayer(0)
     assert lyr.GetFeatureCount() == 10
+
+
+@pytest.mark.require_driver("GPKG")
+def test_gdalalg_vector_pipeline_read_execute_sql(tmp_vsimem):
+
+    tmp_filename = tmp_vsimem / "tmp.gpkg"
+    gdal.alg.vector.convert(input="../ogr/data/poly.shp", output=tmp_filename)
+
+    gdal.alg.vector.pipeline(
+        pipeline=f'read {tmp_filename} --layer poly ! sql --sql "select * from poly group by eas_id" ! write {tmp_vsimem}/out.shp'
+    )
+
+    ds = ogr.Open(f"{tmp_vsimem}/out.shp")
+    lyr = ds.GetLayer(0)
+    assert lyr.GetFeatureCount() == 10