Как объединить растровый фрейм данных с шейп-файлом, чтобы создать карту богатства видов? - PullRequest
Новая
       40

Как объединить растровый фрейм данных с шейп-файлом, чтобы создать карту богатства видов?

0 голосов
/ 03 июля 2019

Я хочу нанести растровый график / фрейм данных поверх шейп-файла, но продолжаю получать различные ошибки в зависимости от того, как я пишу код.

Я пытался использовать + с готовыми объектами, карта местности+ critmapped, я попытался добавить коды вместе и вывести их как данные, но это не сработало:

Я пробовал следующее, а также другие вещи ...

Спасибо за любую помощь / направление. 

critmapped<-ggplot(df, aes(x, y, fill = layer)) + 
geom_raster() +
scale_fill_viridis_c(na.value = "white") +
labs(fill = "Count") + 
theme_minimal() + ggplot() + geom_path(data = land_df, aes(x = long, y = lat, group = group), color = 'black', fill = 'green')

#Error: Don't know how to add ggplot() to a plot

critmapped+landmap2
#Error: Don't know how to add landmap2 to a plot

#The shapefile code

  require(rgdal)
  land <- readOGR(dsn = "C:/Users/tjef631/Desktop/R Stats/Data/NE_10m_full", 
  layer = "ne_10m_land")
  land_df<-fortify(land)
  names(land_df)
  landmap<-ggplot() + geom_path(data = land_df, aes(x = long, y = lat, group 
  = group),
            color = 'black', fill = 'green')
landmap


#the raster/dataframe code

  critmapped<-ggplot(df, aes(x, y, fill = layer)) + 
  geom_raster() +
  scale_fill_viridis_c(na.value = "white") +
  labs(fill = "Count") + 
  theme_minimal()

  critmapped
...