2

我正在尝试从经度、纬度和海拔数据点创建网格,但我遇到了很多困难。这里有更多细节。

我有经度、纬度和海拔数据点,当绘制不同的海拔高度时,这些数据点看起来很闲散:

在此处输入图像描述

我的最终目标是为之前的每个形状创建一个网格,并将这些网格用于一些克里金编码。

我阅读了一堆文档,但我无法弄清楚如何做到这一点。我想我可能必须从经度、纬度和海拔创建一个 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
4

0 回答 0