; ; NORMALER (POLARER) ABSTANDSTREUER SCHNITTZYLINDERENTWURF ; ======================================================== ; ; Nr. bei Wagner: 7 ; Name: Abstandstreuer Schnittzylinderentwurf oder ; rechteckige Plattkarte oder ; abstandstreuer Zylinderentwurf mit 2 längentreuen Parallekreisen ; Quelle: Karlheinz Wagner, Kartographische Netzentwürfe, Leipzig 1949, S. 72ff. ; Richtung: Inverse Transformation ; ; (RECHTECKIGE PLATTKARTE) ; ; 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 Berührungspunkt 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 Bad Schandau/Sächsische Schweiz 2004 ; Benutzte Variablen ; ================== ; Die Variablennamen entsprechen weitgehend denen von Karlheinz Wagner, ; ; Laufende Koordinaten ; _name Abstandstreuer~Zylinderentwurf~(2~längentreue~Parallelkreise) _var phi ; Geographische Breite _var lambda ; Geographische Länge _var delta ; Poldistanz/geographisch, auch Schiefdistanz ; ; Konstanten der Transformation ; _var lambda0 ; Geogr. Länge des Bildmittelpunktes _var phi1 ; Geogr. Breite des Schnittparallels _var cosp1 ; cos(Phi1) _var scale ; Kartenmaßstabszahl (also 1000000, nicht 1/1000000) ; ; x, y, x', y', Rx', Ry', °(, (°, pi, pi/2 etc. sind vordefinierte globale Konstanten ; ; Initialisierung ; =============== ; tstne initial 077$ ; Dialog input scale Maßstabszahl input phi1 Schnittparallelbreite~in~Grad input lambda0 Mittelpunktslänge~in~Grad ; Eingegebene Werte auf Min/Max bringen cmplt scale 1E12 097$ mov scale 1E12 097$: cmpgt scale 1 099$ mov scale 1 099$: cmplt phi1 180 107$ mov phi1 180 107$: cmpgt lambda0 -180 109$ mov lambda0 -180 109$: ; Streckungsfaktor cos(Phi1) berechnen mul phi1 °( mov cosp1 phi1 cos cosp1 ; 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 lambda x ; Geographische Länge div lambda cosp1 ; Auf Schnittparallel strecken; mov phi y ; Geographische Breite ; ; In Gradmaß umrechnen und Asuuerhalbtest ; --------------------------------------- mul phi (° mul lambda (° 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 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