时间戳

首页 > 解决方案 > 从 R 导出克里格值

r - 从 R 导出克里格值

问题描述

有没有人对如何从 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 但我无法将它们加入同一个数据框中。我不知道从这里去哪里。请让我知道您是否希望我添加任何其他内容。目前，我不允许分享我的数据。

标签： rkriging

解决方案

这是一个最小的、可重现的示例gstat。看起来您对 object 感兴趣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)

推荐阅读