gdalalg_vector_pipeline.h

/******************************************************************************
 *
 * Project:  GDAL
 * Purpose:  gdal "vector pipeline" subcommand
 * Author:   Even Rouault <even dot rouault at spatialys.com>
 *
 ******************************************************************************
 * Copyright (c) 2024, Even Rouault <even dot rouault at spatialys.com>
 *
 * SPDX-License-Identifier: MIT
 ****************************************************************************/
 
#ifndef GDALALG_VECTOR_PIPELINE_INCLUDED
#define GDALALG_VECTOR_PIPELINE_INCLUDED
 
#include "gdalalgorithm.h"
#include "gdalalg_abstract_pipeline.h"
 
#include "ogrsf_frmts.h"
#include "ogrlayerwithtranslatefeature.h"
 
#include <map>
#include <tuple>
#include <vector>
 
 
/************************************************************************/
/*                   GDALVectorPipelineStepAlgorithm                    */
/************************************************************************/
 
class GDALRasterAlgorithmStepRegistry;
 
class GDALVectorPipelineStepAlgorithm /* non final */
    : public GDALPipelineStepAlgorithm
{
  public:
    ~GDALVectorPipelineStepAlgorithm() override;
 
  protected:
    GDALVectorPipelineStepAlgorithm(const std::string &name,
                                    const std::string &description,
                                    const std::string &helpURL,
                                    bool standaloneStep);
 
    GDALVectorPipelineStepAlgorithm(const std::string &name,
                                    const std::string &description,
                                    const std::string &helpURL,
                                    const ConstructorOptions &options);
 
    friend class GDALVectorPipelineAlgorithm;
    friend class GDALVectorConcatAlgorithm;
 
    int GetInputType() const override
    {
        return GDAL_OF_VECTOR;
    }
 
    int GetOutputType() const override
    {
        return GDAL_OF_VECTOR;
    }
};
 
/************************************************************************/
/*                   GDALVectorAlgorithmStepRegistry                    */
/************************************************************************/
 
class GDALVectorAlgorithmStepRegistry : public virtual GDALAlgorithmRegistry
{
  public:
    GDALVectorAlgorithmStepRegistry() = default;
    ~GDALVectorAlgorithmStepRegistry() override;
 
    template <class MyAlgorithm>
    bool Register(const std::string &name = std::string())
    {
        static_assert(
            std::is_base_of_v<GDALVectorPipelineStepAlgorithm, MyAlgorithm>,
            "Algorithm is not a GDALVectorPipelineStepAlgorithm");
 
        AlgInfo info;
        info.m_name = name.empty() ? MyAlgorithm::NAME : name;
        info.m_aliases = MyAlgorithm::GetAliasesStatic();
        info.m_creationFunc = []() -> std::unique_ptr<GDALAlgorithm>
        { return std::make_unique<MyAlgorithm>(); };
        return GDALAlgorithmRegistry::Register(info);
    }
};
 
/************************************************************************/
/*                     GDALVectorPipelineAlgorithm                      */
/************************************************************************/
 
class GDALVectorPipelineAlgorithm final : public GDALAbstractPipelineAlgorithm
{
  public:
    static constexpr const char *NAME = "pipeline";
    static constexpr const char *DESCRIPTION =
        "Process a vector dataset applying several steps.";
    static constexpr const char *HELP_URL =
        "/programs/gdal_vector_pipeline.html";
 
    static std::vector<std::string> GetAliasesStatic()
    {
        return {
#ifdef GDAL_PIPELINE_PROJ_NOSTALGIA
            GDALAlgorithmRegistry::HIDDEN_ALIAS_SEPARATOR,
            "+pipeline",
            "+gdal=pipeline",
#endif
        };
    }
 
    GDALVectorPipelineAlgorithm();
 
    std::string GetUsageForCLI(bool shortUsage,
                               const UsageOptions &usageOptions) const override;
 
    static void RegisterAlgorithms(GDALVectorAlgorithmStepRegistry &registry,
                                   bool forMixedPipeline);
 
