Function variogram_directional

pub fn variogram_directional(
    f: ArrayView2<'_, f64>,
    bin_edges: ArrayView1<'_, f64>,
    pos: ArrayView2<'_, f64>,
    direction: ArrayView2<'_, f64>,
    angles_tol: f64,
    bandwidth: f64,
    separate_dirs: bool,
    estimator_type: char,
    num_threads: Option<usize>,
) -> (Array2<f64>, Array2<u64>)

Directional variogram estimation on an unstructured grid.

Calculates the empirical variogram according to the equations shown in the module documentation.

Arguments

• f - the spatially distributed data
       dim = (no. of data fields, no. of spatial data points per field $i$)
• bin_edges - the bins of the variogram
       dim = number of bins j
• pos - the positions of the data f
       dim = (spatial dim. $d$, no. of spatial data points $i$)
• direction - directions in which the variogram will be estimated
       dim = (no. of directions, spatial dim. $d$)
• angles_tol - the tolerance of the angles
• bandwidth - bandwidth to cut off the angular tolerance
• separate_dirs - do the direction bands overlap?
• estimator_type - the estimator function, can be
  • ‘m’ - Matheron, the standard method of moments by Matheron
  • ‘c’ - Cressie, an estimator more robust to outliers
• num_threads - the number of parallel threads used, if None, use rayon’s default