在R中循环浏览nc文件

Question

大家早上好，我目前在这个链接上使用Antonio Olinto Avila da Silva编写的代码：https://oceancolor.gsfc.nasa.gov/forum/oceancolor/topic_show.pl?tid=5954它允许我从nc文件中提取sst/chlor_a类型的数据。它使用一个循环来创建一个包含所有数据的excel文件。不幸的是，我注意到该函数只获取循环中的第一个数据文件。因此，我发现自己在excel文件中一行有20倍相同的数据。有人能解决这个循环正常工作的问题吗？

Answer 1

我首先要检查这两行是否包含您期望的所有文件：

(f <- list.files(".", pattern="*.L3m_MO_SST_sst_9km.nc",full.names=F))
(lf<-length(f))

然后在for循环中出现了一个错误。此行：

data<-nc_open(f)

需要引用迭代器i，所以将其更改为以下内容：

data<-nc_open(f[[i]])

看来两个脚本都有相同的错误。