X3DfoParamEditor.cs 19 KB

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