用R读取栅格数据，格式为jp2格式

Question

我在阅读带有R的jp2格式的sentinel2A数据时收到错误Error in .local(.Object, ...)。我安装了rgdal软件包（rgdal_1.2-8）。

$platform---"x86_64-w64-mingw32"
$system---"x86_64, mingw32"
$version.string---"R version 3.4.1 (2017-06-30)"
gdal version-----"GDAL 2.0.1, released 2015/09/15"

Answer 1

试试raster <- readGDAL('you_file_name.jp2')。您可以在此处详细了解：https://gis.stackexchange.com/questions/211157/open-jpeg2000-sentinel-2-in-r

Answer 2

尝试将jp2文件转换为tif（GeoTif）格式

r代码：

library(rgdal)
gdal_translate("image_band_name.jp2", "new_band_name.tif")
band1 <- readGDAL("new_band_name.tif")