8得票4回答
R:检测“主”路径并使用核函数删除或过滤GPS轨迹?

有没有办法过滤掉不属于主路径的部分?如您在图片中所见,我想删除被划掉的部分,同时保留主路径。我已经尝试使用zoo/rolling median,但没有成功。我认为我可以使用某种核来完成这个任务,但我不确定。我还尝试了不同的平滑方法/函数,但这些方法并没有产生期望的结果,反而使情况更糟。 数据中...

8得票1回答
如何在不使用辅助文件的情况下将栅格数据写入磁盘?

我将使用R语言中的Raster包将数据集写入ERMapper格式(.ers)文件,但我在生成的.aux.xml辅助文件方面遇到了问题(实际上我并不关心它)。 简单的例子: rst <- raster(ncols=15000,nrows=10000) rst[] <- 1.234...

7得票1回答
Rscript: 尝试使用raster::extract时为什么会出现“Error in UseMethod("extract_")”错误?

我尝试使用raster包中的extract方法从Raster*对象中提取值。 RStudioPrompt> jpnpe <- extract(jpnp, jpnb, fun = mean, na.rm = T) jpnp是栅格对象,jpnb是SpatialPolygonsDa...

7得票1回答
在terra包中读取多个栅格图层

我正在改变我的空间工作流程,使用terra包代替raster包。使用raster包时,我曾经直接将多个raster文件读入一个堆栈。 filelist_temp <- list.files(datapath("Climate/World Clim 1 yr Monthly Weathe...

7得票2回答
将边框添加到rasterImage

这是我使用rasterImage创建的渐变色图例: colfunc <- colorRampPalette(c("red", "blue")) legend_image <- as.raster(matrix(colfunc(20), ncol=1)) plot.new() ra...

7得票1回答
为什么terra::cellSize()和raster::area()会给出不同的栅格单元面积估计?

我刚刚注意到 terra::cellSize() 产生的单元格面积估计值与 raster::area() 不匹配。 首先,为什么这两种方法不能提供相同的答案?其次,哪个估算值最准确?请参见以下示例。 library(raster) #> Loading required packag...

7得票1回答
写入栅格输出文件的大小

我有一个函数,它将多波段图像作为raster brick对象读入,遍历各个波段进行各种计算,然后将raster写出为新的.tif文件。所有这些都很好,但新图像文件的文件大小大约是原始图像的四倍(我认为是因为原始图像有4个波段)。我想知道writeRaster()函数中是否有我不知道的参数,或者...