如何使R的“raster”软件包在GeoTIFF中区分正负旋转矩阵？
在 R 中，似乎 raster 包不区分 GeoTIFF 的正旋转和负旋转。我有一种感觉，那就是 R 忽略了旋转矩阵中的负号。我不太懂得深入挖掘 raster 源代码来验证，但我确实创建了一个可重现的示例来演示这个问题:读取 R 徽标并保存为 GeoTIFF。library(raster) b...
使用R raster进行交互式绘图：鼠标悬停显示值
我想用R编写一个小程序，用于交互式可视化和修改一些光栅数据集，这些数据集被视为彩色图像。用户应该打开一个文件（从终端也可以），绘制它，选择要用鼠标单击编辑的点，并插入新值。 到目前为止，我已经轻松实现了这一点。我使用“raster”包中的“plot()”函数来可视化绘图，然后使用“click...
在raster包中处理sf对象
之前我使用raster::crop和raster::mask处理Spatial*类的shapefile文件，使用rgal::readOGR读取。 现在我要“升级”我的脚本，使用sf来读取和操作多边形。 raster::crop raster::crop函数的第二个参数需要一个'exte...
R - 将SpatialLines转换为栅格
在 R 中，我们可以使用函数 rasterToCountour 将一个 raster 转换成一个 SpatialLinesDataFrame。library(raster) f <- system.file("external/test.grd", package="raster") r...
将SpatVector转换为sf或sp
有没有一种简单的方法可以将多边形的SpatVector类对象（来自terra库）转换为简单要素或SpatialPolygonsDataFrames？
如何在R中使用外推法来外推栅格数据
我正在尝试使用this文章中的方法，利用R软件对气候条件进行降尺度处理。我已经接近成功，但还缺少一些步骤。 所需软件包和数据 为了进行这个例子，我上传了一些数据到archive.org网站，使用以下代码加载所需的软件包和数据： library(raster) library(rgdal)...
图像光栅 R 包 - 无背景、边框和图例的光栅打印
我正在使用以下行中显示的栅格函数。 我的最后一行会产生一些输出。 输出中有一行说dimensions : 240, 320, 76800 (nrow, ncol, ncell)。 我想重新打印该图像，但只想显示前200行和前300列。 我该怎么做？以下是倒数第二行绘制整个图像。f <...
在R中，在leaflet地图上识别光栅上的点击位置
我正在使用shinydashboard在R leaflet地图上绘制一个大的经纬度NetCDF raster。当我点击地图时，会弹出一个窗口，显示所点击的栅格点的行、列、经纬度位置和值。（请参见下面的可重现代码） 问题是，如果栅格足够大，我会遇到栅格偏移的问题。例如，在这里我点击了一个应该有...
计算和绘制任意栅格层的矢量场
问题陈述： 使用ggquiver::geom_quiver（），我们可以绘制向量场，前提是我们知道x，y，xend和yend的值。 我如何为任意高程的RasterLayer计算这些参数？ 我如何确保这些箭头的大小显示该向量的坡度，以便箭头在该位置的梯度成比例地不同长度显示（例如下面的第一张图...
在R中对栅格堆栈进行n层求和
我正在使用以.nc文件组织的气候数据进行日常观测。我使用raster软件包的stack命令来读取这些数据。每个文件（对应于一年）都是一个RasterStack元素，具有以下特征：class : RasterStack dimensions : 360, 720, 259200,...