Example of DOV search methods for interpretations (gecodeerde lithologie)

Use cases explained below

• Get ‘gecodeerde lithologie’ in a bounding box

• Get ‘gecodeerde lithologie’ with specific properties within a distance from a point

• Get ‘gecodeerde lithologie’ in a bounding box with specific properties

• Get ‘gecodeerde lithologie’ based on fields not available in the standard output dataframe

• Get ‘gecodeerde lithologie’ data, returning fields not available in the standard output dataframe

%matplotlib inline
import inspect, sys

# check pydov path
import pydov

Get information about the datatype ‘Gecodeerde lithologie’

from pydov.search.interpretaties import GecodeerdeLithologieSearch
itp = GecodeerdeLithologieSearch()

A description is provided for the ‘Gecodeerde lithologie’ datatype:

itp.get_description()

'Een gecodeerde lithologie van een boring is een geologische codering van een lithologische beschrijving op basis van een vaste DOV-standaard per diepte-interval. U vindt een overzicht van de gebruikte standaarden op onze website.'

The different fields that are available for objects of the ‘Gecodeerde lithologie’ datatype can be requested with the get_fields() method:

fields = itp.get_fields()

# print available fields
for f in fields.values():
    print(f['name'])

pkey_interpretatie
Type_proef
Proefnummer
pkey_boring
x
y
start_interpretatie_mtaw
diepte_tot_m
gemeente
Auteurs
Datum
Opdrachten
betrouwbaarheid_interpretatie
Geldig_van
Geldig_tot
Score_categorie
Geschiktheid_boormethode
eerste_invoer
geom
diepte_laag_van
diepte_laag_tot
hoofdnaam1_grondsoort
hoofdnaam2_grondsoort
bijmenging1_plaatselijk
bijmenging1_hoeveelheid
bijmenging1_grondsoort
bijmenging2_plaatselijk
bijmenging2_hoeveelheid
bijmenging2_grondsoort
bijmenging3_plaatselijk
bijmenging3_hoeveelheid
bijmenging3_grondsoort

Alternatively, you can list all the fields and their details by inspecting the get_fields() output or the search instance itself in a notebook:

itp

pydov.search.interpretaties.GecodeerdeLithologieSearch

Een gecodeerde lithologie van een boring is een geologische codering van een lithologische beschrijving op basis van een vaste DOV-standaard per diepte-interval. U vindt een overzicht van de gebruikte standaarden op onze website.

pkey_interpretatie - URL die verwijst naar de gegevens van deze gecodeerde lithologie op de website. Voeg '.xml' toe om een XML voorstelling van deze gegevens te verkrijgen.

• type: string
• notnull: False
• query: True
• cost: 1
• multivalue: False

Type_proef - Het type van de proef waarvan de gecodeerde lithologie een interpretatie is. Bij een gecodeerde lithologie zal dit altijd 'Boring' zijn.

• type: string
• notnull: False
• query: True
• cost: 1
• multivalue: False
• codelist:
pydov.util.codelists.FeatureCatalogueValues

Boring - Boring - De gecodeerde lithologie is een interpretatie van een boring.

Proefnummer - Het proefnummer van de boring waarvan de gecodeerde lithologie een interpretatie is.

• type: string
• notnull: False
• query: True
• cost: 1
• multivalue: False

pkey_boring - URL die verwijst naar de gegevens van de boring op de website. Voeg '.xml' toe om een XML voorstelling van deze gegevens te verkrijgen.

• type: string
• notnull: False
• query: True
• cost: 1
• multivalue: False

x - De x-coördinaat van de boring in het Lambert72 coördinaatsysteem (in meter, EPSG:31370).

• type: float
• notnull: False
• query: True
• cost: 1
• multivalue: False

y - De y-coördinaat van de boring in het Lambert72 coördinaatsysteem (in meter, EPSG:31370).

• type: float
• notnull: False
• query: True
• cost: 1
• multivalue: False

start_interpretatie_mtaw - De hoogte van het aanvangspeil van de boring in het TAW stelsel (in meter).

• type: float
• notnull: False
• query: True
• cost: 1
• multivalue: False

diepte_tot_m - De maximumdiepte van de gecodeerde lithologie in meter. Dit is de 'tot' van de diepste geïnterpreteerde laag.

• type: float
• notnull: False
• query: True
• cost: 1
• multivalue: False

gemeente - De gemeente waarin de boring van deze interpretatie gelegen is.

• type: string
• notnull: False
• query: True
• cost: 1
• multivalue: False

Auteurs - De auteur of auteurs van de gecodeerde lithologie.

• type: string
• notnull: False
• query: True
• cost: 1
• multivalue: False

Datum - De datum waarop de gecodeerde lithologie werd beschreven.

• type: date
• notnull: False
• query: True
• cost: 1
• multivalue: False

Opdrachten - De opdracht(en) waaraan de gecodeerde lithologie gekoppeld is.

• type: string
• notnull: False
• query: True
• cost: 1
• multivalue: False

betrouwbaarheid_interpretatie - De betrouwbaarheid van het geheel van de gecodeerde lithologie (goed, twijfelachtig of onbekend). Deze wordt toegekend door een DOV-medewerker verantwoordelijk voor de ingave.

• type: string
• notnull: False
• query: True
• cost: 1
• multivalue: False

