Fix #708

r-lidar · Aug 29, 2023 · b104500 · b104500
1 parent 107126d
commit b104500
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 0 deletions.
diff --git a/NEWS.md b/NEWS.md
@@ -9,6 +9,7 @@ If you are viewing this file on CRAN, please check [the latest news on GitHub](h
 - New: [#680](https://github.com/r-lidar/lidR/issues/680) it is now possible to use `plot(las, mapview = TRUE)`.
 - Fix: [#702](https://github.com/r-lidar/lidR/issues/702) `plot(las, breaks = <vector>)` now works with a vector of custom break points.
 - Fix: [#701](https://github.com/r-lidar/lidR/issues/701) `crown_metrics` now always remove invalid polygons if any.
+- Fix: [#708](https://github.com/r-lidar/lidR/issues/708) files read with `readMSLAS` can now be written properly with `writeLAS`.
 
 ## lidR v4.0.3 (Release date: 2023-03-09)
 

diff --git a/R/io_readXLAS.R b/R/io_readXLAS.R
@@ -74,6 +74,7 @@ readMSLAS = function(files1, files2, files3, select = "*", filter = "")
     las@data[["ScannerChannel"]][(n1+n2+1):npoints(las)] <- 3L
     las@header@PHB[["Point Data Format ID"]] <- 6L
     las@header@PHB[["Version Minor"]] <- 4L
+    las@header@PHB[["Header Size"]] <- 375L
     message("Multispectral data read from point format < 6 with no ScannerChannel attribute. The LAS object has been upgraded to LAS 1.4 prf 6 with a ScannerChannel.")
   }
   else