;
; KAVRAISKIY VI
; =============
;
; Name:           Kavraiskiy  VI. 
;                 Identisch mit dem Wagner I
;                 Libproj4 (2004), S. 55, 6.2.2
;                 Libproj4 (2006), S. 54, 5.2.4
; Quelle:         Gerald I. Evenden. Libproj4: A Comprehensive Library if Cartographic Projection Functions, March 2004
; Richtung:       Direkttransformation
;
; Maßstab, Berührungsparallel und das Längenverhältnis Pollinie-Äquator werden abgefragt.
; Die Zielbildgeometrie wird dem Sekundäroperanden oder einem Fixbild entnommen.
; Der Berührungspunkt wird bildmittig gesetzt.
;
; Identisch mit Wagner I, hier aber nicht nach Wagner, sondern nach Evenden implementiert.
;
; 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	Kavraiskiy~VI
	_var	phi		; Geographische Breite
	_var	lambda		; Geographische Länge
	_var	psi		; Parameter Psi
	_var	Cx
	_var	Cy
	_var	Cn
;
; Konstanten der Transformation
;
	_var	scale		; Kartenmaßstabszahl (also 1000000, nicht 1/1000000)
	_var	lambda0		; Geogr. Länge des Bildmittelpunktes
;
; x, y, x', y', Cx', Cy', °(, (°, pi, pi/2 etc. sind vordefinierte globale Konstanten
;
; Initialisierung
; ===============
;
	tstne	initial	077$
; Dialog
	input	scale	Maßstabszahl
	input	lambda0	Mittelpunktslänge~in~Grad
088$:
; Eingegebene Werte auf Min/Max bringen
	clip	scale	1	1E12
	clip	lambda0	-180	180
; Konstanten
	mov	Cx	3
	root	Cx	4
	mul	Cx	2
	div	Cx	3

	mov	Cy	3
	root	Cy	4

	mov	Cn	3
	root	Cn	2
	div	Cn	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
; -------------------
;
; x
;
	mov	psi	phi
	sin	psi
	mul	psi	Cn
	asin	psi

	mov	x	psi
	cos	x
	mul	x	lambda
	mul	x	Cx
;
; y
;
	mov	y	psi
	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
;	mov	x'	psi
;	mov	y'	i	
	exit
	_end