将坐标系分配给我的 las 目录，但我的代码不起作用

Question

我有一个没有任何坐标系的 LAS 目录数据。如果没有坐标系，我可以轻松计算像素度量，但是当我分配坐标时它不起作用；我需要 x 和 y 作为结果表。你能帮我解决这个问题吗？

normalized1987<-readLAScatalog("normalized1987")
st_crs(normalized1987)<-4326
heigh<-pixel_metrics(normalized1987,func = .stdmetrics_z,res = 30)

我尝试了 LAS 检查，但收到此错误：

Checking overlapping tiles...Error in wk_handle.wk_wkb(wkb, s2_geography_writer(oriented = oriented,  : 
      Loop 0 is not valid: Edge 0 crosses edge 2
     In addition: Warning messages:
     1: In st_is_longlat(x) :
      bounding box has potentially an invalid value range for longlat data
      2: In st_is_longlat(x) :
       bounding box has potentially an invalid value range for longlat data

但过滤重复项也不起作用。

Answer 1

EPSG:4326

是经度/纬度 CRS。您的点云位于投影的 CRS 中。当尝试使用无效的投影坐标构建长/纬度栅格时，该错误来自

terra

包