如何在 R 中创建一个采用多个栅格并将它们马赛克的循环？并用特定名称保存它们

发布于 2025-01-10 00:08:56 字数 1233 浏览 2 评论 0原文

我有大约 2000 个世界不同地区的光栅文件。其中一些来自同一天，另一些来自不同日期。我想马赛克/合并同一天的所有文件，以按天创建“图片”（光栅）。

所有文件的名称都以“MOD13A3.A2000”开头，然后是捕获图像的年份，然后是其他不重要的内容。我想做的就是获取所有以相同名称开头的文件并将它们合并。

例如，我会在这里获取前 4 个文件，因为它们以 MOD13A3.A2000032... 开头，将它们合并并只有一个名为 MOD13A3.A2000032 的文件。然后获取其他 4 个以下文件，因为它们以 MOD13A3.A2000061... 开头，将它们合并并得到另一个名为 MOD13A3.A2000061 的文件。

MOD13A3.A2000032.h12v13.006.2015138123528 // MOD13A3.A2000

032.h13v13.006.2015138123527 //

MOD13A3.A2000032.h13v14.006.2015138123528 // MOD13A3.A2000

032.h14v14.006.2015138123526 //

MOD13A3.A2000061.h12v13.006.2015136111214 //

MOD13A3.A2000061.h13v13.006.2015136111225 //

MOD13A3.A2000061.h13v14.006.2015136111218 //

MOD13A3.A2000061.h14v14.006.2015136111220 //

.. 我一直在研究其他论坛，

但因为我完全R 编程新手我不太理解发布的代码。我曾考虑过创建一个 for 循环来更改文件名的开头，例如 "MOD13A3.A2000"+i 但我不明白它的语法。我也无法找到清晰的示例来说明如何像这样更改文件名。

如果有人可以解释代码，我将非常感激！

原文

分享到QQ

分享到微博

如果你对这篇内容有疑问，欢迎到本站社区发帖提问参与讨论，获取更多帮助，或者扫码二维码加入 Web 技术交流群。

发布评论

需要登录才能够评论，你可以免费注册一个本站的账号。

橘寄 2025-01-17 00:08:56

我将从 terra 包开始：

library(terra)

r <- list.files(path = "rasters", pattern = "MOD13A3.A2000061.*", full.names = TRUE)
r <- lapply(r, rast)
r <- do.call(merge, r)

第一行创建具有特定模式的光栅文件列表，第二行对它们应用 rast() 函数，第三行进行 merge.

填充 free 为文件名中的模式创建一个循环：list.files 读取文件，然后获取文件名的 substr() 并查找所有 unique() 字符串。

I would start with terra package:

library(terra)

r <- list.files(path = "rasters", pattern = "MOD13A3.A2000061.*", full.names = TRUE)
r <- lapply(r, rast)
r <- do.call(merge, r)

First line creates list of raster files with specific pattern, 2nd applies rast() function to them, and 3rd do merge.

Fill free to create a loop for patterns in your filenames: list.files reads the files, then take a substr() of filenames and find all unique() strings.

回复收藏 0 原文

~没有更多了~