polargrid-mapping 0.1.0__tar.gz

polargrid_mapping-0.1.0/LICENSE.txt

@@ -0,0 +1,19 @@
+Copyright (c) 2018 The Python Packaging Authority
+
+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.

polargrid_mapping-0.1.0/PKG-INFO

@@ -0,0 +1,45 @@
+Metadata-Version: 2.4
+Name: polargrid-mapping
+Version: 0.1.0
+Summary: Generate structured grids from scattered 2D points using convex hull layering
+Author-email: ArmanddeCacqueray <armanddecacqueray@sfr.fr>
+Requires-Python: >=3.8
+Description-Content-Type: text/markdown
+License-File: LICENSE.txt
+Requires-Dist: numpy
+Requires-Dist: scipy
+Dynamic: license-file
+
+\# polargrid
+
+
+
+Generate structured square grids from scattered 2D points using convex hull layering.
+
+
+
+\## Installation
+
+
+
+```bash
+
+pip install polargrid
+
+
+
+import numpy as np
+
+from polargrid import polargrid
+
+
+
+points = np.random.rand(1000, 2)
+
+grid = polargrid(points)
+
+
+
+points are now structured as a perfect grid, such that grid\[i, j] is the point with index i, j
+
+

polargrid_mapping-0.1.0/README.md

@@ -0,0 +1,33 @@
+\# polargrid
+
+
+
+Generate structured square grids from scattered 2D points using convex hull layering.
+
+
+
+\## Installation
+
+
+
+```bash
+
+pip install polargrid
+
+
+
+import numpy as np
+
+from polargrid import polargrid
+
+
+
+points = np.random.rand(1000, 2)
+
+grid = polargrid(points)
+
+
+
+points are now structured as a perfect grid, such that grid\[i, j] is the point with index i, j
+
+

polargrid_mapping-0.1.0/pyproject.toml

@@ -0,0 +1,18 @@
+[build-system]
+requires = ["setuptools>=61"]
+build-backend = "setuptools.build_meta"
+
+[project]
+name = "polargrid-mapping"
+version = "0.1.0"
+description = "Generate structured grids from scattered 2D points using convex hull layering"
+authors = [{ name = "ArmanddeCacqueray", email = "armanddecacqueray@sfr.fr" }]
+readme = "README.md"
+requires-python = ">=3.8"
+dependencies = [
+    "numpy",
+    "scipy"
+]
+
+[tool.setuptools.packages.find]
+where = ["src"]

polargrid_mapping-0.1.0/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

polargrid_mapping-0.1.0/src/polargrid/__init__.py

@@ -0,0 +1,4 @@
+from .core import polargrid
+
+__all__ = ["polargrid"]
+__version__ = "0.1.0"

polargrid_mapping-0.1.0/src/polargrid/core.py

@@ -0,0 +1,83 @@
+import numpy as np
+from scipy.spatial import ConvexHull
+from .utils import select, approx
+
+
+def polargrid(points: np.ndarray) -> np.ndarray:
+    """
+    Generate a structured square grid from scattered 2D points
+    using convex hull layering.
+
+    Parameters
+    ----------
+    points : ndarray of shape (N, 2)
+        Input 2D points.
+
+    Returns
+    -------
+    grid : ndarray of shape (n, n, 2)
+        Structured grid of points.
+    """
+    pts = points - points.mean(axis=0)
+
+    N = len(pts)
+    nn, n, layersizes = approx(N)
+    pts = pts[:nn]
+
+    x, y = pts[:, 0], pts[:, 1]
+    thetas = np.arctan2(y, x)
+    thetas0 = np.abs(thetas)
+
+    grid = np.zeros((n, n, 2))
+    remaining_idx = np.arange(nn)
+
+    for lsize in layersizes:
+        temp_idx = remaining_idx.copy()
+        layer = []
+
+        empty = False
+        while (len(layer) < 1.3 * lsize) and (not empty):
+            if len(temp_idx) > 2:
+                hull = ConvexHull(pts[temp_idx]).vertices
+            else:
+                hull = np.arange(len(temp_idx))
+                empty = True
+
+            hull_global = temp_idx[hull]
+            layer += list(hull_global)
+
+            mask = np.zeros(len(temp_idx), dtype=bool)
+            mask[hull] = True
+            temp_idx = temp_idx[~mask]
+
+        layer = np.array(layer)
+
+        # angular sort
+        layer = layer[np.argsort(thetas[layer])]
+
+        # thinning
+        thinning = select(thetas[layer], lsize)
+        layer = layer[thinning]
+
+        # rotate
+        start_idx = np.argmin(thetas0[layer])
+        layer = np.roll(layer, -start_idx)
+
+        # insert
+        if lsize != 1:
+            d = lsize // 4
+            indices = [i * d for i in range(4)] + [lsize]
+            pl = pts[layer]
+            segments = [pl[slice(indices[i], indices[i + 1])] for i in range(4)]
+
+            L = (n - (d + 1)) // 2
+            R = L + d
+
+            grid[L:R, L] = segments[0]
+            grid[R, L:R] = segments[1]
+            grid[R:L:-1, R] = segments[2]
+            grid[L, R:L:-1] = segments[3]
+
+        remaining_idx = np.setdiff1d(remaining_idx, layer)
+
+    return grid