Geldig_van - Het tijdstip vanaf wanneer de gecodeerde lithologie geldig is.

• type: date
• notnull: False
• query: True
• cost: 1
• multivalue: False

Geldig_tot - Het tijdstip tot wanneer de gecodeerde lithologie geldig is.

• type: date
• notnull: False
• query: True
• cost: 1
• multivalue: False

Score_categorie - De score categorie is een kwalitatieve categorie op basis van een kwantitatieve score, tussen 0 en 10, van de kwaliteit van (gecodeerde) lithologische beschrijvingen. Slecht: < 4,7/10, matig: tussen 4,7/10 en 6,6/10, goed: tussen 6,6/10 en 8,5/10, zeer goed: > 8,5/10. De score is samengesteld op basis van de beschrijving van basiselementen (lithologie, kleur, korrelgrootte), voldoende detail (secundaire kenmerken zoals glauconiet, schelpen, …) en voldoende kleine intervallen. Voor goed hergebruik van de beschrijvingen is minstens categorie ‘matig’ nodig.

• type: string
• notnull: False
• query: True
• cost: 1
• multivalue: False

Geschiktheid_boormethode - De geschiktheid van de boormethode beschrijft of de gebruikte boormethode van de boring waarop de beschrijving gebaseerd is eerder geschikt (hoog; bv. kernboring) of eerder niet geschikt (laag; bv. spoelboring) is voor het opmaken van een gedetailleerde (gecodeerde) lithologische beschrijving. Boormethodes die ongeroerde stalen opleveren zijn doorgaans geschikter dan methodes die geroerde stalen opleveren. Voor beschrijvingen met een hoge geschiktheid van de boormethode mag ook een betere score categorie verwacht worden dan van beschrijvingen met een lage geschiktheid van de boormethode.

• type: string
• notnull: False
• query: True
• cost: 1
• multivalue: False

eerste_invoer - Het tijdstip waarop deze interpretatie voor het eerst in DOV ingevoerd werd.

• type: datetime
• notnull: True
• query: True
• cost: 1
• multivalue: False

geom - None

• type: geometry
• notnull: False
• query: False
• cost: 1
• multivalue: False

diepte_laag_van - Diepte van de bovenkant van de laag gecodeerde lithologie in meter.

• type: float
• notnull: False
• query: False
• cost: 10
• multivalue: False

diepte_laag_tot - Diepte van de onderkant van de laag gecodeerde lithologie in meter.

• type: float
• notnull: False
• query: False
• cost: 10
• multivalue: False

hoofdnaam1_grondsoort - Primaire grondsoort (als code) van de laag gecodeerde lithologie

• type: string
• notnull: False
• query: False
• cost: 10
• multivalue: False
• codelist:
pydov.util.codelists.XsdType

AA - AA - amfiboliet

AB - AB - andesiet

AC - AC - anhydriet

AD - AD - ankeriet

AE - AE - arkose

AF - AF - afval

AG - AG - as

AN - AN - antropogeen niet vloeistofdicht

AR - AR - areniet

AT - AT - asfalt

AV - AV - antropogeen vloeistofdicht

AW - AW - antropogeen weinig vloeistofdicht

BA - BA - basalt

BB - BB - bauxiet

BC - BC - bentoniet

BK - BK - baksteen

BS - BS - blauwshist

BT - BT - beton

CA - CA - carbonaatfragmenten

CB - CB - conglomeraat/breccie

CO - CO - steenkool

CT - CT - chert (gebande)

DA - DA - daciet

DB - DB - diatomiet

DC - DC - dioriet

DD - DD - dolerriet

DE - DE - duniet

DK - DK - dolom kalksteen

DL - DL - dolomiet

EA - EA - evaporiet

EC - EC - eclogiet

FA - FA - fosforiet

FE - FE - ijzer

FG - FG - fijn grind

FO - FO - fossielen

FY - FY - fylliet

FZ - FZ - fijn zand

GA - GA - geroerd / gestoord

GB - GB - gabbro

GC - GC - gips

GD - GD - gneiss

GE - GE - graniet

GF - GF - granodioriet

GG - GG - grof grind

GH - GH - groenshist

GI - GI - glimmers

GL - GL - glauconiet

GM - GM - gespleten metamorf

GS - GS - groensteen

GZ - GZ - grof zand

HA - HA - haliet

HB - HB - hematiet

HC - HC - houtskool

HD - HD - hornfels

HT - HT - houtresten

HU - HU - humus

IA - IA - ijzerooliet

IN - IN - intrusief gesteente

KA - KA - kalk

KB - KB - kalksteenbreccie

KC - KC - kalkschiefer

KL - KL - klei

KO - KO - komatiet

KR - KR - krijt

KS - KS - kalksteen

KT - KT - kleisteen

KW - KW - kwartsiet

LA - LA - lateriet

LE - LE - leem

LG - LG - ligniet

LI - LI - limoniet

LS - LS - leisteen

MA - MA - marmer

MB - MB - micashist

ME - ME - mergel

MG - MG - middelmatig grind

MI - MI - migmatiet

MM - MM - metamorf gesteente

MS - MS - schiefer

MY - MY - myloniet

MZ - MZ - middelmatig zand

NN - NN - nummulieten

