如何使用此R函数拼接相同的HDF文件？

Question

文件夹中有1,000多个MODIS HDF图像：

M:\join

他们的名字告诉我们哪些文件必须拼接在一起。

例如，在以下文件中， 2009090表示这三个图像必须镶嵌在一起：

MOD05_L2.A2009090.0420.051.2010336084010
MOD05_L2.A2009090.0555.051.2010336100338
MOD05_L2.A2009090.0600.051.2010336100514

或这两个日期是同一日期2009091 ：

MOD05_L2.A2009091.0555.051.2010336162871
MOD05_L2.A2009091.0600.051.2010336842395

我将使用此功能（功能来源 ）将它们镶嵌：

mosaicHDF(hdfNames, filename, MRTpath, bands_subset, delete=FALSE)

如何将我的HDF文件引入hdfNames ？

我应该用文件名写什么？

我试图找到此功能的手册，但没有任何内容。

谢谢你的帮助。

Answer 1

幸运的是，我确实找到了答案。 谢谢你的帮助。

hdfs <- c('MOD05_L2.A2009090.0420.051.2010336084010.hdf',
          'MOD05_L2.A2009090.0555.051.2010336100338.hdf',
          'MOD05_L2.A2009090.0600.051.2010336100514.hdf')

mosaicHDF(hdfNames=hdfs, filename='newhdf.hdf', MRTpath='C:/MRT/bin',bands_subset="1 0 0 0", delete=FALSE) 

但是我有一个新问题:-)

由于该文件夹中有数千个HDF文件，

如何编写循环以将所有HDF文件引入该函数？

仅供参考：我对R很陌生。

Answer 2

这个问题已经很老了，但我认为我会张贴要求的R代码Canada2015。 因此，对于2000年至2016年的循环，假设文件名仍为A <YEAR> <MONTH> <DAY>。 格式。 此代码将所有图块拼接在一起，以产生一个用于当前年份的新文件。 如果您必须在一年中拼接许多图块，则可以将pattern =参数更改为类似“ .hdf”的常规名称

for(i in 2000:2016){
  HDFs <- list.files(path = "F:/PATHTOFILES/HDFs/", pattern = paste("A",i,sep = ""))
  mosaicHDF(hdfNames = HDFs, filename = paste('newhdf',i,'.hdf',sep = ""), MRTpath = 'C:/MRT/bin',bands_subset="1 0 0 0", delete=FALSE)
}