Exemplo n.º 1
0
        /// <summary>
        /// 将所有产品归算到一个卫星上面。取平均。
        /// </summary>
        /// <param name="basePrn"></param>
        /// <returns></returns>
        public FcbOfUpd GetWideLaneFcb(SatelliteNumber basePrn)
        {
            Dictionary <SatelliteNumber, RmsedNumeral> result = GetFinalFcbOfBsd(basePrn);

            result[basePrn] = new RmsedNumeral(0, 0);
            var prns = result.Keys.ToList();

            prns.Sort();
            //生成产品
            FcbOfUpd fcb = new FcbOfUpd();

            fcb.BasePrn  = basePrn;
            fcb.WnMarker = FcbOfUpd.WideLaneMaker;
            fcb.Count    = 32;// result.Count;
            fcb.Epoch    = this.First.First.FirstKey.Start;
            //fcb.Prns = prns;
            int i = 0;

            foreach (var prn in prns)
            {
                i++;
                if (i > 32)
                {
                    break;
                }
                fcb.Add(prn, result[prn]);
            }

            return(fcb);
        }
Exemplo n.º 2
0
        /// <summary>
        /// 将所有卫星产品归算到一颗卫星上面
        /// </summary>
        /// <param name="basePrn"></param>
        /// <param name="maxRms"></param>
        /// <param name="epoch"></param>
        /// <returns></returns>
        public FcbOfUpd GetWideLaneFcbOfAllSat(SatelliteNumber basePrn, Time epoch, double maxRms)
        {
            //首先,所有的都计算一次
            var prns = this.GetAllPrns();
            var data = new BaseDictionary <SatelliteNumber, List <RmsedNumeral> >("数据集合", (prn) => new List <RmsedNumeral>());

            foreach (var prn in prns)
            {
                //计算以当前卫星为基准的差分数据
                Dictionary <SatelliteNumber, RmsedNumeral> fcbDic = GetWideLaneFcbDic(prn, epoch, maxRms);
                if (!fcbDic.ContainsKey(basePrn))
                {
                    continue;
                }                                               // 没有基准星,则不考虑

                var baseVal = fcbDic[basePrn];
                //其次,归算到指定的卫星
                foreach (var item in fcbDic)
                {
                    var currentPrn = item.Key;
                    var newVal     = item.Value - baseVal;
                    newVal.Value = Geo.Utils.DoubleUtil.GetRoundFraction(newVal.Value);
                    data.GetOrCreate(currentPrn).Add(newVal);
                }
            }

            //加权平均
            var      ave = Geo.Utils.DoubleUtil.GetRoundFractionAverageValues(data);
            FcbOfUpd fcb = new FcbOfUpd(basePrn, epoch, ave, true);

            return(fcb);
        }
Exemplo n.º 3
0
        private void button_saveWLToGNSSer_Click(object sender, EventArgs e)
        {
            var             service    = GetFcbDataService();
            var             timeperiod = this.timePeriodControl1.TimePeriod;
            SatelliteNumber basePrn    = (SatelliteNumber)this.bindingSource_basePrn.Current;

            var toPath = Path.Combine(Setting.TempDirectory, "WhuWLFcbToGNSSer" + Gnsser.Setting.FcbExtension);

            using (FcbOfUpdWriter writer = new FcbOfUpdWriter(toPath))
            {
                var Interval = TimeSpan.FromDays(1);
                for (var time = timeperiod.Start.Date; time <= timeperiod.End; time += Interval)
                {
                    var bsdFrac = service.GetWLFcbOfBsd(time, basePrn);
                    if (bsdFrac != null && bsdFrac.Count > 0)
                    {
                        FcbOfUpd fcbOfUpd = new FcbOfUpd(basePrn, time, true);
                        foreach (var kv in bsdFrac.KeyValues)
                        {
                            fcbOfUpd.Add(kv.Key, kv.Value);
                        }
                        if (fcbOfUpd.DataCount > 0)
                        {
                            writer.Write(fcbOfUpd);
                        }
                    }
                }
            }
            Geo.Utils.FormUtil.ShowOkAndOpenDirectory(Setting.TempDirectory);
        }
Exemplo n.º 4
0
        private void button_saveToGnsserFcb_Click(object sender, EventArgs e)
        {
            var             service    = GetFcbDataService();
            var             timeperiod = this.timePeriodControl1.TimePeriod;
            SatelliteNumber basePrn    = (SatelliteNumber)this.bindingSource_basePrn.Current;
            var             Interval   = TimeSpan.FromMinutes(15);


            var toPath = Path.Combine(Setting.TempDirectory, "WhuNLFcbToGNSSer" + Gnsser.Setting.FcbExtension);

            using (FcbOfUpdWriter writer = new FcbOfUpdWriter(toPath))
            {
                for (var time = timeperiod.Start; time <= timeperiod.End; time += Interval)
                {
                    FcbOfUpd fcbOfUpd = new FcbOfUpd(basePrn, time, false);
                    foreach (var prn in SatelliteNumber.DefaultGpsPrns)
                    {
                        var bsdFrac = service.GetNLFcbOfBsdValue(time, prn, basePrn);
                        if (bsdFrac != null)
                        {
                            fcbOfUpd.Add(prn, bsdFrac);
                        }
                    }

                    if (fcbOfUpd.DataCount > 0)
                    {
                        writer.Write(fcbOfUpd);
                    }
                }
            }
            Geo.Utils.FormUtil.ShowOkAndOpenDirectory(Setting.TempDirectory);
        }
