preliminary version that uses Mercator projection and passes a modified version of the tests
benign test modifications
had to substantially increase epsilon due to additional error added by reprojection
use package rootSolve for much faster optimization
simplify sanity test: exact test not possible due to Mercator projection