Я пробовал это для извлечения данных из trmm ежедневных файлов nc4.
path<-"F:/TRMM/test"
setwd(path)
## open all `trmm` rainfall data files from folder
trmm_files = list.files(path,pattern = '*.nc', full.names = TRUE)
#I want only HQprecipitation variable from nc4 file.
trmm_stack <- stack(trmm_files, varname = "HQprecipitation")
Но trmm_stack потерял свою систему координат при штабелировании. Следовательно, мне пришлось его назначить. Кроме того, как я могу проверить, отсортирован ли растровый стек по датам и в каком порядке? Мне нужны данные временного ряда.
crs(trmm_stack) <- CRS('+init=EPSG:4326') # also it was rotated
trmm_stack_flip <- t(flip(trmm_stack, direction='y'))
Затем я создаю матрицу координат широты-долготы для местоположений
lat<- c(26.0166666, 25.674444, 25.5666666, 25.6666666, 26.18, 25.57)
lon <- c(91.85, 91.9025, 91.8666666, 91.8833333, 91.75, 91.9)
xym <- cbind(lat, lon)
и извлекли значения ячеек из trmm_stack, используя
data2 <- extract(trmm_stack, xym, method='simple', buffer = NULL, small = FALSE, cellnumbers = FALSE,
fun=NULL, na.rm=TRUE, df = FALSE, factors = FALSE)
Затем я получаю выходную таблицу с именами столбцов, заканчивающимися номером слоя. Но я не знаю, было ли извлечение было выполнено в точном месте и в том порядке дат, как временной ряд, как я хотел. Как это проверить? Любая помощь будет оценена по достоинству.
