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

Question

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

我遇到的问题是：

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

不过，我已经更新了我安装的rgdal包，并使用了支持JPEG2000的更新的gdal二进制文件，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：行情)表示，但似乎还没有明确的答案.

Answer 1

使用这里介绍的命令行实用程序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中做这么多繁重的工作并没有确凿的理由，特别是通过本质上是包装器的东西。