XdCxRhDW/XdCxRhDW.App/UserControl/X3GDOPParam.cs

using DevExpress.XtraMap;
using ExtensionsDev;
using System;
using System.Data;
using System.Linq;
using System.Collections.Generic;
using DxHelper;
using XdCxRhDW.Repostory;
using XdCxRhDW.Entity;
using XdCxRhDW.Api;
using System.Drawing;
using DevExpress.XtraEditors.DXErrorProvider;
using XdCxRhDW.Dto;
using System.Configuration;
namespace XdCxRhDW.App.UserControl
{
    public partial class X3GDOPParam : DevExpress.XtraEditors.XtraUserControl
    {
        public MapControl mapControl1;
        public X3GDOP星地接口 Model => new X3GDOP星地接口()
        {
            EphMain = ucEphXYZMain.EphParam(),
            EphAdja1 = ucEphXYZAdja1.EphParam(),
            EphAdja2 = ucEphXYZAdja2.EphParam(),
            CapTime = txtCapTime.DateTime,
            RefStation = PosResType == EnumPosResType.X3NoRef ? new double[3] { 0, 0, 0 } : txtRefLocation1.GetLonLat(),
            DtousErr = Convert.ToDouble(txtDtousErr1.Text),
            SatLocErr = Convert.ToDouble(txtSatLocErr1.Text),
        };
        private EnumPosResType PosResType;
        public X3GDOPParam(PosRes item)
        {
            InitializeComponent();
            this.layoutControl1.UseDefault();
            txtCapTime.UseDefault();
            txtRefLocation1.UseDoubleClickToSelectAll();
            this.txtCapTime.DateTime = item.SigTime;
            this.txtDtousErr1.EditValue = 1;
            this.txtSatLocErr1.EditValue = 10000;
            PosResType = item.PosResType;
            if (PosResType == EnumPosResType.X3NoRef)
            {
                layoutControlItem14.Visibility = DevExpress.XtraLayout.Utils.LayoutVisibility.Never;
            }

            using (RHDWPartContext db = RHDWPartContext.GetContext(item.SigTime))
            {
                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)
                {
                    /*
                    cg.MainX = 14134595.9416663;
                    cg.MainY = 39728560.9832403;
                    cg.MainZ = -21797.4303592971;

                    cg.Adja1X = 17189204.6561764;
                    cg.Adja1Y = 38508305.0902672;
                    cg.Adja1Z = 12338.5166550109;

                    cg.Adja2X = 25377964.3839722;
                    cg.Adja2Y = 33675447.9787357;
                    cg.Adja2Z = 3106.70347014127;
                    */
                    var mainEph = (cg.MainX, cg.MainY, cg.MainZ, cg.MainVx, cg.MainVy, cg.MainVz);
                    ucEphXYZMain.SetParam($"主星", cg.MainCode, mainEph, Color.Black);

                    var adja1Eph = (cg.Adja1X, cg.Adja1Y, cg.Adja1Z, cg.Adja1Vx, cg.Adja1Vy, cg.Adja1Vz);
                    ucEphXYZAdja1.SetParam($"邻星1", cg.Adja1Code, adja1Eph, Color.Black);

                    var adja2Eph = (cg.Adja2X, cg.Adja2Y, cg.Adja2Z, cg.Adja2Vx, cg.Adja2Vy, cg.Adja2Vz);
                    ucEphXYZAdja2.SetParam($"邻星2", cg.Adja2Code, adja2Eph, Color.Black);
                }
            }
        }

        private void btnClose_Click(object sender, EventArgs e)
        {
            DxHelper.PopupHelper.HidePopup(this);
        }
        private bool ParamValidate()
        {
            if (!ucEphXYZMain.CheckEph(dxErrorProvider))
            {
                return false;
            }
            if (!ucEphXYZAdja1.CheckEph(dxErrorProvider))
            {
                return false;
            }
            if (!ucEphXYZAdja2.CheckEph(dxErrorProvider))
            {
                return false;
            }
            if (PosResType == EnumPosResType.X3 && !txtRefLocation1.CheckLonLat(dxErrorProvider, "参考站"))
            {
                return false;
            }
            if (!txtDtousErr1.CheckDouble(dxErrorProvider, "时差误差"))
            {
                return false;
            }
            if (!txtSatLocErr1.CheckDouble(dxErrorProvider, "星历位置误差"))
            {
                return false;
            }
            return true;
        }
        private void btnX1D2_Click(object sender, EventArgs e)
        {
            if (!ParamValidate()) { return; }
            try
            {
                mapControl1.ClearMap();
                string useNewGDOP = ConfigurationManager.AppSettings["UseNewGDOP"];
                if (useNewGDOP == null || useNewGDOP == "0")
                {
                    var data = GdopHelper.Gdop3SatByXyz(Model.EphMain, Model.EphAdja1, Model.EphAdja2
                                  , Model.DtousErr, Model.SatLocErr, PosResType == EnumPosResType.X3NoRef ? null : Model.RefStation);
                    foreach (var errLins in data)//画GDOP
                    {
                        mapControl1.DrawGdopLineTwo(errLins.ErrDistanceKm, errLins.MapDots, 1);
                    }

                }
                else
                {
                    var data = GdopHelper.Gdop3SatByXyzNew(Model.EphMain, Model.EphAdja1, Model.EphAdja2
                                     , Model.DtousErr, Model.SatLocErr, PosResType == EnumPosResType.X3NoRef ? null : Model.RefStation);

                    List<MapPolyline> polylines = new List<MapPolyline>();
                    for (int i = 0; i < data.Count; i++)
                    {
                        var errLins = data[i];
                        var mapline = mapControl1.GetLine(errLins.ErrDistanceKm, errLins.MapDots);
                        if (mapline == null)
                        {
                            continue;
                        }
                        polylines.Add(mapline);
                    }
                    mapControl1.DrawDtoLine(polylines);
                }
            }
            catch (Exception ex)
            {
                DxHelper.MsgBoxHelper.ShowError($"绘制{PosResType.GetEnumDisplayName()}GDOP失败,{ex.Message}");
            }
        }

