Topo Karten konvertieren

Hallo Thomas,

es schaut so aus als hätten Oliver und ich das Problem bei der Erstellung der GPX Dateien unter WinXP gefunden.

Eine Testversion gibt es zum Download:
http://sourceforge.net/project…_2637_export.exe/download

Bitte teste es und gib uns bescheid, ob es jetzt funktioniert.

Danke,

Helmut

Hallo Oliver und Helmut,

erst einmal vielen Dank für die sehr prompte Umsetzung des patches! Die Erzeugung der Custom Maps funktioniert jetzt. Damit ist die tool chain wieder ein wichtiges Stück weiter. Da ich hier auf meinem Schreibtisch kein Garmin habe, das custom maps beherrscht, muss ich jetzt die erzeugten .kmz files erst mal hohbrugg schicken, damit er sie in sein Gerät laden kann. Ich werd' ihm dann noch ein paar Wegpunkte als .gpx mitschicken, dann können wir gleich sehen, ob die Punkte an der richtigen Stelle auf der Karte liegen. Das wird wahrscheinlich im Lauf des Wochenendes laufen, ich werde Euch auf alle Fälle Rückmeldung geben. Wenn das alles klappt, müsste ich dann nur noch ein script file erstellen, damit wir den Kartenstapel von hohbrugg dann gut verarbeitet bekommen (und vor allem keine Tipfehler die Konvertierung stören...). Aber das ist vertrautes Gebiet.

Was die verzerrte Karte angeht: ich habe noch einmal versucht, die .gif Karte zu laden und zu referenzieren, aber es ist mir nicht mehr gelungen, eine verzerrte geotiff Karte zu erzeugen. Daher vermute ich stark, dass ich wohl bei den Einstellungen für die Referenzierung irgendwo was Falsches oder Unpassendes ausgewählt oder eingegeben habe, was ich jetzt nicht mehr rekonstruiert bekomme. Die beim ersten Versuch entstandene geotiff Karte ist dann wohl nicht nur verzerrt, sondern dürfte auch in der Referenzierung korrupte Daten enthalten, die dann in der Folge QLandkarteGT daran hindern, die Karte in eine Kartensammlung einzubauen. Man sollte halt Experimente immer verifizieren, bevor man darüber berichtet...

Also nochmal vielen Dank für Eure Bemühungen, ich melde mich wieder, wenn ich weiss wie es weitergegangen ist,

Gruß,
Thomas

Hallo GPS Cracks
herzlichen Dank dafür, dass ihr Euch so reinhängt. Ist ja echt fantastisch, was ihr in der Zeit auf die Reihe bekommen habt. Werde nächste Woche sobald ich wieder im Land bin, die Daten testen.
Thomas auch dir herzlichen Dank, für Deine Arbeit.
Gruss
Hohbrugg

Hallo Oliver!

Das Problem mit dem Kartenrand ist mir aufgefallen, als ich zum ersten Mal zwei benachbarte Karten nach QLandkarteGT geladen habe. Nun haben wir ja für jeden Karten-gif-File eine map-File, so dass ich nicht einfach den gif-File kleiner schneiden kann. Ich müsste dann auch die Referenzierungen im map-File korrigieren (da es sich um einige Files handelt, müsste ich dafür wohl ein script schreiben, was zwar machbar, aber nicht ganz einfach ist). Ausserdem ist der Rand bei verschiedenen Karten geringfügig unterschiedlich breit, und esscheint so, als ob die Karten auch minimal verdreht sind, jedenfalls liegen die Referenpunkte im map-File nicht exakt auf dem x/y-Raster des gif-Files.

Zitat

