[FOSSGIS-Talk] gdal_translate fuer GeoJPEG2000

[Markus Neteler]

Hallo,

2018-07-05 17:55 GMT+02:00 Jäger, Frank (KRZ) <F.Jaeger at krz.de>:
> Hallo,
> ich möchte in einem Bash-Script die DOPs vom Land NRW für einen WMS aufbereiten.
>
> Die 2014er hatte ich als RGBI-GeoTIFF bekommen, also 4 Kanäle, und so verarbeiten können.
> Die 2017er stehen nun im Format ".jp2" zum Download. Das Tool "gdalinfo" erzählt mir was über die eingebetteten Koordinaten (geo) und die 4 Kanäle (RGBI).
>
> Aber "gdal_translate" verschluckt sich dran, wenn es RGB daraus ziehen soll.
>
>   error: cannot decode code stream

das liegt am verwendeten Treiber. Ich verwende hier

gdalinfo --version
GDAL 2.2.4, released 2018/03/19

Bei Fedora ist dieser JP2 Treiber eingebunden (kann man sehen mit
"-debug on" flag)

gdalinfo --debug on 15_dop20rgbi_32377_5629_1_nw.jp2
GNM: GNMRegisterAllInternal
GNM: RegisterGNMFile
GNM: RegisterGNMdatabase
OPENJPEG: info: Start to read j2k main header (0).
OPENJPEG: info: Main header has been correctly decoded.
OPENJPEG: Grayscale color space
GDALJP2Metadata: Got projection from GeoJP2 (geotiff) box (0):
PROJCS["ETRS89 / UTM zone
32N",GEOGCS["ETRS89",DATUM["European_Terrestrial_Reference_System_1989",SPHEROID["GRS
1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6258"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4258"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",9],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","25832"]]
GDAL: GDALOpen(15_dop20rgbi_32377_5629_1_nw.jp2, this=0x560b85863480)
succeeds as JP2OpenJPEG.
Driver: JP2OpenJPEG/JPEG-2000 driver based on OpenJPEG library

^---- OpenJPEG library

Files: 15_dop20rgbi_32377_5629_1_nw.jp2
Size is 5000, 5000
Coordinate System is:
PROJCS["ETRS89 / UTM zone 32N",
    GEOGCS["ETRS89",
        DATUM["European_Terrestrial_Reference_System_1989",
            SPHEROID["GRS 1980",6378137,298.257222101,
                AUTHORITY["EPSG","7019"]],
[...]
    AXIS["Northing",NORTH],
    AUTHORITY["EPSG","25832"]]
Origin = (377000.000000000000000,5630000.000000000000000)
Pixel Size = (0.200000000000000,-0.200000000000000)
Image Structure Metadata:
  INTERLEAVE=PIXEL
OGRCT: PROJ >= 4.8.0 features enabled
OGRCT: Using locale-safe proj version
OGRCT: Source: +proj=utm +zone=32 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0
+units=m +no_defs
OGRCT: Target: +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs
Corner Coordinates:
Upper Left  (  377000.000, 5630000.000) (  7d15'15.19"E, 50d48'31.17"N)
Lower Left  (  377000.000, 5629000.000) (  7d15'16.40"E, 50d47'58.82"N)
Upper Right (  378000.000, 5630000.000) (  7d16' 6.26"E, 50d48'31.94"N)
Lower Right (  378000.000, 5629000.000) (  7d16' 7.46"E, 50d47'59.58"N)
Center      (  377500.000, 5629500.000) (  7d15'41.33"E, 50d48'15.38"N)
Band 1 Block=1024x1024 Type=Byte, ColorInterp=Undefined
  Overviews: 2500x2500, 1250x1250, 625x625, 312x312, 156x156, 78x78
  Overviews: arbitrary
  Image Structure Metadata:
    COMPRESSION=JPEG2000
Band 2 Block=1024x1024 Type=Byte, ColorInterp=Undefined
[...]


> Ich versuche nun unter Debian 8 die GDAL-Version 2.3.1 zu compilieren.

Es wird auch mit GDAL 2.2.x gehen (für GDAL 2.1 habe ich mal einen
Patch gemacht, siehe
https://bugzilla.redhat.com/attachment.cgi?id=1437926&action=diff)

> Entweder kann sie keine JP2 (configure .. --with-openjpeg)

Bei Fedora ist es openjpeg2-2.3.0-6.fc28.x86_64. Welcher Fehler tritt auf?

Schöne Grüße,
Markus


-- 
Markus Neteler, PhD
http://www.mundialis.de - free data with free software
http://grass.osgeo.org
http://courses.neteler.org/blog