コード例 #1
0
ファイル: Bearing.java プロジェクト: qlongyinqw/thredds
  /**
   * Calculate the bearing between the 2 points. See calculateBearing below. Uses default Earth
   * object.
   *
   * @param pt1 Point 1
   * @param pt2 Point 2
   * @param result Object to use if non-null
   * @return The bearing
   */
  public static Bearing calculateBearing(LatLonPoint pt1, LatLonPoint pt2, Bearing result) {

    return calculateBearing(
        defaultEarth,
        pt1.getLatitude(),
        pt1.getLongitude(),
        pt2.getLatitude(),
        pt2.getLongitude(),
        result);
  }
コード例 #2
0
ファイル: Mercator.java プロジェクト: kenkyu/thredds
  /**
   * Convert a LatLonPoint to projection coordinates
   *
   * @param latLon convert from these lat, lon coordinates
   * @param result the object to write to
   * @return the given result
   */
  public ProjectionPoint latLonToProj(LatLonPoint latLon, ProjectionPointImpl result) {
    double toX, toY;
    double fromLat = latLon.getLatitude();
    double fromLon = latLon.getLongitude();
    double fromLat_r = Math.toRadians(fromLat);

    // infinite projection
    if ((Math.abs(90.0 - Math.abs(fromLat))) < TOLERANCE) {
      toX = Double.POSITIVE_INFINITY;
      toY = Double.POSITIVE_INFINITY;
    } else {
      toX = A * Math.toRadians(LatLonPointImpl.range180(fromLon - this.lon0));
      toY = A * SpecialMathFunction.atanh(Math.sin(fromLat_r)); // p 41 Snyder
    }

    result.setLocation(toX + falseEasting, toY + falseNorthing);
    return result;
  }
コード例 #3
0
  /**
   * set lat and lon by using UTM coordinates
   *
   * @param strZone UTM-zone, e.g. 32U
   * @param strNorthing Northing component
   * @param strEasting Easting component
   */
  public void set(String strZone, String strNorthing, String strEasting) {
    LatLonPoint ll = new LatLonPoint();

    utm.zone_letter = strZone.charAt(strZone.length() - 1);
    utm.zone_number = Convert.toInt(strZone.substring(0, strZone.length() - 1));
    utm.northing = (float) Common.parseDouble(strNorthing);
    utm.easting = (float) Common.parseDouble(strEasting);

    ll = utm.toLatLonPoint(); // returns null if unvalit UTM-coordinates
    if (ll != null) {
      this.utmValid = true;
      this.latDec = ll.getLatitude();
      this.lonDec = ll.getLongitude();
    } else {
      this.latDec = 91;
      this.lonDec = 361;
    }
  }
コード例 #4
0
ファイル: Orthographic.java プロジェクト: nbald/thredds
  /**
   * Convert a LatLonPoint to projection coordinates
   *
   * @param latLon convert from these lat, lon coordinates
   * @param result the object to write to
   * @return the given result
   */
  public ProjectionPoint latLonToProj(LatLonPoint latLon, ProjectionPointImpl result) {
    double toX, toY;
    double fromLat = latLon.getLatitude();
    double fromLon = latLon.getLongitude();

    fromLat = Math.toRadians(fromLat);
    double lonDiff = Math.toRadians(LatLonPointImpl.lonNormal(fromLon - lon0Degrees));
    double cosc = sinLat0 * Math.sin(fromLat) + cosLat0 * Math.cos(fromLat) * Math.cos(lonDiff);
    if (cosc >= 0) {
      toX = R * Math.cos(fromLat) * Math.sin(lonDiff);
      toY = R * (cosLat0 * Math.sin(fromLat) - sinLat0 * Math.cos(fromLat) * Math.cos(lonDiff));
    } else {
      toX = Double.POSITIVE_INFINITY;
      toY = Double.POSITIVE_INFINITY;
    }

    result.setLocation(toX, toY);
    return result;
  }
コード例 #5
0
  void timeProjection(ProjectionImpl proj) {
    java.util.Random r = new java.util.Random((long) this.hashCode());
    LatLonPointImpl startL = new LatLonPointImpl();

    double[][] from = new double[2][NPTS];
    for (int i = 0; i < NPTS; i++) {
      from[0][i] = (180.0 * (r.nextDouble() - .5)); // random latlon point
      from[1][i] = (360.0 * (r.nextDouble() - .5)); // random latlon point
    }

    int n = REPEAT * NPTS;

    // double[][] result = new double[2][NTRIALS];

    // normal
    long t1 = System.currentTimeMillis();
    for (int k = 0; k < REPEAT; k++) {
      for (int i = 0; i < NPTS; i++) {
        ProjectionPoint p = proj.latLonToProj(from[0][i], from[1][i]);
        LatLonPoint endL = proj.projToLatLon(p);

        if (checkit) {
          assert close(from[0][i], endL.getLatitude())
              : "lat: " + from[0][i] + "!=" + endL.getLatitude();
          assert close(from[1][i], endL.getLongitude())
              : "lon: " + from[1][i] + "!=" + endL.getLongitude();
        }
      }
    }
    long took = System.currentTimeMillis() - t1;
    sumNormal += took;
    System.out.println(
        n
            + " normal "
            + proj.getClassName()
            + " took "
            + took
            + " msecs "); // == "+ .001*took/NTRIALS+" secs/call ");

    // array
    long t2 = System.currentTimeMillis();
    for (int k = 0; k < REPEAT; k++) {
      double[][] to = proj.latLonToProj(from);
      double[][] result2 = proj.projToLatLon(to);

      if (checkit) {
        for (int i = 0; i < NPTS; i++) {
          assert close(from[0][i], result2[0][i]) : "lat: " + from[0][i] + "!=" + result2[0][i];
          assert close(from[1][i], result2[1][i]) : "lon: " + from[1][i] + "!=" + result2[1][i];
        }
      }
    }
    took = System.currentTimeMillis() - t2;
    sumArray += took;
    System.out.println(
        n
            + " array "
            + proj.getClassName()
            + " took "
            + took
            + " msecs "); // == "+ .001*took/NTRIALS+" secs/call ");
    //    System.out.println(NTRIALS +  " array "+ proj.getClassName()+" took "+took+ " msecs == "+
    // .001*took/NTRIALS/REPEAT+" secs/call ");

  }
コード例 #6
0
 /**
  * Create CWPoint by using a LatLonPoint
  *
  * @param CWPoint LatLonPoint
  */
 public CWPoint(LatLonPoint llPoint) {
   super(llPoint.getLatitude(), llPoint.getLongitude());
   this.utmValid = false;
 }
コード例 #7
0
 /**
  * Set CWPoint by using a LatLonPoint
  *
  * @param CWPoint LatLonPoint
  */
 public void set(LatLonPoint llPoint) {
   this.latDec = llPoint.getLatitude();
   this.lonDec = llPoint.getLongitude();
   this.utmValid = false;
 }
コード例 #8
0
ファイル: Bearing.java プロジェクト: qlongyinqw/thredds
 /**
  * Calculate a position given an azimuth and distance from another point. Uses default Earth.
  *
  * @param pt1 Point 1
  * @param az azimuth (degrees)
  * @param dist distance from the point (km)
  * @param result Object to use if non-null
  * @return The LatLonPoint
  * @see #findPoint(double,double,double,double,LatLonPointImpl)
  */
 public static LatLonPointImpl findPoint(
     LatLonPoint pt1, double az, double dist, LatLonPointImpl result) {
   return findPoint(defaultEarth, pt1.getLatitude(), pt1.getLongitude(), az, dist, result);
 }