Map using krigging interpolation

I am trying to do kriging interpolation in the following map.
So, far I have created the map.
Please see the following code.

library(readr)
library(ggplot2)
library(maps)

baseData <- map_data('state')
id <- "10Un9up12GIi_yddG4YQpVhzuZhtYVZmA"
df_obs=read_csv(paste0("https://docs.google.com/uc?id=",id,"&export=download"),
              col_names = TRUE)

colnames(df_obs) <- c('lat', 'lon', 'obs')
limits <- round(range(df_obs$obs, na.rm = TRUE))
#dev.new()

ggplot(data=df_obs, aes(x=lon,y=lat, fill=obs)) +   
  theme_bw() +theme(panel.grid.major = element_blank(),
                    panel.grid.minor = element_blank())+
  theme(plot.title = element_text(hjust = 1, size=15)) +
  geom_tile() +
  geom_polygon(data=baseData, aes(x=long, y=lat, group=group), 
               colour="black", fill="white", alpha=0)+
  coord_fixed(ratio=1, xlim=range(df_obs$lon), ylim=range(df_obs$lat)
              , expand = FALSE)+
  geom_point(data = df_obs,
             mapping = aes(x = lon, y = lat), 
             colour = "red")

This topic was automatically closed 21 days after the last reply. New replies are no longer allowed.

If you have a query related to it or one of the replies, start a new topic and refer back with a link.