Benutzer Diskussion:Spischot/Archiv/2012

Projektionsvorlage für die Alex Lambert azimuthal equal-area projection?

Mal was anderes: Alx hat da mal wieder gezaubert und eine Karte für Amerika geplotted: vgl. hier (Diese Karte meine ich: Datei:Alexrk2 Arbeitsdatei.svg. Das ist laut Alex Lambert azimuthal equal-area projection (centered at 90° West, 20° North). Hast du da auch eine funktionierende Vorlage in peto? Oder müsste man das erst programmieren und würde jeden Rahmen sprengen für (erstmal) nur eine Karte. Oder anders ausgedrückt: könnte ich daraus eine optisch passende Positionskarte fertigen, die dann auch Chancen hätte in einer Poskartenvorlage ansteuerbar ist? --TUBS  11:48, 14. Jul. 2012 (CEST)

Ja, das habe ich heute morgen auch schon mit Besorgnis ;-) gesehen. Nee, hab' nix in peto. Machbar ist es bestimmt; ob für mich (mit vertretbarem Aufwand) machbar kann ich aber noch nicht beurteilen. Vor allem hängt es davon ab, ob ich an eine einigermaßen brauchbare mathematische Umwandlung (lat, long, Projektionsparameter) -> (x, y) komme. Flächentreue Azimutalprojektion gibt nix her. en:Lambert azimuthal equal-area projection oder fr:Projection azimutale équivalente de Lambert scheinen ergiebiger, habe ich aber noch nicht näher betrachtet. Die Angabe (centered at 90° West, 20° North) ist dazu bestimmt unzureichend (es fehlt mindestens was noch zu Ausrichtung und Ausschnitt). --Spischot (Diskussion) 14:02, 14. Jul. 2012 (CEST)
Nagut, dann packe ich die Karte erstmal auch nicht an, da man damit vorläufig nicht viel - zmindest nicht viel in Positionskartenhinsicht - anfangen kann. Grüße. --TUBS  14:13, 14. Jul. 2012 (CEST)
Alex schreibt zur Projektion auch was: [1]--TUBS  14:16, 14. Jul. 2012 (CEST)
Vergiss was ich oben gesagt habe! Ich habe aufgrund der verschiedenen Bezeichnungen etwas den Überblick verloren. Die Notwendige Projektionsvorlage ist {{Positionskarte/Flächentreue Azimutalprojektion}} und existiert schon länger. Neben dem bereits besprochenen Projektionszentrum (hier 90°W, 20°N) werden dazu noch die Position des Projektionszentrum auf der Karte (Standardwert ist die Mitte bei (50%, 50%)) sowie die Skalierungsfaktoren (dx, dy) benötigt. In Alex' Beschreibung sind dies (faktorH,faktorV). Damit ist die Vorlage im Nu erstellt. --Spischot (Diskussion) 14:52, 14. Jul. 2012 (CEST)

File:Americas laea location map.svg. Ich hab Alex auch nochmal gebeten die Parameter in der Dateibeschreibung nachzutragen und dann nochmal hier vorbeizuschauen und dir alles sonst noch benötigte mitzuteilen...Hoffe so kriegen wir eine Poskarte zusammen.--TUBS  16:20, 14. Jul. 2012 (CEST)

Schaffen wir gemeinsam bestimmt. Ich habe schon mal mit grob zusammengeratenen Parametern begonnen: {{Positionskarte Amerika (laea)}}. Mit den genauen Daten von Alex wird es dann besser. --Spischot (Diskussion) 20:25, 14. Jul. 2012 (CEST)
Hallo, die Eckpunkte der Karte sind im Zielkoordinatensystem: Links-Unten -7464417, -8174417; Rechts-Oben 7264417, 8464417 (false_easting und false_northing ist beides 0) ...sorry für die krummen Werte, hätte ich gewusste, dass ihr da jetzt gleich eine Poskarte draus macht, hätte ich sie vorher noch etwas gerundet. Entspr. der besagten Anleitung müssten die benötigten Formalparameter über den Weg 2 zu ermitteln sein (Die Kartenbegrenzung ist in Form von Koordinaten im projizierten Koordinatensystem der Karte bekannt) --alexrk (Diskussion) 23:16, 14. Jul. 2012 (CEST)
Die Positionskartenvorlage {{Positionskarte Amerika (laea)}} ist jetzt fertig. Danke euch beiden! --Spischot (Diskussion) 07:29, 15. Jul. 2012 (CEST)
Hast du auch noch Längen- und Breitengrad der beiden Eckpunkte? Dann könnte ich das Musterbeispiel Vorlage:Positionskarte/Flächentreue Azimutalprojektion#Weg 1: Umrechnung aus Längen- und Breitengrad der Eckpunkte mit dem Kommandozeilenwerkzeug Proj im VNR vervollständigen. --Spischot (Diskussion) 12:46, 15. Jul. 2012 (CEST)

