;
; ECKERT II
; =========
;
; Name:           Eckert II
; Proj4:	  +proj=eck2
;                 Libproj4, S. 63, 6.2.22
; 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 Mittelmeridian 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
; Gerald I. Evenden, USGS: Libproj4. A Comprehensive Library if Cartographic Projection Functions, March 2004

; (C) Rolf Böhm 2004

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

;
; Laufende Koordinaten
;
	_name	Eckert~II
	_var	phi		; Geographische Breite
	_var	lambda		; Geographische Länge
	_var	t1		; temp
	_var	t2		; temp		

;
; Konstanten der Transformation
;
	_var	lambda0		; Geogr. Länge des Bildmittelpunktes
	_var	scale		; Kartenmaßstabszahl (also 1000000, nicht 1/1000000)
	_var	r2pi3		; root(2*pi/3)
	_var	2r6pi		; 2/root(6*pi)
	_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	Eckert~II\\Maßstabszahl
	input	lambda0	Mittelpunktslänge~in~Grad
; Eingegebene Werte auf Min/Max bringen
	clip	scale	1	1E12
	clip	lambda0	-180	180
; Konstanten berechnen
	mov	r2pi3	pi
	mul	r2pi3	2
	div	r2pi3	3
	root	r2pi3	2

	mov	t1	pi	
	mul	t1	6
	root	t1	2
	mov	2r6pi	2
	div	2r6pi	t1
; 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	t2	y		; phi
	abs	t2
	div	t2	r2pi3
	mov	t1	2
	sub	t1	t2
	power	t1	2

	mov	phi	4
	sub	phi	t1
	div	phi	3
	asin	phi
	errjump	out
	abs	phi
	tstgt	y	m11$		; sign from y
	neg	phi
m11$:
	mov	lambda	x		; lambda
	mov	t2	phi	
	abs	t2
	sin	t2
	mul	t2	3
	mov	t1	4
	sub	t1	t2
	root	t1	2
	div	lambda	t1
	div	lambda	2r6pi
	
;
; In Gradmaß umrechnen und Ausserhalbtest
; ---------------------------------------
;
	mul	phi	(°
	mul	lambda	(°
	mul	theta	(°		; Test
	mov	t1	theta
	abs	t1
	cmpgt	t1	90	out	; Das verhindert Gegenwelten infolge Periodizität
	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