;
; CANTERS LOW ERROR POLYCONIC
; ===========================
;
; Name:           Canters Low-error polyconic projection with twofold symmetry, equally spaced
;                 parallels and a correct ratio of the axis
; Quelle:         Eric M. Delmelle. Map Projections Properties. Considerations for
;                 Small-Scale GIS Applications. State Univ. of New York at Buffalo,
;                 Dept. of Geography, September 2001
; Richtung:       Direkt-Transformation
;
; Anmerkung:      Bei Delmelle als Canters Minimum Error bezeichnet. Den angegebenen Namen habe ich
;                 von F. Canters direkt. --  Bei Evenden gibt es weitere 4 Pseudocylindric Canters Low Error Projections.
;
; 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.

; (C) ROlf Böhm 2004

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

;
; Laufende Koordinaten
;
	_name	Canters~Low~Error~Polyconic
	_var	phi		; Geographische Breite
	_var	lambda		; Geographische Länge
	_var	phi2		; Quadrate, Kuben ...
	_var	phi3
	_var	phi4
	_var	lambda2
	_var	lambda3
	_var	lambda4
;
; 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.
	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
; ================

;
; 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 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	lambda2	lambda
	mov	lambda3	lambda
	mov	lambda4	lambda
	power	lambda2	2
	power	lambda3	3
	power	lmabda4	4
	mov	phi2	phi
	mov	phi3	phi
	mov	phi4	phi
	power	phi2	2
	power	phi3	3
	power	phi4	4
;
	clr	x
	mov	t1	lambda	
	mul	t1	0.9305	
	add	x	t1

	mov	t2	lambda
	mul	t2	phi2
	mul	t2	-0.1968
	add	x	t2

	mov	t3	lambda3
	mul	t3	phi2
	mul	t3	-0.0067
	add	x	t3	

	mov	t4	lambda
	mul	t4	phi4
	mul	t4	0.0076
	add	x	t4	
;
	clr	y
	mov	t1	phi
	mul	t1	0.9305
	add	y	t1

	mov	t2	lambda2
	mul	t2	phi
	mul	t2	0.0394
	add	y	t2

	mov	t3	lambda4
	mul	t3	phi
	mul	t3	0.0005
	add	y	t3

	mov	t4	lambda2	
	mul	t4	phi3
	mul	t4	-0.0115
	add	y	t4
;
; 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