首页 > 解决方案 > 大小点 R 图使用分类器的问题

问题描述

  1. 这是普通文本 描述:我尝试使用空间点的类大小制作空间地图 只需重复专着中的工作流程:使用 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 版本中准备情节,我会遇到同样的问题。所以我猜想关于将数据划分为类的代码的问题

sp plot ggplot 版本

标签: sizepoint

解决方案


推荐阅读