- * This class has been derived from the implementation of the Geotools project; + * This class has been derived from the implementation of the Geotools project; * git 8cbf52d, org.geotools.referencing.operation.projection.AlbersEqualArea * at the time of migration. *
-* References: + * References: *
+ * + * This class has been derived from the implementation of the Geotools project; + * git 8cbf52d, org.geotools.referencing.operation.projection.CassiniSoldner + * at the time of migration. + */ +public class CassiniSoldner extends AbstractProj { + + /** + * Meridian distance at the {@code latitudeOfOrigin}. + * Used for calculations for the ellipsoid. + */ + private double ml0; + + /** + * Contants used for the forward and inverse transform for the eliptical + * case of the Cassini-Soldner. + */ + private static final double C1= 0.16666666666666666666, + C2= 0.008333333333333333333, + C3= 0.041666666666666666666, + C4= 0.33333333333333333333, + C5= 0.066666666666666666666; + + @Override + public String getName() { + return tr("Cassini-Soldner"); + } + + @Override + public String getProj4Id() { + return "cass"; + } + + @Override + public void initialize(ProjParameters params) throws ProjectionConfigurationException { + super.initialize(params); + if (params.lat0 == null) + throw new ProjectionConfigurationException(tr("Parameter ''{0}'' required.", "lat_0")); + double latitudeOfOrigin = Math.toRadians(params.lat0); + ml0 = mlfn(latitudeOfOrigin, Math.sin(latitudeOfOrigin), Math.cos(latitudeOfOrigin)); + } + + @Override + public double[] project(double phi, double lam) { + double sinphi = Math.sin(phi); + double cosphi = Math.cos(phi); + + double n = 1.0 / (Math.sqrt(1.0 - e2 * sinphi * sinphi)); + double tn = Math.tan(phi); + double t = tn * tn; + double a1 = lam * cosphi; + double c = cosphi * cosphi * e2 / (1 - e2); + double a2 = a1 * a1; + + double x = n * a1 * (1.0 - a2 * t * (C1 - (8.0 - t + 8.0 * c) * a2 * C2)); + double y = mlfn(phi, sinphi, cosphi) - ml0 + n * tn * a2 * (0.5 + (5.0 - t + 6.0 * c) * a2 * C3); + return new double[] {x, y}; + } + + @Override + public double[] invproject(double x, double y) { + double ph1 = inv_mlfn(ml0 + y); + double tn = Math.tan(ph1); + double t = tn * tn; + double n = Math.sin(ph1); + double r = 1.0 / (1.0 - e2 * n * n); + n = Math.sqrt(r); + r *= (1.0 - e2) * n; + double dd = x / n; + double d2 = dd * dd; + double phi = ph1 - (n * tn / r) * d2 * (0.5 - (1.0 + 3.0 * t) * d2 * C3); + double lam = dd * (1.0 + t * d2 * (-C4 + (1.0 + 3.0 * t) * d2 * C5)) / Math.cos(ph1); + return new double[] {phi, lam}; + } + + @Override + public Bounds getAlgorithmBounds() { + return new Bounds(-89, -1.0, 89, 1.0, false); + } +} Index: /trunk/test/unit/org/openstreetmap/josm/data/projection/ProjectionTest.java =================================================================== --- /trunk/test/unit/org/openstreetmap/josm/data/projection/ProjectionTest.java (revision 9548) +++ /trunk/test/unit/org/openstreetmap/josm/data/projection/ProjectionTest.java (revision 9549) @@ -123,4 +123,5 @@ testProj("stere", 1e-5, "+lat_0=90 +lat_ts=90"); testProj("omerc", EPS, "+lat_0=4 +lonc=115 +alpha=53 +no_uoff +gamma=53.130 +bounds=112,4,116,7"); + testProj("cass", 1e-3, "+lat_0=11 +bounds=-1.0,-89,1.0,89"); if (error2) {