public CassiniProjection() { projectionLatitude = Math.toRadians(0); projectionLongitude = Math.toRadians(0); minLongitude = Math.toRadians(-90); maxLongitude = Math.toRadians(90); initialize(); }
public void initialize() { super.initialize(); if (!spherical) { if ((en = ProjectionMath.enfn(es)) == null) throw new ProjectionException(); m0 = ProjectionMath.mlfn( projectionLatitude, Math.sin(projectionLatitude), Math.cos(projectionLatitude), en); } }
public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate xy) { if (spherical) { xy.x = Math.asin(Math.cos(lpphi) * Math.sin(lplam)); xy.y = Math.atan2(Math.tan(lpphi), Math.cos(lplam)) - projectionLatitude; } else { xy.y = ProjectionMath.mlfn(lpphi, n = Math.sin(lpphi), c = Math.cos(lpphi), en); n = 1. / Math.sqrt(1. - es * n * n); tn = Math.tan(lpphi); t = tn * tn; a1 = lplam * c; c *= es * c / (1 - es); a2 = a1 * a1; xy.x = n * a1 * (1. - a2 * t * (C1 - (8. - t + 8. * c) * a2 * C2)); xy.y -= m0 - n * tn * a2 * (.5 + (5. - t + 6. * c) * a2 * C3); } return xy; }
public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) { if (spherical) { out.y = Math.asin(Math.sin(dd = xyy + projectionLatitude) * Math.cos(xyx)); out.x = Math.atan2(Math.tan(xyx), Math.cos(dd)); } else { double ph1; ph1 = ProjectionMath.inv_mlfn(m0 + xyy, es, en); tn = Math.tan(ph1); t = tn * tn; n = Math.sin(ph1); r = 1. / (1. - es * n * n); n = Math.sqrt(r); r *= (1. - es) * n; dd = xyx / n; d2 = dd * dd; out.y = ph1 - (n * tn / r) * d2 * (.5 - (1. + 3. * t) * d2 * C3); out.x = dd * (1. + t * d2 * (-C4 + (1. + 3. * t) * d2 * C5)) / Math.cos(ph1); } return out; }
public void initialize() { double t; super.initialize(); if (Math.abs((t = Math.abs(projectionLatitude)) - ProjectionMath.HALFPI) < EPS10) mode = projectionLatitude < 0. ? SOUTH_POLE : NORTH_POLE; else mode = t > EPS10 ? OBLIQUE : EQUATOR; trueScaleLatitude = Math.abs(trueScaleLatitude); if (!spherical) { double X; switch (mode) { case NORTH_POLE: case SOUTH_POLE: if (Math.abs(trueScaleLatitude - ProjectionMath.HALFPI) < EPS10) akm1 = 2. * scaleFactor / Math.sqrt(Math.pow(1 + e, 1 + e) * Math.pow(1 - e, 1 - e)); else { akm1 = Math.cos(trueScaleLatitude) / ProjectionMath.tsfn(trueScaleLatitude, t = Math.sin(trueScaleLatitude), e); t *= e; akm1 /= Math.sqrt(1. - t * t); } break; case EQUATOR: akm1 = 2. * scaleFactor; break; case OBLIQUE: t = Math.sin(projectionLatitude); X = 2. * Math.atan(ssfn(projectionLatitude, t, e)) - ProjectionMath.HALFPI; t *= e; akm1 = 2. * scaleFactor * Math.cos(projectionLatitude) / Math.sqrt(1. - t * t); sinphi0 = Math.sin(X); cosphi0 = Math.cos(X); break; } } else { switch (mode) { case OBLIQUE: sinphi0 = Math.sin(projectionLatitude); cosphi0 = Math.cos(projectionLatitude); case EQUATOR: akm1 = 2. * scaleFactor; break; case SOUTH_POLE: case NORTH_POLE: akm1 = Math.abs(trueScaleLatitude - ProjectionMath.HALFPI) >= EPS10 ? Math.cos(trueScaleLatitude) / Math.tan(ProjectionMath.QUARTERPI - .5 * trueScaleLatitude) : 2. * scaleFactor; break; } } }
public StereographicAzimuthalProjection() { this(Math.toRadians(90.0), Math.toRadians(0.0)); }
private double ssfn(double phit, double sinphi, double eccen) { sinphi *= eccen; return Math.tan(.5 * (ProjectionMath.HALFPI + phit)) * Math.pow((1. - sinphi) / (1. + sinphi), .5 * eccen); }
public ProjCoordinate projectInverse(double x, double y, ProjCoordinate lp) { if (spherical) { double c, rh, sinc, cosc; sinc = Math.sin(c = 2. * Math.atan((rh = ProjectionMath.distance(x, y)) / akm1)); cosc = Math.cos(c); lp.x = 0.; switch (mode) { case EQUATOR: if (Math.abs(rh) <= EPS10) lp.y = 0.; else lp.y = Math.asin(y * sinc / rh); if (cosc != 0. || x != 0.) lp.x = Math.atan2(x * sinc, cosc * rh); break; case OBLIQUE: if (Math.abs(rh) <= EPS10) lp.y = projectionLatitude; else lp.y = Math.asin(cosc * sinphi0 + y * sinc * cosphi0 / rh); if ((c = cosc - sinphi0 * Math.sin(lp.y)) != 0. || x != 0.) lp.x = Math.atan2(x * sinc * cosphi0, c * rh); break; case NORTH_POLE: y = -y; case SOUTH_POLE: if (Math.abs(rh) <= EPS10) lp.y = projectionLatitude; else lp.y = Math.asin(mode == SOUTH_POLE ? -cosc : cosc); lp.x = (x == 0. && y == 0.) ? 0. : Math.atan2(x, y); break; } } else { double cosphi, sinphi, tp, phi_l, rho, halfe, halfpi; rho = ProjectionMath.distance(x, y); switch (mode) { case OBLIQUE: case EQUATOR: default: // To prevent the compiler complaining about uninitialized vars. cosphi = Math.cos(tp = 2. * Math.atan2(rho * cosphi0, akm1)); sinphi = Math.sin(tp); phi_l = Math.asin(cosphi * sinphi0 + (y * sinphi * cosphi0 / rho)); tp = Math.tan(.5 * (ProjectionMath.HALFPI + phi_l)); x *= sinphi; y = rho * cosphi0 * cosphi - y * sinphi0 * sinphi; halfpi = ProjectionMath.HALFPI; halfe = .5 * e; break; case NORTH_POLE: y = -y; case SOUTH_POLE: phi_l = ProjectionMath.HALFPI - 2. * Math.atan(tp = -rho / akm1); halfpi = -ProjectionMath.HALFPI; halfe = -.5 * e; break; } for (int i = 8; i-- != 0; phi_l = lp.y) { sinphi = e * Math.sin(phi_l); lp.y = 2. * Math.atan(tp * Math.pow((1. + sinphi) / (1. - sinphi), halfe)) - halfpi; if (Math.abs(phi_l - lp.y) < EPS10) { if (mode == SOUTH_POLE) lp.y = -lp.y; lp.x = (x == 0. && y == 0.) ? 0. : Math.atan2(x, y); return lp; } } throw new ConvergenceFailureException("Iteration didn't converge"); } return lp; }
public ProjCoordinate project(double lam, double phi, ProjCoordinate xy) { double coslam = Math.cos(lam); double sinlam = Math.sin(lam); double sinphi = Math.sin(phi); if (spherical) { double cosphi = Math.cos(phi); switch (mode) { case EQUATOR: xy.y = 1. + cosphi * coslam; if (xy.y <= EPS10) throw new ProjectionException(); xy.x = (xy.y = akm1 / xy.y) * cosphi * sinlam; xy.y *= sinphi; break; case OBLIQUE: xy.y = 1. + sinphi0 * sinphi + cosphi0 * cosphi * coslam; if (xy.y <= EPS10) throw new ProjectionException(); xy.x = (xy.y = akm1 / xy.y) * cosphi * sinlam; xy.y *= cosphi0 * sinphi - sinphi0 * cosphi * coslam; break; case NORTH_POLE: coslam = -coslam; phi = -phi; case SOUTH_POLE: if (Math.abs(phi - ProjectionMath.HALFPI) < TOL) throw new ProjectionException(); xy.x = sinlam * (xy.y = akm1 * Math.tan(ProjectionMath.QUARTERPI + .5 * phi)); xy.y *= coslam; break; } } else { double sinX = 0, cosX = 0, X, A; if (mode == OBLIQUE || mode == EQUATOR) { sinX = Math.sin(X = 2. * Math.atan(ssfn(phi, sinphi, e)) - ProjectionMath.HALFPI); cosX = Math.cos(X); } switch (mode) { case OBLIQUE: A = akm1 / (cosphi0 * (1. + sinphi0 * sinX + cosphi0 * cosX * coslam)); xy.y = A * (cosphi0 * sinX - sinphi0 * cosX * coslam); xy.x = A * cosX; break; case EQUATOR: A = 2. * akm1 / (1. + cosX * coslam); xy.y = A * sinX; xy.x = A * cosX; break; case SOUTH_POLE: phi = -phi; coslam = -coslam; sinphi = -sinphi; case NORTH_POLE: xy.x = akm1 * ProjectionMath.tsfn(phi, sinphi, e); xy.y = -xy.x * coslam; break; } xy.x = xy.x * sinlam; } return xy; }