北京54坐标转换为经纬度

发表于： 2020-10-26 2020-10-26
分类： SRE
标签： GIS

依据：

javascript

<!DOCTYPE html>
<html>
<head>
<meta charset="utf-8">
<title></title>
<script>
function bj54toGPS()
{

    x = 526511.014;     
    y = 3408291.638;
    pi = 3.1415926;
    x = x/1000000.0;
    y = y - 500000.0;

    bf = 9.04353692458*x-0.00001007623*Math.pow(x,2.0)-0.00074438304*Math.pow(x,3.0)-0.00000463064*Math.pow(x,4.0)+0.00000505846*Math.pow(x,5.0)-0.00000016754*Math.pow(x,6.0);
    hbf = bf * pi/ 180.0;

    sa = 6378245.0;
    sb = 6356863.019;
    se2 = 0.006693421623;
    sep2 = 0.006738525415;

    w1 = Math.sin(hbf);
    w2 = 1.0 - se2 * Math.pow(w1,2);
    w = Math.sqrt(w2);
    mf = sa*(1.0-se2)/Math.pow(w,3);
    w3 = Math.cos(hbf);


    w4 = Math.pow(sa,2)*Math.pow(w3,2) + Math.pow(sb,2)*Math.pow(w1,2);
    nf = Math.pow(sa,2) / Math.sqrt(w4);


    ynf = y/nf;
    vf = nf/mf;
    tf = Math.tan(hbf);

    yf2 = sep2 * Math.pow(w3, 2);

    gps_b = bf - 1.0/2.0 * vf * tf * (Math.pow(ynf,2)-1.0/12.0*(5.0+3.0*Math.pow(tf,2)+yf2-9.0*yf2*Math.pow(tf,2))*Math.pow(ynf,4))*180.0/pi;
    gps_l = 1.0/w3*ynf*(1.0-1.0/6.0*(1.0+2.0*Math.pow(tf,2)+yf2)*Math.pow(ynf,2)+1.0/120.0*(5.0+28.0*Math.pow(tf,2)+24.0*Math.pow(tf,2)+6.0*yf2+8.0*yf2*Math.pow(tf,2))*Math.pow(ynf,4))*180.0/pi;

    alert(gps_b);

}
</script>
</head>

<body>
<button onclick="bj54toGPS()">计算</button>
</body>
</html>

c++

#include <iostream>
#include <math.h>
using namespace std;


int main()
{
    double x = 3543664;
    double y = 21310994;
    double bf,vf,nf,ynf,tf,yf2,hbf;
    double sa,sb,se2,sep2,mf;
    double w1,w2,w,w3,w4;
    double pi = 3.1415926;
    double gps_l,gps_b;

    x = x/1000000.0;
    y = y - 500000.0;

    bf = 9.04353692458*x-0.00001007623*pow(x,2.0)-0.00074438304*pow(x,3.0)-0.00000463064*pow(x,4.0)+0.00000505846*pow(x,5.0)-0.00000016754*pow(x,6.0);
    hbf = bf * pi/ 180.0;
    sa = 6378245.0;
    sb = 6356863.019;
    se2 = 0.006693421623;
    sep2 = 0.006738525415;

     w1 = sin(hbf);
     w2 = 1.0 - se2 * pow(w1,(double)2);
     w = sqrt(w2);
     mf = sa*(1.0-se2)/pow(w,(double)3);
     w3 = cos(hbf);

    w4 = pow(sa,(double)2)*pow(w3,(double)2) + pow(sb,(double)2)*pow(w1,(double)2);
    nf = pow(sa,(double)2) / sqrt(w4);

    ynf = y/nf;
    vf = nf/mf;
    tf = tan(hbf);

    yf2 = sep2 * pow(w3, (double)2);

    gps_b = bf - 1.0/2.0 * vf * tf * (pow(ynf,(double)2)-1.0/12.0*(5.0+3.0*pow(tf,(double)2)+yf2-9.0*yf2*pow(tf,(double)2))*pow(ynf,(double)4))*180.0/pi;
    gps_l = 1.0/w3*ynf*(1.0-1.0/6.0*(1.0+2.0*pow(tf,(double)2)+yf2)*pow(ynf,(double)2)+1.0/120.0*(5.0+28.0*pow(tf,(double)2)+24.0*pow(tf,(double)2)+6.0*yf2+8.0*yf2*pow(tf,(double)2))*pow(ynf,(double)4))*180.0/pi;


   cout << gps_b;
   return 0;
}

sre

发表回复取消回复