I am trying to perform stack() in R using enhanced vegetation index (EVI) time-series maps but it doesn't include all the .tiff files in my directory. Here's the list of files when I dir(files.evi)
:
[1] "MOD13Q1_EVI_2019_001.tif" "MOD13Q1_EVI_2019_017.tif"
[3] "MOD13Q1_EVI_2019_033.tif" "MOD13Q1_EVI_2019_049.tif"
[5] "MOD13Q1_EVI_2019_065.tif" "MOD13Q1_EVI_2019_081.tif"
[7] "MOD13Q1_EVI_2019_097.tif" "MOD13Q1_EVI_2019_113.tif"
[9] "MOD13Q1_EVI_2019_129.tif" "MOD13Q1_EVI_2019_145.tif"
[11] "MOD13Q1_EVI_2019_161.tif" "MYD13Q1_EVI_2019_009.tif"
[13] "MYD13Q1_EVI_2019_025.tif" "MYD13Q1_EVI_2019_041.tif"
[15] "MYD13Q1_EVI_2019_057.tif" "MYD13Q1_EVI_2019_073.tif"
[17] "MYD13Q1_EVI_2019_089.tif" "MYD13Q1_EVI_2019_105.tif"
[19] "MYD13Q1_EVI_2019_121.tif" "MYD13Q1_EVI_2019_137.tif"
[21] "MYD13Q1_EVI_2019_153.tif" "MYD13Q1_EVI_2019_169.tif"
However, when I do
evi <- list.files(files.evi, pattern = glob2rx("*MYD13Q1*.tif$"), full.names = TRUE)
all_evi_st <- stack(evi)
it only stack half of the files on my list:
class : RasterStack
dimensions : 668, 1105, 738140, 11 (nrow, ncol, ncell, nlayers)
resolution : 231.6564, 231.6564 (x, y)
extent : 11251086, 11507066, 1348703, 1503450 (xmin, xmax, ymin, ymax)
crs : +proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs
names : MYD13Q1_EVI_2019_009, MYD13Q1_EVI_2019_025, MYD13Q1_EVI_2019_041, MYD13Q1_EVI_2019_057, MYD13Q1_EVI_2019_073, MYD13Q1_EVI_2019_089, MYD13Q1_EVI_2019_105, MYD13Q1_EVI_2019_121, MYD13Q1_EVI_2019_137, MYD13Q1_EVI_2019_153, MYD13Q1_EVI_2019_169
Can anyone tell me what's happening and what's the workaround? I am new in remote sensing and still shabby with R.
Other info I can give about the time-series images is that they are from MODIS TERRA and AQUA time series derived from the 250 m 16-day composite vegetation index products (MOD13Q1/MYD13Q1).
print(files.evi)
andprint(evi)
?stack(dir(files.evi))
? The file list onevi
is clearly not what you want since it has not only half the layers of dir(files.evi) but also has different layers.