; ; SEMICONFORMAL ; ============= ; ; Name: Semiconformal ; Proj4: +proj=near_con ; Libproj4 (2005), S. 69, 5.2.40 ; Quelle: Gerald I. Evenden. Libproj4: A Comprehensive Library if Cartographic Projection Functions, March 2004 ; Richtung: Direkttransformation ; ; Maßstab und Berührungsparallel werden abgefragt. ; 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 ; Fiala: Mathematische Kartographie, Berlin: Verlag Technik 1957 ; (C) Rolf Böhm 2006 ; Benutzte Variablen ; ================== ; ; Die Variablennamen entsprechen weitgehend denen von Karlheinz Wagner, ; ; Laufende Koordinaten ; _name Semiconformal _var phi ; Geographische Breite _var lambda ; Geographische Länge ; ; Konstanten der Transformation ; _var scale ; Kartenmaßstabszahl (also 1000000, nicht 1/1000000) _var lambda0 ; Geogr. Länge des Bildmittelpunktes _var p _var theta _var 2pi ; ; 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 2pi pi mul 2pi 2 ; Programm ist initialisiert mov initial 1 077$: ; ; SIMD-Laufbereich ; ================ ; ; Lambda und Phi übernehmen ; ------------------------- ; 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 ; ------------------- ; ; p ; mov r0 phi abs r0 cmpgt r0 1.5564 130$ mov p phi sin p jump 150$ 130$: mov r0 phi sgn r0 mov p 0.99989 mul p r0 150$: ; ; theta ; mov theta 1 add theta p mov r0 1 sub r0 p div theta r0 log theta div theta 2pi ; ; x ; mov x theta cos x mul x lambda ; ; y ; mov y theta sin y mul y pi ; ; 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