X3GDOPParam.cs 8.6 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213
  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 DxHelper;
  8. using XdCxRhDW.Repostory;
  9. using XdCxRhDW.Entity;
  10. using XdCxRhDW.Api;
  11. using System.Drawing;
  12. using DevExpress.XtraEditors.DXErrorProvider;
  13. using XdCxRhDW.Dto;
  14. namespace XdCxRhDW.App.UserControl
  15. {
  16. public partial class X3GDOPParam : DevExpress.XtraEditors.XtraUserControl
  17. {
  18. public MapControl mapControl1;
  19. public X3GDOP星地接口 Model => new X3GDOP星地接口()
  20. {
  21. EphMain = ucEphXYZMain.EphParam(),
  22. EphAdja1 = ucEphXYZAdja1.EphParam(),
  23. EphAdja2 = ucEphXYZAdja2.EphParam(),
  24. CapTime = txtCapTime.DateTime,
  25. RefStation = PosResType == EnumPosResType.X3NoRef ? new double[3] { 0, 0, 0 } : txtRefLocation1.GetLonLat(),
  26. DtousErr = Convert.ToDouble(txtDtousErr1.Text),
  27. SatLocErr = Convert.ToDouble(txtSatLocErr1.Text),
  28. };
  29. private EnumPosResType PosResType;
  30. public X3GDOPParam(PosRes item)
  31. {
  32. InitializeComponent();
  33. this.layoutControl1.UseDefault();
  34. txtCapTime.UseDefault();
  35. txtRefLocation1.UseDoubleClickToSelectAll();
  36. this.txtCapTime.DateTime = item.SigTime;
  37. this.txtDtousErr1.EditValue = 1;
  38. this.txtSatLocErr1.EditValue = 10000;
  39. PosResType = item.PosResType;
  40. if (PosResType == EnumPosResType.X3NoRef)
  41. {
  42. layoutControlItem14.Visibility = DevExpress.XtraLayout.Utils.LayoutVisibility.Never;
  43. }
  44. using (RHDWPartContext db = RHDWPartContext.GetContext(item.SigTime))
  45. {
  46. var cg = db?.CgRes.Where(m => m.ID == item.CgResID).FirstOrDefault();
  47. var station = db?.StationRes.Where(m => m.ID == item.StationResID).FirstOrDefault();
  48. if (station != null)
  49. {
  50. this.txtRefLocation1.Text = $"{station.RefLon},{station.RefLat}";
  51. }
  52. if (cg != null)
  53. {
  54. var mainEph = (cg.MainX, cg.MainY, cg.MainZ, cg.MainVx, cg.MainVy, cg.MainVz);
  55. ucEphXYZMain.SetParam($"主星", cg.MainCode, mainEph, Color.Black);
  56. var adja1Eph = (cg.Adja1X, cg.Adja1Y, cg.Adja1Z, cg.Adja1Vx, cg.Adja1Vy, cg.Adja1Vz);
  57. ucEphXYZAdja1.SetParam($"邻星1", cg.Adja1Code, adja1Eph, Color.Black);
  58. var adja2Eph = (cg.Adja2X, cg.Adja2Y, cg.Adja2Z, cg.Adja2Vx, cg.Adja2Vy, cg.Adja2Vz);
  59. ucEphXYZAdja2.SetParam($"邻星2", cg.Adja2Code, adja2Eph, Color.Black);
  60. }
  61. }
  62. }
  63. private void btnClose_Click(object sender, EventArgs e)
  64. {
  65. DxHelper.PopupHelper.HidePopup(this);
  66. }
  67. private bool ParamValidate()
  68. {
  69. if (!ucEphXYZMain.CheckEph(dxErrorProvider))
  70. {
  71. return false;
  72. }
  73. if (!ucEphXYZAdja1.CheckEph(dxErrorProvider))
  74. {
  75. return false;
  76. }
  77. if (!ucEphXYZAdja2.CheckEph(dxErrorProvider))
  78. {
  79. return false;
  80. }
  81. if (PosResType == EnumPosResType.X3 && !txtRefLocation1.CheckLonLat(dxErrorProvider, "参考站"))
  82. {
  83. return false;
  84. }
  85. if (!txtDtousErr1.CheckDouble(dxErrorProvider, "时差误差"))
  86. {
  87. return false;
  88. }
  89. if (!txtSatLocErr1.CheckDouble(dxErrorProvider, "星历位置误差"))
  90. {
  91. return false;
  92. }
  93. return true;
  94. }
  95. private void btnX1D2_Click(object sender, EventArgs e)
  96. {
  97. if (!ParamValidate()) { return; }
  98. try
  99. {
  100. mapControl1.ClearMap();
  101. var data = GdopHelper.Gdop3SatByXyz(Model.EphMain, Model.EphAdja1, Model.EphAdja2
  102. , Model.DtousErr, Model.SatLocErr, PosResType == EnumPosResType.X3NoRef ? null : Model.RefStation);
  103. if (data == null)
  104. {
  105. return;
  106. }
  107. foreach (var errLins in data)//画GDOP
  108. {
  109. var mapDots = errLins.MapDots.Select(p => p).Select(p => (p.Lon, p.Lat));
  110. mapControl1.DrawGdopLineTwo(errLins.ErrDistanceKm, mapDots, 1);
  111. }
  112. }
  113. catch (Exception ex)
  114. {
  115. DxHelper.MsgBoxHelper.ShowError($"绘制{PosResType.GetEnumDisplayName()}GDOP失败,{ex.Message}");
  116. }
  117. }
  118. private async void btnEphCalc_Click(object sender, EventArgs e)
  119. {
  120. dxErrorProvider.ClearErrors();
  121. if (txtCapTime.DateTime == DateTime.MinValue)
  122. {
  123. dxErrorProvider.SetError(txtCapTime, "信号时间不能为空!");
  124. return;
  125. }
  126. var sigTime = txtCapTime.DateTime;
  127. try
  128. {
  129. var mainCode = ucEphXYZMain.GetSatCode();
  130. var adja1Code = ucEphXYZAdja1.GetSatCode();
  131. var adja2Code = ucEphXYZAdja2.GetSatCode();
  132. var mainxlInfo = await XlRepository.GetLatestAsync(mainCode, sigTime);
  133. if (mainxlInfo == null)
  134. {
  135. DxHelper.MsgBoxHelper.ShowWarning($"卫星:{mainCode}未找到对应的星历信息,请导入星历");
  136. return;
  137. }
  138. var XlCalcDto = new XlCalcDto() { tleStr = mainxlInfo.TwoLine, SigTime = sigTime };
  139. var maineph = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
  140. ucEphXYZMain.SetParam("主星", mainCode, (maineph.data.X, maineph.data.Y, maineph.data.Z, maineph.data.VX, maineph.data.VY, maineph.data.VZ), Color.Red);
  141. var adja1xlInfo = await XlRepository.GetLatestAsync(adja1Code, sigTime);
  142. if (adja1xlInfo == null)
  143. {
  144. DxHelper.MsgBoxHelper.ShowWarning($"卫星:{adja1Code}未找到对应的星历信息,请导入星历");
  145. return;
  146. }
  147. XlCalcDto = new XlCalcDto() { tleStr = adja1xlInfo.TwoLine, SigTime = sigTime };
  148. var adja1eph = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
  149. ucEphXYZAdja1.SetParam("邻星1", adja1Code, (adja1eph.data.X, adja1eph.data.Y, adja1eph.data.Z, adja1eph.data.VX, adja1eph.data.VY, adja1eph.data.VZ), Color.Red);
  150. var adja2xlInfo = await XlRepository.GetLatestAsync(adja2Code, sigTime);
  151. if (adja2xlInfo == null)
  152. {
  153. DxHelper.MsgBoxHelper.ShowWarning($"卫星:{adja2Code}未找到对应的星历信息,请导入星历");
  154. return;
  155. }
  156. XlCalcDto = new XlCalcDto() { tleStr = adja2xlInfo.TwoLine, SigTime = sigTime };
  157. var adja2eph = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
  158. ucEphXYZAdja2.SetParam("邻星2", adja2Code, (adja2eph.data.X, adja2eph.data.Y, adja2eph.data.Z, adja2eph.data.VX, adja2eph.data.VY, adja2eph.data.VZ), Color.Red);
  159. }
  160. catch (Exception ex)
  161. {
  162. Serilog.Log.Error(ex, $"手动推算{PosResType.GetEnumDisplayName()}星历失败,SigTime={sigTime}");
  163. DxHelper.MsgBoxHelper.ShowError($"手动推算{PosResType.GetEnumDisplayName()}星历失败,{ex.Message}");
  164. }
  165. }
  166. }
  167. public class X3GDOP星地接口
  168. {
  169. /// <summary>
  170. /// 主星星历x y z vx vy vz
  171. /// </summary>
  172. public double[] EphMain { get; set; }
  173. /// <summary>
  174. /// 邻星1星历x y z vx vy vz
  175. /// </summary>
  176. public double[] EphAdja1 { get; set; }
  177. /// <summary>
  178. /// 邻星2星历x y z vx vy vz
  179. /// </summary>
  180. public double[] EphAdja2 { get; set; }
  181. /// <summary>
  182. /// 采集时刻
  183. /// </summary>
  184. public DateTime CapTime { get; set; }
  185. /// <summary>
  186. /// 参考站位置-经度 纬度 高度
  187. /// </summary>
  188. public double[] RefStation { get; set; }
  189. /// <summary>
  190. /// 时差误差(单位us)
  191. /// </summary>
  192. public double DtousErr { get; set; } = 1;
  193. /// <summary>
  194. /// 星历位置误差(单位米)
  195. /// </summary>
  196. public double SatLocErr { get; set; } = 10000;
  197. }
  198. }