update example

IGNF · May 29, 2024 · fb77f91 · fb77f91
1 parent 030654a
commit fb77f91
Showing 1 changed file with 5 additions and 2 deletions.
diff --git 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()