123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441 |
- using DevExpress.XtraMap;
- using DxHelper;
- using ExtensionsDev;
- using System;
- using System.CodeDom;
- using System.Collections.Generic;
- using System.ComponentModel;
- using System.Configuration;
- using System.Data;
- using System.Data.Entity;
- using System.Drawing;
- using System.Linq;
- using System.Text;
- using System.Threading.Tasks;
- using XdCxRhDW.Api;
- using XdCxRhDW.Dto;
- using XdCxRhDW.Entity;
- using XdCxRhDW.Framework;
- using XdCxRhDW.Repostory;
- namespace XdCxRhDW.App.EditForms
- {
- public partial class X2D1DTOParamEditor : DevExpress.XtraEditors.XtraUserControl
- {
- private MapControl mapControl1;
- private PosRes info;
- private CgRes cg;
- private List<SatInfo> listSat;
- public X2D1DTOParamEditor(PosRes info, MapControl mapControl)
- {
- InitializeComponent();
- this.info = info;
- itemSigTime.Text = $"{itemSigTime.Text}({SysConfig.Config.TimeZoneUTC})";
- this.layoutControl1.UseDefault();
- txtSigTime.UseDefault();
- this.Text = $"{info.PosResType.GetEnumDisplayName()}时差参数";
- this.listSat = new List<SatInfo>();
- this.mapControl1 = mapControl;
- }
- private async void X2D1DTOParamEditor_Load(object sender, EventArgs e)
- {
- this.layoutControl1.BestFit();
- StationRes station = null;
- using (RHDWContext db = new RHDWContext())
- {
- listSat = await db.SatInfos.ToListAsync();
- }
- using (RHDWPartContext db = RHDWPartContext.GetContext(info.SigTime))
- {
- cg = await db?.CgRes.Where(m => m.ID == info.CgResID).FirstOrDefaultAsync();
- station = await db?.StationRes.Where(m => m.ID == info.StationResID).FirstOrDefaultAsync();
- }
- if (cg != null)
- {
- //if (info.PosResType == EnumPosResType.X2D1NoRef)
- //{
- // cg.Dto1 = 6348.6030;
- // cg.DtoCdb = 246616.0330;
- // cg.MainX = 11265194.9906;
- // cg.MainY = 40486056.6399;
- // cg.MainZ = -2055810.5412;
- // cg.Adja1X = 25525328.3455;
- // cg.Adja1Y = 33539987.1755;
- // cg.Adja1Z = -49504.0378;
- //}
- this.txtDtoSx.Text = $"{cg.Dto1.Value}";
- this.txtDtoCdb.Text = $"{cg.DtoCdb.Value}";
- if (info.PosResType == EnumPosResType.X2D1)
- {
- this.txtYbMain.Text = $"{cg.YbMainDto.Value}";
- this.txtYbAdja.Text = $"{cg.YbAdja1Dto.Value}";
- }
- this.txtSigTime.EditValue = info.SigTime;
- ucEphXYZMain.SetXYZ("主星", cg.MainCode, (cg.MainX, cg.MainY, cg.MainZ), Color.Black);
- ucEphXYZAdaj.SetXYZ("邻星", cg.Adja1Code, (cg.Adja1X, cg.Adja1Y, cg.Adja1Z), Color.Black);
- }
- if (station != null)
- {
- //if (info.PosResType == EnumPosResType.X2D1NoRef)
- //{
- // station.SatTxLon = 45;
- // station.SatTxLat =20;
- // station.CdbTxLon = 45;
- // station.CdbTxLat = 20;
- //}
- this.txtsatStation.Text = $"{station.SatTxLon},{station.SatTxLat}";
- this.txtcdbStation.Text = $"{station.CdbTxLon},{station.CdbTxLat}";
- this.txtRefLocation.Text = $"{station.RefLon},{station.RefLat}";
- }
- if (info.PosResType == EnumPosResType.X2D1NoRef)
- {
- txtYbMain.Properties.ReadOnly = true;
- txtYbAdja.Properties.ReadOnly = true;
- txtRefLocation.Properties.ReadOnly = true;
- this.txtRefLocation.Text = $"{0},{0}";
- }
- //double tarLon = 122.88, tarLat = 31.29;
- //double recLon = 118.5, recLat = 32.1;
- //double refLon = 121, refLat = 30;
- //double cdbLon = 122, cdbLat = 30;
- //var refEcef = PhysicsHelper.GeoToEcef((refLon, refLat, 0));
- //var tarEcef = PhysicsHelper.GeoToEcef((tarLon, tarLat, 0));
- //var recEcef = PhysicsHelper.GeoToEcef((recLon, recLat, 0));
- //var cdbEcef = PhysicsHelper.GeoToEcef((cdbLon, cdbLat, 0));
- //var xl1 = (-41338971.25, 6141819.55, 5504123.22);
- //var dt1 = PhysicsHelper.Dto(tarEcef, xl1, recEcef);
- //var dt2 = PhysicsHelper.Dto(tarEcef, cdbEcef);
- //var dto1 = dt1 - dt2;
- //var dt3 = PhysicsHelper.Dto(refEcef, xl1, recEcef);
- //var dto2 = dt3 - 0;
- }
- public bool CheckParam()
- {
- dxErrorProvider.ClearErrors();
- if (!txtsatStation.CheckLonLat(dxErrorProvider, "接收站"))
- {
- return false;
- }
- if (!txtcdbStation.CheckLonLat(dxErrorProvider, "超短波"))
- {
- return false;
- }
- if (!string.IsNullOrWhiteSpace(txtRefLocation.Text) && !txtRefLocation.CheckLonLat(dxErrorProvider, "参考站"))
- {
- return false;
- }
- if (!txtDtoSx.CheckDouble(dxErrorProvider, "双星时差"))
- {
- return false;
- }
- if (!txtDtoCdb.CheckDouble(dxErrorProvider, "主星超短时差"))
- {
- return false;
- }
- if (!string.IsNullOrWhiteSpace(txtYbMain.Text) && !txtYbMain.CheckDouble(dxErrorProvider, "样本主星时差"))
- {
- return false;
- }
- if (!string.IsNullOrWhiteSpace(txtYbAdja.Text) && !txtYbAdja.CheckDouble(dxErrorProvider, "样本邻星时差"))
- {
- return false;
- }
- if (!ucEphXYZMain.CheckEphXYZ(dxErrorProvider))
- {
- return false;
- }
- if (!ucEphXYZAdaj.CheckEphXYZ(dxErrorProvider))
- {
- return false;
- }
- return true;
- }
- private void btnDtoLine_Click(object sender, EventArgs e)
- {
- if (!CheckParam()) { return; }
- bool isRef = false;
- if (!string.IsNullOrWhiteSpace(txtRefLocation.Text)
- && !string.IsNullOrWhiteSpace(txtYbMain.Text)
- && !string.IsNullOrWhiteSpace(txtYbAdja.Text))
- {
- isRef = true;
- }
- try
- {
- var MsAnt = txtsatStation.GetLonLat();
- var CDBAnt = txtcdbStation.GetLonLat();
- var RefGeod = isRef ? txtRefLocation.GetLonLat() : new double[3] { 0, 0, 0 };
- var DtoSx = Convert.ToDouble(this.txtDtoSx.Text);
- var DtoCdb = Convert.ToDouble(this.txtDtoCdb.Text);
- var YbMainDto = isRef ? Convert.ToDouble(this.txtYbMain.Text) : 0;
- var YbAdja1Dto = isRef ? Convert.ToDouble(this.txtYbAdja.Text) : 0;
- double[] msEph = ucEphXYZMain.EphXYZ();
- double[] NsEph = ucEphXYZAdaj.EphXYZ();
- DtoLineXdOption dtoLineXd = new DtoLineXdOption();
- dtoLineXd.MsEph = msEph;
- dtoLineXd.MsAnt = MsAnt;
- dtoLineXd.CDBAnt = CDBAnt;
- dtoLineXd.RefGeod = RefGeod;
- dtoLineXd.xdDto = DtoCdb;
- dtoLineXd.RefDto = YbMainDto;
- dtoLineXd.PosLon = info.PosLon;
- dtoLineXd.PosLat = info.PosLat;
- DtoLineTwoStartOption twoStartOption = new DtoLineTwoStartOption();
- twoStartOption.MsEph = msEph;
- twoStartOption.NsEph = NsEph;
- twoStartOption.MsAnt = MsAnt;
- twoStartOption.NsAnt = MsAnt;
- twoStartOption.RefGeod = RefGeod;
- twoStartOption.TargetDto = DtoSx;
- twoStartOption.RefDto = YbMainDto - YbAdja1Dto;
- twoStartOption.PosLon = info.PosLon;
- twoStartOption.PosLat = info.PosLat;
- var msat = listSat.FirstOrDefault(m => m.SatCode == cg.MainCode.Value)?.Sat;
- if (string.IsNullOrWhiteSpace(msat)) msat = cg.MainCode.Value.ToString();
- var nsat = listSat.FirstOrDefault(m => m.SatCode == cg.Adja1Code.Value)?.Sat;
- if (string.IsNullOrWhiteSpace(nsat)) nsat = cg.Adja1Code.Value.ToString();
- if (isRef)
- {
- var xdDtoLine = DrawDtoLineHelper.DtoLineXdNew(dtoLineXd);
- List<MapPolyline> polylines = new List<MapPolyline>();
- foreach (var dtoLine in xdDtoLine)
- {
- var mapline = mapControl1.GetLine($"[{msat},超短{CDBAnt[0]}°]带参时差线", dtoLine.dtoLinePoints, true);
- if (mapline == null)
- {
- continue;
- }
- polylines.Add(mapline);
- }
- mapControl1.DrawDtoLine(polylines);
- var tsDtoLine = DrawDtoLineHelper.DtoLine2XStart(twoStartOption);
- mapControl1.DrawDtoLine($"[{msat},{nsat}]带参时差线", tsDtoLine);
- }
- else
- {
- var xdDtoLine = DrawDtoLineHelper.DtoLineXdNoRefNew(dtoLineXd);
- List<MapPolyline> polylines = new List<MapPolyline>();
- foreach (var dtoLine in xdDtoLine)
- {
- var mapline = mapControl1.GetLine($"[{msat},超短{CDBAnt[0]}°]无参时差线", dtoLine.dtoLinePoints, true);
- if (mapline == null)
- {
- continue;
- }
- polylines.Add(mapline);
- }
- mapControl1.DrawDtoLine(polylines);
- var tsDtoLine = DrawDtoLineHelper.DtoLine2XNoRefStart(twoStartOption);
- mapControl1.DrawDtoLine($"[{msat},{nsat}]无参时差线", tsDtoLine);
- }
- }
- catch (Exception ex)
- {
- XdCxRhDW.Framework.LogHelper.Error($"绘制{info.PosResType.GetEnumDisplayName()}时差线失败.PosID={info.ID},SigTime={info.SigTime}", ex);
- DxHelper.MsgBoxHelper.ShowWarning($"绘制{info.PosResType.GetEnumDisplayName()}时差线失败,{ex.Message}");
- }
- }
- private async void btnEphCalc_Click(object sender, EventArgs e)
- {
- dxErrorProvider.ClearErrors();
- if (this.txtSigTime.DateTime == DateTime.MinValue)
- {
- dxErrorProvider.SetError(txtSigTime, "信号时间不能为空!");
- return;
- }
- var sigTime = this.txtSigTime.DateTime;
- try
- {
- var mainxlInfo = await XlRepository.GetLatestAsync(cg.MainCode.Value, sigTime);
- if (mainxlInfo == null)
- {
- DxHelper.MsgBoxHelper.ShowWarning($"卫星:{cg.MainCode.Value}未找到对应的星历信息,请导入星历");
- return;
- }
- var XlCalcDto = new XlCalcDto() { tleStr = mainxlInfo.TwoLine, SigTime = sigTime };
- var maineph = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
- ucEphXYZMain.SetXYZ("主星", cg.MainCode.Value, (maineph.data.X, maineph.data.Y, maineph.data.Z), Color.Red);
- var adjaxlInfo = await XlRepository.GetLatestAsync(cg.Adja1Code.Value, sigTime);
- if (adjaxlInfo == null)
- {
- DxHelper.MsgBoxHelper.ShowWarning($"卫星:{cg.Adja1Code.Value}未找到对应的星历信息,请导入星历");
- return;
- }
- XlCalcDto = new XlCalcDto() { tleStr = adjaxlInfo.TwoLine, SigTime = sigTime };
- var adjaeph = await HttpHelper.PostRequestAsync<SatEphResDto>(SysConfig.GetUrl("Xl/Calc"), XlCalcDto);
- ucEphXYZAdaj.SetXYZ("邻星", cg.Adja1Code.Value, (adjaeph.data.X, adjaeph.data.Y, adjaeph.data.Z), Color.Red);
- }
- catch (Exception ex)
- {
- XdCxRhDW.Framework.LogHelper.Error($"手动推算{info.PosResType.GetEnumDisplayName()}星历失败,SigTime={sigTime}", ex);
- DxHelper.MsgBoxHelper.ShowError($"手动推算{info.PosResType.GetEnumDisplayName()}星历失败,{ex.Message}");
- }
- }
- private void btnPos_Click(object sender, EventArgs e)
- {
- if (!CheckParam()) { return; }
- txtPosRes.Text = " ";
- bool isRef = false;
- if (!string.IsNullOrWhiteSpace(txtRefLocation.Text)
- && !string.IsNullOrWhiteSpace(txtYbMain.Text)
- && !string.IsNullOrWhiteSpace(txtYbAdja.Text))
- {
- isRef = true;
- }
- try
- {
- var MsAnt = txtsatStation.GetLonLat();
- var CDBAnt = txtcdbStation.GetLonLat();
- var RefGeod = isRef ? txtRefLocation.GetLonLat() : new double[3] { 0, 0, 0 };
- var DtoSx = Convert.ToDouble(this.txtDtoSx.Text);
- var DtoCdb = Convert.ToDouble(this.txtDtoCdb.Text);
- var YbMainDto = isRef ? Convert.ToDouble(this.txtYbMain.Text) : 0;
- var YbAdja1Dto = isRef ? Convert.ToDouble(this.txtYbAdja.Text) : 0;
- double[] msEph = ucEphXYZMain.EphXYZ();
- double[] nsEph = ucEphXYZAdaj.EphXYZ();
- var sigTime = txtSigTime.DateTime;
- var StationRes = new StationRes()
- {
- SatTxLon = MsAnt[0],
- SatTxLat = MsAnt[1],
- CdbTxLon = CDBAnt[0],
- CdbTxLat = CDBAnt[1],
- RefLon = RefGeod[0],
- RefLat = RefGeod[1],
- };
- if (!isRef)
- {
- X2D1NoParPosDto dto = new X2D1NoParPosDto()
- {
- SigTime = sigTime,
- MainCode = cg.MainCode.Value,
- AdjaCode = cg.Adja1Code.Value,
- SxDto = DtoSx,
- XdDto = DtoCdb,
- MainX = msEph[0],
- MainY = msEph[1],
- MainZ = msEph[2],
- AdjaX = nsEph[0],
- AdjaY = nsEph[1],
- AdjaZ = nsEph[2],
- SatTxLon = MsAnt[0],
- SatTxLat = MsAnt[1],
- CdbTxLon = CDBAnt[0],
- CdbTxLat = CDBAnt[1],
- };
- var cgRes = new CgRes()
- {
- SigTime = dto.SigTime,
- Dto1 = dto.SxDto,
- Dfo1 = dto.SxDfo,
- Snr1 = dto.SxSnr,
- DtoCdb = dto.XdDto,
- DfoCdb = dto.XdDfo,
- SnrCdb = dto.XdSnr,
- MainCode = dto.MainCode,
- Adja1Code = dto.AdjaCode,
- MainX = dto.MainX,
- MainY = dto.MainY,
- MainZ = dto.MainZ,
- Adja1X = dto.AdjaX,
- Adja1Y = dto.AdjaY,
- Adja1Z = dto.AdjaZ,
- };
- double[] res;
- //res = PosApi.X2D1_PosNoRef(cgRes, StationRes, dto.CalcConfidence);
- res = PosApi.X2D1_PosNoRef_ZL(cgRes, StationRes, dto.CalcConfidence);
- this.txtPosRes.Text = $"{info.PosResType.GetEnumDisplayName()}定位点:[{res[0]:f4},{res[1]:f4}] 镜像点:[{res[3]:f4},{res[4]:f4}]";
- }
- else
- {
- X2D1PosDto dto = new X2D1PosDto()
- {
- SigTime = sigTime,
- MainCode = cg.MainCode.Value,
- AdjaCode = cg.Adja1Code.Value,
- SxDto = DtoSx,
- XdDto = DtoCdb,
- MainYbDto = YbMainDto,
- AdjaYbDto = YbAdja1Dto,
- MainX = msEph[0],
- MainY = msEph[1],
- MainZ = msEph[2],
- AdjaX = nsEph[0],
- AdjaY = nsEph[1],
- AdjaZ = nsEph[2],
- SatTxLon = MsAnt[0],
- SatTxLat = MsAnt[1],
- CdbTxLon = CDBAnt[0],
- CdbTxLat = CDBAnt[1],
- RefLon = RefGeod[0],
- RefLat = RefGeod[1],
- };
- var cgRes = new CgRes()
- {
- SigTime = dto.SigTime,
- Dto1 = dto.SxDto,
- Dfo1 = dto.SxDfo,
- Snr1 = dto.SxSnr,
- DtoCdb = dto.XdDto,
- DfoCdb = dto.XdDfo,
- SnrCdb = dto.XdSnr,
- YbMainDto = dto.MainYbDto,
- YbAdja1Dto = dto.AdjaYbDto,
- //StationResID = StationRes.ID,
- MainCode = dto.MainCode,
- Adja1Code = dto.AdjaCode,
- // TaskID = runTask.ID,
- MainX = dto.MainX,
- MainY = dto.MainY,
- MainZ = dto.MainZ,
- Adja1X = dto.AdjaX,
- Adja1Y = dto.AdjaY,
- Adja1Z = dto.AdjaZ,
- };
- var res = PosApi.X2D1_Pos(cgRes, StationRes);
- this.txtPosRes.Text = $"{info.PosResType.GetEnumDisplayName()}定位点:[{res[0]:f4},{res[1]:f4}] 镜像点:[{res[3]:f4},{res[4]:f4}]";
- }
- }
- catch (Exception ex)
- {
- XdCxRhDW.Framework.LogHelper.Error($"{info.PosResType.GetEnumDisplayName()}手动定位失败.PosID={info.ID},SigTime={info.SigTime}", ex);
- DxHelper.MsgBoxHelper.ShowWarning($"{info.PosResType.GetEnumDisplayName()}手动定位失败,{ex.Message}");
- }
- }
- }
- }
|