raster相关内容

搜索&terra";最新版本时,R包&raster";无法上载

我在RStudio2021.09.2中使用的是Windows 10和R 4.1.2。我正在处理空间数据,包括矢量和栅格,但三天前命令库(栅格)开始给我这个警告: 错误:加载名称空间(i,c(lib.loc,.libPath()),versionCheck=VI[[i]])中‘raster’的包或命名空间加载失败: 已加载命名空间‘terra’1.4-22,但需要>;=1.5.12 ..
发布时间:2022-07-20 10:41:15 其他开发

使用R处理多个栅格文件时保留原始文件名

我有以下问题:我需要使用R包landscapemetrics中的同一函数处理多个栅格文件。基本上,我的栅格文件是国家地图的一部分,形状和大小都相同(即象限)。我想出了一个文件的代码,但我必须对600多个栅格做同样的事情。因此,手动操作是非常不合理的。我的代码中的步骤如下: # 1. I load "raster" and "landscapemetrics" packages: libra ..
发布时间:2022-07-20 10:37:13 其他开发

R中作为输出文件名的多个输入文件名

我正在尝试使用栅格包中的简单计算操作(将温度从开尔文转换为摄氏度)来处理R中的数千个GeoTIFF图像。 计算正在进行,但我需要原始输入文件名作为输出,以便在以后分配图像。 使用我当前的代码时,只会采用列表中的第一个文件名,并且输出文件会不断被覆盖。 包含要使用‘calc’的路径的完整文件列表: library(raster) files ..
发布时间:2022-05-13 20:26:54 其他开发

在r中更改范围+分辨率后出错

我有一个名为“tiff4”的tiff文件,下面是它的信息: class : RasterLayer dimensions : 49, 89, 4361 (nrow, ncol, ncell) resolution : 0.5, 0.5 (x, y) extent : 60.64135, 105.1413, 14.99839, 39.49839 (xmin, xmax, ymin, ymax) ..
发布时间:2022-05-13 20:18:27 其他开发

使用SF聚合栅格中的值

使用SF聚合栅格中的值 我需要的是为每个栅格聚合一些度量值。 假设我们有一些数据坐标和值,并且我想创建一个热图。 首先,我使用Simple Feature框架创建了一个网格和栅格。 现在我需要获取df中的每个坐标,并检查它们是否在一个栅格中。然后为每个栅格计算任何其他聚合函数的平均值。 # Packages ------------------------------------- ..
发布时间:2022-05-13 20:04:06 其他开发

从R中的netcdf创建栅格最准确的方法是什么?

我处理netCDF数据已多年。NetCDF用于空气污染物数据,纬度和经度作为单独的变量提供,而不是原始网格的一部分。 链接到日期:Sample Netcdf 这些netCDF文件提供2级二氧化氮数据,可从NASA EarthData门户网站下载。卫星是Sentinel-5P,仪器是Tropomi。 因此,在处理这些数据时,您必须为NO2、纬度和经度创建变量。我正在尝试创建栅格图层,然 ..
发布时间:2022-05-13 19:59:28 其他开发

对于给定点,查找具有负值R的最近格网像元

我在R中有一个栅格类型对象,它是一个经度/经度栅格,每个单元格都有一个值,即该单元格的深度或高程(从玩具示例中的marmapgetNOAA.bathy(),bathymap下载)。 然后我就有了一个具有最晚、最长坐标的点列表(points)。对于每个点,我希望找到距离值为负(深度)的点最近的栅格单元。如果有多个单元格,那么我想随机选取其中一个。一旦我有了/一个最近的负网格单元格,我想要查找并保存 ..
发布时间:2022-05-13 19:55:34 其他开发

使用R打开ASCII文件

