我正在尝试使用lidr对树木进行分割。我需要计算每棵树中有多少个点被分割。我有一个las变量,它有一个属性treeID。使用filterpoi,我可以得到单独的树id,看看有多少个点。我想提取点并为所有树id创建一个列表。我的想法是循环遍历treeID并为每个id获取点。但我似乎不能单独提取这些点。当我打印变量时,我得到了一些信息,我只需要"点"。如何提取
las <-segment_trees(las, algo) # segment point cloud
las <- filter_poi(las, treeID == 1)
> las
class : LAS (v1.4 format 6)
memory : 86.8 Kb
extent : 368194, 368203.5, 5807536, 5807544 (xmin, xmax, ymin, ymax)
coord. ref. : WGS 84 / UTM zone 33N
area : 67.9 m²
points : 1.1 thousand points
density : 16.11 points/m²
density : 9.96 pulses/m²
可以聚合data.table
las@data[, length(X), by = treeID]