Retrieving and analysing BHR-GT grainsize samples#

In this example, we will download BHR-GT geotechnical boreholes (BRO BHR-GT objects) and associated grainsize sample data.

ℹ️ This example is work-in-progress. Additional features developed for working with grainsize data will follow in future versions of GeoST

You must download the BHR-GT borehole data (i.e. layer descriptions) and grainsize sample data seperately. In this example we will use a bounding box in WGS84 UTM 31N coordinates (EPSG:32631) that roughly covers the offshore sand mining area L12-2 in the Dutch North Sea

import geost

# BHR-GT cores in sand mining area L12-2, using the object type 'BHR-GT'
bhrgt_cores = geost.bro_api_read(
    "BHR-GT", bbox=(619_000, 592_400_0, 625_000, 593_000_0), epsg=32631
)

# Corresponding samples, using the object type 'BHR-GT-samples'
bhrgt_samples = geost.bro_api_read(
    "BHR-GT-samples", bbox=(619_000, 592_400_0, 625_000, 593_000_0), epsg=32631
)

# Show results of downloaded samples
print("\nThe data tables look like this:")
print(bhrgt_samples.data.head())
The data tables look like this:
                nr              x              y  surface    end   top  \
0  BHR000000451643  119007.510234  610338.277226   -27.13 -32.08  0.20   
1  BHR000000451643  119007.510234  610338.277226   -27.13 -32.08  0.65   
2  BHR000000451643  119007.510234  610338.277226   -27.13 -32.08  1.60   
3  BHR000000451643  119007.510234  610338.277226   -27.13 -32.08  2.20   
4  BHR000000451643  119007.510234  610338.277226   -27.13 -32.08  3.20   

   bottom fractionSmaller63um fractionLarger63um fraction63to90um  ...  \
0    0.40                 0.0              100.0              0.2  ...   
1    0.85                 2.2               97.8              2.6  ...   
2    1.80                 1.7               98.3              1.8  ...   
3    2.40                 1.9               98.1              2.1  ...   
4    3.40                 1.0               99.0              0.4  ...   

  fraction355to500um fraction500to710um fraction710to1000um  \
0                6.6                0.8                 0.3   
1                0.5                0.1                 0.1   
2                0.8                0.3                 0.1   
3                0.9                0.4                 0.1   
4               13.6               49.4                  17   

  fraction1000to1400um fraction1400umto2mm fraction2to4mm fraction4to8mm  \
0                  0.2                 0.2            1.2              0   
1                    0                 0.1            0.2              0   
2                    0                   0              0              0   
3                    0                   0              0              0   
4                    3                 0.7            2.7              0   

  fraction16to31_5mm fraction31_5to63mm fractionLarger63mm  
0                  0                  0                  0  
1                  0                  0                  0  
2                  0                  0                  0  
3                  0                  0                  0  
4                  0                  0                  0  

[5 rows x 24 columns]

For more info, see the API reference of the bro_api_read function.