X3DfoParamEditor.cs 19 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440
  1. using DevExpress.XtraMap;
  2. using ExtensionsDev;
  3. using System;
  4. using System.Collections.Generic;
  5. using System.Data.Entity;
  6. using System.Drawing;
  7. using System.Linq;
  8. using XdCxRhDW.Api;
  9. using XdCxRhDW.Dto;
  10. using XdCxRhDW.Entity;
  11. using XdCxRhDW.Repostory;
  12. namespace XdCxRhDW.App.EditForms
  13. {
  14. public partial class X3DfoParamEditor : DevExpress.XtraEditors.XtraUserControl
  15. {
  16. private PosRes info;
  17. private CgRes cg;
  18. private MapControl mapControl1;
  19. private List<SatInfo> listSat;
  20. StationRes station;
  21. public X3DfoParamEditor(PosRes info, MapControl mapControl)
  22. {
  23. InitializeComponent();
  24. this.info = info;
  25. this.mapControl1 = mapControl;
  26. itemSigTime.Text = $"{itemSigTime.Text}({SysConfig.Config.TimeZoneUTC})";
  27. this.layoutControl1.UseDefault();
  28. sigTime.UseDefault();
  29. this.simpleLabelItem1.Text = $" ";
  30. }
  31. private async void X3DTOParamEditor_Load(object sender, EventArgs e)
  32. {
  33. using (RHDWContext db = new RHDWContext())
  34. {
  35. listSat = await db.SatInfos.ToListAsync();
  36. }
  37. using (RHDWPartContext db = RHDWPartContext.GetContext(info.SigTime))
  38. {
  39. cg = await db?.CgRes.Where(m => m.ID == info.CgResID).FirstOrDefaultAsync();
  40. station = await db?.StationRes.Where(m => m.ID == info.StationResID).FirstOrDefaultAsync();
  41. }
  42. if (cg != null)
  43. {
  44. this.txtDfo1.Text = $"{cg.Dfo1.Value}";
  45. this.txtDfo2.Text = $"{cg.Dfo2.Value}";
  46. this.txtYbDfo1.Text = $"{cg.YbMainDfo.Value}";
  47. this.txtYbDfo2.Text = $"{cg.YbAdja1Dfo.Value}";
  48. this.txtYbDfo3.Text = $"{cg.YbAdja2Dfo.Value}";
  49. this.sigTime.EditValue = info.SigTime;
  50. var mainEph = (cg.MainX, cg.MainY, cg.MainZ, cg.MainVx, cg.MainVy, cg.MainVz);
  51. ucEphXYZMain.SetParam($"主星", cg.MainCode, mainEph, Color.Black);
  52. var adja1Eph = (cg.Adja1X, cg.Adja1Y, cg.Adja1Z, cg.Adja1Vx, cg.Adja1Vy, cg.Adja1Vz);
  53. ucEphXYZAdja1.SetParam($"邻星1", cg.Adja1Code, adja1Eph, Color.Black);
  54. var adja2Eph = (cg.Adja2X, cg.Adja2Y, cg.Adja2Z, cg.Adja2Vx, cg.Adja2Vy, cg.Adja2Vz);
  55. ucEphXYZAdja2.SetParam($"邻星2", cg.Adja2Code, adja2Eph, Color.Black);
  56. this.txtTargetFreq.EditValue = cg.TarFreqUp.HasValue ? cg.TarFreqUp.Value * 1e-6 : 0;
  57. this.txtTargetDFreq.EditValue = cg.TarFreqDown.HasValue ? cg.TarFreqDown.Value * 1e-6 : 0;
  58. this.txtRefFreq.EditValue = cg.RefFreqUp.HasValue ? cg.RefFreqUp.Value * 1e-6 : 0;
  59. this.txtRefDFreq.EditValue = cg.RefFreqDown.HasValue ? cg.RefFreqDown.Value * 1e-6 : 0;
  60. this.txtDto1.Text = cg.Dto1?.ToString();
  61. this.txtDto2.Text = cg.Dto2?.ToString();
  62. this.txtYbDto1.Text = cg.YbMainDto?.ToString();
  63. this.txtYbDto2.Text = cg.YbAdja1Dto?.ToString();
  64. this.txtYbDto3.Text = cg.YbAdja2Dto?.ToString();
  65. }
  66. if (station != null)
  67. {
  68. this.txtsatStation.Text = $"{station.SatTxLon},{station.SatTxLat}";
  69. this.txtRefLocation.Text = $"{station.RefLon},{station.RefLat}";
  70. }
  71. if (info.PosResType == EnumPosResType.X3NoRef)
  72. {
  73. txtYbDfo1.Properties.ReadOnly = true;
  74. txtYbDfo2.Properties.ReadOnly = true;
  75. txtYbDfo3.Properties.ReadOnly = true;
  76. txtRefLocation.Properties.ReadOnly = true;
  77. this.txtRefLocation.Text = $"{0},{0}";
  78. }
  79. }
  80. public bool CheckDtoLineParam()
  81. {
  82. dxErrorProvider.ClearErrors();
  83. if (!txtsatStation.CheckLonLat(dxErrorProvider, "接收站"))
  84. {
  85. return false;
  86. }
  87. if (!txtRefLocation.CheckLonLat(dxErrorProvider, "参考站"))
  88. {
  89. return false;
  90. }
  91. if (!txtDto1.CheckDouble(dxErrorProvider, "主邻1时差"))
  92. {
  93. return false;
  94. }
  95. if (!txtDto2.CheckDouble(dxErrorProvider, "主邻2时差"))
  96. {
  97. return false;
  98. }
  99. if (!txtYbDto1.CheckDouble(dxErrorProvider, "样本主星时差"))
  100. {
  101. return false;
  102. }
  103. if (!txtYbDto2.CheckDouble(dxErrorProvider, "样本邻星1时差"))
  104. {
  105. return false;
  106. }
  107. if (!txtYbDto3.CheckDouble(dxErrorProvider, "样本邻星2时差"))
  108. {
  109. return false;
  110. }
  111. if (!ucEphXYZMain.CheckEphXYZ(dxErrorProvider))
  112. {
  113. return false;
  114. }
  115. if (!ucEphXYZAdja1.CheckEphXYZ(dxErrorProvider))
  116. {
  117. return false;
  118. }
  119. if (!ucEphXYZAdja2.CheckEphXYZ(dxErrorProvider))
  120. {
  121. return false;
  122. }
  123. return true;
  124. }
  125. public bool CheckPosParam()
  126. {
  127. dxErrorProvider.ClearErrors();
  128. if (!txtsatStation.CheckLonLat(dxErrorProvider, "接收站"))
  129. {
  130. return false;
  131. }
  132. if (!txtRefLocation.CheckLonLat(dxErrorProvider, "参考站"))
  133. {
  134. return false;
  135. }
  136. if (!txtDfo1.CheckDouble(dxErrorProvider, "主邻1频差"))
  137. {
  138. return false;
  139. }
  140. if (!txtDfo2.CheckDouble(dxErrorProvider, "主邻2频差"))
  141. {
  142. return false;
  143. }
  144. if (!txtYbDfo1.CheckDouble(dxErrorProvider, "样本主星频差"))
  145. {
  146. return false;
  147. }
  148. if (!txtYbDfo2.CheckDouble(dxErrorProvider, "样本邻星1频差"))
  149. {
  150. return false;
  151. }
  152. if (!txtYbDfo3.CheckDouble(dxErrorProvider, "样本邻星2频差"))
  153. {
  154. return false;
  155. }
  156. if (!txtTargetFreq.CheckDouble(dxErrorProvider, "目标上行频点"))
  157. {
  158. return false;
  159. }
  160. if (!txtTargetDFreq.CheckDouble(dxErrorProvider, "目标下行频点"))
  161. {
  162. return false;
  163. }
  164. if (!txtRefFreq.CheckDouble(dxErrorProvider, "参考上行频点"))
  165. {
  166. return false;
  167. }
  168. if (!txtRefDFreq.CheckDouble(dxErrorProvider, "参考下行频点"))
  169. {
  170. return false;
  171. }
  172. if (!ucEphXYZMain.CheckEph(dxErrorProvider))
  173. {
  174. return false;
  175. }
  176. if (!ucEphXYZAdja1.CheckEph(dxErrorProvider))
  177. {
  178. return false;
  179. }
  180. if (!ucEphXYZAdja2.CheckEph(dxErrorProvider))
  181. {
  182. return false;
  183. }
  184. return true;
  185. }
  186. private void btnOk_Click(object sender, EventArgs e)
  187. {
  188. if (!CheckPosParam()) { return; }
  189. try
  190. {
  191. var MsAnt = txtsatStation.GetLonLat();
  192. var RefGeod = txtRefLocation.GetLonLat();
  193. var sigTime = this.sigTime.DateTime;
  194. var DfoSx = Convert.ToDouble(this.txtDfo1.Text);
  195. var DfoSx1 = Convert.ToDouble(this.txtDfo2.Text);
  196. var YbMainDfo = Convert.ToDouble(this.txtYbDfo1.Text);
  197. var YbAdja1Dfo = Convert.ToDouble(this.txtYbDfo2.Text);
  198. var YbAdja2Dfo = Convert.ToDouble(this.txtYbDfo3.Text);
  199. var tFreq = Convert.ToDouble(this.txtTargetFreq.Text) * 1e6;
  200. var tdFreq = Convert.ToDouble(this.txtTargetDFreq.Text) * 1e6;
  201. var rFreq = Convert.ToDouble(this.txtRefFreq.Text) * 1e6;
  202. var rdFreq = Convert.ToDouble(this.txtRefDFreq.Text) * 1e6;
  203. var msEph = ucEphXYZMain.EphParam();
  204. var ns1Eph = ucEphXYZAdja1.EphParam();
  205. var ns2Eph = ucEphXYZAdja2.EphParam();
  206. var StationRes = new StationRes()
  207. {
  208. SatTxLon = MsAnt[0],
  209. SatTxLat = MsAnt[1],
  210. RefLon = RefGeod[0],
  211. RefLat = RefGeod[1],
  212. };
  213. var cgRes = new CgRes()
  214. {
  215. SigTime = sigTime,
  216. Dfo1 = DfoSx,
  217. Dfo2 = DfoSx1,
  218. YbMainDfo = YbMainDfo,
  219. YbAdja1Dfo = YbAdja1Dfo,
  220. YbAdja2Dfo = YbAdja2Dfo,
  221. TarFreqUp = tFreq,
  222. TarFreqDown = tdFreq,
  223. RefFreqUp = rFreq,
  224. RefFreqDown = rdFreq,
  225. MainX = msEph[0],
  226. MainY = msEph[1],
  227. MainZ = msEph[2],
  228. MainVx = msEph[3],
  229. MainVy = msEph[4],
  230. MainVz = msEph[5],
  231. Adja1X = ns1Eph[0],
  232. Adja1Y = ns1Eph[1],
  233. Adja1Z = ns1Eph[2],
  234. Adja1Vx = ns1Eph[3],
  235. Adja1Vy = ns1Eph[4],
  236. Adja1Vz = ns1Eph[5],
  237. Adja2X = ns2Eph[0],
  238. Adja2Y = ns2Eph[1],
  239. Adja2Z = ns2Eph[2],
  240. Adja2Vx = ns2Eph[3],
  241. Adja2Vy = ns2Eph[4],
  242. Adja2Vz = ns2Eph[5],
  243. };
  244. var res = PosApi.X3_PosTwoDfo(cgRes, StationRes);
  245. this.simpleLabelItem1.Text = $"{info.PosResType.GetEnumDisplayName()}定位点:[{res[0]:f4},{res[1]:f4}] 镜像点:[{res[3]:f4},{res[4]:f4}]";
  246. }
  247. catch (Exception ex)
  248. {
  249. Serilog.Log.Error(ex, $"{info.PosResType.GetEnumDisplayName()}手动定位失败.PosID={info.ID},SigTime={info.SigTime}");
  250. DxHelper.MsgBoxHelper.ShowWarning($"{info.PosResType.GetEnumDisplayName()}手动定位失败,{ex.Message}");
  251. }
  252. }
  253. private async void btnXl_Click(object sender, EventArgs e)
  254. {
  255. dxErrorProvider.ClearErrors();
  256. if (this.sigTime.DateTime == DateTime.MinValue)
  257. {
  258. dxErrorProvider.SetError(this.sigTime, "信号时间不能为空!");
  259. return;
  260. }
  261. var sigTime = this.sigTime.DateTime;
  262. try
  263. {
  264. var xlInfo = await XlRepository.GetLatestAsync(cg.MainCode.Value, sigTime);
  265. if (xlInfo == null)
  266. {
  267. DxHelper.MsgBoxHelper.ShowWarning($"卫星:{cg.MainCode.Value}未找到对应的星历信息,请导入星历");
  268. return;
  269. }
  270. var xlInfo1 = await XlRepository.GetLatestAsync(cg.Adja1Code.Value, sigTime);
  271. if (xlInfo1 == null)
  272. {
  273. DxHelper.MsgBoxHelper.ShowWarning($"卫星:{cg.Adja1Code.Value}未找到对应的星历信息,请导入星历");
  274. return;
  275. }
  276. var xlInfo2 = await XlRepository.GetLatestAsync(cg.Adja2Code.Value, sigTime);
  277. if (xlInfo2 == null)
  278. {
  279. DxHelper.MsgBoxHelper.ShowWarning($"卫星:{cg.Adja2Code.Value}未找到对应的星历信息,请导入星历");
  280. return;
  281. }
  282. var XlCalcDto = new XlCalcDto() { tleStr = xlInfo.TwoLine, SigTime = sigTime };
  283. var mEph = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
  284. XlCalcDto.tleStr = xlInfo1.TwoLine;
  285. var nEph1 = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
  286. XlCalcDto.tleStr = xlInfo2.TwoLine;
  287. var nEph2 = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
  288. var mainEph = (mEph.data.X, mEph.data.Y, mEph.data.Z, mEph.data.VX, mEph.data.VY, mEph.data.VZ);
  289. ucEphXYZMain.SetParam($"主星", cg.MainCode, mainEph, Color.Red);
  290. var adja1Eph = (nEph1.data.X, nEph1.data.Y, nEph1.data.Z, nEph1.data.VX, nEph1.data.VY, nEph1.data.VZ);
  291. ucEphXYZAdja1.SetParam($"邻星1", cg.Adja1Code, adja1Eph, Color.Red);
  292. var adja2Eph = (nEph2.data.X, nEph2.data.Y, nEph2.data.Z, nEph2.data.VX, nEph2.data.VY, nEph2.data.VZ);
  293. ucEphXYZAdja2.SetParam($"邻星2", cg.Adja2Code, adja2Eph, Color.Red);
  294. }
  295. catch (Exception ex)
  296. {
  297. Serilog.Log.Error(ex, $"手动{info.PosResType.GetEnumDisplayName()}推算星历失败.PosID={info.ID},SigTime={info.SigTime}");
  298. DxHelper.MsgBoxHelper.ShowWarning($"手动{info.PosResType.GetEnumDisplayName()}推算星历失败,{ex.Message}");
  299. }
  300. }
  301. private void btnDtoCalc_Click(object sender, EventArgs e)
  302. {
  303. if (!CheckDtoLineParam()) { return; }
  304. try
  305. {
  306. var MsAnt = txtsatStation.GetLonLat();
  307. var RefGeod = txtRefLocation.GetLonLat();
  308. var DtoSx = Convert.ToDouble(this.txtDto1.Text);
  309. var DtoSx1 = Convert.ToDouble(this.txtDto2.Text);
  310. var YbMainDto = Convert.ToDouble(this.txtYbDto1.Text);
  311. var YbAdja1Dto = Convert.ToDouble(this.txtYbDto2.Text);
  312. var YbAdja2Dto = Convert.ToDouble(this.txtYbDto3.Text);
  313. double[] msEph = ucEphXYZMain.EphXYZ();
  314. double[] Ns1Eph = ucEphXYZAdja1.EphXYZ();
  315. double[] Ns2Eph = ucEphXYZAdja2.EphXYZ();
  316. DtoLineTwoStartOption twoStartOption = new DtoLineTwoStartOption();
  317. twoStartOption.MsEph = msEph;
  318. twoStartOption.NsEph = Ns1Eph;
  319. twoStartOption.MsAnt = MsAnt;
  320. twoStartOption.NsAnt = MsAnt;
  321. twoStartOption.RefGeod = RefGeod;
  322. twoStartOption.TargetDto = DtoSx;
  323. twoStartOption.RefDto = YbMainDto - YbAdja1Dto;
  324. twoStartOption.PosLon = info.PosLon;
  325. twoStartOption.PosLat = info.PosLat;
  326. var msat = listSat.FirstOrDefault(m => m.SatCode == cg.MainCode.Value)?.Sat;
  327. if (string.IsNullOrWhiteSpace(msat)) msat = cg.MainCode.Value.ToString();
  328. var nsat = listSat.FirstOrDefault(m => m.SatCode == cg.Adja1Code.Value)?.Sat;
  329. if (string.IsNullOrWhiteSpace(nsat)) nsat = cg.Adja1Code.Value.ToString();
  330. var ad2sat = listSat.FirstOrDefault(m => m.SatCode == cg.Adja2Code.Value)?.Sat;
  331. if (string.IsNullOrWhiteSpace(ad2sat)) ad2sat = cg.Adja2Code.Value.ToString();
  332. var tsDtoLine = DrawDtoLineHelper.DtoLine2XStart(twoStartOption);
  333. mapControl1.DrawDtoLine($"[{msat},{nsat}]带参时差线", tsDtoLine);
  334. twoStartOption.NsEph = Ns2Eph;
  335. twoStartOption.RefDto = YbMainDto - YbAdja2Dto;
  336. twoStartOption.TargetDto = DtoSx1;
  337. var tsDtoLine1 = DrawDtoLineHelper.DtoLine2XStart(twoStartOption);
  338. mapControl1.DrawDtoLine($"[{msat},{ad2sat}]带参时差线", tsDtoLine1);
  339. }
  340. catch (Exception ex)
  341. {
  342. Serilog.Log.Error(ex, $"绘制{info.PosResType.GetEnumDisplayName()}时差线失败.PosID={info.ID},SigTime={info.SigTime}");
  343. DxHelper.MsgBoxHelper.ShowWarning($"绘制{info.PosResType.GetEnumDisplayName()}时差线失败,{ex.Message}");
  344. }
  345. }
  346. private void btnDfoCalc_Click(object sender, EventArgs e)
  347. {
  348. if (!CheckPosParam()) { return; }
  349. try
  350. {
  351. var MsAnt = txtsatStation.GetLonLat();
  352. var RefGeod = txtRefLocation.GetLonLat();
  353. var sigTime = this.sigTime.DateTime;
  354. var DfoSx = Convert.ToDouble(this.txtDfo1.Text);
  355. var DfoSx1 = Convert.ToDouble(this.txtDfo2.Text);
  356. var YbMainDfo = Convert.ToDouble(this.txtYbDfo1.Text);
  357. var YbAdja1Dfo = Convert.ToDouble(this.txtYbDfo2.Text);
  358. var YbAdja2Dfo = Convert.ToDouble(this.txtYbDfo3.Text);
  359. var tFreq = Convert.ToDouble(this.txtTargetFreq.Text) * 1e6;
  360. var tdFreq = Convert.ToDouble(this.txtTargetDFreq.Text) * 1e6;
  361. var rFreq = Convert.ToDouble(this.txtRefFreq.Text) * 1e6;
  362. var rdFreq = Convert.ToDouble(this.txtRefDFreq.Text) * 1e6;
  363. var msEph = ucEphXYZMain.EphParam();
  364. var ns1Eph = ucEphXYZAdja1.EphParam();
  365. var ns2Eph = ucEphXYZAdja2.EphParam();
  366. DfoLineTwoStartOption twoStartOption = new DfoLineTwoStartOption();
  367. twoStartOption.Freq = tFreq;
  368. twoStartOption.RefFreq = rFreq;
  369. twoStartOption.MsEph = msEph;
  370. twoStartOption.NsEph = ns1Eph;
  371. twoStartOption.MsAnt = MsAnt;
  372. twoStartOption.NsAnt = MsAnt;
  373. twoStartOption.RefGeod = RefGeod;
  374. twoStartOption.TargetDfo = DfoSx;
  375. twoStartOption.Turn1 = Math.Abs(twoStartOption.Freq - tdFreq);
  376. twoStartOption.Turn2 = Math.Abs(twoStartOption.RefFreq - rdFreq);
  377. twoStartOption.RefDfo = YbMainDfo - YbAdja1Dfo;
  378. var tsDtoLine = DrawDfoLineHelper.DfoLineTwoStart(twoStartOption);
  379. var msat = listSat?.FirstOrDefault(m => m.SatCode == cg.MainCode.Value)?.Sat; ;
  380. var ad1sat = listSat?.FirstOrDefault(m => m.SatCode == cg.Adja1Code.Value)?.Sat;
  381. if (msat == null) msat = cg.MainCode.ToString();
  382. if (ad1sat == null) ad1sat = cg.Adja1Code.ToString();
  383. mapControl1.DrawDtoLine($"{info.PosResType.GetEnumDisplayName()}[{msat},{ad1sat}]频差线", tsDtoLine);
  384. twoStartOption.NsEph = ns2Eph;
  385. twoStartOption.TargetDfo = DfoSx1;
  386. twoStartOption.RefDfo = YbMainDfo - YbAdja2Dfo;
  387. var tsDtoLine1 = DrawDfoLineHelper.DfoLineTwoStart(twoStartOption);
  388. string ad2sat = string.Empty;
  389. if (cg.Adja2Code.HasValue)
  390. {
  391. ad2sat = listSat.FirstOrDefault(m => m.SatCode == cg.Adja2Code.Value)?.Sat;
  392. }
  393. mapControl1.DrawDtoLine($"{info.PosResType.GetEnumDisplayName()}[{msat},{ad2sat}]频差线", tsDtoLine1);
  394. }
  395. catch (Exception ex)
  396. {
  397. Serilog.Log.Error(ex, $"绘制{info.PosResType.GetEnumDisplayName()}频差线失败.PosID={info.ID},SigTime={sigTime}");
  398. DxHelper.MsgBoxHelper.ShowError($"绘制{info.PosResType.GetEnumDisplayName()}频差线失败,{ex.Message}");
  399. }
  400. }
  401. }
  402. }