有人对如何从R中的克里格数据导出值有任何建议吗?我已经将它们导出为光栅,但我需要实际的单元格值。我正在处理的数据是鱼类密度。我正在对数据进行克里格分析,然后将密度转换为丰度。但是为了转换为丰度,我需要.csv格式的克里格值。
下面是我正在使用的代码:
library(sp)
library(gstat)
library(raster)
library(automap)
library(ggplot2)
library(rgdal)
library(maptools)
fsite.fit=fit.variogram(fsite.vario, vgm(model="Sph",psill=fmy.psill,range=fmy.range,nugget=fmy.nugget),
fit.method=1)
fXloc = data.frame(Data.krig)$fXloc
fYloc = data.frame(Data.krig)$fYloc
fsite.chull = chull(fXloc,fYloc)
plot(fXloc,fYloc)
lines(fXloc[fsite.chull],fYloc[fsite.chull])
fsite.grid = polygrid(
xgrid=seq(min(fXloc),max(fXloc),length=100),
ygrid=seq(min(fYloc),max(fYloc),length=100),
cbind(
fXloc[fsite.chull],
fYloc[fsite.chull]))
names(fsite.grid)=c("fXloc","fYloc")
coordinates(fsite.grid)=c("fXloc","fYloc")
fsite.grid = as(fsite.grid, "SpatialPixels")
plot(fYloc~fXloc,cex=1.2,pch=20,col=2)
points(data.frame(fsite.grid)$fXloc,data.frame(fsite.grid)$fYloc,pch="+")
fsite.ok = krige(fADens~1, Data.krig, fsite.grid, fsite.fit)
现在我想导出fsite.grid中的fsite.ok,但我无法将它们连接到同一数据帧中。我不知道从这里到哪里去。如果你想让我添加其他内容,请好心地告诉我。此时,我不被允许共享我的数据。
这里是gstat
的一个最小、可重复的示例。看起来您对对象fsite.ok
感兴趣。可能是Spatial*DataFrame
。
示例代码
library(gstat)
library(sp)
data(meuse)
coordinates(meuse) = ~x+y
data(meuse.grid)
gridded(meuse.grid) = ~x+y
meuse.gstat <- gstat(id = "zinc", formula = zinc ~ 1, data = meuse,
nmax = 7, set = list(idp = .5))
meuse.gstat
z <- predict(meuse.gstat, meuse.grid)
现在你可以做了
d <- data.frame(z)
write.csv(d, "values.csv", row.names=FALSE)
所以在你的例子中,这将是
write.csv(data.frame(fsite.ok), "values.csv", row.names=FALSE)