using DevExpress.XtraMap; using DxHelper; using ExtensionsDev; using System; using System.Collections.Generic; using System.Data; using System.Drawing; using System.Linq; using XdCxRhDW.Api; using XdCxRhDW.Entity; using XdCxRhDW.Repostory; namespace XdCxRhDW.App.UserControl { public partial class X2DFGDOPParam : DevExpress.XtraEditors.XtraUserControl { public MapControl mapControl1; public X2DFGDOP接口 Model => new X2DFGDOP接口() { EphMain = ucEphXYZMain.EphParam(), EphAdja = ucEphXYZAdja.EphParam(), CapTime = txtCapTime.DateTime, RefLon = Convert.ToDouble(txtRefLocation1.Text.Replace(",", ",").Split(',')[0].Trim()), RefLat = Convert.ToDouble(txtRefLocation1.Text.Replace(",", ",").Split(',')[1].Trim()), DfoErr = Convert.ToDouble(txtDfoErr1.Text), DtousErr = Convert.ToDouble(txtDtousErr1.Text), SatLocErr = Convert.ToDouble(txtSatLocErr1.Text), EphVelErr = Convert.ToDouble(txtEphVelErr1.Text), fu1 = Convert.ToDouble(txtFu1.Text) * 1e6, fu2 = Convert.ToDouble(txtFu2.Text) * 1e6, }; private EnumPosResType PosResType; public X2DFGDOPParam(PosRes item) { InitializeComponent(); this.layoutControl1.UseDefault(); txtCapTime.UseDefault(); txtRefLocation1.UseDoubleClickToSelectAll(); this.txtCapTime.DateTime = item.SigTime; this.txtSatLocErr1.EditValue = 10000; this.txtEphVelErr1.EditValue = 0.1; this.txtDfoErr1.EditValue = 0.01; this.txtDtousErr1.EditValue = 1; PosResType = item.PosResType; using (RHDWPartContext db = RHDWPartContext.GetContext(item.SigTime)) { if (db.DatabaseExist) { var cg = db.CgRes.Where(m => m.ID == item.CgResID).FirstOrDefault(); var station = db.StationRes.Where(m => m.ID == item.StationResID).FirstOrDefault(); if (station != null) { this.txtRefLocation1.Text = $"{station.RefLon},{station.RefLat}"; } if (cg != null) { this.txtFu1.EditValue = cg.TarFreqUp.HasValue ? cg.TarFreqUp.Value * 1e-6 : 950; this.txtFu2.EditValue = cg.RefFreqUp.HasValue ? cg.RefFreqUp.Value * 1e-6 : 950; var mainEph = (cg.MainX, cg.MainY, cg.MainZ, cg.MainVx, cg.MainVy, cg.MainVz); ucEphXYZMain.SetParam($"主星", cg.MainCode, mainEph, Color.Black); var adjaEph = (cg.Adja1X, cg.Adja1Y, cg.Adja1Z, cg.Adja1Vx, cg.Adja1Vy, cg.Adja1Vz); ucEphXYZAdja.SetParam($"邻星", cg.Adja1Code, adjaEph, Color.Black); } } } } private (bool, string) ParamValidate() { var mainparam = ucEphXYZMain.ValidateParam(); if (!mainparam.Item1) { return (false, mainparam.Item2); } var adjaparam = ucEphXYZAdja.ValidateParam(); if (!adjaparam.Item1) { return (false, adjaparam.Item2); } var refsta = txtRefLocation1.CheckLonLat("参考站"); if (!refsta.Item1) { return refsta; } return (true, ""); } private void btnOK_Click(object sender, EventArgs e) { var pv = ParamValidate(); if (!pv.Item1) { DxHelper.MsgBoxHelper.ShowWarning($"{pv.Item2}"); return; } try { mapControl1.ClearMap(); var refs = new double[] { Model.RefLon, Model.RefLat, 0 }; var data = GdopHelper.Gdop2SatDRefByXyz(Model.EphMain, Model.EphAdja, Model.fu1, Model.fu2, Model.DtousErr, Model.DfoErr, Model.SatLocErr, Model.EphVelErr, refs); if (data == null) { return; } foreach (var errLins in data)//画GDOP { var mapDots = errLins.MapDots.Select(p => p).Select(p => (p.Lon, p.Lat)); mapControl1.DrawGdopLineTwo(errLins.ErrDistanceKm, mapDots, mapDots.Count() / 2); } } catch (Exception ex) { DxHelper.MsgBoxHelper.ShowError($"绘制{PosResType.GetEnumDisplayName()}GDOP失败,{ex.Message}"); } } private void btnClose_Click(object sender, EventArgs e) { PopupHelper.HidePopup(this); } private async void btnEphCalc_Click(object sender, EventArgs e) { if (txtCapTime.DateTime == DateTime.MinValue) { DxHelper.MsgBoxHelper.ShowWarning("信号时间不能为空!"); return; } var sigTime = txtCapTime.DateTime; try { var mainCode = ucEphXYZMain.GetSatCode(); var adjaCode = ucEphXYZAdja.GetSatCode(); var mainxlInfo = await XlRepository.GetLatestAsync(mainCode, sigTime); if (mainxlInfo == null) { DxHelper.MsgBoxHelper.ShowWarning($"卫星:{mainCode}未找到对应的星历信息,请导入星历"); return; } var maineph = EphHelper.Calc(mainxlInfo.TwoLine, sigTime); ucEphXYZMain.SetParam("主星", mainCode, (maineph.X, maineph.Y, maineph.Z, maineph.VX, maineph.VY, maineph.VZ), Color.Red); var adjaxlInfo = await XlRepository.GetLatestAsync(adjaCode, sigTime); if (adjaxlInfo == null) { DxHelper.MsgBoxHelper.ShowWarning($"卫星:{adjaCode}未找到对应的星历信息,请导入星历"); return; } var adjaeph = EphHelper.Calc(adjaxlInfo.TwoLine, sigTime); ucEphXYZAdja.SetParam("邻星", adjaCode, (adjaeph.X, adjaeph.Y, adjaeph.Z, adjaeph.VX, adjaeph.VY, adjaeph.VZ), Color.Red); } catch (Exception ex) { Serilog.Log.Error(ex, $"手动推算{PosResType.GetEnumDisplayName()}星历失败,SigTime={sigTime}"); DxHelper.MsgBoxHelper.ShowError($"手动推算{PosResType.GetEnumDisplayName()}星历失败,{ex.Message}"); } } } public class X2DFGDOP接口 { /// /// 主星星历x y z vx vy vz /// public double[] EphMain { get; set; } /// /// 邻星1星历x y z vx vy vz /// public double[] EphAdja { get; set; } /// /// 采集时刻 /// public DateTime CapTime { get; set; } /// /// 参考站位置经度 /// public double RefLon { get; set; } /// /// 参考站位置纬度 /// public double RefLat { get; set; } /// /// 时差误差(单位us) /// public double DtousErr { get; set; } = 1; /// /// 频差误差(Hz) /// public double DfoErr { get; set; } /// /// 星历位置误差(单位米) /// public double SatLocErr { get; set; } = 10000; /// ///星历速度误差 /// public double EphVelErr { get; set; } /// /// 上行频点1(Hz) /// public double fu1 { get; set; } /// /// 上行频点2(Hz) /// public double fu2 { get; set; } } }