polargrid_mapping-0.1.0/src/polargrid/utils.py

@@ -0,0 +1,63 @@
+import numpy as np
+import math
+
+
+def select(input_array, size):
+    input_array = np.asarray(input_array)
+    n = len(input_array)
+
+    order = np.argsort(input_array)
+    sorted_input = input_array[order]
+
+    target = np.linspace(sorted_input[0], sorted_input[-1], size)
+    j = np.searchsorted(sorted_input, target)
+
+    left = np.clip(j - 1, 0, n - 1)
+    right = np.clip(j, 0, n - 1)
+
+    dist_left = np.abs(sorted_input[left] - target)
+    dist_right = np.abs(sorted_input[right] - target)
+
+    selected = np.where(dist_right < dist_left, right, left)
+
+    used = np.zeros(n, dtype=bool)
+    result = np.empty(size, dtype=int)
+
+    for i, idx in enumerate(selected):
+        if not used[idx]:
+            used[idx] = True
+            result[i] = idx
+            continue
+
+        l = idx - 1
+        r = idx + 1
+
+        while True:
+            if l >= 0 and not used[l]:
+                used[l] = True
+                result[i] = l
+                break
+            if r < n and not used[r]:
+                used[r] = True
+                result[i] = r
+                break
+            l -= 1
+            r += 1
+
+    return np.sort(order[result])
+
+
+def approx(N):
+    n = int(math.sqrt(N))
+
+    is_even = (n % 2 == 0)
+    l = n // 2
+
+    idx = np.arange(1, l + 1)
+    if is_even:
+        layersizes = 8 * idx - 4
+    else:
+        layersizes = np.concatenate([[1], 8 * idx])
+
+    nn = n ** 2
+    return nn, n, layersizes[::-1]

polargrid_mapping-0.1.0/src/polargrid_mapping.egg-info/PKG-INFO

@@ -0,0 +1,45 @@
+Metadata-Version: 2.4
+Name: polargrid-mapping
+Version: 0.1.0
+Summary: Generate structured grids from scattered 2D points using convex hull layering
+Author-email: ArmanddeCacqueray <armanddecacqueray@sfr.fr>
+Requires-Python: >=3.8
+Description-Content-Type: text/markdown
+License-File: LICENSE.txt
+Requires-Dist: numpy
+Requires-Dist: scipy
+Dynamic: license-file
+
+\# polargrid
+
+
+
+Generate structured square grids from scattered 2D points using convex hull layering.
+
+
+
+\## Installation
+
+
+
+```bash
+
+pip install polargrid
+
+
+
+import numpy as np
+
+from polargrid import polargrid
+
+
+
+points = np.random.rand(1000, 2)
+
+grid = polargrid(points)
+
+
+
+points are now structured as a perfect grid, such that grid\[i, j] is the point with index i, j
+
+

polargrid_mapping-0.1.0/src/polargrid_mapping.egg-info/SOURCES.txt

@@ -0,0 +1,11 @@
+LICENSE.txt
+README.md
+pyproject.toml
+src/polargrid/__init__.py
+src/polargrid/core.py
+src/polargrid/utils.py
+src/polargrid_mapping.egg-info/PKG-INFO
+src/polargrid_mapping.egg-info/SOURCES.txt
+src/polargrid_mapping.egg-info/dependency_links.txt
+src/polargrid_mapping.egg-info/requires.txt
+src/polargrid_mapping.egg-info/top_level.txt

polargrid_mapping-0.1.0/src/polargrid_mapping.egg-info/dependency_links.txt

@@ -0,0 +1 @@
+

polargrid_mapping-0.1.0/src/polargrid_mapping.egg-info/requires.txt

@@ -0,0 +1,2 @@
+numpy
+scipy

polargrid_mapping-0.1.0/src/polargrid_mapping.egg-info/top_level.txt

@@ -0,0 +1 @@
+polargrid