;
; CRAIG RETROAZIMUTHAL
; ====================
;

; Name:           Craig Retroazimuthal
; Autor:          Craig
; Quelle:         Snyder, Voyland: U. S. G. S. Prof. Paper 1453 1989, S. 231.
; Richtung:       Direkttransformation
;
; Maßstab und Berührungsparallel werden abgefragt.
; Die Zielbildgeometrie wird dem Sekundäroperanden oder einem Fixbild entnommen.
; Der Berührungspunkt wird bildmittig gesetzt.
;
; Auflösung algebraisch nicht möglich. Berechnung mit Newton-Raphson-Iteration
;
; Anmerkung: Da der Bromley bei Evenden nicht abgebildet ist, war keine
; Verifikation möglich.
;
; Literatur:
; Wagner: Kartographische Netzentwürfe, Leipzig: Bibliographisches Institut 1949
; Fiala: Mathematische Kartographie, Berlin: Verlag Technik 1957

; (C) Rolf Böhm 2004

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

;
; Laufende Koordinaten
;
	_name	Craig~Retroazimuthal
	_var	phi		; Geographische Breite
	_var	lambda		; Geographische Länge

;
; Konstanten der Transformation
;
	_var	scale		; Kartenmaßstabszahl (also 1000000, nicht 1/1000000)
	_var	lambda0		; Geogr. Länge des Bildmittelpunktes
	_var	lambda1		; desgl. in Arc		
	_var	lambda1/2	; / 2	 
	_var	phi0		; Geogr. Breite des Bildmittelpuntkes
	_var	phi1		; desgl. in Arc
	_var	tanp1		; tan(Phi1)
;
; 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	phi0	Mittelpunktsbreite~in~Grad
	input	lambda0	Mittelpunktslänge~in~Grad
; Eingegebene Werte auf Min/Max bringen
	clip	scale	1	1E12
	clip	phi0	-90	90	
	mov	phi1	phi0
	mul	phi1	°(
	mov	tanp1	phi1
	tan	tanp1
	clip	lambda0	-180	180
	mov	lambda1	lambda0
	mul	lambda1	°(
	mul	lambda1/2	lambda2
	div	lambda1/2	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

;
; Ausserhalbtest
; --------------
;
	mov	r0	lambda
	sub	r0	lambda0
	abs	r0
	cmpgt	r0	90	$out

	mov	r0	phi
	sub	r0	phi0
	abs	r0
	cmpgt	r0	90	$out

; 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	r0	lambda
	sub	r0	lambda1 ;/2

	mov	x	r0

	mov	r1	r0
	cos	r1	
	mov	r2	phi
	sin	r2
	mul	r1	r2
	
	mov	r3	phi
	cos	r3
	mul	r3	tanp1

	sub	r1	r3

	mov	r4	r0	; "(lambda-lambda0)"
	sin	r4
	tstne	r4	$777
	mov	r4	0.000001	
$777:
	div	r1	r4

	mul	r1	r0
	mov	y	r1

; 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
$out:
	mov	x'	-9999	
	mov	y'	-9999
	exit
	_end