X3DFGDOPParam .cs 9.9 KB

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