|a||Semi-major axis of reference ellipsoid|
||Latitude of computation point|
||Longitude of computation point|
|N||Northing of computation point|
|E||Easting of computation point|
The following equations are divided into three sections:
The conversion of NZGD1949 geographic coordinates (θ, λ) to NZMG projection coordinates ( N , E ) is achieved in several steps.
The conversion requires the latitude and longitude expressed in decimal degrees (north and east are positive). The constants are provided on the NZMG page and coefficients used in these formulae are listed below.
Calculate ΔΨ and Δλ using:
Calculate a complex polynomial function of the complex number θ = ΔΨ + iΔλ using the formula:
Express the complex number z in terms of its real and imaginary parts z = x + iy and calculate the northing N and easting E from x and y as:
The conversion from NZMG to latitude and longitude involves an iterative approximation to reverse step 2 above. This starts with the easting E and northing N. The steps to calculate latitude and longitude are as follows. Note that the constants and coefficients used in these formulae are listed below. Derive the complex number z as the following steps.
Determine the complex number θ as a series of approximations θ0 , θ1 , θ2 and so on. The first approximation is:
Successive approximations are obtained by applying the formula:
Two iterations of this formula will give millimetre transformation accuracy.
Express θ in terms of its real and imaginary parts θ = ΔΨ + iΔλ and calculate the latitude and longitude as follows
|Coefficient||Real Part||Imaginary Part|