        private async void btnEphCalc_Click(object sender, EventArgs e)
        {
            dxErrorProvider.ClearErrors();
            if (txtCapTime.DateTime == DateTime.MinValue)
            {
                dxErrorProvider.SetError(txtCapTime, "信号时间不能为空!");
                return;
            }
            var sigTime = txtCapTime.DateTime;
            try
            {
                var mainCode = ucEphXYZMain.GetSatCode();
                var adja1Code = ucEphXYZAdja1.GetSatCode();
                var adja2Code = ucEphXYZAdja2.GetSatCode();
                var mainxlInfo = await XlRepository.GetLatestAsync(mainCode, sigTime);
                if (mainxlInfo == null)
                {
                    DxHelper.MsgBoxHelper.ShowWarning($"卫星:{mainCode}未找到对应的星历信息,请导入星历");
                    return;
                }
                var XlCalcDto = new XlCalcDto() { tleStr = mainxlInfo.TwoLine, SigTime = sigTime };
                var maineph = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
                ucEphXYZMain.SetParam("主星", mainCode, (maineph.data.X, maineph.data.Y, maineph.data.Z, maineph.data.VX, maineph.data.VY, maineph.data.VZ), Color.Red);

                var adja1xlInfo = await XlRepository.GetLatestAsync(adja1Code, sigTime);
                if (adja1xlInfo == null)
                {
                    DxHelper.MsgBoxHelper.ShowWarning($"卫星:{adja1Code}未找到对应的星历信息,请导入星历");
                    return;
                }
                XlCalcDto = new XlCalcDto() { tleStr = adja1xlInfo.TwoLine, SigTime = sigTime };
                var adja1eph = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
                ucEphXYZAdja1.SetParam("邻星1", adja1Code, (adja1eph.data.X, adja1eph.data.Y, adja1eph.data.Z, adja1eph.data.VX, adja1eph.data.VY, adja1eph.data.VZ), Color.Red);

                var adja2xlInfo = await XlRepository.GetLatestAsync(adja2Code, sigTime);
                if (adja2xlInfo == null)
                {
                    DxHelper.MsgBoxHelper.ShowWarning($"卫星:{adja2Code}未找到对应的星历信息,请导入星历");
                    return;
                }
                XlCalcDto = new XlCalcDto() { tleStr = adja2xlInfo.TwoLine, SigTime = sigTime };
                var adja2eph = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
                ucEphXYZAdja2.SetParam("邻星2", adja2Code, (adja2eph.data.X, adja2eph.data.Y, adja2eph.data.Z, adja2eph.data.VX, adja2eph.data.VY, adja2eph.data.VZ), Color.Red);
            }
            catch (Exception ex)
            {
                Serilog.Log.Error(ex, $"手动推算{PosResType.GetEnumDisplayName()}星历失败,SigTime={sigTime}");
                DxHelper.MsgBoxHelper.ShowError($"手动推算{PosResType.GetEnumDisplayName()}星历失败,{ex.Message}");
            }
        }
    }
    public class X3GDOP星地接口
    {
        /// <summary>
        /// 主星星历x y z vx vy vz
        /// </summary>
        public double[] EphMain { get; set; }

        /// <summary>
        /// 邻星1星历x y z vx vy vz
        /// </summary>
        public double[] EphAdja1 { get; set; }

        /// <summary>
        /// 邻星2星历x y z vx vy vz
        /// </summary>
        public double[] EphAdja2 { get; set; }

        /// <summary>
        /// 采集时刻
        /// </summary>
        public DateTime CapTime { get; set; }

        /// <summary>
        /// 参考站位置-经度 纬度 高度
        /// </summary>
        public double[] RefStation { get; set; }
        /// <summary>
        /// 时差误差(单位us)
        /// </summary>
        public double DtousErr { get; set; } = 1;

        /// <summary>
        /// 星历位置误差(单位米)
        /// </summary>
        public double SatLocErr { get; set; } = 10000;

    }
}