我觉得有人问过许多类似的问题,但是在相关问题中我还没有找到解决问题的方法。
我正在尝试使用spplot和col.region命令更改地图上的颜色。
我尝试使用RColorBrewer软件包设置调色板。我尝试对n使用各种值:5、6、7(我使用5个分位数进行切割),甚至使用1来查看它是否会导致任何可观察到的变化,但无济于事。我也多次更改了颜色选择,但也没有带来任何变化。 我也尝试过设置colorRampPalette函数的颜色函数,但是也没有成功。
我正在使用的数据可以在这里下载(很抱歉,我无法弄清楚如何上传带有问题的数据):http://www.data.drees.sante.gouv.fr/ReportFolders/reportFolders.aspx,方法是在搜索栏中输入“ APL”并取得第一个excel文件(2015);道歉也无法将其正确链接)
我的代码如下:
library(maptools)
library(raster)
library(RColorBrewer)
library(readxl)
library(classInt)
apl15 <- read_excel("DATA USED")
france_com<- getData(name = "GADM",
country="FRA",
level=5)
paca <- subset(fr_com,NAME_1=="Provence-Alpes-Côte d'Azur")
apl15 <- data.frame(apl15)
names(apl15) <- c("code_com","name_com","apl","apl_young","pop_com","dept","apl*pop")
id_paca <- match(paca$NAME_5,apl15$name_com)
key_paca<- apl15[id_paca,"apl"]
paca$apl <- key_paca
apl.qt <- classIntervals(paca$apl,
n=5,
style="quantile")
br <- apl.qt$brks
offs <- 0.0000001
br[1] <- br[1] - offs
br[length(br)] <- br[length(br)] + offs
paca$apl_q <- cut(paca$apl,br)
display.brewer.all(type="seq")
my.palette <- brewer.pal(6, "OrRd")
spplot(paca,"apl_q",
col.region=my.palette,
main="APL en PACA", cex.=.7)
我希望这张地图具有与我选择的调色板相对应的颜色。
我的数据映射到一个区域: 。
在此先感谢您的帮助!