tomvansteijn / xsboringen

plot borehole data in cross-sections
MIT License
20 stars 10 forks source link

Small fixes to raster surfaces and CPT #6

Closed erikgeo closed 3 years ago

erikgeo commented 3 years ago

Three fixes, one being more of a workaround:

  1. It is now possible to have a profile line outside of the raster's spatial extent without getting an IndexError
  2. Fixed check for NaN values in rasters. This went wrong since np.nan == np.nan returns False.
  3. Ignore friction ratio if not available. Some CPT's don't have the luxury of a friction ratio and it would stop the program. If this is the case and translate_cpts = True in the config, a plot will look like this now:

image

Geweldige package overigens. Stond op het punt zelf zoiets te ontwikkelen tot een collega me hier op wees. Zie maar wat je met de fixes doet, de workaround past misschien niet heel mooi in het geheel maar werkt en had ik nodig. Bedankt!

Groeten, Erik van Onselen (Geoloog bij Deltares)

tomvansteijn commented 3 years ago

Hi Erik,

Thanks for your work and appreciation. Merged!

This package has not been updated in years. I do have a local branch with a fix for the correct representation of the blind part of a CPT. And a better sampling method for rasters (using rasterio 's rasterize). I will try to merge soon. The visualisation of well screens is still on my personal wishlist, as well as a better translation method for CPT to lithology, like Robertson.

Feel free to make changes and suggestions for improvement.