Ah, ich sehe schon. hohbrugg hat definitiv auf das richtige Pferd gesetzt [Blockierte Grafik: http://1.1.1.3/bmi/www.naviboa…b/images/smilies/wink.gif] Was sicherlich noch interessant wäre, wenn man die einzelnen Karten zu einer Großen zusammen fügt. Ich kenne jetzt eure Quelldateien nicht, aber poehalis hat immer diesen lästigen Rand. Wenn man den weg lässt und die Karten sich dann noch überschneiden, dann kann man mit:

gdalwarp file1.tif file2.tif (usw.) file_out.tif

alle zusammenfügen. Den Rand entfernen kannst Du entweder in QLGT mit "Karte exportieren -> QLandkarte M" machen (öde) oder, wenn der Rand immer gleich dick ist, mit gdal_translate und einem Skript.

Wenn ich Dich richtig verstehe, baut gdalwarp in Deinen Hinweis auch 'nur' mehrere tif Files zu einem größeren zusammen, wenn die sich überlagern, wird von einem der beiden Files (dem Bild das 'oben' liegt) trotzdem der Rand dargestellt.
Was für unseren Fall eher passen könnte, wäre Folgendes: kann gdal_translate ein referenziertes geotif zuschneiden, ohne dass die Referenzierung verloren geht bzw. ungültig wird? Dann könnte ich die geotif einfach kleiner schneiden, es würde dann zwar ein kleiner Sicherheitsrand stehen bleiben, aber das wäre schon mal wesentlich besser. Ist es das, was Du mit Deinem Hinweisauf gdal_translate meinst?

Übrigens: wo finde ich eine Art Bedienungsanleitung für die gdal-Programme? In http://www.gdal.org/gdal_utilities.html sind zwar die Aufrufe der Programme und einige Parameter kurz beschrieben, aber es ist schwer, damit herauszufinden, was man wo einstellen muss bzw. kann, damit ein Programm das Gewünschte auch tut.

Gruß,
Thomas

Hallo Oliver,

dann werd' ich da mal etwas mit gdal_translate spielen...

Übrigens: warum wird eigentlich der EPSG-Code 3785 benutzt? Laut http://www.epsg-registry.org/ gehört dieser Code zur Region 'New Zealand - South Island - Marlborough meridional circuit area'. Das ist ja nicht gerade Mitteleuropa oder auch nur Nordhalbkugel... Sicher eine theoretische Frage, denn offensichtlich funktioniert die tool chain ja, aber interessant wär's schon.

Gruß,
Thomas

Zitat von thomas_rolf

Übrigens: warum wird eigentlich der EPSG-Code 3785 benutzt?

Die EPSG hat da wohl zwischendurch mal die Nummer geändert.
D.h. um sicher zu gehen müsste man ggf beser 3857 nehmen.
Siehe http://viswaug.wordpress.com/2…ne-crazy-epsg-especially/ oder
http://www.maptiler.org/google…s-tile-bounds-projection/

Bei gdal/fwtools ist offenbar noch die "deprecated" 3875 im Angebot, nicht die neue (New Zealand).

>testepsg EPSG:3785
Validate Succeeds.
WKT[EPSG:3785] =
PROJCS['Popular Visualisation CRS / Mercator (deprecated)',
    GEOGCS['Popular Visualisation CRS',
        DATUM['Popular_Visualisation_Datum',
            SPHEROID['Popular Visualisation Sphere',6378137,0,
                AUTHORITY['EPSG','7059']],
            TOWGS84[0,0,0,0,0,0,0],
            AUTHORITY['EPSG','6055']],
        PRIMEM['Greenwich',0,
            AUTHORITY['EPSG','8901']],
        UNIT['degree',0.01745329251994328,
            AUTHORITY['EPSG','9122']],
        AUTHORITY['EPSG','4055']],
    UNIT['metre',1,
        AUTHORITY['EPSG','9001']],
    PROJECTION['Mercator_1SP'],
    PARAMETER['central_meridian',0],
    PARAMETER['scale_factor',1],
    PARAMETER['false_easting',0],
    PARAMETER['false_northing',0],
    EXTENSION['PROJ4','+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +
x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs'],
    AUTHORITY['EPSG','3785'],
    AXIS['X',EAST],
    AXIS['Y',NORTH]]


Simplified WKT[EPSG:3785] =
PROJCS['Popular Visualisation CRS / Mercator (deprecated)',
    GEOGCS['Popular Visualisation CRS',
        DATUM['Popular_Visualisation_Datum',
            SPHEROID['Popular Visualisation Sphere',6378137,0],
            TOWGS84[0,0,0,0,0,0,0]],
        PRIMEM['Greenwich',0],
        UNIT['degree',0.01745329251994328]],
    UNIT['metre',1],
    PROJECTION['Mercator_1SP'],
    PARAMETER['central_meridian',0],
    PARAMETER['scale_factor',1],
    PARAMETER['false_easting',0],
    PARAMETER['false_northing',0]]


Old Style WKT[EPSG:3785] = PROJCS['Popular Visualisation CRS / Mercator (depreca
ted)',GEOGCS['Popular Visualisation CRS',DATUM['Popular_Visualisation_Datum',SPH
EROID['Popular Visualisation Sphere',6378137,0]],PRIMEM['Greenwich',0],UNIT['deg
ree',0.01745329251994328]],UNIT['metre',1],PROJECTION['Mercator_1SP'],PARAMETER[
'central_meridian',0],PARAMETER['scale_factor',1],PARAMETER['false_easting',0],P
ARAMETER['false_northing',0]]
ESRI'ified WKT[EPSG:3785] =
PROJCS['Popular_Visualisation_CRS_Mercator_deprecated',
    GEOGCS['GCS_Popular Visualisation CRS',
        DATUM['D_',
            SPHEROID['Popular_Visualisation_Sphere',6378137,0]],
        PRIMEM['Greenwich',0],
        UNIT['Degree',0.017453292519943295]],
    PROJECTION['Mercator'],
    PARAMETER['central_meridian',0],
    PARAMETER['false_easting',0],
    PARAMETER['false_northing',0],
    UNIT['Meter',1],
    PARAMETER['standard_parallel_1',0.0]]
