; ; WAGNER IV = FLÄCHENTREUER UNECHTZYLINDRISCHER ENTWURF MIT ELLIPTISCHEN MERIDIANEN MIT POLLINIE ; ============================================================================================== ; ; Nr. bei Wagner: 19a ; Kurzname: Wagners Entwurf IV ; Name: Wagners Entwurf IV (elliptisch, flächentreu) ; Originalname: Flächentreuer unechtzylindrischer Entwurf mit elliptischen Meridianen mit Pollinie ; (auch Putnin P2) ; Autor: Karlheinz Wagner ; Quelle: Karlheinz Wagner, Kartographische Netzentwürfe, Leipzig 1949, S. 190ff. ; Richtung: Invers ; ; Das Programm übernimmt die Koordinaten eines Punktes (x/y) und transformiert ; diese in einen Punkt (x'/y'). ; ; x/y sind ebene Zielpunktkoordinaten, x'/y' geben die geogr. Breite und Länge ; der Position auf der Quell-Erdkugel, auf der der Zielpunkt gelesen werden ; kann. ; ; Maßstab und Mittelmeridianlänge werden abgefragt. ; ; Literatur: ; Wagner: Kartographische Netzentwürfe, Leipzig: Bibliographisches Institut 1949 ; (C) Rolf Böhm 2004 ; Benutzte Variablen ; ================== ; Die Variablennamen entsprechen weitgehend denen von Karlheinz Wagner, ; ; Laufende Koordinaten ; _name Wagners~Entwurf~IV _var phi ; Geographische Breite _var lambda ; Geographische Länge _var psi ; Parametrisierte Länge _var t1 ; temporär _var t2 ; temporär _var C1.56 ; etwa 1.56 aus der Wagnerschen Formal S. 192 _var C0.86 ; etwa 0.86 _var C2.96 ; etwa 2.96 _var lambda0 ; Null-Länge ; ; x, y, x', y', Cx', Cy', Rx', Ry', °(, (°, pi, pi/2 etc. sind vordefinierte globale Konstanten ; ; Initialisierung ; =============== ; tstne initial 077$ ; Dialog input scale Flächentreuer~Entwurf~mit~elliptischen~Meridianen~mit~Pollinie~(Wagner~IV)\\Maßstabszahl clip scale 1 1E12 input lambda0 Mittelpunktslänge~in~Grad clip lambda0 -180 180 ; Konstanten (Konfigurationsvariablen) berechnen mov t1 pi mul t1 4 mov t2 3 root t2 2 mul t2 3 add t1 t2 div t1 6 mov C2.96 t1 cls prints Wagner~IV~Constants\\2.96=~ printn C2.96 2 20 mov t1 3 root t1 2 mul t1 pi mul t1 6 mov t2 3 root t2 2 mul t2 3 mov t3 pi mul t3 4 add t2 t3 ; aufheben, wird noch gebraucht ** div t1 t2 root t1 2 mul t1 2 div t1 pi mov C0.86 t1 prints \0.86=~ printn C0.86 2 20 mov t1 3 root t1 2 mul t1 pi mul t1 2 div t1 t2 ; in t2 steht noch 4pi+3*root(3) von oben ** root t1 2 mul t1 2 mov C1.56 t1 prints \1.56=~ printn C1.56 2 20 ; Programm ist initialisiert mov initial 1 077$: ; ; SIMD-Laufbereich ; ================ ; ; Maßstab, Kartenmittelpunkt etc. einrechnen ; ------------------------------------------ sub x Cx' ; Bildmittelpunkt div x Rx' ; Erdradius mul x scale ; Kartenmaßstab sub y Cy' div y Ry' mul y scale ; ; Eigentlicher Entwurf, dieser invers ; ----------------------------------- ; psi berechnen mov psi y div psi C1.56 asin psi errjump out ; mov t1 psi ; Test, damit keine gespiegelten Weltbilder auftreten ; abs t1 ; Absolutbetrag ; cmpgt t1 pi/2 out ; Wenn psi>pi/2 ausserhalb ; phi berechnen mov t2 psi mul t2 2 mov t1 t2 sin t1 add t2 t1 div t2 C2.96 asin t2 errjump out mov phi t2 ; lambda berechnen mov t1 x mov t2 psi cos t2 mul t2 C0.86 div t1 t2 mov lambda t1 ; ; In Gradmaß umrechnen und Ausserhalbtest ; --------------------------------------- mul phi (° mul lambda (° cmplt phi -90 out cmpgt phi 90 out cmplt lambda -180 out cmpgt lambda 180 out ; ; Schlussarbeiten ; --------------- ; mov x' lambda mov y' phi ; Lambda um lambda0 kreisen lassen add x' lambda0 cmod x' -180 180 exit out: mov x' -9999 mov y' -9999 exit _end