根据快速高斯投影方法，由经纬度计算西安80坐标系坐标

CSharpBoy

1. //根据快速高斯投影方法，由经纬度计算西安80坐标系坐标
   
2. // (double L,double B,double *x,double *y)
   
3. // 经度 100多,纬度 30多, x 七位,y 六位
   
4. ProjGAuss80(double L,double B,double *x,double *y)
   
5. {
   
6.   //得到与中央经线的经差
   
7.   double j1=L-(((int)L/6+1)*6-3);
   
8.   //化为弧度
   
9.   L=(j1/180.0)*PI;
   
10.   B=(B/180.0)*PI;
    
11.   //简化公式，中间参数无物理意义
    
12.   double cb=cos(B);
    
13.   
    
14.   double n;  
    
15.   n=(6399596.652-(21565.045-(108.996-0.603*cb*cb)*cb*cb)*cb*cb)*cb;
    
16.   
    
17.     double P0,P3,P4,P5,P6;
    
18.     P0=(32144.5189-(135.3646-0.7034*cb*cb)*cb*cb)*cb;
    
19.   P3=(0.3333333+0.0011233*cb*cb)*cb*cb-0.1666667;
    
20.     P4=(0.25+0.00253*cb*cb)*cb*cb-0.04167;
    
21.     P5=0.00878-(0.1702-0.20382*cb*cb)*cb*cb;
    
22.     P6=cb*cb*(0.167*cb*cb-0.083);
    
23.   
    
24.   *x=6367452.1328*B-(P0-(0.5+(P4+P6*L*L)*L*L)*L*L*n)*sin(B);
    
25.   *y=(1+(P3+P5*L*L)*L*L)*L*n+500000; //按照国家标准，将Y坐标加500km  
    
26. }

//根据快速高斯投影方法，由经纬度计算西安80坐标系坐标<br /> // (double L,double B,double *x,double *y)<br /> // 经度 100多,纬度 30多, x 七位,y 六位<br /> ProjGAuss80(double L,double B,double *x,double *y)<br /> {<br />   //得到与中央经线的经差<br />   double j1=L-(((int)L/6+1)*6-3);<br />   //化为弧度<br />   L=(j1/180.0)*PI;<br />   B=(B/180.0)*PI;<br />   //简化公式，中间参数无物理意义<br />   double cb=cos(B);<br />   <br />   double n;  <br />   n=(6399596.652-(21565.045-(108.996-0.603*cb*cb)*cb*cb)*cb*cb)*cb;<br />   <br />     double P0,P3,P4,P5,P6;<br />     P0=(32144.5189-(135.3646-0.7034*cb*cb)*cb*cb)*cb;<br />   P3=(0.3333333+0.0011233*cb*cb)*cb*cb-0.1666667;<br />     P4=(0.25+0.00253*cb*cb)*cb*cb-0.04167;<br />     P5=0.00878-(0.1702-0.20382*cb*cb)*cb*cb;<br />     P6=cb*cb*(0.167*cb*cb-0.083);<br />   <br />   *x=6367452.1328*B-(P0-(0.5+(P4+P6*L*L)*L*L)*L*L*n)*sin(B);<br />   *y=(1+(P3+P5*L*L)*L*L)*L*n+500000; //按照国家标准，将Y坐标加500km  <br /> }