ERROR: Couldn't find 'GRIB' or 'TDLP' Is the JPEG2000 driver available?

[nickleb]

Здравствуйте. Пытаюсь
в К 3.4.4 под Winfows 10
прочесть grib-файл:

Код: Выделить всё

library("raster")
r <- raster(FileName)

Говорит:

Код: Выделить всё

ERROR: Couldn't find 'GRIB' or 'TDLP'
There were 4 trailing bytes in the file.
dec_jpeg2000: Unable to open JPEG2000 image within GRIB file.
Is the JPEG2000 driver available?

Как "прорваться"?

[nickleb]

...попытка приложить файл в 18 Mb не удалась -
"Файл слишком большого размера"

[voltusya]

Может rNOMADS ?
А так, по заветам freeExec

Может оказаться, что GDALскомпилирован без его [JPEG2000] поддержки.

Стоит проверить rgdal

[nickleb]

voltusya,

Может rNOMADS ?

да:
по rNOMADS::GribInfo - идёт,
а далее по rNOMADS::ReadGrib
"виснет' без сообщений - вчера заряжал...

Поэтому сегодня и через raster стал пробовать....
Вероятно, надо "не ширше, а глубже"

по заветам freeExec

Спасибо!!!

[voltusya]

Есть у вас тут "JPEG2000" ?

Код: Выделить всё

rgdal::gdalDrivers()$name

У меня есть и так grb2 читается (но может у нас отличается grib как-то принципиально)

Код: Выделить всё

library("raster")
r <- raster(FileName)

[nickleb]

Есть у вас тут "JPEG2000" ?

voltusya, "ужо усё у порядке, Шеф" rgdal переустановил с новыми драйверами и драйверами на новые форматы, в т.ч. и JPEG2000 ! На "испытуемом" файле предупреждение и ошибка:

Код: Выделить всё

Warning: Inside GRIB2Inventory, Message # 171
ERROR: Couldn't find 'GRIB' or 'TDLP'
There were 4 trailing bytes in the file.
Warning: Inside GRIB2Inventory, Message # 171
ERROR: Couldn't find 'GRIB' or 'TDLP'
There were 4 trailing bytes in the file.

- но на др. - "родненьких" NOAA-шных - всё "гладенько".... подозреваю, что"испытуемый" файл какой-то производный "мутно-мутированный" - докопаюсь, поспрашав у авторов его...

[nickleb]

У меня есть и так grb2 читается (но может у нас отличается grib как-то принципиально)

voltusya, - "правильные" файлы посредством raster правильно теперь читаются и у меня! ещё раз спасибо, что надоумили rgdal переустановить! процитированному Вами мудрейшему из мудрейших freeExec - отдельное спасибо!

[nickleb]

КОД: ВЫДЕЛИТЬ ВСЁ

rgdal::gdalDrivers()$name
У меня есть и так grb2 читается (но может у нас отличается grib как-то принципиально)
КОД: ВЫДЕЛИТЬ ВСЁ

library("raster")
r <- raster(FileName)

...да, voltusya, вполне сейчас удобно в R метаданные dataset'ов файлов grib2 извлекать через rNOMADS:

Код: Выделить всё

Grb2Info <- rNOMADS::GribInfo(FileName)

и с самими dataset'ами работать "скормив" их raster:

Код: Выделить всё

Grb2Brick <- raster::brick(FileName)