/******************************************************************************
 *
 * Project:  GDAL
 * Purpose:  GDALAlgorithm class
 * Author:   Even Rouault <even dot rouault at spatialys.com>
 *
 ******************************************************************************
 * Copyright (c) 2024, Even Rouault <even dot rouault at spatialys.com>
 *
 * SPDX-License-Identifier: MIT
 ****************************************************************************/

#include "cpl_port.h"
#include "cpl_conv.h"
#include "cpl_error.h"
#include "cpl_json.h"
#include "cpl_levenshtein.h"
#include "cpl_minixml.h"
#include "cpl_multiproc.h"

#include "gdalalgorithm.h"
#include "gdalalg_abstract_pipeline.h"
#include "gdal_priv.h"
#include "gdal_thread_pool.h"
#include "ogrsf_frmts.h"
#include "ogr_spatialref.h"
#include "vrtdataset.h"

#include <algorithm>
#include <cassert>
#include <cerrno>
#include <cmath>
#include <cstdlib>
#include <limits>
#include <map>
#include <type_traits>
#include <string_view>
#include <regex>

#ifndef _
#define _(x) (x)
#endif

constexpr const char *GDAL_ARG_NAME_OUTPUT_DATA_TYPE = "output-data-type";

constexpr const char *GDAL_ARG_NAME_OUTPUT_OPEN_OPTION = "output-open-option";

constexpr const char *GDAL_ARG_NAME_BAND = "band";

//! @cond Doxygen_Suppress
struct GDALAlgorithmArgHS
{
    GDALAlgorithmArg *ptr = nullptr;

    explicit GDALAlgorithmArgHS(GDALAlgorithmArg *arg) : ptr(arg)
    {
    }
};

//! @endcond

//! @cond Doxygen_Suppress
struct GDALArgDatasetValueHS
{
    GDALArgDatasetValue val{};
    GDALArgDatasetValue *ptr = nullptr;

    GDALArgDatasetValueHS() : ptr(&val)
    {
    }

    explicit GDALArgDatasetValueHS(GDALArgDatasetValue *arg) : ptr(arg)
    {
    }

    GDALArgDatasetValueHS(const GDALArgDatasetValueHS &) = delete;
    GDALArgDatasetValueHS &operator=(const GDALArgDatasetValueHS &) = delete;
};

//! @endcond

/************************************************************************/
/*                     GDALAlgorithmArgTypeIsList()                     */
/************************************************************************/

bool GDALAlgorithmArgTypeIsList(GDALAlgorithmArgType type)
{
    switch (type)
    {
        case GAAT_BOOLEAN:
        case GAAT_STRING:
        case GAAT_INTEGER:
        case GAAT_REAL:
        case GAAT_DATASET:
            break;

        case GAAT_STRING_LIST:
        case GAAT_INTEGER_LIST:
        case GAAT_REAL_LIST:
        case GAAT_DATASET_LIST:
            return true;
    }

    return false;
}

/************************************************************************/
/*                      GDALAlgorithmArgTypeName()                      */
/************************************************************************/

const char *GDALAlgorithmArgTypeName(GDALAlgorithmArgType type)
{
    switch (type)
    {
        case GAAT_BOOLEAN:
            break;
        case GAAT_STRING:
            return "string";
        case GAAT_INTEGER:
            return "integer";
        case GAAT_REAL:
            return "real";
        case GAAT_DATASET:
            return "dataset";
        case GAAT_STRING_LIST:
            return "string_list";
        case GAAT_INTEGER_LIST:
            return "integer_list";
        case GAAT_REAL_LIST:
            return "real_list";
        case GAAT_DATASET_LIST:
            return "dataset_list";
    }

    return "boolean";
}

/************************************************************************/
/*                  GDALAlgorithmArgDatasetTypeName()                   */
/************************************************************************/

std::string GDALAlgorithmArgDatasetTypeName(GDALArgDatasetType type)
{
    std::string ret;
    if ((type & GDAL_OF_RASTER) != 0)
        ret = "raster";
    if ((type & GDAL_OF_VECTOR) != 0)
    {
        if (!ret.empty())
        {
            if ((type & GDAL_OF_MULTIDIM_RASTER) != 0)
                ret += ", ";
            else
                ret += " or ";
        }
        ret += "vector";
    }
    if ((type & GDAL_OF_MULTIDIM_RASTER) != 0)
    {
        if (!ret.empty())
        {
            ret += " or ";
        }
        ret += "multidimensional raster";
    }
    return ret;
}

/************************************************************************/
/*                        GDALAlgorithmArgDecl()                        */
/************************************************************************/

// cppcheck-suppress uninitMemberVar
GDALAlgorithmArgDecl::GDALAlgorithmArgDecl(const std::string &longName,
                                           char chShortName,
                                           const std::string &description,
                                           GDALAlgorithmArgType type)
    : m_longName(longName),
      m_shortName(chShortName ? std::string(&chShortName, 1) : std::string()),
      m_description(description), m_type(type),
      m_metaVar(CPLString(m_type == GAAT_BOOLEAN ? std::string() : longName)
                    .toupper()),
      m_maxCount(GDALAlgorithmArgTypeIsList(type) ? UNBOUNDED : 1)
{
    if (m_type == GAAT_BOOLEAN)
    {
        m_defaultValue = false;
    }
}

/************************************************************************/
/*                 GDALAlgorithmArgDecl::SetMinCount()                  */
/************************************************************************/

GDALAlgorithmArgDecl &GDALAlgorithmArgDecl::SetMinCount(int count)
{
    if (!GDALAlgorithmArgTypeIsList(m_type))
    {
        CPLError(CE_Failure, CPLE_NotSupported,
                 "SetMinCount() illegal on scalar argument '%s'",
                 GetName().c_str());
    }
    else
    {
        m_minCount = count;
    }
    return *this;
}

/************************************************************************/
/*                 GDALAlgorithmArgDecl::SetMaxCount()                  */
/************************************************************************/

GDALAlgorithmArgDecl &GDALAlgorithmArgDecl::SetMaxCount(int count)
{
    if (!GDALAlgorithmArgTypeIsList(m_type))
    {
        CPLError(CE_Failure, CPLE_NotSupported,
                 "SetMaxCount() illegal on scalar argument '%s'",
                 GetName().c_str());
    }
    else
    {
        m_maxCount = count;
    }
    return *this;
}

/************************************************************************/
/*                GDALAlgorithmArg::~GDALAlgorithmArg()                 */
/************************************************************************/

GDALAlgorithmArg::~GDALAlgorithmArg() = default;

/************************************************************************/
/*                       GDALAlgorithmArg::Set()                        */
/************************************************************************/

bool GDALAlgorithmArg::Set(bool value)
{
    if (m_decl.GetType() != GAAT_BOOLEAN)
    {
        CPLError(
            CE_Failure, CPLE_AppDefined,
            "Calling Set(bool) on argument '%s' of type %s is not supported",
            GetName().c_str(), GDALAlgorithmArgTypeName(m_decl.GetType()));
        return false;
    }
    return SetInternal(value);
}

bool GDALAlgorithmArg::ProcessString(std::string &value) const
{
    if (m_decl.IsReadFromFileAtSyntaxAllowed() && !value.empty() &&
        value.front() == '@')
    {
        GByte *pabyData = nullptr;
        if (VSIIngestFile(nullptr, value.c_str() + 1, &pabyData, nullptr,
                          10 * 1024 * 1024))
        {
            // Remove UTF-8 BOM
            size_t offset = 0;
            if (pabyData[0] == 0xEF && pabyData[1] == 0xBB &&
                pabyData[2] == 0xBF)
            {
                offset = 3;
            }
            value = reinterpret_cast<const char *>(pabyData + offset);
            VSIFree(pabyData);
        }
        else
        {
            return false;
        }
    }

    if (m_decl.IsRemoveSQLCommentsEnabled())
        value = CPLRemoveSQLComments(value);

    return true;
}

bool GDALAlgorithmArg::Set(const std::string &value)
{
    switch (m_decl.GetType())
    {
        case GAAT_BOOLEAN:
            if (EQUAL(value.c_str(), "1") || EQUAL(value.c_str(), "TRUE") ||
                EQUAL(value.c_str(), "YES") || EQUAL(value.c_str(), "ON"))
            {
                return Set(true);
            }
            else if (EQUAL(value.c_str(), "0") ||
                     EQUAL(value.c_str(), "FALSE") ||
                     EQUAL(value.c_str(), "NO") || EQUAL(value.c_str(), "OFF"))
            {
                return Set(false);
            }
            break;

        case GAAT_INTEGER:
        case GAAT_INTEGER_LIST:
        {
            errno = 0;
            char *endptr = nullptr;
            const auto v = std::strtoll(value.c_str(), &endptr, 10);
            if (errno == 0 && v >= INT_MIN && v <= INT_MAX &&
                endptr == value.c_str() + value.size())
            {
                if (m_decl.GetType() == GAAT_INTEGER)
                    return Set(static_cast<int>(v));
                else
                    return Set(std::vector<int>{static_cast<int>(v)});
            }
            break;
        }

        case GAAT_REAL:
        case GAAT_REAL_LIST:
        {
            char *endptr = nullptr;
            const double v = CPLStrtod(value.c_str(), &endptr);
            if (endptr == value.c_str() + value.size())
            {
                if (m_decl.GetType() == GAAT_REAL)
                    return Set(v);
                else
                    return Set(std::vector<double>{v});
            }
            break;
        }

        case GAAT_STRING:
            break;

        case GAAT_STRING_LIST:
            return Set(std::vector<std::string>{value});

        case GAAT_DATASET:
            return SetDatasetName(value);

        case GAAT_DATASET_LIST:
        {
            std::vector<GDALArgDatasetValue> v;
            v.resize(1);
            v[0].Set(value);
            return Set(std::move(v));
        }
    }

    if (m_decl.GetType() != GAAT_STRING)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "Calling Set(std::string) on argument '%s' of type %s is not "
                 "supported",
                 GetName().c_str(), GDALAlgorithmArgTypeName(m_decl.GetType()));
        return false;
    }

    std::string newValue(value);
    return ProcessString(newValue) && SetInternal(newValue);
}

bool GDALAlgorithmArg::Set(int value)
{
    if (m_decl.GetType() == GAAT_BOOLEAN)
    {
        if (value == 1)
            return Set(true);
        else if (value == 0)
            return Set(false);
    }
    else if (m_decl.GetType() == GAAT_REAL)
    {
        return Set(static_cast<double>(value));
    }
    else if (m_decl.GetType() == GAAT_STRING)
    {
        return Set(std::to_string(value));
    }
    else if (m_decl.GetType() == GAAT_INTEGER_LIST)
    {
        return Set(std::vector<int>{value});
    }
    else if (m_decl.GetType() == GAAT_REAL_LIST)
    {
        return Set(std::vector<double>{static_cast<double>(value)});
    }
    else if (m_decl.GetType() == GAAT_STRING_LIST)
    {
        return Set(std::vector<std::string>{std::to_string(value)});
    }

    if (m_decl.GetType() != GAAT_INTEGER)
    {
        CPLError(
            CE_Failure, CPLE_AppDefined,
            "Calling Set(int) on argument '%s' of type %s is not supported",
            GetName().c_str(), GDALAlgorithmArgTypeName(m_decl.GetType()));
        return false;
    }
    return SetInternal(value);
}

bool GDALAlgorithmArg::Set(double value)
{
    if (m_decl.GetType() == GAAT_INTEGER && value >= INT_MIN &&
        value <= INT_MAX && static_cast<int>(value) == value)
    {
        return Set(static_cast<int>(value));
    }
    else if (m_decl.GetType() == GAAT_STRING)
    {
        return Set(std::to_string(value));
    }
    else if (m_decl.GetType() == GAAT_INTEGER_LIST && value >= INT_MIN &&
             value <= INT_MAX && static_cast<int>(value) == value)
    {
        return Set(std::vector<int>{static_cast<int>(value)});
    }
    else if (m_decl.GetType() == GAAT_REAL_LIST)
    {
        return Set(std::vector<double>{value});
    }
    else if (m_decl.GetType() == GAAT_STRING_LIST)
    {
        return Set(std::vector<std::string>{std::to_string(value)});
    }
    else if (m_decl.GetType() != GAAT_REAL)
    {
        CPLError(
            CE_Failure, CPLE_AppDefined,
            "Calling Set(double) on argument '%s' of type %s is not supported",
            GetName().c_str(), GDALAlgorithmArgTypeName(m_decl.GetType()));
        return false;
    }
    return SetInternal(value);
}

static bool CheckCanSetDatasetObject(const GDALAlgorithmArg *arg)
{
    if (arg->IsOutput() && arg->GetDatasetInputFlags() == GADV_NAME &&
        arg->GetDatasetOutputFlags() == GADV_OBJECT)
    {
        CPLError(
            CE_Failure, CPLE_AppDefined,
            "Dataset object '%s' is created by algorithm and cannot be set "
            "as an input.",
            arg->GetName().c_str());
        return false;
    }
    else if ((arg->GetDatasetInputFlags() & GADV_OBJECT) == 0)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "Dataset%s '%s' must be provided by name, not as object.",
                 arg->GetMaxCount() > 1 ? "s" : "", arg->GetName().c_str());
        return false;
    }

    return true;
}

bool GDALAlgorithmArg::Set(GDALDataset *ds)
{
    if (m_decl.GetType() != GAAT_DATASET)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "Calling Set(GDALDataset*, bool) on argument '%s' of type %s "
                 "is not supported",
                 GetName().c_str(), GDALAlgorithmArgTypeName(m_decl.GetType()));
        return false;
    }
    if (!CheckCanSetDatasetObject(this))
        return false;
    m_explicitlySet = true;
    auto &val = *std::get<GDALArgDatasetValue *>(m_value);
    val.Set(ds);
    return RunAllActions();
}

bool GDALAlgorithmArg::Set(std::unique_ptr<GDALDataset> ds)
{
    if (m_decl.GetType() != GAAT_DATASET)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "Calling Set(GDALDataset*, bool) on argument '%s' of type %s "
                 "is not supported",
                 GetName().c_str(), GDALAlgorithmArgTypeName(m_decl.GetType()));
        return false;
    }
    if (!CheckCanSetDatasetObject(this))
        return false;
    m_explicitlySet = true;
    auto &val = *std::get<GDALArgDatasetValue *>(m_value);
    val.Set(std::move(ds));
    return RunAllActions();
}

bool GDALAlgorithmArg::SetDatasetName(const std::string &name)
{
    if (m_decl.GetType() != GAAT_DATASET)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "Calling SetDatasetName() on argument '%s' of type %s is "
                 "not supported",
                 GetName().c_str(), GDALAlgorithmArgTypeName(m_decl.GetType()));
        return false;
    }
    m_explicitlySet = true;
    std::get<GDALArgDatasetValue *>(m_value)->Set(name);
    return RunAllActions();
}

bool GDALAlgorithmArg::SetFrom(const GDALArgDatasetValue &other)
{
    if (m_decl.GetType() != GAAT_DATASET)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "Calling SetFrom() on argument '%s' of type %s is "
                 "not supported",
                 GetName().c_str(), GDALAlgorithmArgTypeName(m_decl.GetType()));
        return false;
    }
    if (other.GetDatasetRef() && !CheckCanSetDatasetObject(this))
        return false;
    m_explicitlySet = true;
    std::get<GDALArgDatasetValue *>(m_value)->SetFrom(other);
    return RunAllActions();
}

bool GDALAlgorithmArg::Set(const std::vector<std::string> &value)
{
    if (m_decl.GetType() == GAAT_INTEGER_LIST)
    {
        std::vector<int> v_i;
        for (const std::string &s : value)
        {
            errno = 0;
            char *endptr = nullptr;
            const auto v = std::strtoll(s.c_str(), &endptr, 10);
            if (errno == 0 && v >= INT_MIN && v <= INT_MAX &&
                endptr == s.c_str() + s.size())
            {
                v_i.push_back(static_cast<int>(v));
            }
            else
            {
                break;
            }
        }
        if (v_i.size() == value.size())
            return Set(v_i);
    }
    else if (m_decl.GetType() == GAAT_REAL_LIST)
    {
        std::vector<double> v_d;
        for (const std::string &s : value)
        {
            char *endptr = nullptr;
            const double v = CPLStrtod(s.c_str(), &endptr);
            if (endptr == s.c_str() + s.size())
            {
                v_d.push_back(v);
            }
            else
            {
                break;
            }
        }
        if (v_d.size() == value.size())
            return Set(v_d);
    }
    else if ((m_decl.GetType() == GAAT_INTEGER ||
              m_decl.GetType() == GAAT_REAL ||
              m_decl.GetType() == GAAT_STRING) &&
             value.size() == 1)
    {
        return Set(value[0]);
    }
    else if (m_decl.GetType() == GAAT_DATASET_LIST)
    {
        std::vector<GDALArgDatasetValue> dsVector;
        for (const std::string &s : value)
            dsVector.emplace_back(s);
        return Set(std::move(dsVector));
    }

    if (m_decl.GetType() != GAAT_STRING_LIST)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "Calling Set(const std::vector<std::string> &) on argument "
                 "'%s' of type %s is not supported",
                 GetName().c_str(), GDALAlgorithmArgTypeName(m_decl.GetType()));
        return false;
    }

    if (m_decl.IsReadFromFileAtSyntaxAllowed() ||
        m_decl.IsRemoveSQLCommentsEnabled())
    {
        std::vector<std::string> newValue(value);
        for (auto &s : newValue)
        {
            if (!ProcessString(s))
                return false;
        }
        return SetInternal(newValue);
    }
    else
    {
        return SetInternal(value);
    }
}

bool GDALAlgorithmArg::Set(const std::vector<int> &value)
{
    if (m_decl.GetType() == GAAT_REAL_LIST)
    {
        std::vector<double> v_d;
        for (int i : value)
            v_d.push_back(i);
        return Set(v_d);
    }
    else if (m_decl.GetType() == GAAT_STRING_LIST)
    {
        std::vector<std::string> v_s;
        for (int i : value)
            v_s.push_back(std::to_string(i));
        return Set(v_s);
    }
    else if ((m_decl.GetType() == GAAT_INTEGER ||
              m_decl.GetType() == GAAT_REAL ||
              m_decl.GetType() == GAAT_STRING) &&
             value.size() == 1)
    {
        return Set(value[0]);
    }

    if (m_decl.GetType() != GAAT_INTEGER_LIST)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "Calling Set(const std::vector<int> &) on argument '%s' of "
                 "type %s is not supported",
                 GetName().c_str(), GDALAlgorithmArgTypeName(m_decl.GetType()));
        return false;
    }
    return SetInternal(value);
}

bool GDALAlgorithmArg::Set(const std::vector<double> &value)
{
    if (m_decl.GetType() == GAAT_INTEGER_LIST)
    {
        std::vector<int> v_i;
        for (double d : value)
        {
            if (d >= INT_MIN && d <= INT_MAX && static_cast<int>(d) == d)
            {
                v_i.push_back(static_cast<int>(d));
            }
            else
            {
                break;
            }
        }
        if (v_i.size() == value.size())
            return Set(v_i);
    }
    else if (m_decl.GetType() == GAAT_STRING_LIST)
    {
        std::vector<std::string> v_s;
        for (double d : value)
            v_s.push_back(std::to_string(d));
        return Set(v_s);
    }
    else if ((m_decl.GetType() == GAAT_INTEGER ||
              m_decl.GetType() == GAAT_REAL ||
              m_decl.GetType() == GAAT_STRING) &&
             value.size() == 1)
    {
        return Set(value[0]);
    }

    if (m_decl.GetType() != GAAT_REAL_LIST)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "Calling Set(const std::vector<double> &) on argument '%s' of "
                 "type %s is not supported",
                 GetName().c_str(), GDALAlgorithmArgTypeName(m_decl.GetType()));
        return false;
    }
    return SetInternal(value);
}

bool GDALAlgorithmArg::Set(std::vector<GDALArgDatasetValue> &&value)
{
    if (m_decl.GetType() != GAAT_DATASET_LIST)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "Calling Set(const std::vector<GDALArgDatasetValue> &&) on "
                 "argument '%s' of type %s is not supported",
                 GetName().c_str(), GDALAlgorithmArgTypeName(m_decl.GetType()));
        return false;
    }
    m_explicitlySet = true;
    *std::get<std::vector<GDALArgDatasetValue> *>(m_value) = std::move(value);
    return RunAllActions();
}

GDALAlgorithmArg &
GDALAlgorithmArg::operator=(std::unique_ptr<GDALDataset> value)
{
    Set(std::move(value));
    return *this;
}

bool GDALAlgorithmArg::Set(const OGRSpatialReference &value)
{
    const char *const apszOptions[] = {"FORMAT=WKT2_2019", nullptr};
    return Set(value.exportToWkt(apszOptions));
}

bool GDALAlgorithmArg::SetFrom(const GDALAlgorithmArg &other)
{
    if (m_decl.GetType() != other.GetType())
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "Calling SetFrom() on argument '%s' of type %s whereas "
                 "other argument type is %s is not supported",
                 GetName().c_str(), GDALAlgorithmArgTypeName(m_decl.GetType()),
                 GDALAlgorithmArgTypeName(other.GetType()));
        return false;
    }

    switch (m_decl.GetType())
    {
        case GAAT_BOOLEAN:
            *std::get<bool *>(m_value) = *std::get<bool *>(other.m_value);
            break;
        case GAAT_STRING:
            *std::get<std::string *>(m_value) =
                *std::get<std::string *>(other.m_value);
            break;
        case GAAT_INTEGER:
            *std::get<int *>(m_value) = *std::get<int *>(other.m_value);
            break;
        case GAAT_REAL:
            *std::get<double *>(m_value) = *std::get<double *>(other.m_value);
            break;
        case GAAT_DATASET:
            return SetFrom(other.Get<GDALArgDatasetValue>());
        case GAAT_STRING_LIST:
            *std::get<std::vector<std::string> *>(m_value) =
                *std::get<std::vector<std::string> *>(other.m_value);
            break;
        case GAAT_INTEGER_LIST:
            *std::get<std::vector<int> *>(m_value) =
                *std::get<std::vector<int> *>(other.m_value);
            break;
        case GAAT_REAL_LIST:
            *std::get<std::vector<double> *>(m_value) =
                *std::get<std::vector<double> *>(other.m_value);
            break;
        case GAAT_DATASET_LIST:
        {
            std::get<std::vector<GDALArgDatasetValue> *>(m_value)->clear();
            for (const auto &val :
                 *std::get<std::vector<GDALArgDatasetValue> *>(other.m_value))
            {
                GDALArgDatasetValue v;
                v.SetFrom(val);
                std::get<std::vector<GDALArgDatasetValue> *>(m_value)
                    ->push_back(std::move(v));
            }
            break;
        }
    }
    m_explicitlySet = true;
    return RunAllActions();
}

/************************************************************************/
/*                  GDALAlgorithmArg::RunAllActions()                   */
/************************************************************************/

bool GDALAlgorithmArg::RunAllActions()
{
    if (!RunValidationActions())
        return false;
    RunActions();
    return true;
}

/************************************************************************/
/*                    GDALAlgorithmArg::RunActions()                    */
/************************************************************************/

void GDALAlgorithmArg::RunActions()
{
    for (const auto &f : m_actions)
        f();
}

/************************************************************************/
/*                  GDALAlgorithmArg::ValidateChoice()                  */
/************************************************************************/

// Returns the canonical value if matching a valid choice, or empty string
// otherwise.
std::string GDALAlgorithmArg::ValidateChoice(const std::string &value) const
{
    for (const std::string &choice : GetChoices())
    {
        if (EQUAL(value.c_str(), choice.c_str()))
        {
            return choice;
        }
    }

    for (const std::string &choice : GetHiddenChoices())
    {
        if (EQUAL(value.c_str(), choice.c_str()))
        {
            return choice;
        }
    }

    std::string expected;
    for (const auto &choice : GetChoices())
    {
        if (!expected.empty())
            expected += ", ";
        expected += '\'';
        expected += choice;
        expected += '\'';
    }
    if (m_owner && m_owner->IsCalledFromCommandLine() && value == "?")
    {
        return "?";
    }
    CPLError(CE_Failure, CPLE_IllegalArg,
             "Invalid value '%s' for string argument '%s'. Should be "
             "one among %s.",
             value.c_str(), GetName().c_str(), expected.c_str());
    return std::string();
}

/************************************************************************/
/*                 GDALAlgorithmArg::ValidateIntRange()                 */
/************************************************************************/

bool GDALAlgorithmArg::ValidateIntRange(int val) const
{
    bool ret = true;

    const auto [minVal, minValIsIncluded] = GetMinValue();
    if (!std::isnan(minVal))
    {
        if (minValIsIncluded && val < minVal)
        {
            CPLError(CE_Failure, CPLE_IllegalArg,
                     "Value of argument '%s' is %d, but should be >= %d",
                     GetName().c_str(), val, static_cast<int>(minVal));
            ret = false;
        }
        else if (!minValIsIncluded && val <= minVal)
        {
            CPLError(CE_Failure, CPLE_IllegalArg,
                     "Value of argument '%s' is %d, but should be > %d",
                     GetName().c_str(), val, static_cast<int>(minVal));
            ret = false;
        }
    }

    const auto [maxVal, maxValIsIncluded] = GetMaxValue();
    if (!std::isnan(maxVal))
    {

        if (maxValIsIncluded && val > maxVal)
        {
            CPLError(CE_Failure, CPLE_IllegalArg,
                     "Value of argument '%s' is %d, but should be <= %d",
                     GetName().c_str(), val, static_cast<int>(maxVal));
            ret = false;
        }
        else if (!maxValIsIncluded && val >= maxVal)
        {
            CPLError(CE_Failure, CPLE_IllegalArg,
                     "Value of argument '%s' is %d, but should be < %d",
                     GetName().c_str(), val, static_cast<int>(maxVal));
            ret = false;
        }
    }

    return ret;
}

/************************************************************************/
/*                GDALAlgorithmArg::ValidateRealRange()                 */
/************************************************************************/

bool GDALAlgorithmArg::ValidateRealRange(double val) const
{
    bool ret = true;

    const auto [minVal, minValIsIncluded] = GetMinValue();
    if (!std::isnan(minVal))
    {
        if (minValIsIncluded && !(val >= minVal))
        {
            CPLError(CE_Failure, CPLE_IllegalArg,
                     "Value of argument '%s' is %g, but should be >= %g",
                     GetName().c_str(), val, minVal);
            ret = false;
        }
        else if (!minValIsIncluded && !(val > minVal))
        {
            CPLError(CE_Failure, CPLE_IllegalArg,
                     "Value of argument '%s' is %g, but should be > %g",
                     GetName().c_str(), val, minVal);
            ret = false;
        }
    }

    const auto [maxVal, maxValIsIncluded] = GetMaxValue();
    if (!std::isnan(maxVal))
    {

        if (maxValIsIncluded && !(val <= maxVal))
        {
            CPLError(CE_Failure, CPLE_IllegalArg,
                     "Value of argument '%s' is %g, but should be <= %g",
                     GetName().c_str(), val, maxVal);
            ret = false;
        }
        else if (!maxValIsIncluded && !(val < maxVal))
        {
            CPLError(CE_Failure, CPLE_IllegalArg,
                     "Value of argument '%s' is %g, but should be < %g",
                     GetName().c_str(), val, maxVal);
            ret = false;
        }
    }

    return ret;
}

/************************************************************************/
/*                        CheckDuplicateValues()                        */
/************************************************************************/

template <class T>
static bool CheckDuplicateValues(const GDALAlgorithmArg *arg,
                                 const std::vector<T> &values)
{
    auto tmpValues = values;
    bool bHasDupValues = false;
    if constexpr (std::is_floating_point_v<T>)
    {
        // Avoid undefined behavior with NaN values
        std::sort(tmpValues.begin(), tmpValues.end(),
                  [](T a, T b)
                  {
                      if (std::isnan(a) && !std::isnan(b))
                          return true;
                      if (std::isnan(b))
                          return false;
                      return a < b;
                  });

        bHasDupValues =
            std::adjacent_find(tmpValues.begin(), tmpValues.end(),
                               [](T a, T b)
                               {
                                   if (std::isnan(a) && std::isnan(b))
                                       return true;
                                   return a == b;
                               }) != tmpValues.end();
    }
    else
    {
        std::sort(tmpValues.begin(), tmpValues.end());
        bHasDupValues = std::adjacent_find(tmpValues.begin(),
                                           tmpValues.end()) != tmpValues.end();
    }
    if (bHasDupValues)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "'%s' must be a list of unique values.",
                 arg->GetName().c_str());
        return false;
    }
    return true;
}

/************************************************************************/
/*               GDALAlgorithmArg::RunValidationActions()               */
/************************************************************************/

bool GDALAlgorithmArg::RunValidationActions()
{
    bool ret = true;

    if (GetType() == GAAT_STRING && !GetChoices().empty())
    {
        auto &val = Get<std::string>();
        std::string validVal = ValidateChoice(val);
        if (validVal.empty())
            ret = false;
        else
            val = std::move(validVal);
    }
    else if (GetType() == GAAT_STRING_LIST && !GetChoices().empty())
    {
        auto &values = Get<std::vector<std::string>>();
        for (std::string &val : values)
        {
            std::string validVal = ValidateChoice(val);
            if (validVal.empty())
                ret = false;
            else
                val = std::move(validVal);
        }
    }

    const auto CheckMinCharCount =
        [this, &ret](const std::string &val, int nMinCharCount)
    {
        if (val.size() < static_cast<size_t>(nMinCharCount))
        {
            CPLError(CE_Failure, CPLE_IllegalArg,
                     "Value of argument '%s' is '%s', but should have at least "
                     "%d character%s",
                     GetName().c_str(), val.c_str(), nMinCharCount,
                     nMinCharCount > 1 ? "s" : "");
            ret = false;
        }
    };

    const auto CheckMaxCharCount =
        [this, &ret](const std::string &val, int nMaxCharCount)
    {
        if (val.size() > static_cast<size_t>(nMaxCharCount))
        {
            CPLError(
                CE_Failure, CPLE_IllegalArg,
                "Value of argument '%s' is '%s', but should have no more than "
                "%d character%s",
                GetName().c_str(), val.c_str(), nMaxCharCount,
                nMaxCharCount > 1 ? "s" : "");
            ret = false;
        }
    };

    switch (GetType())
    {
        case GAAT_BOOLEAN:
            break;

        case GAAT_STRING:
        {
            const auto &val = Get<std::string>();
            const int nMinCharCount = GetMinCharCount();
            if (nMinCharCount > 0)
            {
                CheckMinCharCount(val, nMinCharCount);
            }

            const int nMaxCharCount = GetMaxCharCount();
            CheckMaxCharCount(val, nMaxCharCount);
            break;
        }

        case GAAT_STRING_LIST:
        {
            const int nMinCharCount = GetMinCharCount();
            const int nMaxCharCount = GetMaxCharCount();
            const auto &values = Get<std::vector<std::string>>();
            for (const auto &val : values)
            {
                if (nMinCharCount > 0)
                    CheckMinCharCount(val, nMinCharCount);
                CheckMaxCharCount(val, nMaxCharCount);
            }

            if (!GetDuplicateValuesAllowed() &&
                !CheckDuplicateValues(this, values))
                ret = false;
            break;
        }

        case GAAT_INTEGER:
        {
            ret = ValidateIntRange(Get<int>()) && ret;
            break;
        }

        case GAAT_INTEGER_LIST:
        {
            const auto &values = Get<std::vector<int>>();
            for (int v : values)
                ret = ValidateIntRange(v) && ret;

            if (!GetDuplicateValuesAllowed() &&
                !CheckDuplicateValues(this, values))
                ret = false;
            break;
        }

        case GAAT_REAL:
        {
            ret = ValidateRealRange(Get<double>()) && ret;
            break;
        }

        case GAAT_REAL_LIST:
        {
            const auto &values = Get<std::vector<double>>();
            for (double v : values)
                ret = ValidateRealRange(v) && ret;

            if (!GetDuplicateValuesAllowed() &&
                !CheckDuplicateValues(this, values))
                ret = false;
            break;
        }

        case GAAT_DATASET:
            break;

        case GAAT_DATASET_LIST:
        {
            if (!GetDuplicateValuesAllowed())
            {
                const auto &values = Get<std::vector<GDALArgDatasetValue>>();
                std::vector<std::string> aosValues;
                for (const auto &v : values)
                {
                    const GDALDataset *poDS = v.GetDatasetRef();
                    if (poDS)
                    {
                        auto poDriver = poDS->GetDriver();
                        // The dataset name for a MEM driver is not relevant,
                        // so use the pointer address
                        if ((poDriver &&
                             EQUAL(poDriver->GetDescription(), "MEM")) ||
                            poDS->GetDescription()[0] == 0)
                        {
                            aosValues.push_back(CPLSPrintf("%p", poDS));
                        }
                        else
                        {
                            aosValues.push_back(poDS->GetDescription());
                        }
                    }
                    else
                    {
                        aosValues.push_back(v.GetName());
                    }
                }
                if (!CheckDuplicateValues(this, aosValues))
                    ret = false;
            }
            break;
        }
    }

    if (GDALAlgorithmArgTypeIsList(GetType()))
    {
        int valueCount = 0;
        if (GetType() == GAAT_STRING_LIST)
        {
            valueCount =
                static_cast<int>(Get<std::vector<std::string>>().size());
        }
        else if (GetType() == GAAT_INTEGER_LIST)
        {
            valueCount = static_cast<int>(Get<std::vector<int>>().size());
        }
        else if (GetType() == GAAT_REAL_LIST)
        {
            valueCount = static_cast<int>(Get<std::vector<double>>().size());
        }
        else if (GetType() == GAAT_DATASET_LIST)
        {
            valueCount = static_cast<int>(
                Get<std::vector<GDALArgDatasetValue>>().size());
        }

        if (valueCount != GetMinCount() && GetMinCount() == GetMaxCount())
        {
            ReportError(CE_Failure, CPLE_AppDefined,
                        "%d value%s been specified for argument '%s', "
                        "whereas exactly %d %s expected.",
                        valueCount, valueCount > 1 ? "s have" : " has",
                        GetName().c_str(), GetMinCount(),
                        GetMinCount() > 1 ? "were" : "was");
            ret = false;
        }
        else if (valueCount < GetMinCount())
        {
            ReportError(CE_Failure, CPLE_AppDefined,
                        "Only %d value%s been specified for argument '%s', "
                        "whereas at least %d %s expected.",
                        valueCount, valueCount > 1 ? "s have" : " has",
                        GetName().c_str(), GetMinCount(),
                        GetMinCount() > 1 ? "were" : "was");
            ret = false;
        }
        else if (valueCount > GetMaxCount())
        {
            ReportError(CE_Failure, CPLE_AppDefined,
                        "%d value%s been specified for argument '%s', "
                        "whereas at most %d %s expected.",
                        valueCount, valueCount > 1 ? "s have" : " has",
                        GetName().c_str(), GetMaxCount(),
                        GetMaxCount() > 1 ? "were" : "was");
            ret = false;
        }
    }

    if (ret)
    {
        for (const auto &f : m_validationActions)
        {
            if (!f())
                ret = false;
        }
    }

    return ret;
}

