我省分为十二个区,每个区都分为市镇。 我必须按地区映射一些健康数据,但我的shapefile"知道"只有自治市。所以我将我的健康数据与shapefile中的数据合并,并为同一地区的所有城市映射了相同的值。 结果是:
这是代码:
ggplot(df) +
aes(long, lat, group = group, fill = COLOR) +
geom_polygon() +
scale_fill_manual(values = palette) +
geom_path(color=NA) +
coord_equal() +
theme_opts
正如您所看到的,我删除了市政边界,但我必须添加区域边界,因为相邻的地区使用相同的颜色进行映射,但是用户无法区分它们。你有什么建议吗?
这些是要绘制的数据:
district;COLOR
101;-1
102;0
103;1
104;-1
105;0
106;1
109;-1
110;0
111;1
112;-1
113;0
114;1
这是shapefile中数据的提取(最后两列是合并的结果;市政ID是COD_ISTAT):
long;lat;order;hole;piece;id;group;COD_ISTAT;OBJECTID;COD_REG;COD_PRO;PRO_COM;NOME;SHAPE_Leng;SHAPE_Area;district;COLOR
1158442.63774317;4556057.85829892;1;FALSE;1;6445;6445.1;16072001;6608;16;72;72001;Acquaviva delle Fonti;108308.860715;132026431.639;105;0
1158509.6371677;4556058.35816913;2;FALSE;1;6445;6445.1;16072001;6608;16;72;72001;Acquaviva delle Fonti;108308.860715;132026431.639;105;0
1158571.63687283;4556185.35695679;3;FALSE;1;6445;6445.1;16072001;6608;16;72;72001;Acquaviva delle Fonti;108308.860715;132026431.639;105;0
1158659.63610077;4556177.356861;4;FALSE;1;6445;6445.1;16072001;6608;16;72;72001;Acquaviva delle Fonti;108308.860715;132026431.639;105;0
1158679.38614406;4556290.85584435;5;FALSE;1;6445;6445.1;16072001;6608;16;72;72001;Acquaviva delle Fonti;108308.860715;132026431.639;105;0
1158672.88632394;4556356.85528688;6;FALSE;1;6445;6445.1;16072001;6608;16;72;72001;Acquaviva delle Fonti;108308.860715;132026431.639;105;0
1158575.13715917;4556353.85549587;7;FALSE;1;6445;6445.1;16072001;6608;16;72;72001;Acquaviva delle Fonti;108308.860715;132026431.639;105;0
1158567.38710377;4556288.85607141;8;FALSE;1;6445;6445.1;16072001;6608;16;72;72001;Acquaviva delle Fonti;108308.860715;132026431.639;105;0
1158523.63747073;4556283.8561965;9;FALSE;1;6445;6445.1;16072001;6608;16;72;72001;Acquaviva delle Fonti;108308.860715;132026431.639;105;0
1158442.63774317;4556057.85829892;10;FALSE;1;6445;6445.1;16072001;6608;16;72;72001;Acquaviva delle Fonti;108308.860715;132026431.639;105;0
提前致谢。