PROJ.4 rendering of [EPSG:3785] = +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +
lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs




>testepsg EPSG:3857
Validate Succeeds.
WKT[EPSG:3857] =
PROJCS['WGS 84 / Pseudo-Mercator',
    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.01745329251994328,
            AUTHORITY['EPSG','9122']],
        AUTHORITY['EPSG','4326']],
    UNIT['metre',1,
        AUTHORITY['EPSG','9001']],
    PROJECTION['Mercator_1SP'],
    PARAMETER['central_meridian',0],
    PARAMETER['scale_factor',1],
    PARAMETER['false_easting',0],
    PARAMETER['false_northing',0],
    EXTENSION['PROJ4','+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +
x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs'],
    AUTHORITY['EPSG','3857'],
    AXIS['X',EAST],
    AXIS['Y',NORTH]]


Simplified WKT[EPSG:3857] =
PROJCS['WGS 84 / Pseudo-Mercator',
    GEOGCS['WGS 84',
        DATUM['WGS_1984',
            SPHEROID['WGS 84',6378137,298.257223563]],
        PRIMEM['Greenwich',0],
        UNIT['degree',0.01745329251994328]],
    UNIT['metre',1],
    PROJECTION['Mercator_1SP'],
    PARAMETER['central_meridian',0],
    PARAMETER['scale_factor',1],
    PARAMETER['false_easting',0],
    PARAMETER['false_northing',0]]


