Exemple #1
0
        /// <summary>
        /// 目标坐标转源坐标(批量)
        /// </summary>
        /// <param name="targetCoordiantes">目标坐标集合</param>
        /// <returns>源坐标集合</returns>
        public List <ICoordinate> TargetToSourceBatch(IEnumerable <ICoordinate> targetCoordiantes)
        {
            if (targetCoordiantes == null || targetCoordiantes.Count() == 0)
            {
                return(null);
            }
            List <ICoordinate> result = new List <ICoordinate>();
            // 两个中间变量
            SpherePoint spTemp = null;
            PlanePoint  ppTemp = null;
            // 高斯转换
            GaussKrugerTransform gauss_source = new GaussKrugerTransform(sourceCS);
            GaussKrugerTransform gauss_target = new GaussKrugerTransform(targetCS);
            // 大地坐标转换
            GeodeticTransform geo_source = new GeodeticTransform(sourceCS, sourceMeridian);
            GeodeticTransform geo_target = new GeodeticTransform(targetCS, targetMeridian);
            // 七参数模型,对空间直角坐标进行转换,转换后同样是空间直角坐标
            BursaWolfTransform bursa_target = new BursaWolfTransform(sevenParams.Reverse());

            foreach (ICoordinate targetCoordinate in targetCoordiantes)
            {
                // 如果目标是平面坐标
                // 1.平面坐标转高斯坐标;2.高斯反算,转为球面坐标
                if (targetCT == CoordinateType.Plane)
                {
                    ppTemp = gauss_target.PlaneToGauss((PlanePoint)targetCoordinate);
                    spTemp = gauss_target.GaussKrugerReverse(ppTemp, targetMeridian);
                }
                else
                {
                    spTemp = (SpherePoint)targetCoordinate.Clone();
                }

                // 大地坐标转空间直角坐标
                ppTemp = geo_target.GeodeticToThreeDimensions(spTemp);

                // 七参数模型,对空间直角坐标进行转换,转换后同样是空间直角坐标
                ppTemp = (PlanePoint)bursa_target.Transform(ppTemp);

                // 空间直角坐标转大地坐标
                spTemp = geo_source.ThreeDimensionsToGeodetic(ppTemp);

                // 如果源是平面坐标
                // 1.高斯正算,转为高斯坐标;2.高斯坐标转平面坐标
                if (sourceCT == CoordinateType.Plane)
                {
                    ppTemp = gauss_source.GaussKrugerForward(spTemp, sourceMeridian);
                    ppTemp = gauss_source.GaussToPlane(ppTemp);

                    result.Add(ppTemp.Clone());
                }
                else
                {
                    result.Add(spTemp.Clone());
                }
            }

            return(result);
        }
Exemple #2
0
        /// <summary>
        /// 源坐标转目标坐标
        /// </summary>
        /// <param name="sourceCoordinate">源坐标</param>
        /// <returns>目标坐标</returns>
        public ICoordinate SourceToTarget(ICoordinate sourceCoordinate)
        {
            if (sourceCoordinate == null)
            {
                return(null);
            }
            // 两个中间变量
            SpherePoint spTemp = null;
            PlanePoint  ppTemp = null;

            // 如果源是平面坐标
            // 1.平面坐标转高斯坐标;2.高斯反算,转为球面坐标
            if (sourceCT == CoordinateType.Plane)
            {
                GaussKrugerTransform gauss_source = new GaussKrugerTransform(sourceCS);
                ppTemp = gauss_source.PlaneToGauss((PlanePoint)sourceCoordinate);
                spTemp = gauss_source.GaussKrugerReverse(ppTemp, sourceMeridian);
            }
            else
            {
                spTemp = (SpherePoint)sourceCoordinate.Clone();
            }

            // 大地坐标转空间直角坐标
            GeodeticTransform geo_source = new GeodeticTransform(sourceCS, sourceMeridian);

            ppTemp = geo_source.GeodeticToThreeDimensions(spTemp);

            // 七参数模型,对空间直角坐标进行转换,转换后同样是空间直角坐标
            BursaWolfTransform bursa_source = new BursaWolfTransform(sevenParams);

            ppTemp = (PlanePoint)bursa_source.Transform(ppTemp);

            // 空间直角坐标转大地坐标
            GeodeticTransform geo_target = new GeodeticTransform(targetCS, targetMeridian);

            spTemp = geo_target.ThreeDimensionsToGeodetic(ppTemp);


            ICoordinate result = null;

            // 如果目标是平面坐标
            // 1.高斯正算,转为高斯坐标;2.高斯坐标转平面坐标
            if (targetCT == CoordinateType.Plane)
            {
                GaussKrugerTransform gauss_target = new GaussKrugerTransform(targetCS);
                ppTemp = gauss_target.GaussKrugerForward(spTemp, targetMeridian);
                ppTemp = gauss_target.GaussToPlane(ppTemp);

                result = ppTemp.Clone();
            }
            else
            {
                result = spTemp.Clone();
            }

            return(result);
        }