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
..
我有以下问题:我需要使用R包landscapemetrics中的同一函数处理多个栅格文件。基本上,我的栅格文件是国家地图的一部分,形状和大小都相同(即象限)。我想出了一个文件的代码,但我必须对600多个栅格做同样的事情。因此,手动操作是非常不合理的。我的代码中的步骤如下: # 1. I load "raster" and "landscapemetrics" packages: libra
..
我正在尝试使用栅格包中的简单计算操作(将温度从开尔文转换为摄氏度)来处理R中的数千个GeoTIFF图像。 计算正在进行,但我需要原始输入文件名作为输出,以便在以后分配图像。 使用我当前的代码时,只会采用列表中的第一个文件名,并且输出文件会不断被覆盖。 包含要使用‘calc’的路径的完整文件列表: library(raster) files
..
我有一个名为“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)
..
我想把澳大利亚的气候带地图标出来。相关文件(ASCII文件)已从以下位置下载:http://www.bom.gov.au/web01/ncc/www/climatology/climate-classification/kpngrp.zip 原始地图如下所示。 但当我使用ASCII文件重新绘制这张地图时,我无法完成。使用的代码为: library(ggplot2);library
..
使用SF聚合栅格中的值 我需要的是为每个栅格聚合一些度量值。 假设我们有一些数据坐标和值,并且我想创建一个热图。 首先,我使用Simple Feature框架创建了一个网格和栅格。 现在我需要获取df中的每个坐标,并检查它们是否在一个栅格中。然后为每个栅格计算任何其他聚合函数的平均值。 # Packages -------------------------------------
..
我处理netCDF数据已多年。NetCDF用于空气污染物数据,纬度和经度作为单独的变量提供,而不是原始网格的一部分。 链接到日期:Sample Netcdf 这些netCDF文件提供2级二氧化氮数据,可从NASA EarthData门户网站下载。卫星是Sentinel-5P,仪器是Tropomi。 因此,在处理这些数据时,您必须为NO2、纬度和经度创建变量。我正在尝试创建栅格图层,然
..
我在R中有一个栅格类型对象,它是一个经度/经度栅格,每个单元格都有一个值,即该单元格的深度或高程(从玩具示例中的marmapgetNOAA.bathy(),bathymap下载)。 然后我就有了一个具有最晚、最长坐标的点列表(points)。对于每个点,我希望找到距离值为负(深度)的点最近的栅格单元。如果有多个单元格,那么我想随机选取其中一个。一旦我有了/一个最近的负网格单元格,我想要查找并保存
..
我正尝试在R中打开一个ASCII文件。该文件是从热像仪软件导出的。它有农作物的温度值。我正在使用以下代码: library(raster) r = raster("AA092800_1.asc") plot(r) 但每次我都跟错了 “不能识别为支持的文件格式 .rasterObjectFromFile中出错(x,band=band,objecttype=“RasterLayer”,:
..
有人知道为什么rasterize在多边形上使用时会返回以下错误: Error in rv[[ii]] : subscript out of bounds 使用其他常见的R函数(如mean)也会产生相同的错误,但rasterize选项(如'last')可以很好地工作。在以前的raster版本(目前使用的是RASTER_2.3-0)上,我没有遇到过这个问题。我附上了我从?rasterize
..
我正在尝试将一些多波段卫星图像合并到R中(不幸的是,由于版权问题,我不能分享这些图像)。它们具有相同的大小和分辨率,但它们之间略有重叠。当我单独导入每个图像并合并它们时,没有问题: library(raster) raster1
..
我需要在栅格图上绘制一个椭圆。我试着用简单的plot(r1)(r1是栅格层)绘制栅格,然后用add=T绘制椭圆,但不起作用。然后我尝试axes=F绘制栅格,并再次尝试add=T绘制椭圆。它仍然不起作用。 所以我尝试将椭圆数据转换为数据帧,并尝试添加到栅格图。 #Creating a raster r
..
我下载了一个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文件
..
我需要在LatLong中将地图投影到方位角等距投影。 map_proj
..
要提取ENVI.hdr文件的特定信息,我希望使用caTools::Read.ENVI()将其读入R。遗憾的是,R只是抛出错误消息,我不知道如何解决此问题。 到目前为止,我尝试的内容如下: # install.packages("rgdal") # install.packages("raster") # install.packages("caTools") library("rgda
..
我有一个二进制矩阵列表。在每个矩阵中,我希望检测被连接的黑色像素(1)的环(链)包围的白色像素(0)区域。 例如,在下面的矩阵中,有两个白色像素(零)区域,它们都完全被连接的1组成的“链”包围:2x2和3x2组0。 m # [,1] [,2] [,3] [,4] [,5] [,6] [,7] # [1,] 1 1 1 1 0 0
..
所以我将这两个栅格合并为一个包含高程值的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]
..
我有一堆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
..
我尝试在栅格堆栈的层名称中添加空格 names(predstack)[[1]] names(predstack)[[1]] [1] "MSR670.max" 我知道这很愚蠢,但我计算了一个模型,花了7天的时间计算。不幸的是,有一个模型变量的名称中包含空格。 现在我想对堆栈进行预测,但这不起作用,因为出现以下错误: > Prediction2mod
..
我有0.5度分辨率的栅格数据。我想. 使用最近邻插值法将这些值内插到0.125度和1度 如何在R中实现这一点?以下是一些示例数据: library(raster) ras
..