Beispiel #1
0
        /// <summary>
        /// 两个文件合并,本方法只适合计算站点不重复,且不相关的情况。
        /// 信息只是简单的叠加。只有两个文件都有的信息块才会被合并。
        /// </summary>
        /// <param name="fileB"></param>
        /// <param name="eraseNonCoord"> 是否清理非坐标的值和对应矩阵</param>
        public static SinexFile Merge(SinexFile fileA, SinexFile fileB, bool eraseNonCoord = true)
        {
            if (eraseNonCoord)
            {
                fileA.CleanNonCoordSolutionValue();
                fileB.CleanNonCoordSolutionValue();
            }

            SinexFile newFile = EmergeBasic(fileA, fileB);

            //合并解算结果,不更新标准差
            SinexMerger.MergeSolutionValue(newFile.SolutionEstimateBlock, fileA.SolutionEstimateBlock, fileB.SolutionEstimateBlock);
            SinexMerger.MergeSolutionValue(newFile.SolutionAprioriBlock, fileA.SolutionAprioriBlock, fileB.SolutionAprioriBlock);

            //统计数据的合并, 需要做一些叠加工作。
            SinexStatistic statisticA = fileA.GetStatistic();
            SinexStatistic statisticB = fileB.GetStatistic();

            if (fileA.SolutionStattisticsBlock != null && fileB.SolutionStattisticsBlock != null)
            {
                newFile.SolutionStattisticsBlock.Items = fileA.GetStatistic().Merge(fileB.GetStatistic()).GetSolutionStatistics();
            }
            SinexStatistic statisticNew = newFile.GetStatistic();
            double         varFactorA   = statisticNew.VarianceOfUnitWeight / statisticA.VarianceOfUnitWeight;
            double         varFactorB   = statisticNew.VarianceOfUnitWeight / statisticB.VarianceOfUnitWeight;


            //合并矩阵,协方差矩阵,前面已经判断并重新布置matrix了,此处只管合并。
            MergeSolutionMatrix(newFile.SolutionMatrixEstimateCova, fileA.SolutionMatrixEstimateCova, varFactorA, fileB.SolutionMatrixEstimateCova, varFactorB, newFile.SolutionEstimateBlock);
            MergeSolutionMatrix(newFile.SolutionMatrixAprioriCova, fileA.SolutionMatrixAprioriCova, varFactorA, fileB.SolutionMatrixAprioriCova, varFactorB, newFile.SolutionAprioriBlock);

            return(newFile);
        }