gpt4 book ai didi

r - 如何有效地将多个光栅 (.tif) 文件导入 R

转载 作者:行者123 更新时间:2023-12-03 18:45:29 26 4
gpt4 key购买 nike

我是 R 新手,尤其是在空间数据方面。我试图找到一种方法将多个 (~600) 单波段光栅 (.tif) 文件有效地导入到 R 中,所有文件都存储在同一文件夹中。不确定这是否重要,但请注意,在我的 Mac 和 Windows Parallel VM 上的文件夹中查看时,每个 .tif = .TIF 都有以下五 (5) 种文件格式; .tfw; .TIF.aux.xml; .TIF.ovr; .TIF.xml。无论如何,以下代码(以及我尝试过的其他类似变体)似乎不起作用:

library(sp)
library(rgdal)
library(raster)

#path to where all .tif files are located
setwd("/path/to/workingdirectory")

#my attempt to create a list of my .tif files for lapply
temp = list.files(pattern="*.tif")
temp #returns 'character(0)'

#trying to use the raster function to read all .tif files
myfiles = lapply(temp, raster)
myfiles #returns 'list()'

有没有办法使用某种形式的循环来有效地导入所有光栅文件?

最佳答案

我找到了答案,并将发布完整代码以帮助遇到此问题的其他初学者 R 用户。要调用列表元素,请使用双方括号 [[]],如下所示:

#first import all files in a single folder as a list 
rastlist <- list.files(path = "/path/to/wd", pattern='.TIF$',
all.files=TRUE, full.names=FALSE)

#import all raster files in folder using lapply
allrasters <- lapply(rastlist, raster)

#to check the index numbers of all imported raster list elements
allrasters

#call single raster element
allrasters[[1]]

#to run a function on an individual raster e.g., plot
plot(allrasters[[1]])

嘘。感谢冻糕的帮助。

关于r - 如何有效地将多个光栅 (.tif) 文件导入 R,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/52746936/

26 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com