R光栅情节甚至传说

Question

我尝试使用NetCDF文件使用raster和ncdf4包在R中绘制栅格地图。

我的数据范围太大，但几乎我的数据分布在0-2000之间。您可以查看我的数据的直方图：

所以我想绘制一个间隔为200的图，就像seq(0, 2000, 200)一样。 但是当我使用这些间隔时，大于2000的值在绘制图时被R识别为NA。这些大值在绘图时打印为透明。

使用均匀间隔的图：

我试图给情节留出很大的间隔，比如breaks = c(seq(0, 1500, 100), 40000)，但情节的传说看起来很难看。

总之，我想要一个偶数间隔和图例的情节。像开放范围，而不是封闭范围。

下图是我通过其他软件实现的理想情节。如何使用R？

解决我的问题

我的代码在这里：

    library(rgdal)
    library(raster)
    library(ncdf4)
    library(rasterVis)
    library(sp)

    ncname <- "output.nc"
    ncdata <- nc_open(ncname)
    bb <- raster(ncname)

    hist(bb)
    hist(bb, xlim = c(0, 2000), breaks = seq(0, 40000, 500))

    plot(bb)
    plot(bb, xlim = c(25, 53), asp = 1.5,  breaks = c(seq(0, 1500, 100), 40000), col = topo.colors(20)) 

    plot(bb, legend.only = FALSE, zlim = c(0, 2000), col = topo.colors(20), asp = 1.5,
         legend.width = 1, legend.shrink = 1.0,
         axis.args = list(at = seq(0, 2000, 100),
         labels = seq(0, 2000, 100)))