In the vario_estimate function, is there any way to calculate cross variogram:
$$\gamma_{\ell m}(r_k) = \frac{1}{2 N(rk)}\sum{i=1}^{N(r_k)}(z^{(\ell)}(\mathbf{x}_i)-z^{(m)}(\mathbf{x}_i'))^2$$
or
$$\gamma_{\ell m}(r_k) = \frac{1}{2 N(rk)}\sum{i=1}^{N(r_k)}(z^{(\ell)}(\mathbf{x}_i)-z^{(\ell)}(\mathbf{x}_i'))(z^{(m)}(\mathbf{x}_i)-z^{(m)}(\mathbf{x}_i'))$$
instead of semivarogram when two fields $z^{(\ell)}$ and $z^{(m)}$ are given as a list of numpy.ndarray?
Sorry for the newbie question.
In the
vario_estimate
function, is there any way to calculate cross variogram: $$\gamma_{\ell m}(r_k) = \frac{1}{2 N(rk)}\sum{i=1}^{N(r_k)}(z^{(\ell)}(\mathbf{x}_i)-z^{(m)}(\mathbf{x}_i'))^2$$ or$$\gamma_{\ell m}(r_k) = \frac{1}{2 N(rk)}\sum{i=1}^{N(r_k)}(z^{(\ell)}(\mathbf{x}_i)-z^{(\ell)}(\mathbf{x}_i'))(z^{(m)}(\mathbf{x}_i)-z^{(m)}(\mathbf{x}_i'))$$ instead of semivarogram when two fields $z^{(\ell)}$ and $z^{(m)}$ are given as a list of numpy.ndarray?