1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
|
/******************************************************************************
*
* Project: GDAL
* Purpose: "gdal vector simplify"
* Author: Even Rouault <even dot rouault at spatialys.com>
*
******************************************************************************
* Copyright (c) 2025, Even Rouault <even dot rouault at spatialys.com>
*
* SPDX-License-Identifier: MIT
****************************************************************************/
#include "gdalalg_vector_simplify.h"
#include "gdal_priv.h"
#include "ogrsf_frmts.h"
//! @cond Doxygen_Suppress
#ifndef _
#define _(x) (x)
#endif
/************************************************************************/
/* GDALVectorSimplifyAlgorithm() */
/************************************************************************/
GDALVectorSimplifyAlgorithm::GDALVectorSimplifyAlgorithm(bool standaloneStep)
: GDALVectorGeomAbstractAlgorithm(NAME, DESCRIPTION, HELP_URL,
standaloneStep, m_opts)
{
AddArg("tolerance", 0, _("Distance tolerance for simplification."),
&m_opts.m_tolerance)
.SetPositional()
.SetRequired()
.SetMinValueIncluded(0);
}
#ifdef HAVE_GEOS
namespace
{
/************************************************************************/
/* GDALVectorSimplifyAlgorithmLayer */
/************************************************************************/
class GDALVectorSimplifyAlgorithmLayer final
: public GDALVectorGeomOneToOneAlgorithmLayer<GDALVectorSimplifyAlgorithm>
{
protected:
using GDALVectorGeomOneToOneAlgorithmLayer::TranslateFeature;
std::unique_ptr<OGRFeature>
TranslateFeature(std::unique_ptr<OGRFeature> poSrcFeature) const override;
public:
GDALVectorSimplifyAlgorithmLayer(
OGRLayer &oSrcLayer, const GDALVectorSimplifyAlgorithm::Options &opts)
: GDALVectorGeomOneToOneAlgorithmLayer<GDALVectorSimplifyAlgorithm>(
oSrcLayer, opts)
{
}
};
/************************************************************************/
/* TranslateFeature() */
/************************************************************************/
std::unique_ptr<OGRFeature> GDALVectorSimplifyAlgorithmLayer::TranslateFeature(
std::unique_ptr<OGRFeature> poSrcFeature) const
{
const int nGeomFieldCount = poSrcFeature->GetGeomFieldCount();
for (int i = 0; i < nGeomFieldCount; ++i)
{
if (IsSelectedGeomField(i))
{
if (auto poGeom = std::unique_ptr<OGRGeometry>(
poSrcFeature->StealGeometry(i)))
{
poGeom.reset(
poGeom->SimplifyPreserveTopology(m_opts.m_tolerance));
if (poGeom)
{
poGeom->assignSpatialReference(m_srcLayer.GetLayerDefn()
->GetGeomFieldDefn(i)
->GetSpatialRef());
poSrcFeature->SetGeomField(i, std::move(poGeom));
}
}
}
}
return poSrcFeature;
}
} // namespace
#endif // HAVE_GEOS
/************************************************************************/
/* GDALVectorSimplifyAlgorithm::CreateAlgLayer() */
/************************************************************************/
std::unique_ptr<OGRLayerWithTranslateFeature>
GDALVectorSimplifyAlgorithm::CreateAlgLayer([[maybe_unused]] OGRLayer &srcLayer)
{
#ifdef HAVE_GEOS
return std::make_unique<GDALVectorSimplifyAlgorithmLayer>(srcLayer, m_opts);
#else
CPLAssert(false);
return nullptr;
#endif
}
/************************************************************************/
/* GDALVectorSimplifyAlgorithm::RunStep() */
/************************************************************************/
bool GDALVectorSimplifyAlgorithm::RunStep(GDALPipelineStepRunContext &ctxt)
{
#ifdef HAVE_GEOS
return GDALVectorGeomAbstractAlgorithm::RunStep(ctxt);
#else
(void)ctxt;
ReportError(CE_Failure, CPLE_NotSupported,
"This algorithm is only supported for builds against GEOS");
return false;
#endif
}
GDALVectorSimplifyAlgorithmStandalone::
~GDALVectorSimplifyAlgorithmStandalone() = default;
//! @endcond
|
