;
; URMAYEL FLAT-POLAR SINUSOIDAL
; =============================
;
; Name:           Urmayev flat-polar sinusoidal
; Proj4:	  +proj=urmfps
;                 Libproj4, S. 55, 6.2.2
; Quelle:         Gerald I. Evenden. Libproj4: A Comprehensive Library if Cartographic Projection Functions, March 2004
; Richtung:       Inverse Transformation
;
; 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, ausserdem ist das Verhältnis
; Mittelmeridian-Äquator einstellbar.
;
; Die Zielbildgeometrie wird dem Sekundäroperanden oder einem Fixbild entnommen.
; Der Berührungspunkt wird bildmittig gesetzt.

;
; Literatur:
; Wagner: Kartographische Netzentwürfe, Leipzig: Bibliographisches Institut 1949
; Gerald I. Evenden, USGS: Libproj4. A Comprehensive Library if Cartographic Projection Functions, March 2004

; (C) Rolf Böhm 2004

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

;
; Laufende Koordinaten
;
	_name	Urmayev~Flat-Polar~Sinusoidal
	_var	phi		; Geographische Breite
	_var	psi		; Parameter (Bei Evenden Theta)
	_var	lambda		; Geographische Länge
	_var	t1		; temp
	_var	t2		; temp	
;
; Konstanten der Transformation
;
	_var	lambda0		; Geogr. Länge des Bildmittelpunktes
	_var	scale		; Kartenmaßstabszahl (also 1000000, nicht 1/1000000)
	_var	ratio
	_var	m
	_var	n
	_var	initial
	_var	..........	; Ende Symboltabelle

;
; x, y, x', y', Cx', Cy', °(, (°, pi, pi/2 etc. sind vordefinierte globale Konstanten
;
; Initialisierung
; ===============
;
	tstne	initial	077$
; Dialog
	input	scale	Urmayev~Flat-Polar~Sinusoidal\\Maßstabszahl
	input	lambda0	Mittelpunktslänge~in~Grad
	input	ratio	Verhältnis~Mittelmeridian~-~Äquator\\0~=~Standard:\\0.5~(~=~Wagner~1/Kavrayskiy~VI)
	tstne	ratio	088$
	mov	ratio	0.5
088$:
; Eingegebene Werte auf Min/Max bringen
	clip	scale	1	1E12
	clip	lambda0	-180	180
	clip	ratio	0	0.9999
; Konstanten berechnen
	mov	n	1
	mov	t1	ratio
	power	t1	2
	sub	n	t1
	root	n	2

	mov	m	3
	root	m	4
	mul	m	2
	div	m	3
; Programm ist initialisiert
	mov	initial	1
077$:

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

;
; Maßstab, Kartenmittelpunkt etc. einrechnen
; ------------------------------------------

	sub	x	Cx'		; Bildmittelpunkt
	div	x	Rx'		; Erdradius
	mul	x	scale		; Kartenmaßstab

	sub	y	Cy'
	div	y	Ry'
	mul	y	scale
;
; Eigentlicher Entwurf, dieser invers
; -----------------------------------
;
	mov	psi	y		; Geographische Länge
	mul	psi	m
	mul	psi	n
	div	psi	1		; DRUCKFEHLER: HIER IRRT G I EVENDEN!

	mov	t1	psi
	abs	t1
	cmpgt	t1	pi/2	out	; Das verhindert Nebenwelten infolge Periodizität
	
	mov	phi	psi
	sin	phi
	div	phi	n
	asin	phi
	errjump	out

	mov	t1	psi
	cos	t1
	mul	t1	m
	mov	lambda	x
	div	lambda	t1
;
; In Gradmaß umrechnen und Ausserhalbtest
; ---------------------------------------
;
	mul	phi	(°
	mul	lambda	(°
	mul	psi	(°		; Test
	cmplt	phi	-90	out
	cmpgt	phi	90	out
	cmplt	lambda	-180	out
	cmpgt	lambda	180	out

;
; Schlussarbeiten
; ---------------
;
	mov	x'	lambda
	add	x'	lambda0
	mov	y'	phi
	mov	z'	psi		; Test


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

	exit
out:
	mov	x'	-9999
	mov	y'	-9999
	exit
	_end