[GIS] Converting Sentinel 2 jp2 bands to TIFF with gdal

gdalgdal-translate

I want to convert a .jp2 file from the Sentinel 2 project to .tif with gdal. For some bands it works (e.g. B01) and for some bands it always fails (e.g. B04). I use the sentinel public image browser data (e.g this product).

When I run gdal_translate B04.jp2 B04.tif I get the following output:

Input file size is 10980, 10980
0maximum number of samples exceeded (120560400 > 67108864)
error: cannot decode code stream
ERROR 1: B04.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0

My gdal version is 2.3.0, released 2018/05/04

Best Answer

There is no issue in conversion of B04.jp2 (you posted in question) to B04.tif using the above command.

I used gdal 2.3.0 released 2018/05/04. It worked well.

Try opening the data in a viewer. It might be corrupted/not fully download.

The error IReadBlock failed generally occurs when the file is not complete.

Related Solutions

[GIS] Converting large jp2 files using GDAL

Jasper driver is useless for big JPEG2000 files. OpenJPEG driver is OK if source JPEG2000 is tiled, otherwise you will need JP2ECW, JP2KAK, or JP2MrSID driver. Check your drivers with gdalinfo --formats and use --config GDAL_SKIP for skipping the first JP2 drivers from the list for selecting the driver you want to use.

Example output from GDAL compiled with two JP2 drivers:

gdalinfo --formats
...
JP2ECW (rov): ERDAS JPEG2000 (SDK 5.0)
JP2OpenJPEG (rwv): JPEG-2000 driver based on OpenJPEG library

If you wand to use JP2OpenJPEG driver which is the second on the list you must skip the first driver:

gdal_translate -of GTiff input.jp2 output.tif --config GDAL_SKIP JP2ECW

[GIS] Convert TIFF file with multiple bands into a single CSV file

AFAIK gdal_translate is able to convert only one band into XYZ format, so I suggest the following workflow:

use gdal_translate to convert each band into a CSV file;
use ogr2ogr to append each CSV file to an SQLite file;
use ogrinfo to create the index on the ogc_fid field of each SQLite table in order to make the last step fast;
finally, join all the tables and get the desired result.

This is a sample batch file where the above steps are implemented for only three bands (change the FILENAME variable and adapt the join query according to your input data):

set FILENAME=test
if exist %FILENAME%.sqlite del %FILENAME%.sqlite
for /L %%i in (1,1,3) do (
gdal_translate -b %%i -of XYZ %FILENAME%.jpg b%%i.csv -co ADD_HEADER_LINE=YES
ogr2ogr -update -append -f SQLite %FILENAME%.sqlite -nln b%%i b%%i.csv -dsco METADATA=NO -dsco INIT_WITH_EPSG=NO
ogrinfo %FILENAME%.sqlite -sql "CREATE INDEX b%%i_ogc_fid_index ON b%%i(ogc_fid)"
)
ogr2ogr -f CSV -sql "SELECT b1.x AS lon, b1.y AS lat, b1.z AS band1, b2.z AS band2, b3.z AS band3 FROM b1 LEFT JOIN b2 ON b1.ogc_fid = b2.ogc_fid LEFT JOIN b3 ON b1.ogc_fid = b3.ogc_fid" %FILENAME%_XYZ.csv %FILENAME%.sqlite -lco SEPARATOR=SPACE

Best Answer

Related Solutions

[GIS] Converting large jp2 files using GDAL

[GIS] Convert TIFF file with multiple bands into a single CSV file

Related Question