如何在r的活页打印中设置用于打印的插值光栅投影？

我正在尝试将光栅图像添加到传单中。虽然我对插值不太满意，但我确实有一个可以通过正常的

plot

函数进行绘图的插值。我在下面展示的代码应该能够独立运行，因为我为输入数据提供了示例数据

ozone_df<-data.frame("Longitude"=runif(200, -112.245075*10000,-111.455581*10000)/10000)
ozone_df$Latitude<-runif(200, 40.063614*10000,40.827281*10000)/10000
ozone_df$Ozone<-runif(200, 0,115)

#create grid tick marks
small_grid_x = seq(-111.455581,-112.245075,length.out=500)
small_grid_y = seq(40.063614,40.827281,length.out=500)
#create grid nodes
krig_grid_small<-expand.grid(small_grid_x,small_grid_y)
coordinates(krig_grid_small) <- ~ Var1 + Var2
#create kriging fit and apply interpolation to grid
krig_fit_small<-fields::Krig(ozone_df[1:2],ozone_df$Ozone)
ozone_krig_small<-raster::interpolate(grid_raster_small, krig_fit_small)
crs(ozone_krig_small) <-CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")
#plot output raster
plot(ozone_krig_small)

leaflet() %>% 
  addRasterImage(ozone_krig_small, project=T)%>% 
  addTiles() %>% 
  setView(lng = -111.941004,  40.610497, zoom = 10) %>% 
  addMiniMap()

---

如果你更新你的软件包（或者至少是光栅、sp和rgdal），这个错误应该会消失——使用R>=4也很好

这里有一个版本，有一些修正，对我来说是有效的

library(raster)
library(leaflet)
library(fields)
ozone_df<-data.frame("Longitude"=runif(200, -112.245075*10000,-111.455581*10000)/10000)
ozone_df$Latitude<-runif(200, 40.063614*10000,40.827281*10000)/10000
ozone_df$Ozone<-runif(200, 0,115)

#create grid tick marks
small_grid_x = seq(-111.455581,-112.245075,length.out=500)
small_grid_y = seq(40.063614,40.827281,length.out=500)

#create grid nodes
krig_grid_small <- expand.grid(small_grid_x,small_grid_y)
coordinates(krig_grid_small) <- ~ Var1 + Var2

#create kriging fit and apply interpolation to grid
krig_fit_small<-fields::Krig(ozone_df[1:2],ozone_df$Ozone)

grid_raster_small = rasterFromXYZ(krig_grid_small)
ozone_krig_small<-raster::interpolate(grid_raster_small, krig_fit_small)
crs(ozone_krig_small) <-CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")
#plot output raster
plot(ozone_krig_small)

leaflet() %>% 
  setView(lng = -111.941004,  40.610497, zoom = 10) %>% 
  addTiles() %>% 
  addRasterImage(ozone_krig_small, project=T)%>% 
  addMiniMap()
  
  

库（光栅）
图书馆（单张）
图书馆（字段）
ozone_df如果您更新软件包（或者至少是光栅、sp和rgdal），这个错误应该会消失——使用R>=4也很好

这里有一个版本，有一些修正，对我来说是有效的

library(raster)
library(leaflet)
library(fields)
ozone_df<-data.frame("Longitude"=runif(200, -112.245075*10000,-111.455581*10000)/10000)
ozone_df$Latitude<-runif(200, 40.063614*10000,40.827281*10000)/10000
ozone_df$Ozone<-runif(200, 0,115)

#create grid tick marks
small_grid_x = seq(-111.455581,-112.245075,length.out=500)
small_grid_y = seq(40.063614,40.827281,length.out=500)

#create grid nodes
krig_grid_small <- expand.grid(small_grid_x,small_grid_y)
coordinates(krig_grid_small) <- ~ Var1 + Var2

#create kriging fit and apply interpolation to grid
krig_fit_small<-fields::Krig(ozone_df[1:2],ozone_df$Ozone)

grid_raster_small = rasterFromXYZ(krig_grid_small)
ozone_krig_small<-raster::interpolate(grid_raster_small, krig_fit_small)
crs(ozone_krig_small) <-CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")
#plot output raster
plot(ozone_krig_small)

leaflet() %>% 
  setView(lng = -111.941004,  40.610497, zoom = 10) %>% 
  addTiles() %>% 
  addRasterImage(ozone_krig_small, project=T)%>% 
  addMiniMap()
  
  

库（光栅）
图书馆（单张）
图书馆（字段）
臭氧层