坐标光栅图转换表

Question

我是一个完全的初学者，有一个我认为是相当简单的问题，但我有点卡住.

我有一个有三列的数据表：x、y和height值。我导入了csv文件并提取了值：

data <- read.csv("C:/Data/heights.csv")
mydata <- data[, c("x", "y", "height")]

然后，我需要将height值放入具有这些坐标的网格中，这是我作为光栅创建的：

grid <- raster(ncol=2001, nrow=2001, xmn=479975, xmx=500025, ymn=119975, ymx=140025)

Answer 1

我觉得你们很亲密。我建议你读一下光栅包装材料，格列奈特和描述。

你必须注意你的数据和它是如何组织的。一排排？列？栅格单元格值将按行填充。

对help ?raster中提供的示例使用simmilar

library(raster)
vals <- 1:100
r1 <- raster(nrows=10, ncols=10, xmn=479975, xmx=500025, ymn=119975, ymx=140025)
r1[] <- vals
plot(r1)

​

​

但是对于一个空间对象来说，过程可能是不同的。

library(sp)
x <- 1:10
y <- 1:10
df <- data.frame(expand.grid(x=x,y=y), v=1:100)
coordinates(df)=~x+y
gridded(df) = TRUE
r2 <- raster(df)
plot(r2)

​

​

我认为这可能是最好的方法，但最好了解一下您的数据以及它是如何组织的。

Answer 2

你可以试试这个：

library(raster)
x <- rasterFromXYZ(mydata)
extent(x) <- c(479975, 500025, 119975, 140025)