26得票4回答
将ASC文件读入R语言

我目前正在尝试从各种“ASC”文件中提取信息到R中,以便对数据进行分析。 问题是我不确定该如何读入这些文件。我尝试了标准的read.table函数,但所有数字都完全相同(-9999.00)。为了排除数据损坏的可能性,我读入了另一个ASC文件,得到了相同的结果。我唯一确定的是它们之间的文件大小...

23得票6回答
如何在R中创建具有相同颜色刻度的光栅图

我正在使用R语言的"raster"包从栅格文件创建一些地图。我想创建比较地图,将多个地图并排显示。这对于所有地图使用相同的颜色刻度很重要,而不管每个地图中的值如何。例如,如果地图1的值从0到1,地图2的值从0到0.5,则值为0.5的单元格在两个地图上应具有相同的颜色。 例如: - 地图1的...

20得票3回答
PROJ4升级到PROJ6和“废弃数据”警告

Context 这些问题与从PROJ4升级到PROJ6所引起的更改及其对各种R空间包(sp,sf,raster)的影响有关。 我们现在收到了很多关于“已放弃基准”的警告,看起来有点令人担忧,我有点困惑我应该做什么。我可以看到这可能会在某些情况下产生可怕的后果,并且在其他情况下可以忽略它们。 看...

19得票1回答
转换 HTML 到 RMarkdown 时如何处理空格块?

我一直在使用 Raster 和 RMarkdown 记录一些空间数据的可视化过程,但是每个图像上方都有很多负空间,这让我感到困扰。以下是 RMarkdown 代码(有点简化): ```{r global_options, include=FALSE} knitr::opts_chunk$se...

19得票2回答
ggplot2:设置透明度值时,光栅绘图效果不如预期

这是我的第一篇帖子,希望我能够遵守网站礼仪。我在网站上找不到答案,之前我在一个ggplot2专门的群组中发布了这个问题,但至今没有解决方案。 基本上,我正在尝试使用ggplot2叠加两个raster,并要求顶部的raster半透明。我有一个hillShade栅格,它是从高程数据栅格计算出来的...

18得票4回答
在栅格地图上识别线性特征并使用R返回线性形状对象

我希望能够使用R来识别栅格地图上的线性特征,例如道路和河流,并将它们转换为线性空间对象 (SpatialLines类)。可以使用raster和sp包将栅格的特征转换为多边形矢量对象(SpatialPolygons类)。函数rasterToPolygons()将从栅格中提取某个值的单元格并返回一...

18得票5回答
如何在R代码中将小于0的栅格值替换为NA

我正在使用Landsat图像计算净辐射。我的转换后反射率栅格层中有非常低的负值(例如:-0.000003)。为了减少未来计算中的误差,我想确保我的反射率在0到1之间。 在R中,如何将小于0的栅格值替换为“NA”,类似于栅格计算函数。我不确定如何提供示例,但我相信你们中的某个人可以帮助我,是吗?...

18得票2回答
如何在R包“raster”中的“plot.raster”函数中移除盒状框架

我需要移除R包"raster"中图像周围的框架,但是我不知道该改变哪个参数。以下是示例:library(raster) r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) plot(r,...

18得票2回答
R - 将SpatialLines转换为栅格

在 R 中,我们可以使用函数 rasterToCountour 将一个 raster 转换成一个 SpatialLinesDataFrame。library(raster) f <- system.file("external/test.grd", package="raster") r...

17得票2回答
如何使用GDAL Python投影和重采样网格以匹配另一个网格?

澄清:我不知道怎么使用GDAL和Python将NOAA GTX偏移网格的一部分转换为垂直基准面变换。我想将一个网格(在这种情况下是深度和属性值的网格,但它也可以是一个Geotif文件)用作我想要进行操作的模板。如果我能够正确执行此操作,我有信心它将大大帮助人们利用这种类型的数据。 这里是我的...