Old Style WKT[EPSG:3857] = PROJCS['WGS 84 / Pseudo-Mercator',GEOGCS['WGS 84',DAT
UM['WGS_1984',SPHEROID['WGS 84',6378137,298.257223563]],PRIMEM['Greenwich',0],UN
IT['degree',0.01745329251994328]],UNIT['metre',1],PROJECTION['Mercator_1SP'],PAR
AMETER['central_meridian',0],PARAMETER['scale_factor',1],PARAMETER['false_eastin
g',0],PARAMETER['false_northing',0]]
ESRI'ified WKT[EPSG:3857] =
PROJCS['WGS_84_Pseudo_Mercator',
    GEOGCS['GCS_WGS_1984',
        DATUM['D_WGS_1984',
            SPHEROID['WGS_1984',6378137,298.257223563]],
        PRIMEM['Greenwich',0],
        UNIT['Degree',0.017453292519943295]],
    PROJECTION['Mercator'],
    PARAMETER['central_meridian',0],
    PARAMETER['false_easting',0],
    PARAMETER['false_northing',0],
    UNIT['Meter',1],
    PARAMETER['standard_parallel_1',0.0]]
PROJ.4 rendering of [EPSG:3857] = +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +
lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs

Hallo zusammen,

@Helmut:vielen Dank für die Info zu EPSG. Dass die Nummern mal eben um die halbe Welt umziehen können, findet sich auf der EPSG-Seite natürlich nicht so direkt...

@oliver: nachdem ich die Test-kmz-Kacheln an Hohbrugg geschickt hatte, bekam ich heute Antwort: es klappt teilweise. Er kann die Karten wohl sehen, aber sie sind an der falschen Stelle und verzerrt. Man kann das auch in Google Earth sehen, wenn man die kmz-Files dort lädt. Dieser Fehler liegt wahrscheinlich daran, daß ich den Kartenausschnitt etwas speziell gewählt habe (och hab' hohbrugg inzwischen neue Files geschickt, hab' aber noch keine Rückmeldung).
Der Fehler tritt nur auf, wenn man in QLandkarteGT zwei oder mehr Karten (die sich in meinem Fall leicht überlappen, wegen des Kartenrandes) in einer Kartensammlung zusammengefasst hat und dann einen Kartenausschnitt zum Exportieren über den Rand einer der Teilkarten hinweg definiert. QLandkarte GT exportiert den Ausschnitt, legt dabei aber keine Kacheln an, die Teile aus beiden benachbarten Teilkarten enthalten, sondern erzeugt am Übergang kleinere Kacheln, die jeweils am Teilkartenrand enden. Wenn ich die Teilkacheln dann in Google Earth lade, liegen diese nicht mehr nebeneinander, sondern alle übereinander, ausserdem sind sie in der x-Achse verzerrt.
Wenn ich den Ausschnitt nur innerhalb einer Teilkarte wähle, werden die Kacheln ordentlich erzeugt, jedenfalls, sieht in Google Earth alles ordentlich aus.
Wenn ich die Teilkarten mit 'gdalwarp file1.tif file2.tif (usw.) file_out.tif' vorher zu einer Karte zusammenbaue, tritt der Fehler nicht auf, für QLandkarteGT ist das ja dann eine Karte. Beim Kopieren wurde bei mir allerdings nur der erste tif-File korrekt nach file_out übertragen, beim zweiten File war die Auflösung irgendwie deutlich geringer als vorher, ausserdem war die Farb-Map korrupt. Muss ich da noch weitere Parameter setzen?

Gruß,
Thomas

Hallo zusammen,

bin gerade am Ausprobieren, was man mit gdal_translate beim Zuschneiden der Karten so alles anstellen kann. Weiß jemand, in welchem Format beim gdal_translate-Parameter -projwin die Eckpunkte angegeben werden müssen (ulx, uly, lrx, lry)? Vorliegen habe ich die Lat/Lon-Werte direkt von der Karte (sie stehen ja auch im map-File), aber wie muss ich die angeben? Einfach die Dezimalwerte nehmen klappt nicht, gdal_translate vermutet den gesuchten Ausschnitt weit ausserhalb der Karte.

Gruß,
Thomas

Hallo Oliver,

srcwin hab'ich schon ausprobiert, würde schon reichen... Aber projwin wäre eleganter, und ausserdem könnte ich da die Werte aus dem map-File entnehmen (im Idealfall innerhalb eines scriptes). Ich bin halt im Moment am Rumspielen und Testen...

Hier die Antwort von gdalinfo:

Driver: GTiff/GeoTIFF
Files: 200k--m45-28.tif
Size is 3067, 3089
Coordinate System is:
PROJCS['Popular Visualisation CRS / Mercator',
    GEOGCS['Popular Visualisation CRS',
        DATUM['Popular_Visualisation_Datum',
            SPHEROID['Popular Visualisation Sphere',6378137,0,
                AUTHORITY['EPSG','7059']],
            AUTHORITY['EPSG','6055']],
        PRIMEM['Greenwich',0],
        UNIT['degree',0.0174532925199433],
        AUTHORITY['EPSG','4055']],
    PROJECTION['Mercator_1SP'],
    PARAMETER['central_meridian',0],
    PARAMETER['scale_factor',1],
    PARAMETER['false_easting',0],
    PARAMETER['false_northing',0],
    AUTHORITY['EPSG','3785'],
    EXTENSION['PROJ4','+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs'],
    UNIT['metre',1,
        AUTHORITY['EPSG','9001']]]
Origin = (9680403.680058261400000,6335385.307206181800000)
Pixel Size = (39.000290546521789,-39.000290546521789)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  ( 9680403.680, 6335385.307) ( 86d57'37.96"E, 49d 9'54.43"N)
Lower Left  ( 9680403.680, 6214913.410) ( 86d57'37.96"E, 48d27'16.18"N)
Upper Right ( 9800017.571, 6335385.307) ( 88d 2'6.20"E, 49d 9'54.43"N)
Lower Right ( 9800017.571, 6214913.410) ( 88d 2'6.20"E, 48d27'16.18"N)
Center      ( 9740210.626, 6275149.358) ( 87d29'52.08"E, 48d48'39.81"N)
Band 1 Block=3067x2 Type=Byte, ColorInterp=Palette
  Color Table (RGB with 256 entries)
    0: 172,238,254,255
    1: 240,149,116,255

(Color Table gekürzt)

Gruß,
Thomas

Hallo zusammen!

Ich war eine Weile still, aber das lag zum Einen daran, daß ich mit den tools rumexperimentiert habe, zum Anderen an beruflichen Verpflichtungen, die mich auf mehrere Dienstreisen führten...

Der Stand der Dinge: Ich habe eine elegante Methode gefunden, wie ich meine Karten automatisch randbeschneiden kann, so daß ich sie dann anschliessend randfrei in QLandkarteGT darstellen kann, aber auch mit gdal_translate zu großen randfreien Karten zusammenbauen kann. Ein automatischer Ablauf über scripte war nötig, da es sich doch um einige Karten handelt, die bearbeitet werden wollen (und ausserdem war da auch noch etwas persönlicher Ehrgeiz...). Damit auch andere etwas davon haben, hier meine scripte: (Ausgangslage *.gif maps mit *.map ozi-explorer Referenzdaten, die Einzelkarten fangen jeweils mit 100k, 200k oder 500k an, je nach Maßstab der Karte, die Kartenränder folgen den Längen- und Breitengraden).
zunächst convert.bat (die zwischenzeitliche Verdoppelung der Pixel-Zahl sorgt dafür, das sich die Bildqualität nach zweimaligem gdalwarp nicht zu stark verschlechtert (gibt es da noch eie andere Möglichkeit?)):

REM wandelt alle *.gif (mit zugehörigem *.map) jeweils in ein komprimiertes, gekacheltes und randbeschnittenes *.tif geotiff um


REM Achtung: file muss in der directory gestartet werden, in denr die *.gif und *.map liegen!!


REM Einstellungen und Pfade laden
call "C:\Programme\FWTools2.4.7\setfw.bat"


echo on


REM Beginn der Schleife über all Files
for %%i in (*.gif) do (


   for /F "tokens=1 delims=." %%j in ("%%i") do (


REM Erhöhung Pixel-Zahl
   echo start gdal_translate -outsize 200%% 200%% %%j.gif %%j.tifi
   gdal_translate -outsize 200%% 200%% %%j.gif %%j.tifi
   echo ende gdal_translate


REM Umwandlung zunächst nach Lon/Lat, WGS84
   echo start gdalwarp -t_srs "EPSG:4326" %%j.tifi %%j.1.tifi
   gdalwarp -t_srs "EPSG:4326" %%j.tifi %%j.1.tifi
   echo ende gdalwarp


REM Erzeugen info-File *.txt
   echo start gdalinfo %%j.1.tifi > %%j.1.txt
   gdalinfo %%j.1.tifi > %%j.1.txt
   echo ende gdalinfo


REM Extrahieren randbeschnittene Karte
   echo start call "../scripte/extract.pl" %%j.1
   call "../scripte/extract.pl" %%j.1
   echo end extract


REM Umwandlung nach UTM, WGS84
   echo start gdalwarp -t_srs "EPSG:3785" %%j.1.1.tifi %%j.2.tifi
   gdalwarp -t_srs "EPSG:3785" %%j.1.1.tifi %%j.2.tifi
   echo ende gdalwarp


REM Komprimieren, Kacheln und Pixel-Zahl reduzieren
   echo start gdal_translate -outsize 50%% 50%% -co tiled=yes -co compress=LZW %%j.2.tifi %%j.tif
   gdal_translate -outsize 50%% 50%% -co tiled=yes -co compress=LZW %%j.2.tifi %%j.tif
   echo ende gdal_translate
   )
)
pause

Da die Möglichkeiten, innerhalb eines windows command files Berechnungen durchzuführen, doch sehr beschränkt sind, habe ich das Zuschneiden der Karten in ein perl-script ausgelagert, extract.pl (siehe Aufruf oben und Inhalt unten)

use strict;
use warnings;
my $file;
my $in_arg = shift @ARGV;
if ( ! defined $in_arg ) {            # filename holen
    print "Filename (ohne extension):";
    $file = <STDIN>;
    chomp $file;     }
else { $file = $in_arg;  }
my $read_file = "$file.txt";            # oeffnen output von gdalinfo
if ( -e $read_file ) {
    open FILE, '<', $read_file or die "Fehler beim Öffnen von $read_file zum Lesen: $!";
    my $x_min_1 = 99999;
    my $y_min_1 = 99999;
    my $x_max_1 = 0;
    my $y_max_1 = 0;
    my $found_1 = 0;
    my $x_fact  = 1;
    my $y_fact  = 1;
    if ( $file =~ m/100k/ )    {         # festlegen der Rundungsparameter je nach Kartenformat
        $x_fact = 4;            # 1:100k: x-Achse: 15'-Streifen
        $y_fact = 6;       }        #         y-Achse: 10'-Streifen
    elsif ( $file =~ m/200k/ ) {     
        $x_fact = 2;            # 1:200k: x-Achse: 30'-Streifen
        $y_fact = 3;       }        #         y-Achse: 20'-Streifen
    elsif ( $file =~ m/500k/ ) { 
        $x_fact = 2;            # 1:500k: x-Achse: 30'-Streifen (eigentlich 1.5°) 
        $y_fact = 1;       }        #         y-Achse: 1°-Streifen
    else                   { 
        $x_fact = 1;
        $y_fact = 1;       }
    while ( <FILE> )                                                 {
        chomp $_;                # Formatieren und Zerlegen der Zeilen gdalinfo-File
        $_ =~ tr/(/,/;
        $_ =~ tr/)/,/;
        my @array = split ( ",", $_, 4 );
        if ( $array[0] =~ m/Upper Left|Lower Left|Upper Right|Lower Right/ ) {
            $found_1 = 1;
            print "$_\n";
            if ( $array[1] < $x_min_1 ) { $x_min_1 = $array[1] }
            if ( $array[1] > $x_max_1 ) { $x_max_1 = $array[1] }
            if ( $array[2] < $y_min_1 ) { $y_min_1 = $array[2] }
            if ( $array[2] > $y_max_1 ) { $y_max_1 = $array[2] }           }        }
    if ( $found_1 == 1 )        {
        $x_min_1 = $x_min_1 + 0;
        $y_min_1 = $y_min_1 + 0;
        $x_max_1 = $x_max_1 + 0;
        $y_max_1 = $y_max_1 + 0;
        print "Corner Parameter: x_min=$x_min_1 x_max=$x_max_1 y_min=$y_min_1 y_max=$y_max_1 \n"; 
        $x_min_1 = $x_min_1 * $x_fact;        # Runden der gelesenen Corner-Parameter entsprechend Vorgabe
        $x_min_1 = sprintf ("%.0f", $x_min_1); 
        $x_min_1 = $x_min_1 / $x_fact;
        $y_min_1 = $y_min_1 * $y_fact;
        $y_min_1 = sprintf ("%.0f", $y_min_1); 
        $y_min_1 = $y_min_1 / $y_fact;
        $x_max_1 = $x_max_1 * $x_fact;
        $x_max_1 = sprintf ("%.0f", $x_max_1); 
        $x_max_1 = $x_max_1 / $x_fact;
        $y_max_1 = $y_max_1 * $y_fact;
        $y_max_1 = sprintf ("%.0f", $y_max_1); 
        $y_max_1 = $y_max_1 / $y_fact;
        print "korrigierte Corner Parameter: x_min=$x_min_1 x_max=$x_max_1 y_min=$y_min_1 y_max=$y_max_1 \n"; }
    else
        { print "Zeilen mit Corner Parametern wurden nicht gefunden!\n" }
    close FILE;
    print "gdal_translate -projwin $x_min_1 $y_max_1 $x_max_1 $y_min_1 $file.tifi $file.1.tifi\n";
    open FILE, "gdal_translate -projwin $x_min_1 $y_max_1 $x_max_1 $y_min_1 $file.tifi $file.1.tifi |";
    print <FILE>;
    close FILE or die "Fehler bei Ausführung von gdal_translate\n";    }
else    { print "$file.txt existiert nicht!\n" }

Anschliessend kann ich in QLandkarteGT aus den Karten problemlos die *.kmz Kacheln erzeugen.
Uns bleibt jetzt noch ein letztes Problem zu lösen: hohbruck hat ein neues GPSmap62 von Garmin, er kann die Files auch in den richtigen Ordner im Gerät laden, er kann sie auch im Menupunkt '"Karte einrichten/ Karten-informationen Karte Wählen/Benutzer-definierte Karten/aktivieren" aktivieren (zumindest entsprechend markieren und die Aktion starten), aber er schafft es nicht, die Karte dann auch tatsächlich zu sehen. Irgend etwas klemmt da noch...
Hat einer von Euch eine Idee, oder gibt es einen anderen thread, der sich damit/mit solchen Themen beschäftigt?
Gruß,
Thomas

Hallo Thomas,

ein Beispiel-KMZ würde die Frage nach der fehlenden Anzeige im Map62 einfach klären lassen . . .

Grüsse - Anton

Hallo Anton und Oliver!

Ich hatte zunächst versucht, nur mit dem windows command script auszukommen, meistens reichen die Möglichkeiten dieser Sprache. Wenn ich sonst darüber hinaus irgendwelche text-basierten Files zu verarbeiten habe, greife ich ganz gerne zu perl, die hab' ich auch im Büro zur Verfügung. An Python habe ich auch zwischendurch gedacht (würde sicher gut mit gdal zusammenspielen), aber da hätte ich mich erst mal wieder in eine neue Sprache reinfinden müssen...

Zitat

Wenn Du nichts dagegen hast, würde ich das Skript ins QLGT Wiki stellen. Dazu fehlt in der Datei aber noch der Urheber und die Lizenz.

Da hab' ich kein ernsthaftes Problem mit, die scripte passen halt 100% erst mal auf mein Problem und müssten für verwandte Probleme noch etwas umgebaut werden (aber das ist ja meist so...). Ob ich noch (nur aus Interesse...) beide scripte in einem perl-script zusammenfasse, weiß ich noch nicht. Andererseits dürfte es auch einfach sein, das windows-script in ein z.B. unix-shell-script umzuschreiben. Was bräuchtes Du noch, um das script ins QLGT WiKi einzustellen? Muss ich da noch eine Art disclaimer in die scripte einbauen?

Was das Problem der Darstellung im GPSmap62 angeht: ich hatte schon mal einen Wegpunkt unter die Kacheln gelegt, um die Kachel zu finden, aber wir konnten die Karte trotzdem nicht sehen. Wir werden den Versuch aber noch mal wiederholen, da wir damals möglicherweise noch korrupte kmz-FIles hatten...

@Anton: wenn ich da nicht weiterkomme, lade ich mal eine Beispielkachel mit Wegpunkt hier hoch...

Die Referenzierung der kmz-Files sollte stimmen, jedenfalls liegen die Kacheln, wenn man sie nach GoogleEarth lädt, genau an der richtigen Stelle.

Einstweilen herzliche Grüße,
Thomas