source: issm/trunk-jpl/src/m/coordsystems/ll2mercator.m

Last change on this file was 15144, checked in by Mathieu Morlighem, 12 years ago

CHG: added mercator2ll.m for reverse transformation
File size: 381 bytes
Line 
1function [mx my]=ll2mercator(lat, lon),
2%LL2MERCATOR - transform lat long to mercator projection
3%
4% Usage:
5% [mx my]=ll2mercator(lat, lon)
6
7EARTH_RADIUS = 6378137;
8EQUATOR_CIRCUMFERENCE = 2 * pi * EARTH_RADIUS;
9ORIGIN_SHIFT = EQUATOR_CIRCUMFERENCE / 2.0;
10
11mx = (lon * ORIGIN_SHIFT) / 180.0;
12my = log(tan((90 + lat) * pi/360.0))/(pi/180.0);
13my = (my * ORIGIN_SHIFT) /180.0;
Note: See TracBrowser for help on using the repository browser.