Projection and Transformation Calculations
Projection and Transformation Calculations
Projection and Transformation Calculations
33 23-May-2007
This spreadsheet which will carry out common calculations with coordinates. All of the concepts and formulae given in "A guide to co
the user manual for these calculations, are detailed in this spreadsheet. "A guide to coordinate systems in Great Britain" is available fro
(www.ordnancesurvey.co.uk/gps). The guide gives some important warnings about the use of its formulae, which you should be aware
such as scale factors, convergence and t-T in this spreadsheet that are not in the guide.
The Constants sheet is where ellipsoid, projection and Helmert transformation parameters are entered.
The Enter coordinates here sheet is for user input of coordinates to the calculations.
Use the buttons at the bottom left of the screen to scroll through the other sheets.
Each sheet carries out a different calculation. The calculations are shown in "step by step" format where all the intermediate calculatio
results in single cells computed by user defined functions written in Visual Basic. The code for the user defined functions can be view
applications.
Cells for user data entry are coloured Cells giving calculation results are coloured Cells showing intermediate c
Cells which repeat the user inputted data are coloured
and formulae given in "A guide to coordinate systems in Great Britain", which is
ems in Great Britain" is available from the Ordnance Survey GPS web site
ormulae, which you should be aware of. There are also other projection functions
ered.
where all the intermediate calculation results are given and they are also shown as
e user defined functions can be viewed and if required copied to other
84 geodetic coordinates to OSGB36 Eastings and Northings and vice versa can be
e user.
by the user.
er and also the line scale factor between easting, northing and easting2, northing2.
ween easting, northing and easting2, northing2 entered by the user.
ween easting, northing and easting2, northing2 entered by the user.
° ' "
Semi-major axis, a 6377563.396 True origin latitude, j0 N 49 0 0
Semi-minor axis, b 6356256.910 True origin longitude, l0 W 2 0 0
Central Meridan Scale, F0 0.999601271700 a for OSGB36 =
True origin Easting, E0 400000.000 b for OSGB36 =
True origin Northing, N0 -100000.000 a for GRS80 & WGS84 =
aF0 6375020.48098897 b for GRS80 =
bF0 6353722.49048791 b for WGS84 =
All other parameters same
Projection Notes
1. The constants above are correct for working with the Ordnance Survey National Grid.
If you want to work with OS National Grid, do not change anything in the above box.
2. It is also possible to use this spreadsheet to work with Universal Transverse Mercator (UTM) grids.
The constants that must be entered in the yellow boxes if you want to work with UTM are given in Annex A of
'A Guide to coordinate systems in Great Britain', available from the OS GPS website www.gps.gov.uk.
Transformation Notes
1. It is very important to understand the limitations of simple datum transformations.
This method does not take account of the local distortions present in coordinate reference frames (TRFs)
This is especially important when one or both of the coordinate systems involved is based on a traditional triangulation net
For instance, transforming between WGS84 GPS coordinates and OSGB36 National Grid will incur transformation errors
of up to five metres, depending on location.
Also, this method cannot convert ellipsoid heights to heights above mean sea level.
For full details on the limitations of the method coded here, see the OS publication 'A guide to coordinate systems in
Great Britain', available free from our GPS Website (www.gps.gov.uk) or OS customer helpline.
2. The constants in the table above are suitable for transforming from WGS84 GPS coordinates to OSGB36 National Grid
By changing the constants above, this spreadsheet may be used to convert between any two geodetic datums
Some ellipsoid constants are given in 'A guide to coordinate systems in Great Britain'
Page 3
Constants
Ordnance Survey cannot supply suitable transformation parameters for the geodetic systems of other countries -
please make enquiries for these parameters to the national mapping agencies of the country concerned.
Page 4
Constants
n Annex A of
ordinate systems in
Page 5
Constants
her countries -
Page 6
Enter coordinates here
Notes
1. Type your input data in the yellow boxes only on this page.
Enter either Latitude, Longitude and optional Ellipsoidal height, Easting and Northing or Earth centred Cartesian coordina
To convert from decimal latitude and longitude or degrees lat / long and decimal minutes to the degrees, minutes and secon
see the "Lat & Long Format Conversions " sheet.
Ellipsoidal height is only required for the Cartesian XYZ coordiante conversions.
The Earth centred Cartesian coordinates are used as input to the Helmert datum transformation calculations.
Look at the following sheets to see the results of the calculations.
2. Check that the appropriate parameters are set on the "Constants" sheet.
3. For two of the calculations (t-T correction and true azimuth), the eastings and northings of two points are required.
For these calculations only, you must enter Eastings 2 and Northings 2 in addition to Easting and Northing
4. For background information on the calculations contained in this spreadsheet, download 'A Guide to Coordinate System
in Great Britain' from our GPS website www.gps.gov.uk, or see any textbook on map projections or surveying computatio
Page 7
Enter coordinates here
ation calculations.
Page 8
Latitude & Longitude Format Conversions (Enter your coordinates in the yellow boxes)
Degrees, Minutes & Seconds to Decimal Degrees and Degrees & Decimal Minutes
Input Data Output Data
Degrees (°) Minutes (') Seconds (") Decimal Degrees Degrees & Decim
latitude N 52 39 27.253100 52.657570305556 N
longitude E 1 43 4.517700 1.717921583333 E
Decimal Degrees to Degrees, Minutes & Seconds and Degrees & Decimal Minutes
Output Data Input Data Output D
Degrees (°) Minutes (') Seconds (") Decimal Degrees Degrees & Decim
latitude N 52 39 27.253100 52.657570305556 N
longitude E 1 43 4.517700 1.717921583333 E
Degrees & Decimal Minutes to Degrees, Minutes & Seconds and Decimal Degrees
Output Data
Degrees (°) Minutes (') Seconds (") Decimal Degrees Degrees & Decim
latitude N 52 39 27.253100 52.657570305556 N
longitude E 1 43 4.517700 1.717921583333 E
n the yellow boxes)
ecimal Minutes
Output Data
Degrees & Decimal Minutes
52 39.454218333
1 43.075295000
ecimal Minutes
Output Data
Degrees & Decimal Minutes
52 39.454218333
1 43.075295000
Decimal Degrees
Input Data
Degrees & Decimal Minutes
52 39.454218333
1 43.075295000
Helmert Datum Transformation
X 3909833.018
Y -147097.1376
Z 5020322.4777
° ' " Decimal Deg Radians
latitude #VALUE! ### ### #VALUE! #VALUE! #VALUE!
longitude W 2 9 16.510138 -2.154586149489 -0.0376046223
Ellipsoidal height #VALUE! m
p 3912599.110
e2 6.6705397616E-03
Initial j N 52 15 16.6721 52.2546311476693 0.9120153629
n #VALUE!
Page 13
XYZ to lat,long,H
Page 14
lat,long to E,N
Easting #VALUE! m This is the result of the calculation using user defined f
Northing #VALUE! m To see the Visual Basic code of the functions -
Menu = Tools, Macro, Visual Basic Editor,
then view the code in "Module 1".
Page 15
lat,long to E,N
Page 16
E,N to lat,long
latitude #VALUE! ### ### #VALUE! #VALUE! This is the result of the ca
longitude #VALUE! ### ### #VALUE! #VALUE! To see the Visual Basic c
Menu = Tools,
then view the code in "M
Page 17
E,N to lat,long
The radians value of j' and the value of M are computed using user defined functions
To see the Visual Basic code of the functions -
Menu = Tools, Macro, Visual Basic Editor,
then view the code in "Module 1".
Page 18
Latitude, longitude and Ellipsoidal height to Cartesian XYZ
WGS 84 to OSGB36
WGS84 Latitude, longitude and Ellipsoidal height to WGS84 Cartesian XYZ
° ' " Dec Degs
latitude N 52 39 28.723003 52.65797861194 Input
longitude E 1 42 57.787253 1.71605201472
Ellipsoidal height 69.391 m
X #VALUE! m Output
Y #VALUE! m
Z #VALUE! m
X #VALUE! m Input
Y #VALUE! m See the Transformation Notes on the
Z #VALUE! m "Constants" page for information on the
X #VALUE! m Output limitations and accuracy of this
Y #VALUE! m transformation.
Z #VALUE! m
OSGB36 Cartesian XYZ to OSGB36 Latitude, longitude and approx ODN height
X #VALUE! m Input
Y #VALUE! m
Z #VALUE! m
° ' " Dec Degs
latitude #VALUE! ### ### #VALUE! #VALUE! Output
longitude #VALUE! ### ### #VALUE! #VALUE!
~ODN height #VALUE! m
OSGB36 to WGS84
OSGB36 easting and northing to OSGB36 latitude and longitude
OSGB36 latitude, longitude and approx ODN height to OSGB36 cartesian XYZ
° ' " Dec Degs
latitude #VALUE! ### ### #VALUE! #VALUE! Input
longitude #VALUE! ### ### #VALUE! #VALUE!
~ODN height 24.700 m
X #VALUE! m Output
Y #VALUE! m
Z #VALUE! m
X #VALUE! m Input
Y #VALUE! m See the Transformation Notes on the
Z #VALUE! m "Constants" page for information on the
X #VALUE! m Output limitations and accuracy of this
Y #VALUE! m transformation.
Z #VALUE! m
X #VALUE! m Input
Y #VALUE! m
Z #VALUE! m
° ' " Dec Degs
latitude #VALUE! ### ### #VALUE! #VALUE! Output
longitude #VALUE! ### ### #VALUE! #VALUE!
Ellipsoidal height #VALUE! m
sed.
view the code in "Module 1".
d in stages to an OSGB36 easting, northing and approx ODN height.
B36
ght to WGS84 Cartesian XYZ
User inputs DMS and Height. Decimal degrees and
height are the input to the user defined functions.
The DMS values are converted to decimal using
standard Excel functions
Ellipsoid parameters (WGS84)
a= 6378137.000
b= 6356752.314
GS84 to OSGB36)
Transformation (WGS84 to OSGB36)
translation parallel to X = -446.448
translation parallel to Y = 125.157
translation parallel to Z = -542.060
scale change = 20.4894
rotation about X = -0.1502
rotation about Y = -0.2470
rotation about Z = -0.8421
S84
6 latitude and longitude
The eastings, northings and approx ODN height are
input by the user. User defined functions convert
them to dec degrees. The decimal degrees are
converted to DMS using standard Excel functions
Ellipsoid and projection parameters (OSGB36)
a= 6377563.396
b= 6356256.910
f0 = 0.9996012717
e0 = 400000
n0 = -100000
j0 = 49.00000
l0 = -2.00000
SGB36 to WGS84)
Transformation (OSGB36 to WGS84)
translation parallel to X = 446.448
translation parallel to Y = -125.157
translation parallel to Z = 542.060
scale change = -20.4894
rotation about X = 0.1502
rotation about Y = 0.2470
rotation about Z = 0.8421
convergence #VALUE! ### #VALUE! #VALUE! This is the result of the calc
To see the Visual Basic cod
Menu = Tools,
then view the code in "Mod
Notes
Grid convergence is the angle between the north-south grid line and the direction of the meridian.
Page 27
convergence from lat & long
Page 28
convergence from E,N
Easting 651409.903
Northing 313177.270
° ' " Decimal Deg Radians
convergence #VALUE! ### #VALUE! #VALUE! #VALUE!
j' #VALUE! ### #VALUE! #VALUE! #VALUE!
n #VALUE!
r #VALUE!
h2 #VALUE!
Et 2.51409903E+05
XVI #VALUE!
XVII #VALUE!
XVIII #VALUE!
Notes
Grid convergence is the angle between the north-south grid line and the direction of the meridian.
Page 29
convergence from E,N
n of the meridian.
Page 30
local scale factor from lat,lon
factor #VALUE! This is the result of the calculation using a user defined functio
To see the Visual Basic code of the function -
Menu = Tools, Macro, Visual Basic Editor,
then view the code in "Module 1".
Notes
Local scale factor is the scale distortion inherent in the map projection at a point
Page 31
local scale factor from lat,lon
lculation steps
Page 32
local & line sf from E,N
Page 33
local & line sf from E,N
ength this will give a result with an error not exceeding 1 to 2 ppm.
int 2 (East2,North2)
1) and Point 2 (East2,North2)
Page 34
t-T correction from E,N
Easting 651409.903
Northing 313177.270
Easting 2 626238.249
Northing 2 302646.415
° ' " Decimal Deg Radians
(t1 - T1) #VALUE!### #VALUE! #VALUE! #VALUE!
(t2 - T2) #VALUE!### #VALUE! #VALUE! #VALUE!
Nm 307911.8425
M #VALUE!
j' #VALUE!### #VALUE! #VALUE! #VALUE!
n #VALUE!
r #VALUE!
XXIII #VALUE!
Notes
t-T correction is the difference between true direction and projected direction due to the curving of straight lines on the pro
Page 35
t-T correction from E,N
The radians value of j' and the value of M are computed using user defined functions
To see the Visual Basic code of the functions -
Menu = Tools, Macro, Visual Basic Editor,
then view the code in "Module 1".
Page 36
true azimuth from E,N
Easting 651409.903
Northing 313177.270
Easting 2 626238.249
Northing 2 302646.415
° ' " Decimal Deg Radians
True Azimuth(1 to 2) #VALUE! ### #VALUE! #VALUE! #VALUE!
Grid Bearing(1 to 2) 247 17 50.766311 247.29743508633 4.316154474
C1 #VALUE! ### #VALUE! #VALUE! #VALUE!
(t1 - T1) #VALUE! ### #VALUE! #VALUE! #VALUE!
Notes
The true azimuth from one point (easting, northing) to a second point (easting 2, northing 2)
is computed by applying the convergence and t-T corrections to the grid bearing.
Page 37
true azimuth from E,N
Page 38