code: purgatorio

ref: 611dced75f0a6c9fd4b35b88ee0dd9ac5806cb54
dir: /module/math/geodesy.m/

Geodesy: module
{
	PATH: con "/dis/math/geodesy.dis";

	# easting, northing in metres
	Eano: adt{
		e: real;
		n: real;
	};

	# latitude, longitude in radians
	Lalo: adt{
		la: real;
		lo: real;
	};

	# datums
	# WGS84 and ITRS2000 effectively the same
	OSGB36, Ireland65, ED50, WGS84, ITRS2000, ETRS89: con iota;

	# transverse Mercator projections
	Natgrid, IrishNatgrid, UTMEur, UTM: con iota;

	# call first
	# d specifies the datum (default WGS84)
	# t specifies the transverse Mercator projection (default Natgrid)
	# z specifies the UTM zone if relevant (default 30)
	# calls format below
	init: fn(d: int, t: int, z: int);

	# alters the current datum, transverse Mercator projection and UTM zone
	# use a negative value to leave unaltered
	format: fn(d: int, t: int, z: int);

	# OS string to (easting, northing) and back
	# formats XYen, XYeenn, XYeeennn, XYeeeennnn, XYeeeeennnnn or
	# formats eenn, eeennn, eeeennnn, eeeeennnnn, eeeeeennnnnn
	os2en: fn(s: string): (int, Eano);	# returns (0, ...) if bad string format
	en2os: fn(en: Eano): string;

	# latitude/longitude string to (latitude, longitude) and back
	# format latitude longitude
	# formats deg[N|S], deg:min[N|S], deg:min:sec[N|S] for latitude
	# formats deg[E|W], deg:min[E|W], deg:min:sec[E|W] for longitude
	str2lalo: fn(s: string): (int, Lalo);	# returns (0, ...) if bad string format
	lalo2str: fn(lalo: Lalo): string;

	# general string to (easting, northing)
	# OS grid or latitude/longitude format as above
	str2en: fn(s: string): (int, Eano);	# returns (0, ...) if bad string format

	# (easting, northing) to (latitude, longitude) and back
	en2lalo: fn(en: Eano): Lalo;
	lalo2en: fn(lalo: Lalo): Eano;

	# approximate transformations between any of OSGB36, WGS84, ITRS2000, ETRS89
	datum2datum: fn(lalo: Lalo, f: int, t: int): Lalo;
};