double x_PI = 3.14159265358979323846 * 3000.0 / 180.0;
double PI = 3.1415926535897932384626;
double a = 6378245.0;
double ee = 0.00669342162296594323;
POSITION bd_pos;
POSITION gcj_pos;
POSITION wgs_pos;
bool outof_China(double lon, double lat)
{
return(lon<72.004 || lon>137.8374 || lat<0.8293 || lat >55.8271 || false);
}
POSITION bd09togcj02(double bd_lon, double bd_lat)
{
double x = bd_lon - 0.0065;
double y = bd_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);
gcj_pos.longitude = z*cos(theta);
gcj_pos.latitude = z*sin(theta);
return gcj_pos;
}
POSITION gcj02tobd09(double gcj_lon, double gcj_lat)
{
double z = sqrt(gcj_lon*gcj_lon + gcj_lat*gcj_lat) + 0.00002*sin(gcj_lat * x_PI);
double theta = atan2(gcj_lat,gcj_lon) + 0.000003 * cos