RHGDOPParam.cs 8.5 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231
  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.Utils.About;
  13. using DevExpress.XtraEditors.DXErrorProvider;
  14. using XdCxRhDW.Dto;
  15. using System.Configuration;
  16. namespace XdCxRhDW.App.UserControl
  17. {
  18. public partial class RHGDOPParam : DevExpress.XtraEditors.XtraUserControl
  19. {
  20. public MapControl mapControl1;
  21. public GDOP融合接口 Model => new GDOP融合接口()
  22. {
  23. EphMain = ucEphXYZMain.EphParam(),
  24. EphAdja = ucEphXYZAdja.EphParam(),
  25. CapTime = txtCapTime.DateTime,
  26. CDBStation = txtStationLocation1.GetLonLat(),
  27. RefStation = txtRefLocation1.GetLonLat(),
  28. DtousErr = Convert.ToDouble(txtDtousErr1.Text),
  29. SatLocErr = Convert.ToDouble(txtSatLocErr1.Text),
  30. };
  31. private EnumPosResType PosResType;
  32. public RHGDOPParam(PosRes item)
  33. {
  34. InitializeComponent();
  35. this.layoutControl1.UseDefault();
  36. txtCapTime.UseDefault();
  37. txtStationLocation1.UseDoubleClickToSelectAll();
  38. txtRefLocation1.UseDoubleClickToSelectAll();
  39. this.txtCapTime.DateTime = item.SigTime;
  40. this.txtDtousErr1.EditValue = 1;
  41. this.txtSatLocErr1.EditValue = 10000;
  42. PosResType = item.PosResType;
  43. using (RHDWPartContext db = RHDWPartContext.GetContext(item.SigTime))
  44. {
  45. var cg = db?.CgRes.Where(m => m.ID == item.CgResID).FirstOrDefault();
  46. var station = db?.StationRes.Where(m => m.ID == item.StationResID).FirstOrDefault();
  47. if (station != null)
  48. {
  49. this.txtStationLocation1.Text = $"{station.CdbTxLon},{station.CdbTxLat}";
  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 adjaEph = (cg.Adja1X, cg.Adja1Y, cg.Adja1Z, cg.Adja1Vx, cg.Adja1Vy, cg.Adja1Vz);
  57. ucEphXYZAdja.SetParam($"邻星", cg.Adja1Code, adjaEph, Color.Black);
  58. }
  59. }
  60. }
  61. private bool ParamValidate()
  62. {
  63. dxErrorProvider.ClearErrors();
  64. if (!ucEphXYZMain.CheckEph(dxErrorProvider))
  65. {
  66. return false;
  67. }
  68. if (!ucEphXYZAdja.CheckEph(dxErrorProvider))
  69. {
  70. return false;
  71. }
  72. if (!txtStationLocation1.CheckLonLat(dxErrorProvider, "超短波"))
  73. {
  74. return false;
  75. }
  76. if (!txtRefLocation1.CheckLonLat(dxErrorProvider, "参考站"))
  77. {
  78. return false;
  79. }
  80. if (!txtDtousErr1.CheckDouble(dxErrorProvider, "时差误差"))
  81. {
  82. return false;
  83. }
  84. if (!txtSatLocErr1.CheckDouble(dxErrorProvider, "星历位置误差"))
  85. {
  86. return false;
  87. }
  88. return true;
  89. }
  90. private void btnOK_Click(object sender, EventArgs e)
  91. {
  92. if (!ParamValidate()) { return; }
  93. try
  94. {
  95. mapControl1.ClearMap();
  96. string useNewGDOP = ConfigurationManager.AppSettings["UseNewGDOP"];
  97. if (useNewGDOP == null || useNewGDOP == "0")
  98. {
  99. var data = GdopHelper.Gdop2Sat1DByXyz(Model.EphMain, Model.EphAdja, Model.CDBStation
  100. , Model.DtousErr, Model.SatLocErr, Model.RefStation);
  101. foreach (var errLins in data)//画GDOP
  102. {
  103. mapControl1.DrawGdopLineTwo(errLins.ErrDistanceKm, errLins.MapDots, 1);
  104. }
  105. }
  106. else
  107. {
  108. var data = GdopHelper.Gdop2Sat1DByXyzNew(Model.EphMain, Model.EphAdja, Model.CDBStation
  109. , Model.DtousErr, Model.SatLocErr, Model.RefStation);
  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. }
  124. catch (Exception ex)
  125. {
  126. DxHelper.MsgBoxHelper.ShowError($"绘制{PosResType.GetEnumDisplayName()}GDOP失败,{ex.Message}");
  127. }
  128. }
  129. private void btnClose_Click(object sender, EventArgs e)
  130. {
  131. DxHelper.PopupHelper.HidePopup(this);
  132. }
  133. private async void btnEphCalc_Click(object sender, EventArgs e)
  134. {
  135. dxErrorProvider.ClearErrors();
  136. if (txtCapTime.DateTime == DateTime.MinValue)
  137. {
  138. dxErrorProvider.SetError(txtCapTime, "信号时间不能为空!");
  139. return;
  140. }
  141. var sigTime = txtCapTime.DateTime;
  142. try
  143. {
  144. var mainCode = ucEphXYZMain.GetSatCode();
  145. var adjaCode = ucEphXYZAdja.GetSatCode();
  146. var mainxlInfo = await XlRepository.GetLatestAsync(mainCode, sigTime);
  147. if (mainxlInfo == null)
  148. {
  149. DxHelper.MsgBoxHelper.ShowWarning($"卫星:{mainCode}未找到对应的星历信息,请导入星历");
  150. return;
  151. }
  152. var XlCalcDto = new XlCalcDto() { tleStr = mainxlInfo.TwoLine, SigTime = sigTime };
  153. var maineph = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
  154. ucEphXYZMain.SetParam("主星", mainCode, (maineph.data.X, maineph.data.Y, maineph.data.Z, maineph.data.VX, maineph.data.VY, maineph.data.VZ), Color.Red);
  155. var adjaxlInfo = await XlRepository.GetLatestAsync(adjaCode, sigTime);
  156. if (adjaxlInfo == null)
  157. {
  158. DxHelper.MsgBoxHelper.ShowWarning($"卫星:{adjaCode}未找到对应的星历信息,请导入星历");
  159. return;
  160. }
  161. XlCalcDto = new XlCalcDto() { tleStr = adjaxlInfo.TwoLine, SigTime = sigTime };
  162. var adjaeph = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
  163. ucEphXYZAdja.SetParam("邻星", adjaCode, (adjaeph.data.X, adjaeph.data.Y, adjaeph.data.Z, adjaeph.data.VX, adjaeph.data.VY, adjaeph.data.VZ), Color.Red);
  164. }
  165. catch (Exception ex)
  166. {
  167. Serilog.Log.Error(ex, $"手动推算{PosResType.GetEnumDisplayName()}星历失败,SigTime={sigTime}");
  168. DxHelper.MsgBoxHelper.ShowError($"手动推算{PosResType.GetEnumDisplayName()}星历失败,{ex.Message}");
  169. }
  170. }
  171. }
  172. public class GDOP融合接口
  173. {
  174. /// <summary>
  175. /// 主星x y z vx vy vz
  176. /// </summary>
  177. public double[] EphMain { get; set; }
  178. /// <summary>
  179. /// 邻星x y z vx vy vz
  180. /// </summary>
  181. public double[] EphAdja { get; set; }
  182. /// <summary>
  183. /// 采集时刻
  184. /// </summary>
  185. public DateTime CapTime { get; set; }
  186. /// <summary>
  187. /// 超短接收站-经度 纬度 高度
  188. /// </summary>
  189. public double[] CDBStation { get; set; }
  190. /// <summary>
  191. /// 参考站位置-经度 纬度 高度
  192. /// </summary>
  193. public double[] RefStation { get; set; }
  194. /// <summary>
  195. /// 时差误差(单位us)
  196. /// </summary>
  197. public double DtousErr { get; set; } = 1;
  198. /// <summary>
  199. /// 星历位置误差(单位米)
  200. /// </summary>
  201. public double SatLocErr { get; set; } = 10000;
  202. /// <summary>
  203. /// 定位类型
  204. /// </summary>
  205. public EnumPosResType PosResType { get; set; }
  206. }
  207. }