堆栈内存溢出
cost 83 ms
如何读取 GeoTIFF 文件并以 R 语言返回特定数据类型？ - How to read a GeoTIFF file and return in a specific datatype in R language?

我正在尝试通过raster package 读取 R 中无符号字节数据格式的大型 GeoTIFF 文件（18000×18000×40）。我希望将其读取为 R 中的raw数据类型（无符号字节），然后由Rcpp处理。 但是， raster package 的函数似乎只返回numeric数据类型（8 字 ...

Function 用于对焦方法 - Function for the focal method

我正在尝试使用 focal 运行一个简单的 function，它在每个 window 中取 5 个随机值，然后根据这些值计算方差系数 (sd/mean)。 下面的代码当前返回错误：评估错误：第一个参数无效。 提前致谢！ ...

2023-01-30 03:31:14   1   27    r / r-raster  
Sentinel-1 & 2：在不同图块的图像中创建重合点 - Sentinel-1 & 2: create coincident points in images of different tiles

我有一个 Sentinel-1 图像和 Sentinel-2 图像，我想让 Sentinel-1 中的像素与 Sentinel-2 重合。 更简单的解决方案是使用S2 <- resample(S2, S1)将两个图像放在相同的几何形状（空间范围和分辨率）中，但我担心光谱信息的完整性。 更好地重 ...

打开带有光栅的 a.nc 文件时出错 package - Error when opening a .nc file with raster package

我是 R 中光栅 package 的新手，我试图用 package 光栅打开一个 .nc 文件，但弹出了一些错误。 如果您想尝试，我使用的是来自哥白尼的 2018 年和 2019 年每月海水盐度数据集（网格是魁北克圣劳伦斯河口和加斯佩西海岸）。 以前打开过类似的数据文件，但从来没有出现过这个错误， ...

使用边界定义的焦点填充 NA 栅格单元 - fill NA raster cells using focal defined by boundary

赏金将在 2 天后到期。 此问题的答案有资格获得+50声望赏金。 89_Simple想让更多人关注这个问题。 我有一个栅格和一个 shapefile。 栅格包含 NA，我正在使用focal function 填充 NA 但是，还有一些 NA 仍然存在。 所以我这样做： 填满所有 NA 后，我再次 ...

与多个人创建空间线 Object - Creating a Spatial Line Object with multiple individuals

我正在尝试为我的数据集中的多个人创建一条空间线，然后跨线提取栅格值。 我已经尝试了一些事情，但不幸的是我正在使用的代码将所有线路分配给一个人。 WD是其中有37个人的数据集。 x 和 y 是坐标。 如果有人知道如何执行此操作，我将不胜感激。 谢谢！ 我希望能够沿线提取栅格值。 ...

为 R 中的每个 GPS 位置从 RasterBrick 的多层中提取值 - Extract values from multiple layers of a RasterBrick for each GPS location in R

我正在尝试计算 GPS 位置的 USDA CropScape 作物价值，该位置的日历年、上一个日历年和下一个日历年（例如，如果 GPS 位置的日期是 10/1/2016 ，我想提取该位置 2015、2016 和 2017 年的作物价值）。 这是我用来准备 GPS 和 CropScape 数据的代码： ...

2023-01-03 22:37:54   2   102    r / r-raster  
如何使用单层 shapefile 屏蔽 multi-layer.netCDF 或栅格？ - How do I mask a multi-layer netCDF or raster with a single-layer shapefile?

我目前正在处理 .netCDF 格式的每日降水数据。 数据分辨率为 4 公里，覆盖美国。 但是，我想为特定地理区域（大约一个县的大小）使用更高分辨率的 shapefile 来屏蔽/剪辑数据。 最终，我希望 output 成为更小区域的每日降水数据，无论是高分辨率还是原始 4km 分辨率。 我尝试了 ...

尝试使用 raster::projectRaster 重新分配 crs 时出错：rbind 中的错误：矩阵的列数必须匹配 - Error when attempting to reassign crs with raster::projectRaster: error in rbind:number of columns of matrices must match

我有以下名为 pred_th 的 SpatRaster 对象。 我正在尝试将 crs 重新分配给 ESPG 2970。当我尝试使用raster函数projectRaster执行以下操作时： 我收到此错误消息： 我很困惑，如果有人可以提供帮助，我将不胜感激。 ...

使用 R 将栅格 DEM 分辨率从 30m 更改为 500m 时，最佳使用方法是什么？ - What is the best way to use when changing raster DEM resolution from 30m to 500m using R?

我正在尝试使用 R 以 500 米的分辨率使用 SRTM-DEM 计算高程、坡度和坡向。SRTM-DEM 坐标系在 UTM 中，分辨率为 30m。 有人可以告诉我如何将分辨率更改为 500m 并以相同的范围计算高程、坡度和坡向。 我想知道我是否可以使用聚合方法获得准确的 500m 分辨率。 ...

R ggplot 绘图 map 圆形栅格 - 如何删除投影区域外的数据？ - R ggplot plotting map raster with rounded shape - How to remove data outside projected area?

我正在尝试 plot 在协调系统中投影的栅格，该系统像大多数非 WGS84 投影一样遵循地球曲率。 问题是地球环绕数据的地方不应该绘制在地球之外。 我意识到 ggplot 不能做一个圆形/椭圆形 plot 但我如何屏蔽或自动删除地球以外的数据？ 我必须 plot 超过 100 张地图，我不能手动执 ...

R terra：跨层滚动求和 - R terra: rolling sum across layers

拥有大量时间序列栅格，我想创建另一个 rasterstack，它代表图层中每个像素的 5 天滚动总和。 鉴于光栅堆栈的大小， raster::calc会为此任务无休止地运行，因此我想知道是否可以使用terra::lapp更快地实现此操作。 然而： 返回错误： 如果有人可以建议错误出在哪里，或者建议另 ...

使用 R 中的 sp、terra、raster 包将阿尔伯斯等面积转换为纬度/经度 - Using sp, terra, raster packages in R to convert Albers Equal Area to Lat/Long

我有一个关于地球上某个区域的地面温度的.tif 文件，我想找出该区域的坐标。 这是我正在处理的文件的链接： 2017001D.tif 在 R 中使用raster R ，我能够将 .tif 文件加载为RasterLayer class。 我可以将坐标信息提取为 我还使用terra package 将其 ...

如何计算 terra package 中的聚合事实参数？ - How to calculate aggregate fact parameter in terra package?

terra package 具有聚合 function，允许创建具有较低分辨率（较大单元格）但需要事实参数的新 SpatRaster。 转换很多栅格时，每次都需要计算fact，有没有办法根据另一个栅格的目标分辨率传递fact参数？ 其他函数将现有栅格作为输入，例如 function(r1,r2) ...


排序:
质量:
 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM