[FOSSGIS-Talk] gdal_translate fuer GeoJPEG2000
Markus Neteler
neteler at osgeo.org
Fr Jul 6 10:47:54 CEST 2018
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
Mehr Informationen über die Mailingliste FOSSGIS-Talk-Liste