Ich habe ein DOP1,6m meiner Stadt beim LVA Lizenziert und heute per Email zugestellt bekommen.
Inhalt: eine TIFF Datei, tfw, wld und txt. Die letzteren drei enthalten nur paar Zahlenwerte welche die Koordinate des Bildausschnittes wiedergeben.
Das TIFF habe ich in ein JPG umgewandelt und dann mit JOSM und dem PicLayer Modul reingeladen. Das Foto entsprechend skaliert damit es deckungsgleich ist. Soweit alles toll, nur je weiter ich von dem zentrierten Punkt mich in Josm entferne, um so größer werden die Abweichungen von Luftbild und OSM Daten. Das sind dann Abweichungen von bis zu 100m. Das Luftbild kann ich an diesen Abweichstellen deckungsgleich hinschieben und verlagere das Problem dann bei einem anderen Bereich.
Sieht für mich nach dem selben Problem aus, welches damals dieser PlanAT Type hatte, wo auch nichts so 100% richtig passte.
Wie entzerre ich diesen Scheiß in die richtige Projektion, damit das Luftbild wie bei dem WMS aussieht?
einfach: http://warper.geothings.net/ da anmelden, uploaden (als private) kontrollpunkte mit bereits vorhandener OSM karte setzen und die url mit wms-plugin von josm einspielen.
ist zwar ungenau (immerhin verwendest du zur korrektur ungenaue osm kartenpunkte) … aber eben doch am schnellsten.
naja das tiff könntest auch mit gdal (linux) bzw. den fwtools (gdal für windows) umrechnen. dazu brauchst halt deine projektion (dürfte gauss-krüger koordinatensystem sein) und dein datum damit du es auf wgs84+mercator umrechnen kannst. du könntest aus dem tiff (ich nehme an es ist ein geoTiff) genau diese werte auslesen und die projektion umrechnen (alles mit gdal).
Edit: Trotzdem müsstest du glaub ich das geoTiff schlussendlich in einem lokalen WMS-Server verwenden, damit du das genau wie bei JOSM verwenden kannst. Arbeitsaufwand geschätzt: 8h (solltest du nicht wissen wie das alles funktioniert). mit dem map-warper lässt sich das ganze (wenn auch ungenau) in ca. 30 Minuten erledigen.
Ich habe geothing ausprobiert und nachdem ich 30 Punkte gesetzt habe, passt es und alles ist weitestgehend deckungsgleich. Dafür ist die Bildqualität die mir geothing liefert unter aller Sau und in der Auflösung stark verschlechtert.
Das Original hat 3344x2709 Pixel und der Export von Geothing liefert nur 1713x953
Dann rechnen Sie die Längen-/Brei"-tengrad-Rand"-koordinaten der Längen-/Brei"-tengrad-location “`naxosll”’ mit dem Modul $ m.proj in das UTM-System um… Ich verstehe nur Bahnhof
Mir geht’s da auch nicht besser; das grenzt an Arbeit
Was ist mit der anderen Idee, also die Projektionsmethode im JOSM zu verstellen (z.B. von Mercator auf WGS84). Da müssten doch eigentlich die Verzerrungen deutlich geringer werden…
Tim von geothings meinte, ich soll das Luftbild stückeln, dann würde der Server das nicht mehr künstlich verkleinern. Das habe ich nun gemacht und es funktioniert. Das Ergebnis ist optisch immer noch nicht so schön wie das Original, aber man kann damit arbeiten.
Merkator passt im Zentrierten Bereich, darum herum gibt es immer größer werdende Abweichungen. WGS84 bringt nichts brauchbares und das restliche Zeug wie Estland, Frankreich, Schweiz Gitter brauch ich nicht testen. Die UTMs liefern auch nur Müll.
Laut Freund Google soll das Zeug in Gauss-Krüger projeziert sein, was laut Google das Programm Josm nicht kann.
Falls ja: ruf nochmal an, entschuldige Dich für den Fehler, und frage ob sie Dir vielleicht die Bilder nochmal in WGS-84 schicken könnten?
Das ist für die Spezialisten vom LVA-Bayern ganz einfach und da machen sie sicher gern wenn Du nett fragst.
Driver: GTiff/GeoTIFF
Files: fuck.tif
Size is 3344, 2709
Coordinate System is:
PROJCS[“unnamed”,
GEOGCS[“unnamed”,
DATUM[“unknown”,
SPHEROID[“unretrievable - using WGS84”,6378137,298.257223563]],
PRIMEM[“Greenwich”,0],
UNIT[“degree”,0.0174532925199433]],
UNIT[“metre”,1,
AUTHORITY[“EPSG”,“9001”]],
AUTHORITY[“EPSG”,“31468”]]
Origin = (4339136.179999999700000,5278047.320000000300000)
Pixel Size = (1.600000000000000,-1.600000000000000)
Metadata:
AREA_OR_POINT=Area
TIFFTAG_SOFTWARE=Geodaten Bestellung von Matthias Schulze
TIFFTAG_DATETIME=09.06.2006 00:00:00
TIFFTAG_COPYRIGHT=Landesamt für Vermessung und Geoinformation
TIFFTAG_XRESOLUTION=31
TIFFTAG_YRESOLUTION=31
TIFFTAG_RESOLUTIONUNIT=3 (pixels/cm)
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 4339136.180, 5278047.320)
Lower Left ( 4339136.180, 5273712.920)
Upper Right ( 4344486.580, 5278047.320)
Lower Right ( 4344486.580, 5273712.920)
Center ( 4341811.380, 5275880.120)
Band 1 Block=3344x8 Type=Byte, ColorInterp=Red
Band 2 Block=3344x8 Type=Byte, ColorInterp=Green
Band 3 Block=3344x8 Type=Byte, ColorInterp=Blue
Ich hab mal deinen “So oder so ähnlich Code” ausprobiert. Da kommen nur Fehler 1 oder Fehler 4. Hab paar Parameter geändert, die `in " geändert, da epsg:31468 oder nur 31468 eingegeben. Bringt alles nix.
Im Web hab ich dann das da gefunden. gdalwarp -of GTiff -co “TILED=YES” -srcnodata 31468 -dstnodata wgs84 -t_srs “+proj=merc +ellps=sphere +R=6378137 +a=6378137 +units=m” -rcs -order 3 -multi fuck.tif warped.tif
Da arbeitet dann was von 0-100% und erzeugt eine TIF die leicht gedreht ist. Deckungsgleich ist da aber immer noch nichts Was die ganzen Parameter bedeuten, keine Ahnung
Was die ganzen Parameter bedeuten, keine Ahnung