У меня есть фрейм данных, который содержит долготу, широту и осадки во всем мире с 1980 по 2018. Мне нужно показать среднее количество осадков на карте тепла в R. Мои долготы и широты от -180 до 180и от -50 до 50 градусов, и для каждой сетки она изменяется на 0,25 градуса. Моя проблема в том, что я не знаю, как получить только осадки для среднего запада США?
Это заголовок моих данных:
long lat daily_precip_mean(mm)
1 -97.875 35.875 0.0009090909
2 -97.625 35.875 0.0000000000
3 -97.375 35.875 0.0099999994
4 -97.125 35.875 0.0090909085
5 -96.875 35.875 0.0245443763
6 -96.625 35.875 0.0353187964
Это мой код:
# drawing a heatmap
library(ggplot2)
library(maps)
library(mapdata)
states <- map_data("state")
midwest <- subset(states, region %in% c("illinois", "iowa", "indiana", "minnesota", "missouri", "michigan", "wisconsin", "ohio"))
midwestmap = ggplot(data = midwest) +
geom_polygon(aes(x = long, y = lat, group = group), fill = "white", color = "black") +
coord_fixed(1.3)
midwestmap
midwestmap + geom_tile(aes(long, lat, fill=daily_precip_mean), data=data, alpha=0.8)+
geom_polygon(data=midwest,aes(x = long, y = lat, group = group),color = "black", fill = NA)+
scale_fill_gradientn(colours = rev(rainbow(7)))
Я не могу поставить свою карту здесь, но моя карта представляет собой квадрат осадков, но я просто хочу средний запад
