Changeset 8510 in josm for trunk/src/org/openstreetmap/josm/data/projection/Ellipsoid.java
- Timestamp:
- 2015-06-20T23:42:21+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/org/openstreetmap/josm/data/projection/Ellipsoid.java
r8451 r8510 29 29 */ 30 30 public static final Ellipsoid AustSA = Ellipsoid.create_a_rf(6378160.0, 298.25); 31 31 32 32 /** 33 33 * Bessel 1841 ellipsoid 34 34 */ 35 35 public static final Ellipsoid Bessel1841 = Ellipsoid.create_a_rf(6377397.155, 299.1528128); 36 36 37 37 /** 38 38 * Clarke 1866 ellipsoid … … 251 251 double v1 = 1-e*Math.sin(phi); 252 252 double v2 = 1+e*Math.sin(phi); 253 return Math.log(Math.tan(Math.PI/4+phi/2)*Math.pow(v1/v2, e/2));253 return Math.log(Math.tan(Math.PI/4+phi/2)*Math.pow(v1/v2, e/2)); 254 254 } 255 255 … … 262 262 double v1 = 1-e*Math.sin(phi); 263 263 double v2 = 1+e*Math.sin(phi); 264 return Math.log(Math.tan(Math.PI/4+phi/2)*Math.pow(v1/v2, e/2));264 return Math.log(Math.tan(Math.PI/4+phi/2)*Math.pow(v1/v2, e/2)); 265 265 } 266 266 … … 275 275 double lati = lat0; 276 276 double lati1 = 1.0; // random value to start the iterative processus 277 while (Math.abs(lati1-lati)>=epsilon) {277 while (Math.abs(lati1-lati) >= epsilon) { 278 278 lati = lati1; 279 279 double v1 = 1+e*Math.sin(lati); 280 280 double v2 = 1-e*Math.sin(lati); 281 lati1 = 2*Math.atan(Math.pow(v1/v2, e/2)*Math.exp(latIso))-Math.PI/2;281 lati1 = 2*Math.atan(Math.pow(v1/v2, e/2)*Math.exp(latIso))-Math.PI/2; 282 282 } 283 283 return lati1;
Note:
See TracChangeset
for help on using the changeset viewer.