我正在使用shinydashboard在R leaflet地图上绘制一个大的经纬度NetCDF raster。当我点击地图时，会弹出一个窗口，显示所点击的栅格点的行、列、经纬度位置和值。（请参见下面的可重现代码） 问题是，如果栅格足够大，我会遇到栅格偏移的问题。例如，在这里我点击了一个应该有...

我想绘制一个带有两个轴的图例。具体来说，我已经组合了两个已分类的空间对象，第一个显示事件的强度，第二个显示该位置事件发生的概率。我想创建一个图例，显示组合栅格的像素在每个类别中的位置。我想创建的图例将类似于这样：带有两个轴的图例。 普通的分类数据图例如下所示：原始图例。 这是一个可...

我正在使用以下行中显示的栅格函数。 我的最后一行会产生一些输出。 输出中有一行说dimensions : 240, 320, 76800 (nrow, ncol, ncell)。 我想重新打印该图像，但只想显示前200行和前300列。 我该怎么做？以下是倒数第二行绘制整个图像。f <...

我刚刚注意到 terra::cellSize() 产生的单元格面积估计值与 raster::area() 不匹配。 首先，为什么这两种方法不能提供相同的答案？其次，哪个估算值最准确？请参见以下示例。 library(raster) #> Loading required packag...

这是我使用rasterImage创建的渐变色图例： colfunc <- colorRampPalette(c("red", "blue")) legend_image <- as.raster(matrix(colfunc(20), ncol=1)) plot.new() ra...

我有50个geotiff文件在同一个文件夹中。 它们都代表了世界不同地区的高程数据。 我想合并某些geotiff文件，我发现R中的mosaic可能会帮助我们。 我已经将这些geotiff移动到了同一个文件夹中，并编写了下面的R脚本： setwd() a<-list.files(patt...

在 R 中，我们可以使用函数 rasterToCountour 将一个 raster 转换成一个 SpatialLinesDataFrame。library(raster) f <- system.file("external/test.grd", package="raster") r...

我想用R编写一个小程序，用于交互式可视化和修改一些光栅数据集，这些数据集被视为彩色图像。用户应该打开一个文件（从终端也可以），绘制它，选择要用鼠标单击编辑的点，并插入新值。 到目前为止，我已经轻松实现了这一点。我使用“raster”包中的“plot()”函数来可视化绘图，然后使用“click...

问题陈述： 使用ggquiver::geom_quiver（），我们可以绘制向量场，前提是我们知道x，y，xend和yend的值。 我如何为任意高程的RasterLayer计算这些参数？ 我如何确保这些箭头的大小显示该向量的坡度，以便箭头在该位置的梯度成比例地不同长度显示（例如下面的第一张图...

我正在尝试使用this文章中的方法，利用R软件对气候条件进行降尺度处理。我已经接近成功，但还缺少一些步骤。 所需软件包和数据 为了进行这个例子，我上传了一些数据到archive.org网站，使用以下代码加载所需的软件包和数据： library(raster) library(rgdal)...