我是R的初学者,我正在尝试从R语言的薄板样条图中提取网格格式的数据。基本上,我有一个来自英国各地的点数据集,其中包含特定日期的经度、纬度和降雨量。使用以下代码,我可以将此数据绘制到英国地图上:
dat <- read.table("~jan1.csv", header=T, sep=",")
names(dat) <- c("gauge", "date", "station", "mm", "lat", "lon", "location", "county", "days")
library(fields)
quilt.plot(cbind(dat$lon,dat$lat),dat$mm)
world(add=TRUE)
到目前为止,一切都很好。我还可以使用以下方法执行薄板样条插值 (TPS):
fit <- Tps(cbind(dat$lon, dat$lat), dat$mm, scale.type="unscaled")
然后我可以在我选择的网格比例上绘制曲面图,例如:
surface (fit, nx=100, ny=100)
这有效地为我提供了分辨率为 100*100 的网格数据图。所以这是我的问题:
如何从该网格数据集中提取数据(即实际值)以放入 excel 或 .txt 等文件中?
如何更改网格大小,以便网格以我选择的步骤从特定的 x 值(和 y 值)开始?