; ; FOURNIER GLOBULAR I. ; ==================== ; ; Nr. bei Evenden: 8.1.6 (2005) ; Name: Fournier Globular I. ; Quelle: Evenden, G. I.: Libroj4 Documentation ; Richtung: Direkt-Transformation ; ; 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. ; (C) Rolf Böhm 2004 ; Benutzte Variablen ; ================== ; ; Laufende Koordinaten ; _name Fournier~Globular~I _var phi ; Geographische Breite _var lambda ; Geographische Länge _var phi² _var lambda² _var C ; siehe G. I. Evenden _var P _var S _var S² _var A _var t1 ; temporär _var t2 ; temporär ; ; Konstanten der Transformation ; _var lambda0 ; Geogr. Länge des Bildmittelpunktes _var scale ; Kartenmaßstabszahl (also 1000000, nicht 1/1000000) ; ; x, y, x', y', Cx', Cy', °(, (°, pi, pi/2 etc. sind vordefinierte globale Konstanten ; ; Initialisierung ; =============== ; tstne initial 077$ ; Dialog pause Hinweis:~Dieses~Programm~rechnet~eine~Vorwärtstransformation.\\Es~muss~mit~einer~direkt~arbeitenden~Projection~engine~abgearbeitet~werden. input scale Maßstabszahl input lambda0 Mittelpunktslänge~in~Grad ; Eingegebene Werte auf Min/Max bringen clip scale 1 1E12 clip lambda0 -180 180 clip phi0 -90 90 mov C pi power C 2 div C 4 ; Programm ist initialisiert mov initial 1 077$: ; ; SIMD-Laufbereich ; ================ ; ; Eigentlicher Entwurf, dieser direkt! ; ------------------------------------ ; mov lambda x ; Geographische Länge mov phi y ; Geographische Breite sub lambda lambda0 cmpgt lambda -180 10$ add lambda 360 10$: cmpgt lambda -180 30$ add lambda 360 30$: cmplt lambda 180 40$ sub lambda 360 40$: cmplt lambda 180 50$ sub lambda 360 50$: ; ; Umrechnung in Bogenmaß ; ---------------------- ; mul phi °( mul lambda °( ; Netzentwurf rechnen ; ------------------- mov phi² phi power phi² 2 mov lambda² lambda power lambda² 2 mov P phi ; P (8.27) sin P mul P pi abs P mov t2 phi ; S (8.28) abs t2 mul t2 2 mov t1 P sub t1 t2 mov S phi² neg S add S C div S t1 ; main division mov S² S power S² 2 mov A lambda² ; A (8.29) ; power A 2 ; 1st Error "Freias birthday" div A C sub A 1 mov t1 P ; y (8.30) mul t1 S neg t1 sub t1 lambda² add t1 C mul t1 A neg t1 add t1 S² root t1 2 sub t1 S div t1 A mov y phi sgn y mul y t1 ; mov t1 y ; x (8.31) ; 2rd Error power t1 2 ; 3rd Error div t1 C neg t1 add t1 1 root t1 2 mul t1 lambda ; here already the lambda sign ... ; mov x lambda ; ... the sign already correct from the mul?? But here the written formula ... ; sgn x ; mul x t1 mov x t1 ; ; Maßstab, Kartenmittelpunkt etc. einrechnen ; ------------------------------------------ ; mul x Rx' ; Erdradius div x scale ; Kartenmaßstab add x Cx' mul y Ry' div y scale add y Cy' ; ; Schlussarbeiten ; --------------- ; 111$: mov x' x mov y' y exit out: mov x' -9999 mov y' -9999 exit _end