У меня есть временные ряды данных 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)