преобразовать данные modis EVI в дату yy-mm-dd в r - PullRequest
Новая
       20

преобразовать данные modis EVI в дату yy-mm-dd в r

0 голосов
/ 12 декабря 2018

Я работаю с растрами MODIS EVI в 2000 году. У меня есть 6 растров по годам, один растр по месяцам: 

"D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000209_aid0001.tif"
"D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000225_aid0001.tif"
"D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000241_aid0001.tif"
"D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000257_aid0001.tif"
"D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000273_aid0001.tif"
"D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000289_aid0001.tif"
"D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000305_aid0001.tif"
"D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000321_aid0001.tif"
"D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000337_aid0001.tif"
"D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000353_aid0001.tif"

Я хотел бы преобразовать месяц как 2000-02-09, но я нене знаю, как это сделать.

1 Ответ

0 голосов
/ 16 декабря 2018

Похоже, что ваш datetime формат строки соответствует year(4digits)day_of_year(3digit) формату, например, 2000209 означает 2000 год и 209 дней (27 июля 2000 г.).Если это правда, то проблема не сложна:

  • Извлеките эти семизначные числа.(str_extract)
  • Извлечь из него 'datetime'.(вам нужно знать, что j используется для разбора даты из day_of_year.)
  • [: Graph:] отбрасывает все, кроме цифр и знаков препинания.

Код 

dt <- data.frame(string =
     c("D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000209_aid0001.tif",
       "D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000225_aid0001.tif",
       "D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000241_aid0001.tif",
       "D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000257_aid0001.tif",
       "D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000273_aid0001.tif",
       "D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000289_aid0001.tif",
       "D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000305_aid0001.tif",
       "D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000321_aid0001.tif",
       "D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000337_aid0001.tif",
       "D:/Rteledetection/Pivots/MODIS/MOD13Q1.006__250m_16_days_EVI_doy2000353_aid0001.tif"))

dt$string %>% str_extract('\\d{7}') %>% str_replace('2000', '2000-') %>% 
   parse_date_time('y-j') %>% str_subset('[:graph:]')

Выход

[1] "2000-07-27" "2000-08-12""2000-08-28" "2000-09-13" "2000-09-29" "2000-10-15"
[7] "2000-10-31" "2000-11-16" "2000-12-02 "" 2000-12-18 "

...