X3DfoParamEditor.cs 19 KB

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