OK - OK - oölitische kalksteen

OM - OM - organisch materiaal

PA - PA - pegmatiet

PB - PB - peridodiet

PC - PC - poederkool

PD - PD - pyroxeniet

PL - PL - plantenresten

PM - PM - puimsteen

PS - PS - psammiet

PU - PU - puin

PY - PY - pyriet

QU - QU - kwarts

RA - RA - radiolariet

RB - RB - rhyoliet

RO - RO - rots

RV - RV - verweerde rots

SA - SA - sideriet

SC - SC - schelpen

SE - SE - sedimentair

SF - SF - schelpfragmenten

SG - SG - schelpgruis

SH - SH - schalie

SI - SI - silt

SK - SK - schieferige kalkst.

SL - SL - slib

SN - SN - steenfragmenten

SO - SO - scoria

SP - SP - serpentiniet

SS - SS - siltsteen

ST - ST - stenen

SX - SX - silex

SY - SY - shist

TA - TA - teelaarde

TB - TB - tuff

TS - TS - tufsteen (tuffeau)

VA - VA - vulkanisch gesteente

VE - VE - veen

VS - VS - vuursteen

WA - WA - wacksteen

XG - XG - grind

XX - XX - onbekend

XZ - XZ - zand

ZE - ZE - zeepsteen

ZK - ZK - kalkzandsteen

ZL - ZL - zandleem

ZS - ZS - zandsteen

hoofdnaam2_grondsoort - Secundaire grondsoort (als code) van de laag gecodeerde lithologie

• type: string
• notnull: False
• query: False
• cost: 10
• multivalue: False
• codelist:
pydov.util.codelists.XsdType

AA - AA - amfiboliet

AB - AB - andesiet

AC - AC - anhydriet

AD - AD - ankeriet

AE - AE - arkose

AF - AF - afval

AG - AG - as

AN - AN - antropogeen niet vloeistofdicht

AR - AR - areniet

AT - AT - asfalt

AV - AV - antropogeen vloeistofdicht

AW - AW - antropogeen weinig vloeistofdicht

BA - BA - basalt

BB - BB - bauxiet

BC - BC - bentoniet

BK - BK - baksteen

BS - BS - blauwshist

BT - BT - beton

CA - CA - carbonaatfragmenten

CB - CB - conglomeraat/breccie

CO - CO - steenkool

CT - CT - chert (gebande)

DA - DA - daciet

DB - DB - diatomiet

DC - DC - dioriet

DD - DD - dolerriet

DE - DE - duniet

DK - DK - dolom kalksteen

DL - DL - dolomiet

EA - EA - evaporiet

EC - EC - eclogiet

FA - FA - fosforiet

FE - FE - ijzer

FG - FG - fijn grind

FO - FO - fossielen

FY - FY - fylliet

FZ - FZ - fijn zand

GA - GA - geroerd / gestoord

GB - GB - gabbro

GC - GC - gips

GD - GD - gneiss

GE - GE - graniet

GF - GF - granodioriet

GG - GG - grof grind

GH - GH - groenshist

GI - GI - glimmers

GL - GL - glauconiet

GM - GM - gespleten metamorf

GS - GS - groensteen

GZ - GZ - grof zand

HA - HA - haliet

HB - HB - hematiet

HC - HC - houtskool

HD - HD - hornfels

HT - HT - houtresten

HU - HU - humus

IA - IA - ijzerooliet

IN - IN - intrusief gesteente

KA - KA - kalk

KB - KB - kalksteenbreccie

KC - KC - kalkschiefer

KL - KL - klei

KO - KO - komatiet

KR - KR - krijt

KS - KS - kalksteen

KT - KT - kleisteen

KW - KW - kwartsiet

LA - LA - lateriet

LE - LE - leem

LG - LG - ligniet

LI - LI - limoniet

LS - LS - leisteen

MA - MA - marmer

MB - MB - micashist

ME - ME - mergel

MG - MG - middelmatig grind

MI - MI - migmatiet

MM - MM - metamorf gesteente

MS - MS - schiefer

MY - MY - myloniet

MZ - MZ - middelmatig zand

NN - NN - nummulieten

OK - OK - oölitische kalksteen

OM - OM - organisch materiaal

PA - PA - pegmatiet

PB - PB - peridodiet

PC - PC - poederkool

PD - PD - pyroxeniet

PL - PL - plantenresten

PM - PM - puimsteen

PS - PS - psammiet

PU - PU - puin

PY - PY - pyriet

QU - QU - kwarts

RA - RA - radiolariet

RB - RB - rhyoliet

RO - RO - rots

RV - RV - verweerde rots

SA - SA - sideriet

SC - SC - schelpen

SE - SE - sedimentair

SF - SF - schelpfragmenten

SG - SG - schelpgruis

SH - SH - schalie

SI - SI - silt

SK - SK - schieferige kalkst.

SL - SL - slib

SN - SN - steenfragmenten

SO - SO - scoria

SP - SP - serpentiniet

SS - SS - siltsteen

ST - ST - stenen

SX - SX - silex

SY - SY - shist

TA - TA - teelaarde

TB - TB - tuff

TS - TS - tufsteen (tuffeau)

VA - VA - vulkanisch gesteente

VE - VE - veen

VS - VS - vuursteen

