Ignore:
Timestamp:
2018-04-15T18:22:42+02:00 (19 months ago)
Author:
Don-vip
Message:

see #16129 - Transverse Mercator: Minor coefficient error in tmerc inverse transform (see https://github.com/OSGeo/proj.4/issues/174)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/org/openstreetmap/josm/data/projection/proj/TransverseMercator.java

    r13167 r13634  
    8282
    8383    /**
    84      * Contants used for the forward and inverse transform for the eliptical
     84     * Constants used for the forward and inverse transform for the elliptical
    8585     * case of the Transverse Mercator.
    8686     */
     
    101101    /**
    102102     * A derived quantity of excentricity, computed by <code>e'² = (a²-b²)/b² = es/(1-es)</code>
    103      * where <var>a</var> is the semi-major axis length and <var>b</var> is the semi-minor axis
    104      * length.
     103     * where <var>a</var> is the semi-major axis length and <var>b</var> is the semi-minor axis length.
    105104     */
    106105    private double eb2;
     
    109108     * Latitude of origin in <u>radians</u>. Default value is 0, the equator.
    110109     * This is called '<var>phi0</var>' in Snyder.
    111      * <p>
    112      * <strong>Consider this field as final</strong>. It is not final only
    113      * because some classes need to modify it at construction time.
    114110     */
    115111    protected double latitudeOfOrigin;
     
    151147
    152148        if (params.gamma != null) {
    153                 rectifiedGridAngle = Utils.toRadians(params.gamma);
     149            rectifiedGridAngle = Utils.toRadians(params.gamma);
    154150        } else {
    155                 rectifiedGridAngle = 0.0;
     151            rectifiedGridAngle = 0.0;
    156152        }
    157153        sinrot = Math.sin(rectifiedGridAngle);
    158154        cosrot = Math.cos(rectifiedGridAngle);
    159 
    160155    }
    161156
     
    220215                FC4 * (5.0 + t*(3.0 - 9.0*n) + n*(1.0 - 4*n) - ds *
    221216                FC6 * (61.0 + t*(90.0 - 252.0*n + 45.0*t) + 46.0*n - ds *
    222                 FC8 * (1385.0 + t*(3633.0 + t*(4095.0 + 1574.0*t))))));
     217                FC8 * (1385.0 + t*(3633.0 + t*(4095.0 + 1575.0*t))))));
    223218
    224219            x = d*(FC1 - ds * FC3 * (1.0 + 2.0*t + n -
Note: See TracChangeset for help on using the changeset viewer.