DIS3.1 from the Deltares OpenDAP server

DIS (Delfstoffen Informatie Systeem) is a geological voxelmodel created for sand mining within in the 12-mile zone of the Dutch North Sea. See this Deltares webpage for more information on the DIS

ℹ️ This example is work-in-progress. Additional features and examples will follow in future versions of GeoST

In this example we will download a small part of the DIS model using a bounding box in WGS84 UTM 31N coordinates (EPSG:32631) that roughly covers the offshore sand mining areas Q16-6 and Q16-4.

import geost

# Load DIS voxel model from OpenDAP server with bounding box
dis = geost.read_model_from_opendap(
    r"https://opendap.deltares.nl/thredds/dodsC/opendap/rijkswaterstaat/DIS/DIS3.1.nc",
    bbox=(561_000, 5_770_000, 569_000, 5_775_500),
)
dis  # Inspect the loaded DIS voxel model

<xarray.Dataset> Size: 48MB
Dimensions:               (y: 28, x: 40, z: 94)
Coordinates:
  * y                     (y) float64 224B 5.77e+06 5.77e+06 ... 5.776e+06
  * x                     (x) float64 320B 5.611e+05 5.613e+05 ... 5.689e+05
  * z                     (z) float32 376B -55.25 -54.75 -54.25 ... -9.25 -8.75
    spatial_ref           int64 8B 0
Data variables: (12/114)
    lithoklasse           (y, x, z) float32 421kB ...
    slibklasse            (y, x, z) float32 421kB ...
    schelpenklasse        (y, x, z) float32 421kB ...
    stratunit             (y, x, z) float32 421kB ...
    kans_lithoklasse_1    (y, x, z) float32 421kB ...
    kans_lithoklasse_2    (y, x, z) float32 421kB ...
    ...                    ...
    litho_realisatie_096  (y, x, z) float32 421kB ...
    litho_realisatie_097  (y, x, z) float32 421kB ...
    litho_realisatie_098  (y, x, z) float32 421kB ...
    litho_realisatie_099  (y, x, z) float32 421kB ...
    litho_realisatie_100  (y, x, z) float32 421kB ...
    bathy                 (y, x) float32 4kB ...
Attributes:
    epsg:                    32631
    Units:                   Meter
    Resolution:              200 x 200 x 0.5 m (y, x, z)
    Coordinates:             Cellcenters
    Model:                   DIS 3.1
    Source_software:         Isatis
    Organisation:            Rijkswaterstaat, Deltares, TNO
    Conventions:             CF-1.0
    Creator_NetCDF_dataset:  Deltares
    E-mail:                  tommer.vermaas@deltares.nl

xarray.Dataset

• Dimensions:
  • y: 28
  • x: 40
  • z: 94