/************************************************************************/
/*                   GDALAlgorithmArg::ReportError()                    */
/************************************************************************/

void GDALAlgorithmArg::ReportError(CPLErr eErrClass, CPLErrorNum err_no,
                                   const char *fmt, ...) const
{
    va_list args;
    va_start(args, fmt);
    if (m_owner)
    {
        m_owner->ReportError(eErrClass, err_no, "%s",
                             CPLString().vPrintf(fmt, args).c_str());
    }
    else
    {
        CPLError(eErrClass, err_no, "%s",
                 CPLString().vPrintf(fmt, args).c_str());
    }
    va_end(args);
}

/************************************************************************/
/*                 GDALAlgorithmArg::GetEscapedString()                 */
/************************************************************************/

/* static */
std::string GDALAlgorithmArg::GetEscapedString(const std::string &s)
{
    if (s.find_first_of("\" \\,") != std::string::npos &&
        !(s.size() > 4 &&
          s[0] == GDALAbstractPipelineAlgorithm::OPEN_NESTED_PIPELINE[0] &&
          s[1] == ' ' && s[s.size() - 2] == ' ' &&
          s.back() == GDALAbstractPipelineAlgorithm::CLOSE_NESTED_PIPELINE[0]))
    {
        return std::string("\"")
            .append(
                CPLString(s).replaceAll('\\', "\\\\").replaceAll('"', "\\\""))
            .append("\"");
    }
    else
    {
        return s;
    }
}

/************************************************************************/
/*                    GDALAlgorithmArg::Serialize()                     */
/************************************************************************/

bool GDALAlgorithmArg::Serialize(std::string &serializedArg,
                                 bool absolutePath) const
{
    serializedArg.clear();

    if (!IsExplicitlySet())
    {
        return false;
    }

    std::string ret = "--";
    ret += GetName();
    if (GetType() == GAAT_BOOLEAN)
    {
        serializedArg = std::move(ret);
        return true;
    }

    const auto AddListValueSeparator = [this, &ret]()
    {
        if (GetPackedValuesAllowed())
        {
            ret += ',';
        }
        else
        {
            ret += " --";
            ret += GetName();
            ret += ' ';
        }
    };

    const auto MakeAbsolutePath = [](const std::string &filename)
    {
        VSIStatBufL sStat;
        if (VSIStatL(filename.c_str(), &sStat) != 0 ||
            !CPLIsFilenameRelative(filename.c_str()))
            return filename;
        char *pszCWD = CPLGetCurrentDir();
        if (!pszCWD)
            return filename;
        const auto absPath =
            CPLFormFilenameSafe(pszCWD, filename.c_str(), nullptr);
        CPLFree(pszCWD);
        return absPath;
    };

    ret += ' ';
    switch (GetType())
    {
        case GAAT_BOOLEAN:
            break;
        case GAAT_STRING:
        {
            const auto &val = Get<std::string>();
            ret += GetEscapedString(val);
            break;
        }
        case GAAT_INTEGER:
        {
            ret += CPLSPrintf("%d", Get<int>());
            break;
        }
        case GAAT_REAL:
        {
            ret += CPLSPrintf("%.17g", Get<double>());
            break;
        }
        case GAAT_DATASET:
        {
            const auto &val = Get<GDALArgDatasetValue>();
            const auto &str = val.GetName();
            if (str.empty())
            {
                return false;
            }
            ret += GetEscapedString(absolutePath ? MakeAbsolutePath(str) : str);
            break;
        }
        case GAAT_STRING_LIST:
        {
            const auto &vals = Get<std::vector<std::string>>();
            for (size_t i = 0; i < vals.size(); ++i)
            {
                if (i > 0)
                    AddListValueSeparator();
                ret += GetEscapedString(vals[i]);
            }
            break;
        }
        case GAAT_INTEGER_LIST:
        {
            const auto &vals = Get<std::vector<int>>();
            for (size_t i = 0; i < vals.size(); ++i)
            {
                if (i > 0)
                    AddListValueSeparator();
                ret += CPLSPrintf("%d", vals[i]);
            }
            break;
        }
        case GAAT_REAL_LIST:
        {
            const auto &vals = Get<std::vector<double>>();
            for (size_t i = 0; i < vals.size(); ++i)
            {
                if (i > 0)
                    AddListValueSeparator();
                ret += CPLSPrintf("%.17g", vals[i]);
            }
            break;
        }
        case GAAT_DATASET_LIST:
        {
            const auto &vals = Get<std::vector<GDALArgDatasetValue>>();
            for (size_t i = 0; i < vals.size(); ++i)
            {
                if (i > 0)
                    AddListValueSeparator();
                const auto &val = vals[i];
                const auto &str = val.GetName();
                if (str.empty())
                {
                    return false;
                }
                ret += GetEscapedString(absolutePath ? MakeAbsolutePath(str)
                                                     : str);
            }
            break;
        }
    }

    serializedArg = std::move(ret);
    return true;
}

/************************************************************************/
/*                  ~GDALInConstructionAlgorithmArg()                   */
/************************************************************************/

GDALInConstructionAlgorithmArg::~GDALInConstructionAlgorithmArg() = default;

/************************************************************************/
/*              GDALInConstructionAlgorithmArg::AddAlias()              */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALInConstructionAlgorithmArg::AddAlias(const std::string &alias)
{
    m_decl.AddAlias(alias);
    if (m_owner)
        m_owner->AddAliasFor(this, alias);
    return *this;
}

/************************************************************************/
/*           GDALInConstructionAlgorithmArg::AddHiddenAlias()           */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALInConstructionAlgorithmArg::AddHiddenAlias(const std::string &alias)
{
    m_decl.AddHiddenAlias(alias);
    if (m_owner)
        m_owner->AddAliasFor(this, alias);
    return *this;
}

/************************************************************************/
/*         GDALInConstructionAlgorithmArg::AddShortNameAlias()          */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALInConstructionAlgorithmArg::AddShortNameAlias(char shortNameAlias)
{
    m_decl.AddShortNameAlias(shortNameAlias);
    if (m_owner)
        m_owner->AddShortNameAliasFor(this, shortNameAlias);
    return *this;
}

/************************************************************************/
/*           GDALInConstructionAlgorithmArg::SetPositional()            */
/************************************************************************/

GDALInConstructionAlgorithmArg &GDALInConstructionAlgorithmArg::SetPositional()
{
    m_decl.SetPositional();
    if (m_owner)
        m_owner->SetPositional(this);
    return *this;
}

/************************************************************************/
/*              GDALArgDatasetValue::GDALArgDatasetValue()              */
/************************************************************************/

GDALArgDatasetValue::GDALArgDatasetValue(GDALDataset *poDS)
    : m_poDS(poDS), m_name(m_poDS ? m_poDS->GetDescription() : std::string()),
      m_nameSet(true)
{
    if (m_poDS)
        m_poDS->Reference();
}

/************************************************************************/
/*                      GDALArgDatasetValue::Set()                      */
/************************************************************************/

void GDALArgDatasetValue::Set(const std::string &name)
{
    Close();
    m_name = name;
    m_nameSet = true;
    if (m_ownerArg)
        m_ownerArg->NotifyValueSet();
}

/************************************************************************/
/*                      GDALArgDatasetValue::Set()                      */
/************************************************************************/

void GDALArgDatasetValue::Set(std::unique_ptr<GDALDataset> poDS)
{
    Close();
    m_poDS = poDS.release();
    m_name = m_poDS ? m_poDS->GetDescription() : std::string();
    m_nameSet = true;
    if (m_ownerArg)
        m_ownerArg->NotifyValueSet();
}

/************************************************************************/
/*                      GDALArgDatasetValue::Set()                      */
/************************************************************************/

void GDALArgDatasetValue::Set(GDALDataset *poDS)
{
    Close();
    m_poDS = poDS;
    if (m_poDS)
        m_poDS->Reference();
    m_name = m_poDS ? m_poDS->GetDescription() : std::string();
    m_nameSet = true;
    if (m_ownerArg)
        m_ownerArg->NotifyValueSet();
}

/************************************************************************/
/*                    GDALArgDatasetValue::SetFrom()                    */
/************************************************************************/

void GDALArgDatasetValue::SetFrom(const GDALArgDatasetValue &other)
{
    Close();
    m_name = other.m_name;
    m_nameSet = other.m_nameSet;
    m_poDS = other.m_poDS;
    if (m_poDS)
        m_poDS->Reference();
}

/************************************************************************/
/*             GDALArgDatasetValue::~GDALArgDatasetValue()              */
/************************************************************************/

GDALArgDatasetValue::~GDALArgDatasetValue()
{
    Close();
}

/************************************************************************/
/*                     GDALArgDatasetValue::Close()                     */
/************************************************************************/

bool GDALArgDatasetValue::Close()
{
    bool ret = true;
    if (m_poDS && m_poDS->Dereference() == 0)
    {
        ret = m_poDS->Close() == CE_None;
        delete m_poDS;
    }
    m_poDS = nullptr;
    return ret;
}

/************************************************************************/
/*                   GDALArgDatasetValue::operator=()                   */
/************************************************************************/

GDALArgDatasetValue &GDALArgDatasetValue::operator=(GDALArgDatasetValue &&other)
{
    Close();
    m_poDS = other.m_poDS;
    m_name = other.m_name;
    m_nameSet = other.m_nameSet;
    other.m_poDS = nullptr;
    other.m_name.clear();
    other.m_nameSet = false;
    return *this;
}

/************************************************************************/
/*                  GDALArgDatasetValue::GetDataset()                   */
/************************************************************************/

GDALDataset *GDALArgDatasetValue::GetDatasetIncreaseRefCount()
{
    if (m_poDS)
        m_poDS->Reference();
    return m_poDS;
}

/************************************************************************/
/*           GDALArgDatasetValue(GDALArgDatasetValue &&other)           */
/************************************************************************/

GDALArgDatasetValue::GDALArgDatasetValue(GDALArgDatasetValue &&other)
    : m_poDS(other.m_poDS), m_name(other.m_name), m_nameSet(other.m_nameSet)
{
    other.m_poDS = nullptr;
    other.m_name.clear();
}

/************************************************************************/
/*            GDALInConstructionAlgorithmArg::SetIsCRSArg()             */
/************************************************************************/

GDALInConstructionAlgorithmArg &GDALInConstructionAlgorithmArg::SetIsCRSArg(
    bool noneAllowed, const std::vector<std::string> &specialValues)
{
    if (GetType() != GAAT_STRING)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "SetIsCRSArg() can only be called on a String argument");
        return *this;
    }
    AddValidationAction(
        [this, noneAllowed, specialValues]()
        {
            const std::string &osVal =
                static_cast<const GDALInConstructionAlgorithmArg *>(this)
                    ->Get<std::string>();
            if (osVal == "?" && m_owner && m_owner->IsCalledFromCommandLine())
                return true;

            if ((!noneAllowed || (osVal != "none" && osVal != "null")) &&
                std::find(specialValues.begin(), specialValues.end(), osVal) ==
                    specialValues.end())
            {
                OGRSpatialReference oSRS;
                if (oSRS.SetFromUserInput(osVal.c_str()) != OGRERR_NONE)
                {
                    m_owner->ReportError(CE_Failure, CPLE_AppDefined,
                                         "Invalid value for '%s' argument",
                                         GetName().c_str());
                    return false;
                }
            }
            return true;
        });

    SetAutoCompleteFunction(
        [this, noneAllowed, specialValues](const std::string &currentValue)
        {
            bool bIsRaster = false;
            OGREnvelope sDatasetLongLatEnv;
            std::string osCelestialBodyName;
            if (GetName() == GDAL_ARG_NAME_OUTPUT_CRS)
            {
                auto inputArg = m_owner->GetArg(GDAL_ARG_NAME_INPUT);
                if (inputArg && inputArg->GetType() == GAAT_DATASET_LIST)
                {
                    auto &val =
                        inputArg->Get<std::vector<GDALArgDatasetValue>>();
                    if (val.size() == 1)
                    {
                        CPLErrorStateBackuper oBackuper(CPLQuietErrorHandler);
                        auto poDS = std::unique_ptr<GDALDataset>(
                            GDALDataset::Open(val[0].GetName().c_str()));
                        if (poDS)
                        {
                            bIsRaster = poDS->GetRasterCount() != 0;
                            if (auto poCRS = poDS->GetSpatialRef())
                            {
                                const char *pszCelestialBodyName =
                                    poCRS->GetCelestialBodyName();
                                if (pszCelestialBodyName)
                                    osCelestialBodyName = pszCelestialBodyName;

                                if (!pszCelestialBodyName ||
                                    !EQUAL(pszCelestialBodyName, "Earth"))
                                {
                                    OGRSpatialReference oLongLat;
                                    oLongLat.CopyGeogCSFrom(poCRS);
                                    oLongLat.SetAxisMappingStrategy(
                                        OAMS_TRADITIONAL_GIS_ORDER);
                                    poDS->GetExtent(&sDatasetLongLatEnv,
                                                    &oLongLat);
                                }
                                else
                                {
                                    poDS->GetExtentWGS84LongLat(
                                        &sDatasetLongLatEnv);
                                }
                            }
                        }
                    }
                }
            }

            const auto IsCRSCompatible =
                [bIsRaster, &sDatasetLongLatEnv,
                 &osCelestialBodyName](const OSRCRSInfo *crsInfo)
            {
                if (!sDatasetLongLatEnv.IsInit())
                    return true;
                return crsInfo->eType != OSR_CRS_TYPE_VERTICAL &&
                       !(bIsRaster &&
                         crsInfo->eType == OSR_CRS_TYPE_GEOCENTRIC) &&
                       crsInfo->dfWestLongitudeDeg <
                           crsInfo->dfEastLongitudeDeg &&
                       sDatasetLongLatEnv.MinX < crsInfo->dfEastLongitudeDeg &&
                       sDatasetLongLatEnv.MaxX > crsInfo->dfWestLongitudeDeg &&
                       sDatasetLongLatEnv.MinY < crsInfo->dfNorthLatitudeDeg &&
                       sDatasetLongLatEnv.MaxY > crsInfo->dfSouthLatitudeDeg &&
                       ((!osCelestialBodyName.empty() &&
                         crsInfo->pszCelestialBodyName &&
                         osCelestialBodyName ==
                             crsInfo->pszCelestialBodyName) ||
                        (osCelestialBodyName.empty() &&
                         !crsInfo->pszCelestialBodyName));
            };

            std::vector<std::string> oRet;
            if (noneAllowed)
                oRet.push_back("none");
            oRet.insert(oRet.end(), specialValues.begin(), specialValues.end());
            if (!currentValue.empty())
            {
                const CPLStringList aosTokens(
                    CSLTokenizeString2(currentValue.c_str(), ":", 0));
                int nCount = 0;
                std::unique_ptr<OSRCRSInfo *, decltype(&OSRDestroyCRSInfoList)>
                    pCRSList(OSRGetCRSInfoListFromDatabase(aosTokens[0],
                                                           nullptr, &nCount),
                             OSRDestroyCRSInfoList);
                std::string osCode;

                std::vector<const OSRCRSInfo *> candidates;
                for (int i = 0; i < nCount; ++i)
                {
                    const auto *entry = (pCRSList.get())[i];
                    if (!entry->bDeprecated && IsCRSCompatible(entry))
                    {
                        if (aosTokens.size() == 1 ||
                            STARTS_WITH(entry->pszCode, aosTokens[1]))
                        {
                            if (candidates.empty())
                                osCode = entry->pszCode;
                            candidates.push_back(entry);
                        }
                    }
                }
                if (candidates.size() == 1)
                {
                    oRet.push_back(std::move(osCode));
                }
                else
                {
                    if (sDatasetLongLatEnv.IsInit())
                    {
                        std::sort(
                            candidates.begin(), candidates.end(),
                            [](const OSRCRSInfo *a, const OSRCRSInfo *b)
                            {
                                const double dfXa =
                                    a->dfWestLongitudeDeg >
                                            a->dfEastLongitudeDeg
                                        ? a->dfWestLongitudeDeg -
                                              a->dfEastLongitudeDeg
                                        : (180 - a->dfWestLongitudeDeg) +
                                              (a->dfEastLongitudeDeg - -180);
                                const double dfYa = a->dfNorthLatitudeDeg -
                                                    a->dfSouthLatitudeDeg;
                                const double dfXb =
                                    b->dfWestLongitudeDeg >
                                            b->dfEastLongitudeDeg
                                        ? b->dfWestLongitudeDeg -
                                              b->dfEastLongitudeDeg
                                        : (180 - b->dfWestLongitudeDeg) +
                                              (b->dfEastLongitudeDeg - -180);
                                const double dfYb = b->dfNorthLatitudeDeg -
                                                    b->dfSouthLatitudeDeg;
                                const double diffArea =
                                    dfXa * dfYa - dfXb * dfYb;
                                if (diffArea < 0)
                                    return true;
                                if (diffArea == 0)
                                {
                                    if (std::string_view(a->pszName) ==
                                        b->pszName)
                                    {
                                        if (a->eType ==
                                                OSR_CRS_TYPE_GEOGRAPHIC_2D &&
                                            b->eType !=
                                                OSR_CRS_TYPE_GEOGRAPHIC_2D)
                                            return true;
                                        if (a->eType ==
                                                OSR_CRS_TYPE_GEOGRAPHIC_3D &&
                                            b->eType == OSR_CRS_TYPE_GEOCENTRIC)
                                            return true;
                                        return false;
                                    }
                                    return std::string_view(a->pszCode) <
                                           b->pszCode;
                                }
                                return false;
                            });
                    }

                    for (const auto *entry : candidates)
                    {
                        std::string val = std::string(entry->pszCode)
                                              .append(" -- ")
                                              .append(entry->pszName);
                        if (entry->eType == OSR_CRS_TYPE_GEOGRAPHIC_2D)
                            val.append(" (geographic 2D)");
                        else if (entry->eType == OSR_CRS_TYPE_GEOGRAPHIC_3D)
                            val.append(" (geographic 3D)");
                        else if (entry->eType == OSR_CRS_TYPE_GEOCENTRIC)
                            val.append(" (geocentric)");
                        oRet.push_back(std::move(val));
                    }
                }
            }
            if (currentValue.empty() || oRet.empty())
            {
                const CPLStringList aosAuthorities(
                    OSRGetAuthorityListFromDatabase());
                for (const char *pszAuth : cpl::Iterate(aosAuthorities))
                {
                    int nCount = 0;
                    OSRDestroyCRSInfoList(OSRGetCRSInfoListFromDatabase(
                        pszAuth, nullptr, &nCount));
                    if (nCount)
                        oRet.push_back(std::string(pszAuth).append(":"));
                }
            }
            return oRet;
        });

    return *this;
}

/************************************************************************/
/*                    GDALAlgorithm::GDALAlgorithm()                    */
/************************************************************************/

GDALAlgorithm::GDALAlgorithm(const std::string &name,
                             const std::string &description,
                             const std::string &helpURL)
    : m_name(name), m_description(description), m_helpURL(helpURL),
      m_helpFullURL(!m_helpURL.empty() && m_helpURL[0] == '/'
                        ? "https://gdal.org" + m_helpURL
                        : m_helpURL)
{
    auto &helpArg =
        AddArg("help", 'h', _("Display help message and exit"),
               &m_helpRequested)
            .SetHiddenForAPI()
            .SetCategory(GAAC_COMMON)
            .AddAction([this]()
                       { m_specialActionRequested = m_calledFromCommandLine; });
    auto &helpDocArg =
        AddArg("help-doc", 0,
               _("Display help message for use by documentation"),
               &m_helpDocRequested)
            .SetHidden()
            .AddAction([this]()
                       { m_specialActionRequested = m_calledFromCommandLine; });
    auto &jsonUsageArg =
        AddArg("json-usage", 0, _("Display usage as JSON document and exit"),
               &m_JSONUsageRequested)
            .SetHiddenForAPI()
            .SetCategory(GAAC_COMMON)
            .AddAction([this]()
                       { m_specialActionRequested = m_calledFromCommandLine; });
    AddArg("config", 0, _("Configuration option"), &m_dummyConfigOptions)
        .SetMetaVar("<KEY>=<VALUE>")
        .SetHiddenForAPI()
        .SetCategory(GAAC_COMMON)
        .AddAction(
            [this]()
            {
                ReportError(
                    CE_Warning, CPLE_AppDefined,
                    "Configuration options passed with the 'config' argument "
                    "are ignored");
            });

    AddValidationAction(
        [this, &helpArg, &helpDocArg, &jsonUsageArg]()
        {
            if (!m_calledFromCommandLine && m_specialActionRequested)
            {
                for (auto &arg : {&helpArg, &helpDocArg, &jsonUsageArg})
                {
                    if (arg->IsExplicitlySet())
                    {
                        ReportError(CE_Failure, CPLE_AppDefined,
                                    "'%s' argument only available when called "
                                    "from command line",
                                    arg->GetName().c_str());
                        return false;
                    }
                }
            }
            return true;
        });
}

/************************************************************************/
/*                   GDALAlgorithm::~GDALAlgorithm()                    */
/************************************************************************/

GDALAlgorithm::~GDALAlgorithm() = default;

/************************************************************************/
/*                    GDALAlgorithm::ParseArgument()                    */
/************************************************************************/

bool GDALAlgorithm::ParseArgument(
    GDALAlgorithmArg *arg, const std::string &name, const std::string &value,
    std::map<
        GDALAlgorithmArg *,
        std::variant<std::vector<std::string>, std::vector<int>,
                     std::vector<double>, std::vector<GDALArgDatasetValue>>>
        &inConstructionValues)
{
    const bool isListArg =
        GDALAlgorithmArgTypeIsList(arg->GetType()) && arg->GetMaxCount() > 1;
    if (arg->IsExplicitlySet() && !isListArg)
    {
        // Hack for "gdal info" to be able to pass an opened raster dataset
        // by "gdal raster info" to the "gdal vector info" algorithm.
        if (arg->SkipIfAlreadySet())
        {
            arg->SetSkipIfAlreadySet(false);
            return true;
        }

        ReportError(CE_Failure, CPLE_IllegalArg,
                    "Argument '%s' has already been specified.", name.c_str());
        return false;
    }

    if (!arg->GetRepeatedArgAllowed() &&
        cpl::contains(inConstructionValues, arg))
    {
        ReportError(CE_Failure, CPLE_IllegalArg,
                    "Argument '%s' has already been specified.", name.c_str());
        return false;
    }

    switch (arg->GetType())
    {
        case GAAT_BOOLEAN:
        {
            if (value.empty() || value == "true")
                return arg->Set(true);
            else if (value == "false")
                return arg->Set(false);
            else
            {
                ReportError(
                    CE_Failure, CPLE_IllegalArg,
                    "Invalid value '%s' for boolean argument '%s'. Should be "
                    "'true' or 'false'.",
                    value.c_str(), name.c_str());
                return false;
            }
        }

        case GAAT_STRING:
        {
            return arg->Set(value);
        }

        case GAAT_INTEGER:
        {
            errno = 0;
            char *endptr = nullptr;
            const auto val = std::strtol(value.c_str(), &endptr, 10);
            if (errno == 0 && endptr &&
                endptr == value.c_str() + value.size() && val >= INT_MIN &&
                val <= INT_MAX)
            {
                return arg->Set(static_cast<int>(val));
            }
            else
            {
                ReportError(CE_Failure, CPLE_IllegalArg,
                            "Expected integer value for argument '%s', "
                            "but got '%s'.",
                            name.c_str(), value.c_str());
                return false;
            }
        }

        case GAAT_REAL:
        {
            char *endptr = nullptr;
            double dfValue = CPLStrtod(value.c_str(), &endptr);
            if (endptr != value.c_str() + value.size())
            {
                ReportError(
                    CE_Failure, CPLE_IllegalArg,
                    "Expected real value for argument '%s', but got '%s'.",
                    name.c_str(), value.c_str());
                return false;
            }
            return arg->Set(dfValue);
        }

        case GAAT_DATASET:
        {
            return arg->SetDatasetName(value);
        }

        case GAAT_STRING_LIST:
        {
            const CPLStringList aosTokens(
                arg->GetPackedValuesAllowed()
                    ? CSLTokenizeString2(value.c_str(), ",",
                                         CSLT_HONOURSTRINGS |
                                             CSLT_PRESERVEQUOTES)
                    : CSLAddString(nullptr, value.c_str()));
            if (!cpl::contains(inConstructionValues, arg))
            {
                inConstructionValues[arg] = std::vector<std::string>();
            }
            auto &valueVector =
                std::get<std::vector<std::string>>(inConstructionValues[arg]);
            for (const char *v : aosTokens)
            {
                valueVector.push_back(v);
            }
            if (arg->GetMaxCount() == 1)
            {
                bool ret = arg->Set(std::move(valueVector));
                inConstructionValues.erase(inConstructionValues.find(arg));
                return ret;
            }

            break;
        }

        case GAAT_INTEGER_LIST:
        {
            const CPLStringList aosTokens(
                arg->GetPackedValuesAllowed()
                    ? CSLTokenizeString2(
                          value.c_str(), ",",
                          CSLT_HONOURSTRINGS | CSLT_STRIPLEADSPACES |
                              CSLT_STRIPENDSPACES | CSLT_ALLOWEMPTYTOKENS)
                    : CSLAddString(nullptr, value.c_str()));
            if (!cpl::contains(inConstructionValues, arg))
            {
                inConstructionValues[arg] = std::vector<int>();
            }
            auto &valueVector =
                std::get<std::vector<int>>(inConstructionValues[arg]);
            for (const char *v : aosTokens)
            {
                errno = 0;
                char *endptr = nullptr;
                const auto val = std::strtol(v, &endptr, 10);
                if (errno == 0 && endptr && endptr == v + strlen(v) &&
                    val >= INT_MIN && val <= INT_MAX && strlen(v) > 0)
                {
                    valueVector.push_back(static_cast<int>(val));
                }
                else
                {
                    ReportError(
                        CE_Failure, CPLE_IllegalArg,
                        "Expected list of integer value for argument '%s', "
                        "but got '%s'.",
                        name.c_str(), value.c_str());
                    return false;
                }
            }
            if (arg->GetMaxCount() == 1)
            {
                bool ret = arg->Set(std::move(valueVector));
                inConstructionValues.erase(inConstructionValues.find(arg));
                return ret;
            }

            break;
        }

        case GAAT_REAL_LIST:
        {
            const CPLStringList aosTokens(
                arg->GetPackedValuesAllowed()
                    ? CSLTokenizeString2(
                          value.c_str(), ",",
                          CSLT_HONOURSTRINGS | CSLT_STRIPLEADSPACES |
                              CSLT_STRIPENDSPACES | CSLT_ALLOWEMPTYTOKENS)
                    : CSLAddString(nullptr, value.c_str()));
            if (!cpl::contains(inConstructionValues, arg))
            {
                inConstructionValues[arg] = std::vector<double>();
            }
            auto &valueVector =
                std::get<std::vector<double>>(inConstructionValues[arg]);
            for (const char *v : aosTokens)
            {
                char *endptr = nullptr;
                double dfValue = CPLStrtod(v, &endptr);
                if (strlen(v) == 0 || endptr != v + strlen(v))
                {
                    ReportError(
                        CE_Failure, CPLE_IllegalArg,
                        "Expected list of real value for argument '%s', "
                        "but got '%s'.",
                        name.c_str(), value.c_str());
                    return false;
                }
                valueVector.push_back(dfValue);
            }
            if (arg->GetMaxCount() == 1)
            {
                bool ret = arg->Set(std::move(valueVector));
                inConstructionValues.erase(inConstructionValues.find(arg));
                return ret;
            }

            break;
        }

        case GAAT_DATASET_LIST:
        {
            if (!cpl::contains(inConstructionValues, arg))
            {
                inConstructionValues[arg] = std::vector<GDALArgDatasetValue>();
            }
            auto &valueVector = std::get<std::vector<GDALArgDatasetValue>>(
                inConstructionValues[arg]);
            if (!value.empty() && value[0] == '{' && value.back() == '}')
            {
                valueVector.push_back(GDALArgDatasetValue(value));
            }
            else
            {
                const CPLStringList aosTokens(
                    arg->GetPackedValuesAllowed()
                        ? CSLTokenizeString2(value.c_str(), ",",
                                             CSLT_HONOURSTRINGS |
                                                 CSLT_STRIPLEADSPACES)
                        : CSLAddString(nullptr, value.c_str()));
                for (const char *v : aosTokens)
                {
                    valueVector.push_back(GDALArgDatasetValue(v));
                }
            }
            if (arg->GetMaxCount() == 1)
            {
                bool ret = arg->Set(std::move(valueVector));
                inConstructionValues.erase(inConstructionValues.find(arg));
                return ret;
            }

            break;
        }
    }

    return true;
}

/************************************************************************/
/*              GDALAlgorithm::ParseCommandLineArguments()              */
/************************************************************************/

