; ; CANTERS LOW ERROR POLYCONIC C8-4 ; ================================ ; ; Name: Canters Low-error simple oblique polyconic projection with ; pointed meta-pole and constant scale along the axes; ; projection 4 with 8 coefficients. ; Quelle: Canters, Small-scale Map Projection Design, Taylor & Francis, London 2002, ; p. 209f., Table 5.13, Fig. 5.24 ; 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. ; Mit schiefer Quell Erde zu rechnen. Zentrum bei 61°N, 24°E, keine Verdrehung. ; Dies wird erreicht mit einer Metatransformation Erde auf Erde mit dem ; Programm meta_pro1_kugeltransformation.rta und den Parametern ; phiP = -29, lambdaP = 24, lambda0 = 180. ; (C) Rolf Böhm 2006 ; Benutzte Variablen ; ================== ; ; Die Variablennamen entsprechen weitgehend denen von Frank Canters ; ; Laufende Koordinaten ; _name Canters~Low~Error~Oblique~Polyconic~C8-4 ; _var phi _var lambda ; _var c10 _var c01 _var c20 _var c11 _var c02 _var c30 _var c21 _var c12 _var c03 _var c40 _var c31 _var c22 _var c13 _var c04 _var c50 _var c41 _var c32 _var c23 _var c14 _var c05 ; _var c'10 _var c'01 _var c'20 _var c'11 _var c'02 _var c'30 _var c'21 _var c'12 _var c'03 _var c'40 _var c'31 _var c'22 _var c'13 _var c'04 _var c'50 _var c'41 _var c'32 _var c'23 _var c'14 _var c'05 ; _var s0 ; The y-Register _var s1 _var s2 _var s3 _var s4 _var s5 ; ; 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 ; Koeffizienten initialisieren mov c10 +0.9457 mov c12 +0.2962 mov c32 +0.0420 mov c14 -0.1875 ; mov c'01 +0.9836 mov c'21 +0.0472 mov c'41 -0.0106 mov c'23 +0.0929 ; 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 clip 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 r0 1 mov r1 lambda mov r2 lambda mov r3 lambda mov r4 lambda mov s0 1 mov s1 phi mov s2 phi mov s3 phi mov s4 phi mov r6 phi cos r6 power r2 2 power r3 3 power r4 4 power s2 2 power s3 3 power s4 4 ; ; The X-Sum ; clr x mov a c10 mul a r1 add x a mov a c12 mul a r1 mul a s2 add x a mov a c32 mul a r3 mul a s2 add x a mov a c14 mul a r1 mul a s4 add x a mul x r6 ; multiply with cos(phi) -- that's x ; ; The Y-Sum ; clr y mov a c'21 mul a r2 mul a s1 add y a mov a c'41 mul a r4 mul a s1 add y a mov a c'23 mul a r2 mul a s3 add y a mul y r6 ; mul with cos(phi) here! mov a c'01 mul a s1 add y a ; y complete ; ; 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