Как читать JPEG2000 данных Sentinel-2 в R

Я пытаюсь использовать пакет rgdal для открытия нескольких данных Sentinel-2 в формате JPEG2000. .

Проблема, с которой я сталкиваюсь, заключается в следующем:

Error in .local(.Object, ...) : 
d:/data/T35SMD_20161227T091402_B01.jp2' not recognised as a supported file format.

Хотя я обновил свою установку пакета rgdal и использую обновленные двоичные файлы gdal, которые поддерживают JPEG2000, в котором была введена поддержка в версии 1.9.0.

rgdal: version: 1.2-5, (SVN revision 648)
Loaded GDAL runtime: GDAL 2.0.1, released 2015/09/15

Когда я выполняю команду gdalDrivers(), запись для JPG2000 не отображается среди других драйверов.

Может кто-нибудь, пожалуйста, скажите мне, как я могу решить эту проблему? У меня много данных S-2, и конвертировать их в tiff в QGIS или любой подобной программе для меня не вариант.

Я также нашел несколько похожих вопросов по этой теме (см. ), однако прямого ответа пока нет.


r rgdal jpeg2000
person Ursulka Siete    schedule 13.01.2017    source источник

Ответы (1)


arrow_upward
2
arrow_downward

Используйте утилиту командной строки gdal_translate, как описано здесь: http://www.gdal.org/frmt_sentinel2. HTML. У вас будет гораздо больше гибкости и, вероятно, лучшая скорость и обработка памяти. Пример, который дан там:

gdal_translate SENTINEL2_L1C:S2A_OPER_MTD_SAFL1C_PDMC_20150818T101440_R022_V20150813T102406_20150813T102406.xml:10m:EPSG_32632 10m.tif -co TILED=YES --config GDAL_CACHEMAX 1000 --config GDAL_NUM_THREADS 2

На самом деле нет веских аргументов в пользу выполнения такой тяжелой работы в R, особенно с помощью чего-то, что по сути является оболочкой.

person Benjamin    schedule 09.02.2017
comment
Большое спасибо за ваше предложение, я включил gdal_translate() из gdalUtils в свой код, и он отлично работает. - person Ursulka Siete; 17.02.2017