bool GDALAlgorithm::ParseCommandLineArguments(
    const std::vector<std::string> &args)
{
    if (m_parsedSubStringAlreadyCalled)
    {
        ReportError(CE_Failure, CPLE_AppDefined,
                    "ParseCommandLineArguments() can only be called once per "
                    "instance.");
        return false;
    }
    m_parsedSubStringAlreadyCalled = true;

    // AWS like syntax supported too (not advertized)
    if (args.size() == 1 && args[0] == "help")
    {
        auto arg = GetArg("help");
        assert(arg);
        arg->Set(true);
        arg->RunActions();
        return true;
    }

    if (HasSubAlgorithms())
    {
        if (args.empty())
        {
            ReportError(CE_Failure, CPLE_AppDefined, "Missing %s name.",
                        m_callPath.size() == 1 ? "command" : "subcommand");
            return false;
        }
        if (!args[0].empty() && args[0][0] == '-')
        {
            // go on argument parsing
        }
        else
        {
            const auto nCounter = CPLGetErrorCounter();
            m_selectedSubAlgHolder = InstantiateSubAlgorithm(args[0]);
            if (m_selectedSubAlgHolder)
            {
                m_selectedSubAlg = m_selectedSubAlgHolder.get();
                m_selectedSubAlg->SetReferencePathForRelativePaths(
                    m_referencePath);
                m_selectedSubAlg->m_executionForStreamOutput =
                    m_executionForStreamOutput;
                m_selectedSubAlg->m_calledFromCommandLine =
                    m_calledFromCommandLine;
                m_selectedSubAlg->m_skipValidationInParseCommandLine =
                    m_skipValidationInParseCommandLine;
                bool bRet = m_selectedSubAlg->ParseCommandLineArguments(
                    std::vector<std::string>(args.begin() + 1, args.end()));
                m_selectedSubAlg->PropagateSpecialActionTo(this);
                return bRet;
            }
            else
            {
                if (!(CPLGetErrorCounter() == nCounter + 1 &&
                      strstr(CPLGetLastErrorMsg(), "Do you mean")))
                {
                    ReportError(CE_Failure, CPLE_AppDefined,
                                "Unknown command: '%s'", args[0].c_str());
                }
                return false;
            }
        }
    }

    std::map<
        GDALAlgorithmArg *,
        std::variant<std::vector<std::string>, std::vector<int>,
                     std::vector<double>, std::vector<GDALArgDatasetValue>>>
        inConstructionValues;

    std::vector<std::string> lArgs(args);
    bool helpValueRequested = false;
    for (size_t i = 0; i < lArgs.size(); /* incremented in loop */)
    {
        const auto &strArg = lArgs[i];
        GDALAlgorithmArg *arg = nullptr;
        std::string name;
        std::string value;
        bool hasValue = false;
        if (m_calledFromCommandLine && cpl::ends_with(strArg, "=?"))
            helpValueRequested = true;
        if (strArg.size() >= 2 && strArg[0] == '-' && strArg[1] == '-')
        {
            const auto equalPos = strArg.find('=');
            name = (equalPos != std::string::npos) ? strArg.substr(0, equalPos)
                                                   : strArg;
            const std::string nameWithoutDash = name.substr(2);
            auto iterArg = m_mapLongNameToArg.find(nameWithoutDash);
            if (m_arbitraryLongNameArgsAllowed &&
                iterArg == m_mapLongNameToArg.end())
            {
                GetArg(nameWithoutDash);
                iterArg = m_mapLongNameToArg.find(nameWithoutDash);
            }
            if (iterArg == m_mapLongNameToArg.end())
            {
                const std::string bestCandidate =
                    GetSuggestionForArgumentName(nameWithoutDash);
                if (!bestCandidate.empty())
                {
                    ReportError(CE_Failure, CPLE_IllegalArg,
                                "Option '%s' is unknown. Do you mean '--%s'?",
                                name.c_str(), bestCandidate.c_str());
                }
                else
                {
                    ReportError(CE_Failure, CPLE_IllegalArg,
                                "Option '%s' is unknown.", name.c_str());
                }
                return false;
            }
            arg = iterArg->second;
            if (equalPos != std::string::npos)
            {
                hasValue = true;
                value = strArg.substr(equalPos + 1);
            }
        }
        else if (strArg.size() >= 2 && strArg[0] == '-' &&
                 CPLGetValueType(strArg.c_str()) == CPL_VALUE_STRING)
        {
            for (size_t j = 1; j < strArg.size(); ++j)
            {
                name.clear();
                name += strArg[j];
                const auto iterArg = m_mapShortNameToArg.find(name);
                if (iterArg == m_mapShortNameToArg.end())
                {
                    const std::string nameWithoutDash = strArg.substr(1);
                    if (m_mapLongNameToArg.find(nameWithoutDash) !=
                        m_mapLongNameToArg.end())
                    {
                        ReportError(CE_Failure, CPLE_IllegalArg,
                                    "Short name option '%s' is unknown. Do you "
                                    "mean '--%s' (with leading double dash) ?",
                                    name.c_str(), nameWithoutDash.c_str());
                    }
                    else
                    {
                        const std::string bestCandidate =
                            GetSuggestionForArgumentName(nameWithoutDash);
                        if (!bestCandidate.empty())
                        {
                            ReportError(
                                CE_Failure, CPLE_IllegalArg,
                                "Short name option '%s' is unknown. Do you "
                                "mean '--%s' (with leading double dash) ?",
                                name.c_str(), bestCandidate.c_str());
                        }
                        else
                        {
                            ReportError(CE_Failure, CPLE_IllegalArg,
                                        "Short name option '%s' is unknown.",
                                        name.c_str());
                        }
                    }
                    return false;
                }
                arg = iterArg->second;
                if (strArg.size() > 2)
                {
                    if (arg->GetType() != GAAT_BOOLEAN)
                    {
                        ReportError(CE_Failure, CPLE_IllegalArg,
                                    "Invalid argument '%s'. Option '%s' is not "
                                    "a boolean option.",
                                    strArg.c_str(), name.c_str());
                        return false;
                    }

                    if (!ParseArgument(arg, name, "true", inConstructionValues))
                        return false;
                }
            }
            if (strArg.size() > 2)
            {
                lArgs.erase(lArgs.begin() + i);
                continue;
            }
        }
        else
        {
            ++i;
            continue;
        }
        CPLAssert(arg);

        if (arg && arg->GetType() == GAAT_BOOLEAN)
        {
            if (!hasValue)
            {
                hasValue = true;
                value = "true";
            }
        }

        if (!hasValue)
        {
            if (i + 1 == lArgs.size())
            {
                if (m_parseForAutoCompletion)
                {
                    lArgs.erase(lArgs.begin() + i);
                    break;
                }
                ReportError(
                    CE_Failure, CPLE_IllegalArg,
                    "Expected value for argument '%s', but ran short of tokens",
                    name.c_str());
                return false;
            }
            value = lArgs[i + 1];
            lArgs.erase(lArgs.begin() + i + 1);
        }

        if (arg && !ParseArgument(arg, name, value, inConstructionValues))
            return false;

        lArgs.erase(lArgs.begin() + i);
    }

    if (m_specialActionRequested)
    {
        return true;
    }

    const auto ProcessInConstructionValues = [&inConstructionValues]()
    {
        for (auto &[arg, value] : inConstructionValues)
        {
            if (arg->GetType() == GAAT_STRING_LIST)
            {
                if (!arg->Set(std::get<std::vector<std::string>>(
                        inConstructionValues[arg])))
                {
                    return false;
                }
            }
            else if (arg->GetType() == GAAT_INTEGER_LIST)
            {
                if (!arg->Set(
                        std::get<std::vector<int>>(inConstructionValues[arg])))
                {
                    return false;
                }
            }
            else if (arg->GetType() == GAAT_REAL_LIST)
            {
                if (!arg->Set(std::get<std::vector<double>>(
                        inConstructionValues[arg])))
                {
                    return false;
                }
            }
            else if (arg->GetType() == GAAT_DATASET_LIST)
            {
                if (!arg->Set(
                        std::move(std::get<std::vector<GDALArgDatasetValue>>(
                            inConstructionValues[arg]))))
                {
                    return false;
                }
            }
        }
        return true;
    };

    // Process positional arguments that have not been set through their
    // option name.
    size_t i = 0;
    size_t iCurPosArg = 0;

    // Special case for <INPUT> <AUXILIARY>... <OUTPUT>
    if (m_positionalArgs.size() == 3 &&
        (m_positionalArgs[0]->IsRequired() ||
         m_positionalArgs[0]->GetMinCount() == 1) &&
        m_positionalArgs[0]->GetMaxCount() == 1 &&
        (m_positionalArgs[1]->IsRequired() ||
         m_positionalArgs[1]->GetMinCount() == 1) &&
        /* Second argument may have several occurrences */
        m_positionalArgs[1]->GetMaxCount() >= 1 &&
        (m_positionalArgs[2]->IsRequired() ||
         m_positionalArgs[2]->GetMinCount() == 1) &&
        m_positionalArgs[2]->GetMaxCount() == 1 &&
        !m_positionalArgs[0]->IsExplicitlySet() &&
        !m_positionalArgs[1]->IsExplicitlySet() &&
        !m_positionalArgs[2]->IsExplicitlySet())
    {
        if (lArgs.size() - i < 3)
        {
            ReportError(CE_Failure, CPLE_AppDefined,
                        "Not enough positional values.");
            return false;
        }
        bool ok = ParseArgument(m_positionalArgs[0],
                                m_positionalArgs[0]->GetName().c_str(),
                                lArgs[i], inConstructionValues);
        if (ok)
        {
            ++i;
            for (; i + 1 < lArgs.size() && ok; ++i)
            {
                ok = ParseArgument(m_positionalArgs[1],
                                   m_positionalArgs[1]->GetName().c_str(),
                                   lArgs[i], inConstructionValues);
            }
        }
        if (ok)
        {
            ok = ParseArgument(m_positionalArgs[2],
                               m_positionalArgs[2]->GetName().c_str(), lArgs[i],
                               inConstructionValues);
            ++i;
        }
        if (!ok)
        {
            ProcessInConstructionValues();
            return false;
        }
    }

    if (m_inputDatasetCanBeOmitted && m_positionalArgs.size() >= 1 &&
        !m_positionalArgs[0]->IsExplicitlySet() &&
        m_positionalArgs[0]->GetName() == GDAL_ARG_NAME_INPUT &&
        (m_positionalArgs[0]->GetType() == GAAT_DATASET ||
         m_positionalArgs[0]->GetType() == GAAT_DATASET_LIST))
    {
        ++iCurPosArg;
    }

    while (i < lArgs.size() && iCurPosArg < m_positionalArgs.size())
    {
        GDALAlgorithmArg *arg = m_positionalArgs[iCurPosArg];
        while (arg->IsExplicitlySet())
        {
            ++iCurPosArg;
            if (iCurPosArg == m_positionalArgs.size())
                break;
            arg = m_positionalArgs[iCurPosArg];
        }
        if (iCurPosArg == m_positionalArgs.size())
        {
            break;
        }
        if (GDALAlgorithmArgTypeIsList(arg->GetType()) &&
            arg->GetMinCount() != arg->GetMaxCount())
        {
            if (iCurPosArg == 0)
            {
                size_t nCountAtEnd = 0;
                for (size_t j = 1; j < m_positionalArgs.size(); j++)
                {
                    const auto *otherArg = m_positionalArgs[j];
                    if (GDALAlgorithmArgTypeIsList(otherArg->GetType()))
                    {
                        if (otherArg->GetMinCount() != otherArg->GetMaxCount())
                        {
                            ReportError(
                                CE_Failure, CPLE_AppDefined,
                                "Ambiguity in definition of positional "
                                "argument "
                                "'%s' given it has a varying number of values, "
                                "but follows argument '%s' which also has a "
                                "varying number of values",
                                otherArg->GetName().c_str(),
                                arg->GetName().c_str());
                            ProcessInConstructionValues();
                            return false;
                        }
                        nCountAtEnd += otherArg->GetMinCount();
                    }
                    else
                    {
                        if (!otherArg->IsRequired())
                        {
                            ReportError(
                                CE_Failure, CPLE_AppDefined,
                                "Ambiguity in definition of positional "
                                "argument "
                                "'%s', given it is not required but follows "
                                "argument '%s' which has a varying number of "
                                "values",
                                otherArg->GetName().c_str(),
                                arg->GetName().c_str());
                            ProcessInConstructionValues();
                            return false;
                        }
                        nCountAtEnd++;
                    }
                }
                if (lArgs.size() < nCountAtEnd)
                {
                    ReportError(CE_Failure, CPLE_AppDefined,
                                "Not enough positional values.");
                    ProcessInConstructionValues();
                    return false;
                }
                for (; i < lArgs.size() - nCountAtEnd; ++i)
                {
                    if (!ParseArgument(arg, arg->GetName().c_str(), lArgs[i],
                                       inConstructionValues))
                    {
                        ProcessInConstructionValues();
                        return false;
                    }
                }
            }
            else if (iCurPosArg == m_positionalArgs.size() - 1)
            {
                for (; i < lArgs.size(); ++i)
                {
                    if (!ParseArgument(arg, arg->GetName().c_str(), lArgs[i],
                                       inConstructionValues))
                    {
                        ProcessInConstructionValues();
                        return false;
                    }
                }
            }
            else
            {
                ReportError(CE_Failure, CPLE_AppDefined,
                            "Ambiguity in definition of positional arguments: "
                            "arguments with varying number of values must be "
                            "first or last one.");
                return false;
            }
        }
        else
        {
            if (lArgs.size() - i < static_cast<size_t>(arg->GetMaxCount()))
            {
                ReportError(CE_Failure, CPLE_AppDefined,
                            "Not enough positional values.");
                return false;
            }
            const size_t iMax = i + arg->GetMaxCount();
            for (; i < iMax; ++i)
            {
                if (!ParseArgument(arg, arg->GetName().c_str(), lArgs[i],
                                   inConstructionValues))
                {
                    ProcessInConstructionValues();
                    return false;
                }
            }
        }
        ++iCurPosArg;
    }

    if (i < lArgs.size())
    {
        ReportError(CE_Failure, CPLE_AppDefined,
                    "Positional values starting at '%s' are not expected.",
                    lArgs[i].c_str());
        return false;
    }

    if (!ProcessInConstructionValues())
    {
        return false;
    }

    // Skip to first unset positional argument.
    while (iCurPosArg < m_positionalArgs.size() &&
           m_positionalArgs[iCurPosArg]->IsExplicitlySet())
    {
        ++iCurPosArg;
    }
    // Check if this positional argument is required.
    if (iCurPosArg < m_positionalArgs.size() && !helpValueRequested &&
        (GDALAlgorithmArgTypeIsList(m_positionalArgs[iCurPosArg]->GetType())
             ? m_positionalArgs[iCurPosArg]->GetMinCount() > 0
             : m_positionalArgs[iCurPosArg]->IsRequired()))
    {
        ReportError(CE_Failure, CPLE_AppDefined,
                    "Positional arguments starting at '%s' have not been "
                    "specified.",
                    m_positionalArgs[iCurPosArg]->GetMetaVar().c_str());
        return false;
    }

    if (m_calledFromCommandLine)
    {
        for (auto &arg : m_args)
        {
            if (arg->IsExplicitlySet() &&
                ((arg->GetType() == GAAT_STRING &&
                  arg->Get<std::string>() == "?") ||
                 (arg->GetType() == GAAT_STRING_LIST &&
                  arg->Get<std::vector<std::string>>().size() == 1 &&
                  arg->Get<std::vector<std::string>>()[0] == "?")))
            {
                {
                    CPLErrorStateBackuper oBackuper(CPLQuietErrorHandler);
                    ValidateArguments();
                }

                auto choices = arg->GetChoices();
                if (choices.empty())
                    choices = arg->GetAutoCompleteChoices(std::string());
                if (!choices.empty())
                {
                    if (choices.size() == 1)
                    {
                        ReportError(
                            CE_Failure, CPLE_AppDefined,
                            "Single potential value for argument '%s' is '%s'",
                            arg->GetName().c_str(), choices.front().c_str());
                    }
                    else
                    {
                        std::string msg("Potential values for argument '");
                        msg += arg->GetName();
                        msg += "' are:";
                        for (const auto &v : choices)
                        {
                            msg += "\n- ";
                            msg += v;
                        }
                        ReportError(CE_Failure, CPLE_AppDefined, "%s",
                                    msg.c_str());
                    }
                    return false;
                }
            }
        }
    }

    return m_skipValidationInParseCommandLine || ValidateArguments();
}

/************************************************************************/
/*                     GDALAlgorithm::ReportError()                     */
/************************************************************************/

//! @cond Doxygen_Suppress
void GDALAlgorithm::ReportError(CPLErr eErrClass, CPLErrorNum err_no,
                                const char *fmt, ...) const
{
    va_list args;
    va_start(args, fmt);
    CPLError(eErrClass, err_no, "%s",
             std::string(m_name)
                 .append(": ")
                 .append(CPLString().vPrintf(fmt, args))
                 .c_str());
    va_end(args);
}

//! @endcond

/************************************************************************/
/*                  GDALAlgorithm::ProcessDatasetArg()                  */
/************************************************************************/

bool GDALAlgorithm::ProcessDatasetArg(GDALAlgorithmArg *arg,
                                      GDALAlgorithm *algForOutput)
{
    bool ret = true;

    const auto updateArg = algForOutput->GetArg(GDAL_ARG_NAME_UPDATE);
    const bool hasUpdateArg = updateArg && updateArg->GetType() == GAAT_BOOLEAN;
    const bool update = hasUpdateArg && updateArg->Get<bool>();

    const auto appendArg = algForOutput->GetArg(GDAL_ARG_NAME_APPEND);
    const bool hasAppendArg = appendArg && appendArg->GetType() == GAAT_BOOLEAN;
    const bool append = hasAppendArg && appendArg->Get<bool>();

    const auto overwriteArg = algForOutput->GetArg(GDAL_ARG_NAME_OVERWRITE);
    const bool overwrite =
        (arg->IsOutput() && overwriteArg &&
         overwriteArg->GetType() == GAAT_BOOLEAN && overwriteArg->Get<bool>());

    auto outputArg = algForOutput->GetArg(GDAL_ARG_NAME_OUTPUT);
    auto &val = [arg]() -> GDALArgDatasetValue &
    {
        if (arg->GetType() == GAAT_DATASET_LIST)
            return arg->Get<std::vector<GDALArgDatasetValue>>()[0];
        else
            return arg->Get<GDALArgDatasetValue>();
    }();
    const bool onlyInputSpecifiedInUpdateAndOutputNotRequired =
        arg->GetName() == GDAL_ARG_NAME_INPUT && outputArg &&
        !outputArg->IsExplicitlySet() && !outputArg->IsRequired() && update &&
        !overwrite;

    // Used for nested pipelines
    const auto oIterDatasetNameToDataset =
        val.IsNameSet() ? m_oMapDatasetNameToDataset.find(val.GetName())
                        : m_oMapDatasetNameToDataset.end();

    if (!val.GetDatasetRef() && !val.IsNameSet())
    {
        ReportError(CE_Failure, CPLE_AppDefined,
                    "Argument '%s' has no dataset object or dataset name.",
                    arg->GetName().c_str());
        ret = false;
    }
    else if (val.GetDatasetRef() && !CheckCanSetDatasetObject(arg))
    {
        return false;
    }
    else if (m_inputDatasetCanBeOmitted &&
             val.GetName() == GDAL_DATASET_PIPELINE_PLACEHOLDER_VALUE &&
             !arg->IsOutput())
    {
        return true;
    }
    else if (!val.GetDatasetRef() &&
             (arg->AutoOpenDataset() ||
              oIterDatasetNameToDataset != m_oMapDatasetNameToDataset.end()) &&
             (!arg->IsOutput() || (arg == outputArg && update && !overwrite) ||
              onlyInputSpecifiedInUpdateAndOutputNotRequired))
    {
        int flags = arg->GetDatasetType();
        bool assignToOutputArg = false;

        // Check if input and output parameters point to the same
        // filename (for vector datasets)
        if (arg->GetName() == GDAL_ARG_NAME_INPUT && update && !overwrite &&
            outputArg && outputArg->GetType() == GAAT_DATASET)
        {
            auto &outputVal = outputArg->Get<GDALArgDatasetValue>();
            if (!outputVal.GetDatasetRef() &&
                outputVal.GetName() == val.GetName() &&
                (outputArg->GetDatasetInputFlags() & GADV_OBJECT) != 0)
            {
                assignToOutputArg = true;
                flags |= GDAL_OF_UPDATE | GDAL_OF_VERBOSE_ERROR;
            }
            else if (onlyInputSpecifiedInUpdateAndOutputNotRequired)
            {
                flags |= GDAL_OF_UPDATE | GDAL_OF_VERBOSE_ERROR;
            }
        }

        if (!arg->IsOutput() || arg->GetDatasetInputFlags() == GADV_NAME)
            flags |= GDAL_OF_VERBOSE_ERROR;
        if ((arg == outputArg || !outputArg) && update)
        {
            flags |= GDAL_OF_UPDATE;
            if (!append)
                flags |= GDAL_OF_VERBOSE_ERROR;
        }

        const auto readOnlyArg = GetArg(GDAL_ARG_NAME_READ_ONLY);
        const bool readOnly =
            (readOnlyArg && readOnlyArg->GetType() == GAAT_BOOLEAN &&
             readOnlyArg->Get<bool>());
        if (readOnly)
            flags &= ~GDAL_OF_UPDATE;

        CPLStringList aosOpenOptions;
        CPLStringList aosAllowedDrivers;
        if (arg->IsInput())
        {
            if (arg == outputArg)
            {
                if (update && !overwrite)
                {
                    const auto ooArg = GetArg(GDAL_ARG_NAME_OUTPUT_OPEN_OPTION);
                    if (ooArg && ooArg->GetType() == GAAT_STRING_LIST)
                        aosOpenOptions = CPLStringList(
                            ooArg->Get<std::vector<std::string>>());
                }
            }
            else
            {
                const auto ooArg = GetArg(GDAL_ARG_NAME_OPEN_OPTION);
                if (ooArg && ooArg->GetType() == GAAT_STRING_LIST)
                    aosOpenOptions =
                        CPLStringList(ooArg->Get<std::vector<std::string>>());

                const auto ifArg = GetArg(GDAL_ARG_NAME_INPUT_FORMAT);
                if (ifArg && ifArg->GetType() == GAAT_STRING_LIST)
                    aosAllowedDrivers =
                        CPLStringList(ifArg->Get<std::vector<std::string>>());
            }
        }

        std::string osDatasetName = val.GetName();
        if (!m_referencePath.empty())
        {
            osDatasetName = GDALDataset::BuildFilename(
                osDatasetName.c_str(), m_referencePath.c_str(), true);
        }
        if (osDatasetName == "-" && (flags & GDAL_OF_UPDATE) == 0)
            osDatasetName = "/vsistdin/";

        // Handle special case of overview delete in GTiff which would fail
        // if it is COG without IGNORE_COG_LAYOUT_BREAK=YES open option.
        if ((flags & GDAL_OF_UPDATE) != 0 && m_callPath.size() == 4 &&
            m_callPath[2] == "overview" && m_callPath[3] == "delete" &&
            aosOpenOptions.FetchNameValue("IGNORE_COG_LAYOUT_BREAK") == nullptr)
        {
            CPLErrorStateBackuper oBackuper(CPLQuietErrorHandler);
            GDALDriverH hDrv =
                GDALIdentifyDriver(osDatasetName.c_str(), nullptr);
            if (hDrv && EQUAL(GDALGetDescription(hDrv), "GTiff"))
            {
                // Cleaning does not break COG layout
                aosOpenOptions.SetNameValue("IGNORE_COG_LAYOUT_BREAK", "YES");
            }
        }

        GDALDataset *poDS;
        {
            // The PostGISRaster may emit an error message, that is not
            // relevant, if it is the vector driver that was intended
            std::unique_ptr<CPLErrorStateBackuper> poBackuper;
            if (cpl::starts_with(osDatasetName, "PG:") &&
                (flags & (GDAL_OF_RASTER | GDAL_OF_VECTOR)) != 0)
            {
                poBackuper = std::make_unique<CPLErrorStateBackuper>(
                    CPLQuietErrorHandler);
            }

            CPL_IGNORE_RET_VAL(poBackuper);
            poDS = oIterDatasetNameToDataset != m_oMapDatasetNameToDataset.end()
                       ? oIterDatasetNameToDataset->second
                       : GDALDataset::Open(osDatasetName.c_str(), flags,
                                           aosAllowedDrivers.List(),
                                           aosOpenOptions.List());

            // Retry with PostGIS vector driver
            if (!poDS && poBackuper &&
                GetGDALDriverManager()->GetDriverByName("PostGISRaster") &&
                aosAllowedDrivers.empty() && aosOpenOptions.empty())
            {
                poBackuper.reset();
                poDS = GDALDataset::Open(
                    osDatasetName.c_str(), flags & ~GDAL_OF_RASTER,
                    aosAllowedDrivers.List(), aosOpenOptions.List());
            }
        }

        if (poDS)
        {
            if (oIterDatasetNameToDataset != m_oMapDatasetNameToDataset.end())
            {
                if (arg->GetType() == GAAT_DATASET)
                    arg->Get<GDALArgDatasetValue>().Set(poDS->GetDescription());
                poDS->Reference();
                m_oMapDatasetNameToDataset.erase(oIterDatasetNameToDataset);
            }

            // A bit of a hack for situations like 'gdal raster clip --like "PG:..."'
            // where the PG: dataset will be first opened with the PostGISRaster
            // driver whereas the PostgreSQL (vector) one is actually wanted.
            if (poDS->GetRasterCount() == 0 && (flags & GDAL_OF_RASTER) != 0 &&
                (flags & GDAL_OF_VECTOR) != 0 && aosAllowedDrivers.empty() &&
                aosOpenOptions.empty())
            {
                auto poDrv = poDS->GetDriver();
                if (poDrv && EQUAL(poDrv->GetDescription(), "PostGISRaster"))
                {
                    // Retry with PostgreSQL (vector) driver
                    std::unique_ptr<GDALDataset> poTmpDS(GDALDataset::Open(
                        osDatasetName.c_str(), flags & ~GDAL_OF_RASTER));
                    if (poTmpDS)
                    {
                        poDS->ReleaseRef();
                        poDS = poTmpDS.release();
                    }
                }
            }

            if (assignToOutputArg)
            {
                // Avoid opening twice the same datasource if it is both
                // the input and output.
                // Known to cause problems with at least FGdb, SQLite
                // and GPKG drivers. See #4270
                // Restrict to those 3 drivers. For example it is known
                // to break with the PG driver due to the way it
                // manages transactions.
                auto poDriver = poDS->GetDriver();
                if (poDriver && (EQUAL(poDriver->GetDescription(), "FileGDB") ||
                                 EQUAL(poDriver->GetDescription(), "SQLite") ||
                                 EQUAL(poDriver->GetDescription(), "GPKG")))
                {
                    outputArg->Get<GDALArgDatasetValue>().Set(poDS);
                }
            }
            val.SetDatasetOpenedByAlgorithm();
            val.Set(poDS);
            poDS->ReleaseRef();
        }
        else if (!append)
        {
            ret = false;
        }
    }

    // Deal with overwriting the output dataset
    if (ret && arg == outputArg && val.GetDatasetRef() == nullptr)
    {
        if (!append)
        {
            // If outputting to MEM, do not try to erase a real file of the same name!
            const auto outputFormatArg =
                algForOutput->GetArg(GDAL_ARG_NAME_OUTPUT_FORMAT);
            if (!(outputFormatArg &&
                  outputFormatArg->GetType() == GAAT_STRING &&
                  (EQUAL(outputFormatArg->Get<std::string>().c_str(), "MEM") ||
                   EQUAL(outputFormatArg->Get<std::string>().c_str(),
                         "stream") ||
                   EQUAL(outputFormatArg->Get<std::string>().c_str(),
                         "Memory"))))
            {
                const char *pszType = "";
                GDALDriver *poDriver = nullptr;
                if (!val.GetName().empty() &&
                    GDALDoesFileOrDatasetExist(val.GetName().c_str(), &pszType,
                                               &poDriver))
                {
                    if (!overwrite)
                    {
                        std::string options;
                        if (algForOutput->GetArg(GDAL_ARG_NAME_OVERWRITE_LAYER))
                        {
                            options += "--";
                            options += GDAL_ARG_NAME_OVERWRITE_LAYER;
                        }
                        if (hasAppendArg)
                        {
                            if (!options.empty())
                                options += '/';
                            options += "--";
                            options += GDAL_ARG_NAME_APPEND;
                        }
                        if (hasUpdateArg)
                        {
                            if (!options.empty())
                                options += '/';
                            options += "--";
                            options += GDAL_ARG_NAME_UPDATE;
                        }

                        if (poDriver)
                        {
                            const char *pszPrefix = poDriver->GetMetadataItem(
                                GDAL_DMD_CONNECTION_PREFIX);
                            if (pszPrefix &&
                                STARTS_WITH_CI(val.GetName().c_str(),
                                               pszPrefix))
                            {
                                bool bExists = false;
                                {
                                    CPLErrorStateBackuper oBackuper(
                                        CPLQuietErrorHandler);
                                    bExists = std::unique_ptr<GDALDataset>(
                                                  GDALDataset::Open(
                                                      val.GetName().c_str())) !=
                                              nullptr;
                                }
                                if (bExists)
                                {
                                    if (!options.empty())
                                        options = " You may specify the " +
                                                  options + " option.";
                                    ReportError(CE_Failure, CPLE_AppDefined,
                                                "%s '%s' already exists.%s",
                                                pszType, val.GetName().c_str(),
                                                options.c_str());
                                    return false;
                                }

                                return true;
                            }
                        }

                        if (!options.empty())
                            options = '/' + options;
                        ReportError(
                            CE_Failure, CPLE_AppDefined,
                            "%s '%s' already exists. You may specify the "
                            "--overwrite%s option.",
                            pszType, val.GetName().c_str(), options.c_str());
                        return false;
                    }
                    else if (EQUAL(pszType, "File"))
                    {
                        if (VSIUnlink(val.GetName().c_str()) != 0)
                        {
                            ReportError(CE_Failure, CPLE_AppDefined,
                                        "Deleting %s failed: %s",
                                        val.GetName().c_str(),
                                        VSIStrerror(errno));
                            return false;
                        }
                    }
                    else if (EQUAL(pszType, "Directory"))
                    {
                        // We don't want the user to accidentally erase a non-GDAL dataset
                        ReportError(CE_Failure, CPLE_AppDefined,
                                    "Directory '%s' already exists, but is not "
                                    "recognized as a valid GDAL dataset. "
                                    "Please manually delete it before retrying",
                                    val.GetName().c_str());
                        return false;
                    }
                    else if (poDriver)
                    {
                        bool bDeleteOK;
                        {
                            CPLErrorStateBackuper oBackuper(
                                CPLQuietErrorHandler);
                            bDeleteOK = (poDriver->Delete(
                                             val.GetName().c_str()) == CE_None);
                        }
                        VSIStatBufL sStat;
                        if (!bDeleteOK &&
                            VSIStatL(val.GetName().c_str(), &sStat) == 0)
                        {
                            if (VSI_ISDIR(sStat.st_mode))
                            {
                                // We don't want the user to accidentally erase a non-GDAL dataset
                                ReportError(
                                    CE_Failure, CPLE_AppDefined,
                                    "Directory '%s' already exists, but is not "
                                    "recognized as a valid GDAL dataset. "
                                    "Please manually delete it before retrying",
                                    val.GetName().c_str());
                                return false;
                            }
                            else if (VSIUnlink(val.GetName().c_str()) != 0)
                            {
                                ReportError(CE_Failure, CPLE_AppDefined,
                                            "Deleting %s failed: %s",
                                            val.GetName().c_str(),
                                            VSIStrerror(errno));
                                return false;
                            }
                        }
                    }
                }
            }
        }
    }

    // If outputting to stdout, automatically turn off progress bar
    if (arg == outputArg && val.GetName() == "/vsistdout/")
    {
        auto quietArg = GetArg(GDAL_ARG_NAME_QUIET);
        if (quietArg && quietArg->GetType() == GAAT_BOOLEAN)
            quietArg->Set(true);
    }

    return ret;
}

/************************************************************************/
/*                  GDALAlgorithm::ValidateArguments()                  */
/************************************************************************/

bool GDALAlgorithm::ValidateArguments()
{
    if (m_selectedSubAlg)
        return m_selectedSubAlg->ValidateArguments();

    if (m_specialActionRequested)
        return true;

    m_arbitraryLongNameArgsAllowed = false;

    // If only --output=format=MEM/stream is specified and not --output,
    // then set empty name for --output.
    auto outputArg = GetArg(GDAL_ARG_NAME_OUTPUT);
    auto outputFormatArg = GetArg(GDAL_ARG_NAME_OUTPUT_FORMAT);
    if (outputArg && outputFormatArg && outputFormatArg->IsExplicitlySet() &&
        !outputArg->IsExplicitlySet() &&
        outputFormatArg->GetType() == GAAT_STRING &&
        (EQUAL(outputFormatArg->Get<std::string>().c_str(), "MEM") ||
         EQUAL(outputFormatArg->Get<std::string>().c_str(), "stream")) &&
        outputArg->GetType() == GAAT_DATASET &&
        (outputArg->GetDatasetInputFlags() & GADV_NAME))
    {
        outputArg->Get<GDALArgDatasetValue>().Set("");
    }

    // The method may emit several errors if several constraints are not met.
    bool ret = true;
    std::map<std::string, std::string> mutualExclusionGroupUsed;
    std::map<std::string, std::vector<std::string>> mutualDependencyGroupUsed;
    for (auto &arg : m_args)
    {
        // Check mutually exclusive/dependent arguments
        if (arg->IsExplicitlySet())
        {

            const auto &mutualExclusionGroup = arg->GetMutualExclusionGroup();
            if (!mutualExclusionGroup.empty())
            {
                auto oIter =
                    mutualExclusionGroupUsed.find(mutualExclusionGroup);
                if (oIter != mutualExclusionGroupUsed.end())
                {
                    ret = false;
                    ReportError(
                        CE_Failure, CPLE_AppDefined,
                        "Argument '%s' is mutually exclusive with '%s'.",
                        arg->GetName().c_str(), oIter->second.c_str());
                }
                else
                {
                    mutualExclusionGroupUsed[mutualExclusionGroup] =
                        arg->GetName();
                }
            }

            const auto &mutualDependencyGroup = arg->GetMutualDependencyGroup();
            if (!mutualDependencyGroup.empty())
            {
                if (mutualDependencyGroupUsed.find(mutualDependencyGroup) ==
                    mutualDependencyGroupUsed.end())
                {
                    mutualDependencyGroupUsed[mutualDependencyGroup] = {
                        arg->GetName()};
                }
                else
                {
                    mutualDependencyGroupUsed[mutualDependencyGroup].push_back(
                        arg->GetName());
                }
            }

            // Check direct dependencies
            for (const auto &dependency : arg->GetDirectDependencies())
            {
                auto depArg = GetArg(dependency);
                if (!depArg)
                {
                    ret = false;
                    ReportError(CE_Failure, CPLE_AppDefined,
                                "Argument '%s' depends on argument '%s' that "
                                "is not defined.",
                                arg->GetName().c_str(), dependency.c_str());
                }
                else if (!depArg->IsExplicitlySet())
                {
                    ret = false;
                    ReportError(CE_Failure, CPLE_AppDefined,
                                "Argument '%s' depends on argument '%s' that "
                                "has not been specified.",
                                arg->GetName().c_str(),
                                depArg->GetName().c_str());
                }
            }
        }

        if (arg->IsRequired() && !arg->IsExplicitlySet() &&
            !arg->HasDefaultValue())
        {
            bool emitError = true;
            const auto &mutualExclusionGroup = arg->GetMutualExclusionGroup();
            if (!mutualExclusionGroup.empty())
            {
                for (const auto &otherArg : m_args)
                {
                    if (otherArg->GetMutualExclusionGroup() ==
                            mutualExclusionGroup &&
                        otherArg->IsExplicitlySet())
                    {
                        emitError = false;
                        break;
                    }
                }
            }
            if (emitError && !(m_inputDatasetCanBeOmitted &&
                               arg->GetName() == GDAL_ARG_NAME_INPUT &&
                               (arg->GetType() == GAAT_DATASET ||
                                arg->GetType() == GAAT_DATASET_LIST)))
            {
                ReportError(CE_Failure, CPLE_AppDefined,
                            "Required argument '%s' has not been specified.",
                            arg->GetName().c_str());
                ret = false;
            }
        }
        else if (arg->IsExplicitlySet() && arg->GetType() == GAAT_DATASET)
        {
            if (!ProcessDatasetArg(arg.get(), this))
                ret = false;
        }

        if (arg->IsExplicitlySet() && arg->GetType() == GAAT_DATASET_LIST)
        {
            auto &listVal = arg->Get<std::vector<GDALArgDatasetValue>>();
            if (listVal.size() == 1)
            {
                if (!ProcessDatasetArg(arg.get(), this))
                    ret = false;
            }
            else
            {
                for (auto &val : listVal)
                {
                    if (val.GetDatasetRef())
                    {
                        if (!CheckCanSetDatasetObject(arg.get()))
                        {
                            ret = false;
                        }
                        continue;
                    }

                    if (val.GetName().empty())
                    {
                        ReportError(CE_Failure, CPLE_AppDefined,
                                    "Argument '%s' has no dataset object or "
                                    "dataset name.",
                                    arg->GetName().c_str());
                        ret = false;
                        continue;
                    }

                    auto oIter = m_oMapDatasetNameToDataset.find(val.GetName());
                    if (oIter != m_oMapDatasetNameToDataset.end())
                    {
                        auto poDS = oIter->second;
                        val.SetDatasetOpenedByAlgorithm();
                        val.Set(poDS);
                        m_oMapDatasetNameToDataset.erase(oIter);
                        continue;
                    }

                    if (!arg->AutoOpenDataset())
                        continue;

                    int flags = arg->GetDatasetType() | GDAL_OF_VERBOSE_ERROR;

                    CPLStringList aosOpenOptions;
                    CPLStringList aosAllowedDrivers;
                    if (arg->GetName() == GDAL_ARG_NAME_INPUT)
                    {
                        const auto ooArg = GetArg(GDAL_ARG_NAME_OPEN_OPTION);
                        if (ooArg && ooArg->GetType() == GAAT_STRING_LIST)
                        {
                            aosOpenOptions = CPLStringList(
                                ooArg->Get<std::vector<std::string>>());
                        }

                        const auto ifArg = GetArg(GDAL_ARG_NAME_INPUT_FORMAT);
                        if (ifArg && ifArg->GetType() == GAAT_STRING_LIST)
                        {
                            aosAllowedDrivers = CPLStringList(
                                ifArg->Get<std::vector<std::string>>());
                        }

                        const auto updateArg = GetArg(GDAL_ARG_NAME_UPDATE);
                        if (updateArg && updateArg->GetType() == GAAT_BOOLEAN &&
                            updateArg->Get<bool>())
                        {
                            flags |= GDAL_OF_UPDATE;
                        }
                    }

                    auto poDS = std::unique_ptr<GDALDataset>(GDALDataset::Open(
                        val.GetName().c_str(), flags, aosAllowedDrivers.List(),
                        aosOpenOptions.List()));
                    if (poDS)
                    {
                        val.Set(std::move(poDS));
                    }
                    else
                    {
                        ret = false;
                    }
                }
            }
        }

        if (arg->IsExplicitlySet() && !arg->RunValidationActions())
        {
            ret = false;
        }
    }

    // Check mutual dependency groups
    std::vector<std::string> processedGroups;
    // Loop through group map and check there are not required args in the group that are not set
    for (const auto &[groupName, argNames] : mutualDependencyGroupUsed)
    {
        if (std::find(processedGroups.begin(), processedGroups.end(),
                      groupName) != processedGroups.end())
            continue;
        std::vector<std::string> missingArgs;
        for (auto &arg : m_args)
        {
            const auto &mutualDependencyGroup = arg->GetMutualDependencyGroup();
            if (mutualDependencyGroup == groupName &&
                std::find(argNames.begin(), argNames.end(), arg->GetName()) ==
                    argNames.end())
            {
                missingArgs.push_back(arg->GetName());
            }
        }
        if (!missingArgs.empty())
        {
            ret = false;
            std::string missingArgsStr;
            for (const auto &missingArg : missingArgs)
            {
                if (!missingArgsStr.empty())
                    missingArgsStr += ", ";
                missingArgsStr += missingArg;
            }
            std::string givenArgsStr;
            for (const auto &givenArg : argNames)
            {
                if (!givenArgsStr.empty())
                    givenArgsStr += ", ";
                givenArgsStr += givenArg;
            }
            ReportError(CE_Failure, CPLE_AppDefined,
                        "Argument(s) '%s' require(s) that the following "
                        "argument(s) are also specified: %s.",
                        givenArgsStr.c_str(), missingArgsStr.c_str());
        }
        processedGroups.push_back(groupName);
    }

    for (const auto &f : m_validationActions)
    {
        if (!f())
            ret = false;
    }

    return ret;
}

