size - 大小点 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 版本中准备情节,我会遇到同样的问题。所以我猜想关于将数据划分为类的代码的问题
解决方案
推荐阅读
- azure - Azure AD B2C 审核日志 - 图形 API
- elixir - 更改集未在下拉数据更新时更新
- flutter - 导航弹出后Flutter获取参考
- google-apps-script - 如何使用脚本编辑器在单元格中添加今天的日期加上 7 天
- node.js - 为什么请求超时时不抛出错误?
- awk - 通过搜索前缀提取整个单词
- r - 错误使用 gam 模型分析浮游植物丰度和环境参数(mgcv 包)
- jquery - 由于某种原因无法使用更改方法的无线电输入
- c# - MPRemoteCommandCenter 播放/暂停命令未正确更新
- python - Python - 如何从 .csv 文件读取到字典?