~flavour/sahana-eden/trunk

/* Copyright (c) 2006-2012 by OpenLayers Contributors (see authors.txt for 

 * Namespace: OpenLayers.Projection

 * Methods for coordinate transforms between coordinate systems.  By default,

 *     OpenLayers ships with the ability to transform coordinates between

 *     geographic (EPSG:4326) and web or spherical mercator (EPSG:900913 et al.)

 *     coordinate reference systems.  See the <transform> method for details

 *     on usage.

 *

 * Additional transforms may be added by using the <proj4js at http://proj4js.org/>

 *     library.  If the proj4js library is included, the <transform> method 

 *     will work between any two coordinate reference systems with proj4js 

 *     definitions.

 *

 * If the proj4js library is not included, or if you wish to allow transforms

 *     between arbitrary coordinate reference systems, use the <addTransform>

 *     method to register a custom transform method.

        } else {

            var sourceCode = source.getCode();

            var destCode = dest.getCode();

            var transforms = OpenLayers.Projection.transforms;

            if (transforms[sourceCode] && transforms[sourceCode][destCode]) {

                transforms[sourceCode][destCode](point);

            }

        }

 

/**

 * Note: Transforms for web mercator <-> EPSG:4326

 * OpenLayers recognizes EPSG:3857, EPSG:900913, EPSG:102113 and EPSG:102100.

 * OpenLayers originally started referring to EPSG:900913 as web mercator.

 * The EPSG has declared EPSG:3857 to be web mercator.

 * ArcGIS 10 recognizes the EPSG:3857, EPSG:102113, and EPSG:102100 as

 * equivalent.  See http://blogs.esri.com/Dev/blogs/arcgisserver/archive/2009/11/20/ArcGIS-Online-moving-to-Google-_2F00_-Bing-tiling-scheme_3A00_-What-does-this-mean-for-you_3F00_.aspx#12084

 */

(function() {

 

    var pole = 20037508.34;

 

    function inverseMercator(xy) {

        xy.x = 180 * xy.x / pole;

        xy.y = 180 / Math.PI * (2 * Math.atan(Math.exp((xy.y / pole) * Math.PI)) - Math.PI / 2);

        return xy;

    }

 

    function forwardMercator(xy) {

        xy.x = xy.x * pole / 180;

        xy.y = Math.log(Math.tan((90 + xy.y) * Math.PI / 360)) / Math.PI * pole;

        return xy;

    }

 

    // list of equivalent codes for web mercator

    var codes = ["EPSG:900913", "EPSG:3857", "EPSG:102113", "EPSG:102100"];

 

    var add = OpenLayers.Projection.addTransform;

    var same = OpenLayers.Projection.nullTransform;

 

    var i, len, code, other, j;

    for (i=0, len=codes.length; i<len; ++i) {

        code = codes[i];

        add("EPSG:4326", code, forwardMercator);

        add(code, "EPSG:4326", inverseMercator);

        for (j=i+1; j<len; ++j) {

            other = codes[j];

            add(code, other, same);

            add(other, code, same);

        }

    }

 

})();