如何为terra栅格重现raster包的函数getValuesBlock？例如，如何将terra包的getValuesBlock文档中给出的示例翻译成raster包？在terra中有更直接的等价物吗？
浏览 39提问于2021-01-07得票数 0
1回答
几个星期以来，我一直在尝试将NetCDF文件放入数据帧中。尽管我成功地从ncdf文件中提取了变量/尺寸并绘制了切片，但在将其绑定到dataframe中然后绘制时，所有数据都会被隐藏。这里的最终目标是对数据帧进行栅格，以便能够按栅格对天气进行分类。 可以从https://cds.climate.copernicus.eu/cdsapp#!
浏览 32提问于2021-10-22得票数 0
我正在尝试使用terra::extract()来获取地块中的平均NDVI值。我注意到来自terra::extract()的数据帧给出了NA值(实际上，在4519个观察值中，我有2494个NA值)。我也尝试过exactextractr::exact_extract()，但它没有给我提供与raster::extract()相同的值。我想知道使用哪个更好，exactextractr::exact_extract()还是<e
浏览 1提问于2021-03-01得票数 0
1回答
用户应该输入一个n值来选择列和行数(或窗口大小)。将返回一个带有值的光栅。下面是一个可重复的例子：r <- terra::set.ext(rast(volcano), terra::ext(0, nrow(volcano), 0, ncol(volcano为此，我想将min和max范围(使用terra包)作为下面的代码：
w_10 <- c(xmin(terra::ext(r)),
浏览 5提问于2022-03-31得票数 1
回答已采纳
我有一个关于用最大重叠对多边形进行栅格化的问题，即指定与栅格单元重叠面积最高的多边形的值。
现实世界的练习是在R中对土壤is的多边形进行栅格化，以产生相对较低分辨率的土壤性质图作为模型输入。问题是，terra包的rasterize()函数(以及类似恒星的st_rasterize())从包含细胞中点的多边形分配细胞值。如果光栅单元格包含多个多边形，我宁愿选择多边形的值(土壤ID)，它在光栅单元格中具有最高的aerea覆
浏览 18提问于2022-02-10得票数 6
回答已采纳
UPDATE --这个问题与terra 1.4-1有关，现在terra 1.5-12已经过时了(或者更早)。我在terra::gridDistance上找到了一个看上去很旧的网页，来自terra v0.2-8 ，但据我所知，terra::distance是raster::gridDistance ()的最新替代品但是，我不知道如何在omit = my_omit中实现terra::distance (或等效的)。从看，似乎任何非NA都被
浏览 2提问于2021-09-02得票数 0
回答已采纳
1回答
我需要在一个网格单元上总结我的SpatVector，涵盖这10点，并思考地：：栅格化是最适合这样做的功能。栅格化的输出需要考虑SpatVector的两个属性。但是，我只是想不出如何将属性传递给函数。例如，使用这个奇数函数： out <- mean(i_1*i_2, na.rm = T)}
返回具有1行1列(
浏览 4提问于2022-02-28得票数 1
1回答
我是第一次接触terra套餐。我正在尝试aggregate一个只有一个图层的分类栅格(或者更准确地说，是一个SpatRaster )。结果应该是一个具有与原始栅格中的类别一样多的图层的栅格；单元的值应该具有每个类别中的原始(较小)单元的数量。下面的例子展示了我想要实现的目标： library(terra)
set.seed(0)- c(&q
浏览 33提问于2021-10-12得票数 1
回答已采纳
1回答
我有一个由rast()从terra包创建的SpatRaster文件，我想以矩阵的形式读出栅格值。如果底层的栅格有10行10列，我希望这些值的格式相同，但作为一个矩阵。
浏览 5提问于2021-04-15得票数 0
回答已采纳
1回答
我需要通过使用z <- rast(x, vals = y)分配值来创建一些新的栅格。我希望新的光栅(z)在与旧光栅(x)相同的位置上有NA值，但是在非NA单元中有新的赋值值(S) (y)。，什么是创建一个新栅格的最有效的方法？
我没有看到任何参数来处理rast()文档中列出的NA值，所以是否需要执行两个步骤，首先将值分配给所有单元格，然后将NAs从旧的栅格转移到新的单元格？在下面的示例中，我希
浏览 5提问于2022-04-14得票数 0
回答已采纳
2回答
然而，我发现栅格数据集的库/功能非常有限。有几个光栅/火花图书馆，但他们不是很成熟。例如、和。因此，我一直在探索在Databricks平台上高效处理栅格数据的其他方法，该平台利用Spark / Delta表/ Parquet文件。我的一个想法是将栅格数据转储到简单的x，y，value列，并将它们作为表加载。python gdal2xyz.py -band 1 -skipnodata "AR_FLRF_UD_Q1500_RD_02.tif" "AR_FLRF_UD_Q150
浏览 11提问于2022-09-27得票数 0
回答已采纳
1回答
我试着用光栅层计算一个类似于转动惯量的度量，我很难弄清楚每个细胞到一个贴片质心的距离，然后提取这个距离和这个细胞的值。我想要计算转动惯量(求出每个单元到其贴片质心的平方距离，乘以单元格的值，将这些值乘以补丁，然后除以每个块的所有值之和)。下面我提供了一个简化的设置。代码创建了一个简单的栅格层、片状的单元格，并获取它们的质心。我知道接下来要用到的函数可能是terra：：距离(可能与terra：：zonal)结合使用library(<em
浏览 11提问于2022-05-08得票数 1
回答已采纳
我正在尝试使用terra::distance，就像使用raster::distanceFromPoints一样。然而，terra::distance只报告从点到NA单元的距离。这是预期的结果吗？我在解决方案中包括了示例代码。距离计算用带点的栅格图地形图：：点的距离期望输出r <- terra::rast(ncols=10, nrows=10)vecP
浏览 3提问于2021-12-08得票数 1
它使用st_intersection工作，但是这需要很长时间，所以我尝试将多多边形转换为栅格，然后用较小的多边形进行裁剪。
浏览 7提问于2022-10-26得票数 0
1回答
我在一个shapefile中包含了多个多边形，并且我使用rgeos中的gCentroid()函数来创建一堆质心点。简单的例子：library(rgeos)
v <- vect(system.file("ex/lux.shp", package="terra"问题是，我最终将使用这些质心从栅格中导出值，并需要is来区分这些值来自哪个多边形。
浏览 11提问于2021-06-08得票数 0
我有一个NDVI(Normalized植被指数)光栅对象，现在我想在相同的latitude.Like上计算像素值，我想知道30°N latittude.How中的像素值，我能用R做这个吗？但我发现在“栅格”或"terra“包中没有任何功能可以做到这一点。
浏览 5提问于2022-11-08得票数 0
回答已采纳
2回答
我需要替换下面示例中的栅格对象(r)中的NA。- system.file("external/test.grd", package="raster")我还试图删除这些(并将结果放在data.frame中)
浏览 4提问于2012-08-15得票数 14
回答已采纳
问题：如何格式化Terra包中预测函数的“const”参数？我使用Terra的预测功能从模型中创建空间地图。当我为我的三个常量值创建完整的栅格时，这个函数就像预期的那样工作。不幸的是，help文档没有指定如何构造数据以外的数据，而这些数据应该位于dataframe中。 data.frame.可以将常量值作为预测变量添加，这样您就不需要为它
浏览 1提问于2022-02-08得票数 2
我有一个关于app()系列函数的简单问题 我想用一个矢量(n个值)除以一个栅格，得到n个结果栅格 s <- rast(system.file("ex/logo.tif", package="terra"))[[1]]
lapp(s, aa = c(1, 10, 100, 1000), function(x, aa){0.46 * aa/(x + 1)}) 同样的，我也想知道如何使用SpatRasterDataset
浏览 18提问于2021-07-24得票数 1
额外的问题:检查lc_code变量在landcover_classes_cavmVec中是否与原来的1:21 lc_2003级别中的等效数字匹配是一个很好的策略？我不知道如何验证这些代码是否实现了我想要的结果。
浏览 12提问于2021-12-21得票数 2
回答已采纳
云服务器
ICP备案
对象存储
云直播
腾讯会议
腾讯云开发者
