[Proj] Cassini coordinate reprojection problem

[Gerald Evenden]

Just for informational purposes I did another script for Cassini based 
upon
Snyder's Working Manual:

echo test from p. 289-290 pp 1395
lproj +proj=cass +a=6378206.4 +es=0.00676866 +lat_0=40 +lon_0=-75 <<EOF
-73 43
EOF
echo results from Snyder p.289: 163,071.1  335,127.6
echo turn-around test
lproj -o +proj=cass +a=6378206.4 +es=0.00676866 +lat_0=40 +lon_0=-75 
<<EOF | \
lproj -i -I +proj=cass +a=6378206.4 +es=0.00676866 +lat_0=40 +lon_0=-75
-73 43
EOF

The results of the script are:

test from p. 289-290 pp 1395
163071.14       335127.59
results from Snyder p.289: 163,071.1 335,127.6
turn-around test
73dW    43dN

You might put this script through your local copy of proj for test 
purposes.

I did not bother with tmerc as I feel that that its high usage is 
testament to its
proper functioning.

_____________________________________
Jerry and the low riders: Daisy Mae and Joshua.
"The whole religious complexion of the modern world is due to the
absence from Jerusalem of a lunatic asylum." Havelock Ellis, 1914