Tuesday 19 June 2018

python - Installing jp2 driver for GDAL


I cannot open jp2-files with gdal:


This works fine (the Sentinel2 metadata):


~/copernicus/S2A_MSIL2A_20170111T100351_N0204_R122_T33UWP_20170111T100351.SAFE$ gdalinfo MTD_MSIL2A.xml

Driver: SENTINEL2/Sentinel 2
Files: MTD_MSIL2A.xml
Size is 512, 512
Coordinate System is `'

Metadata:
AOT_RETRIEVAL_ACCURACY=0.0
BARE_SOILS_PERCENTAGE=0.000388
CLOUD_COVERAGE_ASSESSMENT=0.1018
CLOUD_COVERAGE_PERCENTAGE=0.001909
CLOUD_SHADOW_PERCENTAGE=0.389902
DARK_FEATURES_PERCENTAGE=6.691170
[...]

But this not (one of the jp2 files):



~/copernicus/S2A_MSIL2A_20170111T100351_N0204_R122_T33UWP_20170111T100351.SAFE/GRANULE/L2A_T33UWP_A008127_20170111T100351/IMG_DATA/R10m$ gdalinfo L2A_T33UWP_20170111T100351_B02_10m.jp2
ERROR 4: `L2A_T33UWP_20170111T100351_B02_10m.jp2' not recognized as a supported file format.

gdalinfo failed - unable to open 'L2A_T33UWP_20170111T100351_B02_10m.jp2'.

I am using Anaconda 4.3.11 and created an environment with Python 2.7.13. I want to open jp2-Files which are part of the Sentinel L2A-Product (converted from L1C with sen2cor).


I've installed openjpeg 2.1.0 and gdal 2.1.0 via "conda install":


> conda list

# packages in environment at /home/roland/anaconda2/envs/copernicus:

#
curl 7.49.0 1
gdal 2.1.0 py27_0
geos 3.5.0 0
geotiff 1.4.1 0
hdf4 4.2.12 0
hdf5 1.8.17 1
jbig 2.1 0
jpeg 8d 2
kealib 1.4.6 0

libgdal 2.1.0 0
libnetcdf 4.4.1 0
libpng 1.6.27 0
libtiff 4.0.6 2
mkl 2017.0.1 0
numpy 1.12.0 py27_0
openjpeg 2.1.0 5 sunpy
openssl 1.0.2k 0
pip 9.0.1 py27_1
proj4 4.9.2 0

python 2.7.13 0
readline 6.2 2
setuptools 27.2.0 py27_0
sqlite 3.13.0 0
tk 8.5.18 0
wheel 0.29.0 py27_0
xerces-c 3.1.4 0
xz 5.2.2 1
zlib 1.2.8 3


The main problem seems to be that I don't manage to install the JPEG2000-driver:


> gdal-config --formats
gxf gtiff hfa aigrid aaigrid ceos ceos2 iso8211 xpm sdts raw dted mem jdem envisat elas fit vrt usgsdem l1b nitf bmp airsar rs2 ilwis rmf leveller sgi srtmhgt idrisi gsg ingr ers jaxapalsar dimap gff cosar pds adrg coasp tsx terragen blx msgn til r northwood saga xyz hf2 kmlsuperoverlay ctg e00grid zmap ngsgeoid iris map cals safe sentinel2 mrf wcs wms plmosaic wmts grib bsb netcdf kea hdf5 hdf4 gif jpeg png pcraster pcidsk rik ozi pdf arg

How can I do this?



Answer



It looks like this is a problem with the conda packaging of gdal. I see that the current gdal version on conda is 2.1.3, not 2.1.0, so hopefully this packaged version works. Consider removing the current gdal and reinstalling.


I looked at the contents of the gdal conda package I mentioned and for 64-bit linux python 2.7 numpy 1.12 the files gdaljp2abstractdataset.h, gdaljp2metadata.h, and gdaljp2metadatagenerator.h are included. Haven't installed it to confirm whether the driver exists.


No comments:

Post a Comment

arcpy - Changing output name when exporting data driven pages to JPG?

Is there a way to save the output JPG, changing the output file name to the page name, instead of page number? I mean changing the script fo...