2

I'm trying to convert Sentinel2 jp2 image to tiff using gdal_translate. When running the code:

Miguel$ gdalinfo /Users/Miguel/Desktop/S2_b1.jp2 
Driver: JP2ECW/ERDAS JPEG2000 (SDK 5.2)
Files: /Users/Miguel/Desktop/S2_b1.jp2
Size is 1830, 1830
Coordinate System is:
PROJCS["WGS 84 / UTM zone 33N",
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AXIS["Latitude",NORTH],
    AXIS["Longitude",EAST],
    AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",15],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
    AUTHORITY["EPSG","9001"]],
AXIS["Easting",EAST],
AXIS["Northing",NORTH],
AUTHORITY["EPSG","32633"]]
Origin = (600000.000000000000000,7000020.000000000000000)
Pixel Size = (60.000000000000000,-60.000000000000000)
Metadata:
ALL_COMMENTS=Kakadu-v7.4
Kdu-Layer-Info: log_2{Delta-D(squared-error)/Delta-L(bytes)}, L(bytes)
-192.0,  3.8e+06

COLORSPACE=GREYSCALE
COMPRESSION_RATE_TARGET=3
Corner Coordinates:
Upper Left  (  600000.000, 7000020.000) ( 16d58'56.33"E, 63d 6'56.41"N)
Lower Left  (  600000.000, 6890220.000) ( 16d55' 3.52"E, 62d 7'50.52"N)
Upper Right (  709800.000, 7000020.000) ( 19d 9'13.82"E, 63d 4' 7.06"N)
Lower Right (  709800.000, 6890220.000) ( 19d 1' 7.24"E, 62d 5' 8.11"N)
Center      (  654900.000, 6945120.000) ( 18d 1' 5.23"E, 62d36'15.28"N)
Band 1 Block=256x256 Type=UInt16, ColorInterp=Gray
Description = Grayscale
Overviews: 915x915, 457x457, 228x228
Image Structure Metadata:
NBITS=15


ns27:~ Miguel$ gdal_translate /Users/Miguel/Desktop/S2_b1.jp2       /Users/Miguel/Desktop/S2_b1tif.tif
Input file size is 1830, 1830
0

The file is created but can´t be opened.

Any suggestion about how to proceed?

  • What means "file is created but can't be opened"? Is it just created so that progress percentage stays at zero? And gdalinfo does not report anything? – user30184 Jan 28 '17 at 12:57
  • When running gdal_translate the percentage stays at zero and a tif file is created in the folder but it is just a 8kb file that can not be opened (in my case I want to import it to R as tif). I'm using macOS Sierra v. 10.12.1. I'm using the Terminal to run the code. When closing Terminal, it seems gdal_translate is still working. I waited for hours but the percentage didn't move from 0 – mcg Jan 28 '17 at 17:51
  • Can you give us your exact GDAL version - as seen in gdalinfo --version? – Kersten Jan 28 '17 at 18:48
  • Which drivers your GDAL has for JPEG2000? Check with gdalinfo --formats and search for JP2 in driver names. – user30184 Jan 28 '17 at 18:50
  • my GDAL version is GDAL 2.1.2. This is what I see when checking the driver for JPEG2000 JPEG2000 -raster,vector- (rwv): JPEG-2000 part 1 (ISO/IEC 15444-1), based on Jasper library If I run this, the drivers are: Miguel:~ Miguel$ gdalinfo --formats | grep JP2 JP2ECW -raster,vector- (rw+v): ERDAS JPEG2000 (SDK 5.2) JP2OpenJPEG -raster,vector- (rwv): JPEG-2000 driver based on OpenJPEG library – mcg Jan 28 '17 at 20:44
5

I can't reproduce your issue with GDAL 2.2-dev (yesterday's trunk version from gisinternals.com).

I downloaded a rando Sentinel2 zip file from https://scihub.copernicus.eu/dhus/#/home. It happened to be "S2A_MSIL1C_20170129T110311_N0204_R094_T31UCT_20170129T110306.zip".

I unzipped the file and changed my working directory into

C:\data\S2A_MSIL1C_20170129T110311_N0204_R094_T31UCT_20170129T110306.SA
FE\GRANULE\L1C_T31UCT_A008385_20170129T110306\IMG_DATA>

I used this command:

C:\data\sentinel\S2A_MSIL1C_20170129T110311_N0204_R094_T31UCT_20170129T110306.SA
FE\GRANULE\L1C_T31UCT_A008385_20170129T110306\IMG_DATA>gdal_translate T31UCT_201
70129T110311_B08.jp2 c:\data\sentinel.tif
Input file size is 10980, 10980
0...10...20...30...40...50...60...70...80...90...100 - done.

Gdalinfo report from the resulting tiff file:

C:\data>gdalinfo sentinel.tif
Driver: GTiff/GeoTIFF
Files: sentinel.tif
Size is 10980, 10980
Coordinate System is:
PROJCS["WGS 84 / UTM zone 31N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",3],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH],
    AUTHORITY["EPSG","32631"]]
Origin = (300000.000000000000000,5800020.000000000000000)
Pixel Size = (10.000000000000000,-10.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  300000.000, 5800020.000) (  0d 3'56.87"E, 52d18'50.48"N)
Lower Left  (  300000.000, 5690220.000) (  0d 7'45.16"E, 51d19'40.99"N)
Upper Right (  409800.000, 5800020.000) (  1d40'33.33"E, 52d20'35.00"N)
Lower Right (  409800.000, 5690220.000) (  1d42'16.48"E, 51d21'21.89"N)
Center      (  354900.000, 5745120.000) (  0d53'37.96"E, 51d50'16.88"N)
Band 1 Block=10980x1 Type=UInt16, ColorInterp=Gray
  Image Structure Metadata:
    NBITS=15
0

Beside the method you have already found, using of="GTiff" might solve the problem

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.