; ; LAMBERT CONFORMAL ; ================= ; ; Name: Lambert Conformal ; Quelle: Eric M. Delmelle: Map Projection Properties, S. 79 ; 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) Dr.-Ing. Rolf Böhm 2005 ; Benutzte Variablen ; ================== ; ; Laufende Koordinaten ; _name Lambert~Conformal _var phi ; Geographische Breite _var lambda ; Geographische Länge _var tanphi/2 ; tan(phi/2) _var sinlambda/2 _var coslambda/2 ; ; Konstanten der Transformation ; _var lambda0 ; Geogr. Länge des Bildmittelpunktes _var scale ; Kartenmaßstabszahl (also 1000000, nicht 1/1000000) _var initial _var -------- ; EofSymbolTable ; ; 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 ; 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$: ; x, y haben eine Doppelbedeutung: ; - einmal die verdefinierten RTA-Input-Koordinaten, also eigentlich Phi und Lambda, ; - dann aber auch die ebenen Kartenkoordinaten, die errechnet werden ; ; Umrechnung in Bogenmaß ; ---------------------- ; mul phi °( mul lambda °( ; Netzentwurf rechnen ; ------------------- ; mov r0 phi div r0 2 tan r0 mov tanphi/2 r0 power r0 2 ; r0: tan²phi/2 mov r1 1 sub r1 r0 root r1 2 ; r1: root(1-tan²(phi/2)) mov r2 lambda div r2 2 mov sinlambda/2 r2 mov coslambda/2 r2 sin sinlambda/2 cos coslambda/2 ; Nun: r7 Zähler, r6: Nenner mov r7 r1 ; Zähler für x mul r7 sinlambda/2 mul r7 2 mov r6 r1 ; Nenner für x mul r6 coslambda/2 add r6 1 mov x r7 ; x Hauptdivision div x r6 mov r7 tanphi/2 ; Zähler für y ; power r7 2 mul r7 2 ; Der Nenner y ist derselbe wie x mov y r7 ; y Hauptdivision div y r6 ; ; 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