绘制栅格数据和点

Question

我有一个点数据框并使用plot 绘制它。现在我想添加一个AnnualTemp 的RasterLayer。能否请您协助完成代码？

library(raster)

MinTemp_plots <- data.frame(
  Station = c( "Labasa", "Laucala", "Lautoka", "Levuka", "Matei", "Matuku",
               "Nabouwalu", "Nacocolevu",  "Nadi", "Nausori",  "Ono-I-Lau",  "Penang",
               "Savusavu",   "UduPoint",   "Viwa", "Yasawa"),
  
  Latitude = c(-16.43, -18.15, -17.6, -17.68, -16.69, -19.15, -16.98,
               -18.1,  -17.75, -18.03, -20.65, -17.37, -16.78, -16.11,
               -17.14, -16.78 ),
  
  Longitude = c(179.36, 178.45, 177.45, 178.83, 180, 179.76, 178.7, 177.55,
                177.45, 178.56, 178.7, 178.15, 179.34, 180, 176.93,
                177.5),
  
  AnnualTempMin = c(1.722, 1.711, 0.042, 0.135, 0.264, 0.276, 0.625, 1.215,  
                    1.522, 0.917,  0.617, 0.072, 0.509, 1.057, 1.201, 0.123))


plot(y= MinTemp_plots$Latitude, x= MinTemp_plots$Longitude)

r <- raster(xmn=1790828.61, xmx=2337149.40, ymn=3577110.39, ymx=4504717.19, res=100000)
crs(r) = crs("+init=epsg:3460")
    

Answer 1

我不完全确定你的问题是什么，但这里是：

为示例 RasterLayer 添加值以进行绘图

values(r) <- 1:ncell(r)

创建一个 SpatialPointsDataFrame

x <- MinTemp_plots
coordinates(x) <- ~ Longitude + Latitude
crs(x) <- "+proj=longlat"

将其转换为与栅格相同的坐标参考系

y <- spTransform(x, crs(r))

现在开始剧情

plot(r)
points(y)