X3DTOParamEditor.cs 13 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301
  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.Entity;
  10. using XdCxRhDW.Repostory;
  11. namespace XdCxRhDW.App.EditForms
  12. {
  13. public partial class X3DTOParamEditor : DevExpress.XtraEditors.XtraUserControl
  14. {
  15. private MapControl mapControl1;
  16. private PosRes info;
  17. private CgRes cg;
  18. private List<SatInfo> listSat;
  19. public X3DTOParamEditor(PosRes info, MapControl mapControl)
  20. {
  21. InitializeComponent();
  22. this.info = info;
  23. itemSigTime.Text = $"{itemSigTime.Text}({SysConfig.Config.TimeZoneUTC})";
  24. this.layoutControl1.UseDefault();
  25. sigTime.UseDefault();
  26. this.Text = $"{info.PosResType.GetEnumDisplayName()}时差参数";
  27. this.mapControl1 = mapControl;
  28. }
  29. private async void X3DTOParamEditor_Load(object sender, EventArgs e)
  30. {
  31. StationRes station = null;
  32. using (RHDWContext db = new RHDWContext())
  33. {
  34. listSat = await db.SatInfos.ToListAsync();
  35. }
  36. using (RHDWPartContext db = RHDWPartContext.GetContext(info.SigTime))
  37. {
  38. cg = await db?.CgRes.Where(m => m.ID == info.CgResID).FirstOrDefaultAsync();
  39. station = await db?.StationRes.Where(m => m.ID == info.StationResID).FirstOrDefaultAsync();
  40. }
  41. if (cg != null)
  42. {
  43. this.txtDtoSx.Text = $"{cg.Dto1.Value}";
  44. this.txtDtoSx1.Text = $"{cg.Dto2.Value}";
  45. if (info.PosResType == EnumPosResType.X3)
  46. {
  47. this.txtYbMain.Text = $"{cg.YbMainDto.Value}";
  48. this.txtYbAdja.Text = $"{cg.YbAdja1Dto.Value}";
  49. this.txtYbAdja1.Text = $"{cg.YbAdja2Dto.Value}";
  50. }
  51. this.sigTime.EditValue = info.SigTime;
  52. var mainEph = (cg.MainX, cg.MainY, cg.MainZ);
  53. ucEphXYZMain.SetXYZ($"主星", cg.MainCode, mainEph, Color.Black);
  54. var adja1Eph = (cg.Adja1X, cg.Adja1Y, cg.Adja1Z);
  55. ucEphXYZAdja1.SetXYZ($"邻星1", cg.Adja1Code, adja1Eph, Color.Black);
  56. var adja2Eph = (cg.Adja2X, cg.Adja2Y, cg.Adja2Z);
  57. ucEphXYZAdja2.SetXYZ($"邻星2", cg.Adja2Code, adja2Eph, Color.Black);
  58. }
  59. if (station != null)
  60. {
  61. this.txtsatStation.Text = $"{station.SatTxLon},{station.SatTxLat}";
  62. this.txtRefLocation.Text = $"{station.RefLon},{station.RefLat}";
  63. }
  64. if (info.PosResType == EnumPosResType.X3NoRef)
  65. {
  66. txtYbMain.Properties.ReadOnly = true;
  67. txtYbAdja.Properties.ReadOnly = true;
  68. txtYbAdja1.Properties.ReadOnly = true;
  69. txtRefLocation.Properties.ReadOnly = true;
  70. this.txtRefLocation.Text = $"{0},{0}";
  71. }
  72. }
  73. public bool CheckParam()
  74. {
  75. dxErrorProvider.ClearErrors();
  76. if (!txtsatStation.CheckLonLat(dxErrorProvider, "接收站"))
  77. {
  78. return false;
  79. }
  80. if (info.PosResType == EnumPosResType.X3 && !txtRefLocation.CheckLonLat(dxErrorProvider, "参考站"))
  81. {
  82. return false;
  83. }
  84. if (!txtDtoSx.CheckDouble(dxErrorProvider, "主邻1时差"))
  85. {
  86. return false;
  87. }
  88. if (!txtDtoSx1.CheckDouble(dxErrorProvider, "主邻2时差"))
  89. {
  90. return false;
  91. }
  92. if (info.PosResType == EnumPosResType.X3 && !txtYbMain.CheckDouble(dxErrorProvider, "样本主星时差"))
  93. {
  94. return false;
  95. }
  96. if (info.PosResType == EnumPosResType.X3 && !txtYbAdja.CheckDouble(dxErrorProvider, "样本邻星1时差"))
  97. {
  98. return false;
  99. }
  100. if (info.PosResType == EnumPosResType.X3 && !txtYbAdja1.CheckDouble(dxErrorProvider, "样本邻星2时差"))
  101. {
  102. return false;
  103. }
  104. if (!ucEphXYZMain.CheckEphXYZ(dxErrorProvider))
  105. {
  106. return false;
  107. }
  108. if (!ucEphXYZAdja1.CheckEphXYZ(dxErrorProvider))
  109. {
  110. return false;
  111. }
  112. if (!ucEphXYZAdja2.CheckEphXYZ(dxErrorProvider))
  113. {
  114. return false;
  115. }
  116. return true;
  117. }
  118. private void btnDtoLine_Click(object sender, EventArgs e)
  119. {
  120. if (!CheckParam()) { return; }
  121. try
  122. {
  123. var MsAnt = txtsatStation.GetLonLat();
  124. var RefGeod = info.PosResType == EnumPosResType.X3 ? txtRefLocation.GetLonLat() : new double[3] { 0, 0, 0 };
  125. var DtoSx = Convert.ToDouble(this.txtDtoSx.Text);
  126. var DtoSx1 = Convert.ToDouble(this.txtDtoSx1.Text);
  127. var YbMainDto = info.PosResType == EnumPosResType.X3 ? Convert.ToDouble(this.txtYbMain.Text) : 0;
  128. var YbAdja1Dto = info.PosResType == EnumPosResType.X3 ? Convert.ToDouble(this.txtYbAdja.Text) : 0;
  129. var YbAdja2Dto = info.PosResType == EnumPosResType.X3 ? Convert.ToDouble(this.txtYbAdja1.Text) : 0;
  130. double[] msEph = ucEphXYZMain.EphXYZ();
  131. double[] Ns1Eph = ucEphXYZAdja1.EphXYZ();
  132. double[] Ns2Eph = ucEphXYZAdja2.EphXYZ();
  133. DtoLineTwoStartOption twoStartOption = new DtoLineTwoStartOption();
  134. twoStartOption.MsEph = msEph;
  135. twoStartOption.NsEph = Ns1Eph;
  136. twoStartOption.MsAnt = MsAnt;
  137. twoStartOption.NsAnt = MsAnt;
  138. twoStartOption.RefGeod = RefGeod;
  139. twoStartOption.TargetDto = DtoSx;
  140. twoStartOption.RefDto = YbMainDto - YbAdja1Dto;
  141. twoStartOption.PosLon = info.PosLon;
  142. twoStartOption.PosLat = info.PosLat;
  143. var msat = listSat.FirstOrDefault(m => m.SatCode == cg.MainCode.Value)?.Sat;
  144. if (string.IsNullOrWhiteSpace(msat)) msat = cg.MainCode.Value.ToString();
  145. var nsat = listSat.FirstOrDefault(m => m.SatCode == cg.Adja1Code.Value)?.Sat;
  146. if (string.IsNullOrWhiteSpace(nsat)) nsat = cg.Adja1Code.Value.ToString();
  147. var ad2sat = listSat.FirstOrDefault(m => m.SatCode == cg.Adja2Code.Value)?.Sat;
  148. if (string.IsNullOrWhiteSpace(ad2sat)) ad2sat = cg.Adja2Code.Value.ToString();
  149. if (info.PosResType == EnumPosResType.X3)
  150. {
  151. var tsDtoLine = DrawDtoLineHelper.DtoLine2XStart(twoStartOption);
  152. mapControl1.DrawDtoLine($"[{msat},{nsat}]带参时差线", tsDtoLine);
  153. twoStartOption.NsEph = Ns2Eph;
  154. twoStartOption.RefDto = YbMainDto - YbAdja2Dto;
  155. twoStartOption.TargetDto = DtoSx1;
  156. var tsDtoLine1 = DrawDtoLineHelper.DtoLine2XStart(twoStartOption);
  157. mapControl1.DrawDtoLine($"[{msat},{ad2sat}]带参时差线", tsDtoLine1);
  158. }
  159. else
  160. {
  161. var tsDtoLine = DrawDtoLineHelper.DtoLine2XNoRefStart(twoStartOption);
  162. mapControl1.DrawDtoLine($"[{msat},{nsat}]无参时差线", tsDtoLine);
  163. twoStartOption.NsEph = Ns2Eph;
  164. twoStartOption.TargetDto = DtoSx1;
  165. var tsDtoLine1 = DrawDtoLineHelper.DtoLine2XNoRefStart(twoStartOption);
  166. mapControl1.DrawDtoLine($"[{msat},{ad2sat}]无参时差线", tsDtoLine1);
  167. }
  168. }
  169. catch (Exception ex)
  170. {
  171. Serilog.Log.Error(ex, $"绘制{info.PosResType.GetEnumDisplayName()}时差线失败.PosID={info.ID},SigTime={info.SigTime}");
  172. DxHelper.MsgBoxHelper.ShowWarning($"绘制{info.PosResType.GetEnumDisplayName()}时差线失败,{ex.Message}");
  173. }
  174. }
  175. private async void btnEphCalc_Click(object sender, EventArgs e)
  176. {
  177. dxErrorProvider.ClearErrors();
  178. if (this.sigTime.DateTime == DateTime.MinValue)
  179. {
  180. dxErrorProvider.SetError(this.sigTime, "信号时间不能为空!");
  181. return;
  182. }
  183. var sigTime = this.sigTime.DateTime;
  184. try
  185. {
  186. var xlInfo = await XlRepository.GetLatestAsync(cg.MainCode.Value, sigTime);
  187. if (xlInfo == null)
  188. {
  189. DxHelper.MsgBoxHelper.ShowWarning($"卫星:{cg.MainCode.Value}未找到对应的星历信息,请导入星历");
  190. return;
  191. }
  192. var xlInfo1 = await XlRepository.GetLatestAsync(cg.Adja1Code.Value, sigTime);
  193. if (xlInfo1 == null)
  194. {
  195. DxHelper.MsgBoxHelper.ShowWarning($"卫星:{cg.Adja1Code.Value}未找到对应的星历信息,请导入星历");
  196. return;
  197. }
  198. var xlInfo2 = await XlRepository.GetLatestAsync(cg.Adja2Code.Value, sigTime);
  199. if (xlInfo2 == null)
  200. {
  201. DxHelper.MsgBoxHelper.ShowWarning($"卫星:{cg.Adja2Code.Value}未找到对应的星历信息,请导入星历");
  202. return;
  203. }
  204. var mEph = EphHelper.Calc(xlInfo.TwoLine, sigTime.ToUtc());
  205. var nEph1 = EphHelper.Calc(xlInfo1.TwoLine, sigTime.ToUtc());
  206. var nEph2 = EphHelper.Calc(xlInfo2.TwoLine, sigTime.ToUtc());
  207. var mainEph = (mEph.X, mEph.Y, mEph.Z);
  208. ucEphXYZMain.SetXYZ($"主星", cg.MainCode, mainEph, Color.Red);
  209. var adja1Eph = (nEph1.X, nEph1.Y, nEph1.Z);
  210. ucEphXYZAdja1.SetXYZ($"邻星1", cg.Adja1Code, adja1Eph, Color.Red);
  211. var adja2Eph = (nEph2.X, nEph2.Y, nEph2.Z);
  212. ucEphXYZAdja2.SetXYZ($"邻星2", cg.Adja2Code, adja2Eph, Color.Red);
  213. }
  214. catch (Exception ex)
  215. {
  216. Serilog.Log.Error(ex, $"手动{info.PosResType.GetEnumDisplayName()}推算星历失败.PosID={info.ID},SigTime={info.SigTime}");
  217. DxHelper.MsgBoxHelper.ShowWarning($"手动{info.PosResType.GetEnumDisplayName()}推算星历失败,{ex.Message}");
  218. }
  219. }
  220. private void btnPos_Click(object sender, EventArgs e)
  221. {
  222. if (!CheckParam()) { return; }
  223. try
  224. {
  225. var MsAnt = txtsatStation.GetLonLat();
  226. var RefGeod = info.PosResType == EnumPosResType.X3 ? txtRefLocation.GetLonLat() : new double[3] { 0, 0, 0 };
  227. var sigTime = this.sigTime.DateTime;
  228. var DtoSx = Convert.ToDouble(this.txtDtoSx.Text);
  229. var DtoSx1 = Convert.ToDouble(this.txtDtoSx1.Text);
  230. var YbMainDto = info.PosResType == EnumPosResType.X3 ? Convert.ToDouble(this.txtYbMain.Text) : 0;
  231. var YbAdja1Dto = info.PosResType == EnumPosResType.X3 ? Convert.ToDouble(this.txtYbAdja.Text) : 0;
  232. var YbAdja2Dto = info.PosResType == EnumPosResType.X3 ? Convert.ToDouble(this.txtYbAdja1.Text) : 0;
  233. double[] msEph = ucEphXYZMain.EphXYZ();
  234. double[] Ns1Eph = ucEphXYZAdja1.EphXYZ();
  235. double[] Ns2Eph = ucEphXYZAdja2.EphXYZ();
  236. var StationRes = new StationRes()
  237. {
  238. SatTxLon = MsAnt[0],
  239. SatTxLat = MsAnt[1],
  240. RefLon = RefGeod[0],
  241. RefLat = RefGeod[1],
  242. };
  243. var cgRes = new CgRes()
  244. {
  245. SigTime = sigTime,
  246. Dto1 = DtoSx,
  247. Dto2 = DtoSx1,
  248. YbMainDto = YbMainDto,
  249. YbAdja1Dto = YbAdja1Dto,
  250. YbAdja2Dto = YbAdja2Dto,
  251. MainX = msEph[0],
  252. MainY = msEph[1],
  253. MainZ = msEph[2],
  254. Adja1X = Ns1Eph[0],
  255. Adja1Y = Ns1Eph[1],
  256. Adja1Z = Ns1Eph[2],
  257. Adja2X = Ns2Eph[0],
  258. Adja2Y = Ns2Eph[1],
  259. Adja2Z = Ns2Eph[2],
  260. };
  261. var res = info.PosResType == EnumPosResType.X3 ? PosApi.X3_Pos(cgRes, StationRes) : PosApi.X3_PosNoRef(cgRes, StationRes);
  262. this.txtPosRes.Text = $"{info.PosResType.GetEnumDisplayName()}定位点:[{res[0]:f4},{res[1]:f4}] 镜像点:[{res[3]:f4},{res[4]:f4}]";
  263. }
  264. catch (Exception ex)
  265. {
  266. Serilog.Log.Error(ex, $"{info.PosResType.GetEnumDisplayName()}手动定位失败.PosID={info.ID},SigTime={info.SigTime}");
  267. DxHelper.MsgBoxHelper.ShowWarning($"{info.PosResType.GetEnumDisplayName()}手动定位失败,{ex.Message}");
  268. }
  269. }
  270. }
  271. }