; ; CANTERS LOW ERROR PSEUDOCYLINDRIC IV ; ==================================== ; ; Name: Canters low-error pseudocylindric IV ; (Pointed pole, correct axis ratio) ; "C6-4" ; Quelle: - G. I. Evenden, Libproj4 (2005) S. 70, 5.2.44 ; - Frank Canters: Small-scale Map Projection Design, London 2002, p. 207, ; Table 5.12, Fig. 5.23 ; 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 2006 ; Benutzte Variablen ; ================== ; ; Laufende Koordinaten ; _name Canters~Low~Error~Pseudocylindric~IV~(C6-4) _var phi ; Geographische Breite _var lambda ; Geographische Länge _var c0 _var c2 _var c4 _var c'1 _var c'3 _var c'5 ; ; 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 ; Konstanten mov c0 0.8333 mov c2 0.3385 ; -0.3385 gives Evendens Fig. 5.11 D mov c4 0.0942 mov c'1 1.0114 mov c'3 0.0243 mov c'5 -0.0391 ; 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 cmod lambda -180 180 ; ; 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 r1 phi mov r2 phi mov r3 phi mov r4 phi mov r5 phi power r2 2 power r3 3 power r4 4 power r5 5 mov r0 c0 mul r1 c'1 mul r2 c2 mul r3 c'3 mul r4 c4 mul r5 c'5 mov x r0 add x r2 add x r4 mul x lambda mov r7 phi ; only canters low error pseudocylindrical iv cos r7 mul x r7 mov y r1 add y r3 add y r5 ; ; 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 _end