X3DFGDOPParam .cs 7.7 KB

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