;
; DER VON CANTERS OPTIMIERTE WAGNER IX (WAGNER-AÏTOFF)
; ====================================================
;
; Kurzname:       Canters Wagner IX
; Originalname:   Optimised Aitoff-Wagner
; Autor:          Karlheinz Wagner 1949, Frank Canters 2002
; Quelle:         Karlheinz Wagner, Kartographische Netzentwürfe, Leipzig 1949, S. 205ff.
;                 Frank Canters, small-scale Map-Projection Design 
; Richtung:       Direkt-Transformation
;
; Das Programm übernimmt die Koordinaten eines Punktes (x/y) und transformiert
; diese in einen Punkt (x'/y').
;
; Literatur:
; Wagner: Kartographische Netzentwürfe, Leipzig: Bibliographisches Institut 1949
; Fiala: Mathematische Kartographie, Berlin: Verlag Technik 1957
; Canters: Small-scale Map Projection Design, London: Taylor & Francis 2002

; (C) Rolf Böhm 2006

; Benutzte Variablen
; ==================
;
; Die Variablennamen entsprechen weitgehend denen von Karlheinz Wagner,

;
; Laufende Koordinaten
;
	_name	Canters~Wagner~IX~(Aïtoff-Wagner)-Optimierung
	_var	phi		; Geographische Breite
	_var	lambda		; Geographische Länge
	_var	psi		; Parameter
	_var	delta
	_var	alpha
	_var	n
	_var	m
	_var	k1
	_var	k2
	_var	Cx
	_var	Cy
;
; 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
; Konstanten
	mov	n	0.4809
	mov	m	0.7459
	mov	k2	1.2664
	mov	k1	1.0226

	mov	r0	m
	mul	r0	n
	root	r0	2
	
	mov	Cx	k1
	div	Cx	r0

	mov	Cy	1
	div	Cy	k2
	div	Cy	r0

; Programm ist initialisiert
	mov	initial	1
077$:

;
; SIMD-Laufbereich
; ================

;
; Vorbereitung
; ------------
;
	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
; -------------------
;
	mov	psi	phi		; psi
	mul	psi	m

	mov	r0	psi		; delta
	cos	r0
	mov	delta	lambda
	mul	delta	n
	cos	delta
	mul	delta	r0
	abs	delta
	acos	delta

	mov	alpha	psi		; alpha
	sin	alpha
	mov	r0	delta
	sin	r0
	div	alpha	r0
	acos	alpha

	mov	x	alpha		; x
	sin	x
	mul	x	delta
	mul	x	Cx

	mov	r7	lambda		; sign of x
	sgn	r7
	mul	x	r7

	mov	y	alpha
	cos	y
	mul	y	delta
	mul	y	Cy
;
; 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