Es gibt an der ganzen Geschichte glaube ich noch einen Schönheitsfehler: wenn ich recht verstehe, verwendet die Poskarten-Vorlage einfachheitshalber sphärische Formeln. Wenn ich die Karte mit einem GIS ausplotte, verwende ich dort hingegen idR. jedoch eine Ellipsoid-Definition. Für die Zwecke der Poskarten scheint das egal zu sein, da bei den kleinen Karten kaum ein Unterschied zu sehen ist. Wenn man aber mal die Zahlen durchrechnet bzw. mit PROJ transformiert, erkennt man halt, dass es doch etwas komplizierter ist und die Rechnungen nicht zu 100% ganz aufgehen. Also in unserem Beispiel (dummerweise habe ich grad noch gesehen, dass ich als PRJ-Datei für die Amerika-Karte einen Hayford-Ellipsoid drin hatte)...

# Ellipsoidisch
proj -I +proj=laea +lat_0=20 +lon_0=-90 +ellps=intl +units=m +no_defs -f "%.6f"
-7464417 -8174417
155.629628      -50.615445
7264417 8464417
52.603980       25.098416
# Sphärisch
proj -I +proj=laea +lat_0=20 +lon_0=-90 +a=6378137 +units=m +no_defs -f "%.6f"
-7464417 -8174417
156.097043      -50.495989
7264417 8464417
52.398066       25.090463

Genau genommen, müsste man also noch beim Ausplotten darauf achten, dass man im GIS als SPHEROID-Parameter eine Kugel in der Projektion verwendet. Aber wie gesagt, ich vermute den Unterschied wird man erst ab Kartengrößen > 2000 Pixel erkennen. --alexrk (Diskussion) 17:33, 15. Jul. 2012 (CEST)

Vielen Dank für die Parameter. Ich habe sie in die Beschreibung aufgenommen und denke, dass sie nun auch für Nichteingeweihte nachvollziehbar ist. Ich habe die Werte der spärischen Projektion verwendet und so gerundet, dass keine „dummen Fragen“ aufkommen. Mach’ deine Karten ruhig weiter mit dem besseren Elispsoid. Ich stimme dir zu, dass man den kleinen Fehler der vereinfachten Positionsanzeige in Kauf nehmen kann. Auf das letzte Quäntchen Präzision kommt es mir nicht an. Mir geht es darum, dass es nun jedem möglich ist, Positionskartenvorlagen für neue Positionskarten selbst zu erstellen. Auf Anfrage werde ich natürlich auch in Zukunft helfen. --Spischot (Diskussion) 20:19, 15. Jul. 2012 (CEST)
Das sieht toll aus. Kann da jmd. von euch beiden in die Dateibeschreibung noch irgendwas sinnvolles rein schreiben? Ich habe da ja nur so Platzhalter reingeschrieben ohne zu wissen ob das die richtigen Parameter sind (hab die glaub ich aus der Russland-laea-Poskarte kopiert, weil ich dachte das man das noch bracuhen könnte). Dann steht da was sinnvolles.--TUBS  23:25, 15. Jul. 2012 (CEST)
Vermutlich weiß Alex am besten, mit welchen Parametern er die Karte erzeugt hat, daher lasse ich ihm mal den Vortritt. Falls er in den kommenden Tagen nicht dazu kommt (oder schlicht keine Lust hat) werde ich ersatzweise die Angaben eintragen, soweit ich sie verstanden habe. --Spischot (Diskussion) 09:08, 16. Jul. 2012 (CEST)
Hab die Werte mal ergänzt. Auf die GMT-Infos würde ich vlt verzichten (wüsst ich jetzt auch nicht, was man da angeben muss) --alexrk (Diskussion) 20:03, 16. Jul. 2012 (CEST)
TipTop. Danke --TUBS  20:06, 16. Jul. 2012 (CEST)

