Changeset 5230 in josm for trunk/src

Timestamp:

2012-05-11T17:11:04+02:00 (13 years ago)

Author:

bastiK

Message:

use standard parameter set for lambert projection

Location:

trunk/src/org/openstreetmap/josm/data/projection

Files:

• Lambert.java (modified) (6 diffs)
• proj/LambertConformalConic.java (modified) (2 diffs)
• proj/ProjParameters.java (modified) (1 diff)

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

@@ -6 +6 @@
 import java.awt.GridBagLayout;
 import java.awt.event.ActionListener;
-import java.io.InputStream;
 import java.util.Collection;
 import java.util.Collections;
@@ -14 +13 @@
 import javax.swing.JPanel;
 
-import org.openstreetmap.josm.Main;
 import org.openstreetmap.josm.data.Bounds;
 import org.openstreetmap.josm.data.coor.LatLon;
 import org.openstreetmap.josm.data.projection.datum.NTV2Datum;
-import org.openstreetmap.josm.data.projection.datum.NTV2GridShiftFile;
 import org.openstreetmap.josm.data.projection.datum.NTV2GridShiftFileWrapper;
 import org.openstreetmap.josm.data.projection.proj.LambertConformalConic;
@@ -29 +26 @@
  * This newer version uses the grid translation NTF<->RGF93 provided by IGN for a submillimetric accuracy.
  * (RGF93 is the French geodetic system similar to WGS84 but not mathematically equal)
+ *
+ * Source: http://professionnels.ign.fr/DISPLAY/000/526/700/5267002/transformation.pdf
  * @author Pieren
  */
 public class Lambert extends AbstractProjection implements ProjectionSubPrefs {
-    /**
-     * Lambert I, II, III, and IV projection exponents
-     */
-    private static final double n[] = { 0.7604059656, 0.7289686274, 0.6959127966, 0.6712679322 };
-
-    /**
-     * Lambert I, II, III, and IV projection constants
-     */
-    private static final double c[] = { 11603796.98, 11745793.39, 11947992.52, 12136281.99 };
+
+    /**
+     * Lambert I, II, III, and IV latitude origin
+     */
+    private static final double lat_0s[] = { 49.5, 46.8, 44.1, 42.165 };
+    
+    /**
+     * Lambert I, II, III, and IV latitude of first standard parallel
+     */
+    private static final double lat_1s[] = { 
+            convertDegreeMinuteSecond(48, 35, 54.682), 
+            convertDegreeMinuteSecond(45, 53, 56.108),
+            convertDegreeMinuteSecond(43, 11, 57.449),
+            convertDegreeMinuteSecond(41, 33, 37.396)};
+    
+    /**
+     * Lambert I, II, III, and IV latitude of second standard parallel
+     */
+    private static final double lat_2s[] = {
+            convertDegreeMinuteSecond(50, 23, 45.282),
+            convertDegreeMinuteSecond(47, 41, 45.652),
+            convertDegreeMinuteSecond(44, 59, 45.938),
+            convertDegreeMinuteSecond(42, 46, 3.588)};
 
     /**
@@ -46 +59 @@
      */
     private static final double x_0s[] = { 600000.0, 600000.0, 600000.0, 234.358 };
-
+    
     /**
      * Lambert I, II, III, and IV false north
      */
-    private static final double y_fs[] = { 5657616.674, 6199695.768, 6791905.085, 7239161.542 };
-
+    private static final double y_0s[] = { 200000.0, 200000.0, 200000.0, 185861.369 };
+    
     /**
      * France is divided in 4 Lambert projection zones (1,2,3 + 4th for Corsica)
@@ -87 +100 @@
         datum = new NTV2Datum("ntf_rgf93Grid", null, ellps, NTV2GridShiftFileWrapper.ntf_rgf93);
         x_0 = x_0s[layoutZone];
+        y_0 = y_0s[layoutZone];
         lon_0 = 2.0 + 20.0 / 60 + 14.025 / 3600; // 0 grade Paris
         if (proj == null) {
@@ -95 +109 @@
             proj.initialize(new ProjParameters() {{
                 ellps = Lambert.this.ellps;
-                lcc_n = n[layoutZone];
-                lcc_F = c[layoutZone] / ellps.a;
-                lcc_r0 = y_fs[layoutZone] / ellps.a;
+                lat_0 = lat_0s[layoutZone];
+                lat_1 = lat_1s[layoutZone];
+                lat_2 = lat_2s[layoutZone];
             }});
         } catch (ProjectionConfigurationException e) {

trunk/src/org/openstreetmap/josm/data/projection/proj/LambertConformalConic.java

@@ -67 +67 @@
         ellps = params.ellps;
         e = ellps.e;
-        if (params.lcc_n != null && params.lcc_F != null && params.lcc_r0 != null) {
-            initializeDirect(params.lcc_n, params.lcc_F, params.lcc_r0);
+        if (params.lat_0 == null)
+            throw new ProjectionConfigurationException(tr("Parameter ''{0}'' required.", "lat_0"));
+        if (params.lat_1 != null && params.lat_2 != null) {
+            initialize2SP(params.lat_0, params.lat_1, params.lat_2);
         } else {
-            if (params.lat_0 == null)
-                throw new ProjectionConfigurationException(tr("Parameter ''{0}'' required.", "lat_0"));
-            if (params.lat_1 != null && params.lat_2 != null) {
-                initialize2SP(params.lat_0, params.lat_1, params.lat_2);
-            } else {
-                initialize1SP(params.lat_0);
-            }
+            initialize1SP(params.lat_0);
         }
     }
@@ -117 +113 @@
         F  = m0 / (n * pow(t0, n));
         r0 = F * pow(t0, n);
-    }
-
-    /**
-     * Initialize LCC by providing the projection parameters directly.
-     *
-     * @param n see field n
-     * @param F see field F
-     * @param r0 see field r0
-     */
-    private void initializeDirect(double n, double F, double r0) {
-        this.params = null;
-        this.n = n;
-        this.F = F;
-        this.r0 = r0;
     }
 

trunk/src/org/openstreetmap/josm/data/projection/proj/ProjParameters.java

@@ -15 +15 @@
     public Double lat_2;
 
-    /* for LambertConformalConic */
-    public Double lcc_n;
-    public Double lcc_F;
-    public Double lcc_r0;
-
 }