~ 1 min read

Variography made easy with Python

Our intern made our clients save time with GeoLime. Should we hire him?

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

Share: