public ProjectionPoint latLonToProj(LatLonPoint latLon, ProjectionPointImpl result) { double fromLat = Math.toRadians(latLon.getLatitude()); double theta = computeTheta(latLon.getLongitude()); double term = isSpherical ? n2 * Math.sin(fromLat) : n * MapMath.qsfn(Math.sin(fromLat), e, one_es); double rho = c - term; if (rho < 0.0) throw new RuntimeException("F"); rho = dd * Math.sqrt(rho); double toX = rho * Math.sin(theta); double toY = rho0 - rho * Math.cos(theta); result.setLocation(totalScale * toX + falseEasting, totalScale * toY + falseNorthing); return result; }
double ml1 = MapMath.qsfn(sinphi, e, one_es); if (secant) { /* secant cone */ sinphi = Math.sin(phi2); cosphi = Math.cos(phi2); double m2 = MapMath.msfn(sinphi, cosphi, es); double ml2 = MapMath.qsfn(sinphi, e, one_es); n = (m1 * m1 - m2 * m2) / (ml2 - ml1); c = m1 * m1 + n * ml1; dd = 1. / n; rho0 = dd * Math.sqrt(c - n * MapMath.qsfn(Math.sin(lat0rad), e, one_es));
public ProjectionPoint latLonToProj(LatLonPoint latLon, ProjectionPointImpl result) { double fromLat = Math.toRadians(latLon.getLatitude()); double theta = computeTheta(latLon.getLongitude()); double term = earth.isSpherical() ? n2 * Math.sin(fromLat) : n * MapMath.qsfn(Math.sin(fromLat), e, one_es); double rho = c - term; if (rho < 0.0) throw new RuntimeException("F"); rho = dd * Math.sqrt(rho); double toX = rho * Math.sin(theta); double toY = rho0 - rho * Math.cos(theta); result.setLocation(totalScale * toX + falseEasting, totalScale * toY + falseNorthing); return result; }
public ProjectionPoint latLonToProj(LatLonPoint latLon, ProjectionPointImpl result) { double fromLat = Math.toRadians(latLon.getLatitude()); double theta = computeTheta(latLon.getLongitude()); double term = earth.isSpherical() ? n2 * Math.sin(fromLat) : n * MapMath.qsfn(Math.sin(fromLat), e, one_es); double rho = c - term; if (rho < 0.0) throw new RuntimeException("F"); rho = dd * Math.sqrt(rho); double toX = rho * Math.sin(theta); double toY = rho0 - rho * Math.cos(theta); result.setLocation(totalScale * toX + falseEasting, totalScale * toY + falseNorthing); return result; }
double ml1 = MapMath.qsfn(sinphi, e, one_es); if (secant) { /* secant cone */ sinphi = Math.sin(phi2); cosphi = Math.cos(phi2); double m2 = MapMath.msfn(sinphi, cosphi, es); double ml2 = MapMath.qsfn(sinphi, e, one_es); n = (m1 * m1 - m2 * m2) / (ml2 - ml1); c = m1 * m1 + n * ml1; dd = 1. / n; rho0 = dd * Math.sqrt(c - n * MapMath.qsfn(Math.sin(lat0rad), e, one_es));
double ml1 = MapMath.qsfn(sinphi, e, one_es); if (secant) { /* secant cone */ sinphi = Math.sin(phi2); cosphi = Math.cos(phi2); double m2 = MapMath.msfn(sinphi, cosphi, es); double ml2 = MapMath.qsfn(sinphi, e, one_es); n = (m1 * m1 - m2 * m2) / (ml2 - ml1); c = m1 * m1 + n * ml1; dd = 1. / n; rho0 = dd * Math.sqrt(c - n * MapMath.qsfn(Math.sin(lat0rad), e, one_es));
public ProjectionPoint latLonToProj(LatLonPoint latLon, ProjectionPointImpl result) { double fromLat = Math.toRadians(latLon.getLatitude()); double theta = computeTheta(latLon.getLongitude()); double term = earth.isSpherical() ? n2 * Math.sin(fromLat) : n * MapMath.qsfn(Math.sin(fromLat), e, one_es); double rho = c - term; if (rho < 0.0) throw new RuntimeException("F"); rho = dd * Math.sqrt(rho); double toX = rho * Math.sin(theta); double toY = rho0 - rho * Math.cos(theta); result.setLocation(totalScale * toX + falseEasting, totalScale * toY + falseNorthing); return result; }
double ml1 = MapMath.qsfn(sinphi, e, one_es); if (secant) { /* secant cone */ sinphi = Math.sin(phi2); cosphi = Math.cos(phi2); double m2 = MapMath.msfn(sinphi, cosphi, es); double ml2 = MapMath.qsfn(sinphi, e, one_es); n = (m1 * m1 - m2 * m2) / (ml2 - ml1); c = m1 * m1 + n * ml1; dd = 1. / n; rho0 = dd * Math.sqrt(c - n * MapMath.qsfn(Math.sin(lat0rad), e, one_es));
public ProjectionPoint latLonToProj(LatLonPoint latlon, ProjectionPointImpl xy) { double lam = Math.toRadians( latlon.getLongitude() - lon0); double phi = Math.toRadians( latlon.getLatitude()); if (earth.isSpherical()) { xy.x = scaleFactor * lam; xy.y = Math.sin(phi) / scaleFactor; } else { xy.x = scaleFactor * lam; xy.y = .5 * MapMath.qsfn(Math.sin(phi), e, one_es) / scaleFactor; } xy.setLocation(totalScale * xy.x + falseEasting, totalScale * xy.y + falseNorthing); return xy; }
public ProjectionPoint latLonToProj(LatLonPoint latlon, ProjectionPointImpl xy) { double lam = Math.toRadians( latlon.getLongitude() - lon0); double phi = Math.toRadians( latlon.getLatitude()); if (earth.isSpherical()) { xy.setLocation(scaleFactor * lam, Math.sin(phi) / scaleFactor); } else { xy.setLocation(scaleFactor * lam, .5 * MapMath.qsfn(Math.sin(phi), e, one_es) / scaleFactor); } xy.setLocation(totalScale * xy.getX() + falseEasting, totalScale * xy.getY() + falseNorthing); return xy; }
public ProjectionPoint latLonToProj(LatLonPoint latlon, ProjectionPointImpl xy) { double lam = Math.toRadians( latlon.getLongitude() - lon0); double phi = Math.toRadians( latlon.getLatitude()); if (earth.isSpherical()) { xy.setLocation(scaleFactor * lam, Math.sin(phi) / scaleFactor); } else { xy.setLocation(scaleFactor * lam, .5 * MapMath.qsfn(Math.sin(phi), e, one_es) / scaleFactor); } xy.setLocation(totalScale * xy.getX() + falseEasting, totalScale * xy.getY() + falseNorthing); return xy; }
scaleFactor /= Math.sqrt(1. - es * t * t); apa = MapMath.authset(es); qp = MapMath.qsfn(1., e, one_es);
scaleFactor /= Math.sqrt(1. - es * t * t); apa = MapMath.authset(es); qp = MapMath.qsfn(1., e, one_es);
scaleFactor /= Math.sqrt(1. - es * t * t); apa = MapMath.authset(es); qp = MapMath.qsfn(1., e, one_es);