我正在尝试从经度、纬度和海拔数据点创建网格,但我遇到了很多困难。这里有更多细节。
我有经度、纬度和海拔数据点,当绘制不同的海拔高度时,这些数据点看起来很闲散:
我的最终目标是为之前的每个形状创建一个网格,并将这些网格用于一些克里金编码。
我阅读了一堆文档,但我无法弄清楚如何做到这一点。我想我可能必须从经度、纬度和海拔创建一个 shapefile,然后从那里创建网格,但我仍然可以弄清楚如何做到这一点。
任何事情都会有所帮助,尤其是一些实现这样的事情的编码示例!
谢谢。(这是上图的代码)
all_pit_data<- read.table("ALL_PIT_DATA.csv",header=TRUE,sep=",")
library(ggplot2)
all_pit_data$intervals <- cut_interval(all_pit_data$z1, n=9)
windows()
ggplot(all_pit_data, aes(x = x1, y = y1, color = z1)) +
geom_point(shape=1) +
facet_wrap ( ~ intervals, nrow=3) +
labs(x = "x", y = "y", colour = "Depth")enter code here