/************************************************************************/
/*                GDALAlgorithm::InstantiateSubAlgorithm                */
/************************************************************************/

std::unique_ptr<GDALAlgorithm>
GDALAlgorithm::InstantiateSubAlgorithm(const std::string &name,
                                       bool suggestionAllowed) const
{
    auto ret = m_subAlgRegistry.Instantiate(name);
    auto childCallPath = m_callPath;
    childCallPath.push_back(name);
    if (!ret)
    {
        ret = GDALGlobalAlgorithmRegistry::GetSingleton()
                  .InstantiateDeclaredSubAlgorithm(childCallPath);
    }
    if (ret)
    {
        ret->SetCallPath(childCallPath);
    }
    else if (suggestionAllowed)
    {
        std::string bestCandidate;
        size_t bestDistance = std::numeric_limits<size_t>::max();
        for (const std::string &candidate : GetSubAlgorithmNames())
        {
            const size_t distance =
                CPLLevenshteinDistance(name.c_str(), candidate.c_str(),
                                       /* transpositionAllowed = */ true);
            if (distance < bestDistance)
            {
                bestCandidate = candidate;
                bestDistance = distance;
            }
            else if (distance == bestDistance)
            {
                bestCandidate.clear();
            }
        }
        if (!bestCandidate.empty() && bestDistance <= 2)
        {
            CPLError(CE_Failure, CPLE_AppDefined,
                     "Algorithm '%s' is unknown. Do you mean '%s'?",
                     name.c_str(), bestCandidate.c_str());
        }
    }
    return ret;
}

/************************************************************************/
/*            GDALAlgorithm::GetSuggestionForArgumentName()             */
/************************************************************************/

std::string
GDALAlgorithm::GetSuggestionForArgumentName(const std::string &osName) const
{
    if (osName.size() >= 3)
    {
        std::string bestCandidate;
        size_t bestDistance = std::numeric_limits<size_t>::max();
        for (const auto &[key, value] : m_mapLongNameToArg)
        {
            CPL_IGNORE_RET_VAL(value);
            const size_t distance = CPLLevenshteinDistance(
                osName.c_str(), key.c_str(), /* transpositionAllowed = */ true);
            if (distance < bestDistance)
            {
                bestCandidate = key;
                bestDistance = distance;
            }
            else if (distance == bestDistance)
            {
                bestCandidate.clear();
            }
        }
        if (!bestCandidate.empty() &&
            bestDistance <= (bestCandidate.size() >= 4U ? 2U : 1U))
        {
            return bestCandidate;
        }
    }
    return std::string();
}

/************************************************************************/
/*         GDALAlgorithm::IsKnownOutputRelatedBooleanArgName()          */
/************************************************************************/

/* static */
bool GDALAlgorithm::IsKnownOutputRelatedBooleanArgName(std::string_view osName)
{
    return osName == GDAL_ARG_NAME_APPEND || osName == GDAL_ARG_NAME_UPDATE ||
           osName == GDAL_ARG_NAME_OVERWRITE ||
           osName == GDAL_ARG_NAME_OVERWRITE_LAYER;
}

/************************************************************************/
/*                   GDALAlgorithm::HasOutputString()                   */
/************************************************************************/

bool GDALAlgorithm::HasOutputString() const
{
    auto outputStringArg = GetArg(GDAL_ARG_NAME_OUTPUT_STRING);
    return outputStringArg && outputStringArg->IsOutput();
}

/************************************************************************/
/*                       GDALAlgorithm::GetArg()                        */
/************************************************************************/

GDALAlgorithmArg *GDALAlgorithm::GetArg(const std::string &osName,
                                        bool suggestionAllowed, bool isConst)
{
    const auto nPos = osName.find_first_not_of('-');
    if (nPos == std::string::npos)
        return nullptr;
    std::string osKey = osName.substr(nPos);
    {
        const auto oIter = m_mapLongNameToArg.find(osKey);
        if (oIter != m_mapLongNameToArg.end())
            return oIter->second;
    }
    {
        const auto oIter = m_mapShortNameToArg.find(osKey);
        if (oIter != m_mapShortNameToArg.end())
            return oIter->second;
    }

    if (!isConst && m_arbitraryLongNameArgsAllowed)
    {
        const auto nDotPos = osKey.find('.');
        const std::string osKeyEnd =
            nDotPos == std::string::npos ? osKey : osKey.substr(nDotPos + 1);
        if (IsKnownOutputRelatedBooleanArgName(osKeyEnd))
        {
            m_arbitraryLongNameArgsValuesBool.emplace_back(
                std::make_unique<bool>());
            AddArg(osKey, 0, std::string("User-provided argument ") + osKey,
                   m_arbitraryLongNameArgsValuesBool.back().get())
                .SetUserProvided();
        }
        else
        {
            const std::string osKeyInit = osKey;
            if (osKey == "oo")
                osKey = GDAL_ARG_NAME_OPEN_OPTION;
            else if (osKey == "co")
                osKey = GDAL_ARG_NAME_CREATION_OPTION;
            else if (osKey == "of")
                osKey = GDAL_ARG_NAME_OUTPUT_FORMAT;
            else if (osKey == "if")
                osKey = GDAL_ARG_NAME_INPUT_FORMAT;
            m_arbitraryLongNameArgsValuesStr.emplace_back(
                std::make_unique<std::string>());
            auto &arg =
                AddArg(osKey, 0, std::string("User-provided argument ") + osKey,
                       m_arbitraryLongNameArgsValuesStr.back().get())
                    .SetUserProvided();
            if (osKey != osKeyInit)
                arg.AddAlias(osKeyInit);
        }
        const auto oIter = m_mapLongNameToArg.find(osKey);
        CPLAssert(oIter != m_mapLongNameToArg.end());
        return oIter->second;
    }

    if (suggestionAllowed)
    {
        const std::string bestCandidate = GetSuggestionForArgumentName(osName);
        if (!bestCandidate.empty())
        {
            CPLError(CE_Failure, CPLE_AppDefined,
                     "Argument '%s' is unknown. Do you mean '%s'?",
                     osName.c_str(), bestCandidate.c_str());
        }
    }

    return nullptr;
}

/************************************************************************/
/*                     GDALAlgorithm::AddAliasFor()                     */
/************************************************************************/

//! @cond Doxygen_Suppress
void GDALAlgorithm::AddAliasFor(GDALInConstructionAlgorithmArg *arg,
                                const std::string &alias)
{
    if (cpl::contains(m_mapLongNameToArg, alias))
    {
        ReportError(CE_Failure, CPLE_AppDefined, "Name '%s' already declared.",
                    alias.c_str());
    }
    else
    {
        m_mapLongNameToArg[alias] = arg;
    }
}

//! @endcond

/************************************************************************/
/*                GDALAlgorithm::AddShortNameAliasFor()                 */
/************************************************************************/

//! @cond Doxygen_Suppress
void GDALAlgorithm::AddShortNameAliasFor(GDALInConstructionAlgorithmArg *arg,
                                         char shortNameAlias)
{
    std::string alias;
    alias += shortNameAlias;
    if (cpl::contains(m_mapShortNameToArg, alias))
    {
        ReportError(CE_Failure, CPLE_AppDefined,
                    "Short name '%s' already declared.", alias.c_str());
    }
    else
    {
        m_mapShortNameToArg[alias] = arg;
    }
}

//! @endcond

/************************************************************************/
/*                    GDALAlgorithm::SetPositional()                    */
/************************************************************************/

//! @cond Doxygen_Suppress
void GDALAlgorithm::SetPositional(GDALInConstructionAlgorithmArg *arg)
{
    CPLAssert(std::find(m_positionalArgs.begin(), m_positionalArgs.end(),
                        arg) == m_positionalArgs.end());
    m_positionalArgs.push_back(arg);
}

//! @endcond

/************************************************************************/
/*                  GDALAlgorithm::HasSubAlgorithms()                   */
/************************************************************************/

bool GDALAlgorithm::HasSubAlgorithms() const
{
    if (!m_subAlgRegistry.empty())
        return true;
    return !GDALGlobalAlgorithmRegistry::GetSingleton()
                .GetDeclaredSubAlgorithmNames(m_callPath)
                .empty();
}

/************************************************************************/
/*                GDALAlgorithm::GetSubAlgorithmNames()                 */
/************************************************************************/

std::vector<std::string> GDALAlgorithm::GetSubAlgorithmNames() const
{
    std::vector<std::string> ret = m_subAlgRegistry.GetNames();
    const auto other = GDALGlobalAlgorithmRegistry::GetSingleton()
                           .GetDeclaredSubAlgorithmNames(m_callPath);
    ret.insert(ret.end(), other.begin(), other.end());
    if (!other.empty())
        std::sort(ret.begin(), ret.end());
    return ret;
}

/************************************************************************/
/*                       GDALAlgorithm::AddArg()                        */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddArg(std::unique_ptr<GDALInConstructionAlgorithmArg> arg)
{
    auto argRaw = arg.get();
    const auto &longName = argRaw->GetName();
    if (!longName.empty())
    {
        if (longName[0] == '-')
        {
            ReportError(CE_Failure, CPLE_AppDefined,
                        "Long name '%s' should not start with '-'",
                        longName.c_str());
        }
        if (longName.find('=') != std::string::npos)
        {
            ReportError(CE_Failure, CPLE_AppDefined,
                        "Long name '%s' should not contain a '=' character",
                        longName.c_str());
        }
        if (cpl::contains(m_mapLongNameToArg, longName))
        {
            ReportError(CE_Failure, CPLE_AppDefined,
                        "Long name '%s' already declared", longName.c_str());
        }
        m_mapLongNameToArg[longName] = argRaw;
    }
    const auto &shortName = argRaw->GetShortName();
    if (!shortName.empty())
    {
        if (shortName.size() != 1 ||
            !((shortName[0] >= 'a' && shortName[0] <= 'z') ||
              (shortName[0] >= 'A' && shortName[0] <= 'Z') ||
              (shortName[0] >= '0' && shortName[0] <= '9')))
        {
            ReportError(CE_Failure, CPLE_AppDefined,
                        "Short name '%s' should be a single letter or digit",
                        shortName.c_str());
        }
        if (cpl::contains(m_mapShortNameToArg, shortName))
        {
            ReportError(CE_Failure, CPLE_AppDefined,
                        "Short name '%s' already declared", shortName.c_str());
        }
        m_mapShortNameToArg[shortName] = argRaw;
    }
    m_args.emplace_back(std::move(arg));
    return *(
        cpl::down_cast<GDALInConstructionAlgorithmArg *>(m_args.back().get()));
}

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddArg(const std::string &longName, char chShortName,
                      const std::string &helpMessage, bool *pValue)
{
    return AddArg(std::make_unique<GDALInConstructionAlgorithmArg>(
        this,
        GDALAlgorithmArgDecl(longName, chShortName, helpMessage, GAAT_BOOLEAN),
        pValue));
}

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddArg(const std::string &longName, char chShortName,
                      const std::string &helpMessage, std::string *pValue)
{
    return AddArg(std::make_unique<GDALInConstructionAlgorithmArg>(
        this,
        GDALAlgorithmArgDecl(longName, chShortName, helpMessage, GAAT_STRING),
        pValue));
}

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddArg(const std::string &longName, char chShortName,
                      const std::string &helpMessage, int *pValue)
{
    return AddArg(std::make_unique<GDALInConstructionAlgorithmArg>(
        this,
        GDALAlgorithmArgDecl(longName, chShortName, helpMessage, GAAT_INTEGER),
        pValue));
}

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddArg(const std::string &longName, char chShortName,
                      const std::string &helpMessage, double *pValue)
{
    return AddArg(std::make_unique<GDALInConstructionAlgorithmArg>(
        this,
        GDALAlgorithmArgDecl(longName, chShortName, helpMessage, GAAT_REAL),
        pValue));
}

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddArg(const std::string &longName, char chShortName,
                      const std::string &helpMessage,
                      GDALArgDatasetValue *pValue, GDALArgDatasetType type)
{
    auto &arg = AddArg(std::make_unique<GDALInConstructionAlgorithmArg>(
                           this,
                           GDALAlgorithmArgDecl(longName, chShortName,
                                                helpMessage, GAAT_DATASET),
                           pValue))
                    .SetDatasetType(type);
    pValue->SetOwnerArgument(&arg);
    return arg;
}

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddArg(const std::string &longName, char chShortName,
                      const std::string &helpMessage,
                      std::vector<std::string> *pValue)
{
    return AddArg(std::make_unique<GDALInConstructionAlgorithmArg>(
        this,
        GDALAlgorithmArgDecl(longName, chShortName, helpMessage,
                             GAAT_STRING_LIST),
        pValue));
}

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddArg(const std::string &longName, char chShortName,
                      const std::string &helpMessage, std::vector<int> *pValue)
{
    return AddArg(std::make_unique<GDALInConstructionAlgorithmArg>(
        this,
        GDALAlgorithmArgDecl(longName, chShortName, helpMessage,
                             GAAT_INTEGER_LIST),
        pValue));
}

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddArg(const std::string &longName, char chShortName,
                      const std::string &helpMessage,
                      std::vector<double> *pValue)
{
    return AddArg(std::make_unique<GDALInConstructionAlgorithmArg>(
        this,
        GDALAlgorithmArgDecl(longName, chShortName, helpMessage,
                             GAAT_REAL_LIST),
        pValue));
}

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddArg(const std::string &longName, char chShortName,
                      const std::string &helpMessage,
                      std::vector<GDALArgDatasetValue> *pValue,
                      GDALArgDatasetType type)
{
    return AddArg(std::make_unique<GDALInConstructionAlgorithmArg>(
                      this,
                      GDALAlgorithmArgDecl(longName, chShortName, helpMessage,
                                           GAAT_DATASET_LIST),
                      pValue))
        .SetDatasetType(type);
}

/************************************************************************/
/*                            MsgOrDefault()                            */
/************************************************************************/

inline const char *MsgOrDefault(const char *helpMessage,
                                const char *defaultMessage)
{
    return helpMessage && helpMessage[0] ? helpMessage : defaultMessage;
}

/************************************************************************/
/*         GDALAlgorithm::SetAutoCompleteFunctionForFilename()          */
/************************************************************************/

/* static */
void GDALAlgorithm::SetAutoCompleteFunctionForFilename(
    GDALInConstructionAlgorithmArg &arg, GDALArgDatasetType type)
{
    arg.SetAutoCompleteFunction(
        [&arg,
         type](const std::string &currentValue) -> std::vector<std::string>
        {
            std::vector<std::string> oRet;

            if (arg.IsHidden())
                return oRet;

            {
                CPLErrorStateBackuper oBackuper(CPLQuietErrorHandler);
                VSIStatBufL sStat;
                if (!currentValue.empty() && currentValue.back() != '/' &&
                    VSIStatL(currentValue.c_str(), &sStat) == 0)
                {
                    return oRet;
                }
            }

            auto poDM = GetGDALDriverManager();
            std::set<std::string> oExtensions;
            if (type)
            {
                for (int i = 0; i < poDM->GetDriverCount(); ++i)
                {
                    auto poDriver = poDM->GetDriver(i);
                    if (((type & GDAL_OF_RASTER) != 0 &&
                         poDriver->GetMetadataItem(GDAL_DCAP_RASTER)) ||
                        ((type & GDAL_OF_VECTOR) != 0 &&
                         poDriver->GetMetadataItem(GDAL_DCAP_VECTOR)) ||
                        ((type & GDAL_OF_MULTIDIM_RASTER) != 0 &&
                         poDriver->GetMetadataItem(GDAL_DCAP_MULTIDIM_RASTER)))
                    {
                        const char *pszExtensions =
                            poDriver->GetMetadataItem(GDAL_DMD_EXTENSIONS);
                        if (pszExtensions)
                        {
                            const CPLStringList aosExts(
                                CSLTokenizeString2(pszExtensions, " ", 0));
                            for (const char *pszExt : cpl::Iterate(aosExts))
                                oExtensions.insert(CPLString(pszExt).tolower());
                        }
                    }
                }
            }

            std::string osDir;
            const CPLStringList aosVSIPrefixes(VSIGetFileSystemsPrefixes());
            std::string osPrefix;
            if (STARTS_WITH(currentValue.c_str(), "/vsi"))
            {
                for (const char *pszPrefix : cpl::Iterate(aosVSIPrefixes))
                {
                    if (STARTS_WITH(currentValue.c_str(), pszPrefix))
                    {
                        osPrefix = pszPrefix;
                        break;
                    }
                }
                if (osPrefix.empty())
                    return aosVSIPrefixes;
                if (currentValue == osPrefix)
                    osDir = osPrefix;
            }
            if (osDir.empty())
            {
                osDir = CPLGetDirnameSafe(currentValue.c_str());
                if (!osPrefix.empty() && osDir.size() < osPrefix.size())
                    osDir = std::move(osPrefix);
            }

            auto psDir = VSIOpenDir(osDir.c_str(), 0, nullptr);
            const std::string osSep = VSIGetDirectorySeparator(osDir.c_str());
            if (currentValue.empty())
                osDir.clear();
            const std::string currentFilename =
                CPLGetFilename(currentValue.c_str());
            if (psDir)
            {
                while (const VSIDIREntry *psEntry = VSIGetNextDirEntry(psDir))
                {
                    if ((currentFilename.empty() ||
                         STARTS_WITH(psEntry->pszName,
                                     currentFilename.c_str())) &&
                        strcmp(psEntry->pszName, ".") != 0 &&
                        strcmp(psEntry->pszName, "..") != 0 &&
                        (oExtensions.empty() ||
                         !strstr(psEntry->pszName, ".aux.xml")))
                    {
                        if (oExtensions.empty() ||
                            cpl::contains(
                                oExtensions,
                                CPLString(CPLGetExtensionSafe(psEntry->pszName))
                                    .tolower()) ||
                            VSI_ISDIR(psEntry->nMode))
                        {
                            std::string osVal;
                            if (osDir.empty() || osDir == ".")
                                osVal = psEntry->pszName;
                            else
                                osVal = CPLFormFilenameSafe(
                                    osDir.c_str(), psEntry->pszName, nullptr);
                            if (VSI_ISDIR(psEntry->nMode))
                                osVal += osSep;
                            oRet.push_back(std::move(osVal));
                        }
                    }
                }
                VSICloseDir(psDir);
            }
            return oRet;
        });
}

/************************************************************************/
/*                 GDALAlgorithm::AddInputDatasetArg()                  */
/************************************************************************/

GDALInConstructionAlgorithmArg &GDALAlgorithm::AddInputDatasetArg(
    GDALArgDatasetValue *pValue, GDALArgDatasetType type,
    bool positionalAndRequired, const char *helpMessage)
{
    auto &arg = AddArg(
        GDAL_ARG_NAME_INPUT, 'i',
        MsgOrDefault(helpMessage,
                     CPLSPrintf("Input %s dataset",
                                GDALAlgorithmArgDatasetTypeName(type).c_str())),
        pValue, type);
    if (positionalAndRequired)
        arg.SetPositional().SetRequired();

    SetAutoCompleteFunctionForFilename(arg, type);

    AddValidationAction(
        [pValue]()
        {
            if (pValue->GetName() == "-")
                pValue->Set("/vsistdin/");
            return true;
        });

    return arg;
}

/************************************************************************/
/*                 GDALAlgorithm::AddInputDatasetArg()                  */
/************************************************************************/

GDALInConstructionAlgorithmArg &GDALAlgorithm::AddInputDatasetArg(
    std::vector<GDALArgDatasetValue> *pValue, GDALArgDatasetType type,
    bool positionalAndRequired, const char *helpMessage)
{
    auto &arg =
        AddArg(GDAL_ARG_NAME_INPUT, 'i',
               MsgOrDefault(
                   helpMessage,
                   CPLSPrintf("Input %s datasets",
                              GDALAlgorithmArgDatasetTypeName(type).c_str())),
               pValue, type)
            .SetPackedValuesAllowed(false);
    if (positionalAndRequired)
        arg.SetPositional().SetRequired();

    SetAutoCompleteFunctionForFilename(arg, type);

    AddValidationAction(
        [pValue]()
        {
            for (auto &val : *pValue)
            {
                if (val.GetName() == "-")
                    val.Set("/vsistdin/");
            }
            return true;
        });
    return arg;
}

/************************************************************************/
/*                 GDALAlgorithm::AddOutputDatasetArg()                 */
/************************************************************************/

GDALInConstructionAlgorithmArg &GDALAlgorithm::AddOutputDatasetArg(
    GDALArgDatasetValue *pValue, GDALArgDatasetType type,
    bool positionalAndRequired, const char *helpMessage)
{
    auto &arg =
        AddArg(GDAL_ARG_NAME_OUTPUT, 'o',
               MsgOrDefault(
                   helpMessage,
                   CPLSPrintf("Output %s dataset",
                              GDALAlgorithmArgDatasetTypeName(type).c_str())),
               pValue, type)
            .SetIsInput(true)
            .SetIsOutput(true)
            .SetDatasetInputFlags(GADV_NAME)
            .SetDatasetOutputFlags(GADV_OBJECT);
    if (positionalAndRequired)
        arg.SetPositional().SetRequired();

    AddValidationAction(
        [this, &arg, pValue]()
        {
            if (pValue->GetName() == "-")
                pValue->Set("/vsistdout/");

            auto outputFormatArg = GetArg(GDAL_ARG_NAME_OUTPUT_FORMAT);
            if (outputFormatArg && outputFormatArg->GetType() == GAAT_STRING &&
                (!outputFormatArg->IsExplicitlySet() ||
                 outputFormatArg->Get<std::string>().empty()) &&
                arg.IsExplicitlySet())
            {
                const auto vrtCompatible =
                    outputFormatArg->GetMetadataItem(GAAMDI_VRT_COMPATIBLE);
                if (vrtCompatible && !vrtCompatible->empty() &&
                    vrtCompatible->front() == "false" &&
                    EQUAL(
                        CPLGetExtensionSafe(pValue->GetName().c_str()).c_str(),
                        "VRT"))
                {
                    ReportError(
                        CE_Failure, CPLE_NotSupported,
                        "VRT output is not supported.%s",
                        outputFormatArg->GetDescription().find("GDALG") !=
                                std::string::npos
                            ? " Consider using the GDALG driver instead (files "
                              "with .gdalg.json extension)"
                            : "");
                    return false;
                }
                else if (pValue->GetName().size() > strlen(".gdalg.json") &&
                         EQUAL(pValue->GetName()
                                   .substr(pValue->GetName().size() -
                                           strlen(".gdalg.json"))
                                   .c_str(),
                               ".gdalg.json") &&
                         outputFormatArg->GetDescription().find("GDALG") ==
                             std::string::npos)
                {
                    ReportError(CE_Failure, CPLE_NotSupported,
                                "GDALG output is not supported");
                    return false;
                }
            }
            return true;
        });

    return arg;
}

/************************************************************************/
/*                   GDALAlgorithm::AddOverwriteArg()                   */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddOverwriteArg(bool *pValue, const char *helpMessage)
{
    return AddArg(
               GDAL_ARG_NAME_OVERWRITE, 0,
               MsgOrDefault(
                   helpMessage,
                   _("Whether overwriting existing output dataset is allowed")),
               pValue)
        .SetDefault(false);
}

/************************************************************************/
/*                GDALAlgorithm::AddOverwriteLayerArg()                 */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddOverwriteLayerArg(bool *pValue, const char *helpMessage)
{
    AddValidationAction(
        [this]
        {
            auto updateArg = GetArg(GDAL_ARG_NAME_UPDATE);
            if (!(updateArg && updateArg->GetType() == GAAT_BOOLEAN))
            {
                ReportError(CE_Failure, CPLE_AppDefined,
                            "--update argument must exist for "
                            "--overwrite-layer, even if hidden");
                return false;
            }
            return true;
        });
    return AddArg(
               GDAL_ARG_NAME_OVERWRITE_LAYER, 0,
               MsgOrDefault(
                   helpMessage,
                   _("Whether overwriting existing output layer is allowed")),
               pValue)
        .SetDefault(false)
        .AddAction(
            [this]
            {
                auto updateArg = GetArg(GDAL_ARG_NAME_UPDATE);
                if (updateArg && updateArg->GetType() == GAAT_BOOLEAN)
                {
                    updateArg->Set(true);
                }
            });
}

/************************************************************************/
/*                    GDALAlgorithm::AddUpdateArg()                     */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddUpdateArg(bool *pValue, const char *helpMessage)
{
    return AddArg(GDAL_ARG_NAME_UPDATE, 0,
                  MsgOrDefault(
                      helpMessage,
                      _("Whether to open existing dataset in update mode")),
                  pValue)
        .SetDefault(false);
}

/************************************************************************/
/*                  GDALAlgorithm::AddAppendLayerArg()                  */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddAppendLayerArg(bool *pValue, const char *helpMessage)
{
    AddValidationAction(
        [this]
        {
            auto updateArg = GetArg(GDAL_ARG_NAME_UPDATE);
            if (!(updateArg && updateArg->GetType() == GAAT_BOOLEAN))
            {
                ReportError(CE_Failure, CPLE_AppDefined,
                            "--update argument must exist for --append, even "
                            "if hidden");
                return false;
            }
            return true;
        });
    return AddArg(GDAL_ARG_NAME_APPEND, 0,
                  MsgOrDefault(
                      helpMessage,
                      _("Whether appending to existing layer is allowed")),
                  pValue)
        .SetDefault(false)
        .AddAction(
            [this]
            {
                auto updateArg = GetArg(GDAL_ARG_NAME_UPDATE);
                if (updateArg && updateArg->GetType() == GAAT_BOOLEAN)
                {
                    updateArg->Set(true);
                }
            });
}

/************************************************************************/
/*                GDALAlgorithm::AddOptionsSuggestions()                */
/************************************************************************/

/* static */
bool GDALAlgorithm::AddOptionsSuggestions(const char *pszXML, int datasetType,
                                          const std::string &currentValue,
                                          std::vector<std::string> &oRet)
{
    if (!pszXML)
        return false;
    CPLXMLTreeCloser poTree(CPLParseXMLString(pszXML));
    if (!poTree)
        return false;

    std::string typedOptionName = currentValue;
    const auto posEqual = typedOptionName.find('=');
    std::string typedValue;
    if (posEqual != 0 && posEqual != std::string::npos)
    {
        typedValue = currentValue.substr(posEqual + 1);
        typedOptionName.resize(posEqual);
    }

    for (const CPLXMLNode *psChild = poTree.get()->psChild; psChild;
         psChild = psChild->psNext)
    {
        const char *pszName = CPLGetXMLValue(psChild, "name", nullptr);
        if (pszName && typedOptionName == pszName &&
            (strcmp(psChild->pszValue, "Option") == 0 ||
             strcmp(psChild->pszValue, "Argument") == 0))
        {
            const char *pszType = CPLGetXMLValue(psChild, "type", "");
            const char *pszMin = CPLGetXMLValue(psChild, "min", nullptr);
            const char *pszMax = CPLGetXMLValue(psChild, "max", nullptr);
            if (EQUAL(pszType, "string-select"))
            {
                for (const CPLXMLNode *psChild2 = psChild->psChild; psChild2;
                     psChild2 = psChild2->psNext)
                {
                    if (EQUAL(psChild2->pszValue, "Value"))
                    {
                        oRet.push_back(CPLGetXMLValue(psChild2, "", ""));
                    }
                }
            }
            else if (EQUAL(pszType, "boolean"))
            {
                if (typedValue == "YES" || typedValue == "NO")
                {
                    oRet.push_back(currentValue);
                    return true;
                }
                oRet.push_back("NO");
                oRet.push_back("YES");
            }
            else if (EQUAL(pszType, "int"))
            {
                if (pszMin && pszMax && atoi(pszMax) - atoi(pszMin) > 0 &&
                    atoi(pszMax) - atoi(pszMin) < 25)
                {
                    const int nMax = atoi(pszMax);
                    for (int i = atoi(pszMin); i <= nMax; ++i)
                        oRet.push_back(std::to_string(i));
                }
            }

            if (oRet.empty())
            {
                if (pszMin && pszMax)
                {
                    oRet.push_back(std::string("##"));
                    oRet.push_back(std::string("validity range: [")
                                       .append(pszMin)
                                       .append(",")
                                       .append(pszMax)
                                       .append("]"));
                }
                else if (pszMin)
                {
                    oRet.push_back(std::string("##"));
                    oRet.push_back(
                        std::string("validity range: >= ").append(pszMin));
                }
                else if (pszMax)
                {
                    oRet.push_back(std::string("##"));
                    oRet.push_back(
                        std::string("validity range: <= ").append(pszMax));
                }
                else if (const char *pszDescription =
                             CPLGetXMLValue(psChild, "description", nullptr))
                {
                    oRet.push_back(std::string("##"));
                    oRet.push_back(std::string("type: ")
                                       .append(pszType)
                                       .append(", description: ")
                                       .append(pszDescription));
                }
            }

            return true;
        }
    }

    for (const CPLXMLNode *psChild = poTree.get()->psChild; psChild;
         psChild = psChild->psNext)
    {
        const char *pszName = CPLGetXMLValue(psChild, "name", nullptr);
        if (pszName && (strcmp(psChild->pszValue, "Option") == 0 ||
                        strcmp(psChild->pszValue, "Argument") == 0))
        {
            const char *pszScope = CPLGetXMLValue(psChild, "scope", nullptr);
            if (!pszScope ||
                (EQUAL(pszScope, "raster") &&
                 (datasetType & GDAL_OF_RASTER) != 0) ||
                (EQUAL(pszScope, "vector") &&
                 (datasetType & GDAL_OF_VECTOR) != 0))
            {
                oRet.push_back(std::string(pszName).append("="));
            }
        }
    }

    return false;
}

/************************************************************************/
/*             GDALAlgorithm::OpenOptionCompleteFunction()              */
/************************************************************************/

//! @cond Doxygen_Suppress
std::vector<std::string>
GDALAlgorithm::OpenOptionCompleteFunction(const std::string &currentValue) const
{
    std::vector<std::string> oRet;

    int datasetType = GDAL_OF_RASTER | GDAL_OF_VECTOR | GDAL_OF_MULTIDIM_RASTER;
    auto inputArg = GetArg(GDAL_ARG_NAME_INPUT);
    if (inputArg && (inputArg->GetType() == GAAT_DATASET ||
                     inputArg->GetType() == GAAT_DATASET_LIST))
    {
        datasetType = inputArg->GetDatasetType();
    }

    auto inputFormat = GetArg(GDAL_ARG_NAME_INPUT_FORMAT);
    if (inputFormat && inputFormat->GetType() == GAAT_STRING_LIST &&
        inputFormat->IsExplicitlySet())
    {
        const auto &aosAllowedDrivers =
            inputFormat->Get<std::vector<std::string>>();
        if (aosAllowedDrivers.size() == 1)
        {
            auto poDriver = GetGDALDriverManager()->GetDriverByName(
                aosAllowedDrivers[0].c_str());
            if (poDriver)
            {
                AddOptionsSuggestions(
                    poDriver->GetMetadataItem(GDAL_DMD_OPENOPTIONLIST),
                    datasetType, currentValue, oRet);
            }
            return oRet;
        }
    }

    const auto AddSuggestions = [datasetType, &currentValue,
                                 &oRet](const GDALArgDatasetValue &datasetValue)
    {
        auto poDM = GetGDALDriverManager();

        const auto &osDSName = datasetValue.GetName();
        const std::string osExt = CPLGetExtensionSafe(osDSName.c_str());
        if (!osExt.empty())
        {
            std::set<std::string> oVisitedExtensions;
            for (int i = 0; i < poDM->GetDriverCount(); ++i)
            {
                auto poDriver = poDM->GetDriver(i);
                if (((datasetType & GDAL_OF_RASTER) != 0 &&
                     poDriver->GetMetadataItem(GDAL_DCAP_RASTER)) ||
                    ((datasetType & GDAL_OF_VECTOR) != 0 &&
                     poDriver->GetMetadataItem(GDAL_DCAP_VECTOR)) ||
                    ((datasetType & GDAL_OF_MULTIDIM_RASTER) != 0 &&
                     poDriver->GetMetadataItem(GDAL_DCAP_MULTIDIM_RASTER)))
                {
                    const char *pszExtensions =
                        poDriver->GetMetadataItem(GDAL_DMD_EXTENSIONS);
                    if (pszExtensions)
                    {
                        const CPLStringList aosExts(
                            CSLTokenizeString2(pszExtensions, " ", 0));
                        for (const char *pszExt : cpl::Iterate(aosExts))
                        {
                            if (EQUAL(pszExt, osExt.c_str()) &&
                                !cpl::contains(oVisitedExtensions, pszExt))
                            {
                                oVisitedExtensions.insert(pszExt);
                                if (AddOptionsSuggestions(
                                        poDriver->GetMetadataItem(
                                            GDAL_DMD_OPENOPTIONLIST),
                                        datasetType, currentValue, oRet))
                                {
                                    return;
                                }
                                break;
                            }
                        }
                    }
                }
            }
        }
    };

    if (inputArg && inputArg->GetType() == GAAT_DATASET)
    {
        auto &datasetValue = inputArg->Get<GDALArgDatasetValue>();
        AddSuggestions(datasetValue);
    }
    else if (inputArg && inputArg->GetType() == GAAT_DATASET_LIST)
    {
        auto &datasetValues = inputArg->Get<std::vector<GDALArgDatasetValue>>();
        if (datasetValues.size() == 1)
            AddSuggestions(datasetValues[0]);
    }

    return oRet;
}

//! @endcond

