[GIS] Converting latitude, longitude (epsg:4326) into EPSG:3857

cconvertcoordinate systemmathematicspython

I'm looking to study the math involved to see how a given latitude and longitude in decimal degrees can be converted to EPSG:3857. Can anybody point to a good reference or a possible Ppython/C/C++/Java library source code for checking the math?

Best Answer

from pyproj import Proj, transform

P3857 = Proj(init='epsg:3857')
P4326 = Proj(init='epsg:4326')

x,y = transform(P4326, P3857, lon, lat)

Related Solutions

[GIS] Converting latitude longitude coordinate to x y coordinate

I assume your data are stored inside a PostGIS database. You can use PostGIS functions to do the transformations directly. No need for external libraries.

    select 
        st_x(st_transform(geom,<SRID>)) as X, 
        st_y(st_transform(geom,<srid>)) as Y 
    from 
        mTable 
    where 
        ...

[GIS] the formula for converting from WGS84 to EPSG:3071 (Wisconsin Transverse Mercator)

as whuber said, there are several pages of explanations to get the exact formulas. Though you can find some simplified series on wikipedia. For Wisconsin Transverse Mercator, please also note that implementing the projection is not enough, because you also need to code the datum transform between WGS 84 and NAD 83.

My suggestion is to take advantage of open source and pick the transverse Mercator projection in one of the open source libraries (e.g. proj.4).

Best Answer

Related Solutions

[GIS] Converting latitude longitude coordinate to x y coordinate

[GIS] the formula for converting from WGS84 to EPSG:3071 (Wisconsin Transverse Mercator)

Related Question