Извлечение данных из NetCDF

Question

Я скачал температуру поверхности моря за январь отсюда https://oceancolor.gsfc.nasa.gov/l3/

и импортировал его в R.

Я знаю, как обрезать, используя экстент (ymax, ymin, xmax, xmin), но я не могу понять, как сделать это только для одной станции (53.9S, 174,1W) или ближайшей к этой координате. Можно ли обрезать данные только для одной станции?

val <- выдержка (174,1,53,9) Ошибка в .local (x, y, ...): экстенты не перекрываются </p>

SST_Jan <- brick("~https://oceandata.sci.gsfc.nasa.gov/cgi/getfile/A20021822018212.L3m_MC_SST_sst_9km.nc", stopIfNotEqualSpaced = FALSE, varname = "sst")

print(SST_Jan)

val<-extract(174.1, 53.9)

SST_Jan_station <- extract(SST_Jan, val)

Я хотел бы иметь возможность составить график изменений SST в этом конкретном месте за 12 месяцев

Спасибо,

Answer 1

Функция extract не работает с числовым вектором.

Вы можете поместить координаты в matrix -

pnt = matrix(c(174.1, 53.9), ncol = 2)
pnt
##       [,1] [,2]
## [1,] 174.1 53.9

И тогда extract будет работать -

extract(SST_Jan, pnt)
##      layer
## [1,]  8.24