    int GetInputType() const override
    {
        return GDAL_OF_VECTOR;
    }
 
    int GetOutputType() const override
    {
        return GDAL_OF_VECTOR;
    }
 
  protected:
    GDALVectorAlgorithmStepRegistry m_stepRegistry{};
 
    GDALAlgorithmRegistry &GetStepRegistry() override
    {
        return m_stepRegistry;
    }
 
    const GDALAlgorithmRegistry &GetStepRegistry() const override
    {
        return m_stepRegistry;
    }
 
  private:
    std::unique_ptr<GDALAbstractPipelineAlgorithm>
    CreateNestedPipeline() const override
    {
        auto pipeline = std::make_unique<GDALVectorPipelineAlgorithm>();
        pipeline->m_bInnerPipeline = true;
        return pipeline;
    }
};
 
/************************************************************************/
/*                      GDALVectorDecoratedDataset                      */
/************************************************************************/
 
class GDALVectorDecoratedDataset /* non final */
    : public GDALDataset
{
  public:
    ~GDALVectorDecoratedDataset() override;
 
    CSLConstList GetMetadata(const char *pszDomain) override
    {
        return m_srcDS.GetMetadata(pszDomain);
    }
 
    const char *GetMetadataItem(const char *pszName,
                                const char *pszDomain) override
    {
        return m_srcDS.GetMetadataItem(pszName, pszDomain);
    }
 
    std::vector<std::string>
    GetFieldDomainNames(CSLConstList papszOptions) const override
    {
        return m_srcDS.GetFieldDomainNames(papszOptions);
    }
 
    const OGRFieldDomain *GetFieldDomain(const std::string &name) const override
    {
        return m_srcDS.GetFieldDomain(name);
    }
 
    std::vector<std::string>
    GetRelationshipNames(CSLConstList papszOptions) const override
    {
        return m_srcDS.GetRelationshipNames(papszOptions);
    }
 
    const GDALRelationship *
    GetRelationship(const std::string &name) const override
    {
        return m_srcDS.GetRelationship(name);
    }
 
  private:
    std::unique_ptr<GDALDataset> m_dummySrcDS{};
 
  protected:
    explicit GDALVectorDecoratedDataset(GDALDataset *poSrcDS);
 
    GDALDataset &m_srcDS;
};
 
/************************************************************************/
/*                       GDALVectorOutputDataset                        */
/************************************************************************/
 
class GDALVectorOutputDataset : public GDALVectorDecoratedDataset
{
 
  public:
    explicit GDALVectorOutputDataset(GDALDataset *poSrcDS);
 
    int GetLayerCount() const override
    {
        return static_cast<int>(m_layers.size());
    }
 
    const OGRLayer *GetLayer(int idx) const override
    {
        return idx >= 0 && idx < GetLayerCount() ? m_layers[idx].get()
                                                 : nullptr;
    }
 
    int TestCapability(const char *) const override;
 
    void AddLayer(std::unique_ptr<OGRLayer> layer)
    {
        m_layers.emplace_back(std::move(layer));
    }
 
  private:
    std::vector<std::unique_ptr<OGRLayer>> m_layers{};
};
 
/************************************************************************/
/*                GDALVectorAlgorithmLayerProgressHelper                */
/************************************************************************/
 
class GDALVectorAlgorithmLayerProgressHelper
{
  public:
    GDALVectorAlgorithmLayerProgressHelper(GDALProgressFunc pfnProgress,
                                           void *pProgressData);
    explicit GDALVectorAlgorithmLayerProgressHelper(
        const GDALPipelineStepRunContext &ctxt);
 
    void AddProcessedLayer(OGRLayer &srcLayer);
 
    void AddPassThroughLayer(OGRLayer &srcLayer);
 
    class iterator
    {
      public:
        explicit iterator(const GDALVectorAlgorithmLayerProgressHelper &helper,
                          bool start)
            : m_helper(helper),
              m_nLayerIdx(start ? 0 : m_helper.m_apoSrcLayers.size())
        {
        }
 
