矩形像素阵列。栅格定义占据平面的特定矩形区域的像素的值,不一定包括(0,0)。
矢量化函数,调用数据帧以在 raster::overlay() 中使用
我有一个函数,可以根据数据帧中指定的给定条件更改 RasterLayer。我想使用overlay()函数来调用这个函数来改变RasterStack,但是遇到了麻烦
有聪明的人知道如何在栅格堆栈而不是矩阵或 data.frame 上编写 quantForestError() 预测吗?例如... 图书馆(护林员) 图书馆(tidyverse) 库(光栅)
我想在栅格中插入缺失值: <- rast(xmin=0, xmax=10, ymin=0, ymax=10, nrows=50, ncols=50) e = rast(xmin=2, xmax=4, ymin=2, ymax=4, nrows=10, ncols=10) values(r) <...
GDAL 关于像素行/列到地图坐标之间的地理变换的用户指南说使用以下参数数组: GT[0] 左上角像素的 x 坐标。 GT[1] w-e...
如何在 R terra 中裁剪、投影和合并来自反子午线两侧的栅格数据?
我正在从全球latlon栅格(EPSG:4326)中裁剪数据,然后希望将该数据投影到合适的局部坐标参考系中。我正在裁剪的区域跨越了反子午线,所以我......
我想将名为“havza”的shapefile重叠到名为“gio_ara”的栅格数据上,但shapefile文件出现在绘图中栅格的后面,我该如何解决这个问题? 黑线位于...
我正在尝试使用光栅包中的堆栈功能将多个光栅文件(.tif)放入一个堆栈中。当我这样做时,我遇到了两个错误,即: > covStack <- stack(paste0(getwd(),&q...
terra::predict() 直接应用于 spatraster
我在R中创建了一个随机森林模型,并希望使用terra中的predict()函数将其应用到spatraster对象。我可以将栅格转换为数据帧,而不是直接预测...
我刚开始在 R 中使用 terra 包。我正在尝试使用包的绘图函数来绘制 SpatRaster。但是,我想更改多个布局特性但没有找到新的...
如何在没有 RGDAL 的情况下在 R 中打开 .tif 文件(因为现在未使用)?
我有一个 .tif 文件,但无法在 R 中打开它。 我无法安装 rgdal,因为如果我是正确的,它就不再使用了。 我没有找到任何其他有效的方法来打开它。 我尝试了以下脚本: 我...
我有两个光栅堆栈,每个都有大约。 80000层。这些堆栈具有相同的范围、相同的 CRS、相同的分辨率和相同的层数。我想添加堆栈元素的层...
我正在研究 R 中的物种分布模型,并试图在此过程的早期解决代码中的一些问题。我使用 WorldClim 的气候数据栅格创建了一个掩模。然后我创建
图书馆(terra) 数据 <- structure(list(lon = c(73.0754166666667, 73.0754166666667, 73.0754166666667, 73.0754166666667, 73.07625, 73.07625, 73.0770833333333, ...
我正在按纬度、经度和图层聚合一个大型 RasterBrick。我原来的 RasterBrick 对象有一个设定的 Z 维度(即为每个图层分配了日期),但是当我聚合时我丢失了这些日期...
我在很长一段时间内有一大堆海面温度栅格层(每小时一个),但是我想找到这些层在六个小时周期 00:00 - 0600、0600 - 12 的平均值。 ..
如何在rasterio中保存(作为PNG)带有假色的单波段热图像geotiff栅格
我是新手,请原谅我的无知。 我有一个带有一个波段的 GEOTIFF 栅格(这是芝加哥的热图像)。 这是元数据: 司机: GTiff 数据类型:uint8 无数据:无 宽度:7921 身高:
我想在南美降水栅格上添加世界地图形状文件。我使用了下面的代码 # 调用变量 预 <- raster("C:/Users/DELL/OneDrive/Desktop/test1/wc2.
我想并行执行以下向量*光栅乘法,因为我的向量很长: 图书馆(土地) # 创建一个栅格示例 设置.种子(123) 光栅_示例 <- rast(nro...
如何选择在使用Python Rioxarray重新投影时不更改图像大小?
我在坐标系 EPSG:4326 中有一个大小为 18346x10218 的 tif 格式的图像,我想将其转换为新的坐标系,我使用 Python 中的 roxarray 使用以下方法执行此操作。 我...
给定这两个栅格: 库(光栅) rr1 <- raster(ncol=10, nrow=10) rr1[] <- runif(ncell(rr1)) rr1[1]=2 rr1[10]=2 rr1[50]=2 ...