43 std::size_t nLags = matrix.size1();
46 boost::numeric::ublas::matrix<double> m(nLags, 2);
49 std::vector<double> sphericalValues;
51 for(std::size_t t= 0; t < nLags; ++t)
53 if(matrix(t, 0) == 0.)
55 sphericalValues.push_back(m_nugget);
59 sphericalValues.push_back(m_nugget + (m_sill - m_nugget) * (1.5 * (std::abs(matrix(t, 0)) / m_range)) - 0.5 * (pow(std::abs(matrix(t, 0)) / m_range, 3)));
61 m(t, 0) = matrix(t, 0);
62 m(t, 1) = sphericalValues[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.
This file contains a class that represents the geostatistical spherical model.
GeostatisticalModelSpherical()
Default constructor.
virtual ~GeostatisticalModelSpherical()
Virtual destructor.
te::sa::GeostatisticalModelType m_type
Attribute used to define the model type used.