        inline bool operator==(const iterator &other) const
        {
            return m_nLayerIdx == other.m_nLayerIdx;
        }
 
        inline bool operator!=(const iterator &other) const
        {
            return m_nLayerIdx != other.m_nLayerIdx;
        }
 
        inline iterator &operator++()
        {
            if (!m_helper.m_anFeatures.empty())
                m_nFeatureIdx += m_helper.m_anFeatures[m_nProcessedLayerIdx];
            if (m_helper.m_apoSrcLayers[m_nLayerIdx].second)
                ++m_nProcessedLayerIdx;
            ++m_nLayerIdx;
            return *this;
        }
 
        using progress_data_unique_ptr =
            std::unique_ptr<void, decltype(&GDALDestroyScaledProgress)>;
        using value_type = std::tuple<OGRLayer *, bool, GDALProgressFunc,
                                      progress_data_unique_ptr>;
 
        value_type operator*() const;
 
      private:
        const GDALVectorAlgorithmLayerProgressHelper &m_helper;
        size_t m_nLayerIdx = 0;
        size_t m_nProcessedLayerIdx = 0;
        GIntBig m_nFeatureIdx = 0;
    };
 
 
    iterator begin() const
    {
        return iterator(*this, true);
    }
 
    iterator end() const
    {
        return iterator(*this, false);
    }
 
    bool HasProcessedLayers() const
    {
        return !m_anFeatures.empty();
    }
 
  private:
    GDALProgressFunc m_pfnProgress = nullptr;
    void *m_pProgressData = nullptr;
    int64_t m_nTotalFeatures = 0;
    std::vector<std::pair<OGRLayer *, bool>> m_apoSrcLayers{};
    std::vector<int64_t> m_anFeatures{};
 
    CPL_DISALLOW_COPY_ASSIGN(GDALVectorAlgorithmLayerProgressHelper)
};
 
/************************************************************************/
/*                    GDALVectorPipelineOutputLayer                     */
/************************************************************************/
 
class GDALVectorPipelineOutputLayer /* non final */
    : public OGRLayerWithTranslateFeature
{
  protected:
    explicit GDALVectorPipelineOutputLayer(OGRLayer &oSrcLayer);
    ~GDALVectorPipelineOutputLayer() override;
 
    OGRLayer &m_srcLayer;
 
    void FailTranslation()
    {
        m_translateError = true;
    }
 
    bool PassesFilters(OGRFeature *poFeature);
 
  public:
    void ResetReading() override;
 
    OGRFeature *GetNextFeature() override final;
 
  private:
    std::vector<std::unique_ptr<OGRFeature>> m_pendingFeatures{};
    size_t m_idxInPendingFeatures = 0;
    bool m_translateError = false;
};
 
/************************************************************************/
/*                  GDALVectorPipelinePassthroughLayer                  */
/************************************************************************/
 
class GDALVectorPipelinePassthroughLayer /* non final */
    : public GDALVectorPipelineOutputLayer
{
  public:
    explicit GDALVectorPipelinePassthroughLayer(OGRLayer &oSrcLayer)
        : GDALVectorPipelineOutputLayer(oSrcLayer)
    {
        SetDescription(oSrcLayer.GetDescription());
    }
 
    const OGRFeatureDefn *GetLayerDefn() const override;
 
    int TestCapability(const char *pszCap) const override
    {
        if (EQUAL(pszCap, OLCFastFeatureCount))
            return false;
        return m_srcLayer.TestCapability(pszCap);
    }
 
    OGRErr IGetExtent(int iGeomField, OGREnvelope *psExtent,
                      bool bForce) override
    {
        return m_srcLayer.GetExtent(iGeomField, psExtent, bForce);
    }
 
    OGRErr IGetExtent3D(int iGeomField, OGREnvelope3D *psExtent,
                        bool bForce) override
    {
        return m_srcLayer.GetExtent3D(iGeomField, psExtent, bForce);
    }
 
    void TranslateFeature(
        std::unique_ptr<OGRFeature> poSrcFeature,
        std::vector<std::unique_ptr<OGRFeature>> &apoOutFeatures) override;
};
 
