X3DFGDOPParam .cs 9.8 KB

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