swcgeom 0.13.2__py3-none-any.whl → 0.15.0__py3-none-any.whl

{swcgeom-0.13.2.dist-info → swcgeom-0.15.0.dist-info}/METADATA

@@ -1,9 +1,9 @@
 Metadata-Version: 2.1
 Name: swcgeom
-Version: 0.13.2
+Version: 0.15.0
 Summary: Neuron geometry library for swc format
 Author-email: yzx9 <yuan.zx@outlook.com>
-License: CC4.0 BY-NC-SA
+License: Apache-2.0
 Project-URL: repository, https://github.com/yzx9/swcgeom
 Keywords: neuronscience,neuron,neuroanatomy,neuron-morphology
 Requires-Python: >=3.10
@@ -15,6 +15,7 @@ Requires-Dist: numpy >=1.22.3
 Requires-Dist: pandas >=1.4.2
 Requires-Dist: pynrrd >=1.0.0
 Requires-Dist: scipy >=1.9.1
+Requires-Dist: sdflit >=0.2.1
 Requires-Dist: seaborn >=0.12.0
 Requires-Dist: tifffile >=2022.8.12
 Requires-Dist: typing-extensions >=4.4.0
@@ -30,8 +31,6 @@ Requires-Dist: urllib3 >=1.26.0 ; extra == 'all'
 
 # SWCGEOM
 
