;
; ECKERT-GREIFFENDORFS POLYKONISCHE PROJEKTION
; ============================================
;

; Name:           Eckert~Greiffendorfs polykonische Projektion (Eckert VII von 1935)
; Autor:          Max Eckert-Greifendorf (1935)
; Quelle:         Canters & Decleire 1989, S. 65
; Richtung:       Direkttransformation
;
; Der Maßstab und die Mittelmeridianlänge werden abgefragt.
; Die Zielbildgeometrie wird dem Sekundäroperanden oder einem Fixbild entnommen.
; Der Berührungspunkt wird bildmittig gesetzt.
;
; Literatur:
; Canters & Decleire 1989


; (C) Rolf Böhm 2006

; Benutzte Variablen
; ==================
;

;
; Laufende Koordinaten
;
	_name	Eckert-Greiffendorfs~polykonische~Projektion~(Eckert~VII)
	_var	phi		; Geographische Breite
	_var	lambda		; Geographische Länge
	_var	lambda4		; Lambda/4
	_var 	root2		; Root(2)

;
; 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	root2	2
	root	root2	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
; -------------------
;

; Lambda/4

	mov	lambda4	lambda
	div 	lambda4	4

; x
	
	mov	r1	phi
	cos	r1
	mov	r2	lambda4
	sin	r2
	mul	r1	r2

	mov	r3	phi
	cos	r3
	mov	r4	lambda4
	cos	r4	
	mul	r3	r4
	add	r3	1
	root	r3	2

	mov	x	r1
	div	x	r3
	mul	x	root2
	mul	x	4	

; y

	mov	r1	phi
	sin	r1

	mov	y	r1
	div	y	r3
	mul	y	root2
;
; 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