;
; FURUTIS VERALLGEMEINERTER ABSTANDSGLEICHER POLYKONISCHER ENTWURF
; ================================================================
;
; Name:           Generalized Transversal Equidistant Polyconic Projection
; Autor:          Carlos A. Furuti
; Richtung:       Direkt-Transformation
;
; Anmerkung:      Furuti rechnet eine transversale Erde mit Metapol 0° N, 95° E und 
;                 95° Verdrehung nach der Metatransformation. Dies ist vor der Projektion
;                 mit einer sphärischen Transformation Kugel auf Kugel auszuführen.
;                 Anschließend erst Projektion mit diesem Programm, wobei lambda0 auf -25 und
;                 lambdaAS auf -30 zu setzen ist. Natürlich lassen sich auch beliebig andere 
;                 Parametrierungen rechnen.
;
; Note:		  The Original Furuti Transformation uses a transversal Earth sphere with
;                 meta-pole 0° N, 95° E and 95° post-sherical-transformation rotation. This must be
;                 done before run of this projection program. In these programm are the Original Furuti
;                 parameters zero-longitude lambda0 -25 and longitude assymetrical shift -30.
;        	  But any other parametrizations are also possible.
;
; (C) Rolf Böhm 2004, 2008

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

;
; Laufende Koordinaten
;
	_name	Äquidistanter~polykonischer~Entwurf,~Furutis~Verallgemeinerung
	_var	phi		; Geographische Breite
	_var	lambda		; Geographische Länge

	_var	r		; Polarkoordinaten: Radiua
	_var	alpha		; Polarkoordinaten: Azimut
	_var	c		; Der Polparameter c der polykonischen Entwürfe

	_var	t1	
	_var	t2
	_var	t3
	_var	lambdaAS	;Lambda Assymetrical Shift
;
; 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.\\¶
		Die~Original~Furuti-Projektion~muss~mit~einer~transversalen~Erdlage,~Metapol~¶
		0°~N,~95°~E~und~95°~Längenverdrehung~(Postrotation)~gerechnet~werden.
	input	scale	Maßstabszahl
	mov	lambda0	-25
	input	lambda0	Mittelpunktslänge~in~Grad~(Original-Furuti:~-25) 
	mov	lambdaAS	-30
	input	lambdaAS	Assymetrischer~Längendefekt~in~Grad~(Original-Furuti:~-30)
; Eingegebene Werte auf Min/Max bringen
	clip	scale	1	1E12
	clip	lambda0	-180	180
	clip	lambdaAS	-180	180
	mul	lambdaAS	°(
; Programm ist initialisiert
	mov	initial	1
077$:

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

;
; Eigentlicher Entwurf, dieser direkt!
; ------------------------------------
;
	mov	lambda	x		; Geographische Länge
	mov	phi	y		; Geographische Breite

	sub	lambda	lambda0	

	cmpgt	lambda	-180	10$
	add	lambda	360
10$:
	cmpgt	lambda	-180	30$
	add	lambda	360
30$:
	cmplt	lambda	180	40$
	sub	lambda	360
40$:
	cmplt	lambda	180	50$
	sub	lambda	360
50$:

; x, y haben eine Doppelbedeutung:

; - einmal die vordefinierten 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
; -------------------
					; Nur dieser eine Befehl ist Furutis Verallgemeinerung! 
	add	lambda	lambdaAS	; Only that instruction is the whole Furuti generalization! 

	mov	r	phi
	cot	r

	mov	alpha	phi
	
	sin	alpha
	mul	alpha	lambda
	

	mov	c	phi
	add	c	r

	mov	t1	alpha
	sin	t1
	mul	t1	r
	mov	x	t1

	mov	t2	alpha
	cos	t2
	mul	t2	r
	
	mov	y	c
	sub	y	t2
;	mov	x'	alpha
;	mov	y'	r
;	exit
		
;
; 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