Back to Paul's main page.

Or Go to Paul's GPS links page.

In 1946 New Zealand introduced a Transverse Mercator mapping system based on the International Spheroid as recommended at the Madrid meeting of the International Geodetic and Geophysical Union in 1924.

The fundamental elements of that spheroid are

compression, f = 1/297

The projection was base on the Sears(1927) metre-foot relationship

( 12 inches = 1 foot, 3 feet = 1 yard, 1 link = 7.92 inches

and 100 links = 1 chain)

This all seems very odd today but New Zealand was surveyed in links , chains and acres.

A further complication arose in that New Zealand is a country with a comparatively small extent in longitude and more in latitude. It also consists basically of two large islands. To accommodate this factor the first New Zealand projection was split into two separate Transverse Mercator systems. These were the

and South Island National Yards Grid.

Both projections though were Transverse Mercator. However as time evolved it became obvious that there were some serious drawbacks.

- The North Island grid did not reach the Three Kings Islands
- The range of scale error was found to be unsatisfactory in Fiordland
- It was not possible to use the National Grids to provide a single unified map of New Zealand
- Computation of interisland air routes became vastly more difficult
- The requirements of off-shore prospecting created significant difficulties.

At the present the datum which should be selected is Geodetic Datum 49. This is equivalent to the spheroidal datum used for the New Zealand National Yards Grids. Unless a GPS unit has a built-in conversion it is not possible to output east and north coordinates in terms of the New Zealand Map Grid.

Plans are in train to establish a new datum for mapping in New Zealand.

See the links on my Surveying/Mapping for details. All the information on this page relates solely to Geodetic Datum 1949.

Latitudes and longitudes obtained by using Geodetic Datum 49 are plottable directly on to maps produced in NZ Map Grid but only if the graticule of lats/longs. is shown. Various approximations have been made to enable 'user grids' to be input in a GPS unit. There are all though, just that - approximations. I make no further comment on this.

Further information could probably be found from the popular

news group 'sci.geo.satellite-nav'

My solution has been to write my own moving map software (for PC laptop or desktop) which accepts scanned maps in NZ Map Grid. It converts NMEA information from the GPS to NZMG coordinates. The software though will only work in New Zealand. Anyone interested in the software can Email me.
Other software replicates the Garmin 12XL satellites in view/signal strength's screen and a lat/long. to NZMG converter.

It is not my intention to make these available on this page but anyone interested can Email me.

New Zealand lies between 166 degrees East and 180 degrees East in longitude and between 34 degrees South and 48 degrees South latitude.

To achieve the small range of scale variation over the land area of New Zealand it has been necessary to abandon the orderly arrangement of scale curves. The result is a unique projection known simply as the

- The true origin of the projection is placed at
**latitude 41 degrees South**, this being the latitude nearest to the middle latitude of the country. - The true origin is placed at
**longitude 173 degrees East**, this being the whole degree nearest to the meridian to which equal perpendiculars can be drawn from the easternmost and westernmost points of the country. - The meridian of longitude 173 degrees East, which is not represented by a straight line, is oriented so that its tangent at the origin is the north-south axis of coordinates.
- The true origin is assigned arbitrary coordinates sufficiently large to render all coordinates positive, or east and north of a so-called "false origin".
- In a metre coordinate system sufficiently extensive to cover the whole country, the northing must at some stage reach seven integral figures. In the New Zealand Map Grid the coordinates have seven integral figures in all cases.
- The Easting is always less than 5,000,000 metres, the Northing always greater than 5,000,000 metres, so that no confusion can arise between Easting and Northing whichever is stated first.
- The coordinates assigned to the true origin are sufficiently large for the grid to be extended a considerable distance out to sea without departing from the characteristics stated in (5) and (6).
- Thus, the true origin is place at latitude 41 degrees South, longitude 173 degrees East, and the coordinates of this point are

**2,510,000 metres East, 6,023150 metres North.**

The land area of New Zealand is thus fitted into a rectangle 1,000,000 metres wide by 1,500,000 metres from north to south. Coordinate values have been alloted to the nominated origin so that the easting ranges from 2,000,000 to 3,000,000 metres, and the Northing ranges from 5,300,000 to 6,800,000 metres. The seemingly strange values assigned to the origin have also been chosen to fit into a scheme of map sheets at a scale of 1:50,000 with boundaries at 10,000 metre values.

Latitudes and longitudes in other datums (unless already in WGS84) must first be converted to WGS84 datum and then further converted to Geodetic Datum 1949.

Transformation parameters: WGS84 - Geo49

diff Y = -22

diff Z = 209

Diff.major semi-axis(a) = -251

Diff.flattening * 10^4 =-0.14192702

Formulae for those conversions are not currently provided here.

- Both '*' and '.' are used to indicate multiplication (for clarity)
- 'd' is used to denote difference (or delta).
- 'u' is used to denote difference of isometric latitude.
- '^' is used to denote 'to the power of'
- 'i' is the imaginary part of a complex number
- Phi is the normally used symbol for latitude
- Lambda is the normally used symbol for longitude

d.u = 0.6399175073 * (d.phi * 10^-5) - 0.1358797613 * (d.phi * 10^-5)^2 +0.063294409 * (d.phi * 10^-5)^3 -0.025268530 * (d.phi * 10^-5)^4 +0.0117879 * (d.phi * 10^-5)^5 -0.0055161 * (d.phi * 10^-5)^6 +0.0026906 * (d.phi * 10^-5)^7 -0.001333 * (d.phi * 10^-5)^8 +0.00067 * (d.phi * 10^-5)^9 -0.00034 * (d.phi * 10^-5)^10These coefficients give 10-figure accuracy for the range of +/-7 degrees from 41 degrees south.

d.lambda being the difference of longitude

That is, coordinates from the true origin are given by

z = a(B1 * zeta + B2 * zeta^2 + B3 * zeta^3 + ...).

For the NZ Map Grid , the series is terminated after the sixth term, and constants are added to the coordinates, which are then given by

For the International Spheroid , the major semi-axis is given by

and the coefficients of the series defining the NZ Map Grid are as follows, all except B1 being complex numbers. i.e the right hand figures are the imaginery (i) parts.

B1 = 0.7557853228

B2 = +0.249204646 and + i(0.003371507)

B3 = -0.001541739 and + i(0.041058560)

B4 = -0.10162907 and + i(0.01727609)

B5 = -0.26623489 and + i(-0.36249218)

B6 = -0.6870983 and + i(-1.1651967)

then

The coefficients of the first six terms are

c2 = -0.577245789 and i(-0.007809598)

c3 = +0.508307513 and i(-0.1152208952)

c4 = -0.150947620 and i(+0.18200602)

c5 = +1.014181790 and i(+1.64497696)

c6 = +1.966054900 and i(+2.5127645)

zeta = __________________________________________________________________

B1 + 2B2.zeta + 3b3.zeta^2 + 4B4.zeta^3 + 5B5.zeta^4 + 6B6.zeta^5

A second application of this formula gives sufficient accuracy at any point within the land area of New Zealand.

When the final value of zeta is obtained, the latitude can be derived by from the following inverse series.

d(phi) * 10^-5 = 1.5627014243 * d.u +0.5185406398 * d.u^2 -0.03333098 * d.u^3 -0.1052906 * d.u^4 -0.0368594 * d.u^5 +0.007317 * d.u^6 +0.01220 * d.u^7 +0.00394 * d.u^8 -0.0013 * d.u^9

Scale coefficient and convergence formulae are not included here.

And, if I ever get around to it, some examples will be here.

Back to last page.

Back to Paul's main page.

Or Go to Paul's GPS links page.