WA - WA - wacksteen

XG - XG - grind

XX - XX - onbekend

XZ - XZ - zand

ZE - ZE - zeepsteen

ZK - ZK - kalkzandsteen

ZL - ZL - zandleem

ZS - ZS - zandsteen

bijmenging1_plaatselijk - plaatselijk of niet-plaatselijk

• type: boolean
• notnull: False
• query: False
• cost: 10
• multivalue: False

bijmenging1_hoeveelheid - aanduiding van de hoeveelheid bijmenging

• type: string
• notnull: False
• query: False
• cost: 10
• multivalue: False
• codelist:
pydov.util.codelists.XsdType

M - M - Xhoudend

N - N - geen aanduiding van de hoeveelheid

S - S - sterk Xhoudend

V - V - veel X

W - W - weinig X

Z - Z - weinig Xhoudend

bijmenging1_grondsoort - type grondsoort (als code) van de laag gecodeerde lithologie of geotechnische codering

• type: string
• notnull: False
• query: False
• cost: 10
• multivalue: False
• codelist:
pydov.util.codelists.XsdType

AA - AA - amfiboliet

AB - AB - andesiet

AC - AC - anhydriet

AD - AD - ankeriet

AE - AE - arkose

AF - AF - afval

AG - AG - as

AN - AN - antropogeen niet vloeistofdicht

AR - AR - areniet

AT - AT - asfalt

AV - AV - antropogeen vloeistofdicht

AW - AW - antropogeen weinig vloeistofdicht

BA - BA - basalt

BB - BB - bauxiet

BC - BC - bentoniet

BK - BK - baksteen

BS - BS - blauwshist

BT - BT - beton

CA - CA - carbonaatfragmenten

CB - CB - conglomeraat/breccie

CO - CO - steenkool

CT - CT - chert (gebande)

DA - DA - daciet

DB - DB - diatomiet

DC - DC - dioriet

DD - DD - dolerriet

DE - DE - duniet

DK - DK - dolom kalksteen

DL - DL - dolomiet

EA - EA - evaporiet

EC - EC - eclogiet

FA - FA - fosforiet

FE - FE - ijzer

FG - FG - fijn grind

FO - FO - fossielen

FY - FY - fylliet

FZ - FZ - fijn zand

GA - GA - geroerd / gestoord

GB - GB - gabbro

GC - GC - gips

GD - GD - gneiss

GE - GE - graniet

GF - GF - granodioriet

GG - GG - grof grind

GH - GH - groenshist

GI - GI - glimmers

GL - GL - glauconiet

GM - GM - gespleten metamorf

GS - GS - groensteen

GZ - GZ - grof zand

HA - HA - haliet

HB - HB - hematiet

HC - HC - houtskool

HD - HD - hornfels

HT - HT - houtresten

HU - HU - humus

IA - IA - ijzerooliet

IN - IN - intrusief gesteente

KA - KA - kalk

KB - KB - kalksteenbreccie

KC - KC - kalkschiefer

KL - KL - klei

KO - KO - komatiet

KR - KR - krijt

KS - KS - kalksteen

KT - KT - kleisteen

KW - KW - kwartsiet

LA - LA - lateriet

LE - LE - leem

LG - LG - ligniet

LI - LI - limoniet

LS - LS - leisteen

MA - MA - marmer

MB - MB - micashist

ME - ME - mergel

MG - MG - middelmatig grind

MI - MI - migmatiet

MM - MM - metamorf gesteente

MS - MS - schiefer

MY - MY - myloniet

MZ - MZ - middelmatig zand

NN - NN - nummulieten

OK - OK - oölitische kalksteen

OM - OM - organisch materiaal

PA - PA - pegmatiet

PB - PB - peridodiet

PC - PC - poederkool

PD - PD - pyroxeniet

PL - PL - plantenresten

PM - PM - puimsteen

PS - PS - psammiet

PU - PU - puin

PY - PY - pyriet

QU - QU - kwarts

RA - RA - radiolariet

RB - RB - rhyoliet

RO - RO - rots

RV - RV - verweerde rots

SA - SA - sideriet

SC - SC - schelpen

SE - SE - sedimentair

SF - SF - schelpfragmenten

SG - SG - schelpgruis

SH - SH - schalie

SI - SI - silt

SK - SK - schieferige kalkst.

SL - SL - slib

SN - SN - steenfragmenten

SO - SO - scoria

SP - SP - serpentiniet

SS - SS - siltsteen

ST - ST - stenen

SX - SX - silex

SY - SY - shist

TA - TA - teelaarde

TB - TB - tuff

TS - TS - tufsteen (tuffeau)

VA - VA - vulkanisch gesteente

VE - VE - veen

VS - VS - vuursteen

WA - WA - wacksteen

XG - XG - grind

XX - XX - onbekend

XZ - XZ - zand

ZE - ZE - zeepsteen

ZK - ZK - kalkzandsteen

ZL - ZL - zandleem

ZS - ZS - zandsteen

bijmenging2_plaatselijk - plaatselijk of niet-plaatselijk

• type: boolean
• notnull: False
• query: False
• cost: 10
• multivalue: False

bijmenging2_hoeveelheid - aanduiding van de hoeveelheid bijmenging

