Thank, but didn't work.

> test = readGDAL("C:/Users/   ...  M_B03.jp2")

      Error in .local(.Object, ...) :

     `C:\Users\  ...  B03.jp2' not recognised as a supported file format.

you can follow this link

http://gis.stackexchange.com/questions/73849/jp2-fails-to-open-in-qgis-and-gdal

Hi Sayan, thanks but it didn't worked because there is no plugins like that for R. =(

I'd suggest you a two steps process.

First with gdal_translate, to obtain tif from jp2. USGS provides a number of solutions to batch convert these file formats. See here:

https://viewer.nationalmap.gov/tools/rasterconversion/Convert_JP2_Data_to_GeoTIFF_with_GDAL.html

Once with the tif files, stack() them directly into R.

I've successfully converted sentinel-2 jp2 to tif using GDAL_JP2_2_GeoTIFF.bat 

https://viewer.nationalmap.gov/tools/rasterconversion/Convert_JP2_Data_to_GeoTIFF_with_GDAL.html

Thank you Paulo.

Jonathan, gdal_translate is failing to retrieve crs from jp2 files. Please let me know it you found a solution for this. (This could be solved from R side but anyway...).

All the best

Hi Paulo, I tried the bat file, it converted all jp2 into tif, however i cannot use stack() function since sentinel2 bands have different resolution. Do I need to resample them first before doing stack()? any suggestions?

Thanks

Hi Yask.

To build a rasterstack all layers must share the same extent, resolution, origin, crs.

see ?stack ?compareRaster

Did you try to resample those with higher resolution?

It is easier in R when you convert the ".jp2" files to ".tif"

Hi all

Just wanted to get back to this thread.

Does anyone know how do we atmospherically correct sentinel 2 data processing level 1C data?

I downloaded this from sentinel hub and it is at processing level 1C.

Can we do this in R?

Hi Paulo Eduardo, the link that you mentioned leads to an "unavailable page". Could you possibly share the script please?

Julius, I still being able to access https://viewer.nationalmap.gov/tools/rasterconversion/Convert_JP2_Data_to_GeoTIFF_with_GDAL.html

You could write a batch (.bat) file with something like:

for %%N in (path/to/your/sentinel2/*.jp2) DO C:/Program Files/QGIS 2.14.5/bin/gdal_translate.exe -of GTiff C:/temp/sentinel2/%%~nN.tif

@ECHO OFF

SETLOCAL EnableDelayedExpansion

SET mypath=C:\temp\sentinel2\

FOR /F %%i IN ('DIR /B %mypath%*.jp2') DO (

SET infile=%%i

SET outfile=!infile:.jp2=.tif!

C:/Program Files/QGIS 2.14.5/bin/gdal_translate.exe -of "GTiff" %mypath%!infile! %mypath%!outfile!

)

You'll need to adjust it and redirect to your gdal_translate.exe and jp2 folders.

Thank you Paulo!

I'm use follows steps (https://stackoverflow.com/questions/40044403/gdal-jp2-driver-missing-in-python)

require(rgdal)

require(gdalUtils)

gdal_chooseInstallation('JP2OpenJPEG')

its works

Thank you Masita! That solves all I wanted!

Dear try to use Land viewer to get sentinel image with all bands not JPG

Hi dear , am using the atmospherically corrected sentinel2 level 2A image with .tif format , and when I tried to run the plotRGB function using R , I got a black image ,,,, any one with solution in this regard??

thank you

Hi, I just had the same problem and this was the easy way to solve it:

- download the free software multispec (I am using it for basic remote sensing stuff): https://engineering.purdue.edu/~biehl/MultiSpec/

- open each jp2 image in the software

- save as geotiff

then:

library(raster)

# Blue2, Green(3), Red(4), NIR(8), 10m res, Sentinel-2

b2