; ; LASKOWSKIS MINIMUM ERROR PROJECTION (DELMELLES ZEICHNUNG) ; =================================== ; ; Name: Laskowskis Minimum Error Projection ; Quelle: Eric M. Delmelle. Map Projections Properties. Considerations for ; Small-Scale GIS Applications. State Univ. of New York at Buffalo, ; Dept. of Geography, September 2001 ; Richtung: Direkt-Transformation ; ; Note: There is an picture in ; MAP PROJECTION PROPERTIES: CONSIDERATIONS FOR SMALL-SCALE GIS APPLICATIONS ; by Eric M. Delmelle (2001) on p. 90. ; follows by formulas. But there is an difference between the picture and ; the formulas. I believe the formulas are the Original Laskowski. ; ; This programs instruction #2 is an nop-instruction. Not an power-instruction. ; ; And this program generates the picture. Not computed by the formula. ; ; 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) Dr.-Ing. Rolf Böhm 2005 ; Benutzte Variablen ; ================== ; ; Laufende Koordinaten ; _name Laskowski~Minimum~Error~(Delmelles~Zeichung) _var phi ; Geographische Breite _var lambda ; Geographische Länge _var phi2 ; Quadrate, Kuben ... _var phi3 _var phi4 _var phi5 _var lambda2 _var lambda3 _var lambda4 ; ; 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 ; 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 cmpgt lambda -180 10$ add lambda 360 10$: cmpgt lambda -180 30$ add lambda 360 30$: cmplt lambda 180 40$ sub lambda 360 40$: cmplt lambda 180 50$ sub lambda 360 50$: ; 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 lambda2 lambda mov lambda3 lambda mov lambda4 lambda power lambda2 2 power lambda3 3 ; power lambda4 4 ; SO ERGIBT ES DIE ABBILDUNG IM DELLMELLE (2001) S. 90 !!!! nop mov phi2 phi mov phi3 phi mov phi4 phi mov phi5 phi power phi2 2 power phi3 3 power phi4 4 power phi5 5 ; clr x mov t1 lambda mul t1 0.975534 add x t1 mov t2 lambda mul t2 phi2 mul t2 -0.119161 add x t2 mov t3 lambda3 mul t3 phi2 mul t3 -0.0143059 add x t3 mov t4 lambda mul t4 phi4 mul t4 -0.0547009 add x t4 ; clr y mov t1 phi mul t1 1.00384 add y t1 mov t2 lambda2 mul t2 phi mul t2 0.080294 add y t2 mov t3 phi3 mul t3 0.0998909 add y t3 mov t4 lambda4 mul t4 phi mul t4 0.00199025 add y t4 mov t5 lambda2 mul t5 phi3 mul t5 -0.0285500 add y t5 mov t6 phi5 mul t6 -0.0490132 add y t6 ; ; 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