子午线弧长计算公式X=α[B+βsin2B+γsin4B+δsin6B+εsin8B]a=椭球的长半径b=椭球的短半径n=(a-b)/(a+b)n2=n*nn3=n*n*nα=（a+b）(1+n2/4+n4/64)/2β=-3n/2+9nnn/16-3nnnnn/32γ=15nn/16-15nnnn/32δ=-35nnn/48+105nnnnn/256ε=315nnnn/512B=在G13、G14、G15(或H13、H14、H15;I13、I14、I15)中输入纬度的度、分、秒值B=纬度值（化为弧度值）X=弧长的值高斯投影正算e2=(a2-b2)/a2e'2(a2-b2)/b2N=卯酉圈曲率半径:a/√（1-e2sin2B）t=tgBt2=cos2Bη2=e'2cos2BX+t*N*cos2B*l2/2+t*N*cos4B(5-t2+9η2+4η4)*l4/24+t*N*cos6B(61-58t2+t4+270η2-330x=3111t2+543t4-t6)*l8/40320L=1232.146754980L=126213.13622.19976032l=;l2=;l3=;l4=0.053005340.0028095660.0001489227.89366E-06cosB;cos2B;co0.6205248450.3850510830.2389337640.148264337s3B;cos4Bx=X+t*N*cos2B*l2/2+t*N*cos4B(5-t2+9η2+4η4)*l4/24+t*N*cos6B(61-58t2+t4+270η2-33111t2+543t4-t6)*l8/40320t*N*cos2B*l2/2x的计算t*N*cos4B(5-t2+9η2+4η4)*l4/24t*N*cos6B(61-58t2+t4+270η2-330t2η2)l6/720t*Ncos8B(1385-3111t2+543t4-t6)*l8/40320x=y=N*cosB*l+N*cos3B(1-t2+η2)*l3/6+N*cos5B(5-18t2+t4+14η2-58t2η2)l5/120+Ncos7B(N*cosB*ly的计算N*cos3B(1-t2+η2)*l3/6N*cos5B(5-18t2+t4+14η2-58t2η2)l5/120Ncos7B(61-479t2+179t4-t6)*l7/5040y=WGS84椭球54椭球80椭球6378137637824563781406356752.3146356863.0196356755.2880.001679220.0016789790.0016792222.81978E-062.81897E-062.81979E-064.73503E-094.73299E-094.73504E-096367449.1466367558.4976367452.133-0.002518828-0.002518466-0.002518832.64354E-062.64278E-062.64355E-06-3.45262E-09-3.45114E-09-3.45263E-094.89183E-124.88902E-124.89185E-1251515138383843.902343.902343.90230.9013845150.9013845150.9013845155723903.635724004.1735723906.3030.7841867870.006694380.006693420.006694380.006739500.006738530.006739506391306.1186391412.456391309.1341.2637476051.2637476051.2637476051.5970580091.5970580091.5970580090.3850510830.3850510830.3850510830.0025950510.0025946770.002595052*cos6B(61-58t2+t4+270η2-330t2η2)l6/720+t*Ncos8B(1385-*N*cos6B(61-58t2+t4+270η2-330t2η2)l6/720+t*Ncos8B(1385-4368.9549744369.027664368.9570361.3495295541.3495506771.349530198-0.000422489-0.000422495-0.000422489-6.04345E-07-6