У меня есть следующие данные:
https://drive.google.com/file/d/1AZph5WjBjq_e_N1QqoTHzvUNuX4J-SUh/view?usp=sharing
Я хочу создать пространственный_пункт данныхDataframe, состоящий из широты, долготы, переменной времени и значений орошения на единицу,
library(ncdf4)
ncin <- nc_open("cons_irr_h08.nc")
lon= ncvar_get(ncin,'lon')
lat= ncvar_get(ncin,'lat')
head(cbind(lat,lon))
cons_irr <- ncvar_get(ncin, 'cons_irr')
plot(lon,lat,type='p',pch='.')
ncin <- nc_open("cons_irr_h08.nc")
tvar= ncvar_get(ncin,'month')
add.months= function(date,n) seq(date, by = paste (n, "months"), length = 2)[2]
times = as.Date(sapply(tvar-1,add.months,date=as.Date('1971-01-15')),origin='1970-01-01')
С помощью этого кода я могу получить все необходимые мне компоненты, но теперь я не знаю, как объединить их вместе в пространственный объект PointDataframe.
https://gis.stackexchange.com/questions/303855/find-time-variable-in-raster-data/303941?noredirect=1#comment490105_303941 в этом сообщении Дэйв помог мне с кодом, который я публикую здесь.