Трудно дать количественную оценку «лучше» - но описанная вами работа, безусловно, будет проще и быстрее, если вы используете формат пакета sf
, а не пространственные полигоны из пакета sp
.Он также доступен на сайте GADM.
Пространственные объекты sf
представляют собой модифицированные фреймы данных, что означает, что вы можете применять к ним стандартные методы обработки данных - например, dplyr::mutate()
для манипулирования именами и т. Д.
Что касается объединения регионов, вы можете использовать рабочий процесс на основе sf::st_union()
, который я описал некоторое время назад в сообщении в блоге: https://www.jla -data.net / eng / disolving-polygons-in-sf-environment / - он даже использует tmap
в качестве инструмента для построения графиков, поэтому вы должны найти рабочий процесс знакомым.
И, наконец, предупреждение относительно пятилетнего вопроса о построении картс ggplot вы понимаете: подход, использующий fortify()
, к настоящему времени устарел, вы получите гораздо лучшие результаты, используя функцию geom_sf()
текущего ggplot *1019*.