X2DFGDOPParam.cs 7.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207
  1. 
  2. using DevExpress.XtraMap;
  3. using DxHelper;
  4. using ExtensionsDev;
  5. using System;
  6. using System.Collections.Generic;
  7. using System.Data;
  8. using System.Linq;
  9. using XdCxRhDW.Core.Api;
  10. using XdCxRhDW.Repostory.EFContext;
  11. using XdCxRhDW.Repostory.Model;
  12. namespace XdCxRhDW.App.UserControl
  13. {
  14. public partial class X2DFGDOPParam : DevExpress.XtraEditors.XtraUserControl
  15. {
  16. public MapControl mapControl1;
  17. public X2DFGDOP接口 Model => new X2DFGDOP接口()
  18. {
  19. TleLeo1 = txtTleLeo1.Text.Trim(),
  20. TleLeo2 = txtTleLeo2.Text.Trim(),
  21. CapTime = txtCapTime.DateTime,
  22. RefLon = Convert.ToDouble(txtRefLocation1.Text.Replace(",", ",").Split(',')[0].Trim()),
  23. RefLat = Convert.ToDouble(txtRefLocation1.Text.Replace(",", ",").Split(',')[1].Trim()),
  24. DfoErr = Convert.ToDouble(txtDfoErr1.Text),
  25. DtousErr = Convert.ToDouble(txtDtousErr1.Text),
  26. SatLocErr = Convert.ToDouble(txtSatLocErr1.Text),
  27. EphVelErr = Convert.ToDouble(txtEphVelErr1.Text),
  28. fu1 = Convert.ToDouble(txtFu1.Text) * 1e6,
  29. fu2 = Convert.ToDouble(txtFu2.Text) * 1e6,
  30. };
  31. public X2DFGDOPParam(PosRes item)
  32. {
  33. InitializeComponent();
  34. txtCapTime.UseDefault();
  35. txtTleLeo1.UseDoubleClickToSelectAll();
  36. txtTleLeo2.UseDoubleClickToSelectAll();
  37. txtRefLocation1.UseDoubleClickToSelectAll();
  38. this.txtCapTime.DateTime = item.SigTime;
  39. this.txtSatLocErr1.EditValue = 10000;
  40. this.txtEphVelErr1.EditValue = 0.1;
  41. this.txtDfoErr1.EditValue = 0.01;
  42. this.txtDtousErr1.EditValue = 1;
  43. List<string> xlall = new List<string>();
  44. string mainTle = string.Empty;
  45. string adjaTle1 = string.Empty;
  46. using (RHDWContext db = new RHDWContext())
  47. {
  48. var sats = db.SatInfos.ToList();
  49. var cg = db.CgRes.Where(m => m.ID == item.CgResID).FirstOrDefault();
  50. var station = db.StationRes.Where(m => m.ID == item.StationResID).FirstOrDefault();
  51. if (station != null)
  52. {
  53. this.txtRefLocation1.Text = $"{station.RefLon},{station.RefLat}";
  54. }
  55. this.txtFu1.EditValue = cg.TarFreqUp.HasValue ? cg.TarFreqUp.Value * 1e-6 : 950;
  56. this.txtFu2.EditValue = cg.RefFreqUp.HasValue ? cg.RefFreqUp.Value * 1e-6 : 950;
  57. var xlList = db.XlInfos.OrderBy(p => p.SatName).OrderByDescending(p => p.TimeBJ).ToList();
  58. xlall.AddRange(xlList.Select(m => m.TwoLine));
  59. if (xlall.Count() == 0) return;
  60. if (cg.MainCode.HasValue && xlList.Any(m => m.SatCode == cg.MainCode.Value))
  61. {
  62. mainTle = xlList.First(m => m.SatCode == cg.MainCode.Value).TwoLine;
  63. }
  64. else
  65. {
  66. mainTle = xlList.First().TwoLine;
  67. }
  68. if (cg.Adja1Code.HasValue && xlList.Any(m => m.SatCode == cg.Adja1Code.Value))
  69. {
  70. adjaTle1 = xlList.First(m => m.SatCode == cg.Adja1Code.Value).TwoLine;
  71. }
  72. else
  73. {
  74. adjaTle1 = xlList.First().TwoLine;
  75. }
  76. }
  77. txtTleLeo1.UseDefault().SetStringData(xlall).Text = mainTle;
  78. txtTleLeo2.UseDefault().SetStringData(xlall).Text = adjaTle1;
  79. }
  80. private (bool, string) ParamValidate()
  81. {
  82. if (string.IsNullOrWhiteSpace(txtTleLeo1.Text.Trim()))
  83. {
  84. return (false, "主星星历不能为空!");
  85. }
  86. if (string.IsNullOrWhiteSpace(txtTleLeo2.Text.Trim()))
  87. {
  88. return (false, "邻星星历不能为空!");
  89. }
  90. if (txtTleLeo1.Text.Trim() == txtTleLeo2.Text.Trim())
  91. {
  92. return (false, "主邻星历不能相同!");
  93. }
  94. if (txtCapTime.DateTime == DateTime.MinValue)
  95. {
  96. return (false, "采集时刻不能为空!");
  97. }
  98. var refsta = txtRefLocation1.CheckLonLat("参考站");
  99. if ( !refsta.Item1)
  100. {
  101. return refsta;
  102. }
  103. return (true, "");
  104. }
  105. private void btnOK_Click(object sender, EventArgs e)
  106. {
  107. var pv = ParamValidate();
  108. if (!pv.Item1)
  109. {
  110. DxHelper.MsgBoxHelper.ShowWarning($"{pv.Item2}");
  111. return;
  112. }
  113. try
  114. {
  115. mapControl1.ClearMap();
  116. var refs = new double[] { Model.RefLon, Model.RefLat, 0 };
  117. var (listSat, data) = GdopHelper.Gdop2SatDRef(Model.TleLeo1, Model.TleLeo2,
  118. Model.CapTime, Model.fu1, Model.fu2, Model.DtousErr, Model.DfoErr, Model.SatLocErr, Model.EphVelErr, refs);
  119. if (data == null)
  120. {
  121. return;
  122. }
  123. foreach (var errLins in data)//画GDOP
  124. {
  125. var mapDots = errLins.MapDots.Select(p => p).Select(p => (p.Lon, p.Lat));
  126. mapControl1.DrawGdopLineTwo(errLins.ErrDistanceKm, mapDots, mapDots.Count() / 2);
  127. }
  128. }
  129. catch (Exception ex)
  130. {
  131. DxHelper.MsgBoxHelper.ShowError($"绘制GDOP失败,{ex.Message}");
  132. }
  133. }
  134. private void btnClose_Click(object sender, EventArgs e)
  135. {
  136. PopupHelper.HidePopup(this);
  137. }
  138. }
  139. public class X2DFGDOP接口
  140. {
  141. /// <summary>
  142. /// 主星星历(Tle)
  143. /// </summary>
  144. public string TleLeo1 { get; set; }
  145. /// <summary>
  146. /// 邻星1星历(Tle)
  147. /// </summary>
  148. public string TleLeo2 { get; set; }
  149. /// <summary>
  150. /// 采集时刻
  151. /// </summary>
  152. public DateTime CapTime { get; set; }
  153. /// <summary>
  154. /// 参考站位置经度
  155. /// </summary>
  156. public double RefLon { get; set; }
  157. /// <summary>
  158. /// 参考站位置纬度
  159. /// </summary>
  160. public double RefLat { get; set; }
  161. /// <summary>
  162. /// 时差误差(单位us)
  163. /// </summary>
  164. public double DtousErr { get; set; } = 1;
  165. /// <summary>
  166. /// 频差误差(Hz)
  167. /// </summary>
  168. public double DfoErr { get; set; }
  169. /// <summary>
  170. /// 星历位置误差(单位米)
  171. /// </summary>
  172. public double SatLocErr { get; set; } = 10000;
  173. /// <summary>
  174. ///星历速度误差
  175. /// </summary>
  176. public double EphVelErr { get; set; }
  177. /// <summary>
  178. /// 上行频点1(Hz)
  179. /// </summary>
  180. public double fu1 { get; set; }
  181. /// <summary>
  182. /// 上行频点2(Hz)
  183. /// </summary>
  184. public double fu2 { get; set; }
  185. }
  186. }