我有四点意见,代表四个空气质素监测站的位置。这四个观测站的参数如下。
Station Lat Long
1 Toronto Downtown 43.7 -79.4
2 Toronto East 43.7 -79.3
3 Toronto North 43.8 -79.5
4 Toronto West 43.7 -79.5 我想要创建一个具有以下参数的光栅:
dimensions : 14, 26, 364 (nrow, ncol, ncell)
resolution : 0.02, 0.02 (x, y)
extent : -79.64, -79.12, 43.58, 43.86 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +no_defs 我想要创建一个栅格,其中四个站被绘制在最能代表他们的位置的单元格中。但是光栅有与我提供的格式相匹配的。我该怎么做?
发布于 2022-01-20 16:16:35
还不清楚你的问题是什么,但我猜你在映射这方面有困难,因为你的数据是"lat/lon",而不是"lon/lat“。你可以这样解决:
library(terra)
r <- rast(nrow=24, ncol=26, xmin=-79.64, xmax=-79.12, ymin=43.58, ymax=43.86)
values(r) <- 1:ncell(r)
d <- data.frame(Station=c("D", "E", "N", "W"), Lat=c(43.7,43.7,43.8,43.7), Lon=c(-79.4,-79.3,-79.5,-79.5))
plot(r)
points(d[, c("Lon", "Lat")])https://stackoverflow.com/questions/70776533
复制相似问题