py-pilecore 0.4.0__py3-none-any.whl → 0.4.1__py3-none-any.whl

{py_pilecore-0.4.0.dist-info → py_pilecore-0.4.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: py-pilecore
-Version: 0.4.0
+Version: 0.4.1
 Summary: Public python SDK for the CEMS PileCore web-API.
 License: MIT License
         
@@ -37,6 +37,7 @@ Requires-Dist: matplotlib <4,>=3.8
 Requires-Dist: tqdm[notebook] <5,>4
 Requires-Dist: natsort <9,>8
 Requires-Dist: shapely <3,>=2
+Requires-Dist: scipy <2,>=1
 Provides-Extra: docs
 Requires-Dist: Sphinx ==6.1.3 ; extra == 'docs'
 Requires-Dist: sphinx-autodoc-typehints ==1.22 ; extra == 'docs'

{py_pilecore-0.4.0.dist-info → py_pilecore-0.4.1.dist-info}/RECORD

@@ -1,24 +1,24 @@
 pypilecore/__init__.py,sha256=T2Uuao6fboVhrzKI_Sa7jXbcNAUl_PKk9BFdsozOG98,78
-pypilecore/_version.py,sha256=c8cd26QQ6zmNw4FZZlzQGhINwVjZ_hWaE2xv6ELt5Nw,175
+pypilecore/_version.py,sha256=0_HYvCjqXJ8bLPFLe6E2tkYlW15w7EOTtTD9FQ5jQGM,175
 pypilecore/api.py,sha256=dXgjfttEUcQHUCMkTX9eYip0dd3bf_ppE6cdx87ouHY,6822
 pypilecore/exceptions.py,sha256=-MZOfsxyHLCI0k1-wZFfVsMxc1lya5buuhLks5rxlCo,89
 pypilecore/plot_utils.py,sha256=rK5_067-4-x7LzZgt_t6ahcGrZInxNrqHqsy0RzCnq8,954
 pypilecore/utils.py,sha256=ib9LgJBIgWukL7zd_Zk1LP27UTMIZTRJ4RBB6ubn97o,1186
 pypilecore/input/__init__.py,sha256=tlmThdPtO8e6L6pqxuRQ7EOHRxYwuIcaNNGlZyAnzig,606
-pypilecore/input/grouper_properties.py,sha256=AqZHvg-CXcoTW7u9wgyqSm1oZLBc7sSoTawxGK69144,10895
-pypilecore/input/multi_cpt.py,sha256=1gEK-Imdhr0xOfcn50lq1OwYvjFs2NMf51LDaVddGas,19539
-pypilecore/input/pile_properties.py,sha256=mK8_VRy-tHZCKFHcv8ftMlCEs8KEsZpG9eGLq13hYPA,8400
-pypilecore/input/soil_properties.py,sha256=PlAutpueCzd5WHkEau16gVai7CT2j3kEAKVHxllMRz8,7688
+pypilecore/input/grouper_properties.py,sha256=qOc-SNYSPdOcZZ9DsRE-sBNVOKHbUD044EWRsWDu_kc,10893
+pypilecore/input/multi_cpt.py,sha256=Ua5IN9-m1jgCwa5xFNAYrMcxkVe2L1LIW907CWLJBW0,19577
+pypilecore/input/pile_properties.py,sha256=6O7iH-f7_G9PZbBsbzSWT_zUF2RtwYLfF5GRtUFvhHk,8404
+pypilecore/input/soil_properties.py,sha256=oqZWeSvhGyIKl4j0ZlGr923qyxkBjYxbVKBckG2rLAw,7687
 pypilecore/results/__init__.py,sha256=ggTI2QzILhX_oNx1YMOih6IVCkBFg8I5-Jyn9Sw3_h0,389
-pypilecore/results/grouper_result.py,sha256=tJl6NEdaFQKgx4W62P3dO5dNlOtf5OqZISTfNY35kSs,30903
+pypilecore/results/grouper_result.py,sha256=wQxWSJRrgh7aO5pglrCNTPvCNdERBMDXAlm6YZE3v-s,30927
 pypilecore/results/load_settlement.py,sha256=EbfTrSvH_g96KE-x8ZjmO8D0mt5KFaQ_-AR8u4blLsU,9752
 pypilecore/results/multi_cpt_results.py,sha256=dzNpP3hu0Sl00gNPHT37DcQqLbfp8IYci46QLSycyqA,25398
 pypilecore/results/pile_properties.py,sha256=z1R5UNoYjBem2rS8LMJ_ye2J5ejDoSilU5fCjUrNdUg,28592
-pypilecore/results/post_processing.py,sha256=ddIcbb3j5tva17F7hnkspDB6nd-D2_qlNmuJP_0vq10,15580
+pypilecore/results/post_processing.py,sha256=UWXcdff5dhPFDwzKbVIayEIp3HX6pxW8oQR3Z7AHPn0,22262
 pypilecore/results/single_cpt_results.py,sha256=irNqsL_cBKUh-QEmWC216KZjkv9MAZkaiLTvDDlEt20,17448
 pypilecore/results/soil_properties.py,sha256=Mcwz23PcaIW1xKLabMUSfIeKCTBli7ouuZPfT5dfrb4,20722
-py_pilecore-0.4.0.dist-info/LICENSE,sha256=3OCAZXffN0Bettjeya8uF_ZYegyvvCfH1WUt6CrHb_0,1061
-py_pilecore-0.4.0.dist-info/METADATA,sha256=QJEV5dvj3JqCrudtj9XSefcu0wGqr4TXjIjsJJXgtz8,5806
-py_pilecore-0.4.0.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
-py_pilecore-0.4.0.dist-info/top_level.txt,sha256=7BKIWZuSkbQtJ0ho5P1JvcaEbHzqADCcBuOduZmIaiI,11
-py_pilecore-0.4.0.dist-info/RECORD,,
+py_pilecore-0.4.1.dist-info/LICENSE,sha256=3OCAZXffN0Bettjeya8uF_ZYegyvvCfH1WUt6CrHb_0,1061
+py_pilecore-0.4.1.dist-info/METADATA,sha256=oQzHxa-c_UofjtUdg5FCaFWsqvanEgxCxuAkn3M_A3g,5834
+py_pilecore-0.4.1.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
+py_pilecore-0.4.1.dist-info/top_level.txt,sha256=7BKIWZuSkbQtJ0ho5P1JvcaEbHzqADCcBuOduZmIaiI,11
+py_pilecore-0.4.1.dist-info/RECORD,,

pypilecore/_version.py

@@ -4,4 +4,4 @@ try:
     __version__ = version("py-pilecore")
 # during CI
 except PackageNotFoundError:
-    __version__ = "0.4.0"
+    __version__ = "0.4.1"

pypilecore/input/grouper_properties.py

@@ -43,7 +43,7 @@ def create_grouper_payload(
     "/grouper/group_cpts"
 
 
-    Notes
+    Note
     ------
     The grouper uses pile bearing capacity results calculated by PileCore or other software to
     form  subgroups of the total group of CPT’s belonging to this project.
@@ -56,7 +56,6 @@ def create_grouper_payload(
         - is spatially coherent, which means there are no other CPTs in between the members
           of the subgroup. (Spatial check)
 
-
     Additionally, centre to centre validation (include_centre_to_centre_check; NEN9997-1 3.2.3) can be
     added to the cluster method. This check adds restrictions to the maximum allowable R;c;cal outliers
     and makes sure that the suitable data density requirements for the subgroup are met, by checking the

pypilecore/input/multi_cpt.py

@@ -344,7 +344,7 @@ def create_multi_cpt_payload(
         pile_head_level_nap=pile_head_level_nap,
         stiff_construction=stiff_construction,
         rel_pile_load=relative_pile_load,
-        soil_load=soil_load_sls,
+        soil_load=soil_load_sls if soil_load_sls is not None else 0.0,
         excavation_param_t=excavation_param_t,
         use_almere_rules=use_almere_rules,
         gamma_f_nk=gamma_f_nk,

pypilecore/input/pile_properties.py

@@ -29,8 +29,9 @@ def create_pile_properties_payload(
     Creates a dictionary with the `pile_properties` payload content for the PileCore
     endpoints.
 
-    Note that
-    the dictionary should be converted to a jsonifyable message before it can be passed
+    Note
+    ----
+    The dictionary should be converted to a jsonifyable message before it can be passed
     to a `requests` call directly, for instance with
     `nuclei.client.utils.python_types_to_message()`.
 

pypilecore/input/soil_properties.py

@@ -35,7 +35,7 @@ def create_soil_properties_payload(
     Creates a dictionary with the `soil_properties` payload content for the PileCore
     endpoints.
 
-    Notes
+    Note
     ------
     the dictionary should be converted to a jsonifyable message before it can be passed
     to a `requests` call directly, for instance with

pypilecore/results/grouper_result.py

@@ -211,7 +211,7 @@ class SingleClusterData:
         """
         Plot the bearing capacity and variation coefficient in a subplot
 
-        Notes
+        Note
         ------
         For the `Net bearing capacity` subplot there are two colors plotted:
          - orange:  conservative bearing capacity
@@ -251,7 +251,7 @@ class SingleClusterData:
         """
         Plot the spacing checks in a subplot
 
-        Notes
+        Note
         ------
         For the `spacing` subplot there are two colors plotted:
          - red:     invalid spacing
@@ -308,7 +308,7 @@ class SingleClusterData:
         """
         Plot the xi factor in a subplot
 
-        Notes
+        Note
         ------
         For the `xi factor` subplot there are two colors plotted:
          - olive:   xi3
@@ -499,7 +499,7 @@ class SingleClusterResult:
             - xi factor
             - centre to centre validation
 
-        Notes
+        Note
         ------
         For the `Net bearing capacity` subplot there are two colors plotted:
          - orange:  conservative bearing capacity
@@ -775,6 +775,9 @@ class GrouperResults:
     ) -> plt.Figure:
         """
         Plot a summary of the valid subgroups.
+
+        Note
+        -----
         Plot contains the:
 
             - cpts within a subgroup

pypilecore/results/post_processing.py

@@ -9,8 +9,10 @@ import numpy as np
 import pandas as pd
 from matplotlib import pyplot as plt
 from matplotlib.axes import Axes
+from matplotlib.collections import PatchCollection
 from matplotlib.figure import Figure
 from numpy.typing import NDArray
+from scipy.spatial import Delaunay, Voronoi, voronoi_plot_2d
 
 from pypilecore.results.soil_properties import SoilProperties, get_soil_layer_handles
 
@@ -94,7 +96,7 @@ class MaxBearingResult:
         The object with soil properties
     pile_head_level_nap
         The elevation of the pile-head, in [m] w.r.t. NAP.
-    results_table
+    table
         The object with CPT results.
     """
 
@@ -371,6 +373,59 @@ class MaxBearingResults:
 
         return cpt_results_df
 
+    @lru_cache()
+    def triangulation(self, pile_tip_level_nap: float) -> List[Dict[str, list]]:
+        """
+        Delaunay tessellation based on the CPT location
+
+        Returns
+        -------
+        collection: List
+            A list of dictionaries containing the tessellation
+            geometry and corresponding cpt names:
+
+                - geometry: List[Tuple[float, float]]
+                - test_id: List[str]
+
+        """
+        _lookup = {
+            (point.soil_properties.x, point.soil_properties.y): key
+            for key, point in self.cpt_results_dict.items()
+        }
+        # select point with valid bearing capacity at pile tip level
+        _points = (
+            self.to_pandas()
+            .loc[
+                (self.to_pandas()["pile_tip_level_nap"] == pile_tip_level_nap)
+                & (~pd.isna(self.to_pandas()["R_c_d_net"])),
+                ["x", "y"],
+            ]
+            .to_numpy()
+            .tolist()
+        )
+
+        # check if enough points Delaunay
+        if len(_points) < 4:
+            raise ValueError(
+                "Not enough points at this pile tip level to construct "
+                "the delaunay tessellation based on the CPT location."
+            )
+        tri = Delaunay(
+            _points,
+            incremental=False,
+            furthest_site=False,
+            qhull_options="Qbb",
+        )
+        geometries = np.array(_points)[tri.simplices]
+
+        return [
+            {
+                "geometry": geometry.tolist(),
+                "test_id": [_lookup[(xy[0], xy[1])] for xy in geometry],
+            }
+            for geometry in geometries
+        ]
+
     def plot(
         self,
         projection: Optional[Literal["3d"]] = "3d",
@@ -413,20 +468,20 @@ class MaxBearingResults:
         kwargs_subplot.update(kwargs)
         fig = plt.figure(**kwargs_subplot)
         axes = fig.add_subplot(projection=projection)
+        df = self.to_pandas().dropna()
         # create color list based on hue option
         if hue == "category":
             colors = [
-                "red" if var < pile_load_uls else "green"
-                for var in self.to_pandas()["R_c_d_net"]
+                "red" if var < pile_load_uls else "green" for var in df["R_c_d_net"]
             ]
         else:
-            colors = self.to_pandas()["R_c_d_net"].tolist()
+            colors = df["R_c_d_net"].tolist()
         # create scatter plot
         if projection == "3d":
             cmap = axes.scatter(
-                self.to_pandas()["x"],
-                self.to_pandas()["y"],
-                self.to_pandas()["pile_tip_level_nap"],
+                df["x"],
+                df["y"],
+                df["pile_tip_level_nap"],
                 c=colors,
             )
             axes.set_xlabel("X")
@@ -444,12 +499,13 @@ class MaxBearingResults:
                 )
         else:
             cmap = axes.scatter(
-                self.to_pandas()["test_id"],
-                self.to_pandas()["pile_tip_level_nap"],
+                df["test_id"],
+                df["pile_tip_level_nap"],
                 c=colors,
             )
             axes.set_ylabel("Z [m w.r.t NAP]")
             axes.tick_params(axis="x", labelrotation=90)
+            axes.grid()
 
         if hue == "category":
             fig.legend(
@@ -475,3 +531,149 @@ class MaxBearingResults:
             fig.colorbar(cmap, orientation="vertical", label="$R_{c;d;net}$ [kN]")
 
         return fig
+
+    def map(
+        self,
+        pile_tip_level_nap: float,
+        pile_load_uls: float = 100,
+        show_delaunay_vertices: bool = True,
+        show_voronoi_vertices: bool = False,
+        figsize: Tuple[int, int] | None = None,
+        **kwargs: Any,
+    ) -> plt.Figure:
+        """
+        Plot a map of the valid ULS load for a given depth.
+
+        Note
+        ------
+        Based on the Delaunay methode a tessellation is created with
+        the location of the CPT's. Each triangle is then colored according to
+        the bearing capacity of the CPT its based on. If any of the CPT does
+        not meet the required capacity the triangle becomes also invalid.
+
+        Warnings
+        --------
+        Please note that this map indication of valid ULS zones is intended as a visual aid to help
+        the geotechnical engineer. It does not necessarily comply with the NEN 9997-1+C2:2017 since the NEN is open
+        to interpretation. It is therefore that the interpretation provided by this methode must be carefully
+        validated by a geotechnical engineer.
+
+        Parameters
+        ----------
+        pile_tip_level_nap:
+            Pile tip level to generate map.
+        pile_load_uls
+            default is 100 kN
+            ULS load in kN. Used to determine if a pile tip level configuration is valid.
+        show_delaunay_vertices
+            default is True
+            Add delaunay vertices to the figure
+        show_voronoi_vertices
+            default is False
+            Add voronoi vertices to the figure
+        figsize:
+            Size of the activate figure, as the `plt.figure()` argument.
+        **kwargs:
+            All additional keyword arguments are passed to the `pyplot.subplots()` call.
+
+        Returns
+        -------
+        figure:
+            The `Figure` object where the data was plotted on.
+        """
+        kwargs_subplot = {
+            "figsize": figsize,
+            "tight_layout": True,
+        }
+
+        kwargs_subplot.update(kwargs)
+        fig, axes = plt.subplots(**kwargs_subplot)
+
+        # filter data
+        df = (
+            self.to_pandas()
+            .loc[self.to_pandas()["pile_tip_level_nap"] == pile_tip_level_nap]
+            .dropna()
+        )
+
+        if df.empty:
+            raise ValueError(
+                "Pile tip level is not valid pile tip level. "
+                "Please select one of the following pile tip level: "
+                f"[{(self.to_pandas()['pile_tip_level_nap']).unique()}]"
+            )
+
+        df["valid"] = [
+            False if var < pile_load_uls else True for var in df["R_c_d_net"]
+        ]
+
+        # iterate over geometry
+        if show_delaunay_vertices:
+            _patches = []
+            for tri in self.triangulation(pile_tip_level_nap):
+                color = (
+                    "green"
+                    if all(
+                        df.where(df["test_id"].isin(tri["test_id"])).dropna()["valid"]
+                    )
+                    else "red"
+                )
+                _patches.append(
+                    patches.Polygon(
+                        np.array(tri["geometry"]), facecolor=color, edgecolor="grey"
+                    )
+                )
+
+            collection = PatchCollection(_patches, match_original=True)
+            axes.add_collection(collection)
+
+        if show_voronoi_vertices:
+            points = [
+                (point.soil_properties.x, point.soil_properties.y)
+                for point in self.cpt_results_dict.values()
+            ]
+            vor = Voronoi(points)
+            voronoi_plot_2d(
+                vor,
+                show_vertices=False,
+                show_points=False,
+                ax=axes,
+                line_colors="black",
+                line_alpha=0.7,
+                line_width=0.1,
+                point_size=0.0,
+            )
+
+        # add the cpt names
+        axes.scatter(
+            df["x"],
+            df["y"],
+            c=["green" if val else "red" for val in df["valid"]],
+        )
+        for label, x, y in zip(df["test_id"], df["x"], df["y"]):
+            axes.annotate(label, xy=(x, y), xytext=(3, 3), textcoords="offset points")
+        axes.set_xlabel("X")
+        axes.set_ylabel("Y")
+        axes.ticklabel_format(useOffset=False)
+        fig.legend(
+            title="$R_{c;d;net}$ [kN]",
+            title_fontsize=18,
+            fontsize=15,
+            loc="lower right",
+            handles=[
+                patches.Patch(
+                    facecolor=color,
+                    label=label,
+                    alpha=0.9,
+                    linewidth=2,
+                    edgecolor="black",
+                )
+                for label, color in zip(
+                    [f">= {pile_load_uls}", f"< {pile_load_uls}"],
+                    ["green", "red"],
+                )
+            ],
+        )
+        axes.set_title(f"Pile tip level at: {pile_tip_level_nap} [m w.r.t NAP]")
+
+        return fig