• type: string
• notnull: False
• query: False
• cost: 10
• multivalue: False
• codelist:
pydov.util.codelists.XsdType

M - M - Xhoudend

N - N - geen aanduiding van de hoeveelheid

S - S - sterk Xhoudend

V - V - veel X

W - W - weinig X

Z - Z - weinig Xhoudend

bijmenging2_grondsoort - type grondsoort (als code) van de laag gecodeerde lithologie of geotechnische codering

• type: string
• notnull: False
• query: False
• cost: 10
• multivalue: False
• codelist:
pydov.util.codelists.XsdType

AA - AA - amfiboliet

AB - AB - andesiet

AC - AC - anhydriet

AD - AD - ankeriet

AE - AE - arkose

AF - AF - afval

AG - AG - as

AN - AN - antropogeen niet vloeistofdicht

AR - AR - areniet

AT - AT - asfalt

AV - AV - antropogeen vloeistofdicht

AW - AW - antropogeen weinig vloeistofdicht

BA - BA - basalt

BB - BB - bauxiet

BC - BC - bentoniet

BK - BK - baksteen

BS - BS - blauwshist

BT - BT - beton

CA - CA - carbonaatfragmenten

CB - CB - conglomeraat/breccie

CO - CO - steenkool

CT - CT - chert (gebande)

DA - DA - daciet

DB - DB - diatomiet

DC - DC - dioriet

DD - DD - dolerriet

DE - DE - duniet

DK - DK - dolom kalksteen

DL - DL - dolomiet

EA - EA - evaporiet

EC - EC - eclogiet

FA - FA - fosforiet

FE - FE - ijzer

FG - FG - fijn grind

FO - FO - fossielen

FY - FY - fylliet

FZ - FZ - fijn zand

GA - GA - geroerd / gestoord

GB - GB - gabbro

GC - GC - gips

GD - GD - gneiss

GE - GE - graniet

GF - GF - granodioriet

GG - GG - grof grind

GH - GH - groenshist

GI - GI - glimmers

GL - GL - glauconiet

GM - GM - gespleten metamorf

GS - GS - groensteen

GZ - GZ - grof zand

HA - HA - haliet

HB - HB - hematiet

HC - HC - houtskool

HD - HD - hornfels

HT - HT - houtresten

HU - HU - humus

IA - IA - ijzerooliet

IN - IN - intrusief gesteente

KA - KA - kalk

KB - KB - kalksteenbreccie

KC - KC - kalkschiefer

KL - KL - klei

KO - KO - komatiet

KR - KR - krijt

KS - KS - kalksteen

KT - KT - kleisteen

KW - KW - kwartsiet

LA - LA - lateriet

LE - LE - leem

LG - LG - ligniet

LI - LI - limoniet

LS - LS - leisteen

MA - MA - marmer

MB - MB - micashist

ME - ME - mergel

MG - MG - middelmatig grind

MI - MI - migmatiet

MM - MM - metamorf gesteente

MS - MS - schiefer

MY - MY - myloniet

MZ - MZ - middelmatig zand

NN - NN - nummulieten

OK - OK - oölitische kalksteen

OM - OM - organisch materiaal

PA - PA - pegmatiet

PB - PB - peridodiet

PC - PC - poederkool

PD - PD - pyroxeniet

PL - PL - plantenresten

PM - PM - puimsteen

PS - PS - psammiet

PU - PU - puin

PY - PY - pyriet

QU - QU - kwarts

RA - RA - radiolariet

RB - RB - rhyoliet

RO - RO - rots

RV - RV - verweerde rots

SA - SA - sideriet

SC - SC - schelpen

SE - SE - sedimentair

SF - SF - schelpfragmenten

SG - SG - schelpgruis

SH - SH - schalie

SI - SI - silt

SK - SK - schieferige kalkst.

SL - SL - slib

SN - SN - steenfragmenten

SO - SO - scoria

SP - SP - serpentiniet

SS - SS - siltsteen

ST - ST - stenen

SX - SX - silex

SY - SY - shist

TA - TA - teelaarde

TB - TB - tuff

TS - TS - tufsteen (tuffeau)

VA - VA - vulkanisch gesteente

VE - VE - veen

VS - VS - vuursteen

WA - WA - wacksteen

XG - XG - grind

XX - XX - onbekend

XZ - XZ - zand

ZE - ZE - zeepsteen

ZK - ZK - kalkzandsteen

ZL - ZL - zandleem

ZS - ZS - zandsteen

bijmenging3_plaatselijk - plaatselijk of niet-plaatselijk

• type: boolean
• notnull: False
• query: False
• cost: 10
• multivalue: False

bijmenging3_hoeveelheid - aanduiding van de hoeveelheid bijmenging

• type: string
• notnull: False
• query: False
• cost: 10
• multivalue: False
• codelist:
pydov.util.codelists.XsdType

M - M - Xhoudend

N - N - geen aanduiding van de hoeveelheid

S - S - sterk Xhoudend

V - V - veel X

W - W - weinig X

Z - Z - weinig Xhoudend

bijmenging3_grondsoort - type grondsoort (als code) van de laag gecodeerde lithologie of geotechnische codering

• type: string
• notnull: False
• query: False
• cost: 10
• multivalue: False
• codelist:
pydov.util.codelists.XsdType

