43 std::size_t nLags = matrix.size1();
46 boost::numeric::ublas::matrix<double> m(nLags, 2);
49 std::vector<double> exponentialValues;
51 for(std::size_t t= 0; t < nLags; ++t)
53 if(matrix(t, 0) == 0.)
55 exponentialValues.push_back(m_nugget);
59 exponentialValues.push_back(m_nugget + (m_sill - m_nugget) * (1 - exp(- (3. * matrix(t, 0)) / m_range)));
62 m(t, 0) = matrix(t, 0);
63 m(t, 1) = exponentialValues[t];
The empirical variogram cannot be computed at every lag distance h and due to variation in the estima...
virtual boost::numeric::ublas::matrix< double > calculate(boost::numeric::ublas::matrix< double > matrix)
Function to calculate the geostatistical information from model.
virtual ~GeostatisticalModelExponential()
Virtual destructor.
GeostatisticalModelExponential()
Default constructor.
This file contains a class that represents the geostatistical exponential model.
te::sa::GeostatisticalModelType m_type
Attribute used to define the model type used.