wamblee.org Git - utils/blob - gps/src/org/wamblee/gpx/SphericalCoordinateSystem.java

   1 /*
   2  * Copyright 2006 the original author or authors.
   3  *
   4  * Licensed under the Apache License, Version 2.0 (the "License");
   5  * you may not use this file except in compliance with the License.
   6  * You may obtain a copy of the License at
   7  *
   8  *      http://www.apache.org/licenses/LICENSE-2.0
   9  *
  10  * Unless required by applicable law or agreed to in writing, software
  11  * distributed under the License is distributed on an "AS IS" BASIS,
  12  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
  13  * See the License for the specific language governing permissions and
  14  * limitations under the License.
  15  */
  16
  17 package org.wamblee.gpx;
  18
  19 import static java.lang.Math.PI;
  20 import static java.lang.Math.cos;
  21 import static java.lang.Math.sin;
  22
  23 /**
  24  * Represents the coordinate system for a GPS measurement identified by
  25  * <ul>
  26  *   <li> x1: latitude in degrees </li>
  27  *   <li> x2: longitude in degrees </li>
  28  *   <li> x3: elevation in meters </li>
  29  * </ul>
  30  * This coordinate system models the earth as a sphere of a specific radius.
  31  */
  32 public class SphericalCoordinateSystem implements CoordinateSystem {
  33     /**
  34      * Earth radius in meters.
  35      */
  36     private static final double EARTH_RADIUS = 6371000;
  37
  38
  39     /* (non-Javadoc)
  40      * @see org.wamblee.gpx.CoordinateSystem#toReferenceSystem(org.wamblee.gpx.Coordinates)
  41      */
  42     public Coordinates toReferenceSystem(Coordinates aCoordinates) {
  43         double latrad = radians(aCoordinates.getX1());
  44         double lonrad = radians(aCoordinates.getX2());
  45         double coslat = cos(latrad);
  46         double sinlat = sin(latrad);
  47         double coslon = cos(lonrad);
  48         double sinlon = sin(lonrad);
  49
  50         double trueElevation = EARTH_RADIUS + aCoordinates.getX3();
  51         return new Coordinates(trueElevation*coslat*coslon,
  52                                trueElevation*coslat*sinlon,
  53                                trueElevation*sinlat);
  54
  55     }
  56
  57     private double radians(double aDegrees) {
  58         return aDegrees/180.0*PI;
  59     }
  60 }