AA - AA - amfiboliet

AB - AB - andesiet

AC - AC - anhydriet

AD - AD - ankeriet

AE - AE - arkose

AF - AF - afval

AG - AG - as

AN - AN - antropogeen niet vloeistofdicht

AR - AR - areniet

AT - AT - asfalt

AV - AV - antropogeen vloeistofdicht

AW - AW - antropogeen weinig vloeistofdicht

BA - BA - basalt

BB - BB - bauxiet

BC - BC - bentoniet

BK - BK - baksteen

BS - BS - blauwshist

BT - BT - beton

CA - CA - carbonaatfragmenten

CB - CB - conglomeraat/breccie

CO - CO - steenkool

CT - CT - chert (gebande)

DA - DA - daciet

DB - DB - diatomiet

DC - DC - dioriet

DD - DD - dolerriet

DE - DE - duniet

DK - DK - dolom kalksteen

DL - DL - dolomiet

EA - EA - evaporiet

EC - EC - eclogiet

FA - FA - fosforiet

FE - FE - ijzer

FG - FG - fijn grind

FO - FO - fossielen

FY - FY - fylliet

FZ - FZ - fijn zand

GA - GA - geroerd / gestoord

GB - GB - gabbro

GC - GC - gips

GD - GD - gneiss

GE - GE - graniet

GF - GF - granodioriet

GG - GG - grof grind

GH - GH - groenshist

GI - GI - glimmers

GL - GL - glauconiet

GM - GM - gespleten metamorf

GS - GS - groensteen

GZ - GZ - grof zand

HA - HA - haliet

HB - HB - hematiet

HC - HC - houtskool

HD - HD - hornfels

HT - HT - houtresten

HU - HU - humus

IA - IA - ijzerooliet

IN - IN - intrusief gesteente

KA - KA - kalk

KB - KB - kalksteenbreccie

KC - KC - kalkschiefer

KL - KL - klei

KO - KO - komatiet

KR - KR - krijt

KS - KS - kalksteen

KT - KT - kleisteen

KW - KW - kwartsiet

LA - LA - lateriet

LE - LE - leem

LG - LG - ligniet

LI - LI - limoniet

LS - LS - leisteen

MA - MA - marmer

MB - MB - micashist

ME - ME - mergel

MG - MG - middelmatig grind

MI - MI - migmatiet

MM - MM - metamorf gesteente

MS - MS - schiefer

MY - MY - myloniet

MZ - MZ - middelmatig zand

NN - NN - nummulieten

OK - OK - oölitische kalksteen

OM - OM - organisch materiaal

PA - PA - pegmatiet

PB - PB - peridodiet

PC - PC - poederkool

PD - PD - pyroxeniet

PL - PL - plantenresten

PM - PM - puimsteen

PS - PS - psammiet

PU - PU - puin

PY - PY - pyriet

QU - QU - kwarts

RA - RA - radiolariet

RB - RB - rhyoliet

RO - RO - rots

RV - RV - verweerde rots

SA - SA - sideriet

SC - SC - schelpen

SE - SE - sedimentair

SF - SF - schelpfragmenten

SG - SG - schelpgruis

SH - SH - schalie

SI - SI - silt

SK - SK - schieferige kalkst.

SL - SL - slib

SN - SN - steenfragmenten

SO - SO - scoria

SP - SP - serpentiniet

SS - SS - siltsteen

ST - ST - stenen

SX - SX - silex

SY - SY - shist

TA - TA - teelaarde

TB - TB - tuff

TS - TS - tufsteen (tuffeau)

VA - VA - vulkanisch gesteente

VE - VE - veen

VS - VS - vuursteen

WA - WA - wacksteen

XG - XG - grind

XX - XX - onbekend

XZ - XZ - zand

ZE - ZE - zeepsteen

ZK - ZK - kalkzandsteen

ZL - ZL - zandleem

ZS - ZS - zandsteen

Example use cases

Get ‘Gecodeerde lithologie’ in a bounding box

Get data for all the ‘Gecodeerde lithologie’ interpretations that are geographically located within the bounds of the specified box.

The coordinates are in the Belgian Lambert72 (EPSG:31370) coordinate system and are given in the order of lower left x, lower left y, upper right x, upper right y.

from pydov.util.location import Within, Box

df = itp.search(location=Within(Box(153145, 206930, 153150, 206935, epsg=31370)))
df.head()

[000/001] .
[000/001] .

	pkey_interpretatie	pkey_boring	betrouwbaarheid_interpretatie	x	y	start_interpretatie_mtaw	diepte_laag_van	diepte_laag_tot	hoofdnaam1_grondsoort	hoofdnaam2_grondsoort	bijmenging1_plaatselijk	bijmenging1_hoeveelheid	bijmenging1_grondsoort	bijmenging2_plaatselijk	bijmenging2_hoeveelheid	bijmenging2_grondsoort	bijmenging3_plaatselijk	bijmenging3_hoeveelheid	bijmenging3_grondsoort
