X2DFGDOPParam.cs 9.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246
  1. 
  2. using DevExpress.XtraEditors.DXErrorProvider;
  3. using DevExpress.XtraMap;
  4. using DxHelper;
  5. using ExtensionsDev;
  6. using System;
  7. using System.Collections.Generic;
  8. using System.Configuration;
  9. using System.Data;
  10. using System.Drawing;
  11. using System.Linq;
  12. using XdCxRhDW.Api;
  13. using XdCxRhDW.Dto;
  14. using XdCxRhDW.Entity;
  15. using XdCxRhDW.Framework;
  16. using XdCxRhDW.Repostory;
  17. namespace XdCxRhDW.App.UserControl
  18. {
  19. public partial class X2DFGDOPParam : DevExpress.XtraEditors.XtraUserControl
  20. {
  21. public MapControl mapControl1;
  22. public X2DFGDOP接口 Model => new X2DFGDOP接口()
  23. {
  24. EphMain = ucEphXYZMain.EphParam(),
  25. EphAdja = ucEphXYZAdja.EphParam(),
  26. CapTime = txtCapTime.DateTime,
  27. RefStation = txtRefLocation1.GetLonLat(),
  28. DfoErr = Convert.ToDouble(txtDfoErr1.Text),
  29. DtousErr = Convert.ToDouble(txtDtousErr1.Text),
  30. SatLocErr = Convert.ToDouble(txtSatLocErr1.Text),
  31. EphVelErr = Convert.ToDouble(txtEphVelErr1.Text),
  32. fu1 = Convert.ToDouble(txtFu1.Text) * 1e6,
  33. fu2 = Convert.ToDouble(txtFu2.Text) * 1e6,
  34. };
  35. private EnumPosResType PosResType;
  36. public X2DFGDOPParam(PosRes item)
  37. {
  38. InitializeComponent();
  39. this.layoutControl1.UseDefault();
  40. txtCapTime.UseDefault();
  41. txtRefLocation1.UseDoubleClickToSelectAll();
  42. this.txtCapTime.DateTime = item.SigTime;
  43. this.txtSatLocErr1.EditValue = 10000;
  44. this.txtEphVelErr1.EditValue = 0.1;
  45. this.txtDfoErr1.EditValue = 0.01;
  46. this.txtDtousErr1.EditValue = 1;
  47. PosResType = item.PosResType;
  48. using (RHDWPartContext db = RHDWPartContext.GetContext(item.SigTime))
  49. {
  50. var cg = db?.CgRes.Where(m => m.ID == item.CgResID).FirstOrDefault();
  51. var station = db?.StationRes.Where(m => m.ID == item.StationResID).FirstOrDefault();
  52. if (station != null)
  53. {
  54. this.txtRefLocation1.Text = $"{station.RefLon},{station.RefLat}";
  55. }
  56. if (cg != null)
  57. {
  58. this.txtFu1.EditValue = cg.TarFreqUp.HasValue ? cg.TarFreqUp.Value * 1e-6 : 0;
  59. this.txtFu2.EditValue = cg.RefFreqUp.HasValue ? cg.RefFreqUp.Value * 1e-6 : 0;
  60. var mainEph = (cg.MainX, cg.MainY, cg.MainZ, cg.MainVx, cg.MainVy, cg.MainVz);
  61. ucEphXYZMain.SetParam($"主星", cg.MainCode, mainEph, Color.Black);
  62. var adjaEph = (cg.Adja1X, cg.Adja1Y, cg.Adja1Z, cg.Adja1Vx, cg.Adja1Vy, cg.Adja1Vz);
  63. ucEphXYZAdja.SetParam($"邻星", cg.Adja1Code, adjaEph, Color.Black);
  64. }
  65. }
  66. }
  67. private bool ParamValidate()
  68. {
  69. dxErrorProvider.ClearErrors();
  70. if (!ucEphXYZMain.CheckEph(dxErrorProvider))
  71. {
  72. return false;
  73. }
  74. if (!ucEphXYZAdja.CheckEph(dxErrorProvider))
  75. {
  76. return false;
  77. }
  78. if (!txtRefLocation1.CheckLonLat(dxErrorProvider, "参考站"))
  79. {
  80. return false;
  81. }
  82. if (!txtDtousErr1.CheckDouble(dxErrorProvider, "时差误差"))
  83. {
  84. return false;
  85. }
  86. if (!txtDfoErr1.CheckDouble(dxErrorProvider, "频差误差"))
  87. {
  88. return false;
  89. }
  90. if (!txtSatLocErr1.CheckDouble(dxErrorProvider, "星历位置误差"))
  91. {
  92. return false;
  93. }
  94. if (!txtEphVelErr1.CheckDouble(dxErrorProvider, "星历速度误差"))
  95. {
  96. return false;
  97. }
  98. if (!txtFu1.CheckDouble(dxErrorProvider, "目标上行频点"))
  99. {
  100. return false;
  101. }
  102. if (!txtFu2.CheckDouble(dxErrorProvider, "参考上行频点"))
  103. {
  104. return false;
  105. }
  106. return true;
  107. }
  108. private void btnOK_Click(object sender, EventArgs e)
  109. {
  110. if (!ParamValidate()) { return; }
  111. try
  112. {
  113. mapControl1.ClearMap();
  114. var data = GdopHelper.Gdop2SatDRefByXyzNew(Model.EphMain, Model.EphAdja, Model.fu1, Model.fu2, Model.DtousErr, Model.DfoErr, Model.SatLocErr, Model.EphVelErr, Model.RefStation);
  115. if (data.Count == 0)
  116. {
  117. DxHelper.MsgBoxHelper.ShowInfo($"绘制{PosResType.GetEnumDisplayName()}GDOP没有{GdopParam.误差配置.误差距离km.Max()}公里以内的误差!");
  118. return;
  119. }
  120. List<MapPolyline> polylines = new List<MapPolyline>();
  121. for (int i = 0; i < data.Count; i++)
  122. {
  123. var errLins = data[i];
  124. var mapline = mapControl1.GetLine(errLins.ErrDistanceKm, errLins.MapDots);
  125. if (mapline == null)
  126. {
  127. continue;
  128. }
  129. polylines.Add(mapline);
  130. }
  131. mapControl1.DrawDtoLine(polylines);
  132. }
  133. catch (Exception ex)
  134. {
  135. DxHelper.MsgBoxHelper.ShowError($"绘制{PosResType.GetEnumDisplayName()}GDOP失败,{ex.Message}");
  136. }
  137. }
  138. private void btnClose_Click(object sender, EventArgs e)
  139. {
  140. PopupHelper.HidePopup(this);
  141. }
  142. private async void btnEphCalc_Click(object sender, EventArgs e)
  143. {
  144. dxErrorProvider.ClearErrors();
  145. if (txtCapTime.DateTime == DateTime.MinValue)
  146. {
  147. dxErrorProvider.SetError(txtCapTime, "信号时间不能为空!");
  148. return;
  149. }
  150. var sigTime = txtCapTime.DateTime;
  151. try
  152. {
  153. var mainCode = ucEphXYZMain.GetSatCode();
  154. var adjaCode = ucEphXYZAdja.GetSatCode();
  155. var mainxlInfo = await XlRepository.GetLatestAsync(mainCode, sigTime);
  156. if (mainxlInfo == null)
  157. {
  158. DxHelper.MsgBoxHelper.ShowWarning($"卫星:{mainCode}未找到对应的星历信息,请导入星历");
  159. return;
  160. }
  161. var XlCalcDto = new XlCalcDto() { tleStr = mainxlInfo.TwoLine, SigTime = sigTime };
  162. var maineph = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
  163. ucEphXYZMain.SetParam("主星", mainCode, (maineph.data.X, maineph.data.Y, maineph.data.Z, maineph.data.VX, maineph.data.VY, maineph.data.VZ), Color.Red);
  164. var adjaxlInfo = await XlRepository.GetLatestAsync(adjaCode, sigTime);
  165. if (adjaxlInfo == null)
  166. {
  167. DxHelper.MsgBoxHelper.ShowWarning($"卫星:{adjaCode}未找到对应的星历信息,请导入星历");
  168. return;
  169. }
  170. XlCalcDto = new XlCalcDto() { tleStr = adjaxlInfo.TwoLine, SigTime = sigTime };
  171. var adjaeph = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
  172. ucEphXYZAdja.SetParam("邻星", adjaCode, (adjaeph.data.X, adjaeph.data.Y, adjaeph.data.Z, adjaeph.data.VX, adjaeph.data.VY, adjaeph.data.VZ), Color.Red);
  173. }
  174. catch (Exception ex)
  175. {
  176. Serilog.Log.Error(ex, $"手动推算{PosResType.GetEnumDisplayName()}星历失败,SigTime={sigTime}");
  177. DxHelper.MsgBoxHelper.ShowError($"手动推算{PosResType.GetEnumDisplayName()}星历失败,{ex.Message}");
  178. }
  179. }
  180. }
  181. public class X2DFGDOP接口
  182. {
  183. /// <summary>
  184. /// 主星星历x y z vx vy vz
  185. /// </summary>
  186. public double[] EphMain { get; set; }
  187. /// <summary>
  188. /// 邻星1星历x y z vx vy vz
  189. /// </summary>
  190. public double[] EphAdja { 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. /// 频差误差(Hz)
  205. /// </summary>
  206. public double DfoErr { get; set; }
  207. /// <summary>
  208. /// 星历位置误差(单位米)
  209. /// </summary>
  210. public double SatLocErr { get; set; } = 10000;
  211. /// <summary>
  212. ///星历速度误差
  213. /// </summary>
  214. public double EphVelErr { get; set; }
  215. /// <summary>
  216. /// 上行频点1(Hz)
  217. /// </summary>
  218. public double fu1 { get; set; }
  219. /// <summary>
  220. /// 上行频点2(Hz)
  221. /// </summary>
  222. public double fu2 { get; set; }
  223. }
  224. }