我目前正在尝试从各种“ASC”文件中提取信息到R中,以便对数据进行分析。 问题是我不确定该如何读入这些文件。我尝试了标准的read.table函数,但所有数字都完全相同(-9999.00)。为了排除数据损坏的可能性,我读入了另一个ASC文件,得到了相同的结果。我唯一确定的是它们之间的文件大小...
我正在使用R语言的"raster"包从栅格文件创建一些地图。我想创建比较地图,将多个地图并排显示。这对于所有地图使用相同的颜色刻度很重要,而不管每个地图中的值如何。例如,如果地图1的值从0到1,地图2的值从0到0.5,则值为0.5的单元格在两个地图上应具有相同的颜色。 例如: - 地图1的...
我一直在使用 Raster 和 RMarkdown 记录一些空间数据的可视化过程,但是每个图像上方都有很多负空间,这让我感到困扰。以下是 RMarkdown 代码(有点简化): ```{r global_options, include=FALSE} knitr::opts_chunk$se...
我希望能够使用R来识别栅格地图上的线性特征,例如道路和河流,并将它们转换为线性空间对象 (SpatialLines类)。可以使用raster和sp包将栅格的特征转换为多边形矢量对象(SpatialPolygons类)。函数rasterToPolygons()将从栅格中提取某个值的单元格并返回一...
我正在使用Landsat图像计算净辐射。我的转换后反射率栅格层中有非常低的负值(例如:-0.000003)。为了减少未来计算中的误差,我想确保我的反射率在0到1之间。 在R中,如何将小于0的栅格值替换为“NA”,类似于栅格计算函数。我不确定如何提供示例,但我相信你们中的某个人可以帮助我,是吗?...
我需要移除R包"raster"中图像周围的框架,但是我不知道该改变哪个参数。以下是示例:library(raster) r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) plot(r,...
在 R 中,我们可以使用函数 rasterToCountour 将一个 raster 转换成一个 SpatialLinesDataFrame。library(raster) f <- system.file("external/test.grd", package="raster") r...