如何在R中读取Sentinel-2数据的JPEG2000

时间:2017-01-13 15:22:48

标签: r rgdal jpeg2000

我正在尝试使用rgdal包打开JPEG2000格式的多个Sentinel-2数据。

我遇到的问题是:

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

虽然,我已更新了rgdal软件包的安装,并使用了支持JPEG2000的更新的gdal二进制文件,该文件支持introduced in v1.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数据,并在QGIS中将它们转换为tiff,或者任何类似的程序都不适合我。

我也发现很少有关于这个主题的类似问题(c.f。open jpeg2000 sentinel 2 in r),但似乎还没有直接的答案。

1 个答案:

答案 0 :(得分:2)

使用命令行实用程序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中做那么多繁重的工作并不是一个坚实的论据,特别是通过一些本质上是包装的东西。