Как открыть Himawari-8 ежечасно в режиме реального времени набор данных SST? - PullRequest
0 голосов
/ 09 марта 2020

У меня есть временные ряды данных SST, почасовые данные которых получены с помощью спутника Himawari-8 ( здесь ). К сожалению, у меня есть выданный при открытии этого набора данных в растровый объект. мой код выдает сообщение об ошибке:

r<-raster('d:/sst_nrt/20200301040000-JAXA-L3C_GHRSST-SSTskin-H08_AHI_NRT-v1.2-v02.0-fv01.0.nc', varname='sea_surface_temperature')
Error in .rasterObjectFromCDF(x, type = objecttype, band = band, ...) : 
cells are not equally spaced; you should extract values as points

Затем я пишу дополнительный аргумент stopIfNotEqualSpaced = FALSE и выдаю вывод сетки с разборкой:

> r
class      : RasterLayer 
dimensions : 6001, 6001, 36012001  (nrow, ncol, ncell)
resolution : 0.05999667, 0.02  (x, y)
extent     : -180.01, 180.03, -60.01, 60.01  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : d:/sst_nrt/20200301040000-JAXA-L3C_GHRSST-SSTskin-H08_AHI_NRT-v1.2-v02.0-fv01.0.nc 
names      : sea.surface.skin.temperature 
z-value    : 1235880017 
zvar       : sea_surface_temperature 

здесь график ( рис. 1), и дать ложное / другое значение долготы сравнить с официальным сайтом Химавари (рис. 2)

fig. 1Fig. 2

...