Karte mit map Datei laden

Hallo,

GDAL kann ja Ozi map Dateien in gewissen Grenzen lesen. Ich schaffe es aber nicht, sowas in QLandkarte zu laden, weder direkt noch über F1, F6. Wird das geblockt oder wie muß man da vorgehen?

Was ich in Erinnerung habe, kann gdal eine map+ozfx öffnen. Man muß dann die *.map angeben und die ozfx muß in der map eingetragen sein. Weiters kann gdal gif+map oder jpg+map öffnen, dann muß man aber die *.gif angeben und die *.map muß den gleichen Namen haben. Das liegt daran, daß das von verschiedenen Leuten völlig unkoordiniert da reingehackt wurde. tif+map geht merkwürdigerweise nicht.

Ich hab jetzt mal folgendes ausprobiert: gif+map (russische Genstabkarte) mit gleichem Namen: gdalinfo xy.gif liefert was vernünftiges. gdalinfo xy.map liefert not supported file format. QLandkarte Karte laden xy.gif liefert "keine Georefernzinformationen gefunden"
F1 F6 hinzufügen xy.gif liefert "Laden der Datei xy.gif fehlgeschlagen"

Nachtrag: für map+ozfx wurde ein eigener Treiber geschrieben, das ist also die "ordentliche" Variante, und nachdem sie von Des Newman die Erlaubnis bekommen haben, die Ozi Datumsdefinitionen zu verwenden, wird das wohl da eingebaut werden.

Für gif+map wurde einfach nur der gif Treiber gehackt und deshalb funktioniert das wohl nur mit einigen ausgewählten Datums und Projektionen und wird vermutlich nicht weiter ausgebaut werden.

Das liefert das gleiche wie mit xy.gif. Das wird auch gar nicht erst in die Kartensammlung übernommen - weder in eine neue noch in eine bestehende.

Dieses F1 F6 ist mir eh suspekt, weil da dann steht "Erzeuge Kartensammlung aus Geotiff" Das sind aber keine Geotiffs.

Ja, eine ozf kann man nicht direkt laden, sondern nur die map, in welcher dann in der 2. oder 3.Zeile der Name der ozf stehen muß. Der Name kann dann auch anders sein. Für die ozf ist der Ozi Treiber zuständig, das ist ein völlig anderer Teil der Software.

Mir ist es jetzt auch gelungen eine ozf2 zu laden - einfach über Datei|Karte laden und dann die map anklicken. Nicht geklappt hat es aber mit einer ozfx3. Da kommt die Meldung "keine Georeferenzinformationen gefunden", obwohl gdalinfo diese ausspuckt und gdal_translate sogar in der Lage ist, das in ein tif umzuwandeln.

Nachtrag der Vollständigkeit halber: hier steht das mit den Ozi Datums:
http://trac.osgeo.org/gdal/ticket/3929

Das neue GDAL 1.10 bietet jetzt erweitere Möglichkeiten mit .map Dateien:

http://trac.osgeo.org/gdal/wiki/Release/1.10.0-News

Zitat

ozi2vrt.py a script to convert .MAP file into .VRT
ozi2gdal.py an enhanced version of the script to convert .MAP file into .VRT
gdal-ozimap.patch.gz OziExplorer? .MAP driver

sowie die original Oziexplorer Datumsdefinitionen:

Zitat

data/ozi_datum.csv
data/ozi_ellips.csv

In GDAL 1.10 schreibt man jetzt immer die *.map in die Befehlszeile. Damit ist das endlich vereinheitlicht und funktionioniert auch mit anderen als ozf Dateien. Der Treiber heißt map driver.

Ich habe nun einen Grund gefunden, warum QLandkarte GT manche dieser Ozidateien nicht akzeptiert. Das hat mit Ozi erstmal garnichts zu tun, sondern damit, daß es 2 verschiedene Arten gibt, eine Georeferenzierung in eine GeoTiff zu schreiben. Bei der ersten Variante werden Koordinaten im Stil eines Wordfiles eingetragen. Das können die meisten Programme inklusive QLGT lesen.

Wenn man jedoch von einer map Datei ausgeht, dann trägt gdal_translate einfach die GCPs aus der *.map in die Tiff Datei ein. Das ist offenbar auch eine zulässige Variante, wird aber von QLGT nicht akzeptiert: keine Georeferenzinformationen gefunden. TTQV4 und QV6 interpretieren das auch nicht richtig. So eine GeoTiff hab ich mal angehängt; in gdalinfo sieht das dann so aus:

Zitat

