不使用渐变在R中制作合唱曲
我对R很陌生,我正试图用adm1(或州)分区绘制阿根廷蚂蚁物种的分布图 我已经从GADM网站下载了数据,我创建了一个csv文件,其中包含说明每个adm1中是否存在该物种的信息。 即使我没有梯度,我还能做一个合唱团吗?如果没有,我可以使用哪些其他类型的地图 我看过几个网站,包括,和Choropleth地图挑战赛,它们确实很有帮助,但它们都有数字数据,我使用的是当前(1)或不存在(0)列。我尝试过的不同软件包有不使用渐变在R中制作合唱曲,r,map,R,Map,我对R很陌生,我正试图用adm1(或州)分区绘制阿根廷蚂蚁物种的分布图 我已经从GADM网站下载了数据,我创建了一个csv文件,其中包含说明每个adm1中是否存在该物种的信息。 即使我没有梯度,我还能做一个合唱团吗?如果没有,我可以使用哪些其他类型的地图 我看过几个网站,包括,和Choropleth地图挑战赛,它们确实很有帮助,但它们都有数字数据,我使用的是当前(1)或不存在(0)列。我尝试过的不同软件包有sp(使用RColorBrewer)、ggplot2、rgeos、和maptools 以下
sp
(使用RColorBrewer
)、ggplot2
、rgeos
、和maptools
以下是我目前掌握的代码:
library(sp)
library(RColorBrewer)
write.csv(atr, "atr_data.csv")
atr_data<-read.csv("atr_data.csv", header=TRUE)
spcode country_code adm1_code newcol
1 atr VEN VE.AR 0
2 atr PRY PY.CE 0
3 atr PAN PA.CL 0
4 atr PAN PA.CL 0
5 atr PAN PA.PN 0
6 atr PAN PA.PN 0
库(sp)
图书馆(RColorBrewer)
write.csv(atr,“atr_data.csv”)
atr_数据在您提供的示例中出现了几个小问题
首先,可以使用gadm@data$COLU no
而不是gadm$COLU no
。填写完状态/缺勤表后,通过直接访问插槽或使用maptools
包中的spCbind
可以将状态/缺勤数据添加到空间多边形数据框中
第二,如果您的col_no
系数中只有2个级别,则必须将mypalete
子集为2种颜色,因为Brewer调色板只能使用至少3个级别
library(sp)
library(RColorBrewer)
con <- url("http://gadm.org/data/rda/ARG_adm1.RData")
print(load(con))
close(con)
# Randomly assigning presence/absence data for display purposes only
gadm@data$col_no <- as.factor(rbinom(n = 24, size = 1, prob = 0.5))
myPalette <- brewer.pal(3, "Purples")
# col.regions is limited to 2 colors below with the middle color dropped.
spplot(gadm, zcol = "col_no", colorkey = TRUE, col.regions = myPalette[-2],
main="Distribution of Atratus in Argentina")
库(sp)
图书馆(RColorBrewer)
con在您提供的示例中,出现了几个小问题
首先,可以使用gadm@data$COLU no
而不是gadm$COLU no
。填写完状态/缺勤表后,通过直接访问插槽或使用maptools
包中的spCbind
可以将状态/缺勤数据添加到空间多边形数据框中
第二,如果您的col_no
系数中只有2个级别,则必须将mypalete
子集为2种颜色,因为Brewer调色板只能使用至少3个级别
library(sp)
library(RColorBrewer)
con <- url("http://gadm.org/data/rda/ARG_adm1.RData")
print(load(con))
close(con)
# Randomly assigning presence/absence data for display purposes only
gadm@data$col_no <- as.factor(rbinom(n = 24, size = 1, prob = 0.5))
myPalette <- brewer.pal(3, "Purples")
# col.regions is limited to 2 colors below with the middle color dropped.
spplot(gadm, zcol = "col_no", colorkey = TRUE, col.regions = myPalette[-2],
main="Distribution of Atratus in Argentina")
库(sp)
图书馆(RColorBrewer)
骗局