/************************************************************************/
/*                  GDALAlgorithm::AddOpenOptionsArg()                  */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddOpenOptionsArg(std::vector<std::string> *pValue,
                                 const char *helpMessage)
{
    auto &arg = AddArg(GDAL_ARG_NAME_OPEN_OPTION, 0,
                       MsgOrDefault(helpMessage, _("Open options")), pValue)
                    .AddAlias("oo")
                    .SetMetaVar("<KEY>=<VALUE>")
                    .SetPackedValuesAllowed(false)
                    .SetCategory(GAAC_ADVANCED);

    arg.AddValidationAction([this, &arg]()
                            { return ParseAndValidateKeyValue(arg); });

    arg.SetAutoCompleteFunction(
        [this](const std::string &currentValue)
        { return OpenOptionCompleteFunction(currentValue); });

    return arg;
}

/************************************************************************/
/*               GDALAlgorithm::AddOutputOpenOptionsArg()               */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddOutputOpenOptionsArg(std::vector<std::string> *pValue,
                                       const char *helpMessage)
{
    auto &arg =
        AddArg(GDAL_ARG_NAME_OUTPUT_OPEN_OPTION, 0,
               MsgOrDefault(helpMessage, _("Output open options")), pValue)
            .AddAlias("output-oo")
            .SetMetaVar("<KEY>=<VALUE>")
            .SetPackedValuesAllowed(false)
            .SetCategory(GAAC_ADVANCED);

    arg.AddValidationAction([this, &arg]()
                            { return ParseAndValidateKeyValue(arg); });

    arg.SetAutoCompleteFunction(
        [this](const std::string &currentValue)
        { return OpenOptionCompleteFunction(currentValue); });

    return arg;
}

/************************************************************************/
/*                           ValidateFormat()                           */
/************************************************************************/

bool GDALAlgorithm::ValidateFormat(const GDALAlgorithmArg &arg,
                                   bool bStreamAllowed,
                                   bool bGDALGAllowed) const
{
    if (arg.GetChoices().empty())
    {
        const auto Validate =
            [this, &arg, bStreamAllowed, bGDALGAllowed](const std::string &val)
        {
            if (const auto extraFormats =
                    arg.GetMetadataItem(GAAMDI_EXTRA_FORMATS))
            {
                for (const auto &extraFormat : *extraFormats)
                {
                    if (EQUAL(val.c_str(), extraFormat.c_str()))
                        return true;
                }
            }

            if (bStreamAllowed && EQUAL(val.c_str(), "stream"))
                return true;

            if (EQUAL(val.c_str(), "GDALG") &&
                arg.GetName() == GDAL_ARG_NAME_OUTPUT_FORMAT)
            {
                if (bGDALGAllowed)
                {
                    return true;
                }
                else
                {
                    ReportError(CE_Failure, CPLE_NotSupported,
                                "GDALG output is not supported.");
                    return false;
                }
            }

            const auto vrtCompatible =
                arg.GetMetadataItem(GAAMDI_VRT_COMPATIBLE);
            if (vrtCompatible && !vrtCompatible->empty() &&
                vrtCompatible->front() == "false" && EQUAL(val.c_str(), "VRT"))
            {
                ReportError(CE_Failure, CPLE_NotSupported,
                            "VRT output is not supported.%s",
                            bGDALGAllowed
                                ? " Consider using the GDALG driver instead "
                                  "(files with .gdalg.json extension)."
                                : "");
                return false;
            }

            const auto allowedFormats =
                arg.GetMetadataItem(GAAMDI_ALLOWED_FORMATS);
            if (allowedFormats && !allowedFormats->empty() &&
                std::find(allowedFormats->begin(), allowedFormats->end(),
                          val) != allowedFormats->end())
            {
                return true;
            }

            const auto excludedFormats =
                arg.GetMetadataItem(GAAMDI_EXCLUDED_FORMATS);
            if (excludedFormats && !excludedFormats->empty() &&
                std::find(excludedFormats->begin(), excludedFormats->end(),
                          val) != excludedFormats->end())
            {
                ReportError(CE_Failure, CPLE_NotSupported,
                            "%s output is not supported.", val.c_str());
                return false;
            }

            auto hDriver = GDALGetDriverByName(val.c_str());
            if (!hDriver)
            {
                auto poMissingDriver =
                    GetGDALDriverManager()->GetHiddenDriverByName(val.c_str());
                if (poMissingDriver)
                {
                    const std::string msg =
                        GDALGetMessageAboutMissingPluginDriver(poMissingDriver);
                    ReportError(CE_Failure, CPLE_AppDefined,
                                "Invalid value for argument '%s'. Driver '%s' "
                                "not found but is known. However plugin %s",
                                arg.GetName().c_str(), val.c_str(),
                                msg.c_str());
                }
                else
                {
                    ReportError(CE_Failure, CPLE_AppDefined,
                                "Invalid value for argument '%s'. Driver '%s' "
                                "does not exist.",
                                arg.GetName().c_str(), val.c_str());
                }
                return false;
            }

            const auto caps = arg.GetMetadataItem(GAAMDI_REQUIRED_CAPABILITIES);
            if (caps)
            {
                for (const std::string &cap : *caps)
                {
                    const char *pszVal =
                        GDALGetMetadataItem(hDriver, cap.c_str(), nullptr);
                    if (!(pszVal && pszVal[0]))
                    {
                        if (cap == GDAL_DCAP_CREATECOPY &&
                            std::find(caps->begin(), caps->end(),
                                      GDAL_DCAP_RASTER) != caps->end() &&
                            GDALGetMetadataItem(hDriver, GDAL_DCAP_RASTER,
                                                nullptr) &&
                            GDALGetMetadataItem(hDriver, GDAL_DCAP_CREATE,
                                                nullptr))
                        {
                            // if it supports Create, it supports CreateCopy
                        }
                        else if (cap == GDAL_DMD_EXTENSIONS)
                        {
                            ReportError(
                                CE_Failure, CPLE_AppDefined,
                                "Invalid value for argument '%s'. Driver '%s' "
                                "does "
                                "not advertise any file format extension.",
                                arg.GetName().c_str(), val.c_str());
                            return false;
                        }
                        else
                        {
                            if (cap == GDAL_DCAP_CREATE)
                            {
                                auto updateArg = GetArg(GDAL_ARG_NAME_UPDATE);
                                if (updateArg &&
                                    updateArg->GetType() == GAAT_BOOLEAN &&
                                    updateArg->IsExplicitlySet())
                                {
                                    continue;
                                }

                                ReportError(
                                    CE_Failure, CPLE_AppDefined,
                                    "Invalid value for argument '%s'. "
                                    "Driver '%s' does not have write support.",
                                    arg.GetName().c_str(), val.c_str());
                                return false;
                            }
                            else
                            {
                                ReportError(
                                    CE_Failure, CPLE_AppDefined,
                                    "Invalid value for argument '%s'. Driver "
                                    "'%s' "
                                    "does "
                                    "not expose the required '%s' capability.",
                                    arg.GetName().c_str(), val.c_str(),
                                    cap.c_str());
                                return false;
                            }
                        }
                    }
                }
            }
            return true;
        };

        if (arg.GetType() == GAAT_STRING)
        {
            return Validate(arg.Get<std::string>());
        }
        else if (arg.GetType() == GAAT_STRING_LIST)
        {
            for (const auto &val : arg.Get<std::vector<std::string>>())
            {
                if (!Validate(val))
                    return false;
            }
        }
    }

    return true;
}

/************************************************************************/
/*                     FormatAutoCompleteFunction()                     */
/************************************************************************/

/* static */
std::vector<std::string> GDALAlgorithm::FormatAutoCompleteFunction(
    const GDALAlgorithmArg &arg, bool /* bStreamAllowed */, bool bGDALGAllowed)
{
    std::vector<std::string> res;
    auto poDM = GetGDALDriverManager();
    const auto vrtCompatible = arg.GetMetadataItem(GAAMDI_VRT_COMPATIBLE);
    const auto allowedFormats = arg.GetMetadataItem(GAAMDI_ALLOWED_FORMATS);
    const auto excludedFormats = arg.GetMetadataItem(GAAMDI_EXCLUDED_FORMATS);
    const auto caps = arg.GetMetadataItem(GAAMDI_REQUIRED_CAPABILITIES);
    if (auto extraFormats = arg.GetMetadataItem(GAAMDI_EXTRA_FORMATS))
        res = std::move(*extraFormats);
    for (int i = 0; i < poDM->GetDriverCount(); ++i)
    {
        auto poDriver = poDM->GetDriver(i);

        if (vrtCompatible && !vrtCompatible->empty() &&
            vrtCompatible->front() == "false" &&
            EQUAL(poDriver->GetDescription(), "VRT"))
        {
            // do nothing
        }
        else if (allowedFormats && !allowedFormats->empty() &&
                 std::find(allowedFormats->begin(), allowedFormats->end(),
                           poDriver->GetDescription()) != allowedFormats->end())
        {
            res.push_back(poDriver->GetDescription());
        }
        else if (excludedFormats && !excludedFormats->empty() &&
                 std::find(excludedFormats->begin(), excludedFormats->end(),
                           poDriver->GetDescription()) !=
                     excludedFormats->end())
        {
            continue;
        }
        else if (caps)
        {
            bool ok = true;
            for (const std::string &cap : *caps)
            {
                if (cap == GDAL_ALG_DCAP_RASTER_OR_MULTIDIM_RASTER)
                {
                    if (!poDriver->GetMetadataItem(GDAL_DCAP_RASTER) &&
                        !poDriver->GetMetadataItem(GDAL_DCAP_MULTIDIM_RASTER))
                    {
                        ok = false;
                        break;
                    }
                }
                else if (const char *pszVal =
                             poDriver->GetMetadataItem(cap.c_str());
                         pszVal && pszVal[0])
                {
                }
                else if (cap == GDAL_DCAP_CREATECOPY &&
                         (std::find(caps->begin(), caps->end(),
                                    GDAL_DCAP_RASTER) != caps->end() &&
                          poDriver->GetMetadataItem(GDAL_DCAP_RASTER)) &&
                         poDriver->GetMetadataItem(GDAL_DCAP_CREATE))
                {
                    // if it supports Create, it supports CreateCopy
                }
                else
                {
                    ok = false;
                    break;
                }
            }
            if (ok)
            {
                res.push_back(poDriver->GetDescription());
            }
        }
    }
    if (bGDALGAllowed)
        res.push_back("GDALG");
    return res;
}

/************************************************************************/
/*                 GDALAlgorithm::AddInputFormatsArg()                  */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddInputFormatsArg(std::vector<std::string> *pValue,
                                  const char *helpMessage)
{
    auto &arg = AddArg(GDAL_ARG_NAME_INPUT_FORMAT, 0,
                       MsgOrDefault(helpMessage, _("Input formats")), pValue)
                    .AddAlias("if")
                    .SetCategory(GAAC_ADVANCED);
    arg.AddValidationAction([this, &arg]()
                            { return ValidateFormat(arg, false, false); });
    arg.SetAutoCompleteFunction(
        [&arg](const std::string &)
        { return FormatAutoCompleteFunction(arg, false, false); });
    return arg;
}

/************************************************************************/
/*                 GDALAlgorithm::AddOutputFormatArg()                  */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddOutputFormatArg(std::string *pValue, bool bStreamAllowed,
                                  bool bGDALGAllowed, const char *helpMessage)
{
    auto &arg = AddArg(GDAL_ARG_NAME_OUTPUT_FORMAT, 'f',
                       MsgOrDefault(helpMessage,
                                    bGDALGAllowed
                                        ? _("Output format (\"GDALG\" allowed)")
                                        : _("Output format")),
                       pValue)
                    .AddAlias("of")
                    .AddAlias("format");
    arg.AddValidationAction(
        [this, &arg, bStreamAllowed, bGDALGAllowed]()
        { return ValidateFormat(arg, bStreamAllowed, bGDALGAllowed); });
    arg.SetAutoCompleteFunction(
        [&arg, bStreamAllowed, bGDALGAllowed](const std::string &)
        {
            return FormatAutoCompleteFunction(arg, bStreamAllowed,
                                              bGDALGAllowed);
        });
    return arg;
}

/************************************************************************/
/*                GDALAlgorithm::AddOutputDataTypeArg()                 */
/************************************************************************/
GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddOutputDataTypeArg(std::string *pValue,
                                    const char *helpMessage)
{
    auto &arg =
        AddArg(GDAL_ARG_NAME_OUTPUT_DATA_TYPE, 0,
               MsgOrDefault(helpMessage, _("Output data type")), pValue)
            .AddAlias("ot")
            .AddAlias("datatype")
            .AddMetadataItem("type", {"GDALDataType"})
            .SetChoices("UInt8", "Int8", "UInt16", "Int16", "UInt32", "Int32",
                        "UInt64", "Int64", "CInt16", "CInt32", "Float16",
                        "Float32", "Float64", "CFloat32", "CFloat64")
            .SetHiddenChoices("Byte");
    return arg;
}

/************************************************************************/
/*                    GDALAlgorithm::AddNodataArg()                     */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddNodataArg(std::string *pValue, bool noneAllowed,
                            const std::string &optionName,
                            const char *helpMessage)
{
    auto &arg = AddArg(
        optionName, 0,
        MsgOrDefault(helpMessage,
                     noneAllowed
                         ? _("Assign a specified nodata value to output bands "
                             "('none', numeric value, 'nan', 'inf', '-inf')")
                         : _("Assign a specified nodata value to output bands "
                             "(numeric value, 'nan', 'inf', '-inf')")),
        pValue);
    arg.AddValidationAction(
        [this, pValue, noneAllowed, optionName]()
        {
            if (!(noneAllowed && EQUAL(pValue->c_str(), "none")))
            {
                char *endptr = nullptr;
                CPLStrtod(pValue->c_str(), &endptr);
                if (endptr != pValue->c_str() + pValue->size())
                {
                    ReportError(CE_Failure, CPLE_IllegalArg,
                                "Value of '%s' should be %sa "
                                "numeric value, 'nan', 'inf' or '-inf'",
                                optionName.c_str(),
                                noneAllowed ? "'none', " : "");
                    return false;
                }
            }
            return true;
        });
    return arg;
}

/************************************************************************/
/*                 GDALAlgorithm::AddOutputStringArg()                  */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddOutputStringArg(std::string *pValue, const char *helpMessage)
{
    return AddArg(
               GDAL_ARG_NAME_OUTPUT_STRING, 0,
               MsgOrDefault(helpMessage,
                            _("Output string, in which the result is placed")),
               pValue)
        .SetHiddenForCLI()
        .SetIsInput(false)
        .SetIsOutput(true);
}

/************************************************************************/
/*                    GDALAlgorithm::AddStdoutArg()                     */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddStdoutArg(bool *pValue, const char *helpMessage)
{
    return AddArg(GDAL_ARG_NAME_STDOUT, 0,
                  MsgOrDefault(helpMessage,
                               _("Directly output on stdout. If enabled, "
                                 "output-string will be empty")),
                  pValue)
        .SetHidden();
}

/************************************************************************/
/*                   GDALAlgorithm::AddLayerNameArg()                   */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddLayerNameArg(std::string *pValue, const char *helpMessage)
{
    return AddArg(GDAL_ARG_NAME_INPUT_LAYER, 'l',
                  MsgOrDefault(helpMessage, _("Input layer name")), pValue);
}

/************************************************************************/
/*                   GDALAlgorithm::AddArrayNameArg()                   */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddArrayNameArg(std::string *pValue, const char *helpMessage)
{
    return AddArg("array", 0, MsgOrDefault(helpMessage, _("Array name")),
                  pValue)
        .SetAutoCompleteFunction([this](const std::string &)
                                 { return AutoCompleteArrayName(); });
}

/************************************************************************/
/*                   GDALAlgorithm::AddArrayNameArg()                   */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddArrayNameArg(std::vector<std::string> *pValue,
                               const char *helpMessage)
{
    return AddArg("array", 0, MsgOrDefault(helpMessage, _("Array name(s)")),
                  pValue)
        .SetAutoCompleteFunction([this](const std::string &)
                                 { return AutoCompleteArrayName(); });
}

/************************************************************************/
/*                GDALAlgorithm::AutoCompleteArrayName()                */
/************************************************************************/

std::vector<std::string> GDALAlgorithm::AutoCompleteArrayName() const
{
    std::vector<std::string> ret;
    std::string osDSName;
    auto inputArg = GetArg(GDAL_ARG_NAME_INPUT);
    if (inputArg && inputArg->GetType() == GAAT_DATASET_LIST)
    {
        auto &inputDatasets = inputArg->Get<std::vector<GDALArgDatasetValue>>();
        if (!inputDatasets.empty())
        {
            osDSName = inputDatasets[0].GetName();
        }
    }
    else if (inputArg && inputArg->GetType() == GAAT_DATASET)
    {
        auto &inputDataset = inputArg->Get<GDALArgDatasetValue>();
        osDSName = inputDataset.GetName();
    }

    if (!osDSName.empty())
    {
        CPLStringList aosAllowedDrivers;
        const auto ifArg = GetArg(GDAL_ARG_NAME_INPUT_FORMAT);
        if (ifArg && ifArg->GetType() == GAAT_STRING_LIST)
            aosAllowedDrivers =
                CPLStringList(ifArg->Get<std::vector<std::string>>());

        CPLStringList aosOpenOptions;
        const auto ooArg = GetArg(GDAL_ARG_NAME_OPEN_OPTION);
        if (ooArg && ooArg->GetType() == GAAT_STRING_LIST)
            aosOpenOptions =
                CPLStringList(ooArg->Get<std::vector<std::string>>());

        if (auto poDS = std::unique_ptr<GDALDataset>(GDALDataset::Open(
                osDSName.c_str(), GDAL_OF_MULTIDIM_RASTER,
                aosAllowedDrivers.List(), aosOpenOptions.List(), nullptr)))
        {
            if (auto poRG = poDS->GetRootGroup())
            {
                ret = poRG->GetMDArrayFullNamesRecursive();
            }
        }
    }

    return ret;
}

/************************************************************************/
/*                  GDALAlgorithm::AddMemorySizeArg()                   */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddMemorySizeArg(size_t *pValue, std::string *pStrValue,
                                const std::string &optionName,
                                const char *helpMessage)
{
    return AddArg(optionName, 0, helpMessage, pStrValue)
        .SetDefault(*pStrValue)
        .AddValidationAction(
            [this, pValue, pStrValue]()
            {
                CPLDebug("GDAL", "StrValue `%s`", pStrValue->c_str());
                GIntBig nBytes;
                bool bUnitSpecified;
                if (CPLParseMemorySize(pStrValue->c_str(), &nBytes,
                                       &bUnitSpecified) != CE_None)
                {
                    return false;
                }
                if (!bUnitSpecified)
                {
                    ReportError(CE_Failure, CPLE_AppDefined,
                                "Memory size must have a unit or be a "
                                "percentage of usable RAM (2GB, 5%%, etc.)");
                    return false;
                }
                if constexpr (sizeof(std::uint64_t) > sizeof(size_t))
                {
                    // -1 to please CoverityScan
                    if (static_cast<std::uint64_t>(nBytes) >
                        std::numeric_limits<size_t>::max() - 1U)
                    {
                        ReportError(CE_Failure, CPLE_AppDefined,
                                    "Memory size %s is too large.",
                                    pStrValue->c_str());
                        return false;
                    }
                }

                *pValue = static_cast<size_t>(nBytes);
                return true;
            });
}

/************************************************************************/
/*                GDALAlgorithm::AddOutputLayerNameArg()                */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddOutputLayerNameArg(std::string *pValue,
                                     const char *helpMessage)
{
    return AddArg(GDAL_ARG_NAME_OUTPUT_LAYER, 0,
                  MsgOrDefault(helpMessage, _("Output layer name")), pValue);
}

/************************************************************************/
/*                   GDALAlgorithm::AddLayerNameArg()                   */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddLayerNameArg(std::vector<std::string> *pValue,
                               const char *helpMessage)
{
    return AddArg(GDAL_ARG_NAME_INPUT_LAYER, 'l',
                  MsgOrDefault(helpMessage, _("Input layer name")), pValue);
}

/************************************************************************/
/*                 GDALAlgorithm::AddGeometryTypeArg()                  */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddGeometryTypeArg(std::string *pValue, const char *helpMessage)
{
    return AddArg("geometry-type", 0,
                  MsgOrDefault(helpMessage, _("Geometry type")), pValue)
        .SetAutoCompleteFunction(
            [](const std::string &currentValue)
            {
                std::vector<std::string> oRet;
                for (const char *type :
                     {"GEOMETRY", "POINT", "LINESTRING", "POLYGON",
                      "MULTIPOINT", "MULTILINESTRING", "MULTIPOLYGON",
                      "GEOMETRYCOLLECTION", "CURVE", "CIRCULARSTRING",
                      "COMPOUNDCURVE", "SURFACE", "CURVEPOLYGON", "MULTICURVE",
                      "MULTISURFACE", "POLYHEDRALSURFACE", "TIN"})
                {
                    if (currentValue.empty() ||
                        STARTS_WITH(type, currentValue.c_str()))
                    {
                        oRet.push_back(type);
                        oRet.push_back(std::string(type).append("Z"));
                        oRet.push_back(std::string(type).append("M"));
                        oRet.push_back(std::string(type).append("ZM"));
                    }
                }
                return oRet;
            })
        .AddValidationAction(
            [this, pValue]()
            {
                if (wkbFlatten(OGRFromOGCGeomType(pValue->c_str())) ==
                        wkbUnknown &&
                    !STARTS_WITH_CI(pValue->c_str(), "GEOMETRY"))
                {
                    ReportError(CE_Failure, CPLE_AppDefined,
                                "Invalid geometry type '%s'", pValue->c_str());
                    return false;
                }
                return true;
            });
}

/************************************************************************/
/*         GDALAlgorithm::SetAutoCompleteFunctionForLayerName()         */
/************************************************************************/

/* static */
void GDALAlgorithm::SetAutoCompleteFunctionForLayerName(
    GDALInConstructionAlgorithmArg &layerArg, GDALAlgorithmArg &datasetArg)
{
    CPLAssert(datasetArg.GetType() == GAAT_DATASET ||
              datasetArg.GetType() == GAAT_DATASET_LIST);

    layerArg.SetAutoCompleteFunction(
        [&datasetArg](const std::string &currentValue)
        {
            std::vector<std::string> ret;
            CPLErrorStateBackuper oBackuper(CPLQuietErrorHandler);
            GDALArgDatasetValue *dsVal = nullptr;
            if (datasetArg.GetType() == GAAT_DATASET)
            {
                dsVal = &(datasetArg.Get<GDALArgDatasetValue>());
            }
            else
            {
                auto &val = datasetArg.Get<std::vector<GDALArgDatasetValue>>();
                if (val.size() == 1)
                {
                    dsVal = &val[0];
                }
            }
            if (dsVal && !dsVal->GetName().empty())
            {
                auto poDS = std::unique_ptr<GDALDataset>(GDALDataset::Open(
                    dsVal->GetName().c_str(), GDAL_OF_VECTOR));
                if (poDS)
                {
                    for (auto &&poLayer : poDS->GetLayers())
                    {
                        if (currentValue == poLayer->GetDescription())
                        {
                            ret.clear();
                            ret.push_back(poLayer->GetDescription());
                            break;
                        }
                        ret.push_back(poLayer->GetDescription());
                    }
                }
            }
            return ret;
        });
}

/************************************************************************/
/*         GDALAlgorithm::SetAutoCompleteFunctionForFieldName()         */
/************************************************************************/

void GDALAlgorithm::SetAutoCompleteFunctionForFieldName(
    GDALInConstructionAlgorithmArg &fieldArg,
    const GDALAlgorithmArg *layerNameArg, bool attributeFields,
    bool geometryFields, std::vector<GDALArgDatasetValue> &datasetArg)
{

    fieldArg.SetAutoCompleteFunction(
        [&datasetArg, layerNameArg, attributeFields,
         geometryFields](const std::string &currentValue)
        {
            std::set<std::string> ret;
            if (!datasetArg.empty())
            {
                CPLErrorStateBackuper oBackuper(CPLQuietErrorHandler);

                const auto getLayerFields =
                    [&ret, &currentValue, attributeFields,
                     geometryFields](const OGRLayer *poLayer)
                {
                    const auto poDefn = poLayer->GetLayerDefn();
                    if (attributeFields)
                    {
                        for (const auto poFieldDefn : poDefn->GetFields())
                        {
                            const char *fieldName = poFieldDefn->GetNameRef();
                            if (currentValue == fieldName)
                            {
                                ret.clear();
                                ret.insert(fieldName);
                                break;
                            }
                            ret.insert(fieldName);
                        }
                    }
                    if (geometryFields)
                    {
                        for (const auto poFieldDefn : poDefn->GetGeomFields())
                        {
                            const char *fieldName = poFieldDefn->GetNameRef();
                            if (fieldName[0] == 0)
                                fieldName = "OGR_GEOMETRY";
                            if (currentValue == fieldName)
                            {
                                ret.clear();
                                ret.insert(fieldName);
                                break;
                            }
                            ret.insert(fieldName);
                        }
                    }
                };

                const GDALArgDatasetValue &dsVal = datasetArg[0];

                if (!dsVal.GetName().empty())
                {
                    auto poDS = std::unique_ptr<GDALDataset>(
                        GDALDataset::Open(dsVal.GetName().c_str(),
                                          GDAL_OF_VECTOR | GDAL_OF_READONLY));
                    if (poDS)
                    {
                        std::vector<std::string> layerNames;
                        if (layerNameArg && layerNameArg->IsExplicitlySet())
                        {
                            if (layerNameArg->GetType() == GAAT_STRING_LIST)
                            {
                                layerNames =
                                    layerNameArg
                                        ->Get<std::vector<std::string>>();
                            }
                            else if (layerNameArg->GetType() == GAAT_STRING)
                            {
                                layerNames.push_back(
                                    layerNameArg->Get<std::string>());
                            }
                        }
                        if (layerNames.empty())
                        {
                            // Loop through all layers
                            for (const auto *poLayer : poDS->GetLayers())
                            {
                                getLayerFields(poLayer);
                            }
                        }
                        else
                        {
                            for (const std::string &layerName : layerNames)
                            {
                                const auto poLayer =
                                    poDS->GetLayerByName(layerName.c_str());
                                if (poLayer)
                                {
                                    getLayerFields(poLayer);
                                }
                            }
                        }
                    }
                }
            }
            std::vector<std::string> retVector(ret.begin(), ret.end());
            return retVector;
        });
}

/************************************************************************/
/*                   GDALAlgorithm::AddFieldNameArg()                   */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddFieldNameArg(std::string *pValue, const char *helpMessage)
{
    return AddArg("field-name", 0, MsgOrDefault(helpMessage, _("Field name")),
                  pValue);
}

/************************************************************************/
/*                GDALAlgorithm::ParseFieldDefinition()                 */
/************************************************************************/
bool GDALAlgorithm::ParseFieldDefinition(const std::string &posStrDef,
                                         OGRFieldDefn *poFieldDefn,
                                         std::string *posError)
{
    static const std::regex re(
        R"(^([^:]+):([^(\s]+)(?:\((\d+)(?:,(\d+))?\))?$)");
    std::smatch match;
    if (std::regex_match(posStrDef, match, re))
    {
        const std::string name = match[1];
        const std::string type = match[2];
        const int width = match[3].matched ? std::stoi(match[3]) : 0;
        const int precision = match[4].matched ? std::stoi(match[4]) : 0;
        poFieldDefn->SetName(name.c_str());

        const auto typeEnum{OGRFieldDefn::GetFieldTypeByName(type.c_str())};
        if (typeEnum == OFTString && !EQUAL(type.c_str(), "String"))
        {
            if (posError)
                *posError = "Unsupported field type: " + type;

            return false;
        }
        poFieldDefn->SetType(typeEnum);
        poFieldDefn->SetWidth(width);
        poFieldDefn->SetPrecision(precision);
        return true;
    }

    if (posError)
        *posError = "Invalid field definition format. Expected "
                    "<NAME>:<TYPE>[(<WIDTH>[,<PRECISION>])]";

    return false;
}

/************************************************************************/
/*                GDALAlgorithm::AddFieldDefinitionArg()                */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddFieldDefinitionArg(std::vector<std::string> *pValues,
                                     std::vector<OGRFieldDefn> *pFieldDefns,
                                     const char *helpMessage)
{
    auto &arg =
        AddArg("field", 0, MsgOrDefault(helpMessage, _("Field definition")),
               pValues)
            .SetMetaVar("<NAME>:<TYPE>[(<WIDTH>[,<PRECISION>])]")
            .SetPackedValuesAllowed(true)
            .SetRepeatedArgAllowed(true);

    auto validationFunction = [this, pFieldDefns, pValues]()
    {
        pFieldDefns->clear();
        for (const auto &strValue : *pValues)
        {
            OGRFieldDefn fieldDefn("", OFTString);
            std::string error;
            if (!GDALAlgorithm::ParseFieldDefinition(strValue, &fieldDefn,
                                                     &error))
            {
                ReportError(CE_Failure, CPLE_AppDefined, "%s", error.c_str());
                return false;
            }
            // Check uniqueness of field names
            for (const auto &existingFieldDefn : *pFieldDefns)
            {
                if (EQUAL(existingFieldDefn.GetNameRef(),
                          fieldDefn.GetNameRef()))
                {
                    ReportError(CE_Failure, CPLE_AppDefined,
                                "Duplicate field name: '%s'",
                                fieldDefn.GetNameRef());
                    return false;
                }
            }
            pFieldDefns->push_back(fieldDefn);
        }
        return true;
    };

    arg.AddValidationAction(std::move(validationFunction));

    return arg;
}

/************************************************************************/
/*               GDALAlgorithm::AddFieldTypeSubtypeArg()                */
/************************************************************************/

GDALInConstructionAlgorithmArg &GDALAlgorithm::AddFieldTypeSubtypeArg(
    OGRFieldType *pTypeValue, OGRFieldSubType *pSubtypeValue,
    std::string *pStrValue, const std::string &argName, const char *helpMessage)
{
    auto &arg =
        AddArg(argName.empty() ? std::string("field-type") : argName, 0,
               MsgOrDefault(helpMessage, _("Field type or subtype")), pStrValue)
            .SetAutoCompleteFunction(
                [](const std::string &currentValue)
                {
                    std::vector<std::string> oRet;
                    for (int i = 1; i <= OGRFieldSubType::OFSTMaxSubType; i++)
                    {
                        const char *pszSubType =
                            OGRFieldDefn::GetFieldSubTypeName(
                                static_cast<OGRFieldSubType>(i));
                        if (pszSubType != nullptr)
                        {
                            if (currentValue.empty() ||
                                STARTS_WITH(pszSubType, currentValue.c_str()))
                            {
                                oRet.push_back(pszSubType);
                            }
                        }
                    }

                    for (int i = 0; i <= OGRFieldType::OFTMaxType; i++)
                    {
                        // Skip deprecated
                        if (static_cast<OGRFieldType>(i) ==
                                OGRFieldType::OFTWideString ||
                            static_cast<OGRFieldType>(i) ==
                                OGRFieldType::OFTWideStringList)
                            continue;
                        const char *pszType = OGRFieldDefn::GetFieldTypeName(
                            static_cast<OGRFieldType>(i));
                        if (pszType != nullptr)
                        {
                            if (currentValue.empty() ||
                                STARTS_WITH(pszType, currentValue.c_str()))
                            {
                                oRet.push_back(pszType);
                            }
                        }
                    }
                    return oRet;
                });

    auto validationFunction =
        [this, &arg, pTypeValue, pSubtypeValue, pStrValue]()
    {
        bool isValid{true};
        *pTypeValue = OGRFieldDefn::GetFieldTypeByName(pStrValue->c_str());

        // String is returned for unknown types
        if (!EQUAL(pStrValue->c_str(), "String") && *pTypeValue == OFTString)
        {
            isValid = false;
        }

        *pSubtypeValue =
            OGRFieldDefn::GetFieldSubTypeByName(pStrValue->c_str());

        if (*pSubtypeValue != OFSTNone)
        {
            isValid = true;
            switch (*pSubtypeValue)
            {
                case OFSTBoolean:
                case OFSTInt16:
                {
                    *pTypeValue = OFTInteger;
                    break;
                }
                case OFSTFloat32:
                {
                    *pTypeValue = OFTReal;
                    break;
                }
                default:
                {
                    *pTypeValue = OFTString;
                    break;
                }
            }
        }

        if (!isValid)
        {
            ReportError(CE_Failure, CPLE_AppDefined,
                        "Invalid value for argument '%s': '%s'",
                        arg.GetName().c_str(), pStrValue->c_str());
        }

        return isValid;
    };

    if (!pStrValue->empty())
    {
        arg.SetDefault(*pStrValue);
        validationFunction();
    }

    arg.AddValidationAction(std::move(validationFunction));

    return arg;
}

/************************************************************************/
/*                   GDALAlgorithm::ValidateBandArg()                   */
/************************************************************************/

bool GDALAlgorithm::ValidateBandArg() const
{
    bool ret = true;
    const auto bandArg = GetArg(GDAL_ARG_NAME_BAND);
    const auto inputDatasetArg = GetArg(GDAL_ARG_NAME_INPUT);
    if (bandArg && bandArg->IsExplicitlySet() && inputDatasetArg &&
        (inputDatasetArg->GetType() == GAAT_DATASET ||
         inputDatasetArg->GetType() == GAAT_DATASET_LIST) &&
        (inputDatasetArg->GetDatasetType() & GDAL_OF_RASTER) != 0)
    {
        const auto CheckBand = [this](const GDALDataset *poDS, int nBand)
        {
            if (nBand > poDS->GetRasterCount())
            {
                ReportError(CE_Failure, CPLE_AppDefined,
                            "Value of 'band' should be greater or equal than "
                            "1 and less or equal than %d.",
                            poDS->GetRasterCount());
                return false;
            }
            return true;
        };

        const auto ValidateForOneDataset =
            [&bandArg, &CheckBand](const GDALDataset *poDS)
        {
            bool l_ret = true;
            if (bandArg->GetType() == GAAT_INTEGER)
            {
                l_ret = CheckBand(poDS, bandArg->Get<int>());
            }
            else if (bandArg->GetType() == GAAT_INTEGER_LIST)
            {
                for (int nBand : bandArg->Get<std::vector<int>>())
                {
                    l_ret = l_ret && CheckBand(poDS, nBand);
                }
            }
            return l_ret;
        };

        if (inputDatasetArg->GetType() == GAAT_DATASET)
        {
            auto poDS =
                inputDatasetArg->Get<GDALArgDatasetValue>().GetDatasetRef();
            if (poDS && !ValidateForOneDataset(poDS))
                ret = false;
        }
        else
        {
            CPLAssert(inputDatasetArg->GetType() == GAAT_DATASET_LIST);
            for (auto &datasetValue :
                 inputDatasetArg->Get<std::vector<GDALArgDatasetValue>>())
            {
                auto poDS = datasetValue.GetDatasetRef();
                if (poDS && !ValidateForOneDataset(poDS))
                    ret = false;
            }
        }
    }
    return ret;
}

