Context Navigation

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

Rev	Line
[15139]	1	function [mx my]=ll2mercator(lat, lon),
	2	%LL2MERCATOR - transform lat long to mercator projection
	3	%
	4	% Usage:
	5	% [mx my]=ll2mercator(lat, lon)
	6
	7	EARTH_RADIUS = 6378137;
	8	EQUATOR_CIRCUMFERENCE = 2 * pi * EARTH_RADIUS;
	9	ORIGIN_SHIFT = EQUATOR_CIRCUMFERENCE / 2.0;
	10
	11	mx = (lon * ORIGIN_SHIFT) / 180.0;
	12	my = log(tan((90 + lat) * pi/360.0))/(pi/180.0);
	13	my = (my * ORIGIN_SHIFT) /180.0;

Note: See TracBrowser for help on using the repository browser.