如何解决使用分类器用于大小点 R 图的问题
- 这是普通文本 描述:我尝试使用空间点的类大小制作空间地图 我只是从专着中重复工作流程:显示时间序列、空间和时空 数据与 R / Oscar Perpiñán Lamigueiro。见章节 8.4 优化分类以改进 歧视。脚本开始:
library(sp)
library(rgdal)
Issano_Gap_Rocks <- read_excel("Issano Gap Rocks Original File_for_plot.xls")
xy <- Issano_Gap_Rocks[,c('EASTING','NORTHING')]
Issano_spdf_Rocks <- SpatialPointsDataFrame(coords = xy,data = Issano_Gap_Rocks,proj4string = CRS("+init=epsg:24821"))
library(classInt)
airPal_Rock <- colorRampPalette(c("green","yellow","red"))(5)
op <- options(digits=4)
tab <- print(intervals)
options(op)
nClasses <- 5
intervals <- classIntervals(Issano_spdf_Rocks$Sc_ppm_INAA,n = nClasses,style = "quantile")
tab <- print(intervals)
nClasses <- length(intervals$brks) -1
dent <- c(0.64,1.14,1.65,2.79,4.32,6.22,9.65,12.95,15.11) #point size
dentAQ <- dent[seq_len(nClasses)]
idx <- findCols(intervals)
cexNO2 <- dentAQ[idx]
Issano_spdf_Rocks$classNO2 <- factor(names(tab)[idx])
2.这是普通文本
但毕竟我的点类不是按值 Sc_ppm_INAA 定期排列的。看起来
不匹配:放置在错误位置的最大间隔。我不明白是我的
错了吗?
如果在 GGPLOT 版本中准备绘图,我会遇到同样的问题。所以我猜有问题
将数据分到类的代码
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。