/************************************************************************/
/*            GDALAlgorithm::RunPreStepPipelineValidations()            */
/************************************************************************/

bool GDALAlgorithm::RunPreStepPipelineValidations() const
{
    return ValidateBandArg();
}

/************************************************************************/
/*                     GDALAlgorithm::AddBandArg()                      */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddBandArg(int *pValue, const char *helpMessage)
{
    AddValidationAction([this]() { return ValidateBandArg(); });

    return AddArg(GDAL_ARG_NAME_BAND, 'b',
                  MsgOrDefault(helpMessage, _("Input band (1-based index)")),
                  pValue)
        .AddValidationAction(
            [pValue]()
            {
                if (*pValue <= 0)
                {
                    CPLError(CE_Failure, CPLE_AppDefined,
                             "Value of 'band' should greater or equal to 1.");
                    return false;
                }
                return true;
            });
}

/************************************************************************/
/*                     GDALAlgorithm::AddBandArg()                      */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddBandArg(std::vector<int> *pValue, const char *helpMessage)
{
    AddValidationAction([this]() { return ValidateBandArg(); });

    return AddArg(GDAL_ARG_NAME_BAND, 'b',
                  MsgOrDefault(helpMessage, _("Input band(s) (1-based index)")),
                  pValue)
        .AddValidationAction(
            [pValue]()
            {
                for (int val : *pValue)
                {
                    if (val <= 0)
                    {
                        CPLError(CE_Failure, CPLE_AppDefined,
                                 "Value of 'band' should greater or equal "
                                 "to 1.");
                        return false;
                    }
                }
                return true;
            });
}

/************************************************************************/
/*                      ParseAndValidateKeyValue()                      */
/************************************************************************/

bool GDALAlgorithm::ParseAndValidateKeyValue(GDALAlgorithmArg &arg)
{
    const auto Validate = [this, &arg](const std::string &val)
    {
        if (val.find('=') == std::string::npos)
        {
            ReportError(
                CE_Failure, CPLE_AppDefined,
                "Invalid value for argument '%s'. <KEY>=<VALUE> expected",
                arg.GetName().c_str());
            return false;
        }

        return true;
    };

    if (arg.GetType() == GAAT_STRING)
    {
        return Validate(arg.Get<std::string>());
    }
    else if (arg.GetType() == GAAT_STRING_LIST)
    {
        std::vector<std::string> &vals = arg.Get<std::vector<std::string>>();
        if (vals.size() == 1)
        {
            // Try to split A=B,C=D into A=B and C=D if there is no ambiguity
            std::vector<std::string> newVals;
            std::string curToken;
            bool canSplitOnComma = true;
            char lastSep = 0;
            bool inString = false;
            bool equalFoundInLastToken = false;
            for (char c : vals[0])
            {
                if (!inString && c == ',')
                {
                    if (lastSep != '=' || !equalFoundInLastToken)
                    {
                        canSplitOnComma = false;
                        break;
                    }
                    lastSep = c;
                    newVals.push_back(curToken);
                    curToken.clear();
                    equalFoundInLastToken = false;
                }
                else if (!inString && c == '=')
                {
                    if (lastSep == '=')
                    {
                        canSplitOnComma = false;
                        break;
                    }
                    equalFoundInLastToken = true;
                    lastSep = c;
                    curToken += c;
                }
                else if (c == '"')
                {
                    inString = !inString;
                    curToken += c;
                }
                else
                {
                    curToken += c;
                }
            }
            if (canSplitOnComma && !inString && equalFoundInLastToken)
            {
                if (!curToken.empty())
                    newVals.emplace_back(std::move(curToken));
                vals = std::move(newVals);
            }
        }

        for (const auto &val : vals)
        {
            if (!Validate(val))
                return false;
        }
    }

    return true;
}

/************************************************************************/
/*                           IsGDALGOutput()                            */
/************************************************************************/

bool GDALAlgorithm::IsGDALGOutput() const
{
    bool isGDALGOutput = false;
    const auto outputFormatArg = GetArg(GDAL_ARG_NAME_OUTPUT_FORMAT);
    const auto outputArg = GetArg(GDAL_ARG_NAME_OUTPUT);
    if (outputArg && outputArg->GetType() == GAAT_DATASET &&
        outputArg->IsExplicitlySet())
    {
        if (outputFormatArg && outputFormatArg->GetType() == GAAT_STRING &&
            outputFormatArg->IsExplicitlySet())
        {
            const auto &val =
                outputFormatArg->GDALAlgorithmArg::Get<std::string>();
            isGDALGOutput = EQUAL(val.c_str(), "GDALG");
        }
        else
        {
            const auto &filename =
                outputArg->GDALAlgorithmArg::Get<GDALArgDatasetValue>();
            isGDALGOutput =
                filename.GetName().size() > strlen(".gdalg.json") &&
                EQUAL(filename.GetName().c_str() + filename.GetName().size() -
                          strlen(".gdalg.json"),
                      ".gdalg.json");
        }
    }
    return isGDALGOutput;
}

/************************************************************************/
/*                         ProcessGDALGOutput()                         */
/************************************************************************/

GDALAlgorithm::ProcessGDALGOutputRet GDALAlgorithm::ProcessGDALGOutput()
{
    if (!SupportsStreamedOutput())
        return ProcessGDALGOutputRet::NOT_GDALG;

    if (IsGDALGOutput())
    {
        const auto outputArg = GetArg(GDAL_ARG_NAME_OUTPUT);
        const auto &filename =
            outputArg->GDALAlgorithmArg::Get<GDALArgDatasetValue>().GetName();
        VSIStatBufL sStat;
        if (VSIStatL(filename.c_str(), &sStat) == 0)
        {
            const auto overwriteArg = GetArg(GDAL_ARG_NAME_OVERWRITE);
            if (overwriteArg && overwriteArg->GetType() == GAAT_BOOLEAN)
            {
                if (!overwriteArg->GDALAlgorithmArg::Get<bool>())
                {
                    CPLError(CE_Failure, CPLE_AppDefined,
                             "File '%s' already exists. Specify the "
                             "--overwrite option to overwrite it.",
                             filename.c_str());
                    return ProcessGDALGOutputRet::GDALG_ERROR;
                }
            }
        }

        std::string osCommandLine;

        for (const auto &path : GDALAlgorithm::m_callPath)
        {
            if (!osCommandLine.empty())
                osCommandLine += ' ';
            osCommandLine += path;
        }

        for (const auto &arg : GetArgs())
        {
            if (arg->IsExplicitlySet() &&
                arg->GetName() != GDAL_ARG_NAME_OUTPUT &&
                arg->GetName() != GDAL_ARG_NAME_OUTPUT_FORMAT &&
                arg->GetName() != GDAL_ARG_NAME_UPDATE &&
                arg->GetName() != GDAL_ARG_NAME_OVERWRITE)
            {
                osCommandLine += ' ';
                std::string strArg;
                if (!arg->Serialize(strArg))
                {
                    CPLError(CE_Failure, CPLE_AppDefined,
                             "Cannot serialize argument %s",
                             arg->GetName().c_str());
                    return ProcessGDALGOutputRet::GDALG_ERROR;
                }
                osCommandLine += strArg;
            }
        }

        osCommandLine += " --output-format stream --output streamed_dataset";

        std::string outStringUnused;
        return SaveGDALG(filename, outStringUnused, osCommandLine)
                   ? ProcessGDALGOutputRet::GDALG_OK
                   : ProcessGDALGOutputRet::GDALG_ERROR;
    }

    return ProcessGDALGOutputRet::NOT_GDALG;
}

/************************************************************************/
/*                      GDALAlgorithm::SaveGDALG()                      */
/************************************************************************/

/* static */ bool GDALAlgorithm::SaveGDALG(const std::string &filename,
                                           std::string &outString,
                                           const std::string &commandLine)
{
    CPLJSONDocument oDoc;
    oDoc.GetRoot().Add("type", "gdal_streamed_alg");
    oDoc.GetRoot().Add("command_line", commandLine);
    oDoc.GetRoot().Add("gdal_version", GDALVersionInfo("VERSION_NUM"));

    if (!filename.empty())
        return oDoc.Save(filename);

    outString = oDoc.GetRoot().Format(CPLJSONObject::PrettyFormat::Pretty);
    return true;
}

/************************************************************************/
/*                GDALAlgorithm::AddCreationOptionsArg()                */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddCreationOptionsArg(std::vector<std::string> *pValue,
                                     const char *helpMessage)
{
    auto &arg = AddArg(GDAL_ARG_NAME_CREATION_OPTION, 0,
                       MsgOrDefault(helpMessage, _("Creation option")), pValue)
                    .AddAlias("co")
                    .SetMetaVar("<KEY>=<VALUE>")
                    .SetPackedValuesAllowed(false);
    arg.AddValidationAction([this, &arg]()
                            { return ParseAndValidateKeyValue(arg); });

    arg.SetAutoCompleteFunction(
        [this](const std::string &currentValue)
        {
            std::vector<std::string> oRet;

            int datasetType =
                GDAL_OF_RASTER | GDAL_OF_VECTOR | GDAL_OF_MULTIDIM_RASTER;
            auto outputArg = GetArg(GDAL_ARG_NAME_OUTPUT);
            if (outputArg && (outputArg->GetType() == GAAT_DATASET ||
                              outputArg->GetType() == GAAT_DATASET_LIST))
            {
                datasetType = outputArg->GetDatasetType();
            }

            auto outputFormat = GetArg(GDAL_ARG_NAME_OUTPUT_FORMAT);
            if (outputFormat && outputFormat->GetType() == GAAT_STRING &&
                outputFormat->IsExplicitlySet())
            {
                auto poDriver = GetGDALDriverManager()->GetDriverByName(
                    outputFormat->Get<std::string>().c_str());
                if (poDriver)
                {
                    AddOptionsSuggestions(
                        poDriver->GetMetadataItem(GDAL_DMD_CREATIONOPTIONLIST),
                        datasetType, currentValue, oRet);
                }
                return oRet;
            }

            if (outputArg && outputArg->GetType() == GAAT_DATASET)
            {
                auto poDM = GetGDALDriverManager();
                auto &datasetValue = outputArg->Get<GDALArgDatasetValue>();
                const auto &osDSName = datasetValue.GetName();
                const std::string osExt = CPLGetExtensionSafe(osDSName.c_str());
                if (!osExt.empty())
                {
                    std::set<std::string> oVisitedExtensions;
                    for (int i = 0; i < poDM->GetDriverCount(); ++i)
                    {
                        auto poDriver = poDM->GetDriver(i);
                        if (((datasetType & GDAL_OF_RASTER) != 0 &&
                             poDriver->GetMetadataItem(GDAL_DCAP_RASTER)) ||
                            ((datasetType & GDAL_OF_VECTOR) != 0 &&
                             poDriver->GetMetadataItem(GDAL_DCAP_VECTOR)) ||
                            ((datasetType & GDAL_OF_MULTIDIM_RASTER) != 0 &&
                             poDriver->GetMetadataItem(
                                 GDAL_DCAP_MULTIDIM_RASTER)))
                        {
                            const char *pszExtensions =
                                poDriver->GetMetadataItem(GDAL_DMD_EXTENSIONS);
                            if (pszExtensions)
                            {
                                const CPLStringList aosExts(
                                    CSLTokenizeString2(pszExtensions, " ", 0));
                                for (const char *pszExt : cpl::Iterate(aosExts))
                                {
                                    if (EQUAL(pszExt, osExt.c_str()) &&
                                        !cpl::contains(oVisitedExtensions,
                                                       pszExt))
                                    {
                                        oVisitedExtensions.insert(pszExt);
                                        if (AddOptionsSuggestions(
                                                poDriver->GetMetadataItem(
                                                    GDAL_DMD_CREATIONOPTIONLIST),
                                                datasetType, currentValue,
                                                oRet))
                                        {
                                            return oRet;
                                        }
                                        break;
                                    }
                                }
                            }
                        }
                    }
                }
            }

            return oRet;
        });

    return arg;
}

/************************************************************************/
/*             GDALAlgorithm::AddLayerCreationOptionsArg()              */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddLayerCreationOptionsArg(std::vector<std::string> *pValue,
                                          const char *helpMessage)
{
    auto &arg =
        AddArg(GDAL_ARG_NAME_LAYER_CREATION_OPTION, 0,
               MsgOrDefault(helpMessage, _("Layer creation option")), pValue)
            .AddAlias("lco")
            .SetMetaVar("<KEY>=<VALUE>")
            .SetPackedValuesAllowed(false);
    arg.AddValidationAction([this, &arg]()
                            { return ParseAndValidateKeyValue(arg); });

    arg.SetAutoCompleteFunction(
        [this](const std::string &currentValue)
        {
            std::vector<std::string> oRet;

            auto outputFormat = GetArg(GDAL_ARG_NAME_OUTPUT_FORMAT);
            if (outputFormat && outputFormat->GetType() == GAAT_STRING &&
                outputFormat->IsExplicitlySet())
            {
                auto poDriver = GetGDALDriverManager()->GetDriverByName(
                    outputFormat->Get<std::string>().c_str());
                if (poDriver)
                {
                    AddOptionsSuggestions(poDriver->GetMetadataItem(
                                              GDAL_DS_LAYER_CREATIONOPTIONLIST),
                                          GDAL_OF_VECTOR, currentValue, oRet);
                }
                return oRet;
            }

            auto outputArg = GetArg(GDAL_ARG_NAME_OUTPUT);
            if (outputArg && outputArg->GetType() == GAAT_DATASET)
            {
                auto poDM = GetGDALDriverManager();
                auto &datasetValue = outputArg->Get<GDALArgDatasetValue>();
                const auto &osDSName = datasetValue.GetName();
                const std::string osExt = CPLGetExtensionSafe(osDSName.c_str());
                if (!osExt.empty())
                {
                    std::set<std::string> oVisitedExtensions;
                    for (int i = 0; i < poDM->GetDriverCount(); ++i)
                    {
                        auto poDriver = poDM->GetDriver(i);
                        if (poDriver->GetMetadataItem(GDAL_DCAP_VECTOR))
                        {
                            const char *pszExtensions =
                                poDriver->GetMetadataItem(GDAL_DMD_EXTENSIONS);
                            if (pszExtensions)
                            {
                                const CPLStringList aosExts(
                                    CSLTokenizeString2(pszExtensions, " ", 0));
                                for (const char *pszExt : cpl::Iterate(aosExts))
                                {
                                    if (EQUAL(pszExt, osExt.c_str()) &&
                                        !cpl::contains(oVisitedExtensions,
                                                       pszExt))
                                    {
                                        oVisitedExtensions.insert(pszExt);
                                        if (AddOptionsSuggestions(
                                                poDriver->GetMetadataItem(
                                                    GDAL_DS_LAYER_CREATIONOPTIONLIST),
                                                GDAL_OF_VECTOR, currentValue,
                                                oRet))
                                        {
                                            return oRet;
                                        }
                                        break;
                                    }
                                }
                            }
                        }
                    }
                }
            }

            return oRet;
        });

    return arg;
}

/************************************************************************/
/*                     GDALAlgorithm::AddBBOXArg()                      */
/************************************************************************/

/** Add bbox=xmin,ymin,xmax,ymax argument. */
GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddBBOXArg(std::vector<double> *pValue, const char *helpMessage)
{
    auto &arg = AddArg("bbox", 0,
                       MsgOrDefault(helpMessage,
                                    _("Bounding box as xmin,ymin,xmax,ymax")),
                       pValue)
                    .SetRepeatedArgAllowed(false)
                    .SetMinCount(4)
                    .SetMaxCount(4)
                    .SetDisplayHintAboutRepetition(false);
    arg.AddValidationAction(
        [&arg]()
        {
            const auto &val = arg.Get<std::vector<double>>();
            CPLAssert(val.size() == 4);
            if (!(val[0] <= val[2]) || !(val[1] <= val[3]))
            {
                CPLError(CE_Failure, CPLE_AppDefined,
                         "Value of 'bbox' should be xmin,ymin,xmax,ymax with "
                         "xmin <= xmax and ymin <= ymax");
                return false;
            }
            return true;
        });
    return arg;
}

/************************************************************************/
/*                  GDALAlgorithm::AddActiveLayerArg()                  */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddActiveLayerArg(std::string *pValue, const char *helpMessage)
{
    return AddArg("active-layer", 0,
                  MsgOrDefault(helpMessage,
                               _("Set active layer (if not specified, all)")),
                  pValue);
}

/************************************************************************/
/*                  GDALAlgorithm::AddNumThreadsArg()                   */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddNumThreadsArg(int *pValue, std::string *pStrValue,
                                const char *helpMessage)
{
    auto &arg =
        AddArg(GDAL_ARG_NAME_NUM_THREADS, 'j',
               MsgOrDefault(helpMessage, _("Number of jobs (or ALL_CPUS)")),
               pStrValue);

    AddArg(GDAL_ARG_NAME_NUM_THREADS_INT_HIDDEN, 0,
           _("Number of jobs (read-only, hidden argument)"), pValue)
        .SetHidden();

    auto lambda = [this, &arg, pValue, pStrValue]
    {
        bool bOK = false;
        const char *pszVal = CPLGetConfigOption("GDAL_NUM_THREADS", nullptr);
        const int nLimit = std::clamp(
            pszVal && !EQUAL(pszVal, "ALL_CPUS") ? atoi(pszVal) : INT_MAX, 1,
            CPLGetNumCPUs());
        const int nNumThreads =
            GDALGetNumThreads(pStrValue->c_str(), nLimit,
                              /* bDefaultToAllCPUs = */ false, nullptr, &bOK);
        if (bOK)
        {
            *pValue = nNumThreads;
        }
        else
        {
            ReportError(CE_Failure, CPLE_IllegalArg,
                        "Invalid value for '%s' argument",
                        arg.GetName().c_str());
        }
        return bOK;
    };
    if (!pStrValue->empty())
    {
        arg.SetDefault(*pStrValue);
        lambda();
    }
    arg.AddValidationAction(std::move(lambda));
    return arg;
}

/************************************************************************/
/*                 GDALAlgorithm::AddAbsolutePathArg()                  */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddAbsolutePathArg(bool *pValue, const char *helpMessage)
{
    return AddArg(
        "absolute-path", 0,
        MsgOrDefault(helpMessage, _("Whether the path to the input dataset "
                                    "should be stored as an absolute path")),
        pValue);
}

/************************************************************************/
/*               GDALAlgorithm::AddPixelFunctionNameArg()               */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddPixelFunctionNameArg(std::string *pValue,
                                       const char *helpMessage)
{

    const auto pixelFunctionNames =
        VRTDerivedRasterBand::GetPixelFunctionNames();
    return AddArg(
               "pixel-function", 0,
               MsgOrDefault(
                   helpMessage,
                   _("Specify a pixel function to calculate output value from "
                     "overlapping inputs")),
               pValue)
        .SetChoices(pixelFunctionNames);
}

/************************************************************************/
/*               GDALAlgorithm::AddPixelFunctionArgsArg()               */
/************************************************************************/

GDALInConstructionAlgorithmArg &
GDALAlgorithm::AddPixelFunctionArgsArg(std::vector<std::string> *pValue,
                                       const char *helpMessage)
{
    auto &pixelFunctionArgArg =
        AddArg("pixel-function-arg", 0,
               MsgOrDefault(
                   helpMessage,
                   _("Specify argument(s) to pass to the pixel function")),
               pValue)
            .SetMetaVar("<NAME>=<VALUE>")
            .SetRepeatedArgAllowed(true);
    pixelFunctionArgArg.AddValidationAction(
        [this, &pixelFunctionArgArg]()
        { return ParseAndValidateKeyValue(pixelFunctionArgArg); });

    pixelFunctionArgArg.SetAutoCompleteFunction(
        [this](const std::string &currentValue)
        {
            std::string pixelFunction;
            const auto pixelFunctionArg = GetArg("pixel-function");
            if (pixelFunctionArg && pixelFunctionArg->GetType() == GAAT_STRING)
            {
                pixelFunction = pixelFunctionArg->Get<std::string>();
            }

            std::vector<std::string> ret;

            if (!pixelFunction.empty())
            {
                const auto *pair = VRTDerivedRasterBand::GetPixelFunction(
                    pixelFunction.c_str());
                if (!pair)
                {
                    ret.push_back("**");
                    // Non printable UTF-8 space, to avoid autocompletion to pickup on 'd'
                    ret.push_back(std::string("\xC2\xA0"
                                              "Invalid pixel function name"));
                }
                else if (pair->second.find("Argument name=") ==
                         std::string::npos)
                {
                    ret.push_back("**");
                    // Non printable UTF-8 space, to avoid autocompletion to pickup on 'd'
                    ret.push_back(
                        std::string(
                            "\xC2\xA0"
                            "No pixel function arguments for pixel function '")
                            .append(pixelFunction)
                            .append("'"));
                }
                else
                {
                    AddOptionsSuggestions(pair->second.c_str(), 0, currentValue,
                                          ret);
                }
            }

            return ret;
        });

    return pixelFunctionArgArg;
}

/************************************************************************/
/*                   GDALAlgorithm::AddProgressArg()                    */
/************************************************************************/

void GDALAlgorithm::AddProgressArg()
{
    AddArg(GDAL_ARG_NAME_QUIET, 'q',
           _("Quiet mode (no progress bar or warning message)"), &m_quiet)
        .SetAvailableInPipelineStep(false)
        .SetCategory(GAAC_COMMON)
        .AddAction([this]() { m_progressBarRequested = false; });

    AddArg("progress", 0, _("Display progress bar"), &m_progressBarRequested)
        .SetAvailableInPipelineStep(false)
        .SetHidden();
}

/************************************************************************/
/*                         GDALAlgorithm::Run()                         */
/************************************************************************/

bool GDALAlgorithm::Run(GDALProgressFunc pfnProgress, void *pProgressData)
{
    WarnIfDeprecated();

    if (m_selectedSubAlg)
    {
        if (m_calledFromCommandLine)
            m_selectedSubAlg->m_calledFromCommandLine = true;
        return m_selectedSubAlg->Run(pfnProgress, pProgressData);
    }

    if (m_helpRequested || m_helpDocRequested)
    {
        if (m_calledFromCommandLine)
            printf("%s", GetUsageForCLI(false).c_str()); /*ok*/
        return true;
    }

    if (m_JSONUsageRequested)
    {
        if (m_calledFromCommandLine)
            printf("%s", GetUsageAsJSON().c_str()); /*ok*/
        return true;
    }

    if (!ValidateArguments())
        return false;

    if (m_alreadyRun)
    {
        ReportError(CE_Failure, CPLE_AppDefined,
                    "Run() can be called only once per algorithm instance");
        return false;
    }
    m_alreadyRun = true;

    switch (ProcessGDALGOutput())
    {
        case ProcessGDALGOutputRet::GDALG_ERROR:
            return false;

        case ProcessGDALGOutputRet::GDALG_OK:
            return true;

        case ProcessGDALGOutputRet::NOT_GDALG:
            break;
    }

    if (m_executionForStreamOutput)
    {
        if (!CheckSafeForStreamOutput())
        {
            return false;
        }
    }

    return RunImpl(pfnProgress, pProgressData);
}

/************************************************************************/
/*              GDALAlgorithm::CheckSafeForStreamOutput()               */
/************************************************************************/

bool GDALAlgorithm::CheckSafeForStreamOutput()
{
    const auto outputFormatArg = GetArg(GDAL_ARG_NAME_OUTPUT_FORMAT);
    if (outputFormatArg && outputFormatArg->GetType() == GAAT_STRING)
    {
        const auto &val = outputFormatArg->GDALAlgorithmArg::Get<std::string>();
        if (!EQUAL(val.c_str(), "stream"))
        {
            // For security reasons, to avoid that reading a .gdalg.json file
            // writes a file on the file system.
            ReportError(
                CE_Failure, CPLE_NotSupported,
                "in streamed execution, --format stream should be used");
            return false;
        }
    }
    return true;
}

/************************************************************************/
/*                      GDALAlgorithm::Finalize()                       */
/************************************************************************/

bool GDALAlgorithm::Finalize()
{
    bool ret = true;
    if (m_selectedSubAlg)
        ret = m_selectedSubAlg->Finalize();

    for (auto &arg : m_args)
    {
        if (arg->GetType() == GAAT_DATASET)
        {
            ret = arg->Get<GDALArgDatasetValue>().Close() && ret;
        }
        else if (arg->GetType() == GAAT_DATASET_LIST)
        {
            for (auto &ds : arg->Get<std::vector<GDALArgDatasetValue>>())
            {
                ret = ds.Close() && ret;
            }
        }
    }
    return ret;
}

/************************************************************************/
/*                  GDALAlgorithm::GetArgNamesForCLI()                  */
/************************************************************************/

std::pair<std::vector<std::pair<GDALAlgorithmArg *, std::string>>, size_t>
GDALAlgorithm::GetArgNamesForCLI() const
{
    std::vector<std::pair<GDALAlgorithmArg *, std::string>> options;

    size_t maxOptLen = 0;
    for (const auto &arg : m_args)
    {
        if (arg->IsHidden() || arg->IsHiddenForCLI())
            continue;
        std::string opt;
        bool addComma = false;
        if (!arg->GetShortName().empty())
        {
            opt += '-';
            opt += arg->GetShortName();
            addComma = true;
        }
        for (char alias : arg->GetShortNameAliases())
        {
            if (addComma)
                opt += ", ";
            opt += "-";
            opt += alias;
            addComma = true;
        }
        for (const std::string &alias : arg->GetAliases())
        {
            if (addComma)
                opt += ", ";
            opt += "--";
            opt += alias;
            addComma = true;
        }
        if (!arg->GetName().empty())
        {
            if (addComma)
                opt += ", ";
            opt += "--";
            opt += arg->GetName();
        }
        const auto &metaVar = arg->GetMetaVar();
        if (!metaVar.empty())
        {
            opt += ' ';
            if (metaVar.front() != '<')
                opt += '<';
            opt += metaVar;
            if (metaVar.back() != '>')
                opt += '>';
        }
        maxOptLen = std::max(maxOptLen, opt.size());
        options.emplace_back(arg.get(), opt);
    }

    return std::make_pair(std::move(options), maxOptLen);
}

/************************************************************************/
/*                   GDALAlgorithm::GetUsageForCLI()                    */
/************************************************************************/

std::string
GDALAlgorithm::GetUsageForCLI(bool shortUsage,
                              const UsageOptions &usageOptions) const
{
    if (m_selectedSubAlg)
        return m_selectedSubAlg->GetUsageForCLI(shortUsage, usageOptions);

    std::string osRet(usageOptions.isPipelineStep ? "*" : "Usage:");
    std::string osPath;
    for (const std::string &s : m_callPath)
    {
        if (!osPath.empty())
            osPath += ' ';
        osPath += s;
    }
    osRet += ' ';
    osRet += osPath;

    bool hasNonPositionals = false;
    for (const auto &arg : m_args)
    {
        if (!arg->IsHidden() && !arg->IsHiddenForCLI() && !arg->IsPositional())
            hasNonPositionals = true;
    }

    if (HasSubAlgorithms())
    {
        if (m_callPath.size() == 1)
        {
            osRet += " <COMMAND>";
            if (hasNonPositionals)
                osRet += " [OPTIONS]";
            if (usageOptions.isPipelineStep)
            {
                const size_t nLenFirstLine = osRet.size();
                osRet += '\n';
                osRet.append(nLenFirstLine, '-');
                osRet += '\n';
            }
            osRet += "\nwhere <COMMAND> is one of:\n";
        }
        else
        {
            osRet += " <SUBCOMMAND>";
            if (hasNonPositionals)
                osRet += " [OPTIONS]";
            if (usageOptions.isPipelineStep)
            {
                const size_t nLenFirstLine = osRet.size();
                osRet += '\n';
                osRet.append(nLenFirstLine, '-');
                osRet += '\n';
            }
            osRet += "\nwhere <SUBCOMMAND> is one of:\n";
        }
        size_t maxNameLen = 0;
        for (const auto &subAlgName : GetSubAlgorithmNames())
        {
            maxNameLen = std::max(maxNameLen, subAlgName.size());
        }
        for (const auto &subAlgName : GetSubAlgorithmNames())
        {
            auto subAlg = InstantiateSubAlgorithm(subAlgName);
            if (subAlg && !subAlg->IsHidden())
            {
                const std::string &name(subAlg->GetName());
                osRet += "  - ";
                osRet += name;
                osRet += ": ";
                osRet.append(maxNameLen - name.size(), ' ');
                osRet += subAlg->GetDescription();
                if (!subAlg->m_aliases.empty())
                {
                    bool first = true;
                    for (const auto &alias : subAlg->GetAliases())
                    {
                        if (alias ==
                            GDALAlgorithmRegistry::HIDDEN_ALIAS_SEPARATOR)
                            break;
                        if (first)
                            osRet += " (alias: ";
                        else
                            osRet += ", ";
                        osRet += alias;
                        first = false;
                    }
                    if (!first)
                    {
                        osRet += ')';
                    }
                }
                osRet += '\n';
            }
        }

        if (shortUsage && hasNonPositionals)
        {
            osRet += "\nTry '";
            osRet += osPath;
            osRet += " --help' for help.\n";
        }
    }
    else
    {
        if (!m_args.empty())
        {
            if (hasNonPositionals)
                osRet += " [OPTIONS]";
            for (const auto *arg : m_positionalArgs)
            {
                if ((!arg->IsHidden() && !arg->IsHiddenForCLI()) ||
                    (GetName() == "pipeline" && arg->GetName() == "pipeline"))
                {
                    const bool optional =
                        (!arg->IsRequired() && !(GetName() == "pipeline" &&
                                                 arg->GetName() == "pipeline"));
                    osRet += ' ';
                    if (optional)
                        osRet += '[';
                    const std::string &metavar = arg->GetMetaVar();
                    if (!metavar.empty() && metavar[0] == '<')
                    {
                        osRet += metavar;
                    }
                    else
                    {
                        osRet += '<';
                        osRet += metavar;
                        osRet += '>';
                    }
                    if (arg->GetType() == GAAT_DATASET_LIST &&
                        arg->GetMaxCount() > 1)
                    {
                        osRet += "...";
                    }
                    if (optional)
                        osRet += ']';
                }
            }
        }

        const size_t nLenFirstLine = osRet.size();
        osRet += '\n';
        if (usageOptions.isPipelineStep)
        {
            osRet.append(nLenFirstLine, '-');
            osRet += '\n';
        }

        if (shortUsage)
        {
            osRet += "Try '";
            osRet += osPath;
            osRet += " --help' for help.\n";
            return osRet;
        }

        osRet += '\n';
        osRet += m_description;
        osRet += '\n';
    }

    if (!m_args.empty() && !shortUsage)
    {
        std::vector<std::pair<GDALAlgorithmArg *, std::string>> options;
        size_t maxOptLen;
        std::tie(options, maxOptLen) = GetArgNamesForCLI();
        if (usageOptions.maxOptLen)
            maxOptLen = usageOptions.maxOptLen;

        const std::string userProvidedOpt = "--<user-provided-option>=<value>";
        if (m_arbitraryLongNameArgsAllowed)
            maxOptLen = std::max(maxOptLen, userProvidedOpt.size());

        const auto OutputArg =
            [this, maxOptLen, &osRet,
             &usageOptions](const GDALAlgorithmArg *arg, const std::string &opt)
        {
            osRet += "  ";
            osRet += opt;
            osRet += "  ";
            osRet.append(maxOptLen - opt.size(), ' ');
            osRet += arg->GetDescription();

            const auto &choices = arg->GetChoices();
            if (!choices.empty())
            {
                osRet += ". ";
                osRet += arg->GetMetaVar();
                osRet += '=';
                bool firstChoice = true;
                for (const auto &choice : choices)
                {
                    if (!firstChoice)
                        osRet += '|';
                    osRet += choice;
                    firstChoice = false;
                }
            }

            if (arg->GetType() == GAAT_DATASET ||
                arg->GetType() == GAAT_DATASET_LIST)
            {
                if (arg->IsOutput() &&
                    arg->GetDatasetInputFlags() == GADV_NAME &&
                    arg->GetDatasetOutputFlags() == GADV_OBJECT)
                {
                    osRet += " (created by algorithm)";
                }
            }

            if (arg->GetType() == GAAT_STRING && arg->HasDefaultValue())
            {
                osRet += " (default: ";
                osRet += arg->GetDefault<std::string>();
                osRet += ')';
            }
            else if (arg->GetType() == GAAT_BOOLEAN && arg->HasDefaultValue())
            {
                if (arg->GetDefault<bool>())
                    osRet += " (default: true)";
            }
            else if (arg->GetType() == GAAT_INTEGER && arg->HasDefaultValue())
            {
                osRet += " (default: ";
                osRet += CPLSPrintf("%d", arg->GetDefault<int>());
                osRet += ')';
            }
            else if (arg->GetType() == GAAT_REAL && arg->HasDefaultValue())
            {
                osRet += " (default: ";
                osRet += CPLSPrintf("%g", arg->GetDefault<double>());
                osRet += ')';
            }
            else if (arg->GetType() == GAAT_STRING_LIST &&
                     arg->HasDefaultValue())
            {
                const auto &defaultVal =
                    arg->GetDefault<std::vector<std::string>>();
                if (defaultVal.size() == 1)
                {
                    osRet += " (default: ";
                    osRet += defaultVal[0];
                    osRet += ')';
                }
            }
            else if (arg->GetType() == GAAT_INTEGER_LIST &&
                     arg->HasDefaultValue())
            {
                const auto &defaultVal = arg->GetDefault<std::vector<int>>();
                if (defaultVal.size() == 1)
                {
                    osRet += " (default: ";
                    osRet += CPLSPrintf("%d", defaultVal[0]);
                    osRet += ')';
                }
            }
            else if (arg->GetType() == GAAT_REAL_LIST && arg->HasDefaultValue())
            {
                const auto &defaultVal = arg->GetDefault<std::vector<double>>();
                if (defaultVal.size() == 1)
                {
                    osRet += " (default: ";
                    osRet += CPLSPrintf("%g", defaultVal[0]);
                    osRet += ')';
                }
            }

            if (arg->GetDisplayHintAboutRepetition())
            {
                if (arg->GetMinCount() > 0 &&
                    arg->GetMinCount() == arg->GetMaxCount())
                {
                    if (arg->GetMinCount() != 1)
                        osRet += CPLSPrintf(" [%d values]", arg->GetMaxCount());
                }
                else if (arg->GetMinCount() > 0 &&
                         arg->GetMaxCount() < GDALAlgorithmArgDecl::UNBOUNDED)
                {
                    osRet += CPLSPrintf(" [%d..%d values]", arg->GetMinCount(),
                                        arg->GetMaxCount());
                }
                else if (arg->GetMinCount() > 0)
                {
                    osRet += CPLSPrintf(" [%d.. values]", arg->GetMinCount());
                }
                else if (arg->GetMaxCount() > 1)
                {
                    osRet += " [may be repeated]";
                }
            }

            if (arg->IsRequired())
            {
                osRet += " [required]";
            }

            if (!arg->IsAvailableInPipelineStep() &&
                !usageOptions.isPipelineStep)
            {
                osRet += " [not available in pipelines]";
            }

            osRet += '\n';

            const auto &mutualExclusionGroup = arg->GetMutualExclusionGroup();
            if (!mutualExclusionGroup.empty())
            {
                std::string otherArgs;
                for (const auto &otherArg : m_args)
                {
                    if (otherArg->IsHidden() || otherArg->IsHiddenForCLI() ||
                        otherArg.get() == arg)
                        continue;
                    if (otherArg->GetMutualExclusionGroup() ==
                        mutualExclusionGroup)
                    {
                        if (!otherArgs.empty())
                            otherArgs += ", ";
                        otherArgs += "--";
                        otherArgs += otherArg->GetName();
                    }
                }
                if (!otherArgs.empty())
                {
                    osRet += "  ";
                    osRet += "  ";
                    osRet.append(maxOptLen, ' ');
                    osRet += "Mutually exclusive with ";
                    osRet += otherArgs;
                    osRet += '\n';
                }
            }

            // Check dependency
            std::string dependencyArgs;

            for (const auto &dependencyArgumentName :
                 GetArgDependencies(arg->GetName()))
            {
                const auto otherArg{GetArg(dependencyArgumentName)};
                if (otherArg != nullptr)
                {
                    if (otherArg->IsHidden() || otherArg->IsHiddenForCLI() ||
                        otherArg == arg)
                    {
                        continue;
                    }

                    if (!dependencyArgs.empty())
                    {
                        dependencyArgs += ", ";
                    }

                    dependencyArgs += "--";
                    dependencyArgs += otherArg->GetName();
                }
                else
                {
                    CPLError(CE_Warning, CPLE_AppDefined,
                             "Argument '%s' depends on unknown argument '%s'",
                             arg->GetName().c_str(),
                             dependencyArgumentName.c_str());
                }
            }

            if (!dependencyArgs.empty())
            {
                osRet += "  ";
                osRet += "  ";
                osRet.append(maxOptLen, ' ');
                osRet += "Depends on ";
                osRet += dependencyArgs;
                osRet += '\n';
            }
        };

        if (!m_positionalArgs.empty())
        {
            osRet += "\nPositional arguments:\n";
            for (const auto &[arg, opt] : options)
            {
                if (arg->IsPositional())
                    OutputArg(arg, opt);
            }
        }

        if (hasNonPositionals)
        {
            bool hasCommon = false;
            bool hasBase = false;
            bool hasAdvanced = false;
            bool hasEsoteric = false;
            std::vector<std::string> categories;
            for (const auto &iter : options)
            {
                const auto &arg = iter.first;
                if (!arg->IsPositional())
                {
                    const auto &category = arg->GetCategory();
                    if (category == GAAC_COMMON)
                    {
                        hasCommon = true;
                    }
                    else if (category == GAAC_BASE)
                    {
                        hasBase = true;
                    }
                    else if (category == GAAC_ADVANCED)
                    {
                        hasAdvanced = true;
                    }
                    else if (category == GAAC_ESOTERIC)
                    {
                        hasEsoteric = true;
                    }
                    else if (std::find(categories.begin(), categories.end(),
                                       category) == categories.end())
                    {
                        categories.push_back(category);
                    }
                }
            }
            if (hasAdvanced || m_arbitraryLongNameArgsAllowed)
                categories.insert(categories.begin(), GAAC_ADVANCED);
            if (hasBase)
                categories.insert(categories.begin(), GAAC_BASE);
            if (hasCommon && !usageOptions.isPipelineStep)
                categories.insert(categories.begin(), GAAC_COMMON);
            if (hasEsoteric)
                categories.push_back(GAAC_ESOTERIC);

            for (const auto &category : categories)
            {
                osRet += "\n";
                if (category != GAAC_BASE)
                {
                    osRet += category;
                    osRet += ' ';
                }
                osRet += "Options:\n";
                for (const auto &[arg, opt] : options)
                {
                    if (!arg->IsPositional() && arg->GetCategory() == category)
                        OutputArg(arg, opt);
                }
                if (m_arbitraryLongNameArgsAllowed && category == GAAC_ADVANCED)
                {
                    osRet += "  ";
                    osRet += userProvidedOpt;
                    osRet += "  ";
                    if (userProvidedOpt.size() < maxOptLen)
                        osRet.append(maxOptLen - userProvidedOpt.size(), ' ');
                    osRet += "Argument provided by user";
                    osRet += '\n';
                }
            }
        }
    }

    if (!m_longDescription.empty())
    {
        osRet += '\n';
        osRet += m_longDescription;
        osRet += '\n';
    }

    if (!m_helpDocRequested && !usageOptions.isPipelineMain)
    {
        if (!m_helpURL.empty())
        {
            osRet += "\nFor more details, consult ";
            osRet += GetHelpFullURL();
            osRet += '\n';
        }
        osRet += GetUsageForCLIEnd();
    }

    return osRet;
}

