正算经度0.001m,反算经度0.0001“
正算公式:
N_metre = X + N * tB * Math.pow(cB * l, 2) / 2 + N * tB * (5 - Math.pow(tB, 2) + 9 * itaf + 4 * Math.pow(itaf, 2)) * Math.pow(cB * l, 4) / 24 + N * tB * (61 - 58 * Math.pow(tB, 2) + Math.pow(tB, 4)) * Math.pow(cB * l, 6) / 720;
E_metre = N * cB * l + N * (1 - Math.pow(tB, 2) + itaf) * Math.pow(cB * l, 3) / 6 + N * (5 - 18 * Math.pow(tB, 2) + Math.pow(tB, 4) + 14 * itaf - 58 * itaf * Math.pow(tB, 2)) * Math.pow(cB * l, 5) / 120 + 500000;
反算公式:
B_degree = Bf - tBf * Math.pow(y, 2) / (2 * Mf * Nf) + tBf * (5 + 3 * Math.pow(tBf, 2) + itaff - 9 * itaff * Math.pow(tBf, 2)) * Math.pow(y, 4) / (24 * Mf * Math.pow(Nf, 3)) - tBf * (61 + 90 * Math.pow(tBf, 2) + 45 * Math.pow(tBf, 4)) * Math.pow(y, 6) / (720 * Mf * Math.pow(Nf, 5));
L_degree = y / (Nf * cBf) - (1 + 2 * Math.pow(tBf, 2) + itaff) * Math.pow(y, 3) / (6 * Math.pow(Nf, 3) * cBf) + (5 + 28 * Math.pow(tBf, 2) + 24 * Math.pow(tBf, 4) + 6 * itaff + 8 * itaff * Math.pow(tBf, 2)) * Math.pow(y, 5) / (120 * Math.pow(Nf, 5) * cBf);
参考资料:孔祥元,郭际明,刘宗泉.大地测量学基础