使用R语言对netCDF栅格数据集进行时间和地理子集的处理
对于以下2016年全球每日海表面温度的netcdf文件，我试图 (i) 在时间上进行子集，(ii) 在地理上进行子集，(iii) 然后对每个像素取长期平均值并创建一个基本图。 文件链接: 这里 library(raster) library(ncdf4) 设置工作目录后打开netcdf...
在raster包中处理sf对象
之前我使用raster::crop和raster::mask处理Spatial*类的shapefile文件，使用rgal::readOGR读取。 现在我要“升级”我的脚本，使用sf来读取和操作多边形。 raster::crop raster::crop函数的第二个参数需要一个'exte...
图像光栅 R 包 - 无背景、边框和图例的光栅打印
我正在使用以下行中显示的栅格函数。 我的最后一行会产生一些输出。 输出中有一行说dimensions : 240, 320, 76800 (nrow, ncol, ncell)。 我想重新打印该图像，但只想显示前200行和前300列。 我该怎么做？以下是倒数第二行绘制整个图像。f <...
如果栅格值为 NA，则搜索并提取最近的非 NA 像素。
当我将栅格的值提取到点时，我发现有几个 NA，而我不想使用extract函数中的buffer和fun参数，相反，我想提取最近的非NA像素到与NA重叠的点。 我正在使用基本的extract函数： data.extr<-extract(loc.thr, data[,11:10])
使用raster包进行空间自相关性分析
亲爱的众人 问题 我尝试使用nfc、pgirmess、SpatialPack和spdep包计算空间自相关图。然而，我在定义距离的起点和终点时遇到了困难。我只对小距离下的空间自相关性感兴趣，但是这里要使用更小的区间。此外，由于栅格数据集相当大（1.8百万像素），使用这些包但SpatialPa...
如何在R中使用外推法来外推栅格数据
我正在尝试使用this文章中的方法，利用R软件对气候条件进行降尺度处理。我已经接近成功，但还缺少一些步骤。 所需软件包和数据 为了进行这个例子，我上传了一些数据到archive.org网站，使用以下代码加载所需的软件包和数据： library(raster) library(rgdal)...
如何使R的“raster”软件包在GeoTIFF中区分正负旋转矩阵？
在 R 中，似乎 raster 包不区分 GeoTIFF 的正旋转和负旋转。我有一种感觉，那就是 R 忽略了旋转矩阵中的负号。我不太懂得深入挖掘 raster 源代码来验证，但我确实创建了一个可重现的示例来演示这个问题:读取 R 徽标并保存为 GeoTIFF。library(raster) b...
俄罗斯的“球形”地图
我使用GADM数据绘制了俄罗斯地区的地图： setwd("~/Desktop/Master thesis/") library(sp) library(RColorBrewer) library(raster) data <- getData('GADM', country='RU...
将SpatVector转换为sf或sp
有没有一种简单的方法可以将多边形的SpatVector类对象（来自terra库）转换为简单要素或SpatialPolygonsDataFrames？