20060303-partial

1 files changed, 58 insertions, 0 deletions

diff --git a/module/math/geodesy.m b/module/math/geodesy.m
new file mode 100644
index 00000000..77ae1af0
--- /dev/null
+++ b/module/math/geodesy.m
@@ -0,0 +1,58 @@
+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;
+};