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.