geopandas 1.2.0 | pyproj 3.8.0 | python-ags4 1.2.0

=== Step 1: LOCA to points in British National Grid ===
112 holes | CRS OSGB36 / British National Grid (EPSG:27700)
extent E 392,081 to 394,938, N 213,164 to 216,419
centre in WGS84: 51.83674 N, 2.09630 W (Crickley Hill and Birdlip)

=== Step 2: what goes wrong with the wrong CRS ===
CP102 with easting and northing swapped lands at 53.393 N, -4.773 W: in the sea off north Wales, 249 km away
tagged as EPSG:4326 instead, x runs to 394,938 'degrees'; valid longitudes stop at 180, so nothing can be drawn

=== Step 3: summarise GEOL for each hole ===
type_name
Rotary Cored                          71
Trial Pit/trench                      30
Cone Penetration Test                  6
Rotary open hole                       3
Cable Percussion (shell and auger)     2

90 holes reach bedrock; depth to rock 0.10 to 20.70 m (median 1.20 m)
unit at rockhead: {'FE': 20, 'BLPL': 17, 'WHM': 15, 'SALS': 9, 'ASLS': 5, 'GOG': 5, 'DYS': 5, 'HMB': 5, 'WHL': 4, 'CHAM': 3, 'MRB': 1, 'IOG': 1}
LOCA_ID    type_name  LOCA_GL  LOCA_FDEP  rockhead_depth  rockhead_mOD             rockhead_name first_water_strike  n_samples
  CP217 Rotary Cored   194.95       33.5            9.50        185.45 Whitby Mudstone Formation               2.25         61
DSRC415 Rotary Cored   287.20       51.0            0.75        286.45  Fuller's Earth Formation                NaN         53
  RC516 Rotary Cored   277.35       22.2            1.20        276.15  Fuller's Earth Formation                NaN         35

=== Step 4: check LOCA_GL against Environment Agency LIDAR ===
DTM 2000 x 2000 cells of 2 m | 81.7 to 302.9 m OD | no-data cells: 0
LOCA_GL minus LIDAR: median -0.13 m | middle half -0.23 to -0.04 m | within 0.5 m: 103 of 112 | within 1 m: 109
more than 1 m apart:
LOCA_ID        type_name  LOCA_CLST  LOCA_GL   lidar_gl  gl_minus_lidar
  TP602 Trial Pit/trench PHASE 2A.1   274.40 275.630005           -1.23
  CP106     Rotary Cored PHASE 2A.1   186.25 187.350006           -1.10
  CP204     Rotary Cored PHASE 2A.1   178.90 177.770004            1.13

=== Step 5: maps ===

=== Step 6: write a GeoPackage for QGIS ===
wrote output/a417_boreholes.gpkg: layer 'boreholes', 112 points, 19 fields, CRS 27700
rasters for the QGIS basemap: a417_dtm_2m.tif, a417_hillshade.tif
figures written: 01-ground-level-check.png, 02-borehole-map.png, 03-rockhead-unit-map.png, 04-qgis-borehole-map.png