/************************************************************************/
/*                  GDALAlgorithm::GetUsageForCLIEnd()                  */
/************************************************************************/

//! @cond Doxygen_Suppress
std::string GDALAlgorithm::GetUsageForCLIEnd() const
{
    std::string osRet;

    if (!m_callPath.empty() && m_callPath[0] == "gdal")
    {
        osRet += "\nWARNING: the gdal command is provisionally provided as an "
                 "alternative interface to GDAL and OGR command line "
                 "utilities.\nThe project reserves the right to modify, "
                 "rename, reorganize, and change the behavior of the utility\n"
                 "until it is officially frozen in a future feature release of "
                 "GDAL.\n";
    }
    return osRet;
}

//! @endcond

/************************************************************************/
/*                   GDALAlgorithm::GetUsageAsJSON()                    */
/************************************************************************/

std::string GDALAlgorithm::GetUsageAsJSON() const
{
    CPLJSONDocument oDoc;
    auto oRoot = oDoc.GetRoot();

    if (m_displayInJSONUsage)
    {
        oRoot.Add("name", m_name);
        CPLJSONArray jFullPath;
        for (const std::string &s : m_callPath)
        {
            jFullPath.Add(s);
        }
        oRoot.Add("full_path", jFullPath);
    }

    oRoot.Add("description", m_description);
    if (!m_helpURL.empty())
    {
        oRoot.Add("short_url", m_helpURL);
        oRoot.Add("url", GetHelpFullURL());
    }

    CPLJSONArray jSubAlgorithms;
    for (const auto &subAlgName : GetSubAlgorithmNames())
    {
        auto subAlg = InstantiateSubAlgorithm(subAlgName);
        if (subAlg && subAlg->m_displayInJSONUsage && !subAlg->IsHidden())
        {
            CPLJSONDocument oSubDoc;
            CPL_IGNORE_RET_VAL(oSubDoc.LoadMemory(subAlg->GetUsageAsJSON()));
            jSubAlgorithms.Add(oSubDoc.GetRoot());
        }
    }
    oRoot.Add("sub_algorithms", jSubAlgorithms);

    if (m_arbitraryLongNameArgsAllowed)
    {
        oRoot.Add("user_provided_arguments_allowed", true);
    }

    const auto ProcessArg = [this](const GDALAlgorithmArg *arg)
    {
        CPLJSONObject jArg;
        jArg.Add("name", arg->GetName());
        jArg.Add("type", GDALAlgorithmArgTypeName(arg->GetType()));
        jArg.Add("description", arg->GetDescription());

        const auto &metaVar = arg->GetMetaVar();
        if (!metaVar.empty() && metaVar != CPLString(arg->GetName()).toupper())
        {
            if (metaVar.front() == '<' && metaVar.back() == '>' &&
                metaVar.substr(1, metaVar.size() - 2).find('>') ==
                    std::string::npos)
                jArg.Add("metavar", metaVar.substr(1, metaVar.size() - 2));
            else
                jArg.Add("metavar", metaVar);
        }

        if (!arg->IsAvailableInPipelineStep())
        {
            jArg.Add("available_in_pipeline_step", false);
        }

        const auto &choices = arg->GetChoices();
        if (!choices.empty())
        {
            CPLJSONArray jChoices;
            for (const auto &choice : choices)
                jChoices.Add(choice);
            jArg.Add("choices", jChoices);
        }
        if (arg->HasDefaultValue())
        {
            switch (arg->GetType())
            {
                case GAAT_BOOLEAN:
                    jArg.Add("default", arg->GetDefault<bool>());
                    break;
                case GAAT_STRING:
                    jArg.Add("default", arg->GetDefault<std::string>());
                    break;
                case GAAT_INTEGER:
                    jArg.Add("default", arg->GetDefault<int>());
                    break;
                case GAAT_REAL:
                    jArg.Add("default", arg->GetDefault<double>());
                    break;
                case GAAT_STRING_LIST:
                {
                    const auto &val =
                        arg->GetDefault<std::vector<std::string>>();
                    if (val.size() == 1)
                    {
                        jArg.Add("default", val[0]);
                    }
                    else
                    {
                        CPLJSONArray jArr;
                        for (const auto &s : val)
                        {
                            jArr.Add(s);
                        }
                        jArg.Add("default", jArr);
                    }
                    break;
                }
                case GAAT_INTEGER_LIST:
                {
                    const auto &val = arg->GetDefault<std::vector<int>>();
                    if (val.size() == 1)
                    {
                        jArg.Add("default", val[0]);
                    }
                    else
                    {
                        CPLJSONArray jArr;
                        for (int i : val)
                        {
                            jArr.Add(i);
                        }
                        jArg.Add("default", jArr);
                    }
                    break;
                }
                case GAAT_REAL_LIST:
                {
                    const auto &val = arg->GetDefault<std::vector<double>>();
                    if (val.size() == 1)
                    {
                        jArg.Add("default", val[0]);
                    }
                    else
                    {
                        CPLJSONArray jArr;
                        for (double d : val)
                        {
                            jArr.Add(d);
                        }
                        jArg.Add("default", jArr);
                    }
                    break;
                }
                case GAAT_DATASET:
                case GAAT_DATASET_LIST:
                    CPLError(CE_Warning, CPLE_AppDefined,
                             "Unhandled default value for arg %s",
                             arg->GetName().c_str());
                    break;
            }
        }

        const auto [minVal, minValIsIncluded] = arg->GetMinValue();
        if (!std::isnan(minVal))
        {
            if (arg->GetType() == GAAT_INTEGER ||
                arg->GetType() == GAAT_INTEGER_LIST)
                jArg.Add("min_value", static_cast<int>(minVal));
            else
                jArg.Add("min_value", minVal);
            jArg.Add("min_value_is_included", minValIsIncluded);
        }

        const auto [maxVal, maxValIsIncluded] = arg->GetMaxValue();
        if (!std::isnan(maxVal))
        {
            if (arg->GetType() == GAAT_INTEGER ||
                arg->GetType() == GAAT_INTEGER_LIST)
                jArg.Add("max_value", static_cast<int>(maxVal));
            else
                jArg.Add("max_value", maxVal);
            jArg.Add("max_value_is_included", maxValIsIncluded);
        }

        jArg.Add("required", arg->IsRequired());
        if (GDALAlgorithmArgTypeIsList(arg->GetType()))
        {
            jArg.Add("packed_values_allowed", arg->GetPackedValuesAllowed());
            jArg.Add("repeated_arg_allowed", arg->GetRepeatedArgAllowed());
            jArg.Add("min_count", arg->GetMinCount());
            jArg.Add("max_count", arg->GetMaxCount());
        }

        // Process dependencies
        const auto &mutualDependencyGroup = arg->GetMutualDependencyGroup();
        if (!mutualDependencyGroup.empty())
        {
            jArg.Add("mutual_dependency_group", mutualDependencyGroup);
        }

        CPLJSONArray jDependencies;
        for (const auto &dependencyArgumentName :
             GetArgDependencies(arg->GetName()))
        {
            jDependencies.Add(dependencyArgumentName);
        }

        if (jDependencies.Size() > 0)
        {
            jArg.Add("depends_on", jDependencies);
        }

        jArg.Add("category", arg->GetCategory());

        if (arg->GetType() == GAAT_DATASET ||
            arg->GetType() == GAAT_DATASET_LIST)
        {
            {
                CPLJSONArray jAr;
                if (arg->GetDatasetType() & GDAL_OF_RASTER)
                    jAr.Add("raster");
                if (arg->GetDatasetType() & GDAL_OF_VECTOR)
                    jAr.Add("vector");
                if (arg->GetDatasetType() & GDAL_OF_MULTIDIM_RASTER)
                    jAr.Add("multidim_raster");
                jArg.Add("dataset_type", jAr);
            }

            const auto GetFlags = [](int flags)
            {
                CPLJSONArray jAr;
                if (flags & GADV_NAME)
                    jAr.Add("name");
                if (flags & GADV_OBJECT)
                    jAr.Add("dataset");
                return jAr;
            };

            if (arg->IsInput())
            {
                jArg.Add("input_flags", GetFlags(arg->GetDatasetInputFlags()));
            }
            if (arg->IsOutput())
            {
                jArg.Add("output_flags",
                         GetFlags(arg->GetDatasetOutputFlags()));
            }
        }

        const auto &mutualExclusionGroup = arg->GetMutualExclusionGroup();
        if (!mutualExclusionGroup.empty())
        {
            jArg.Add("mutual_exclusion_group", mutualExclusionGroup);
        }

        const auto &metadata = arg->GetMetadata();
        if (!metadata.empty())
        {
            CPLJSONObject jMetadata;
            for (const auto &[key, values] : metadata)
            {
                CPLJSONArray jValue;
                for (const auto &value : values)
                    jValue.Add(value);
                jMetadata.Add(key, jValue);
            }
            jArg.Add("metadata", jMetadata);
        }

        return jArg;
    };

    {
        CPLJSONArray jArgs;
        for (const auto &arg : m_args)
        {
            if (!arg->IsHiddenForAPI() && arg->IsInput() && !arg->IsOutput())
                jArgs.Add(ProcessArg(arg.get()));
        }
        oRoot.Add("input_arguments", jArgs);
    }

    {
        CPLJSONArray jArgs;
        for (const auto &arg : m_args)
        {
            if (!arg->IsHiddenForAPI() && !arg->IsInput() && arg->IsOutput())
                jArgs.Add(ProcessArg(arg.get()));
        }
        oRoot.Add("output_arguments", jArgs);
    }

    {
        CPLJSONArray jArgs;
        for (const auto &arg : m_args)
        {
            if (!arg->IsHiddenForAPI() && arg->IsInput() && arg->IsOutput())
                jArgs.Add(ProcessArg(arg.get()));
        }
        oRoot.Add("input_output_arguments", jArgs);
    }

    if (m_supportsStreamedOutput)
    {
        oRoot.Add("supports_streamed_output", true);
    }

    return oDoc.SaveAsString();
}

/************************************************************************/
/*                   GDALAlgorithm::GetAutoComplete()                   */
/************************************************************************/

std::vector<std::string>
GDALAlgorithm::GetAutoComplete(std::vector<std::string> &args,
                               bool lastWordIsComplete, bool showAllOptions)
{
    std::vector<std::string> ret;

    // Get inner-most algorithm
    std::unique_ptr<GDALAlgorithm> curAlgHolder;
    GDALAlgorithm *curAlg = this;
    while (!args.empty() && !args.front().empty() && args.front()[0] != '-')
    {
        auto subAlg = curAlg->InstantiateSubAlgorithm(
            args.front(), /* suggestionAllowed = */ false);
        if (!subAlg)
            break;
        if (args.size() == 1 && !lastWordIsComplete)
        {
            int nCount = 0;
            for (const auto &subAlgName : curAlg->GetSubAlgorithmNames())
            {
                if (STARTS_WITH(subAlgName.c_str(), args.front().c_str()))
                    nCount++;
            }
            if (nCount >= 2)
            {
                for (const std::string &subAlgName :
                     curAlg->GetSubAlgorithmNames())
                {
                    subAlg = curAlg->InstantiateSubAlgorithm(subAlgName);
                    if (subAlg && !subAlg->IsHidden())
                        ret.push_back(subAlg->GetName());
                }
                return ret;
            }
        }
        showAllOptions = false;
        args.erase(args.begin());
        curAlgHolder = std::move(subAlg);
        curAlg = curAlgHolder.get();
    }
    if (curAlg != this)
    {
        curAlg->m_calledFromCommandLine = m_calledFromCommandLine;
        return curAlg->GetAutoComplete(args, lastWordIsComplete,
                                       /* showAllOptions = */ false);
    }

    std::string option;
    std::string value;
    ExtractLastOptionAndValue(args, option, value);

    if (option.empty() && !args.empty() && !args.back().empty() &&
        args.back()[0] == '-')
    {
        const auto &lastArg = args.back();
        // List available options
        for (const auto &arg : GetArgs())
        {
            if (arg->IsHidden() || arg->IsHiddenForCLI() ||
                (!showAllOptions &&
                 (arg->GetName() == "help" || arg->GetName() == "config" ||
                  arg->GetName() == "version" ||
                  arg->GetName() == "json-usage")))
            {
                continue;
            }
            if (!arg->GetShortName().empty())
            {
                std::string str = std::string("-").append(arg->GetShortName());
                if (lastArg == str)
                    ret.push_back(std::move(str));
            }
            if (lastArg != "-" && lastArg != "--")
            {
                for (const std::string &alias : arg->GetAliases())
                {
                    std::string str = std::string("--").append(alias);
                    if (cpl::starts_with(str, lastArg))
                        ret.push_back(std::move(str));
                }
            }
            if (!arg->GetName().empty())
            {
                std::string str = std::string("--").append(arg->GetName());
                if (cpl::starts_with(str, lastArg))
                    ret.push_back(std::move(str));
            }
        }
        std::sort(ret.begin(), ret.end());
    }
    else if (!option.empty())
    {
        // List possible choices for current option
        auto arg = GetArg(option);
        if (arg && arg->GetType() != GAAT_BOOLEAN)
        {
            ret = arg->GetChoices();
            if (ret.empty())
            {
                {
                    CPLErrorStateBackuper oErrorQuieter(CPLQuietErrorHandler);
                    SetParseForAutoCompletion();
                    CPL_IGNORE_RET_VAL(ParseCommandLineArguments(args));
                }
                ret = arg->GetAutoCompleteChoices(value);
            }
            else
            {
                std::sort(ret.begin(), ret.end());
            }
            if (!ret.empty() && ret.back() == value)
            {
                ret.clear();
            }
            else if (ret.empty())
            {
                ret.push_back("**");
                // Non printable UTF-8 space, to avoid autocompletion to pickup on 'd'
                ret.push_back(std::string("\xC2\xA0"
                                          "description: ")
                                  .append(arg->GetDescription()));
            }
        }
    }
    else
    {
        // List possible sub-algorithms
        for (const std::string &subAlgName : GetSubAlgorithmNames())
        {
            auto subAlg = InstantiateSubAlgorithm(subAlgName);
            if (subAlg && !subAlg->IsHidden())
                ret.push_back(subAlg->GetName());
        }
        if (!ret.empty())
        {
            std::sort(ret.begin(), ret.end());
        }

        // Try filenames
        if (ret.empty() && !args.empty())
        {
            {
                CPLErrorStateBackuper oErrorQuieter(CPLQuietErrorHandler);
                SetParseForAutoCompletion();
                CPL_IGNORE_RET_VAL(ParseCommandLineArguments(args));
            }

            const std::string &lastArg = args.back();
            GDALAlgorithmArg *arg = nullptr;
            for (const char *name : {GDAL_ARG_NAME_INPUT, "dataset", "filename",
                                     "like", "source", "destination"})
            {
                if (!arg)
                {
                    auto newArg = GetArg(name);
                    if (newArg)
                    {
                        if (!newArg->IsExplicitlySet())
                        {
                            arg = newArg;
                        }
                        else if (newArg->GetType() == GAAT_STRING ||
                                 newArg->GetType() == GAAT_STRING_LIST ||
                                 newArg->GetType() == GAAT_DATASET ||
                                 newArg->GetType() == GAAT_DATASET_LIST)
                        {
                            VSIStatBufL sStat;
                            if ((!lastArg.empty() && lastArg.back() == '/') ||
                                VSIStatL(lastArg.c_str(), &sStat) != 0)
                            {
                                arg = newArg;
                            }
                        }
                    }
                }
            }
            if (arg)
            {
                ret = arg->GetAutoCompleteChoices(lastArg);
            }
        }
    }

    return ret;
}

/************************************************************************/
/*                   GDALAlgorithm::GetFieldIndices()                   */
/************************************************************************/

bool GDALAlgorithm::GetFieldIndices(const std::vector<std::string> &names,
                                    OGRLayerH hLayer, std::vector<int> &indices)
{
    VALIDATE_POINTER1(hLayer, __func__, false);

    const OGRLayer &layer = *OGRLayer::FromHandle(hLayer);

    if (names.size() == 1 && names[0] == "ALL")
    {
        const int nSrcFieldCount = layer.GetLayerDefn()->GetFieldCount();
        for (int i = 0; i < nSrcFieldCount; ++i)
        {
            indices.push_back(i);
        }
    }
    else if (!names.empty() && !(names.size() == 1 && names[0] == "NONE"))
    {
        std::set<int> fieldsAdded;
        for (const std::string &osFieldName : names)
        {

            const int nIdx =
                layer.GetLayerDefn()->GetFieldIndex(osFieldName.c_str());

            if (nIdx < 0)
            {
                CPLError(CE_Failure, CPLE_AppDefined,
                         "Field '%s' does not exist in layer '%s'",
                         osFieldName.c_str(), layer.GetName());
                return false;
            }

            if (fieldsAdded.insert(nIdx).second)
            {
                indices.push_back(nIdx);
            }
        }
    }

    return true;
}

/************************************************************************/
/*              GDALAlgorithm::ExtractLastOptionAndValue()              */
/************************************************************************/

void GDALAlgorithm::ExtractLastOptionAndValue(std::vector<std::string> &args,
                                              std::string &option,
                                              std::string &value) const
{
    if (!args.empty() && !args.back().empty() && args.back()[0] == '-')
    {
        const auto nPosEqual = args.back().find('=');
        if (nPosEqual == std::string::npos)
        {
            // Deal with "gdal ... --option"
            if (GetArg(args.back()))
            {
                option = args.back();
                args.pop_back();
            }
        }
        else
        {
            // Deal with "gdal ... --option=<value>"
            if (GetArg(args.back().substr(0, nPosEqual)))
            {
                option = args.back().substr(0, nPosEqual);
                value = args.back().substr(nPosEqual + 1);
                args.pop_back();
            }
        }
    }
    else if (args.size() >= 2 && !args[args.size() - 2].empty() &&
             args[args.size() - 2][0] == '-')
    {
        // Deal with "gdal ... --option <value>"
        auto arg = GetArg(args[args.size() - 2]);
        if (arg && arg->GetType() != GAAT_BOOLEAN)
        {
            option = args[args.size() - 2];
            value = args.back();
            args.pop_back();
        }
    }

    const auto IsKeyValueOption = [](const std::string &osStr)
    {
        return osStr == "--co" || osStr == "--creation-option" ||
               osStr == "--lco" || osStr == "--layer-creation-option" ||
               osStr == "--oo" || osStr == "--open-option";
    };

    if (IsKeyValueOption(option))
    {
        const auto nPosEqual = value.find('=');
        if (nPosEqual != std::string::npos)
        {
            value.resize(nPosEqual);
        }
    }
}

/************************************************************************/
/*                 GDALAlgorithm::GetArgDependencies()                  */
/************************************************************************/

std::vector<std::string>
GDALAlgorithm::GetArgDependencies(const std::string &osName) const
{
    const auto arg = GetArg(osName, false);
    if (!arg)
    {
        ReportError(CE_Failure, CPLE_AppDefined, "Argument '%s' does not exist",
                    osName.c_str());
        return {};
    }
    std::vector<std::string> dependencies = arg->GetDirectDependencies();
    if (const auto &mutualDependencyGroup = arg->GetMutualDependencyGroup();
        !mutualDependencyGroup.empty())
    {
        for (const auto &otherArg : m_args)
        {
            if (otherArg.get() == arg ||
                mutualDependencyGroup.compare(
                    otherArg->GetMutualDependencyGroup()) != 0)
                continue;
            dependencies.push_back(otherArg->GetName());
        }
    }
    return dependencies;
}

//! @cond Doxygen_Suppress

/************************************************************************/
/*                  GDALContainerAlgorithm::RunImpl()                   */
/************************************************************************/

bool GDALContainerAlgorithm::RunImpl(GDALProgressFunc, void *)
{
    return false;
}

//! @endcond

/************************************************************************/
/*                        GDALAlgorithmRelease()                        */
/************************************************************************/

/** Release a handle to an algorithm.
 *
 * @since 3.11
 */
void GDALAlgorithmRelease(GDALAlgorithmH hAlg)
{
    delete hAlg;
}

/************************************************************************/
/*                        GDALAlgorithmGetName()                        */
/************************************************************************/

/** Return the algorithm name.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @return algorithm name whose lifetime is bound to hAlg and which must not
 * be freed.
 * @since 3.11
 */
const char *GDALAlgorithmGetName(GDALAlgorithmH hAlg)
{
    VALIDATE_POINTER1(hAlg, __func__, nullptr);
    return hAlg->ptr->GetName().c_str();
}

/************************************************************************/
/*                    GDALAlgorithmGetDescription()                     */
/************************************************************************/

/** Return the algorithm (short) description.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @return algorithm description whose lifetime is bound to hAlg and which must
 * not be freed.
 * @since 3.11
 */
const char *GDALAlgorithmGetDescription(GDALAlgorithmH hAlg)
{
    VALIDATE_POINTER1(hAlg, __func__, nullptr);
    return hAlg->ptr->GetDescription().c_str();
}

/************************************************************************/
/*                  GDALAlgorithmGetLongDescription()                   */
/************************************************************************/

/** Return the algorithm (longer) description.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @return algorithm description whose lifetime is bound to hAlg and which must
 * not be freed.
 * @since 3.11
 */
const char *GDALAlgorithmGetLongDescription(GDALAlgorithmH hAlg)
{
    VALIDATE_POINTER1(hAlg, __func__, nullptr);
    return hAlg->ptr->GetLongDescription().c_str();
}

/************************************************************************/
/*                    GDALAlgorithmGetHelpFullURL()                     */
/************************************************************************/

/** Return the algorithm full URL.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @return algorithm URL whose lifetime is bound to hAlg and which must
 * not be freed.
 * @since 3.11
 */
const char *GDALAlgorithmGetHelpFullURL(GDALAlgorithmH hAlg)
{
    VALIDATE_POINTER1(hAlg, __func__, nullptr);
    return hAlg->ptr->GetHelpFullURL().c_str();
}

/************************************************************************/
/*                   GDALAlgorithmHasSubAlgorithms()                    */
/************************************************************************/

/** Return whether the algorithm has sub-algorithms.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @since 3.11
 */
bool GDALAlgorithmHasSubAlgorithms(GDALAlgorithmH hAlg)
{
    VALIDATE_POINTER1(hAlg, __func__, false);
    return hAlg->ptr->HasSubAlgorithms();
}

/************************************************************************/
/*                 GDALAlgorithmGetSubAlgorithmNames()                  */
/************************************************************************/

/** Get the names of registered algorithms.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @return a NULL terminated list of names, which must be destroyed with
 * CSLDestroy()
 * @since 3.11
 */
char **GDALAlgorithmGetSubAlgorithmNames(GDALAlgorithmH hAlg)
{
    VALIDATE_POINTER1(hAlg, __func__, nullptr);
    return CPLStringList(hAlg->ptr->GetSubAlgorithmNames()).StealList();
}

/************************************************************************/
/*                GDALAlgorithmInstantiateSubAlgorithm()                */
/************************************************************************/

/** Instantiate an algorithm by its name (or its alias).
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @param pszSubAlgName Algorithm name. Must NOT be null.
 * @return an handle to the algorithm (to be freed with GDALAlgorithmRelease),
 * or NULL if the algorithm does not exist or another error occurred.
 * @since 3.11
 */
GDALAlgorithmH GDALAlgorithmInstantiateSubAlgorithm(GDALAlgorithmH hAlg,
                                                    const char *pszSubAlgName)
{
    VALIDATE_POINTER1(hAlg, __func__, nullptr);
    VALIDATE_POINTER1(pszSubAlgName, __func__, nullptr);
    auto subAlg = hAlg->ptr->InstantiateSubAlgorithm(pszSubAlgName);
    return subAlg
               ? std::make_unique<GDALAlgorithmHS>(std::move(subAlg)).release()
               : nullptr;
}

/************************************************************************/
/*               GDALAlgorithmParseCommandLineArguments()               */
/************************************************************************/

/** Parse a command line argument, which does not include the algorithm
 * name, to set the value of corresponding arguments.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @param papszArgs NULL-terminated list of arguments, not including the algorithm name.
 * @return true if successful, false otherwise
 * @since 3.11
 */

bool GDALAlgorithmParseCommandLineArguments(GDALAlgorithmH hAlg,
                                            CSLConstList papszArgs)
{
    VALIDATE_POINTER1(hAlg, __func__, false);
    return hAlg->ptr->ParseCommandLineArguments(CPLStringList(papszArgs));
}

/************************************************************************/
/*                  GDALAlgorithmGetActualAlgorithm()                   */
/************************************************************************/

/** Return the actual algorithm that is going to be invoked, when the
 * current algorithm has sub-algorithms.
 *
 * Only valid after GDALAlgorithmParseCommandLineArguments() has been called.
 *
 * Note that the lifetime of the returned algorithm does not exceed the one of
 * the hAlg instance that owns it.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @return an handle to the algorithm (to be freed with GDALAlgorithmRelease).
 * @since 3.11
 */
GDALAlgorithmH GDALAlgorithmGetActualAlgorithm(GDALAlgorithmH hAlg)
{
    VALIDATE_POINTER1(hAlg, __func__, nullptr);
    return GDALAlgorithmHS::FromRef(hAlg->ptr->GetActualAlgorithm()).release();
}

/************************************************************************/
/*                          GDALAlgorithmRun()                          */
/************************************************************************/

/** Execute the algorithm, starting with ValidateArguments() and then
 * calling RunImpl().
 *
 * This function must be called at most once per instance.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @param pfnProgress Progress callback. May be null.
 * @param pProgressData Progress callback user data. May be null.
 * @return true if successful, false otherwise
 * @since 3.11
 */

bool GDALAlgorithmRun(GDALAlgorithmH hAlg, GDALProgressFunc pfnProgress,
                      void *pProgressData)
{
    VALIDATE_POINTER1(hAlg, __func__, false);
    return hAlg->ptr->Run(pfnProgress, pProgressData);
}

/************************************************************************/
/*                       GDALAlgorithmFinalize()                        */
/************************************************************************/

/** Complete any pending actions, and return the final status.
 * This is typically useful for algorithm that generate an output dataset.
 *
 * Note that this function does *NOT* release memory associated with the
 * algorithm. GDALAlgorithmRelease() must still be called afterwards.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @return true if successful, false otherwise
 * @since 3.11
 */

bool GDALAlgorithmFinalize(GDALAlgorithmH hAlg)
{
    VALIDATE_POINTER1(hAlg, __func__, false);
    return hAlg->ptr->Finalize();
}

/************************************************************************/
/*                    GDALAlgorithmGetUsageAsJSON()                     */
/************************************************************************/

/** Return the usage of the algorithm as a JSON-serialized string.
 *
 * This can be used to dynamically generate interfaces to algorithms.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @return a string that must be freed with CPLFree()
 * @since 3.11
 */
char *GDALAlgorithmGetUsageAsJSON(GDALAlgorithmH hAlg)
{
    VALIDATE_POINTER1(hAlg, __func__, nullptr);
    return CPLStrdup(hAlg->ptr->GetUsageAsJSON().c_str());
}

/************************************************************************/
/*                      GDALAlgorithmGetArgNames()                      */
/************************************************************************/

/** Return the list of available argument names.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @return a NULL terminated list of names, which must be destroyed with
 * CSLDestroy()
 * @since 3.11
 */
char **GDALAlgorithmGetArgNames(GDALAlgorithmH hAlg)
{
    VALIDATE_POINTER1(hAlg, __func__, nullptr);
    CPLStringList list;
    for (const auto &arg : hAlg->ptr->GetArgs())
        list.AddString(arg->GetName().c_str());
    return list.StealList();
}

/************************************************************************/
/*                        GDALAlgorithmGetArg()                         */
/************************************************************************/

/** Return an argument from its name.
 *
 * The lifetime of the returned object does not exceed the one of hAlg.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @param pszArgName Argument name. Must NOT be null.
 * @return an argument that must be released with GDALAlgorithmArgRelease(),
 * or nullptr in case of error
 * @since 3.11
 */
GDALAlgorithmArgH GDALAlgorithmGetArg(GDALAlgorithmH hAlg,
                                      const char *pszArgName)
{
    VALIDATE_POINTER1(hAlg, __func__, nullptr);
    VALIDATE_POINTER1(pszArgName, __func__, nullptr);
    auto arg = hAlg->ptr->GetArg(pszArgName, /* suggestionAllowed = */ true,
                                 /* isConst = */ true);
    if (!arg)
        return nullptr;
    return std::make_unique<GDALAlgorithmArgHS>(arg).release();
}

/************************************************************************/
/*                    GDALAlgorithmGetArgNonConst()                     */
/************************************************************************/

/** Return an argument from its name, possibly allowing creation of user-provided
 * argument if the algorithm allow it.
 *
 * The lifetime of the returned object does not exceed the one of hAlg.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @param pszArgName Argument name. Must NOT be null.
 * @return an argument that must be released with GDALAlgorithmArgRelease(),
 * or nullptr in case of error
 * @since 3.12
 */
GDALAlgorithmArgH GDALAlgorithmGetArgNonConst(GDALAlgorithmH hAlg,
                                              const char *pszArgName)
{
    VALIDATE_POINTER1(hAlg, __func__, nullptr);
    VALIDATE_POINTER1(pszArgName, __func__, nullptr);
    auto arg = hAlg->ptr->GetArg(pszArgName, /* suggestionAllowed = */ true,
                                 /* isConst = */ false);
    if (!arg)
        return nullptr;
    return std::make_unique<GDALAlgorithmArgHS>(arg).release();
}

