X3GDOPParam.cs 9.0 KB

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