ign-pdal-tools 1.7.1__py3-none-any.whl → 1.7.3__py3-none-any.whl

{ign_pdal_tools-1.7.1.dist-info → ign_pdal_tools-1.7.3.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: ign-pdal-tools
-Version: 1.7.1
+Version: 1.7.3
 Summary: Library for common LAS files manipulation with PDAL
 Author-email: Guillaume Liegard <guillaume.liegard@ign.fr>
 Description-Content-Type: text/markdown
@@ -50,6 +50,18 @@ filename = ...
 las_infos.las_info_metadata(filename)
 ```
 
+## Point cloud infos
+
+Misc tools to get information on a point cloud (numpy array). Eg. get expected origin of a point cloud based on a square tiling:
+
+```python
+from pdaltools import pcd_infos
+
+points = ...
+pcd_infos.get_pointcloud_origin_from_tile_width(points, tile_width=1000)
+```
+
+
 ## Stitching
 
 * [las_clip.py](pdaltools/las_clip.py): crop a LAS file using 2d bounding box

{ign_pdal_tools-1.7.1.dist-info → ign_pdal_tools-1.7.3.dist-info}/RECORD

@@ -1,15 +1,16 @@
-pdaltools/_version.py,sha256=xwBg7KNoeyzwyonEgM9Mm-9CRDfrRb3V1tY-iPeL3Os,74
+pdaltools/_version.py,sha256=CIj6_M9SXht6CUH70oDuQQl4_AP2fggn-8tFZbGL0AY,74
 pdaltools/color.py,sha256=7U-SThIKqrfE1xXXnFqpbIhmZEqna29nRiyLW8l8Y1c,8075
 pdaltools/las_add_buffer.py,sha256=sBpTywlfsHHS8KuCUa-eydB2hylshEvjrMQt5TrqXb8,11275
 pdaltools/las_clip.py,sha256=GvEOYu8RXV68e35kU8i42GwSkbo4P9TvmS6rkrdPmFM,1034
 pdaltools/las_info.py,sha256=RE-UBdEUXqKvSrMV3mOlvE_16mhum7bw-p-ERu5bGOc,6979
 pdaltools/las_merge.py,sha256=tcFVueV9X9nNEaoAl5zCduY5DETlBg63MAgP2SuKiNo,4121
 pdaltools/las_remove_dimensions.py,sha256=0zhv9LBvlL69TLmXTJlRQcUBOaBmCRZEQU2Qadx27aM,1805
+pdaltools/pcd_info.py,sha256=d7ukCEdOFBQqaFp3HYM2MNbKDl1hdxUVHI9vSu9jN-M,1897
 pdaltools/replace_attribute_in_las.py,sha256=po1F-fi8s7iilqKWaryW4JRbsmdMOUe0yGvG3AEKxtk,4771
-pdaltools/standardize_format.py,sha256=KM_jC_aC9yLD5rrSUGgTwfyakbh86FXsAI-y8gokF4M,2883
+pdaltools/standardize_format.py,sha256=lUzrjWY22WfxT_2zLtxuXbmWgiAvjy9o7fBXWsAZmDg,3240
 pdaltools/unlock_file.py,sha256=pIThdWMNkTph0xgJVVRaM1o9aUMQhM6804PscScB3JI,1963
-ign_pdal_tools-1.7.1.dist-info/LICENSE.md,sha256=iVzCFZTUXeiqP8bP474iuWZiWO_kDCD4SPh1Wiw125Y,1120
-ign_pdal_tools-1.7.1.dist-info/METADATA,sha256=4xHIuVdU1-CJuSZZZhUIWUGg1rLzTtq9JV_VRsvrBPQ,4825
-ign_pdal_tools-1.7.1.dist-info/WHEEL,sha256=Z4pYXqR_rTB7OWNDYFOm1qRk0RX6GFP2o8LgvP453Hk,91
-ign_pdal_tools-1.7.1.dist-info/top_level.txt,sha256=KvGW0ZzqQbhCKzB5_Tp_buWMZyIgiO2M2krWF_ecOZc,10
-ign_pdal_tools-1.7.1.dist-info/RECORD,,
+ign_pdal_tools-1.7.3.dist-info/LICENSE.md,sha256=iVzCFZTUXeiqP8bP474iuWZiWO_kDCD4SPh1Wiw125Y,1120
+ign_pdal_tools-1.7.3.dist-info/METADATA,sha256=PfhYtt73bGGZxehInuRGdzwfhfW_1UxMO7K3CrPlDL0,5112
+ign_pdal_tools-1.7.3.dist-info/WHEEL,sha256=UvcQYKBHoFqaQd6LKyqHw9fxEolWLQnlzP0h_LgJAfI,91
+ign_pdal_tools-1.7.3.dist-info/top_level.txt,sha256=KvGW0ZzqQbhCKzB5_Tp_buWMZyIgiO2M2krWF_ecOZc,10
+ign_pdal_tools-1.7.3.dist-info/RECORD,,

{ign_pdal_tools-1.7.1.dist-info → ign_pdal_tools-1.7.3.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (70.3.0)
+Generator: setuptools (74.0.0)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

pdaltools/_version.py

@@ -1,4 +1,4 @@
-__version__ = "1.7.1"
+__version__ = "1.7.3"
 
 
 if __name__ == "__main__":

pdaltools/pcd_info.py

@@ -0,0 +1,46 @@
+"""Tools to get information from a point cloud (points as a numpy array)"""
+
+from typing import Tuple
+
+import numpy as np
+
+
+def get_pointcloud_origin_from_tile_width(
+    points: np.ndarray, tile_width: int = 1000, buffer_size: float = 0
+) -> Tuple[int, int]:
+    """Get point cloud theoretical origin (xmin, ymax) for a data that originates from a square tesselation/tiling
+    using the tesselation tile width only.
+
+    Edge values are supposed to be included in the tile
+
+
+    Args:
+        points (np.ndarray): numpy array with the tile points
+        tile_width (int, optional): Edge size of the square used for tiling. Defaults to 1000.
+        buffer_size (float, optional): Optional buffer around the tile. Defaults to 0.
+
+    Raises:
+        ValueError: Raise an error when the bounding box of the tile is not included in a tile
+
+    Returns:
+        Tuple[int, int]: (origin_x, origin_y) origin coordinates
+    """
+    # Extract coordinates xmin, xmax, ymin and ymax of the original tile without buffer
+    x_min, y_min = np.min(points[:, :2], axis=0) + buffer_size
+    x_max, y_max = np.max(points[:, :2], axis=0) - buffer_size
+
+    # Calculate the tiles to which x, y bounds belong
+    tile_x_min = np.floor(x_min / tile_width)
+    tile_x_max = np.floor(x_max / tile_width) if x_max % tile_width != 0 else np.floor(x_max / tile_width) - 1
+    tile_y_min = np.ceil(y_min / tile_width) if y_min % tile_width != 0 else np.floor(y_min / tile_width) + 1
+    tile_y_max = np.ceil(y_max / tile_width)
+
+    if not (tile_x_max - tile_x_min) and not (tile_y_max - tile_y_min):
+        origin_x = tile_x_min * tile_width
+        origin_y = tile_y_max * tile_width
+        return origin_x, origin_y
+    else:
+        raise ValueError(
+            f"Min values (x={x_min} and y={y_min}) do not belong to the same theoretical tile as"
+            f"max values (x={x_max} and y={y_max})."
+        )

pdaltools/standardize_format.py

@@ -7,6 +7,7 @@
     - precision
     - no extra-dims
 """
+
 import argparse
 import os
 import subprocess as sp
@@ -42,6 +43,14 @@ def parse_args():
         "--record_format", choices=[6, 8], type=int, help="Record format: 6 (no color) or 8 (4 color channels)"
     )
     parser.add_argument("--projection", default="EPSG:2154", type=str, help="Projection, eg. EPSG:2154")
+    parser.add_argument(
+        "--extra_dims",
+        default=[],
+        nargs="*",
+        type=str,
+        help="List of extra dims to keep in the output (default=[], use 'all' to keep all extra dims), "
+        "extra_dims must be specified with their type (see pdal.writers.las documentation, eg 'dim1=double')",
+    )
 
     return parser.parse_args()
 
@@ -86,5 +95,5 @@ def standardize(input_file: str, output_file: str, params_from_parser: Dict) ->
 
 if __name__ == "__main__":
     args = parse_args()
-    params_from_parser = dict(dataformat_id=args.record_format, a_srs=args.projection)
+    params_from_parser = dict(dataformat_id=args.record_format, a_srs=args.projection, extra_dims=args.extra_dims)
     standardize(args.input_file, args.output_file, params_from_parser)