GeoStat-Framework / GSTools

GSTools - A geostatistical toolbox: random fields, variogram estimation, covariance models, kriging and much more
https://geostat-framework.org
GNU Lesser General Public License v3.0
567 stars 74 forks source link

What is the difference between "direction" and "angles" pars in vario_estimate() function? #265

Closed flydream0428 closed 2 years ago

flydream0428 commented 2 years ago

I have tried to estimate a directional variogram by the 'direction' and 'angles' parameters. It turns out it gives very similar results. For 'direction', I am using,

bin_center, dir_vario_ani = gs.vario_estimate((x_con, y_con), z_con, direction=gs.rotated_main_axes(dim=2, angles=1/8*np.pi))

For 'angles', I am using, bin_center, dir_vario_ani = gs.vario_estimate((x_con, y_con), z_con, angles=[1/8*np.pi, 5/8*np.pi]

What is the right one to use here? I have tried to read the function docstring, but I am still confused.