Positionskarten

Moin! Es gibt einige Vorlagen mit „Test“ im Namen, die von dir angelegt wurden. Werden die noch gebraucht?

Viele Grüße, NNW 17:43, 27. Jan. 2012 (CET)

Moin! Überwiegend ja, weil sie derzeit die einzigen Beispiele für die Projektionsvorlagen
sind. Ich habe mich nicht getraut, existierende, ausverformelte Positionskartenvorlagen auf die obigen Projektionsorlagen umzustellen - daher habe ich mir eigene Beispiele geschaffen. Sobald ein echtes Beispiel vorhanden ist, kann je eine der Test-Positionskarten weg. {{Positionskarte Arktischer Ozean (Test)}} ist dagegen jetzt schon überflüssig; ich stelle sofort einen sla.
Der Vollständigkeit halber zur Info:
Mit deiner Frage erinnerst du mich aber daran, dass ich die Verwendung noch mal ordentlich dokumentieren wollte. Mach ich dieser Tage mal. --Spischot 19:15, 27. Jan. 2012 (CET)
Die Verwendung ist jetzt dokumentiert. --Spischot 12:25, 29. Jan. 2012 (CET)
Nachtrag: Durch {{Positionskarte Amerika (laea)}} wurden nun auch die beiden Beispiele {{Positionskarte Afrika (Test)}} und {{Positionskarte Europa (Test)}} überflüssig und zum Löschen vorgeschlagen. --Spischot (Diskussion) 13:23, 15. Jul. 2012 (CEST)
Nachtrag 2: Durch {{Positionskarte Brasilien Pernambuco}} wurde nun auch das Beispiel {{Positionskarte Japan (Test)}} überflüssig und zum Löschen vorgeschlagen. --Spischot (Diskussion) 12:38, 21. Jul. 2012 (CEST)

Neue Positionskarten

Moin! Die Franzosen haben fr:Modèle:Géolocalisation/Yukon, fr:Modèle:Géolocalisation/Territoires du Nord-Ouest und fr:Modèle:Géolocalisation/Nunavut neu im Angebot. Ich habe mich bisher nicht wirklich mit den Vorlagen für Kegelprojektionen auseinandergesetzt. Könntest du Vorlage:Positionskarte Kanada Yukon, Vorlage:Positionskarte Kanada Nunavut und Vorlage:Positionskarte Kanada Nordwest-Territorien anlegen? Viele Grüße, NNW 20:05, 16. Aug. 2012 (CEST)

Au fein! Die drei Karten verwenden die Projektionsart fr:Modèle:Géolocalisation/Projection conique avec DL, zu der ich vermutlich zuerst hier hoch die passende Projektionsvorlage erstellen muss. Ich schau mir das am Wochende mal an. --Spischot (Diskussion) 06:57, 17. Aug. 2012 (CEST)
Klasse! Ich werde dann alternative Kartenversionen erstellen, die den hiesigen Farbvorgaben entsprechen. NNW 09:51, 17. Aug. 2012 (CEST)
So, die Projektionsvorlage {{Positionskarte/Längentreue Kegelprojektion A}} und die drei oben erbetenen Vorlagen sowie zusätzlich noch {{Positionskarte Kanada Québec (Kegelprojektion)}} sind fertig. Würdest du bitte abschließend noch mal hier deine Meinung beisteuern? --Spischot (Diskussion) 16:34, 18. Aug. 2012 (CEST)
Wunderbar, vielen Dank! NNW 10:26, 19. Aug. 2012 (CEST)

Umwandlungen in Positionskarten

Tachauch Spischot,

