Changeset 8346 in josm for trunk/src/org/openstreetmap/josm/data/projection/proj/LambertConformalConic.java
- Timestamp:
- 2015-05-11T13:34:53+02:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/org/openstreetmap/josm/data/projection/proj/LambertConformalConic.java
r8345 r8346 15 15 import static org.openstreetmap.josm.tools.I18n.tr; 16 16 17 import org.openstreetmap.josm.data.projection.CustomProjection.Param; 17 18 import org.openstreetmap.josm.data.projection.Ellipsoid; 18 19 import org.openstreetmap.josm.data.projection.ProjectionConfigurationException; … … 60 61 * projection factor 61 62 */ 62 protected double F;63 protected double f; 63 64 /** 64 65 * radius of the parallel of latitude of the false origin (2SP) or at … … 76 77 ellps = params.ellps; 77 78 e = ellps.e; 78 if (params.lat _0 == null)79 throw new ProjectionConfigurationException(tr("Parameter ''{0}'' required.", "lat_0"));80 if (params.lat _1 != null && params.lat_2 != null) {81 initialize2SP(params.lat _0, params.lat_1, params.lat_2);79 if (params.lat0 == null) 80 throw new ProjectionConfigurationException(tr("Parameter ''{0}'' required.", Param.lat_0.key)); 81 if (params.lat1 != null && params.lat2 != null) { 82 initialize2SP(params.lat0, params.lat1, params.lat2); 82 83 } else { 83 initialize1SP(params.lat _0);84 initialize1SP(params.lat0); 84 85 } 85 86 } … … 103 104 104 105 n = (log(m1) - log(m2)) / (log(t1) - log(t2)); 105 F= m1 / (n * pow(t1, n));106 r0 = F* pow(tf, n);106 f = m1 / (n * pow(t1, n)); 107 r0 = f * pow(tf, n); 107 108 } 108 109 … … 120 121 121 122 n = sin(lat_0_rad); 122 F= m0 / (n * pow(t0, n));123 r0 = F* pow(t0, n);123 f = m0 / (n * pow(t0, n)); 124 r0 = f * pow(t0, n); 124 125 } 125 126 … … 153 154 double sinphi = sin(phi); 154 155 double L = (0.5*log((1+sinphi)/(1-sinphi))) - e/2*log((1+e*sinphi)/(1-e*sinphi)); 155 double r = F*exp(-n*L);156 double r = f*exp(-n*L); 156 157 double gamma = n*lambda; 157 158 double X = r*sin(gamma); … … 165 166 double gamma = atan(east / (r0-north)); 166 167 double lambda = gamma/n; 167 double latIso = (-1/n) * log(abs(r/ F));168 double latIso = (-1/n) * log(abs(r/f)); 168 169 double phi = ellps.latitude(latIso, e, epsilon); 169 170 return new double[] { phi, lambda };
Note:
See TracChangeset
for help on using the changeset viewer.