-<a rel="license" href="http://creativecommons.org/licenses/by-nc-sa/4.0/"><img alt="Creative Commons License" style="border-width:0" src="https://i.creativecommons.org/l/by-nc-sa/4.0/88x31.png" /></a><br />
-
 [![Run tests](https://github.com/yzx9/swcgeom/actions/workflows/test.yml/badge.svg)](https://github.com/yzx9/swcgeom/actions/workflows/test.yml)
 [![Release to GitHub](https://github.com/yzx9/swcgeom/actions/workflows/github-publish.yml/badge.svg)](https://github.com/yzx9/swcgeom/releases)
 [![Release to PyPI](https://github.com/yzx9/swcgeom/actions/workflows/pypi-publish.yml/badge.svg)](https://pypi.org/project/swcgeom/)
@@ -67,4 +66,6 @@ Static analysis don't support import hook used in editable install for [PEP660](
 
 ## LICENSE
 
-This work is licensed under a <a rel="license" href="http://creativecommons.org/licenses/by-nc-sa/4.0/">Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License</a>.
+This work is licensed under a <a rel="license" href="https://www.apache.org/licenses/">Apache-2.0</a>.
+
+Copyright (c) 2022-present, Zexin Yuan

{swcgeom-0.13.2.dist-info → swcgeom-0.15.0.dist-info}/RECORD

@@ -1,5 +1,5 @@
 swcgeom/__init__.py,sha256=z88Zwcjv-ii7c7dYd9QPg9XrUVorQjtrgGbQCsEnQhc,265
-swcgeom/_version.py,sha256=ByxP_AGaDBx1nhN9W4fte2GtMXAPh9eMZ4pSGlztQPc,413
+swcgeom/_version.py,sha256=oHv-EAjiXbJma3jZ0Tq6UPimiWYyyw2Ao9S8zdq9uWs,413
 swcgeom/analysis/__init__.py,sha256=esL_poW8u-_bmp7vR9ldcumX3_xodtaVfM1USqxQo5w,377
 swcgeom/analysis/branch_features.py,sha256=s6PTMwwvxrVtZXRZlQbUSIw4M9-1IG63kf-Nxc0tMB0,1958
 swcgeom/analysis/feature_extractor.py,sha256=coe07_bJau96BkimcnXzuf4KqjY5_QRLwqaumFsu_tQ,14031
@@ -8,55 +8,55 @@ swcgeom/analysis/path_features.py,sha256=iE21HBxAoGLxk_qK7MBwQhyUOBqNPcnk4urVHr9
 swcgeom/analysis/sholl.py,sha256=9hSW8rZm1gvSIgcEZg8IVPT8kzBgBfwqbwP4E8R7L44,6390
 swcgeom/analysis/trunk.py,sha256=L2tjUIUmrRQpah_W3ZETGWd16bDXJ5F8Sk2XBNGms0Q,5558
 swcgeom/analysis/visualization.py,sha256=mKOpzTPkLpr1ggGL1MZPZRTG92GEg4idLT4eN5z5KOs,5654
-swcgeom/analysis/volume.py,sha256=jC_HFpZ41D28JnSI5tjp5_zVkKAG47A7Z7WU9xmsezo,1918
+swcgeom/analysis/volume.py,sha256=nWPR7wGOk3Wl5eh97YMws0X-2jk8K7lmFp4-03wL3lY,4628
 swcgeom/core/__init__.py,sha256=ZUudZavxAIUU6Q0lBHrQ4ybmL5lBfvzyYsTtpuih9wg,332
 swcgeom/core/branch.py,sha256=uuJCxaByRu-OdDZVWEffSFcmZWY-6ZWUhHN1M2Awj1s,3980
-swcgeom/core/branch_tree.py,sha256=sN0viBVg5A4r9dMCkGNAaVttrdR4bEoZZBbHZFKdXj4,1892
+swcgeom/core/branch_tree.py,sha256=Ece6q1VNCRLLMj29N_MjXmmlHT8h4tpWCuDE0uSgKJo,1873
 swcgeom/core/node.py,sha256=HvfgsW4WU01hkRIPci8KF4bQMAkwtAxOGfUL4yUbuBs,3623
 swcgeom/core/path.py,sha256=CsEelHiDR0JPBP1dTvoCSRvX3kBlZxkQilosnncV4nQ,4188
-swcgeom/core/population.py,sha256=YZLAtkZKJeErOsE5MwhqieGjoMa7sWjFl5xxmXVPTco,8786
+swcgeom/core/population.py,sha256=I9xSeGUveLhxkOg_5FWLcvf4yJVJ6j9_n03Hna3y_6w,8790
 swcgeom/core/segment.py,sha256=yabRdFj7KlkJP4V67jAlCIRzpHduNnq3bRBIRMuANfA,3158
 swcgeom/core/swc.py,sha256=lSYxAa25l6O8WZ9JtSSET-RZMr6EA1Tq_aXL_x0H9Rc,6795
-swcgeom/core/tree.py,sha256=oAx31r2jFzxSW7vnM17m8sMVVI3_xzxRyDQh19MPTlQ,12126
-swcgeom/core/tree_utils.py,sha256=xPy9b3MO43QR7VSvJvrwioXyQLGAPLotQjekn_-UiLg,7308
-swcgeom/core/tree_utils_impl.py,sha256=5Cb63ziVVLADnkjfuq1T-ePw2TQQ5TKk4gcPZR6f_wY,1123
+swcgeom/core/tree.py,sha256=K2k7o8OZ9r2YeMmIebT8_0qvkwJYtLLzFHECz3fmmts,12365
+swcgeom/core/tree_utils.py,sha256=3aCHghny5Z727sxkt6P8E2MMr34vK6irfPCelMn3vk4,7681
+swcgeom/core/tree_utils_impl.py,sha256=kN2ByjqqQtZUfmC_ac25tXOaE-CMiV2lP58VxFphLEU,1616
 swcgeom/core/swc_utils/__init__.py,sha256=qghRxjtzvq5KKfN4HhvLpZNsGPfZQu-Jj2x62_5-TbQ,575
-swcgeom/core/swc_utils/assembler.py,sha256=RoMJ3RjLC4O7mk62QxXVTQ5SUHagrFmpEw3nOnnqeJo,4563
+swcgeom/core/swc_utils/assembler.py,sha256=_ByaVFc61rfCS2p9QUw4g40uF4pZ6NJaDc--TcV4jWo,642
 swcgeom/core/swc_utils/base.py,sha256=huVxjuMLlTHbEb-KSEFDLgU0Ss3723t2Gr4Z_gQtl00,4737
 swcgeom/core/swc_utils/checker.py,sha256=E72GtLZ_1IqQQ7aWQGs0dZ3Z609__bw3EYQqeWrk-EI,2657
 swcgeom/core/swc_utils/io.py,sha256=6_--Qoe8kDja4PWsjwqRAvPJZNMFILFgauHaeWeGikU,6444
 swcgeom/core/swc_utils/normalizer.py,sha256=_Ysi8bSJ2JBnIGB8o6BvAg2mcz6xuJp9rgNLZqpLuR8,5083
-swcgeom/core/swc_utils/subtree.py,sha256=bd4XOLmRDfQSn_ktfQM3Hn8ONpCuZ_TdTWhE9-7QXW4,1999
+swcgeom/core/swc_utils/subtree.py,sha256=43QITYvgXu3b_kfIod2Irrj3dSfrA-gTFev5VxzRafI,1995
 swcgeom/images/__init__.py,sha256=QBP1ZGGo2nWAcV7Krz-vbvW_jN4ChqXrrpoScXcUURs,96
 swcgeom/images/augmentation.py,sha256=v9zluYXmBEbafaDBTpvJovi4_KWJmHZZSvcYHzG0oWo,4099
-swcgeom/images/folder.py,sha256=urh60ITreGgEwSbCbgexJFM8_-C9WLAQ9jUN50sox10,4034
-swcgeom/images/io.py,sha256=VZWBq-_TMrKKdMaqpbZWZP6fCGIxWdhdIfb0ePcDy-4,20272
+swcgeom/images/folder.py,sha256=2REkrdNghLm1z8kZ2PDVvtsupzog8kARkeMjLuLiLFo,4955
+swcgeom/images/io.py,sha256=jUyKjtau8_5V-PN1kRsPHtP1OtueJi1zEf5-7PZ_tG8,21226
 swcgeom/transforms/__init__.py,sha256=Mi2mOgkQ50JbZ9LmgXgJIuAA5eio67oe2AOs2CCCxTQ,463
-swcgeom/transforms/base.py,sha256=22kkpQ11YCrH1SWWMHVqYjHp6cAlhslhdi1wKgAnT_s,3416
-swcgeom/transforms/branch.py,sha256=mdrTC9T4zENS9DdwlFaoFHcWDuFAgKgiNicJZ1gLsOc,5445
-swcgeom/transforms/geometry.py,sha256=zPmEf6ApJKzeZBbZ7qJYZXpY3-mrABN2dr7AUNqENiQ,6619
-swcgeom/transforms/image_stack.py,sha256=DSIkJzK_SxLB4bvwYBt0ACW5pbB_NJLT8raq2tjOE6E,7235
-swcgeom/transforms/images.py,sha256=tUl3dtqWfrhWR7WrOUc6LzWHnOzz3fwEgV3mMLRhn0c,908
-swcgeom/transforms/mst.py,sha256=B6ZRqeEfCHLHFu0tdOsOgplbYHra4VPYjnneXMU-cNY,6413
+swcgeom/transforms/base.py,sha256=gN5Iqi-OHkYrsjllSOdxI6Yzav3jJGoi6kUPy-38FAs,4101
+swcgeom/transforms/branch.py,sha256=R0rVti--u70IiUKyHSx6MsDYJyy6zSCf18Uia2Cmh28,5410
+swcgeom/transforms/geometry.py,sha256=XR73fO_8T7otUFIllqKOWW0OnrsXBc7yA01oDT99yMc,7385
+swcgeom/transforms/image_stack.py,sha256=KhZ26Ps88jk_7NkI9dkS2E7NZXUvMaN3Ln9WzJ-vPu0,5823
+swcgeom/transforms/images.py,sha256=l5Hx8x27zoClUz3s11j2oj0Lt9ROh5JJfjoIU_vd3H8,898
+swcgeom/transforms/mst.py,sha256=ceL_EWpCtoSy9zYD6wPP5zO68eNTrDl9JExoPLcE_Lw,6236
 swcgeom/transforms/path.py,sha256=Gk2iunGQMX7vE83bdo8xoDO-KAT1Vvep0iZs7oFLzFQ,1089
-swcgeom/transforms/population.py,sha256=ZrKfMAMx4l729f-JLgw0dnGIPtPUoV0ZZoNNyA5cBw8,826
-swcgeom/transforms/tree.py,sha256=Q5OnSti0ZeTNb-WpA_UZsDN7dhLN8YweEF_Siyrj66c,6420
-swcgeom/transforms/tree_assembler.py,sha256=UZ9OUg1bQNecYIY_7ippg3S8gpuoi617ZUE0jg6BrQE,3177
-swcgeom/utils/__init__.py,sha256=VeHspOo13k0Nxr4sZY5lT0amyaAE6Zfz96TF7sWEHp4,425
+swcgeom/transforms/population.py,sha256=EmZ6ntuOKe8mXJxMW7nCUA-w2DVlEVe2n0IOVz49tCY,833
+swcgeom/transforms/tree.py,sha256=Q45sVA7yi8pINV7ENOWPKz6DnEM3U5uJxwkFKNFS5VE,6262
+swcgeom/transforms/tree_assembler.py,sha256=vi_X9CNo5IxHP5J7bRl2z91PWufU6HmYlz1iyfdPUxE,5121
+swcgeom/utils/__init__.py,sha256=QfezYuQzgmPziNiWz2T1h77V-gjpuBoUi3mC-K-PZlI,427
 swcgeom/utils/debug.py,sha256=qay2qJpViLX82mzxdndxQFn-pi1vaEj9CbLGuGt8Y9k,465
 swcgeom/utils/download.py,sha256=By2qZezo6h1Ke_4YpSIhDgcisOrpjVqRmNzbhynC2xs,2834
 swcgeom/utils/dsu.py,sha256=3aCbtpnl_D0OXnowTS8-kuwnCS4BKBYL5ECiFQ1fUW8,1435
 swcgeom/utils/ellipse.py,sha256=LB3q5CIy75GEUdTauIpKySwIHaDrwXzzkBhOCnjJ8Vw,3259
 swcgeom/utils/file.py,sha256=1hchQDsPgn-i-Vz5OQtcogxav_ajCQ_OaEZCLmqczRg,2515
-swcgeom/utils/geometry_object.py,sha256=ayw9hzrxJXSIcEfYVSOXvPaQl8J9HYCGotGOfuH5RKA,7436
 swcgeom/utils/neuromorpho.py,sha256=CDK2tUM2pNwHv_lEserHhQs_VlY3Rn557-jtV63EFlk,14420
 swcgeom/utils/numpy_helper.py,sha256=A-F-eFdGktCHVAQ_HcXiFB3Y1YhhSNfAmtOl8483Dvo,1292
 swcgeom/utils/renderer.py,sha256=xHVZ06Z1MeKBPC3nKzuwA1HryzR0ga79y6johZA4-q0,7290
-swcgeom/utils/sdf.py,sha256=cmkHwdnh33u1gAXlzol5u9ZI3_SdMoMiY4pIkOgEcK0,5206
-swcgeom/utils/solid_geometry.py,sha256=8ao2GAme8kX-gCjB7JSGe2iGhbkGIy3T4Bg4TxsQNhI,1758
+swcgeom/utils/sdf.py,sha256=K-LSnwwNsGF85GW1kNBaqXAAVesxZAEsrn0nvOA2LcA,10614
+swcgeom/utils/solid_geometry.py,sha256=TV02jhcoCLCqtYA9hfE50LFD_VRfixMiOSiHB5Jb2_U,2431
 swcgeom/utils/transforms.py,sha256=PmP5fL_iVguq4GR2aqXhM0TeCsiFVnrPZMZG6zLohrE,6983
-swcgeom-0.13.2.dist-info/LICENSE,sha256=aiTJ_1to1Xx6PaByy-pSXg42VYzE4FvF6GIt69WSDDI,202
-swcgeom-0.13.2.dist-info/METADATA,sha256=feoVNhkkkzcDdqP1bgZd5IGiSN6dC54Tug3yPSDz2Bg,2575
-swcgeom-0.13.2.dist-info/WHEEL,sha256=oiQVh_5PnQM0E3gPdiz09WCNmwiHDMaGer_elqB3coM,92
-swcgeom-0.13.2.dist-info/top_level.txt,sha256=hmLyUXWS61Gxl07haswFEKKefYPBVJYlUlol8ghNkjY,8
-swcgeom-0.13.2.dist-info/RECORD,,
+swcgeom/utils/volumetric_object.py,sha256=DVRGGmQrAL0oaW6hbNtp5TStbic9DfyJdCzsv2FNw2c,15134
+swcgeom-0.15.0.dist-info/LICENSE,sha256=JPtohhZ4XURqoKI0ZqnMYb7dobCOoZR_n5EpnaLTp3E,11344
+swcgeom-0.15.0.dist-info/METADATA,sha256=VJX8KlMJWpYCZ-bMUYy_kL05sAHXmvHLMqUepQaKNTQ,2347
+swcgeom-0.15.0.dist-info/WHEEL,sha256=oiQVh_5PnQM0E3gPdiz09WCNmwiHDMaGer_elqB3coM,92
+swcgeom-0.15.0.dist-info/top_level.txt,sha256=hmLyUXWS61Gxl07haswFEKKefYPBVJYlUlol8ghNkjY,8
+swcgeom-0.15.0.dist-info/RECORD,,

swcgeom/utils/geometry_object.py

@@ -1,255 +0,0 @@
-"""Geometry object."""
-
-from abc import ABC, abstractmethod
-from functools import lru_cache
-
-import numpy as np
-import numpy.typing as npt
-
-from swcgeom.utils.solid_geometry import (
-    find_sphere_line_intersection,
-    find_unit_vector_on_plane,
-    project_point_on_line,
-)
-
-__all__ = ["GeomObject", "GeomSphere", "GeomFrustumCone"]
-
-eps = 1e-6
-
-
-class GeomObject(ABC):
-    """Geometry object."""
-
-    @abstractmethod
-    def get_volume(self) -> float:
-        """Get volume."""
-        raise NotImplementedError()
-
-    @abstractmethod
-    def get_intersect_volume(self, obj: "GeomObject") -> float:
-        """Get intersect volume.
-
-        Parameters
-        ----------
-        obj : GeometryObject
-            Another geometry object.
-
-        Returns
-        -------
-        volume : float
-            Intersect volume.
-        """
-        raise NotImplementedError()
-
-
-class GeomSphere(GeomObject):
-    """Geometry Sphere."""
-
-    def __init__(self, center: npt.ArrayLike, radius: float):
-        super().__init__()
-
-        self.center = np.array(center)
-        assert len(self.center) == 3
-
-        self.radius = radius
-
-    def get_volume(self) -> float:
-        return self.calc_volume(self.radius)
-
-    def get_volume_spherical_cap(self, h: float) -> float:
-        return self.calc_volume_spherical_cap(self.radius, h)
-
-    def get_intersect_volume_sphere(self, obj: "GeomSphere") -> float:
-        return self.calc_intersect_volume_sphere(self, obj)
-
-    def get_intersect_volume_sphere_frustum_cone(
-        self, frustum_cone: "GeomFrustumCone"
-    ) -> float:
-        return calc_intersect_volume_sphere_frustum_cone(self, frustum_cone)
-
-    def get_intersect_volume(self, obj: GeomObject) -> float:
-        if isinstance(obj, GeomSphere):
-            return self.get_intersect_volume_sphere(obj)
-
-        if isinstance(obj, GeomFrustumCone):
-            return self.get_intersect_volume_sphere_frustum_cone(obj)
-
-        classname = obj.__class__.__name__
-        raise NotImplementedError(f"unsupported geometry object: {classname}")
-
-    @staticmethod
-    def calc_volume(radius: float) -> float:
-        r"""Calculate volume of sphere.
-
-        \being{equation}
-        V = \frac{4}{3} * π * r^3
-        \end{equation}
-
-        Returns
-        -------
-        volume : float
-            Volume.
-        """
-        return 4 / 3 * np.pi * radius**3
-
-    @staticmethod
-    def calc_volume_spherical_cap(r: float, h: float) -> float:
-        r"""Calculate the volume of a spherical cap.
-
-        \being{equation}
-        V = π * h^2 * (3r - h) / 3
-        \end{equation}
-
-        Parameters
-        ----------
-        r : float
-            radius of the sphere
-        h : float
-            height of the spherical cap
-
-        Returns
-        -------
-        volume : float
-            volume of the spherical cap
-        """
-        return np.pi * h**2 * (3 * r - h) / 3
-
-    @classmethod
-    def calc_intersect_volume_sphere(
-        cls, obj1: "GeomSphere", obj2: "GeomSphere"
-    ) -> float:
-        r"""Calculate intersect volume of two spheres.
-
-        \being{equation}
-        V = \frac{\pi}{12d} * (r_1 + r_2 - d)^2 (d^2 + 2d r_1 - 3r_1^2 + 2d r_2 - 3r_2^2 + 6 r_1r_2)
-        \end{equation}
-
-        Returns
-        -------
-        volume : float
-            Intersect volume.
-        """
-
-        r1, r2 = obj1.radius, obj2.radius
-        d = np.linalg.norm(obj1.center - obj2.center).item()
-        if d > r1 + r2:
-            return 0
-
-        if d <= abs(r1 - r2):
-            return cls.calc_volume(min(r1, r2))
-
-        part1 = (np.pi / (12 * d)) * (r1 + r2 - d) ** 2
-        part2 = (
-            d**2 + 2 * d * r1 - 3 * r1**2 + 2 * d * r2 - 3 * r2**2 + 6 * r1 * r2
-        )
-        return part1 * part2
-
-
-class GeomFrustumCone(GeomObject):
-    """Geometry Frustum."""
-
-    def __init__(self, c1: npt.ArrayLike, r1: float, c2: npt.ArrayLike, r2: float):
-        super().__init__()
-
-        self.c1 = np.array(c1)
-        assert len(self.c1) == 3
-
-        self.c2 = np.array(c2)
-        assert len(self.c2) == 3
-
-        self.r1 = r1
-        self.r2 = r2
-
-    def height(self) -> float:
-        """Get height of frustum."""
-        return np.linalg.norm(self.c1 - self.c2).item()
-
-    def get_volume(self) -> float:
-        return self.calc_volume(self.r1, self.r2, self.height())
-
-    def get_intersect_volume_sphere(self, sphere: GeomSphere) -> float:
-        return calc_intersect_volume_sphere_frustum_cone(sphere, self)
-
-    def get_intersect_volume(self, obj: GeomObject) -> float:
-        if isinstance(obj, GeomSphere):
-            return self.get_intersect_volume_sphere(obj)
-
-        classname = obj.__class__.__name__
-        raise NotImplementedError(f"unsupported geometry object: {classname}")
-
-    @staticmethod
-    def calc_volume(r1: float, r2: float, height: float) -> float:
-        r"""Calculate volume of frustum.
-
-        \being{equation}
-        V = \frac{1}{3} * π * h * (r^2 + r * R + R^2)
-        \end{equation}
-
-        Returns
-        -------
-        volume : float
-            Volume.
-        """
-        return (1 / 3) * np.pi * height * (r1**2 + r1 * r2 + r2**2)
-
-
-@lru_cache
-def calc_intersect_volume_sphere_frustum_cone(
-    sphere: GeomSphere, frustum_cone: GeomFrustumCone
-) -> float:
-    r"""Calculate intersect volume of sphere and frustum cone.
-
-    Returns
-    -------
-    volume : float
-        Intersect volume.
-    """
-    h = frustum_cone.height()
-    c1, r1 = sphere.center, sphere.radius
-    if np.allclose(c1, frustum_cone.c1) and np.allclose(r1, frustum_cone.r1):
-        c2, r2 = frustum_cone.c2, frustum_cone.r2
-    elif np.allclose(c1, frustum_cone.c2) and np.allclose(r1, frustum_cone.r2):
-        c2, r2 = frustum_cone.c1, frustum_cone.r1
-    else:
-        raise NotImplementedError("unsupported to calculate intersect volume")
-
-    # Fast-Path: The surface of the frustum concentric with the sphere
-    # is the surface with smaller radius
-    if r2 - r1 >= -eps:  # r2 >= r1:
-        v_himisphere = GeomSphere.calc_volume_spherical_cap(r1, r1)
-        if h >= r1:
-            # The hemisphere is completely inside the frustum cone
-            return v_himisphere
-
-        # The frustum cone is lower than the hemisphere
-        v_cap = GeomSphere.calc_volume_spherical_cap(r1, r1 - h)
-        return v_himisphere - v_cap
-
-    up = (c2 - c1) / np.linalg.norm(c2 - c1)
-    v = find_unit_vector_on_plane(up)
-
-    intersections = find_sphere_line_intersection(c1, r1, c1 + r1 * v, c2 + r2 * v)
-    if len(intersections) == 0:
-        # Tricky case: Since the intersection point not found with
-        # numerical precision, we can simply assume that there are two
-        # intersection points and at the same position
-        intersections = [(0, c1 + r1 * v), (0, c1 + r1 * v)]
-    assert len(intersections) == 2
-    t, p = max(intersections, key=lambda x: x[0])
-
-    # Fast-Path: The frustum cone is completely inside the sphere
-    if t > 1 + eps:
-        return frustum_cone.get_volume()
-
-    M = project_point_on_line(c1, up, p)
-    h1 = np.linalg.norm(M - c1).item()
-    r3 = np.linalg.norm(M - p).item()
-    v_cap1 = GeomSphere.calc_volume_spherical_cap(r1, r1 - h1)
-    v_frustum = GeomFrustumCone.calc_volume(r1, r3, h1)
-
-    # Fast-Path: The frustum cone is higher than the sphere
-    if h >= r1:
-        return v_cap1 + v_frustum
-
-    v_cap2 = GeomSphere.calc_volume_spherical_cap(r1, r1 - h)
-    return v_cap1 + v_frustum - v_cap2

swcgeom-0.13.2.dist-info/LICENSE

@@ -1,3 +0,0 @@
-The code is released under the Creative Commons Attribution-NonCommercial-
-ShareAlike 4.0 International Public License for Noncommercial use only. Any
-commercial use should get formal permission first.