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
|
/******************************************************************************
*
* Project: GDAL
* Purpose: gdal "vector grid invdist" subcommand
* 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_grid_invdist.h"
#include <limits>
//! @cond Doxygen_Suppress
#ifndef _
#define _(x) (x)
#endif
/************************************************************************/
/* GDALVectorGridInvdistAlgorithm::GDALVectorGridInvdistAlgorithm() */
/************************************************************************/
GDALVectorGridInvdistAlgorithm::GDALVectorGridInvdistAlgorithm(
bool standaloneStep)
: GDALVectorGridAbstractAlgorithm(NAME, DESCRIPTION, HELP_URL,
standaloneStep)
{
AddArg("power", 0, _("Weighting power"), &m_power).SetDefault(m_power);
AddArg("smoothing", 0, _("Smoothing parameter"), &m_smoothing)
.SetDefault(m_smoothing);
AddRadiusArg();
AddRadius1AndRadius2Arg();
AddAngleArg();
AddMinPointsArg();
AddMaxPointsArg();
AddMinMaxPointsPerQuadrantArg();
AddNodataArg();
AddValidationAction(
[this]()
{
bool ret = true;
if (m_minPoints > 0 && m_radius == 0 && m_radius1 == 0)
{
ReportError(CE_Failure, CPLE_AppDefined,
"'radius' or 'radius1' and 'radius2' should be "
"defined when 'min-points' is.");
ret = false;
}
if (m_maxPoints < std::numeric_limits<int>::max() &&
m_radius == 0 && m_radius1 == 0)
{
ReportError(CE_Failure, CPLE_AppDefined,
"'radius' or 'radius1' and 'radius2' should be "
"defined when 'max-points' is.");
ret = false;
}
return ret;
});
}
/************************************************************************/
/* GDALVectorGridInvdistAlgorithm::RunImpl() */
/************************************************************************/
std::string GDALVectorGridInvdistAlgorithm::GetGridAlgorithm() const
{
std::string ret = CPLSPrintf(
"invdist:power=%.17g:smoothing=%.17g:angle=%.17g:nodata=%.17g", m_power,
m_smoothing, m_angle, m_nodata);
if (m_radius > 0)
{
ret += CPLSPrintf(":radius=%.17g", m_radius);
}
else
{
if (m_radius1 > 0)
ret += CPLSPrintf(":radius1=%.17g", m_radius1);
if (m_radius2 > 0)
ret += CPLSPrintf(":radius2=%.17g", m_radius2);
}
if (m_minPoints > 0)
ret += CPLSPrintf(":min_points=%d", m_minPoints);
if (m_maxPoints < std::numeric_limits<int>::max())
ret += CPLSPrintf(":max_points=%d", m_maxPoints);
if (m_minPointsPerQuadrant > 0)
ret +=
CPLSPrintf(":min_points_per_quadrant=%d", m_minPointsPerQuadrant);
if (m_maxPointsPerQuadrant < std::numeric_limits<int>::max())
ret +=
CPLSPrintf(":max_points_per_quadrant=%d", m_maxPointsPerQuadrant);
return ret;
}
GDALVectorGridInvdistAlgorithmStandalone::
~GDALVectorGridInvdistAlgorithmStandalone() = default;
//! @endcond
|
