Я использовал Растровый пакет, чтобы загрузить карту Испании по провинциям (уровень 2), и я хотел бы заполнить их цветами в соответствии с theincome per capita. Здесь есть файл с доходом на акцию и номер, присвоенный каждой провинции. Обратите внимание, что ID_2 - это номер, назначенный пакетом Raster, а переменная PROV является официальным номером, присвоенным испанским правительством.Цветные полигоны карты в R?
library(raster)
esp<-getData('GADM', country="ESP", level=2)
espPols <- unionSpatialPolygons(esp, esp$ID_2)
renta <- read.table("renta.csv",sep = ";", header=TRUE)
espMapRenta <- SpatialPolygonsDataFrame(espPols, renta)
plot(espMapRenta)
Первая проблема, я сталкиваюсь, что есть некоторые провинции повторяются в пакете, а второй является то, что я не знаю, как заполнять каждую область в цвете градиента по уровню дохода.
Спасибо вам большое за помощь! PS. Ссылка на данные о доходе на душу населения приведена здесь: https://www.dropbox.com/s/si6zpv7p2nap9zg/renta.csv?dl=0
да! definitiley да! Тем не менее, я только что понял, что есть ошибка, когда я объединяю оба набора данных. Альмерия, желтая провинция на юге, имеет доход на душу населения в 20828, поэтому я думаю, что данные не сливаются должным образом. – adrian1121
если ответ сделал трюк, пожалуйста, не забудьте принять ответ :-) – MLavoie