From fb77f919da89050f9263522ea82c4ce00939c25c Mon Sep 17 00:00:00 2001 From: Antoine Lavenant Date: Wed, 29 May 2024 17:23:40 +0200 Subject: [PATCH] update example --- macro/ex_filtering_points_with_add_attributs.py | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) diff --git a/macro/ex_filtering_points_with_add_attributs.py b/macro/ex_filtering_points_with_add_attributs.py index 8d9a9a4..c5700dc 100755 --- a/macro/ex_filtering_points_with_add_attributs.py +++ b/macro/ex_filtering_points_with_add_attributs.py @@ -65,11 +65,14 @@ def parse_args(): pipeline |= pdal.Filter.assign(value=["PT_GRID_DSM = 0 WHERE " + macro.build_condition("Classification", [2,3,4,5,9]) + " && PT_ON_BRIDGE==1"]) - # 4 - export du nuage et des DSM + ## 4 - point pour DTM servent au DSM également + pipeline |= pdal.Filter.assign(value=["PT_GRID_DSM = 1 WHERE PT_GRID_DTM==1"]) + + ## 5 - export du nuage et des DSM pipeline |= pdal.Writer.las(extra_dims="all", minor_version=4, dataformat_id=6, filename=args.output_las) pipeline |= pdal.Writer.gdal(gdaldriver="GTiff", output_type="max", resolution=2.0, filename=args.output_dtm, where="PT_GRID_DTM==1") - pipeline |= pdal.Writer.gdal(gdaldriver="GTiff", output_type="max", resolution=2.0, filename=args.output_dsm, where="PT_GRID_DTM==1 || PT_GRID_DSM==1") + pipeline |= pdal.Writer.gdal(gdaldriver="GTiff", output_type="max", resolution=2.0, filename=args.output_dsm, where="PT_GRID_DSM==1") pipeline.execute()