问题标签 [r-raster]

我在 R 中的一个循环中苦苦挣扎，我必须使用动态变量名（我被告知这是其他关于动态变量名的帖子的坏主意，但我很确定我需要基于我的文件结构）。循环进入的每个文件夹都有不同数量的文件。

动态变量名称包含矩阵，我需要查看矩阵的每一行/列并输出一个新矩阵。

精简示例：

我想从变量 3 的列表中获取第 f 个变量矩阵，并为 [g,h] 的位置分配一个值，我以前用真实的变量名做过这个，但我在努力添加动态元素。这就是它的样子和我得到的错误。

我尝试了 eval 和 assign 的各种组合来组织最后一个语句，但我得到了诸如“无效的第一个赋值”、不正确的维数和赋值目标扩展到非语言对象等错误。

谢谢你的帮助！

R 版本 3.1.1 “Sock it to Me”与库（光栅）

我想为我的地图制作一个自定义颜色键。我的数据涵盖负值和正值。这是我尝试做的事情：

结果：

问题：如何将值 0 与黄色匹配，然后将颜色键标签放置在中断的中点而不是边缘？我需要削减= 10。

谢谢，AT。

此处可重现的示例

我有从某人那里得到的这个 mat 文件，我想将它更改为 R 中的光栅砖。

导入没问题，我可以使用R.matlab包导入

列表中有纬度、经度和 Psim。

我曾经将光栅砖功能转换为光栅砖。但是，我不知道matlab中栅格的结构，它应该是24年的时间序列每日降雨网格。但我不知道为什么，我无法从这个 mat 文件中识别出层数。我试着像阵列中的砖一样，比如

但是，我发现了这个错误。

无论如何可以帮助我吗？最好的祝福

我有一个要拼接在一起的栅格列表。除了 utm 区域外，投影是相同的。这是 UTM 区域 50 中一个 rasterStack 的坐标线：

这是另一个具有相同投影的 rasterStack，但它是 UTM 区域 51：

当我尝试在不先重新投影的情况下进行镶嵌时（就像在其他 GIS 中一样），我收到以下错误：

这让我相信我需要重新投影所有的 rasterStacks。这在计算上很耗时projectRaster，但我相信我可以这样做：

我是否正确假设我必须先重新投影？感谢您的任何指导。

我一直在尝试开发一种更好的解决方案来使用 R 可视化我的数据。但是，到目前为止，我的努力还没有取得成果。

假设您有 6 个栅格（数据是针对在不同时间测量的单个变量），其中栅格内变化不大，但栅格间变化较大。

6 个栅格需要一个颜色条。我所追求的是显示 6 个栅格，仍然能够看到栅格内部和栅格之间的变化，并创建一个带有中断标签的单一颜色键。

一个建议是重新调整数据并创建一个通用比例。就像是：

重新缩放数据后，栅格之间的变化减少并且看起来更好。

此外，如何生成使用未缩放数据的彩zlim=(min(s),max(s))条图例？s

大型可重现示例：

r1

r2

r3

r4

r5

r6

我有光栅图像，我想应用一个函数，该函数将根据相邻单元格的值给出一个值。在以该单元格为中心的 3x3 窗口中，我想给一个单元格一个介于 0 和 8/8（=1）之间的值，具体取决于多少个单元格（来自该窗口中的其他 8 个单元格）。例如，如果其他 5 个单元格的值与中心单元格不同，则函数必须在该单元格中分配 5/8。我正在尝试使用光栅包中的焦点来做到这一点。

首先，我编写了一个处理 3x3 窗口的函数：

然后我尝试像这样使用焦点：

r2<-focal(r, w=matrix(1,3,3), fun=mix(w))

但是：“混合错误（w）：找不到对象'v'”

我认为我遗漏了一些东西，也许方法不正确。

任何帮助将不胜感激。

提前谢谢约翰

有很好的示例如何在栅格图层上叠加矢量/栅格地图，例如在 R 中绘制叠加在栅格背景上的矢量地图或在 R中的背景 GM 图层上叠加栅格图层。但是合并的图形中缺少图例。在示例rproblv中有一个图例但migmaplv + rproblv没有。

有没有一种简单的方法可以组合两个levelplot或spplot对象并在结果中保留一个或两个图例？

从在 R 中的背景 GM 图层上覆盖您的栅格图层

我已经使用包 raster 和 rgdal 从 CSV 制作了 shapefile，但我想将这些 shapefile 剪辑到某个区域，我可以这样做。
但是，当我尝试 write.csv 时它不起作用。

另外，在 writeOGR 之后我丢失了 X 和 Y。

当我将 X、Y 返回到 shapefile 属性表中时，我想将栅格值提取到点，然后将点的坐标和新属性写入 .csv。

如何在最终产品中获得 X、Y？

我有一个 202 行和 201 列的栅格地图 在这个地图中有一些像素值为 0 的网格 我想编写一个函数来返回所有像素值 0 网格的坐标 我该怎么做 我试图使用 if loop 和 while 循环，但它总是说 TRUE/FALSE 需要这里是我的示例代码

我之前曾将 INLA 与 SpatialPolygonsDataFrame 一起使用过，但我从未将它与光栅一起使用过。我已经在互联网上搜索了好几个小时，阅读了小插图和我能找到的所有内容，但仍然无法弄清楚如何将 INLA 与光栅一起使用。下面只是我的数据的一个简单示例：

于是，我就站在你们面前，以这种方式乞求怜悯和帮助。