33 #if defined(ATVu2v0_RECURSES) 34 #error Recursive header files inclusion detected in ATVu2v0.h 35 #else // defined(ATVu2v0_RECURSES) 37 #define ATVu2v0_RECURSES 39 #if !defined ATVu2v0_h 47 #include "DGtal/math/linalg/EigenSupport.h" 48 #include "DGtal/base/Common.h" 49 #include "DGtal/dec/DiscreteExteriorCalculus.h" 50 #include "DGtal/dec/DiscreteExteriorCalculusSolver.h" 51 #include "DGtal/dec/DiscreteExteriorCalculusFactory.h" 52 #include "DECImageHelpers.h" 75 template <
typename TKSpace,
76 typename TLinearAlgebra = EigenLinearAlgebraBackend >
204 template <
typename Image,
typename Function>
207 bool perfect_data =
false );
364 std::vector< PrimalForm2 >
g2;
366 std::vector< PrimalForm2 >
i2;
368 std::vector< PrimalForm2 >
u2;
419 template <
typename TKSpace,
typename TLinearAlgebra>
421 operator<< ( std::ostream & out, const ATVu2v0<TKSpace, TLinearAlgebra> & object );
429 #include "ATVu2v0.ih" 434 #endif // !defined ATVu2v0_h 436 #undef ATVu2v0_RECURSES 437 #endif // else defined(ATVu2v0_RECURSES) void setLambda(Scalar _lambda)
Calculus::PrimalHodge0 PrimalHodge0
PrimalIdentity0 left_V0
(lambda / (4*epsilon)) * Id0 + lambda epsilon D0^t D0
PrimalForm0 l_over_e_1
lambda/epsilon 1
void setEpsilon(Scalar _epsilon)
PrimalIdentity0 primal_L0
primal laplace operator on 0-forms
Calculus::DualHodge0 DualHodge0
std::vector< PrimalForm2 > i2
The ideal input 2-forms (for snr computation).
const PrimalForm2 & getG(int i) const
Calculus::PrimalDerivative0 PrimalDerivative0
double alpha
Smoothness parameter alpha of AT (in 1/area unit)
DualDerivative0 dual_D0
dual derivative dual 0-form -> dual 1-form
const PrimalForm2 & getU(int i) const
Calculus::DualDerivative1 DualDerivative1
double lambda
Amount of discontinuity parameter lambda (in 1/length unit).
DECImage2D< TKSpace, TLinearAlgebra > Base
Calculus::PrimalAntiderivative1 PrimalAntiderivative1
std::vector< PrimalForm2 > alpha_g2
alpha g2
PrimalForm0 v0
The v 0-form.
Aim: This class simplifies the development of 2D image processing tools using discrete exterior calcu...
Calculus::PrimalHodge1 PrimalHodge1
SolverV solver_v
The solver for 0-form v.
PrimalAntiderivative2 primal_AD2
Antiderivative 2-form -> to 1-form.
Scalar computeSNR() const
Computes the SNR of u wrt ideal input (should have been given.
DiscreteExteriorCalculusSolver< Calculus, LinearAlgebraSolver, 0, PRIMAL, 0, PRIMAL > SolverV
PrimalIdentity0 l_over_e_Id0
lambda/epsilon Id0
PrimalHodge2 primal_h2
hodge star: 2-form -> dual 2-form
DualDerivative1 dual_D1
dual derivative dual 1-form -> dual 2-form
Aim: This class solves a variant of Ambrosio-Tortorelli functional in a plane for u a (vector of) 2-f...
PrimalForm0 former_v0
The v 0-form at the previous iteration.
Scalar delta_v_loo
The Linfinity-norm of variation of v.
DualHodge0 dual_h0
hodge star: dual 0-form -> 0-form
bool metric_average
When 'true', use metric average, otherwise use identity.
void selfDisplay(std::ostream &out) const
Calculus::PrimalHodge2 PrimalHodge2
std::vector< PrimalForm2 > u2
The u 2-forms.
double epsilon
Thickness of discontinuity set (in length unit).
Calculus::DualAntiderivative2 DualAntiderivative2
Calculus::PrimalForm2 PrimalForm2
Calculus::PrimalIdentity2 PrimalIdentity2
ATVu2v0 & operator=(const ATVu2v0 &other)=delete
LinearAlgebra::SolverSimplicialLLT LinearAlgebraSolver
Calculus::DualAntiderivative1 DualAntiderivative1
TLinearAlgebra LinearAlgebra
Space::RealVector RealVector
Calculus::PrimalForm1 PrimalForm1
DiscreteExteriorCalculus< 2, 2, LinearAlgebra > Calculus
DualHodge1 dual_h1
hodge star: dual 1-form -> 1-form
Calculus::DualHodge2 DualHodge2
Calculus::PrimalIdentity1 PrimalIdentity1
Scalar delta_v_l1
The L1-norm of variation of v.
PrimalDerivative1 M12
edge_to_face average operator
Calculus::PrimalDerivative1 PrimalDerivative1
const PrimalForm0 & getV() const
Calculus::DualHodge1 DualHodge1
BOOST_STATIC_ASSERT((KSpace::dimension==2))
int verbose
The verbose level (0: silent).
void init(Clone< KSpace > K)
Scalar delta_v_l2
The L2-norm of variation of v.
Calculus::PrimalIdentity0 PrimalIdentity0
SolverU solver_u
The solver for every 2-form u[i].
Calculus calculus
The discrete exterior calculus instance.
void setUFromInputAndMask()
PrimalHodge0 primal_h0
hodge star: 0-form -> dual 0-form
Calculus::DualDerivative0 DualDerivative0
Calculus::PrimalAntiderivative2 PrimalAntiderivative2
DualHodge2 dual_h2
hodge star: dual 2-form -> 2-form
unsigned int size0() const
void setAlpha(Scalar _alpha)
PrimalDerivative0 D0
primal derivative: 0-form -> 1-form
std::vector< PrimalForm2 > g2
The g 2-forms.
Scalar getEpsilon() const
DiscreteExteriorCalculusSolver< Calculus, LinearAlgebraSolver, 2, PRIMAL, 2, PRIMAL > SolverU
static const constexpr Dimension dimension
void addInput(const Image &image, const Function &f, bool perfect_data=false)
PrimalDerivative1 D1
primal derivative: 1-form -> 2-form
PrimalHodge1 primal_h1
hodge star: 1-form -> dual 1-form
Scalar computeVariation()
void setMetricAverage(bool average)
Use metric average to smooth L1-metric effects.
Calculus::PrimalForm0 PrimalForm0
HyperRectDomain< Space > Domain
PrimalIdentity2 alpha_Id2
alpha Id2
RealVector::Component Scalar
PrimalDerivative0 M01
point_to_edge average operator