/************************************************************************/
/*                  GDALAlgorithmGetArgDependencies()                   */
/************************************************************************/

/** Return the list of argument names the specified argument depends on.
 *
 *  This includes both regular dependencies and mutual dependencies.
 *
 * @param hAlg Handle to an algorithm. Must NOT be null.
 * @param pszArgName Argument name. Must NOT be null.
 * @return a NULL terminated list of names, which must be destroyed with
 * CSLDestroy()
 * @since 3.11
 */
char **GDALAlgorithmGetArgDependencies(GDALAlgorithmH hAlg,
                                       const char *pszArgName)
{
    VALIDATE_POINTER1(hAlg, __func__, nullptr);
    VALIDATE_POINTER1(pszArgName, __func__, nullptr);
    return CPLStringList(hAlg->ptr->GetArgDependencies(pszArgName)).StealList();
}

/************************************************************************/
/*                      GDALAlgorithmArgRelease()                       */
/************************************************************************/

/** Release a handle to an argument.
 *
 * @since 3.11
 */
void GDALAlgorithmArgRelease(GDALAlgorithmArgH hArg)
{
    delete hArg;
}

/************************************************************************/
/*                      GDALAlgorithmArgGetName()                       */
/************************************************************************/

/** Return the name of an argument.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return argument name whose lifetime is bound to hArg and which must not
 * be freed.
 * @since 3.11
 */
const char *GDALAlgorithmArgGetName(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    return hArg->ptr->GetName().c_str();
}

/************************************************************************/
/*                      GDALAlgorithmArgGetType()                       */
/************************************************************************/

/** Get the type of an argument
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
GDALAlgorithmArgType GDALAlgorithmArgGetType(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, GAAT_STRING);
    return hArg->ptr->GetType();
}

/************************************************************************/
/*                   GDALAlgorithmArgGetDescription()                   */
/************************************************************************/

/** Return the description of an argument.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return argument description whose lifetime is bound to hArg and which must not
 * be freed.
 * @since 3.11
 */
const char *GDALAlgorithmArgGetDescription(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    return hArg->ptr->GetDescription().c_str();
}

/************************************************************************/
/*                    GDALAlgorithmArgGetShortName()                    */
/************************************************************************/

/** Return the short name, or empty string if there is none
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return short name whose lifetime is bound to hArg and which must not
 * be freed.
 * @since 3.11
 */
const char *GDALAlgorithmArgGetShortName(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    return hArg->ptr->GetShortName().c_str();
}

/************************************************************************/
/*                     GDALAlgorithmArgGetAliases()                     */
/************************************************************************/

/** Return the aliases (potentially none)
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return a NULL terminated list of names, which must be destroyed with
 * CSLDestroy()

 * @since 3.11
 */
char **GDALAlgorithmArgGetAliases(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    return CPLStringList(hArg->ptr->GetAliases()).StealList();
}

/************************************************************************/
/*                     GDALAlgorithmArgGetMetaVar()                     */
/************************************************************************/

/** Return the "meta-var" hint.
 *
 * By default, the meta-var value is the long name of the argument in
 * upper case.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return meta-var hint whose lifetime is bound to hArg and which must not
 * be freed.
 * @since 3.11
 */
const char *GDALAlgorithmArgGetMetaVar(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    return hArg->ptr->GetMetaVar().c_str();
}

/************************************************************************/
/*                    GDALAlgorithmArgGetCategory()                     */
/************************************************************************/

/** Return the argument category
 *
 * GAAC_COMMON, GAAC_BASE, GAAC_ADVANCED, GAAC_ESOTERIC or a custom category.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return category whose lifetime is bound to hArg and which must not
 * be freed.
 * @since 3.11
 */
const char *GDALAlgorithmArgGetCategory(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    return hArg->ptr->GetCategory().c_str();
}

/************************************************************************/
/*                    GDALAlgorithmArgIsPositional()                    */
/************************************************************************/

/** Return if the argument is a positional one.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
bool GDALAlgorithmArgIsPositional(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->IsPositional();
}

/************************************************************************/
/*                     GDALAlgorithmArgIsRequired()                     */
/************************************************************************/

/** Return whether the argument is required. Defaults to false.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
bool GDALAlgorithmArgIsRequired(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->IsRequired();
}

/************************************************************************/
/*                    GDALAlgorithmArgGetMinCount()                     */
/************************************************************************/

/** Return the minimum number of values for the argument.
 *
 * Defaults to 0.
 * Only applies to list type of arguments.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
int GDALAlgorithmArgGetMinCount(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, 0);
    return hArg->ptr->GetMinCount();
}

/************************************************************************/
/*                    GDALAlgorithmArgGetMaxCount()                     */
/************************************************************************/

/** Return the maximum number of values for the argument.
 *
 * Defaults to 1 for scalar types, and INT_MAX for list types.
 * Only applies to list type of arguments.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
int GDALAlgorithmArgGetMaxCount(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, 0);
    return hArg->ptr->GetMaxCount();
}

/************************************************************************/
/*               GDALAlgorithmArgGetPackedValuesAllowed()               */
/************************************************************************/

/** Return whether, for list type of arguments, several values, space
 * separated, may be specified. That is "--foo=bar,baz".
 * The default is true.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
bool GDALAlgorithmArgGetPackedValuesAllowed(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->GetPackedValuesAllowed();
}

/************************************************************************/
/*               GDALAlgorithmArgGetRepeatedArgAllowed()                */
/************************************************************************/

/** Return whether, for list type of arguments, the argument may be
 * repeated. That is "--foo=bar --foo=baz".
 * The default is true.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
bool GDALAlgorithmArgGetRepeatedArgAllowed(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->GetRepeatedArgAllowed();
}

/************************************************************************/
/*                     GDALAlgorithmArgGetChoices()                     */
/************************************************************************/

/** Return the allowed values (as strings) for the argument.
 *
 * Only honored for GAAT_STRING and GAAT_STRING_LIST types.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return a NULL terminated list of names, which must be destroyed with
 * CSLDestroy()

 * @since 3.11
 */
char **GDALAlgorithmArgGetChoices(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    return CPLStringList(hArg->ptr->GetChoices()).StealList();
}

/************************************************************************/
/*                  GDALAlgorithmArgGetMetadataItem()                   */
/************************************************************************/

/** Return the values of the metadata item of an argument.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param pszItem Name of the item. Must NOT be null.
 * @return a NULL terminated list of values, which must be destroyed with
 * CSLDestroy()

 * @since 3.11
 */
char **GDALAlgorithmArgGetMetadataItem(GDALAlgorithmArgH hArg,
                                       const char *pszItem)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    VALIDATE_POINTER1(pszItem, __func__, nullptr);
    const auto pVecOfStrings = hArg->ptr->GetMetadataItem(pszItem);
    return pVecOfStrings ? CPLStringList(*pVecOfStrings).StealList() : nullptr;
}

/************************************************************************/
/*                  GDALAlgorithmArgIsExplicitlySet()                   */
/************************************************************************/

/** Return whether the argument value has been explicitly set with Set()
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
bool GDALAlgorithmArgIsExplicitlySet(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->IsExplicitlySet();
}

/************************************************************************/
/*                  GDALAlgorithmArgHasDefaultValue()                   */
/************************************************************************/

/** Return if the argument has a declared default value.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
bool GDALAlgorithmArgHasDefaultValue(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->HasDefaultValue();
}

/************************************************************************/
/*                GDALAlgorithmArgGetDefaultAsBoolean()                 */
/************************************************************************/

/** Return the argument default value as a integer.
 *
 * Must only be called on arguments whose type is GAAT_BOOLEAN
 *
 * GDALAlgorithmArgHasDefaultValue() must be called to determine if the
 * argument has a default value.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.12
 */
bool GDALAlgorithmArgGetDefaultAsBoolean(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    if (hArg->ptr->GetType() != GAAT_BOOLEAN)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "%s must only be called on arguments of type GAAT_BOOLEAN",
                 __func__);
        return false;
    }
    return hArg->ptr->GetDefault<bool>();
}

/************************************************************************/
/*                 GDALAlgorithmArgGetDefaultAsString()                 */
/************************************************************************/

/** Return the argument default value as a string.
 *
 * Must only be called on arguments whose type is GAAT_STRING.
 *
 * GDALAlgorithmArgHasDefaultValue() must be called to determine if the
 * argument has a default value.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return string whose lifetime is bound to hArg and which must not
 * be freed.
 * @since 3.11
 */
const char *GDALAlgorithmArgGetDefaultAsString(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    if (hArg->ptr->GetType() != GAAT_STRING)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "%s must only be called on arguments of type GAAT_STRING",
                 __func__);
        return nullptr;
    }
    return hArg->ptr->GetDefault<std::string>().c_str();
}

/************************************************************************/
/*                GDALAlgorithmArgGetDefaultAsInteger()                 */
/************************************************************************/

/** Return the argument default value as a integer.
 *
 * Must only be called on arguments whose type is GAAT_INTEGER
 *
 * GDALAlgorithmArgHasDefaultValue() must be called to determine if the
 * argument has a default value.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.12
 */
int GDALAlgorithmArgGetDefaultAsInteger(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, 0);
    if (hArg->ptr->GetType() != GAAT_INTEGER)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "%s must only be called on arguments of type GAAT_INTEGER",
                 __func__);
        return 0;
    }
    return hArg->ptr->GetDefault<int>();
}

/************************************************************************/
/*                 GDALAlgorithmArgGetDefaultAsDouble()                 */
/************************************************************************/

/** Return the argument default value as a double.
 *
 * Must only be called on arguments whose type is GAAT_REAL
 *
 * GDALAlgorithmArgHasDefaultValue() must be called to determine if the
 * argument has a default value.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.12
 */
double GDALAlgorithmArgGetDefaultAsDouble(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, 0);
    if (hArg->ptr->GetType() != GAAT_REAL)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "%s must only be called on arguments of type GAAT_REAL",
                 __func__);
        return 0;
    }
    return hArg->ptr->GetDefault<double>();
}

/************************************************************************/
/*               GDALAlgorithmArgGetDefaultAsStringList()               */
/************************************************************************/

/** Return the argument default value as a string list.
 *
 * Must only be called on arguments whose type is GAAT_STRING_LIST.
 *
 * GDALAlgorithmArgHasDefaultValue() must be called to determine if the
 * argument has a default value.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return a NULL terminated list of names, which must be destroyed with
 * CSLDestroy()

 * @since 3.12
 */
char **GDALAlgorithmArgGetDefaultAsStringList(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    if (hArg->ptr->GetType() != GAAT_STRING_LIST)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "%s must only be called on arguments of type GAAT_STRING_LIST",
                 __func__);
        return nullptr;
    }
    return CPLStringList(hArg->ptr->GetDefault<std::vector<std::string>>())
        .StealList();
}

/************************************************************************/
/*              GDALAlgorithmArgGetDefaultAsIntegerList()               */
/************************************************************************/

/** Return the argument default value as a integer list.
 *
 * Must only be called on arguments whose type is GAAT_INTEGER_LIST.
 *
 * GDALAlgorithmArgHasDefaultValue() must be called to determine if the
 * argument has a default value.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param[out] pnCount Pointer to the number of values in the list. Must NOT be null.
 * @since 3.12
 */
const int *GDALAlgorithmArgGetDefaultAsIntegerList(GDALAlgorithmArgH hArg,
                                                   size_t *pnCount)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    VALIDATE_POINTER1(pnCount, __func__, nullptr);
    if (hArg->ptr->GetType() != GAAT_INTEGER_LIST)
    {
        CPLError(
            CE_Failure, CPLE_AppDefined,
            "%s must only be called on arguments of type GAAT_INTEGER_LIST",
            __func__);
        *pnCount = 0;
        return nullptr;
    }
    const auto &val = hArg->ptr->GetDefault<std::vector<int>>();
    *pnCount = val.size();
    return val.data();
}

/************************************************************************/
/*               GDALAlgorithmArgGetDefaultAsDoubleList()               */
/************************************************************************/

/** Return the argument default value as a real list.
 *
 * Must only be called on arguments whose type is GAAT_REAL_LIST.
 *
 * GDALAlgorithmArgHasDefaultValue() must be called to determine if the
 * argument has a default value.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param[out] pnCount Pointer to the number of values in the list. Must NOT be null.
 * @since 3.12
 */
const double *GDALAlgorithmArgGetDefaultAsDoubleList(GDALAlgorithmArgH hArg,
                                                     size_t *pnCount)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    VALIDATE_POINTER1(pnCount, __func__, nullptr);
    if (hArg->ptr->GetType() != GAAT_REAL_LIST)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "%s must only be called on arguments of type GAAT_REAL_LIST",
                 __func__);
        *pnCount = 0;
        return nullptr;
    }
    const auto &val = hArg->ptr->GetDefault<std::vector<double>>();
    *pnCount = val.size();
    return val.data();
}

/************************************************************************/
/*                      GDALAlgorithmArgIsHidden()                      */
/************************************************************************/

/** Return whether the argument is hidden (for GDAL internal use)
 *
 * This is an alias for GDALAlgorithmArgIsHiddenForCLI() &&
 * GDALAlgorithmArgIsHiddenForAPI().
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.12
 */
bool GDALAlgorithmArgIsHidden(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->IsHidden();
}

/************************************************************************/
/*                   GDALAlgorithmArgIsHiddenForCLI()                   */
/************************************************************************/

/** Return whether the argument must not be mentioned in CLI usage.
 *
 * For example, "output-value" for "gdal raster info", which is only
 * meant when the algorithm is used from a non-CLI context.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
bool GDALAlgorithmArgIsHiddenForCLI(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->IsHiddenForCLI();
}

/************************************************************************/
/*                   GDALAlgorithmArgIsHiddenForAPI()                   */
/************************************************************************/

/** Return whether the argument must not be mentioned in the context of an
 * API use.
 * Said otherwise, if it is only for CLI usage.
 *
 * For example "--help"
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.12
 */
bool GDALAlgorithmArgIsHiddenForAPI(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->IsHiddenForAPI();
}

/************************************************************************/
/*                    GDALAlgorithmArgIsOnlyForCLI()                    */
/************************************************************************/

/** Return whether the argument must not be mentioned in the context of an
 * API use.
 * Said otherwise, if it is only for CLI usage.
 *
 * For example "--help"
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 * @deprecated Use GDALAlgorithmArgIsHiddenForAPI() instead.
 */
bool GDALAlgorithmArgIsOnlyForCLI(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->IsHiddenForAPI();
}

/************************************************************************/
/*             GDALAlgorithmArgIsAvailableInPipelineStep()              */
/************************************************************************/

/** Return whether the argument is available in a pipeline step.
 *
 * If false, it is only available in standalone mode.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.13
 */
bool GDALAlgorithmArgIsAvailableInPipelineStep(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->IsAvailableInPipelineStep();
}

/************************************************************************/
/*                      GDALAlgorithmArgIsInput()                       */
/************************************************************************/

/** Indicate whether the value of the argument is read-only during the
 * execution of the algorithm.
 *
 * Default is true.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
bool GDALAlgorithmArgIsInput(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->IsInput();
}

/************************************************************************/
/*                      GDALAlgorithmArgIsOutput()                      */
/************************************************************************/

/** Return whether (at least part of) the value of the argument is set
 * during the execution of the algorithm.
 *
 * For example, "output-value" for "gdal raster info"
 * Default is false.
 * An argument may return both IsInput() and IsOutput() as true.
 * For example the "gdal raster convert" algorithm consumes the dataset
 * name of its "output" argument, and sets the dataset object during its
 * execution.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
bool GDALAlgorithmArgIsOutput(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->IsOutput();
}

/************************************************************************/
/*                   GDALAlgorithmArgGetDatasetType()                   */
/************************************************************************/

/** Get which type of dataset is allowed / generated.
 *
 * Binary-or combination of GDAL_OF_RASTER, GDAL_OF_VECTOR and
 * GDAL_OF_MULTIDIM_RASTER.
 * Only applies to arguments of type GAAT_DATASET or GAAT_DATASET_LIST.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
GDALArgDatasetType GDALAlgorithmArgGetDatasetType(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, 0);
    return hArg->ptr->GetDatasetType();
}

/************************************************************************/
/*                GDALAlgorithmArgGetDatasetInputFlags()                */
/************************************************************************/

/** Indicates which components among name and dataset are accepted as
 * input, when this argument serves as an input.
 *
 * If the GADV_NAME bit is set, it indicates a dataset name is accepted as
 * input.
 * If the GADV_OBJECT bit is set, it indicates a dataset object is
 * accepted as input.
 * If both bits are set, the algorithm can accept either a name or a dataset
 * object.
 * Only applies to arguments of type GAAT_DATASET or GAAT_DATASET_LIST.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return string whose lifetime is bound to hAlg and which must not
 * be freed.
 * @since 3.11
 */
int GDALAlgorithmArgGetDatasetInputFlags(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, 0);
    return hArg->ptr->GetDatasetInputFlags();
}

/************************************************************************/
/*               GDALAlgorithmArgGetDatasetOutputFlags()                */
/************************************************************************/

/** Indicates which components among name and dataset are modified,
 * when this argument serves as an output.
 *
 * If the GADV_NAME bit is set, it indicates a dataset name is generated as
 * output (that is the algorithm will generate the name. Rarely used).
 * If the GADV_OBJECT bit is set, it indicates a dataset object is
 * generated as output, and available for use after the algorithm has
 * completed.
 * Only applies to arguments of type GAAT_DATASET or GAAT_DATASET_LIST.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return string whose lifetime is bound to hAlg and which must not
 * be freed.
 * @since 3.11
 */
int GDALAlgorithmArgGetDatasetOutputFlags(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, 0);
    return hArg->ptr->GetDatasetOutputFlags();
}

/************************************************************************/
/*              GDALAlgorithmArgGetMutualExclusionGroup()               */
/************************************************************************/

/** Return the name of the mutual exclusion group to which this argument
 * belongs to.
 *
 * Or empty string if it does not belong to any exclusion group.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return string whose lifetime is bound to hArg and which must not
 * be freed.
 * @since 3.11
 */
const char *GDALAlgorithmArgGetMutualExclusionGroup(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    return hArg->ptr->GetMutualExclusionGroup().c_str();
}

/************************************************************************/
/*              GDALAlgorithmArgGetMutualDependencyGroup()              */
/************************************************************************/

/** Return the name of the mutual dependency group to which this argument
 * belongs to.
 *
 * Or empty string if it does not belong to any dependency group.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return string whose lifetime is bound to hArg and which must not
 * be freed.
 * @since 3.13
 */
const char *GDALAlgorithmArgGetMutualDependencyGroup(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    return hArg->ptr->GetMutualDependencyGroup().c_str();
}

/************************************************************************/
/*               GDALAlgorithmArgGetDirectDependencies()                */
/************************************************************************/

/** Return the list of names of arguments that this argument depends on.
 *
 *  This is not necessarily a symmetric relationship.
 *  If argument A depends on argument B, it doesn't mean that B depends on A.
 *  Mutual dependency groups are a special case of dependencies,
 *  where all arguments of the group depend on each other and are not
 *  returned by this method.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return a NULL terminated list of names, which must be destroyed with
 * CSLDestroy()
 * @since 3.13
 */
char **GDALAlgorithmArgGetDirectDependencies(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    return CPLStringList(hArg->ptr->GetDirectDependencies()).StealList();
}

/************************************************************************/
/*                    GDALAlgorithmArgGetAsBoolean()                    */
/************************************************************************/

/** Return the argument value as a boolean.
 *
 * Must only be called on arguments whose type is GAAT_BOOLEAN.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
bool GDALAlgorithmArgGetAsBoolean(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    if (hArg->ptr->GetType() != GAAT_BOOLEAN)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "%s must only be called on arguments of type GAAT_BOOLEAN",
                 __func__);
        return false;
    }
    return hArg->ptr->Get<bool>();
}

/************************************************************************/
/*                    GDALAlgorithmArgGetAsString()                     */
/************************************************************************/

/** Return the argument value as a string.
 *
 * Must only be called on arguments whose type is GAAT_STRING.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return string whose lifetime is bound to hArg and which must not
 * be freed.
 * @since 3.11
 */
const char *GDALAlgorithmArgGetAsString(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    if (hArg->ptr->GetType() != GAAT_STRING)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "%s must only be called on arguments of type GAAT_STRING",
                 __func__);
        return nullptr;
    }
    return hArg->ptr->Get<std::string>().c_str();
}

/************************************************************************/
/*                 GDALAlgorithmArgGetAsDatasetValue()                  */
/************************************************************************/

/** Return the argument value as a GDALArgDatasetValueH.
 *
 * Must only be called on arguments whose type is GAAT_DATASET
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return handle to a GDALArgDatasetValue that must be released with
 * GDALArgDatasetValueRelease(). The lifetime of that handle does not exceed
 * the one of hArg.
 * @since 3.11
 */
GDALArgDatasetValueH GDALAlgorithmArgGetAsDatasetValue(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    if (hArg->ptr->GetType() != GAAT_DATASET)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "%s must only be called on arguments of type GAAT_DATASET",
                 __func__);
        return nullptr;
    }
    return std::make_unique<GDALArgDatasetValueHS>(
               &(hArg->ptr->Get<GDALArgDatasetValue>()))
        .release();
}

/************************************************************************/
/*                    GDALAlgorithmArgGetAsInteger()                    */
/************************************************************************/

/** Return the argument value as a integer.
 *
 * Must only be called on arguments whose type is GAAT_INTEGER
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
int GDALAlgorithmArgGetAsInteger(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, 0);
    if (hArg->ptr->GetType() != GAAT_INTEGER)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "%s must only be called on arguments of type GAAT_INTEGER",
                 __func__);
        return 0;
    }
    return hArg->ptr->Get<int>();
}

/************************************************************************/
/*                    GDALAlgorithmArgGetAsDouble()                     */
/************************************************************************/

/** Return the argument value as a double.
 *
 * Must only be called on arguments whose type is GAAT_REAL
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @since 3.11
 */
double GDALAlgorithmArgGetAsDouble(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, 0);
    if (hArg->ptr->GetType() != GAAT_REAL)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "%s must only be called on arguments of type GAAT_REAL",
                 __func__);
        return 0;
    }
    return hArg->ptr->Get<double>();
}

/************************************************************************/
/*                  GDALAlgorithmArgGetAsStringList()                   */
/************************************************************************/

/** Return the argument value as a string list.
 *
 * Must only be called on arguments whose type is GAAT_STRING_LIST.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @return a NULL terminated list of names, which must be destroyed with
 * CSLDestroy()

 * @since 3.11
 */
char **GDALAlgorithmArgGetAsStringList(GDALAlgorithmArgH hArg)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    if (hArg->ptr->GetType() != GAAT_STRING_LIST)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "%s must only be called on arguments of type GAAT_STRING_LIST",
                 __func__);
        return nullptr;
    }
    return CPLStringList(hArg->ptr->Get<std::vector<std::string>>())
        .StealList();
}

/************************************************************************/
/*                  GDALAlgorithmArgGetAsIntegerList()                  */
/************************************************************************/

/** Return the argument value as a integer list.
 *
 * Must only be called on arguments whose type is GAAT_INTEGER_LIST.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param[out] pnCount Pointer to the number of values in the list. Must NOT be null.
 * @since 3.11
 */
const int *GDALAlgorithmArgGetAsIntegerList(GDALAlgorithmArgH hArg,
                                            size_t *pnCount)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    VALIDATE_POINTER1(pnCount, __func__, nullptr);
    if (hArg->ptr->GetType() != GAAT_INTEGER_LIST)
    {
        CPLError(
            CE_Failure, CPLE_AppDefined,
            "%s must only be called on arguments of type GAAT_INTEGER_LIST",
            __func__);
        *pnCount = 0;
        return nullptr;
    }
    const auto &val = hArg->ptr->Get<std::vector<int>>();
    *pnCount = val.size();
    return val.data();
}

/************************************************************************/
/*                  GDALAlgorithmArgGetAsDoubleList()                   */
/************************************************************************/

/** Return the argument value as a real list.
 *
 * Must only be called on arguments whose type is GAAT_REAL_LIST.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param[out] pnCount Pointer to the number of values in the list. Must NOT be null.
 * @since 3.11
 */
const double *GDALAlgorithmArgGetAsDoubleList(GDALAlgorithmArgH hArg,
                                              size_t *pnCount)
{
    VALIDATE_POINTER1(hArg, __func__, nullptr);
    VALIDATE_POINTER1(pnCount, __func__, nullptr);
    if (hArg->ptr->GetType() != GAAT_REAL_LIST)
    {
        CPLError(CE_Failure, CPLE_AppDefined,
                 "%s must only be called on arguments of type GAAT_REAL_LIST",
                 __func__);
        *pnCount = 0;
        return nullptr;
    }
    const auto &val = hArg->ptr->Get<std::vector<double>>();
    *pnCount = val.size();
    return val.data();
}

/************************************************************************/
/*                    GDALAlgorithmArgSetAsBoolean()                    */
/************************************************************************/

/** Set the value for a GAAT_BOOLEAN argument.
 *
 * It cannot be called several times for a given argument.
 * Validation checks and other actions are run.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param value value.
 * @return true if success.
 * @since 3.11
 */

bool GDALAlgorithmArgSetAsBoolean(GDALAlgorithmArgH hArg, bool value)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->Set(value);
}

/************************************************************************/
/*                    GDALAlgorithmArgSetAsString()                     */
/************************************************************************/

/** Set the value for a GAAT_STRING argument.
 *
 * It cannot be called several times for a given argument.
 * Validation checks and other actions are run.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param value value (may be null)
 * @return true if success.
 * @since 3.11
 */

bool GDALAlgorithmArgSetAsString(GDALAlgorithmArgH hArg, const char *value)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->Set(value ? value : "");
}

/************************************************************************/
/*                    GDALAlgorithmArgSetAsInteger()                    */
/************************************************************************/

/** Set the value for a GAAT_INTEGER (or GAAT_REAL) argument.
 *
 * It cannot be called several times for a given argument.
 * Validation checks and other actions are run.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param value value.
 * @return true if success.
 * @since 3.11
 */

bool GDALAlgorithmArgSetAsInteger(GDALAlgorithmArgH hArg, int value)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->Set(value);
}

/************************************************************************/
/*                    GDALAlgorithmArgSetAsDouble()                     */
/************************************************************************/

/** Set the value for a GAAT_REAL argument.
 *
 * It cannot be called several times for a given argument.
 * Validation checks and other actions are run.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param value value.
 * @return true if success.
 * @since 3.11
 */

bool GDALAlgorithmArgSetAsDouble(GDALAlgorithmArgH hArg, double value)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->Set(value);
}

/************************************************************************/
/*                 GDALAlgorithmArgSetAsDatasetValue()                  */
/************************************************************************/

/** Set the value for a GAAT_DATASET argument.
 *
 * It cannot be called several times for a given argument.
 * Validation checks and other actions are run.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param value Handle to a GDALArgDatasetValue. Must NOT be null.
 * @return true if success.
 * @since 3.11
 */
bool GDALAlgorithmArgSetAsDatasetValue(GDALAlgorithmArgH hArg,
                                       GDALArgDatasetValueH value)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    VALIDATE_POINTER1(value, __func__, false);
    return hArg->ptr->SetFrom(*(value->ptr));
}

/************************************************************************/
/*                     GDALAlgorithmArgSetDataset()                     */
/************************************************************************/

/** Set dataset object, increasing its reference counter.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param hDS Dataset object. May be null.
 * @return true if success.
 * @since 3.11
 */

bool GDALAlgorithmArgSetDataset(GDALAlgorithmArgH hArg, GDALDatasetH hDS)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->Set(GDALDataset::FromHandle(hDS));
}

/************************************************************************/
/*                  GDALAlgorithmArgSetAsStringList()                   */
/************************************************************************/

/** Set the value for a GAAT_STRING_LIST argument.
 *
 * It cannot be called several times for a given argument.
 * Validation checks and other actions are run.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param value value as a NULL terminated list (may be null)
 * @return true if success.
 * @since 3.11
 */

bool GDALAlgorithmArgSetAsStringList(GDALAlgorithmArgH hArg, CSLConstList value)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->Set(
        static_cast<std::vector<std::string>>(CPLStringList(value)));
}

/************************************************************************/
/*                  GDALAlgorithmArgSetAsIntegerList()                  */
/************************************************************************/

/** Set the value for a GAAT_INTEGER_LIST argument.
 *
 * It cannot be called several times for a given argument.
 * Validation checks and other actions are run.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param nCount Number of values in pnValues.
 * @param pnValues Pointer to an array of integer values of size nCount.
 * @return true if success.
 * @since 3.11
 */
bool GDALAlgorithmArgSetAsIntegerList(GDALAlgorithmArgH hArg, size_t nCount,
                                      const int *pnValues)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->Set(std::vector<int>(pnValues, pnValues + nCount));
}

/************************************************************************/
/*                  GDALAlgorithmArgSetAsDoubleList()                   */
/************************************************************************/

/** Set the value for a GAAT_REAL_LIST argument.
 *
 * It cannot be called several times for a given argument.
 * Validation checks and other actions are run.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param nCount Number of values in pnValues.
 * @param pnValues Pointer to an array of double values of size nCount.
 * @return true if success.
 * @since 3.11
 */
bool GDALAlgorithmArgSetAsDoubleList(GDALAlgorithmArgH hArg, size_t nCount,
                                     const double *pnValues)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    return hArg->ptr->Set(std::vector<double>(pnValues, pnValues + nCount));
}

/************************************************************************/
/*                    GDALAlgorithmArgSetDatasets()                     */
/************************************************************************/

/** Set dataset objects to a GAAT_DATASET_LIST argument, increasing their reference counter.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param nCount Number of values in pnValues.
 * @param pahDS Pointer to an array of dataset of size nCount.
 * @return true if success.
 * @since 3.11
 */

bool GDALAlgorithmArgSetDatasets(GDALAlgorithmArgH hArg, size_t nCount,
                                 GDALDatasetH *pahDS)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    std::vector<GDALArgDatasetValue> values;
    for (size_t i = 0; i < nCount; ++i)
    {
        values.emplace_back(GDALDataset::FromHandle(pahDS[i]));
    }
    return hArg->ptr->Set(std::move(values));
}

/************************************************************************/
/*                  GDALAlgorithmArgSetDatasetNames()                   */
/************************************************************************/

/** Set dataset names to a GAAT_DATASET_LIST argument.
 *
 * @param hArg Handle to an argument. Must NOT be null.
 * @param names Dataset names as a NULL terminated list (may be null)
 * @return true if success.
 * @since 3.11
 */

bool GDALAlgorithmArgSetDatasetNames(GDALAlgorithmArgH hArg, CSLConstList names)
{
    VALIDATE_POINTER1(hArg, __func__, false);
    std::vector<GDALArgDatasetValue> values;
    for (size_t i = 0; names[i]; ++i)
    {
        values.emplace_back(names[i]);
    }
    return hArg->ptr->Set(std::move(values));
}

/************************************************************************/
/*                     GDALArgDatasetValueCreate()                      */
/************************************************************************/

/** Instantiate an empty GDALArgDatasetValue
 *
 * @return new handle to free with GDALArgDatasetValueRelease()
 * @since 3.11
 */
GDALArgDatasetValueH GDALArgDatasetValueCreate()
{
    return std::make_unique<GDALArgDatasetValueHS>().release();
}

/************************************************************************/
/*                     GDALArgDatasetValueRelease()                     */
/************************************************************************/

/** Release a handle to a GDALArgDatasetValue
 *
 * @since 3.11
 */
void GDALArgDatasetValueRelease(GDALArgDatasetValueH hValue)
{
    delete hValue;
}

/************************************************************************/
/*                     GDALArgDatasetValueGetName()                     */
/************************************************************************/

/** Return the name component of the GDALArgDatasetValue
 *
 * @param hValue Handle to a GDALArgDatasetValue. Must NOT be null.
 * @return string whose lifetime is bound to hAlg and which must not
 * be freed.
 * @since 3.11
 */
const char *GDALArgDatasetValueGetName(GDALArgDatasetValueH hValue)
{
    VALIDATE_POINTER1(hValue, __func__, nullptr);
    return hValue->ptr->GetName().c_str();
}

/************************************************************************/
/*                  GDALArgDatasetValueGetDatasetRef()                  */
/************************************************************************/

/** Return the dataset component of the GDALArgDatasetValue.
 *
 * This does not modify the reference counter, hence the lifetime of the
 * returned object is not guaranteed to exceed the one of hValue.
 *
 * @param hValue Handle to a GDALArgDatasetValue. Must NOT be null.
 * @since 3.11
 */
GDALDatasetH GDALArgDatasetValueGetDatasetRef(GDALArgDatasetValueH hValue)
{
    VALIDATE_POINTER1(hValue, __func__, nullptr);
    return GDALDataset::ToHandle(hValue->ptr->GetDatasetRef());
}

/************************************************************************/
/*           GDALArgDatasetValueGetDatasetIncreaseRefCount()            */
/************************************************************************/

/** Return the dataset component of the GDALArgDatasetValue, and increase its
 * reference count if not null. Once done with the dataset, the caller should
 * call GDALReleaseDataset().
 *
 * @param hValue Handle to a GDALArgDatasetValue. Must NOT be null.
 * @since 3.11
 */
GDALDatasetH
GDALArgDatasetValueGetDatasetIncreaseRefCount(GDALArgDatasetValueH hValue)
{
    VALIDATE_POINTER1(hValue, __func__, nullptr);
    return GDALDataset::ToHandle(hValue->ptr->GetDatasetIncreaseRefCount());
}

/************************************************************************/
/*                     GDALArgDatasetValueSetName()                     */
/************************************************************************/

/** Set dataset name
 *
 * @param hValue Handle to a GDALArgDatasetValue. Must NOT be null.
 * @param pszName Dataset name. May be null.
 * @since 3.11
 */

void GDALArgDatasetValueSetName(GDALArgDatasetValueH hValue,
                                const char *pszName)
{
    VALIDATE_POINTER0(hValue, __func__);
    hValue->ptr->Set(pszName ? pszName : "");
}

/************************************************************************/
/*                   GDALArgDatasetValueSetDataset()                    */
/************************************************************************/

/** Set dataset object, increasing its reference counter.
 *
 * @param hValue Handle to a GDALArgDatasetValue. Must NOT be null.
 * @param hDS Dataset object. May be null.
 * @since 3.11
 */

void GDALArgDatasetValueSetDataset(GDALArgDatasetValueH hValue,
                                   GDALDatasetH hDS)
{
    VALIDATE_POINTER0(hValue, __func__);
    hValue->ptr->Set(GDALDataset::FromHandle(hDS));
}
