2016-03-13 3 views
3

Я использовал Растровый пакет, чтобы загрузить карту Испании по провинциям (уровень 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

ответ

1

Вы ищете что-то в этом роде?

library(sp) 
spplot(espMapRenta, zcol="Renta.per.Capita") 

enter image description here

+0

да! definitiley да! Тем не менее, я только что понял, что есть ошибка, когда я объединяю оба набора данных. Альмерия, желтая провинция на юге, имеет доход на душу населения в 20828, поэтому я думаю, что данные не сливаются должным образом. – adrian1121

+0

если ответ сделал трюк, пожалуйста, не забудьте принять ответ :-) – MLavoie

Смежные вопросы