Estimación rápida de distancia / área desde coordenadas geodésicas

Dado un par de coordenadas (lon, lat) en grados, obtenemos las coordenadas en Mercator Esférico con: DEGREES_PER_RADIAN = 180.0 / Math.PI; MERCATOR_EARTH_RADIUS = 6378137.0; METERS_PER_EQUATOR_DEGREE = Math.PI * MERCATOR_EARTH_RADIUS / 180.0; rlat = lat / DEGREES_PER_RADIAN; y = 0.5 * Math.log((1 + Math.sin(rlat)) / (1 – Math.sin(rlat))); merx = METERS_PER_EQUATOR_DEGREE * lon; mery = METERS_PER_EQUATOR_DEGREE…