0	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1974...	goed	153147.0	206931.0	14.12	0.0	0.5	LE	NaN	False	N	XZ	NaN	NaN	NaN	NaN	NaN	NaN
1	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1974...	goed	153147.0	206931.0	14.12	0.5	1.0	LE	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
2	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1974...	goed	153147.0	206931.0	14.12	1.0	2.5	KL	NaN	False	N	XZ	NaN	NaN	NaN	NaN	NaN	NaN
3	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1974...	goed	153147.0	206931.0	14.12	2.5	3.5	FZ	NaN	False	N	GL	False	N	KL	NaN	NaN	NaN
4	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1974...	goed	153147.0	206931.0	14.12	3.5	7.0	FZ	NaN	False	N	GL	NaN	NaN	NaN	NaN	NaN	NaN

The dataframe contains one ‘Gecodeerde lithologie’ interpretation where five layers (‘laag’) were identified. The available data are flattened to represent unique attributes per row of the dataframe.

Using the pkey_interpretatie field one can request the details of this interpretation in a webbrowser:

for pkey_interpretatie in set(df.pkey_interpretatie):
    print(pkey_interpretatie)

https://www.dov.vlaanderen.be/data/interpretatie/2001-186544

Get ‘Gecodeerde lithologie’ with specific properties within a distance from a point

Next to querying interpretations based on their geographic location within a bounding box, we can also search for interpretations matching a specific set of properties. For this we can build a query using a combination of the ‘Gecodeerde lithologie’ fields and operators provided by the WFS protocol.

A list of possible operators can be found below:

[i for i,j in inspect.getmembers(sys.modules['owslib.fes2'], inspect.isclass) if 'Property' in i]

['PropertyIsBetween',
 'PropertyIsEqualTo',
 'PropertyIsGreaterThan',
 'PropertyIsGreaterThanOrEqualTo',
 'PropertyIsLessThan',
 'PropertyIsLessThanOrEqualTo',
 'PropertyIsLike',
 'PropertyIsNotEqualTo',
 'PropertyIsNull',
 'SortProperty']

In this example we build a query using the PropertyIsGreaterThan and PropertyIsEqualTo operators to find all interpretations that are at least 20 m deep, that are deemed appropriate for a range of 1 km from a defined point:

from owslib.fes2 import And, PropertyIsGreaterThan, PropertyIsEqualTo
from pydov.util.location import WithinDistance, Point

query = And([PropertyIsEqualTo(propertyname='Betrouwbaarheid',
                              literal='goed'),
            PropertyIsGreaterThan(propertyname='diepte_tot_m',
                                 literal='20'),
           ])

df = itp.search(query=query,
                location=WithinDistance(Point(153145, 206930, epsg=31370), 1000))

df.head()

[000/001] .
[000/002] ..

	pkey_interpretatie	pkey_boring	betrouwbaarheid_interpretatie	x	y	start_interpretatie_mtaw	diepte_laag_van	diepte_laag_tot	hoofdnaam1_grondsoort	hoofdnaam2_grondsoort	bijmenging1_plaatselijk	bijmenging1_hoeveelheid	bijmenging1_grondsoort	bijmenging2_plaatselijk	bijmenging2_hoeveelheid	bijmenging2_grondsoort	bijmenging3_plaatselijk	bijmenging3_hoeveelheid	bijmenging3_grondsoort
0	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1994...	goed	152566.5	207059.57	18.71	0.0	1.0	FZ	NaN	False	M	KA	NaN	NaN	NaN	NaN	NaN	NaN
1	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1994...	goed	152566.5	207059.57	18.71	1.0	2.0	FZ	NaN	False	M	SF	False	M	ST	NaN	NaN	NaN
2	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1994...	goed	152566.5	207059.57	18.71	2.0	4.0	FZ	NaN	False	M	SF	False	M	ST	NaN	NaN	NaN
3	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1994...	goed	152566.5	207059.57	18.71	4.0	18.0	FZ	NaN	False	M	GL	False	M	SF	False	M	ST
4	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1994...	goed	152566.5	207059.57	18.71	18.0	19.0	FZ	NaN	False	S	LE	NaN	NaN	NaN	NaN	NaN	NaN

Once again we can use the pkey_interpretatie as a permanent link to the information of these interpretations:

for pkey_interpretatie in set(df.pkey_interpretatie):
    print(pkey_interpretatie)

https://www.dov.vlaanderen.be/data/interpretatie/2016-294473
https://www.dov.vlaanderen.be/data/interpretatie/2016-294475

Get ‘Gecodeerde lithologie’ in a bounding box based on specific properties

We can combine a query on attributes with a query on geographic location to get the interpretations within a bounding box that have specific properties.

The following example requests the interpretations of boreholes only, within the given bounding box.

(Note that the datatype of the literal parameter should be a string, regardless of the datatype of this field in the output dataframe.)

from owslib.fes2 import PropertyIsEqualTo

query = PropertyIsEqualTo(
            propertyname='Type_proef',
            literal='Boring')

df = itp.search(
    location=Within(Box(153145, 206930, 154145, 207930, epsg=31370)),
    query=query
    )

df.head()

[000/001] .
[000/003] c..

	pkey_interpretatie	pkey_boring	betrouwbaarheid_interpretatie	x	y	start_interpretatie_mtaw	diepte_laag_van	diepte_laag_tot	hoofdnaam1_grondsoort	hoofdnaam2_grondsoort	bijmenging1_plaatselijk	bijmenging1_hoeveelheid	bijmenging1_grondsoort	bijmenging2_plaatselijk	bijmenging2_hoeveelheid	bijmenging2_grondsoort	bijmenging3_plaatselijk	bijmenging3_hoeveelheid	bijmenging3_grondsoort
