; ; LOXIMUTHAL ; ========== ; ; Name: Loximuthal ; Proj4: +proj=loxim ; Libproj4, S. 66, 6.2.33 ; Quelle: Gerald I. Evenden. Libproj4: A Comprehensive Library if Cartographic Projection Functions, March 2004 ; Richtung: Inverse 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 Mittellänge werden abgefragt. ; ; Weiterhin wird eine Bezugsbreite abgefragt. Alle Geraden, die durch den Schnittpunkt ; Mittelmeridian-Bezugsparallel verlaufen, sind längentreue Loxodromen. ; ; Die Zielbildgeometrie wird dem Sekundäroperanden oder einem Fixbild entnommen. ; Der Berührungspunkt wird bildmittig gesetzt. ; ; Literatur: ; Wagner: Kartographische Netzentwürfe, Leipzig: Bibliographisches Institut 1949 ; Gerald I. Evenden, USGS: Libproj4. A Comprehensive Library if Cartographic Projection Functions, March 2004 ; (C) Rolf Böhm 2004 ; Benutzte Variablen ; ================== ; ; Laufende Koordinaten ; _name Loximuthal _var phi ; Geographische Breite _var theta ; Parameter _var lambda ; Geographische Länge _var t1 ; temp _var t2 ; temp _var t3 ; temp ; ; Konstanten der Transformation ; _var scale ; Kartenmaßstabszahl (also 1000000, nicht 1/1000000) _var phi0 ; Geogr. Breite des Bildmittelpunktes _var lambda0 ; Geogr. Länge des Bildmittelpunktes _var lntan.. ; ln(tan(pi/4+phi0/2)) _var initial _var .......... ; Ende Symboltabelle ; ; x, y, x', y', Cx', Cy', °(, (°, pi, pi/2 etc. sind vordefinierte globale Konstanten ; ; Initialisierung ; =============== ; tstne initial 077$ ; Dialog input scale Loximuthal~(Alle~Geraden~durch~den~Mittelpunkt~werden~längentreue~Loxodromen)\\Maßstabszahl input phi0 Mittelpunktsbreite~in~Grad input lambda0 Mittelpunktslänge~in~Grad ; Eingegebene Werte auf Min/Max bringen clip scale 1 1E12 clip phi0 -89.9999 89.9999 clip lambda0 -180 180 ; Konstanten mul phi0 °( mov lntan.. phi0 div lntan.. 2 add lntan.. pi/4 tan lntan.. log lntan.. ; Programm ist initialisiert mov initial 1 077$: ; ; SIMD-Laufbereich ; ================ ; ; Maßstab, Kartenmittelpunkt etc. einrechnen ; ------------------------------------------ sub x Cx' ; Bildmittelpunkt div x Rx' ; Erdradius mul x scale ; Kartenmaßstab sub y Cy' div y Ry' mul y scale ; ; Eigentlicher Entwurf, dieser invers ; ----------------------------------- ; mov phi y add phi phi0 mov t1 phi sub t1 phi0 mov t2 phi div t2 2 add t2 pi/4 tan t2 log t2 sub t2 lntan.. mov lambda x mul lambda t2 div lambda t1 ; ; In Gradmaß umrechnen und Ausserhalbtest ; --------------------------------------- ; mul phi (° mul lambda (° mul theta (° ; Test cmplt phi -90 out cmpgt phi 90 out cmplt lambda -180 out cmpgt lambda 180 out ; ; Schlussarbeiten ; --------------- ; mov x' lambda add x' lambda0 mov y' phi mov z' theta ; Test cmpgt x' -180 10$ add x' 360 10$: cmpgt x' -180 30$ add x' 360 30$: cmplt x' 180 40$ sub x' 360 40$: cmplt x' 180 50$ sub x' 360 50$: exit out: mov x' -9999 mov y' -9999 exit _end