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

Question

我有一个没有坐标系的 las 目录数据，没有坐标系我可以轻松计算像素指标，但是当我分配坐标时它不起作用，我需要 x 和 y 作为结果表，你可以吗帮我解决问题

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

我尝试了 las check，但收到此错误 检查重叠的图块...wk_handle.wk_wkb（wkb，s2_geography_writer（orient = oriented，： 循环 0 无效：边沿 0 与边沿 2 交叉 另外： 警告消息： 1：在 st_is_longlat(x) 中： 对于 longlat 数据，边界框可能具有无效的值范围 2：在 st_is_longlat(x) 中： 对于 longlat 数据，边界框可能具有无效的值范围 但过滤重复也不起作用

Answer 1

EPSG:4326

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

terra

包