Exemplo n.º 5
0
        private void button_convertToGnsserFcb_Click(object sender, EventArgs e)
        {
            var basePrn = baseSatSelectingControl1.SelectedPrn;

            if (basePrn.SatelliteType != SatelliteType.G)
            {
                MessageBox.Show("请选择GPS系统!目前似乎还不支持其它系统,如果支持了请尝试高版本先,若还没有请 Email To: [email protected]");
                return;
            }
            var            toPath = Path.Combine(Setting.TempDirectory, "FcbOfDcb" + Gnsser.Setting.FcbExtension);
            FcbOfUpdWriter writer = new FcbOfUpdWriter(toPath);
            var            period = this.timePeriodControl1.TimePeriod;

            var path = this.fileOpenControlOpath.FilePath;

            WideLaneBiasService = new WideLaneBiasService(path);

            for (var time = period.Start.Date; time <= period.End; time += TimeSpan.FromDays(1))
            {
                var data = WideLaneBiasService.Get(time);
                if (data == null)
                {
                    continue;
                }
                FcbOfUpd fcb = new FcbOfUpd(basePrn, data);
                writer.Write(fcb);
            }
            writer.Dispose();

            Geo.Utils.FormUtil.ShowOkAndOpenDirectory(Setting.TempDirectory);
        }
Exemplo n.º 6
0
        private FcbOfUpdWriter WriteFinal(FcbOfUpd bsdProduct)
        {
            //写最终FCB 产品
            var            fcbPath = Path.Combine(this.OutputDirectory, "FcbProduct" + Setting.FcbExtension);
            FcbOfUpdWriter writer  = new FcbOfUpdWriter(fcbPath);

            writer.Write(bsdProduct);
            writer.Dispose();
            return(writer);
        }
Exemplo n.º 7
0
        /// <summary>
        /// 当前历元产品,为所在时段的平均。
        /// </summary>
        /// <param name="basePrn">基准星</param>
        /// <param name="epoch">历元</param>
        /// <returns></returns>
        public FcbOfUpd GetWideLaneFcb(SatelliteNumber basePrn, Time epoch, double maxRms = 0.6)
        {
            //获取各测站内部差值。
            Dictionary <SatelliteNumber, RmsedNumeral> result = GetWideLaneFcbDic(basePrn, epoch, maxRms);

            if (result.Count == 0)
            {
                return(null);
            }
            FcbOfUpd fcb = new FcbOfUpd(basePrn, epoch, result);

            return(fcb);
        }
Exemplo n.º 8
0
        /// <summary>
        /// 将所有卫星产品归算到一颗卫星上面。
        /// </summary>
        /// <param name="basePrn"></param>
        /// <param name="intervalSeconds"></param>
        /// <param name="maxRms"></param>
        /// <returns></returns>
        public List <FcbOfUpd> GetWideLaneFcbOfAllSatAverage(SatelliteNumber basePrn, double intervalSeconds = 30, double maxRms = 0.6)
        {
            var timePeriod = this.TimePeriod;
            var list       = new List <FcbOfUpd>();

            for (Time epoch = timePeriod.Start; epoch <= timePeriod.End; epoch += intervalSeconds)
            {
                FcbOfUpd fcb = GetWideLaneFcbOfAllSat(basePrn, epoch, maxRms);

                list.Add(fcb);
            }
            return(list);
        }
Exemplo n.º 9
0
        /// <summary>
        /// 获取产品
        /// </summary>
        /// <param name="basePrn">基准卫星,并不做计算,只是赋值</param>
        /// <returns></returns>
        public List <FcbOfUpd> GetFcbProduct(SatelliteNumber basePrn, double maxRms = 0.6)
        {
            List <FcbOfUpd> result = new List <FcbOfUpd>();

            foreach (var epcochRow in this.Data)
            {
                var resultItem = new FcbOfUpd(basePrn, epcochRow.Key, false);
                result.Add(resultItem);

                foreach (var item in epcochRow.Value.KeyValues)
                {
                    if (item.Value.Rms > maxRms)
                    {
                        continue;
                    }

                    resultItem.Add(item.Key, item.Value);
                }
            }
            return(result);
        }