示例#1
0
        public ge_coords Transform(ge_coords _from)
        {
            ge_coords _to = new ge_coords();

            _to.coord_x = ge_UTMconversion.Helmert_X(_from.coord_x,
                                                     _from.coord_y,
                                                     _from.coord_z,
                                                     _t.trans_x,
                                                     _t.rot_y,
                                                     _t.rot_z,
                                                     _t.scale_change);

            _to.coord_y = ge_UTMconversion.Helmert_Y(_from.coord_x,
                                                     _from.coord_y,
                                                     _from.coord_z,
                                                     _t.trans_y,
                                                     _t.rot_x,
                                                     _t.rot_z,
                                                     _t.scale_change);

            _to.coord_z = ge_UTMconversion.Helmert_Z(_from.coord_x,
                                                     _from.coord_y,
                                                     _from.coord_z,
                                                     _t.trans_z,
                                                     _t.rot_x,
                                                     _t.rot_y,
                                                     _t.scale_change);
            return(_to);
        }
示例#2
0
        public ge_data WGS84_Lat_Long_Height_to_OSGB_East_North(ge_data _wg84)
        {
            // Step 1 Convert WGS84 Latitude, longitude and Ellipsoidal height to WGS84 Cartesian XYZ

            wg84 = _wg84;

            ge_projectionWGS84 pwg84 = new ge_projectionWGS84(wg84);

            if (!pwg84.calcXYZ_fromLatLong())
            {
                errorMsg = "Error pwg84.calcXYZ_fromLatLong()";
                return(null);
            }
            ;


            // Step 2 Apply Helmert Datum Transformation (WGS84 to OSGB36)

            ge_projection_transform t  = new WGS84_to_OSGB();
            HelmertTransform        ht = new HelmertTransform(t);
            ge_coords wg84_coords      = new ge_coords(wg84);
            ge_coords osgb_coords      = ht.Transform(wg84_coords);

            // Step 3 Convert OSGB36 Cartesian XYZ to OSGB36 Latitude, longitude and approx ODN height

            osg36      = new ge_data();
            osg36.locX = osgb_coords.coord_x;
            osg36.locY = osgb_coords.coord_y;
            osg36.locZ = osgb_coords.coord_z;

            ge_projectionOSGB36 posgb36 = new ge_projectionOSGB36(osg36);

            if (!posgb36.calcLatLongH_fromXYZ())
            {
                errorMsg = "Error posgb36.calcLatLongH_fromXYZ()";
                return(null);
            }
            ;

            // Step 4 Convert OSGB36 Latitude and longitude to OSGB36 easting and northing

            if (!posgb36.calcEN_fromLatLong())
            {
                errorMsg = "Error posgb36.calcEN_fromLatLong()";
                return(null);
            }

            return(osg36);
        }
示例#3
0
        public ge_data OSGB_East_North_to_WGS84_Lat_Long_Height(ge_data _osg36)
        {
            ge_data wg84 = new ge_data();

            // Step 1 Convert OSGB36 easting, northing and height to OSGB36 latitude and longitude

            osg36 = _osg36;

            ge_projectionOSGB36 posgb36 = new ge_projectionOSGB36(osg36);

            if (!posgb36.calcLatLong_fromEN())
            {
                errorMsg = "Error: posgb36.calcLatLong_fromEN()";
                return(null);
            }

            // Step 2 Convert OSGB36 latitude, longitude and approx ODN height to OSGB36 cartesian XYZ

            if (!posgb36.calcXYZ_fromLatLong())
            {
                errorMsg = "Error: posgb36.calcXYZ_fromLatLong()";
                return(null);
            }

            // Step 3 Apply Helmert Datum Transformation (OSGB36 to WGS84)
            ge_projection_transform t  = new OSGB_to_WGS84();
            HelmertTransform        ht = new HelmertTransform(t);
            ge_coords osg36_coords     = new ge_coords(osg36);
            ge_coords wg84_coords      = ht.Transform(osg36_coords);

            wg84.locX = wg84_coords.coord_x;
            wg84.locY = wg84_coords.coord_y;
            wg84.locZ = wg84_coords.coord_z;

            // Step Convert WGS84 Cartesian XYZ to WGS84 Latitude, longitude and Ellipsoidal height

            ge_projectionWGS84 pwg84 = new ge_projectionWGS84(wg84);

            if (!pwg84.calcLatLongH_fromXYZ())
            {
                errorMsg = "Error: pwg84.calcLatLongH_fromXYZ()";
                return(null);
            }
            ;

            return(wg84);
        }