我正尝试在R中打开一个ASCII文件。该文件是从热像仪软件导出的。它有农作物的温度值。我正在使用以下代码: library(raster) r = raster("AA092800_1.asc") plot(r) 但每次我都跟错了 “不能识别为支持的文件格式 .rasterObjectFromFile中出错(x,band=band,objecttype=“RasterLayer”,: ..
发布时间:2022-05-13 19:46:18 其他开发

R:在多边形上使用光栅化时出错

有人知道为什么rasterize在多边形上使用时会返回以下错误: Error in rv[[ii]] : subscript out of bounds 使用其他常见的R函数(如mean)也会产生相同的错误,但rasterize选项(如'last')可以很好地工作。在以前的raster版本(目前使用的是RASTER_2.3-0)上,我没有遇到过这个问题。我附上了我从?rasterize ..
发布时间:2022-05-13 19:35:11 其他开发

在R中的栅格地块上添加椭圆

我需要在栅格图上绘制一个椭圆。我试着用简单的plot(r1)(r1是栅格层)绘制栅格,然后用add=T绘制椭圆,但不起作用。然后我尝试axes=F绘制栅格,并再次尝试add=T绘制椭圆。它仍然不起作用。 所以我尝试将椭圆数据转换为数据帧,并尝试添加到栅格图。 #Creating a raster r ..
发布时间:2022-05-13 19:27:59 其他开发

R在打开netcdf文件时崩溃

我下载了一个netcdf文件,并尝试在R中打开它。以下是我的代码 download.file("https://data.giss.nasa.gov/impacts/agmipcf/agmerra/AgMERRA_1980_prate.nc4",destfile = "AgMERRA_1980_prate.nc4", method="libcurl") 我要使用R 打开netcdf文件 ..
发布时间:2022-05-13 19:21:26 其他开发

R如何读取ENVI.hdr-文件?

要提取ENVI.hdr文件的特定信息,我希望使用caTools::Read.ENVI()将其读入R。遗憾的是,R只是抛出错误消息,我不知道如何解决此问题。 到目前为止,我尝试的内容如下: # install.packages("rgdal") # install.packages("raster") # install.packages("caTools") library("rgda ..
发布时间:2022-05-13 19:05:15 其他开发

标识矩阵中被一包围的零区域

我有一个二进制矩阵列表。在每个矩阵中,我希望检测被连接的黑色像素(1)的环(链)包围的白色像素(0)区域。 例如,在下面的矩阵中,有两个白色像素(零)区域,它们都完全被连接的1组成的“链”包围:2x2和3x2组0。 m # [,1] [,2] [,3] [,4] [,5] [,6] [,7] # [1,] 1 1 1 1 0 0 ..
发布时间:2022-05-13 19:01:57 其他开发

寻找一座山的顶峰

所以我将这两个栅格合并为一个包含高程值的DEM栅格: dem1 = read_stars("srtm_43_06.tif") dem2 = read_stars("srtm_44_06.tif") pol = st_read("israel_borders.shp") dem = st_mosaic(dem1, dem2) dem = dem[, 5687:6287, 2348:2948] ..
发布时间:2022-05-13 18:56:02 其他开发

计算相邻像素随时间变化的平均相关性

我有一堆4个栅格。我想要一个像素与其8个相邻像素之间的平均时间相关性。 部分数据: library(raster) r1=raster(matrix(runif(25),nrow=5)) r2=raster(matrix(runif(25),nrow=5)) r3=raster(matrix(runif(25),nrow=5)) r4=raster(matrix(runif(25 ..
发布时间:2022-05-13 18:44:25 其他开发

如何在栅格名称中插入空格?R

我尝试在栅格堆栈的层名称中添加空格 names(predstack)[[1]] names(predstack)[[1]] [1] "MSR670.max" 我知道这很愚蠢,但我计算了一个模型,花了7天的时间计算。不幸的是,有一个模型变量的名称中包含空格。 现在我想对堆栈进行预测,但这不起作用,因为出现以下错误: > Prediction2mod ..
发布时间:2022-04-22 21:51:49 其他开发