[英]Plotting raster data and points
我有一個 dataframe 點並使用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")
我不完全確定你的問題是什么,但這里是:
將值添加到示例 RasterLayer 以進行繪圖
values(r) <- 1:ncell(r)
創建一個 SpatialPointsDataFrame
x <- MinTemp_plots
coordinates(x) <- ~ Longitude + Latitude
crs(x) <- "+proj=longlat"
將其轉換為與柵格相同的坐標參考系
y <- spTransform(x, crs(r))
現在 plot
plot(r)
points(y)
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.