X1D1GDOPParam.cs 7.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206
  1. using DevExpress.XtraMap;
  2. using ExtensionsDev;
  3. using System;
  4. using System.Data;
  5. using System.Linq;
  6. using System.Collections.Generic;
  7. using XdCxRhDW.Repostory;
  8. using XdCxRhDW.Entity;
  9. using XdCxRhDW.Api;
  10. using System.Drawing;
  11. using DevExpress.XtraCharts.Design;
  12. using System.Runtime.ConstrainedExecution;
  13. namespace XdCxRhDW.App.UserControl
  14. {
  15. public partial class X1D1GDOPParam : DevExpress.XtraEditors.XtraUserControl
  16. {
  17. public MapControl mapControl1;
  18. public GDOP星地一星一地接口 Model => new GDOP星地一星一地接口()
  19. {
  20. EphMain = ucEphXYZMain.EphParam(),
  21. CapTime = txtCapTime.DateTime,
  22. CDBStation = txtStationLocation1.GetLonLat(),
  23. CXStation = txtCXStation.GetLonLat(),
  24. RefStation = txtRefLocation1.GetLonLat(),
  25. DtousErr = Convert.ToDouble(txtDtousErr1.Text),
  26. SatLocErr = Convert.ToDouble(txtSatLocErr1.Text),
  27. CXErr = Convert.ToDouble(txtCxErr.Text),
  28. };
  29. private EnumPosResType PosResType;
  30. public X1D1GDOPParam(PosRes item)
  31. {
  32. InitializeComponent();
  33. this.layoutControl1.UseDefault();
  34. txtCapTime.UseDefault();
  35. txtStationLocation1.UseDoubleClickToSelectAll();
  36. txtRefLocation1.UseDoubleClickToSelectAll();
  37. txtCXStation.UseDoubleClickToSelectAll();
  38. this.txtCapTime.DateTime = item.SigTime;
  39. this.txtDtousErr1.EditValue = 1;
  40. this.txtSatLocErr1.EditValue = 10000;
  41. this.txtCxErr.EditValue = 0.001;
  42. PosResType = item.PosResType;
  43. using (RHDWPartContext db = RHDWPartContext.GetContext(item.SigTime))
  44. {
  45. var cg = db?.CgRes.Where(m => m.ID == item.CgResID).FirstOrDefault();
  46. var station = db?.StationRes.Where(m => m.ID == item.StationResID).FirstOrDefault();
  47. if (station != null)
  48. {
  49. this.txtStationLocation1.Text = $"{station.CdbTxLon},{station.CdbTxLat}";
  50. this.txtRefLocation1.Text = $"{station.RefLon},{station.RefLat}";
  51. this.txtCXStation.Text = $"{station.CxLon},{station.CxLat}";
  52. }
  53. if (cg != null)
  54. {
  55. var mainEph = (cg.MainX, cg.MainY, cg.MainZ, cg.MainVx, cg.MainVy, cg.MainVz);
  56. ucEphXYZMain.SetParam($"主星", cg.MainCode, mainEph, Color.Black);
  57. }
  58. }
  59. }
  60. private bool ParamValidate()
  61. {
  62. dxErrorProvider.ClearErrors();
  63. if (!ucEphXYZMain.CheckEph(dxErrorProvider))
  64. {
  65. return false;
  66. }
  67. if (!txtStationLocation1.CheckLonLat(dxErrorProvider,"超短波"))
  68. {
  69. return false;
  70. }
  71. if (!txtCXStation.CheckLonLat(dxErrorProvider, "测向站"))
  72. {
  73. return false;
  74. }
  75. if (!txtRefLocation1.CheckLonLat(dxErrorProvider,"参考站"))
  76. {
  77. return false;
  78. }
  79. if (!txtDtousErr1.CheckDouble(dxErrorProvider, "时差误差"))
  80. {
  81. return false;
  82. }
  83. if (!txtSatLocErr1.CheckDouble(dxErrorProvider, "星历位置误差"))
  84. {
  85. return false;
  86. }
  87. if (!txtCxErr.CheckDouble(dxErrorProvider, "测向误差"))
  88. {
  89. return false;
  90. }
  91. return true;
  92. }
  93. private void btnOK_Click(object sender, EventArgs e)
  94. {
  95. if (!ParamValidate()) { return; }
  96. try
  97. {
  98. mapControl1.ClearMap();
  99. var data = GdopHelper.Gdop1Sat1DByXyz(Model.EphMain, Model.CDBStation, Model.CXStation,
  100. Model.DtousErr, Model.CXErr, Model.SatLocErr, Model.RefStation);
  101. if (data == null)
  102. {
  103. return;
  104. }
  105. foreach (var errLins in data)//画GDOP
  106. {
  107. var mapDots = errLins.MapDots.Select(p => p).Select(p => (p.Lon, p.Lat));
  108. mapControl1.DrawGdopLineTwo(errLins.ErrDistanceKm, mapDots, 1);
  109. }
  110. }
  111. catch (Exception ex)
  112. {
  113. DxHelper.MsgBoxHelper.ShowError($"绘制{PosResType.GetEnumDisplayName()}GDOP失败,{ex.Message}");
  114. }
  115. }
  116. private void btnClose_Click(object sender, EventArgs e)
  117. {
  118. DxHelper.PopupHelper.HidePopup(this);
  119. }
  120. private async void btnEphCalc_Click(object sender, EventArgs e)
  121. {
  122. dxErrorProvider.ClearErrors();
  123. if (txtCapTime.DateTime == DateTime.MinValue)
  124. {
  125. DxHelper.MsgBoxHelper.ShowWarning("信号时间不能为空!");
  126. return;
  127. }
  128. var sigTime = txtCapTime.DateTime;
  129. try
  130. {
  131. var mainCode = ucEphXYZMain.GetSatCode();
  132. var mainxlInfo = await XlRepository.GetLatestAsync(mainCode, sigTime);
  133. if (mainxlInfo == null)
  134. {
  135. DxHelper.MsgBoxHelper.ShowWarning($"卫星:{mainCode}未找到对应的星历信息,请导入星历");
  136. return;
  137. }
  138. var maineph = EphHelper.Calc(mainxlInfo.TwoLine, sigTime.ToUtc());
  139. ucEphXYZMain.SetParam("主星", mainCode, (maineph.X, maineph.Y, maineph.Z, maineph.VX, maineph.VY, maineph.VZ), Color.Red);
  140. }
  141. catch (Exception ex)
  142. {
  143. Serilog.Log.Error(ex, $"手动推算{PosResType.GetEnumDisplayName()}星历失败,SigTime={sigTime}");
  144. DxHelper.MsgBoxHelper.ShowError($"手动推算{PosResType.GetEnumDisplayName()}星历失败,{ex.Message}");
  145. }
  146. }
  147. }
  148. public class GDOP星地一星一地接口
  149. {
  150. /// <summary>
  151. /// 主星星历 x y z vx vy vz
  152. /// </summary>
  153. public double[] EphMain { get; set; }
  154. /// <summary>
  155. /// 采集时刻
  156. /// </summary>
  157. public DateTime CapTime { get; set; }
  158. /// <summary>
  159. /// 超短接收站-经度 纬度 高度
  160. /// </summary>
  161. public double[] CDBStation { get; set; }
  162. /// <summary>
  163. /// 参考站位置-经度 纬度 高度
  164. /// </summary>
  165. public double[] RefStation { get; set; }
  166. /// <summary>
  167. /// 测向站位置-经度 纬度 高度
  168. /// </summary>
  169. public double[] CXStation { get; set; }
  170. /// <summary>
  171. /// 时差误差(单位us)
  172. /// </summary>
  173. public double DtousErr { get; set; } = 1;
  174. /// <summary>
  175. /// 星历位置误差(单位米)
  176. /// </summary>
  177. public double SatLocErr { get; set; } = 10000;
  178. /// <summary>
  179. /// 测向误差(单位°)
  180. /// </summary>
  181. public double CXErr { get; set; } = 0.01;
  182. }
  183. }