Tôi có thể chiếu lại tệp LAS trong LidR không

Question

Tôi đang làm việc với 2 tập dữ liệu LAS 2014 liền kề được cung cấp trong 2 dự báo:

Tuy nhiên, khi tôi nhập espg (las) thì câu trả lời của tôi là 0

Nếu tôi nhập tóm tắt (las) coord của mình. Ref: NA

Tôi muốn chiếu lại cả hai điều này lên NAD83 / UTM18 và xử lý chúng cùng nhau thông qua LiDR. Tôi không rõ liệu điều này có thể thực hiện được trong LiDR không? Làm thế nào tốt nhất để tiến hành? Tôi có 1000 ô vuông 1km để giải quyết.

Answer 1

Bạn có thể gán CRS cho đối tượng LAS của mình bằng epsg()<-

epsg(las) <- 12345

Bạn có thể chiếu lại với spTransform()nhưng chức năng này không tối ưu trong cách triển khai hiện tại. Nó thực sự biến đổi LASthành một SpatialPoints, sử dụng sp::spTransform()và cập nhật lại LAS. Tôi khuyên bạn nên sử dụng las2lasfrom LAStoolscho loại tác vụ xử lý tệp LAS này nhưng đối với tập dữ liệu nhỏ spTransformcó thể thuận tiện.

library(lidR)
LASfile <- system.file("extdata", "Megaplot.laz", package="lidR")
las = readLAS(LASfile)
las
#> class        : LAS (v1.2 format 1)
#> memory       : 6.2 Mb 
#> extent       : 684766.4, 684993.3, 5017773, 5018007 (xmin, xmax, ymin, ymax)
#> coord. ref.  : +proj=utm +zone=17 +datum=NAD83 +units=m +no_defs 
#> area         : 53112.69 m²
#> points       : 81.6 thousand points
#> density      : 1.54 points/m²
las2 = spTransform(las, sp::CRS(SRS_string = "EPSG:26918"))
las2
#> class        : LAS (v1.2 format 1)
#> memory       : 6.2 Mb 
#> extent       : 214261.7, 214504.9, 5021517, 5021767 (xmin, xmax, ymin, ymax)
#> coord. ref.  : +proj=utm +zone=18 +datum=NAD83 +units=m +no_defs 
#> area         : 53174.62 m²
#> points       : 81.6 thousand points
#> density      : 1.53 points/m²

---

Lưu ý - trước khi lidR 3.0.4chức năng có một chút lỗi. Nếu bạn kiểm tra, las2bạn sẽ thấy nó không còn tuân thủ nghiêm ngặt các thông số kỹ thuật. Nó đã được cố định trong lidR 3.0.4.

las_check(las2)