Changeset 5073 in josm for trunk/src/org/openstreetmap/josm/data/projection/AbstractProjection.java

Timestamp:

2012-03-11T21:19:59+01:00 (13 years ago)

Author:

bastiK

Message:

move ntv2 to datum (see #7494)

File:

• trunk/src/org/openstreetmap/josm/data/projection/AbstractProjection.java (modified) (6 diffs)

trunk/src/org/openstreetmap/josm/data/projection/AbstractProjection.java

@@ -20 +20 @@
  * Subclasses of AbstractProjection must set ellps and proj to a non-null value.
  * In addition, either datum or nadgrid has to be initialized to some value.
- *
- * FIXME: nadgrids should probably be implemented as a Datum
  */
 abstract public class AbstractProjection implements Projection {
@@ -32 +30 @@
     protected double lon_0 = 0.0;   /* central meridian */
     protected double k_0 = 1.0;     /* general scale factor */
-    protected NTV2GridShiftFile nadgrids = null;
-    
+
     public final Ellipsoid getEllipsoid() {
         return ellps;
@@ -64 +61 @@
     @Override
     public EastNorth latlon2eastNorth(LatLon ll) {
-        if (nadgrids != null) {
-            NTV2GridShift gs = new NTV2GridShift(ll);
-            nadgrids.gridShiftReverse(gs);
-            ll = new LatLon(ll.lat()+gs.getLatShiftDegrees(), ll.lon()+gs.getLonShiftPositiveEastDegrees());
-        } else {
-            ll = datum.fromWGS84(ll);
-        }
+        ll = datum.fromWGS84(ll);
         double[] en = proj.project(Math.toRadians(ll.lat()), Math.toRadians(ll.lon() - lon_0));
         return new EastNorth(ellps.a * k_0 * en[0] + x_0, ellps.a * k_0 * en[1] + y_0);
@@ -79 +70 @@
         double[] latlon_rad = proj.invproject((en.east() - x_0) / ellps.a / k_0, (en.north() - y_0) / ellps.a / k_0);
         LatLon ll = new LatLon(Math.toDegrees(latlon_rad[0]), Math.toDegrees(latlon_rad[1]) + lon_0);
-        if (nadgrids != null) {
-            NTV2GridShift gs = new NTV2GridShift(ll);
-            nadgrids.gridShiftForward(gs);
-            ll = new LatLon(ll.lat()+gs.getLatShiftDegrees(), ll.lon()+gs.getLonShiftPositiveEastDegrees());
-        } else {
-            ll = datum.toWGS84(ll);
-        }
-        return ll;
+        return datum.toWGS84(ll);
     }
 
@@ -94 +78 @@
         return 10;
     }
-    
+
     /**
      * @return The EPSG Code of this CRS, null if it doesn't have one.
@@ -108 +92 @@
         return "EPSG:" + getEpsgCode();
     }
-    
+
     protected static final double convertMinuteSecond(double minute, double second) {
         return (minute/60.0) + (second/3600.0);
     }
-    
+
     protected static final double convertDegreeMinuteSecond(double degree, double minute, double second) {
         return degree + convertMinuteSecond(minute, second);