Utility functions to convert terrestrial map coordinates from latitude-longitude pairs (in degrees) to eastings-northings pairs (in metres), and vice-versa. The routines can handle Transverse Mercator (TM) projections, including the Ordnance Survey of Great Britain (OSGB) projection and a number of Universal Transverse Mercator-type projections including the UTM projections based on the WGS84 and Hayford (International 1924) ellipsoids and the Gauss-Krüger projection using the Pulkovo 1942 ellipsoid. Other TM projections can be added to the code straightforwardly if the projection parameters are known (instructions are given within the code). The functions are based on the OSGB implementation of the Redfearn series and are accurate to around 1 metre. They are written in the GBU Octave language and should therefore also work correctly in MATLAB.

See: https://en.wikipedia.org/wiki/Transverse_Mercator:_Redfearn_series https://www.ordnancesurvey.co.uk/documents/resources/guide-coordinate-systems-great-britain.pdf