/************************************************************************/
/*                 GDALVectorNonStreamingAlgorithmLayer                 */
/************************************************************************/
 
class GDALVectorNonStreamingAlgorithmLayer
    : public OGRLayer,
      public OGRGetNextFeatureThroughRaw<GDALVectorNonStreamingAlgorithmLayer>
{
  public:
    GDALVectorNonStreamingAlgorithmLayer(OGRLayer &srcLayer,
                                         int geomFieldIndex);
 
    ~GDALVectorNonStreamingAlgorithmLayer() override;
 
    const char *GetDescription() const override
    {
        return GetName();
    }
 
    virtual bool Process(GDALProgressFunc pfnProgress, void *pProgressData) = 0;
 
    virtual std::unique_ptr<OGRFeature> GetNextProcessedFeature() = 0;
 
    OGRFeature *GetNextRawFeature();
 
    DEFINE_GET_NEXT_FEATURE_THROUGH_RAW(GDALVectorNonStreamingAlgorithmLayer)
 
  protected:
    OGRLayer &m_srcLayer;
    int m_geomFieldIndex{0};
 
  private:
    CPL_DISALLOW_COPY_ASSIGN(GDALVectorNonStreamingAlgorithmLayer)
};
 
/************************************************************************/
/*                GDALVectorNonStreamingAlgorithmDataset                */
/************************************************************************/
 
class GDALVectorNonStreamingAlgorithmDataset /* non final */
    : public GDALVectorDecoratedDataset
{
  public:
    explicit GDALVectorNonStreamingAlgorithmDataset(GDALDataset &oSrcDS);
    ~GDALVectorNonStreamingAlgorithmDataset() override;
 
    bool AddProcessedLayer(
        std::unique_ptr<GDALVectorNonStreamingAlgorithmLayer> srcLayer,
        GDALProgressFunc progressFunc, void *progressData);
 
    void AddPassThroughLayer(OGRLayer &oLayer);
 
    int GetLayerCount() const final override;
    OGRLayer *GetLayer(int idx) const final override;
    int TestCapability(const char *pszCap) const override;
 
  private:
    std::vector<std::unique_ptr<OGRLayer>> m_owned_layers{};
    std::vector<OGRLayer *> m_layers{};
};
 
/************************************************************************/
/*                   GDALVectorPipelineOutputDataset                    */
/************************************************************************/
 
class GDALVectorPipelineOutputDataset /* non final */
    : public GDALVectorDecoratedDataset
{
    std::map<OGRLayer *, OGRLayerWithTranslateFeature *>
        m_mapSrcLayerToNewLayer{};
    std::vector<std::unique_ptr<OGRLayerWithTranslateFeature>>
        m_layersToDestroy{};
    std::vector<OGRLayerWithTranslateFeature *> m_layers{};
 
    OGRLayerWithTranslateFeature *m_belongingLayer = nullptr;
    std::vector<std::unique_ptr<OGRFeature>> m_pendingFeatures{};
    size_t m_idxInPendingFeatures = 0;
 
    CPL_DISALLOW_COPY_ASSIGN(GDALVectorPipelineOutputDataset)
 
  public:
    explicit GDALVectorPipelineOutputDataset(GDALDataset &oSrcDS);
 
    void AddLayer(OGRLayer &oSrcLayer,
                  std::unique_ptr<OGRLayerWithTranslateFeature> poNewLayer);
 
    int GetLayerCount() const override;
 
    OGRLayer *GetLayer(int idx) const override;
 
    int TestCapability(const char *pszCap) const override;
 
    void ResetReading() override;
 
    OGRFeature *GetNextFeature(OGRLayer **ppoBelongingLayer,
                               double *pdfProgressPct,
                               GDALProgressFunc pfnProgress,
                               void *pProgressData) override;
};
 
 
#endif