给出LiDAR点的X Y Z坐标(取自itcSegment文档):
data(lasData)
检索最高点可以通过以下方式完成:
max(lasData$Z)
现在,我正在寻找一种方法来获得这个最高点的xy坐标。
这可能不是什么大问题,不过我试了很长时间......
答案 0 :(得分:0)
以下是使用矩阵的示例:
> m <- matrix(c(1:9),3,3)
> colnames(m) <- c('X','Y','Z')
> m
X Y Z
[1,] 1 4 7
[2,] 2 5 8
[3,] 3 6 9
> max(m[,'Z'])
[1] 9
> which(m[,'Z']==max(m[,'Z']))
[1] 3
> m[which(m[,'Z']==max(m[,'Z'])),]
X Y Z
3 6 9
答案 1 :(得分:0)
> data<-data.frame(x=c(2,5,4,6,2),y=c(34,55,23,3,6),z=c(54,23,6,11,93))
> data
x y z
1 2 34 54
2 5 55 23
3 4 23 6
4 6 3 11
5 2 6 93
> max(data$z)
[1] 93
> which(data$z == max(data$z))
[1] 5
答案 2 :(得分:0)
这可能会有所帮助
# data.frame
data(meuse.grid)
coordinates(meuse.grid) <- ~x+y
gridded(meuse.grid) <- TRUE
class(meuse.grid)
bbox(meuse.grid)
data(meuse)
meuse.xy = meuse[c("x", "y")]
coordinates(meuse.xy) <- ~x+y
class(meuse.xy)