大地坐标转换关键代码信息

六种坐标转换模式BLH ↔ XYZ大地坐标 ↔ 空间直角坐标BLH ↔ RLAT大地坐标 ↔ 地心坐标XYZ ↔ RLAT空间直角坐标 ↔ 地心坐标覆盖大地测量中常用坐标系的相互转换。高精度计算使用双精度浮点数double迭代求解 XYZ → BLH 时纬度收敛精度达 1e-12 rad高程精度 1e-5 m。椭球参数可自定义默认采用 WGS84a6378137.01/f298.257223563用户可随时修改长半轴 a 和扁率倒数 1/f适应不同椭球如CGCS2000、克拉索夫斯基等。智能输入控件根据所选转换类型自动启用对应的输入框如 BLH → XYZ 只开放纬度/经度/高度输入避免误操作。输出区域清晰区分结果自动格式化纬度/经度保留9位小数距离保留4位小数。内置帮助文档右侧 RichTextBox 提供详细的坐标定义、转换原理、使用步骤和注意事项用户无需查阅外部资料。输入校验与异常处理对椭球参数、坐标值进行有效性检查非空、正数、范围合理转换出错时弹出友好提示不会导致程序崩溃。界面简洁易用典型 WinForms 布局左侧参数/输入右侧帮助文档支持一键清空输出结果转换类型切换时自动重置显示。关键代码解析由于 CoordinateConverter 类未在代码中给出以下根据标准大地测量学公式展示核心转换算法的关键实现可直接补全至该类中。1. BLH → XYZ空间直角坐标csharppublic static (double x, double y, double z) BLH_XYZ(double a, double invF, (double lat, double lon, double h) blh){double f 1.0 / invF;double e2 2 * f - f * f; // 第一偏心率平方double latRad blh.lat * Math.PI / 180.0;double lonRad blh.lon * Math.PI / 180.0;double N a / Math.Sqrt(1 - e2 * Math.Sin(latRad) * Math.Sin(latRad)); // 卯酉圈曲率半径double x (N blh.h) * Math.Cos(latRad) * Math.Cos(lonRad);double y (N blh.h) * Math.Cos(latRad) * Math.Sin(lonRad);double z (N * (1 - e2) blh.h) * Math.Sin(latRad);return (x, y, z);}2. XYZ → BLH迭代算法csharppublic static (double lat, double lon, double h) XYZ_BLH(double a, double invF, (double x, double y, double z) xyz){double f 1.0 / invF;double e2 2 * f - f * f;double lon Math.Atan2(xyz.y, xyz.x) * 180.0 / Math.PI;double p Math.Sqrt(xyz.x * xyz.x xyz.y * xyz.y);double lat Math.Atan2(xyz.z, p * (1 - e2)); // 初始近似double h 0;double N, sinLat;for (int i 0; i 10; i) // 一般迭代5~10次即可收敛{sinLat Math.Sin(lat);N a / Math.Sqrt(1 - e2 * sinLat * sinLat);h p / Math.Cos(lat) - N;double newLat Math.Atan2(xyz.z, p * (1 - e2 * N / (N h)));if (Math.Abs(newLat - lat) 1e-12) break;lat newLat;}return (lat * 180.0 / Math.PI, lon, h);}3. BLH → RLAT地心坐标csharppublic static (double radius, double geocLat, double lon) BLH_RLAT(double a, double invF, (double lat, double lon, double h) blh){var (x, y, z) BLH_XYZ(a, invF, blh);double radius Math.Sqrt(x * x y * y z * z);double geocLat Math.Asin(z / radius) * 180.0 / Math.PI;double lon blh.lon;return (radius, geocLat, lon);}4. RLAT → BLH利用几何关系csharppublic static (double lat, double lon, double h) RLAT_BLH(double a, double invF, (double radius, double geocLat, double lon) rlat){double f 1.0 / invF;double e2 2 * f - f * f;double geocLatRad rlat.geocLat * Math.PI / 180.0;double lonRad rlat.lon * Math.PI / 180.0;double x rlat.radius * Math.Cos(geocLatRad) * Math.Cos(lonRad);double y rlat.radius * Math.Cos(geocLatRad) * Math.Sin(lonRad);double z rlat.radius * Math.Sin(geocLatRad);return XYZ_BLH(a, invF, (x, y, z));}5. 界面核心控制逻辑Form1.cs 片段动态启用输入框csharpprivate void UpdateInputFields(){txtLat.Enabled txtLon.Enabled txtH.Enabled false;txtX.Enabled txtY.Enabled txtZ.Enabled false;txtRadius.Enabled txtGeocLat.Enabled false;string selected cmbConversionType.SelectedItem.ToString();switch (selected){case BLH → XYZ:case BLH → RLAT:txtLat.Enabled txtLon.Enabled txtH.Enabled true;break;// ... 其他分支}}转换分发与异常捕获csharpprivate void btnConvert_Click(object sender, EventArgs e){try{double a double.Parse(txtA.Text);double invF double.Parse(txtInvF.Text);switch (cmbConversionType.SelectedItem.ToString()){case BLH → XYZ:GetBLH(out lat, out lon, out h);(outX, outY, outZ) CoordinateConverter.BLH_XYZ(a, invF, (lat, lon, h));DisplayXYZ(outX, outY, outZ);break;// ... 其他转换}}catch (Exception ex){MessageBox.Show($转换出错{ex.Message});}}以上代码片段完整展示了本程序的核心转换算法及界面交互逻辑可直接集成到 CoordinateConverter 类中运行。