原始到火星
void wgs2gcj(double &lat, double &lon)
{
double old_lon = lon - 105.0;
double old_lat = lat - 35.0;
double dLat = transformLat(old_lon, old_lat);
double dLon = transformLon(old_lon, old_lat);
double radLat = lat / 180.0 * pi;
double magic = sin(radLat);
magic = 1 - ee * magic * magic;
double sqrtMagic = sqrt(magic);
dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * pi);
dLon = (dLon * 180.0) / (a / sqrtMagic * cos(radLat) * pi);
lat = lat + dLat;
lon = lon + dLon;
}
火星到百度
void gcj2bd(double &lat, double &lon)
{
double x = lon, y = lat;
double z = sqrt(x * x + y * y) + 0.00002 * sin(y * x_pi);
double theta = atan2(y, x) + 0.000003 * cos(x * x_pi);
lon = z * cos(theta) + 0.0065;
lat = z * sin(theta) + 0.006;
}
原始到百度
void wgs2bd(double &lat, double &lon)
{
wgs2gcj(lat, lon);
gcj2bd (lat, lon);
}
百度到火星
void bd2gcj(double &lat, double &lon)
{
double x = lon - 0.0065, y = lat - 0.006;
double z = sqrt(x * x + y * y) - 0.00002 * sin(y * x_pi);
double theta = atan2(y, x) - 0.000003 * cos(x * x_pi);
lon = z * cos(theta);
lat = z * sin(theta);
}
额,没找到百度直接到原始坐标的,可以逆运算回去吧。。。