geost 0.2.0__tar.gz

geost-0.2.0/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2022 Deltares / TGG projects / Subsurface Toolbox
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

geost-0.2.0/PKG-INFO

@@ -0,0 +1,142 @@
+Metadata-Version: 2.1
+Name: geost
+Version: 0.2.0
+Summary: The Python interface for geological point data in The Netherlands
+Author-email: Erik van Onselen <erik.vanonselen@deltares.nl>, Bas Knaake <bas.knaake@deltares.nl>, Tommer Vermaas <tommer.vermaas@deltares.nl>
+License: MIT License
+        
+        Copyright (c) 2022 Deltares / TGG projects / Subsurface Toolbox
+        
+        Permission is hereby granted, free of charge, to any person obtaining a copy
+        of this software and associated documentation files (the "Software"), to deal
+        in the Software without restriction, including without limitation the rights
+        to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+        copies of the Software, and to permit persons to whom the Software is
+        furnished to do so, subject to the following conditions:
+        
+        The above copyright notice and this permission notice shall be included in all
+        copies or substantial portions of the Software.
+        
+        THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+        IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+        FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+        AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+        LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+        OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+        SOFTWARE.
+        
+Project-URL: Homepage, https://github.com/Deltares-research/geost
+Project-URL: Repository, https://github.com/Deltares-research/geost
+Project-URL: Issues, https://github.com/Deltares-research/geost/issues
+Keywords: geost,bro,borehole,well log,cpt,geology,subsurface toolbox,deltares
+Classifier: Development Status :: 2 - Pre-Alpha
+Classifier: Intended Audience :: Science/Research
+Classifier: Topic :: Scientific/Engineering
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Operating System :: OS Independent
+Classifier: Programming Language :: Python
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Topic :: Scientific/Engineering
+Requires-Python: >=3.12
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: numpy
+Requires-Dist: pandas>=2.0.1
+Requires-Dist: pyarrow
+Requires-Dist: lxml
+Requires-Dist: pyproj
+Requires-Dist: geopandas
+Requires-Dist: xlrd
+Requires-Dist: openpyxl
+Requires-Dist: rioxarray>=0.15.3
+Requires-Dist: shapely>=2.0.1
+Requires-Dist: pyvista>=0.36.1
+Provides-Extra: dev
+Requires-Dist: pyvista; extra == "dev"
+Requires-Dist: black; extra == "dev"
+Requires-Dist: pytest; extra == "dev"
+Requires-Dist: pytest-cov; extra == "dev"
+Requires-Dist: ruff; extra == "dev"
+Requires-Dist: sphinx; extra == "dev"
+Requires-Dist: sphinx-rtd-theme; extra == "dev"
+
+# GeoST - Geological Subsurface Toolbox
+[![License: MIT](https://img.shields.io/pypi/l/imod)](https://choosealicense.com/licenses/mit)
+[![Lifecycle: experimental](https://lifecycle.r-lib.org/articles/figures/lifecycle-experimental.svg)](https://lifecycle.r-lib.org/articles/stages.html)
+[![Build: status](https://img.shields.io/github/actions/workflow/status/deltares-research/geost/ci.yml)](https://github.com/Deltares-research/geost/actions)
+[![codecov](https://codecov.io/gh/Deltares-research/geost/graph/badge.svg?token=HCNGLWTQ2H)](https://codecov.io/gh/Deltares-research/geost)
+[![Formatting: ruff](https://camo.githubusercontent.com/18c26428c337f9d641fa09b629a3a03b514e8ac84b57974a0ed7d1b38e14e060/68747470733a2f2f696d672e736869656c64732e696f2f656e64706f696e743f75726c3d68747470733a2f2f7261772e67697468756275736572636f6e74656e742e636f6d2f61737472616c2d73682f727566662f6d61696e2f6173736574732f62616467652f76322e6a736f6e)](https://github.com/charliermarsh/ruff)
+
+The Geological Subsurface Toolbox (geost) package is designed to be an easy-to-use Python interface for working with subsurface point data in The Netherlands (boreholes, well logs and CPT's). It provides selection, analysis and export methods that can be applied generically to the loaded data. It is designed to connect with other Deltares developments such as [iMod](https://gitlab.com/deltares/imod) and [DataFusionTools](https://bitbucket.org/DeltaresGEO/datafusiontools/src/master/).
+
+The internal BoreholeCollection, LogCollection and CptCollection classes use [Pandas](https://pandas.pydata.org/) for storing data and header information. It utilizes a custom, lightweight validation module inspired by the [Pandera](https://pandera.readthedocs.io/en/stable/) API. For spatial functions [Geopandas](https://geopandas.org/en/stable/) is used. The package also supports reading/writing parquet and geoparquet files through Pandas and Geopandas respectively. 
+
+GeoST is a work-in-progress and currently support a limited number of data sources.
+
+## Installation (user)
+In a Python > 3.12 environment, install the latest stable release using pip:
+
+    pip install geost
+
+Or the latest (experimental) version of the main branch directly from GitHub using:
+
+    pip install git+https://github.com/Deltares-research/geost.git
+
+## Installation (developer)
+We use [Pixi](https://github.com/prefix-dev/pixi) for package management and workflows.
+
+With pixi installed, navigate to the folder of the cloned repository and run:
+
+    pixi install
+
+To install all GeoST dependencies. Next install GeoST in editable mode by running the 
+pixi task 'install':
+
+    pixi run install
+
+See the [Pixi documentation](https://pixi.sh/latest/) for more information.
+
+## Supported borehole and CPT formats
+- From local files
+    - geost .parquet file (complete)
+    - Dino csv geological boreholes (complete)
+    - Dino XML geological boreholes (planned)
+    - BRO XML geotechnical boreholes (planned)
+    - BRO XML soil boreholes (planned)
+    - GEF boreholes (planned)
+    - GEF CPT's (complete)
+    - BRO XML CPT's (planned)
+    - BRO geopackage CPT's (planned)
+    - Well log LAS files (planned)
+    - Well log ASCII files (planned)
+- Directly from the BRO (REST API) (all planned)
+    - CPT
+    - BHR-P
+    - BHR-GT
+    - BHR-G
+
+## Features
+After loading data from one of the supported formats it will automatically be validated. If the validation is succesful, a Collection object will be returned depending on your input data type (mixed CPT/well log/borehole collections are not allowed). A collection object consists of two main attributes: the **header** and **data**. The header contains a table with one entry per object and provides information about the name, location, surface level, and borehole/log/cpt start and end depths. The data attribute is a table that includes the data for every described layer (boreholes) or measurement (well logs, cpt's).
+
+The collection object comes with a comprehensive set of methods that can be applied generically while ensuring that the header and data remain synchronized:
+
+- Selection/slicing methods (e.g., objects within bounding box, within or close to geometries, based on depth and other conditions)   
+- Export methods (e.g. to csv, parquet, geopackage, VTK, DataFusionTools, Kingdom* , etc)
+- Datafusion methods (e.g. combining collections*, combining with data from maps, conversion of description protocols*)
+- Miscellaneous methods (e.g. changing vertical/horizontal position reference system)
+
+For a better overview of basic functionality, see the [Basics Tutorial](https://gitlab.com/deltares/tgg-projects/subsurface-toolbox/pysst/-/blob/main/tutorials/basics.ipynb?ref_type=heads).
+
+## Contributing
+
+You can contribute by testing, raising issues and making pull requests. Some general guidelines:
+
+- Use new branches for developing new features or bugfixes. Use prefixes such as feature/ bugfix/ experimental/ to indicate the type of branch
+- Add unit tests (and test data) for new methods and functions. We use pytest.
+- Add Numpy-style docstrings
+- Use Black formatting with default line lenght (88 characters)
+- Update requirement.txt en environment.yml files if required
+
+## License
+MIT license (Note: may change to a copyleft license in the future, depending on Deltares management decisions)
+

geost-0.2.0/README.md

@@ -0,0 +1,80 @@
+# GeoST - Geological Subsurface Toolbox
+[![License: MIT](https://img.shields.io/pypi/l/imod)](https://choosealicense.com/licenses/mit)
+[![Lifecycle: experimental](https://lifecycle.r-lib.org/articles/figures/lifecycle-experimental.svg)](https://lifecycle.r-lib.org/articles/stages.html)
+[![Build: status](https://img.shields.io/github/actions/workflow/status/deltares-research/geost/ci.yml)](https://github.com/Deltares-research/geost/actions)
+[![codecov](https://codecov.io/gh/Deltares-research/geost/graph/badge.svg?token=HCNGLWTQ2H)](https://codecov.io/gh/Deltares-research/geost)
+[![Formatting: ruff](https://camo.githubusercontent.com/18c26428c337f9d641fa09b629a3a03b514e8ac84b57974a0ed7d1b38e14e060/68747470733a2f2f696d672e736869656c64732e696f2f656e64706f696e743f75726c3d68747470733a2f2f7261772e67697468756275736572636f6e74656e742e636f6d2f61737472616c2d73682f727566662f6d61696e2f6173736574732f62616467652f76322e6a736f6e)](https://github.com/charliermarsh/ruff)
+
+The Geological Subsurface Toolbox (geost) package is designed to be an easy-to-use Python interface for working with subsurface point data in The Netherlands (boreholes, well logs and CPT's). It provides selection, analysis and export methods that can be applied generically to the loaded data. It is designed to connect with other Deltares developments such as [iMod](https://gitlab.com/deltares/imod) and [DataFusionTools](https://bitbucket.org/DeltaresGEO/datafusiontools/src/master/).
+
+The internal BoreholeCollection, LogCollection and CptCollection classes use [Pandas](https://pandas.pydata.org/) for storing data and header information. It utilizes a custom, lightweight validation module inspired by the [Pandera](https://pandera.readthedocs.io/en/stable/) API. For spatial functions [Geopandas](https://geopandas.org/en/stable/) is used. The package also supports reading/writing parquet and geoparquet files through Pandas and Geopandas respectively. 
+
+GeoST is a work-in-progress and currently support a limited number of data sources.
+
+## Installation (user)
+In a Python > 3.12 environment, install the latest stable release using pip:
+
+    pip install geost
+
+Or the latest (experimental) version of the main branch directly from GitHub using:
+
+    pip install git+https://github.com/Deltares-research/geost.git
+
+## Installation (developer)
+We use [Pixi](https://github.com/prefix-dev/pixi) for package management and workflows.
+
+With pixi installed, navigate to the folder of the cloned repository and run:
+
+    pixi install
+
+To install all GeoST dependencies. Next install GeoST in editable mode by running the 
+pixi task 'install':
+
+    pixi run install
+
+See the [Pixi documentation](https://pixi.sh/latest/) for more information.
+
+## Supported borehole and CPT formats
+- From local files
+    - geost .parquet file (complete)
+    - Dino csv geological boreholes (complete)
+    - Dino XML geological boreholes (planned)
+    - BRO XML geotechnical boreholes (planned)
+    - BRO XML soil boreholes (planned)
+    - GEF boreholes (planned)
+    - GEF CPT's (complete)
+    - BRO XML CPT's (planned)
+    - BRO geopackage CPT's (planned)
+    - Well log LAS files (planned)
+    - Well log ASCII files (planned)
+- Directly from the BRO (REST API) (all planned)
+    - CPT
+    - BHR-P
+    - BHR-GT
+    - BHR-G
+
+## Features
+After loading data from one of the supported formats it will automatically be validated. If the validation is succesful, a Collection object will be returned depending on your input data type (mixed CPT/well log/borehole collections are not allowed). A collection object consists of two main attributes: the **header** and **data**. The header contains a table with one entry per object and provides information about the name, location, surface level, and borehole/log/cpt start and end depths. The data attribute is a table that includes the data for every described layer (boreholes) or measurement (well logs, cpt's).
+
+The collection object comes with a comprehensive set of methods that can be applied generically while ensuring that the header and data remain synchronized:
+
+- Selection/slicing methods (e.g., objects within bounding box, within or close to geometries, based on depth and other conditions)   
+- Export methods (e.g. to csv, parquet, geopackage, VTK, DataFusionTools, Kingdom* , etc)
+- Datafusion methods (e.g. combining collections*, combining with data from maps, conversion of description protocols*)
+- Miscellaneous methods (e.g. changing vertical/horizontal position reference system)
+
+For a better overview of basic functionality, see the [Basics Tutorial](https://gitlab.com/deltares/tgg-projects/subsurface-toolbox/pysst/-/blob/main/tutorials/basics.ipynb?ref_type=heads).
+
+## Contributing
+
+You can contribute by testing, raising issues and making pull requests. Some general guidelines:
+
+- Use new branches for developing new features or bugfixes. Use prefixes such as feature/ bugfix/ experimental/ to indicate the type of branch
+- Add unit tests (and test data) for new methods and functions. We use pytest.
+- Add Numpy-style docstrings
+- Use Black formatting with default line lenght (88 characters)
+- Update requirement.txt en environment.yml files if required
+
+## License
+MIT license (Note: may change to a copyleft license in the future, depending on Deltares management decisions)
+

geost-0.2.0/geost/__init__.py

@@ -0,0 +1,14 @@
+from geost.read import (
+    get_bro_objects_from_bbox,
+    get_bro_objects_from_geometry,
+    read_borehole_table,
+    read_gef_cores,
+    read_gef_cpts,
+    read_nlog_cores,
+    read_sst_cpts,
+    read_xml_cpts,
+    read_xml_geological_cores,
+    read_xml_geotechnical_cores,
+    read_xml_soil_cores,
+)
+from geost.utils import csv_to_parquet, excel_to_parquet

geost-0.2.0/geost/_borehole.py

@@ -0,0 +1,86 @@
+# class CptCollection(PointDataCollection):
+#     """
+#     Class for collections of CPT data.
+
+#     Users must use the reader functions in
+#     :py:mod:`~geost.read` to create collections. The following readers generate CPT
+#     objects:
+
+#     :func:`~geost.read.read_sst_cpts`, :func:`~geost.read.read_gef_cpts`
+
+#     Args:
+#         data (pd.DataFrame): Dataframe containing borehole/CPT data.
+
+#         vertical_reference (str): Vertical reference, see
+#          :py:attr:`~geost.base.PointDataCollection.vertical_reference`
+
+#         horizontal_reference (int): Horizontal reference, see
+#          :py:attr:`~geost.base.PointDataCollection.horizontal_reference`
+
+#         header (pd.DataFrame): Header used for construction. see
+#          :py:attr:`~geost.base.PointDataCollection.header`
+#     """
+
+#     def __init__(
+#         self,
+#         data: pd.DataFrame,
+#         vertical_reference: str = "NAP",
+#         horizontal_reference: int = 28992,
+#         header: Optional[pd.DataFrame] = None,
+#         is_inclined: bool = False,
+#     ):
+#         super().__init__(
+#             data,
+#             vertical_reference,
+#             horizontal_reference,
+#             header=header,
+#             is_inclined=is_inclined,
+#         )
+
+#     def add_ic(
+#         self,
+#     ):  # Move to cpt analysis functions, use something like 'apply' function in classes
+#         """
+#         Calculate soil behaviour type index (Ic) for all CPT's in the collection.
+
+#         The data is added to :py:attr:`~geost.base.PointDataCollection.header`.
+#         """
+#         self.data["ic"] = calc_ic(self.data["qc"], self.data["friction_number"])
+
+#     def add_lithology(
+#         self,
+#     ):  # Move to cpt analysis functions, use something like 'apply' function in classes
+#         """
+#         Interpret lithoclass for all CPT's in the collection.
+
+#         The data is added to :py:attr:`~geost.base.PointDataCollection.header`.
+#         """
+#         if "ic" not in self.data.columns:
+#             self.add_ic()
+#         self.data["lith"] = calc_lithology(
+#             self.data["ic"], self.data["qc"], self.data["friction_number"]
+#         )
+
+#     def as_boreholecollection(self):  # No change
+#         """
+#         Export CptCollection to BoreholeCollection. Requires the "lith" column to be
+#         present. Use the method :py:meth:`~geost.borehole.CptCollection.add_lithology`
+
+#         Returns
+#         -------
+#         Instance of :class:`~geost.borehole.BoreholeCollection`
+#         """
+#         if "lith" not in self.data.columns:
+#             raise IndexError(
+#                 r"The column \"lith\" is required to convert to BoreholeCollection"
+#             )
+
+#         borehole_converted_dataframe = self.data[
+#             ["nr", "x", "y", "surface", "end", "top", "bottom", "lith"]
+#         ]
+#         cptcollection_as_bhcollection = BoreholeCollection(
+#             borehole_converted_dataframe,
+#             vertical_reference=self.vertical_reference,
+#             header=self.header,
+#         )
+#         return cptcollection_as_bhcollection

geost-0.2.0/geost/abstract_classes.py

@@ -0,0 +1,251 @@
+from abc import ABC, abstractmethod
+
+
+class AbstractHeader(ABC):
+    @property
+    @abstractmethod
+    def gdf(self):
+        pass
+
+    @property
+    @abstractmethod
+    def horizontal_reference(self):
+        pass
+
+    @property
+    @abstractmethod
+    def vertical_reference(self):
+        pass
+
+    @gdf.setter
+    @abstractmethod
+    def gdf(self, gdf):
+        pass
+
+    @abstractmethod
+    def change_horizontal_reference(self):
+        pass
+
+    @abstractmethod
+    def change_vertical_reference(self):
+        pass
+
+    @abstractmethod
+    def get(self):
+        pass
+
+    @abstractmethod
+    def select_within_bbox(self):
+        pass
+
+    @abstractmethod
+    def select_with_points(self):
+        pass
+
+    @abstractmethod
+    def select_with_lines(self):
+        pass
+
+    @abstractmethod
+    def select_within_polygons(self):
+        pass
+
+    @abstractmethod
+    def select_by_depth(self):
+        pass
+
+    @abstractmethod
+    def select_by_length(self):
+        pass
+
+    @abstractmethod
+    def get_area_labels(self):
+        pass
+
+
+class AbstractData(ABC):
+    @property
+    @abstractmethod
+    def df(self):
+        pass
+
+    @property
+    @abstractmethod
+    def datatype(self):
+        pass
+
+    @df.setter
+    @abstractmethod
+    def df(self, df):
+        pass
+
+    @datatype.setter
+    @abstractmethod
+    def datatype(self):
+        pass
+
+    @abstractmethod
+    def to_header(self):
+        pass
+
+    @abstractmethod
+    def to_collection(self):
+        pass
+
+    @abstractmethod
+    def select_by_values(self):
+        pass
+
+    @abstractmethod
+    def slice_depth_interval(self):
+        pass
+
+    @abstractmethod
+    def slice_by_values(self):
+        pass
+
+    @abstractmethod
+    def get_cumulative_layer_thickness(self):
+        # Not sure if this should be here, potentially unsuitable with DiscreteData
+        pass
+
+    @abstractmethod
+    def get_layer_top(self):
+        pass
+
+    @abstractmethod
+    def to_multiblock(self):
+        pass
+
+    @abstractmethod
+    def to_vtm(self):
+        pass
+
+    @abstractmethod
+    def to_datafusiontools(self):
+        # supporting this is low priority, perhaps even deprecate
+        pass
+
+
+class AbstractCollection(ABC):
+    @property
+    @abstractmethod
+    def header(self):
+        pass
+
+    @property
+    @abstractmethod
+    def data(self):
+        pass
+
+    @property
+    def n_points(self):
+        pass
+
+    @property
+    @abstractmethod
+    def horizontal_reference(self):  # Move to header class in future refactor
+        pass
+
+    @property
+    @abstractmethod
+    def vertical_reference(self):  # move to data class in future refactor
+        pass
+
+    @header.setter
+    @abstractmethod
+    def header(self, header):
+        pass
+
+    @data.setter
+    @abstractmethod
+    def data(self, data):
+        pass
+
+    @horizontal_reference.setter
+    @abstractmethod
+    def horizontal_reference(self, to_epsg: int):
+        pass
+
+    @vertical_reference.setter
+    @abstractmethod
+    def vertical_reference(self, to_epsg: str):  # will use epsg after refactor
+        pass
+
+    @abstractmethod
+    def get(self):
+        pass
+
+    @abstractmethod
+    def reset_header(self):
+        pass
+
+    @abstractmethod
+    def check_header_to_data_alignment(self):
+        pass
+
+    @abstractmethod
+    def select_within_bbox(self):
+        pass
+
+    @abstractmethod
+    def select_with_points(self):
+        pass
+
+    @abstractmethod
+    def select_with_lines(self):
+        pass
+
+    @abstractmethod
+    def select_within_polygons(self):
+        pass
+
+    @abstractmethod
+    def select_by_depth(self):
+        pass
+
+    @abstractmethod
+    def select_by_length(self):
+        pass
+
+    @abstractmethod
+    def get_area_labels(self):
+        pass
+
+    @abstractmethod
+    def select_by_values(self):
+        pass
+
+    @abstractmethod
+    def slice_depth_interval(self):
+        pass
+
+    @abstractmethod
+    def slice_by_values(self):
+        pass
+
+    @abstractmethod
+    def get_cumulative_layer_thickness(self):
+        # Not sure if this should be here, potentially unsuitable with DiscreteData
+        # These kind of methods should go to a seperate layer_analysis module with
+        # functions to cover such analyses
+        pass
+
+    @abstractmethod
+    def get_layer_top(self):
+        # These kind of methods should go to a seperate layer_analysis module with
+        # functions to cover such analyses
+        pass
+
+    @abstractmethod
+    def to_multiblock(self):
+        pass
+
+    @abstractmethod
+    def to_vtm(self):
+        pass
+
+    @abstractmethod
+    def to_datafusiontools(self):
+        # supporting this is low priority, perhaps even deprecate
+        pass

geost-0.2.0/geost/analysis/__init__.py

@@ -0,0 +1,6 @@
+from geost.analysis.layer_analysis import (
+    cumulative_thickness,
+    find_top_sand,
+    layer_top,
+    top_of_sand,
+)

geost-0.2.0/geost/analysis/interpret_cpt.py

@@ -0,0 +1,50 @@
+import numpy as np
+
+# NOTE: temporary functions ripped from old code. To be made compatible with Geolib
+# functions
+
+
+def calc_ic(qc, rf) -> np.ndarray:
+    """
+    Calculate non-normalized IC values (I_SBT in Robertson 2010). The non-normalized
+    variant does not require calculations of stresses to normalize with and hence no PWP
+    data is required.
+
+    Please note the following when using non-normalized IC values:
+
+    "The non-normalized SBT index (ISBT) is essentially the same as the normalized SBTn
+    index (Ic) but only uses the basic CPT measurements. In general, the normalized Ic
+    provides more reliable identification of SBT than the non-normalized ISBT, but when
+    the insitu vertical effective stress is between 50 kPa to 150 kPa there is often
+    little difference between normalized and non-normalized SBT."
+
+    Parameters
+    ----------
+    qc : np.ndarray
+        Cone resistance values
+    rf : np.ndarray
+        Friction number
+
+    Returns
+    -------
+    np.ndarray
+        Non-normalized IC
+    """
+    ic = np.sqrt((3.47 - np.log10(qc / 0.1)) ** 2 + (np.log10(rf) + 1.22) ** 2)
+    return ic
+
+
+# TODO numpy searchsort
+def calc_lithology(ic, qc, rf) -> np.ndarray:
+    boundaries = [1.6, 2.0, 2.2, 2.6, 2.95, 3.6]
+    lith = np.full_like(ic, "NBE", dtype="<U3")
+    lith[ic < boundaries[0]] = "Z"
+    lith[(ic >= boundaries[0]) & (ic < boundaries[1])] = "Z"
+    lith[(ic >= boundaries[1]) & (ic < boundaries[2])] = "Z"
+    lith[(ic >= boundaries[2]) & (ic < boundaries[3])] = "Z"
+    lith[(ic >= boundaries[3]) & (ic < boundaries[4])] = "Kz"
+    lith[(ic >= boundaries[4]) & (ic < boundaries[5])] = "K"
+    lith[(ic >= boundaries[5]) & (rf > 8)] = "V"
+    lith[((rf > 5) & (qc < 1.5)) | (rf > 6)] = "V"
+    lith[(ic >= boundaries[5]) & (rf <= 8)] = "Kh"
+    return lith