Geo::HelmertTransform (3pm)

1.14

    use Geo::HelmertTransform;

    my ($lat, $lon, $h) = ...; # from OS map
    my $airy1830 = Geo::HelmertTransform::datum('Airy1830');
    my $wgs84 = Geo::HelmertTransform::datum('WGS84');

    ($lat, $lon, $h)
        = Geo::HelmertTransform::convert_datum($airy1830, $wgs84,
                                                $lat, $lon, $h);

Perform transformations between geographical coordinates in different datums.

It is usual to describe geographical points in terms of their polar coordinates (latitude, longitude and altitude) referenced to a \*(L"datum ellipsoid\*(R", which is used to approximate the Earth's geoid. The latitude, longitude and altitude of a given physical point vary depending on which datum ellipsoid is in use. Unfortunately, a number of ellipsoids are in everyday use, and so it is often necessary to transform geographical coordinates between different datum ellipsoids.

Two different datum ellipsoids may differ in the locations of their centers, or in their shape; and there may be an angle between their equatorial planes or the meridians relative to which longitude is measured. The Helmert Transform, which this module implements, is a linear transformation of coordinates between pairs of datum ellipsoids in the limit of small angles of deviation between them.

Latitude is expressed in degrees, positive-north; longitude in degrees, positive-east. Heights (ellipsoid) and cartesian coordinates are in meters.

rad_to_deg \s-1RADIANS\s0

Convert \s-1RADIANS\s0 to degrees.

deg_to_rad \s-1DEGREES\s0

Convert \s-1DEGREES\s0 to radians.

geo_to_xyz \s-1DATUM\s0 \s-1LAT\s0 \s-1LON\s0 H

Return the Cartesian (X, Y, Z) coordinates for the geographical coordinates (\s-1LAT\s0, \s-1LON\s0, H) in the given \s-1DATUM\s0.

xyz_to_geo \s-1DATUM\s0 X Y Z

Return the geographical (\s-1LAT\s0, \s-1LON\s0, H) coordinates for the Cartesian coordinates (X, Y, Z) in the given \s-1DATUM\s0. This is an iterative procedure.

convert_datum D1 D2 \s-1LAT\s0 \s-1LON\s0 H

Given geographical coordinates (\s-1LAT\s0, \s-1LON\s0, H) in datum D1, return the corresponding coordinates in datum D2. This assumes that the transformations are small, and always converts via \s-1WGS84\s0.

datum \s-1NAME\s0

Return the datum of the given \s-1NAME\s0. Currently implemented are:

Airy1830

The 1830 Airy ellipsoid to which the British Ordnance Survey's National Grid is referenced.

Airy1830Modified

The modified 1830 Airy ellipsoid to which the Irish Grid (as used by Ordnance Survey Ireland and Ordnance Survey Northern Ireland); also known as the Ireland 1975 datum.

\s-1WGS84\s0

The global datum used for \s-1GPS\s0.

A guide to coordinate systems in Great Britain, http://www.gps.gov.uk/guidecontents.asp

Making maps compatible with \s-1GPS\s0, http://www.osni.gov.uk/downloads/Making%20maps%20GPS%20compatible.pdf

Written by Chris Lightfoot, [email protected]

Copyright (c) \s-1UK\s0 Citizens Online Democracy.

This module is free software; you can redistribute it and/or modify it under the same terms as Perl itself.