POSC Specifications
Version 2.2
Epicentre Usage Guide
Projections and Projected Coordinate Systems

Previous Section
Projection Formula
Table of Contents
Next Section

3.4.4.3 Cassini-Soldner

The Cassini-Soldner projection is the ellipsoidal version of the Cassini projection for the sphere. It is not conformal but as it is relatively simple to construct it was extensively used in the last century and is still useful for mapping areas with limited longitudinal extent. It has now largely been replaced by the conformal Transverse Mercator which it resembles. Like this, it has a straight central meridian along which the scale is true, all other meridians and parallels are curved, and the scale distortion increases rapidly with increasing distance from the central meridian.

The formulas to derive projected Easting and Northing coordinates are:

Easting:

Northing:

where

and M, the distance along the meridian from equator to latitude j, is given by:

with j in radians.

M0 is the value of M calculated for the latitude of the chosen origin. This may not necessarily be chosen as the equator.

To compute latitude and longitude from Easting and Northing the reverse formulas are:

where r1 is r calculated at j = j1, and j1 is the latitude of the point on the central meridian which has the same Northing as the point whose coordinates are sought, and is found from:

where:

Example:

For Projected Coordinate System Trinidad 1903 / Trinidad Grid

Parameters:

Ellipsoid

Clarke 1858

a = 20926348 ft = 31706587.879 lkCla

 

 

 

b = 20855233 ft

then 1/f = 294.26068

 

 

e = 0.082371998

 

 

e2 = 0.006785146

Latitude of Natural Origin

j0

10°26' 30"N =

0.182241463 rad

Longitude of Natural Origin

l0

61°20' 00"W =

= -1.070468608 rad1

False Eastings

FE

430000.00 lkCla

 

False Northings

FN

325000.00 lkCla

 

Note that the semimajor axis is given in links (Clarke's ratio). This will also be true of the False Easting and False Northing below. The final Easting and Northing values are also in this unit.

Forward calculation for:

Latitude

j

10°00' 00.00" N =

0.174532925rad

Longitude

l

62°00' 00.00" W =

= -1.082104136 rad

 

A =

-0.01145876

C =

0.00662550

T =

0.03109120

M =

5496860.24

n =

31709831.92

M0 =

5739691.12

gives the following easting and northing:

Easting

E = 66644.94 lkCla

Northing

N = 82536.22 lkCla

Reverse the calculations. Use the same parameters and use the E and N values to calculate the latitude and longitude:

D =

-0.01145875

e1 =

0.00170207

T1 =

0.03109544

M1 =

5497227.34

n1 =

31709832.34

m1 =

0.17367306

r1 =

31501122.40

j1 =

0.17454458

Then

Latitude

j =

10°00' 00.00" N

Longitude

l =

62°00' 00.00" W


Previous Section
Projection Formula
Table of Contents
Next Section

1Value corrected from 1.0704468608 to 1.070468608 in 2006-01


Last modified: 12 July 2000
© Copyright 1997-2000 POSC. All rights reserved.