| The constants in CassiniSoldner.java do not match the corresponding code in proj.4. The geotools class contains:
private static final double C1= 0.16666666666666666666,
C2= 0.08333333333333333333,
C3= 0.41666666666666666666,
C4= 0.33333333333333333333,
C5= 0.66666666666666666666;
compare this to
# define C1 .16666666666666666666
# define C2 .00833333333333333333
# define C3 .04166666666666666666
# define C4 .33333333333333333333
# define C5 .06666666666666666666
in PJ_cass.c. The constants C2, C3 and C5 are too large by a factor of 10. This seems to be a copy&paste error - as far as I can tell, the discrepancy is not compensated in the code below. The references (1),(2) confirm that C3 should be 1/24=0.041666. (1) Snyder: Map Projections - A Working Manual
(2) EPSG Guidence Note |
