; ; RAISZ-ARMADILLO ; =============== ; ; Nr. bei Evenden: 8.1.3 (2005) ; Name: Raisz Armadillo ; 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. ; ; Maßstab und Berührungspunkt werden abgefragt. ; Die Zielbildgeometrie wird dem Sekundäroperanden oder einem Fixbild entnommen. ; Der Berührungspunkt wird bildmittig gesetzt. ; ; (C) Rolf Böhm 2004 ; Benutzte Variablen ; ================== ; ; Laufende Koordinaten ; _name Raisz~Armadillo _var phi ; Geographische Breite _var lambda ; Geographische Länge _var phi-s ; Phi(s) (Formel (8.5)) _var phi0 ; Berührungsbreite; im Original 20° _var sinphi0 ; Sinus davon _var cosphi0 ; Cosinus davon _var tanphi0 ; Tangens davon _var lambda/2 ; Lambda/2 _var coslambda/2 ; Cos(Lambda/2) _var cosphi ; Cos(Phi) _var sinphi ; Sin(Phi) _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 phi0 Zentralbreite~in~Grad~(Raisz'~Originalwert:~20) input lambda0 Mittelpunktslänge~in~Grad ; Eingegebene Werte auf Min/Max bringen clip scale 1 1E12 clip lambda0 -180 180 clip phi0 -90 90 mul phi0 °( mov tanphi0 phi0 tan tanphi0 mov sinphi0 phi0 sin sinphi0 mov cosphi0 phi0 cos cosphi0 ; 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 ; ------------------- ; Phi(s) mov lambda/2 lambda div lambda/2 2 mov coslambda/2 lambda/2 cos coslambda/2 mov phi-s coslambda/2 div phi-s tanphi0 atan phi-s neg phi-s tsteq phi0 057$ ; eq: Äquator tstlt phi0 055$ ; lt: Südhalbkugel cmpgt phi-s phi out ; Out-test Nordhalbkugel jump 057$ 055$: cmplt phi-s phi out ; Out-test Südhalbkugel 057$: ; sin(Phi), Cos(Phi) mov sinphi phi sin sinphi mov cosphi phi cos cosphi ; x mov t1 lambda/2 sin t1 mov x 1 add x cosphi mul x t1 ; y mov y 1 ; 1. Summand add y sinphi0 sub y cosphi0 div y 2 mov t1 sinphi ; 2. Summand mul t1 cosphi0 add y t1 mov t1 1 ; 3. Summand add t1 cosphi mul t1 sinphi0 mul t1 coslambda/2 sub y 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