[英]R spplot label overlapping
我正在处理的数据对于每个级别(东南5个州的每日降雨量)具有很多相同的值,例如:
values(r1) <- c(rep(0.1,25575),rep(0.2,25575),rep(0.3,25575),rep(0.3,25575))
因此,在输出图像中,映射对于每个级别都充满了相同的标签/值,是否存在这样的方式,即每个级别仅将标签显示一次或两次,例如:标签0.1对于级别1仅显示一次,而不是25575次?
这是一个测试代码:
library(sp)
library(rgdal)
library(raster)
library(maptools)
setwd("C:/Users/rosaR/Data/PRISM")
r1 <- raster(xmn=-88.77083, xmx=-75.02083, ymn=24.0625,ymx=36.97917, nrow=310, ncol=330,crs='+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0')
#values(r1) <- c(rep(0.1,25575),rep(0.2,25575),rep(0.3,25575),rep(0.4,25575))
values(r1) <- c(rep(seq(from=0.1,to=3,by=0.1),3410))
newproj <- "+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs"
r1 <- projectRaster(r1, crs=newproj)
png(file=paste("test.png",sep=""), width=1350, height=700, units="px")
at=c(0,0.15,0.3,0.45,0.5,0.65,0.7,0.85,0.9,1.05,1.1,1.25,1.3,1.45,1.5,1.65,1.7,1.85,1.9,2.05,2.1,2.25,2.3,2.45,2.5,2.65,2.7,2.85,2.9,3)
spplot(r1,
#sp.layout=list(
#list('sp.polygons', spTransform(SE_sea1, CRS(newproj)),fill="white", first=F, col=NA),
#list('sp.polygons', spTransform(SE1, CRS(newproj)), alpha=1, lwd=2, first=F)),
xlim=c(613190.6,2173441),ylim=c(121319.4,1774699),col.regions=rainbow(29),colorkey=F,at=at+.Machine$double.eps,contour=T,labels=list(labels=format(at)),col=grey.colors(1),maxpixels=102300)
dev.off()
如果我将分界点设置为= c(0,0.15,0.25,0.35,0.5),地图将看起来很完美,但是由于我的原始数据集非常庞大,并且不可避免地存在与级别值相同的值,可能不超过25575 ,但要20-100倍,我希望减少到图像中显示的少于5倍。
我试图添加:1)method =“ flattest”; 2)list('sp.pointLabel',r1,label,col ='grey'); 3)list(“ sp.text”,坐标(raster_data),raster_data); 但似乎没有一个适用于我的RasterLayer,
任何建议将不胜感激,谢谢!
您可以at
值上添加.Machine$double.eps
,以便截止值仅会错过at值。 然后,您可以使用labels
参数调整显示的标签(这样您将获得“ 0.1”而不是“ 0.100000e + 02”)
spplot
的labels
参数传递给levelplot
( ?sp
说...
传递给levelplot
(网格,多边形)或xyplot
(点),但是raster
是“网格”类型的)。 阅读?levelplot
表示要查看?panel.levelplot
以获取有关labels
参数的更多详细信息。 最后, ?panel.levelplot
说, labels
可以是带有元素“标签”的列表,该元素描述放置在轮廓上的标签。
所以结果是您的代码,除了:
at
以添加.Machine$double.eps
(少量) label
,使标签为“好”数字。 赠送:
at <- c(0, 0.1, 0.15, 0.25, 0.35, 0.45)
spplot(r1, xlim=c(613190.6,2173441), ylim=c(121319.4,1774699),
col.regions=rainbow(4), colorkey=F, contour=T,
# add a small number to the colour cutoffs to prevent ties
at=at + .Machine$double.eps,
# change the labels to be "pretty". See `?format`
labels=list(labels=format(at)),
col=grey.colors(1), maxpixels=102300)
(颜色看起来很糟糕,但我认为您可以从这里修复它:P)
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.