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

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

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

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

val ‹- extract (174.1,53.9) Ошибка в .local (x, y, ...): экстенты не перекрываются

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 месяцев.

Спасибо,


extract raster r-raster ncdf4
person Community    schedule 02.05.2019    source источник

Ответы (1)


arrow_upward
1
arrow_downward

Функция 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
person Michael Dorman    schedule 02.05.2019