; ; EISENLOHR ; ========= ; ; Nr. bei Evenden: 8.1.5 (2005) ; Name: Eisenlohr ; 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 2005 ; Benutzte Variablen ; ================== ; ; Laufende Koordinaten ; _name Eisenlohr _var phi ; Geographische Breite _var lambda ; Geographische Länge _var S1 ; siehe G. I. Evenden _var C1 _var Q _var T _var C _var P _var V _var 1/V _var 3+root8 _var t1 ; temporär _var t2 ; temporär _var ------- ; End of Symbol Table ; ; 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 3+root8 8 ; very sophisticated. 3+root8 is a variable. Thats theis initialization root 3+root8 2 add 3+root8 3 ; 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 S1 lambda ; Evendens Formuls (8.14) div S1 2 mov C1 S1 ; see later sin S1 cos C1 ; (8.15) mov Q phi ; (8.16) div Q 2 mov T Q ; see later cos Q sin T ; (8.17) mov t1 phi cos t1 mul t1 2 root t1 2 mul t1 C1 add t1 Q div T t1 mov t1 T ; (8.18) power t1 2 add t1 1 mov C 2 div C t1 root C 2 mov P phi ; (8.19) cos P div P 2 root P 2 mov V C1 ; (8.20) add V S1 mul V P add V Q mov t1 C1 sub t1 S1 mul t1 P add t1 Q div V t1 ; main division root V 2 mov 1/V 1 ; V reciproke div 1/V V ; x (Formula 8.21) mov t1 V sub t1 1/V mul t1 C mov x V log x mul x -2 add x t1 mul x 3+root8 ; y (Formula 8.22) mov t1 V add t1 1/V mul t1 C mul t1 T mov y T atan y mul y -2 add y t1 mul y 3+root8 ; ; 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