osgearth经纬度转墨卡托坐标

一、代码

#include <iostream>
#include <Windows.h>
#include <osgEarth/MapNode>
#include <osgEarthUtil/EarthManipulator>
#include <osgEarthUtil/ExampleResources>
#include <osgViewer/Viewer>
#include <osgDB/WriteFile>
#include <osg/LineWidth>
using namespace std;

const double epsilon = 0.000000000000001;
const double pi = 3.14159265358979323846;
const double d2r = pi / 180;
const double r2d = 180 / pi;

const double a = 6378137.0;		//椭球长半轴
const double f_inverse = 298.257223563;			//扁率倒数
const double b = a - a / f_inverse;
//const double b = 6356752.314245;			//椭球短半轴
const double e = sqrt(a * a - b * b) / a;

osg::Node* aix = osgDB::readNodeFile("axes.osgt");