da Du offenbar der Hauptersteller der ganzen Kegelprojektionsvorlagen bist, dürftest Du der sinnvollste Ansprechpartner für mein Anliegen sein.

Mein Ziel ist es, die Karte Naturraeumliche Grossregionen Deutschlands plus.png und Ausschnitte daraus als Positionskarten verfügbar zu machen. Da die Karte naturgemäß nicht wenige Längen- und Breitengrade umfaßt (und definitiv nicht winkeltreu ist), dürfte sie mit Kegelprojektionen am besten zur Positionskarte umfunktionierbar sein. Eine Annäherung als Plattkarte findet sich hier. Dort zeigt sich auch - die eingezeichneten Punkte sind allesamt solche, deren Soll-Position auf der Karte eindeutig sein dürfte - zeigt die Plattkartenproblematik ganz gut. Und bei Ausschittskarten, die den linken oder rechten Rand zeigen - z. B. bei dieser (alle Punkte links unten müßten deutlich nach links oben wandern - zeigt sich das ebenfalls eindringlich.

Ich bin aber bislang zu blöd, Deine Vorlagen zu verstehen. Wie genau muß ich - ausgehend von einer grob hinhauenden Plattkarte - die Parameter wählen, um möglichst in endlicher Zeit eine möglichst deutlich bessere Annäherung per Kegelprojektion hinzukriegen, die man dann auch auf Ausschnittskarten anwenden könnte? Wenn Du mir da Seiten empfehlen könntest, läse ich mich auch gerne näher ein - die WP-Artikel helfen mir da aber nicht wirklich weiter. Als gelernter Mathefritze müßte ich das eigentlich verstehen können, indes habe ich mich mit Kartographie bislang null beschäftigt (und bin auch kein Fan von Klammerbergen, weshalb mir Deine Vorlagen eigentlich helfen können müßten).

LieGrü, --Elop 22:41, 23. Aug. 2012 (CEST)

Inzwischen stelle ich mich doch nicht mehr so blöd an - ich mußte halt beim n-Wert experimentell vorgehen. Bis auf den Osten ist jetzt eigentlich alles da, wo es sein soll - siehe Positionskarte Deutschland Naturraum. Aber vielleicht hast Du ja trotzdem noch ein paar Tipps ....
LieGrü --Elop 20:07, 24. Aug. 2012 (CEST)
Naja, die Methode, des unbekümmerten Herumprobierens will mir nicht so recht gefallen. Von den Experten der französischen Kartenwerkstatt wird zwar dies im Regelfall so gemacht aber die teilweise unpräsisen Ergebnisse zeigen, dass das nicht der Weisheit letzter Schluss ist.
Methode ganz grob:
  1. Bestimme die Projektion (so gibt es beispielsweise gleich mehrere Kegelprojektionen)
  2. Bestimme die Parameter der Projektion (geht leicht, wenn man weiß, wie die Karte erstellt wurde, schwieriger, wenn die Karte ein Gradnetz hat, bestimmt mühsam, wenn man nur eine Reihe von Stützpunkten hat, und nahezu aussichtslos, wenn die erste Aufgabe nicht erfüllt wurde.
  3. Erstelle die Positionskartenvorlage, ggfs. müssen die Parameter dazu etwas umgerechnet werden.
  4. Teste anhand von Stichproben.
Im obigen Fall haben wir aber den glücklichen Fall, dass die Erzeugung der Originalvorlage vom Autor dankenswerterweise präzise dokumentiert wurde und damit Projektionsart und Parameter gut bestimmbar sind. Demnach ist die Projektion
-JL10/51/42/54/18c
was gemäß der Dokumentation von pscoast
  • eine Lambert Conic Conformal-Projektion
  • mit longitude0=Zentralmeridian=10° E,
  • latitude0=51° N,
  • latitude1=Erstes Standardparallel=42° N,
  • latitude2=Zweites Standardparallel=54° N und
  • scale=18 Zentimeter(?) pro Grad
bedeutet. Der Kartenausschnitt ist
-R5/47/16/55r
und reicht damit von
  • 47° N, 5° O in der linken unteren Ecke bis
  • 55° N 16° O in der rechten oberen Ecke.
Deine Ausschnitte haben die natürlich gleiche Projektion; du müstest allerdings jeweils die beiden Ecken der Ausschnitt neu bestimmen.
Soweit die guten Nachrichten.
Die schlechte Nachricht ist, dass es unter den Kegel-Projektionsvorlagen bislang noch keine {{Positionskarte/Lambertsche Schnittkegelprojektion}} mit den obigen Parametern gibt. --Spischot (Diskussion) 21:44, 24. Aug. 2012 (CEST)
Die korrekte Projektionsvorlage ist in Arbeit. Ich melde hier wenn sie fertig ist. Vielleicht gelingt es mir am Ende auch noch, die lästige Umrechnung der linken-unteren und rechten-oberen Ecke in die Paramter px0, py0, dx und dy zu automatisieren – hierauf bitte zunächst nicht so viel Zeit investieren.
Ich habe gesehen, dass du zwischenzeitlich mit der {{Positionskarte/Lineare Kegelprojektion}} begonnen hast. Diese bildet eine Längentreue Kegelprojektion ab, während deine Karten eine Lambertsche Schnittkegelprojektion haben. Die Abbildung der Positionen wird also niemals ganz stimmen können. Wie groß der Fehler am Ende ist, kann ich allerdings nicht sagen. --Spischot (Diskussion) 07:20, 26. Aug. 2012 (CEST)
Ich glaube, das Hauptproblem ist, daß unsere Kartenbastler in der Regel nicht viel mit Mathematik am Hut haben. Die meisten unserer Positionskartenvorlagen wurden ja von anderen Leuten erstellt.
Überhaupt ist es merkwürdig, daß die meisten unserer großflächigen Positionskarten Plattkarten sind - was ja die flächenmäßigen Proportionen zwischen Norden und Süden erheblich verzerrt!
Ich bin eigentlich nur auf die Idee gekommen, die Naturraumkarte und Ausschnitte zu Poskarten zu machen, weil Morty die Vorlage:Infobox Naturraum in Deutschland erstellt hat, die eine Poskarte vorsieht - was bei den Hauptgruppen es nahe legt, da je gleich breite Ausschnitte der Deutschlandkarte zu verwenden. Dann hat er die Vorlage:Positionskarte Deutschland Naturraum Alpenvorland und Alpen als Plattkarte erstellt, was in dem Bereich (längengradtechnische Mitte der Karte, wenig Nord-Süd-Ausdehnung) auch nicht so fehlerbehaftet ist. Aber schon im Westen ist das nicht sinnvoll übernehmbar, wie ich feststellte.
Daß es mehrere Kegelprojektionen gibt, ist schon klar, aber ich habe von vornherein nur mit denen experimentiert, für die Du vereinfachende Vorlagen geschrieben hattest. Blöd ist eben auch, daß die WP-Artikel nicht viel über die konkreten Formeln verraten - die müßte man sich ergo selber herleiten.
Was ich mir eigentlich erhofft hatte, war, daß, sobald die Deutschlandkarte positioniert wäre, die Ausschnittskarten problemlos per Dreisatz vom Himmel fielen. Das ging aber mit der linearen Kegelprojektion auf jeden Fall nicht! Da muß mindestens das n erheblich angepaßt werden.
Zum Punkt 4. Teste anhand von Stichproben.:
Ich habe Dir mal zu Vorlage Diskussion:Positionskarte Deutschland Naturraum Analoges dorthin übertragen.
Danke für Deine Mühen bis hierher schon mal! --Elop
Da laut Test Die von Dir ermittelten Werte hinzuhauen scheinen, habe ich sie mal in die aktuellere Karte übertragen.
Jetzt würde sich die Frage nach der Rückwärtsformel stellen, damit man z. B. "links unten, 1000 Pixel breit und 1016 hoch, also die hier, direkt umrechnen könnte (jene stimmt ja, siehe Test, an der Luxemburger Grenze noch sehr schlecht mit der hier untauglichen Plattkarte). Hast Du dafür die Formeln parat? --Elop 13:52, 8. Sep. 2012 (CEST)