• Coordinates: (4)
  • y
    (y)
    float64
    5.77e+06 5.77e+06 ... 5.776e+06

    axis :

    Y

    long_name :

    y coordinate of projection

    standard_name :

    projection_y_coordinate

    units :

    metre

    array([5770100., 5770300., 5770500., 5770700., 5770900., 5771100., 5771300.,
           5771500., 5771700., 5771900., 5772100., 5772300., 5772500., 5772700.,
           5772900., 5773100., 5773300., 5773500., 5773700., 5773900., 5774100.,
           5774300., 5774500., 5774700., 5774900., 5775100., 5775300., 5775500.])

  • x
    (x)
    float64
    5.611e+05 5.613e+05 ... 5.689e+05

    axis :

    X

    long_name :

    x coordinate of projection

    standard_name :

    projection_x_coordinate

    units :

    metre

    array([561100., 561300., 561500., 561700., 561900., 562100., 562300., 562500.,
           562700., 562900., 563100., 563300., 563500., 563700., 563900., 564100.,
           564300., 564500., 564700., 564900., 565100., 565300., 565500., 565700.,
           565900., 566100., 566300., 566500., 566700., 566900., 567100., 567300.,
           567500., 567700., 567900., 568100., 568300., 568500., 568700., 568900.])

  • z
    (z)
    float32
    -55.25 -54.75 ... -9.25 -8.75

    reference :

    MSL

    units :

    m

    array([-55.25, -54.75, -54.25, -53.75, -53.25, -52.75, -52.25, -51.75, -51.25,
           -50.75, -50.25, -49.75, -49.25, -48.75, -48.25, -47.75, -47.25, -46.75,
           -46.25, -45.75, -45.25, -44.75, -44.25, -43.75, -43.25, -42.75, -42.25,
           -41.75, -41.25, -40.75, -40.25, -39.75, -39.25, -38.75, -38.25, -37.75,
           -37.25, -36.75, -36.25, -35.75, -35.25, -34.75, -34.25, -33.75, -33.25,
           -32.75, -32.25, -31.75, -31.25, -30.75, -30.25, -29.75, -29.25, -28.75,
           -28.25, -27.75, -27.25, -26.75, -26.25, -25.75, -25.25, -24.75, -24.25,
           -23.75, -23.25, -22.75, -22.25, -21.75, -21.25, -20.75, -20.25, -19.75,
           -19.25, -18.75, -18.25, -17.75, -17.25, -16.75, -16.25, -15.75, -15.25,
           -14.75, -14.25, -13.75, -13.25, -12.75, -12.25, -11.75, -11.25, -10.75,
           -10.25,  -9.75,  -9.25,  -8.75], dtype=float32)

  • spatial_ref
    ()
    int64
    0

    crs_wkt :

    PROJCS["unnamed",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",3],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["METERS",1],AXIS["Easting",EAST],AXIS["Northing",NORTH]]

    semi_major_axis :

    6378137.0

    semi_minor_axis :

    6356752.314245179

    inverse_flattening :

    298.257223563

    reference_ellipsoid_name :

    WGS 84

    longitude_of_prime_meridian :

    0.0

    prime_meridian_name :

    Greenwich

    geographic_crs_name :

    WGS 84

    horizontal_datum_name :

    World Geodetic System 1984

    projected_crs_name :

    unnamed

    grid_mapping_name :

    transverse_mercator

    latitude_of_projection_origin :

    0.0

    longitude_of_central_meridian :

    3.0

    false_easting :

    500000.0

    false_northing :

    0.0

    scale_factor_at_central_meridian :

    0.9996

    spatial_ref :

    PROJCS["unnamed",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",3],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["METERS",1],AXIS["Easting",EAST],AXIS["Northing",NORTH]]

    GeoTransform :

    561000.0 200.0 0.0 5770000.0 0.0 200.0

    array(0)

• Data variables: (114)
  • lithoklasse
    (y, x, z)
    float32
    ...

    1 :

    Zand uiterst fijn (63-105)

    2 :

    Zand matig fijn-zeer fijn (105-210)

    3 :

    Zand matig grof (210-300)

    4 :

    Zand zeer grof (300-420)

    5 :

    Zand uiterst grof (420-2000)

    6 :

    Zand (ongedifferentieerd)

    7 :

    Klei/Leem

    8 :

    Veen

    9 :

    Schelpen

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • slibklasse
    (y, x, z)
    float32
    ...

    1 :

    Zwak slibhoudend (0 - 2%)

    2 :

    Matig zwak slibhoudend (2 - 4%)

    3 :

    Matig sterk slibhoudend (4 - 10%)

    4 :

    Sterk slibhoudend (10 - 30%)

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • schelpenklasse
    (y, x, z)
    float32
    ...

    1 :

    Geen schelpen (0%)

    2 :

    Spoor schelpen (0 - 1%)

    3 :

    Weinig schelpen (1 - 10%)

    4 :

    Veel schelpen (10 - 30%)

    5 :

    Schelpen % onbekend

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • stratunit
    (y, x, z)
    float32
    ...

    0 :

    Southern Bight

    1 :

    Naaldwijk 1

    2 :

    Naaldwijk 2

    3 :

    Naaldwijk 3

    4 :

    Vroeg Holoceen

    5 :

    Pleistoceen

    6 :

    Paleogeen

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • kans_lithoklasse_1
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • kans_lithoklasse_2
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • kans_lithoklasse_3
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • kans_lithoklasse_4
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • kans_lithoklasse_5
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • kans_lithoklasse_6
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • kans_lithoklasse_7
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • kans_lithoklasse_8
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • kans_lithoklasse_9
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_001
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_002
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_003
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_004
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_005
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_006
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_007
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_008
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_009
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_010
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_011
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_012
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_013
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_014
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_015
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_016
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_017
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_018
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_019
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_020
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_021
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_022
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_023
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_024
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_025
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_026
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_027
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_028
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_029
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_030
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_031
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_032
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_033
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_034
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_035
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_036
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_037
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_038
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_039
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_040
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_041
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_042
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_043
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_044
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_045
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_046
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_047
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_048
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_049
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_050
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_051
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_052
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_053
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_054
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_055
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_056
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_057
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_058
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_059
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_060
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_061
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_062
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_063
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_064
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_065
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_066
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_067
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_068
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_069
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_070
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_071
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_072
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_073
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_074
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_075
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_076
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_077
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_078
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_079
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_080
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_081
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_082
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_083
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_084
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_085
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_086
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_087
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_088
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_089
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_090
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_091
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_092
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_093
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_094
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_095
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_096
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_097
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_098
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_099
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • litho_realisatie_100
    (y, x, z)
    float32
    ...

    _ChunkSizes :

    [30 60 12]

    [105280 values with dtype=float32]

  • bathy
    (y, x)
    float32
    ...

    reference :

    MSL

    units :

    m

    _ChunkSizes :

    [59 60]

    [1120 values with dtype=float32]

• Attributes: (10)

  epsg :

  32631

  Units :

  Meter

  Resolution :

  200 x 200 x 0.5 m (y, x, z)

  Coordinates :

  Cellcenters

  Model :

  DIS 3.1

  Source_software :

  Isatis

  Organisation :

  Rijkswaterstaat, Deltares, TNO

  Conventions :

  CF-1.0

  Creator_NetCDF_dataset :

  Deltares

  E-mail :

  tommer.vermaas@deltares.nl

We can apply any method available from GeoST’s extension of Xarray to further select and analyse the downloaded DIS model.

Let’s e.g. find the total thickness of fine sand and very fine sand between the -20 and -25 m MSL depth interval. For this we will use the ‘lithoklasse’ data variable which classifies lithology. First, let’s do a quick inspection of the “lithoklasse” variable.

print(dis["lithoklasse"].attrs)

{'1': 'Zand uiterst fijn (63-105)', '2': 'Zand matig fijn-zeer fijn (105-210)', '3': 'Zand matig grof (210-300)', '4': 'Zand zeer grof (300-420)', '5': 'Zand uiterst grof (420-2000)', '6': 'Zand (ongedifferentieerd)', '7': 'Klei/Leem', '8': 'Veen', '9': 'Schelpen', '_ChunkSizes': array([30, 60, 12], dtype=int32)}

We can see that lithoclasses 1 and 2 represent the required sand fractions so we use these to determine the total thickness of (very) fine sand within our depth interval.

dis_selected = dis.gst.slice_depth_interval(upper=-20, lower=-25)
dis_sand = dis_selected.gst.get_thickness(dis_selected["lithoklasse"] <= 2)

dis_sand.plot.imshow(
    figsize=(8, 6),
    cmap="YlOrBr",
    vmin=0,
    vmax=4,
    cbar_kwargs={"label": "Total thickness of fine sand (m)"},
)

<matplotlib.image.AxesImage at 0x7b55387f5550>