· Posts · 1 min read
Variography made easy with Python
lags, tol = geo.generate_lags(
lags=20,
plag=50,
nlags=10
)
geo.vario_contour(
composite,
attribute='Fe203',
region='HighGradeZone',
lags=lags,
tol=tol,
n_az=20,
atol=20
)
vario_exp_160_major = geo.variogram(
composite,
attribute='Fe203',
region='HighGradeZone',
geographic_azimuth=160,
dip=0,
pitch=0,
lags=lags,
tol=tol,
atol=20
)
vario_exp_160_minor = geo.variogram(
composite,
attribute='Fe203',
region='HighGradeZone',
geographic_azimuth=70,
dip=0,
pitch=0,
lags=lags,
tol=tol,
atol=45
)
geo.plot_semivariogram(
variograms=[
vario_exp_160_major,
vario_exp_160_minor
],
diplay_npairs=True
)
cov_model = geo.Nugget() + geo.Spherical()
geo.model_fit(
[
vario_exp_160_major,
vario_exp_160_minor
],
cov_model
)
geo.plot_semivariogram(
variograms=[
vario_exp_160_major,
vario_exp_160_minor
],
model=cov_model,
model_angles=[
{
"azi":160,
"dip":0,
"pitch":0,
"label":"N160"
},
{
"azi":70,
"dip":0,
"pitch":0,
"label":"N070"
}
]
)