C:\maps\50000>gdalinfo -noct test.tif
Driver: GTiff/GeoTIFF
Files: test.tif
Size is 3403, 3117
Coordinate System is `'
GCP Projection =
PROJCS['unnamed',
GEOGCS['Pulkovo 1942',
DATUM['Pulkovo_1942',
SPHEROID['Krassowsky 1940',6378245,298.2999999999985,
AUTHORITY['EPSG','7024']],
AUTHORITY['EPSG','6284']],
PRIMEM['Greenwich',0],
UNIT['degree',0.0174532925199433],
AUTHORITY['EPSG','4284']],
PROJECTION['Transverse_Mercator'],
PARAMETER['latitude_of_origin',0],
PARAMETER['central_meridian',45],
PARAMETER['scale_factor',1],
PARAMETER['false_easting',500000],
PARAMETER['false_northing',0],
UNIT['metre',1,
AUTHORITY['EPSG','9001']]]
GCP[ 0]: Id=1, Info=
(100,106) -> (357685.766463653,4782807.11213897,0)
GCP[ 1]: Id=2, Info=
(1647,100) -> (367851.254555488,4782602.22720938,0)
GCP[ 2]: Id=3, Info=
(3297,106) -> (378016.70146789,4782412.52758289,0)
GCP[ 3]: Id=4, Info=
(98,1564) -> (357492.398996449,4773548.73195437,0)
GCP[ 4]: Id=5, Info=
(1698,1564) -> (367671.707560777,4773343.8868211,0)
GCP[ 5]: Id=6, Info=
(3300,1561) -> (377850.973054314,4773154.22407581,0)
GCP[ 6]: Id=7, Info=
(98,3019) -> (357299.334569795,4764290.47450573,0)
GCP[ 7]: Id=8, Info=
(1647,3017) -> (367492.44198133,4764085.67090163,0)
GCP[ 8]: Id=9, Info=
(3301,3020) -> (377685.504426164,4763896.04664196,0)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 3117.0)
Upper Right ( 3403.0, 0.0)
Lower Right ( 3403.0, 3117.0)
Center ( 1701.5, 1558.5)
Band 1 Block=3403x2 Type=Byte, ColorInterp=Palette
Metadata:
GIF_BACKGROUND=0
Color Table (RGB with 256 entries)

Alles anzeigen

Globalmapper liest sowas richtig, QLGT aber nicht! Erst wenn man das dann nochmal durch gdalwarp wurschtelt, dann wird die Referenzierung im Stil eines Worldfiles eingetragen und dann liest das auch QLGT. Ich vermute mal, daß gdal_translate nicht die Fähigkeit hat, sowas umzurechnen.

Hallo Martin,

dieses Skript scheint mir überflüssig. GDAL kann die GCPs aus der .map wunderbar extrahieren und damit ein GeoTiff erzeugen. Nur QLGT kann das dann nicht lesen. Das ist meiner Meinung nach der eigentliche Knackpunkt bei der Geschichte.

Zitat von teleskopix

...und wenn qlgt zu einer Bilddatei eine *.gcp hat, macht qlgt im Handumdrehen daraus eine Geotiff

Besteht das Handumdrehen dann in einer Reprojektion mit gdalwarp? Dann kann man alles allein mit GDAL machen.

Ich hab mich noch ein bischen mit GDAL beschäftigt:

die Befehlszeile zur Umwandlung lautet:

gdal_translate -co COPY_SRC_OVERVIEWS=YES -co TILED=YES -co COMPRESS=LZW name.map name.tif

In der name.map muß die zugehörige Bilddatei angegeben sein. Das kann jetzt im Prinzip alles sein, was GDAL verarbeiten kann, insbesondere kann es auch eine tif Datei sein. Man kann jetzt also map+tif in GeoTiff umwandeln.

Die Projektionen UTM und Lambert Conformal Conic funktionieren.

Mit "-co COPY_SRC_OVERVIEWS=YES" übernimmt man die in der ozf Datei enthaltenen Overviews direkt in die GeoTiff. gdaladdo ist dann nicht mehr nötig.

Im gdal-data Verzeichnis befinden sich eine ozi_datum.csv und eine ozi_ellips.csv, die man manuell editieren und so weitere Datums hinzufügen kann.

Es gibt eine Umgebungsvariable OZI_APPROX_GEOTRANSFORM. Wenn man die auf YES setzt, dann zwingt man GDAL, die Geotransformation durchzuführen.

Und noch was sehr interessantes, was aber nichts mit Ozi zu tun hat: man kann jetzt JPEG Bilder ohne Rekompression - also verlustfrei - in Tiff mit JPEG Kompression umwandeln!

Zitat von Even Rouault

Note: in GDAL 1.10, I've introduced in the GeoTIFF driver an optimization that
will take the JPEG coefficients of a JPEG dataset and put them as a JPEG-in-TIFF
without doing any decompression/decompression, so preserving exactly the
quality of the source JPEG (and being much faster)

This is accomplished by : gdal_translate my.jpeg my.tif -co COMPRESS=JPEG
(if you add -co QUALITY=xxx , or some other options, then the optimization
will not be triggered)

Note that it currently only works with a single JPEG image, not a VRT of
several JPEG for example.

Alles anzeigen

http://osgeo-org.1560.x6.nabbl…o-geotiffs-td5039714.html