;
; WIECHELS PSEUDOAZIMUTALE PROJEKTION
; ===================================
;
; Name:           Wiechels pseudoazimutale Projektion
; Autor:          H. Wiechel, 1879      
; Quelle:         Bugayevskij, Snyder: Map Projections. London: Taylor 1998
; Richtung:       Direkt-Transformation
;
; Flächentreu. Wohl von akademischer Bedeutung.
;
; 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:
; Wiechel, H.: Rationelle Gradnetzprojektionen, Civilingenieur, new series, 25. Jg., 
; cols 401-22 (1879). 
; Bugayevskij, Snyder: Map Projections. London: Taylor 1998
;

; (C) Rolf Böhm 2007

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

;
; Laufende Koordinaten
;
	_name	Wiechels~Projektion
	_var	phi		; Geographische Breite
	_var	lambda		; Geographische Länge
;
; Konstanten der Transformation
;
	_var	lambda0		; Geogr. Länge des Bildmittelpunktes
	_var	scale		; Kartenmaßstabszahl (also 1000000, nicht 1/1000000)
; Sonstiges
	_var	sinphi
	_var	cosphi
	_var	sinlambda
	_var	coslambda		
;
; 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
; ================

;
; Eingabekoordinaten übernehmen
; -----------------------------
;
; 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
;
	mov	lambda	x		; Geographische Länge
	mov	phi	y		; Geographische Breite
;
; Um lamda0 verdrehen
; -------------------
;
	sub	lambda	lambda0
	cmod	lambda	-180	180
;
; Umrechnung in Bogenmaß
; ----------------------
;
	mul	phi	°(
	mul	lambda	°(
;
; Netz rechnen
; ------------
;
	mov	coslambda	lambda
	cos	coslambda
	mov	sinlambda	lambda
	sin	sinlambda
	mov	cosphi		phi
	cos	cosphi
	mov	sinphi		phi
	sin	sinphi
;
	mov	r1	cosphi
	mul	r1	sinlambda
;
	mov	r2	1
	sub	r2	sinphi
	mul	r2	coslambda
;
	mov	r3	cosphi
	mul	r3	coslambda
;
	mov	r4	1
	sub	r4	sinphi
	mul	r4	sinlambda
;
	mov	x	r1	
	sub	x	r2
;
	mov	y	r3
	add	y	r4
	neg	y
	



;
; 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