· Posts  · 1 min read

Variography made easy with Python

This browser does not support PDFs. Please download the PDF to view it: Download PDF.

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
)

image

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
)

image

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"
        }
    ]
)

image

View on LinkedIn

Back to Blog