0	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1974...	goed	153147.0	206931.0	14.12	0.0	0.5	LE	NaN	False	N	XZ	NaN	NaN	NaN	NaN	NaN	NaN
1	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1974...	goed	153147.0	206931.0	14.12	0.5	1.0	LE	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
2	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1974...	goed	153147.0	206931.0	14.12	1.0	2.5	KL	NaN	False	N	XZ	NaN	NaN	NaN	NaN	NaN	NaN
3	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1974...	goed	153147.0	206931.0	14.12	2.5	3.5	FZ	NaN	False	N	GL	False	N	KL	NaN	NaN	NaN
4	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1974...	goed	153147.0	206931.0	14.12	3.5	7.0	FZ	NaN	False	N	GL	NaN	NaN	NaN	NaN	NaN	NaN

We can look at one of the interpretations in a webbrowser using its pkey_interpretatie:

for pkey_interpretatie in set(df.pkey_interpretatie):
    print(pkey_interpretatie)

https://www.dov.vlaanderen.be/data/interpretatie/2001-186544
https://www.dov.vlaanderen.be/data/interpretatie/2024-379798
https://www.dov.vlaanderen.be/data/interpretatie/2024-379795

Get ‘Gecodeerde lithologie’ based on fields not available in the standard output dataframe

To keep the output dataframe size acceptable, not all available WFS fields are included in the standard output. However, one can use this information to select interpretations as illustrated below.

For example, make a selection of the interpretations in municipality the of Antwerp, before 1/1/1990:

!remark: mind that the municipality attribute is merely an attribute that is defined by the person entering the data. It can be ok, empty, outdated or wrong!

from owslib.fes2 import And, PropertyIsEqualTo, PropertyIsLessThan

query = And([PropertyIsEqualTo(propertyname='gemeente',
                               literal='Antwerpen'),
             PropertyIsLessThan(propertyname='Datum',
                                 literal='2010-01-01')]
            )
df = itp.search(query=query,
                return_fields=('pkey_interpretatie', 'Datum'))
df.head()

[000/001] .

	pkey_interpretatie	Datum
0	https://www.dov.vlaanderen.be/data/interpretat...	2008-03-28
1	https://www.dov.vlaanderen.be/data/interpretat...	2008-03-28
2	https://www.dov.vlaanderen.be/data/interpretat...	2001-07-04
3	https://www.dov.vlaanderen.be/data/interpretat...	2001-10-18
4	https://www.dov.vlaanderen.be/data/interpretat...	2001-07-09

Get ‘Gecodeerde lithologie’ data, returning fields not available in the standard output dataframe

As denoted in the previous example, not all available fields are available in the default output frame to keep its size limited. However, you can request any available field by including it in the return_fields parameter of the search:

query = PropertyIsEqualTo(
            propertyname='gemeente',
            literal='Leuven')

df = itp.search(query=query,
                return_fields=('pkey_interpretatie', 'pkey_boring',
                               'x', 'y', 'start_interpretatie_mtaw', 'gemeente', 'Auteurs', 'Proefnummer'))

df.head()

[000/001] .

	pkey_interpretatie	pkey_boring	x	y	start_interpretatie_mtaw	gemeente	Auteurs	Proefnummer
0	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1932...	173252.0	179257.0	17.0	Leuven	Vanguestaine, M. - onbekend	kb24d74e-B75
1	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1875...	174974.0	174852.0	30.0	Leuven	Vanguestaine, M. - onbekend	kb32d89e-B1
2	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1932...	173252.0	179257.0	17.0	Leuven	Buffel, Ph. - Matthijs-Buffel	kb24d74e-B75
3	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1966...	173909.0	180204.0	11.0	Leuven	Buffel, Ph. - Matthijs-Buffel	kb24d74e-B114
4	https://www.dov.vlaanderen.be/data/interpretat...	https://www.dov.vlaanderen.be/data/boring/1970...	172553.0	180685.0	16.0	Leuven	Buffel, Ph. - Matthijs-Buffel	kb24d74e-B115

Visualize results

Using Folium, we can display the results of our search on a map.

# import the necessary modules (not included in the requirements of pydov!)
import folium
from folium.plugins import MarkerCluster
from pyproj import Transformer

# convert the coordinates to lat/lon for folium
def convert_latlon(x1, y1):
    transformer = Transformer.from_crs("epsg:31370", "epsg:4326", always_xy=True)
    x2,y2 = transformer.transform(x1, y1)
    return x2, y2

df['lon'], df['lat'] = zip(*map(convert_latlon, df['x'], df['y']))
# convert to list
loclist = df[['lat', 'lon']].values.tolist()

# initialize the Folium map on the centre of the selected locations, play with the zoom until ok
fmap = folium.Map(location=[df['lat'].mean(), df['lon'].mean()], zoom_start=12)
marker_cluster = MarkerCluster().add_to(fmap)
for loc in range(0, len(loclist)):
    folium.Marker(loclist[loc], popup=df['Proefnummer'][loc]).add_to(marker_cluster)
fmap

Make this Notebook Trusted to load map: File -> Trust Notebook