Exemple #1
0
        public static void SetDistanceAndAngleGeocacheFromLocation(GAPPSF.Core.Data.Geocache gc, GAPPSF.Core.Data.Location loc)
        {
            GeodeticMeasurement gm = CalculateDistance(loc.Lat, loc.Lon, gc.Lat, gc.Lon);

            gc.DistanceToCenter = (long)gm.EllipsoidalDistance;
            gc.AngleToCenter    = (int)gm.Azimuth.Degrees;
        }
Exemple #2
0
        public static GeodeticMeasurement CalculateDistance(double lat1, double lon1, double lat2, double lon2)
        {
            GlobalCoordinates   p1  = new GlobalCoordinates(new Angle(lat1), new Angle(lon1));
            GlobalCoordinates   p2  = new GlobalCoordinates(new Angle(lat2), new Angle(lon2));
            GeodeticCalculator  gc  = new GeodeticCalculator();
            GlobalPosition      gp1 = new GlobalPosition(p1);
            GlobalPosition      gp2 = new GlobalPosition(p2);
            GeodeticMeasurement gm  = gc.CalculateGeodeticMeasurement(Ellipsoid.WGS84, gp1, gp2);

            return(gm);
        }
Exemple #3
0
        public override object Execute(object[] args, ExecutionContext ctx)
        {
            string          res     = "";
            ArgumentChecker checker = new ArgumentChecker(this.GetType().Name);

            checker.CheckForNumberOfArguments(ref args, 2, null);
            Utils.Location ll1 = Utils.Conversion.StringToLocation(args[0].ToString());
            Utils.Location ll2 = Utils.Conversion.StringToLocation(args[1].ToString());
            if ((ll1 != null) && (ll2 != null))
            {
                GeodeticMeasurement gm = Utils.Calculus.CalculateDistance(ll1.Lat, ll1.Lon, ll2.Lat, ll2.Lon);
                res = gm.Azimuth.Degrees.ToString("0");
            }
            return(res);
        }
Exemple #4
0
        public static double ComputeDistance(Waypoint w1, Waypoint w2, out double elevationChange)
        {
            GlobalCoordinates startCoords = new GlobalCoordinates(w1.Latitude, w1.Longitude);
            GlobalCoordinates endCoords   = new GlobalCoordinates(w2.Latitude, w2.Longitude);
            GeodeticCurve     curve       = geoCal.CalculateGeodeticCurve(Ellipsoid.WGS84, startCoords, endCoords);

            elevationChange = 0;
            if (!double.IsNaN(w1.Elevation) && !double.IsNaN(w2.Elevation))
            {
                elevationChange = w2.Elevation - w1.Elevation;
            }
            GeodeticMeasurement measurement = new GeodeticMeasurement(curve, elevationChange);

            return(measurement.PointToPointDistance);
        }
Exemple #5
0
        public IList <DerivedPointDto> CreateArc(double radius, SuppliedPointDto centerPoint, SuppliedPointDto startPoint, SuppliedPointDto endPoint)
        {
            Angle centerLatitude  = new Angle(centerPoint.Latitude);
            Angle centerLongitude = new Angle(centerPoint.Longitude);
            Angle startLatitude   = new Angle(startPoint.Latitude);
            Angle startLongitude  = new Angle(startPoint.Longitude);
            Angle endLatitude     = new Angle(endPoint.Latitude);
            Angle endLongitude    = new Angle(endPoint.Longitude);

            GlobalCoordinates centerCoordinate = new GlobalCoordinates(centerLatitude, centerLongitude);
            GlobalCoordinates startCoordinate  = new GlobalCoordinates(startLatitude, startLongitude);
            GlobalCoordinates endCoordinate    = new GlobalCoordinates(endLatitude, endLongitude);

            GlobalPosition centerPosition = new GlobalPosition(centerCoordinate);
            GlobalPosition startPosition  = new GlobalPosition(startCoordinate);
            GlobalPosition endPosition    = new GlobalPosition(endCoordinate);

            GeodeticMeasurement centerToStartMeasurement = calculator.CalculateGeodeticMeasurement(centerPosition, startPosition);
            GeodeticMeasurement centerToEndMeasurement   = calculator.CalculateGeodeticMeasurement(centerPosition, endPosition);

            // Radius should ensure that it is tolerable (we are going to ingore it in this implementation)
            double startDistance   = centerToStartMeasurement.PointToPointDistance;
            double endDistance     = centerToEndMeasurement.PointToPointDistance;
            double distanceDelta   = Math.Abs(startDistance - endDistance);
            double distanceEpsilon = distanceDelta / 10;

            Angle startBearing = centerToStartMeasurement.Azimuth;
            Angle endBearing   = centerToEndMeasurement.Azimuth;

            double degreesDelta   = Math.Abs(startBearing.Degrees - endBearing.Degrees);
            double degreesEpsilon = degreesDelta / 10;  // Our angular distribution

            IList <DerivedPointDto> arcPoints = new List <DerivedPointDto>();

            // We are going to assume that arc derivation is done in a clockwise position.
            for (int i = 0; i < 10; i++) // Such that the 10th increment should equate to the end point and the initial the start
            {
                Angle             nextBearing   = new Angle(startBearing.Degrees + (degreesEpsilon * (i + 1)));
                GlobalCoordinates thetaLocation = calculator.CalculateEndingGlobalCoordinates(centerCoordinate, nextBearing, startDistance + (distanceEpsilon * (i + 1)));
                DerivedPointDto   thetaPoint    = new DerivedPointDto(thetaLocation.Latitude.Degrees, thetaLocation.Longitude.Degrees);

                arcPoints.Add(thetaPoint);
            }

            return(arcPoints);
        }
Exemple #6
0
 private void calculateAngleAndDegrees()
 {
     Framework.Data.Location ll1 = Utils.Conversion.StringToLocation(textBox5.Text);
     Framework.Data.Location ll2 = Utils.Conversion.StringToLocation(textBox6.Text);
     if (ll1 != null && ll2 != null)
     {
         GeodeticMeasurement gm = Utils.Calculus.CalculateDistance(ll1.Lat, ll1.Lon, ll2.Lat, ll2.Lon);
         textBox7.Text = gm.EllipsoidalDistance.ToString("0");
         textBox8.Text = (0.0006214 * gm.EllipsoidalDistance).ToString("0.000");
         textBox9.Text = gm.Azimuth.Degrees.ToString("0");
     }
     else
     {
         textBox7.Text = "";
         textBox8.Text = "";
         textBox9.Text = "";
     }
 }
        public void TestCalculateGeodeticMeasurement()
        {
            // instantiate the calculator
            GeodeticCalculator geoCalc = new GeodeticCalculator();

            // select a reference elllipsoid
            Ellipsoid reference = Ellipsoid.WGS84;

            // set Pike's Peak position
            GlobalPosition pikesPeak = new GlobalPosition(new GlobalCoordinates(Angle.FromDegrees(38.840511), Angle.FromDegrees(-105.0445896)), 4301);

            // set Alcatraz Island coordinates
            GlobalPosition alcatrazIsland = new GlobalPosition(new GlobalCoordinates(Angle.FromDegrees(37.826389), Angle.FromDegrees(-122.4225)), 0);

            // calculate the geodetic measurement
            GeodeticMeasurement geoMeasurement = geoCalc.CalculateGeodeticMeasurement(reference, pikesPeak, alcatrazIsland);

            Assert.AreEqual(-4301, geoMeasurement.ElevationChangeMeters, 0.001);
            Assert.AreEqual(1521788.826, geoMeasurement.PointToPointDistanceMeters, 0.001);
            Assert.AreEqual(1521782.748, geoMeasurement.EllipsoidalDistanceMeters, 0.001);
            Assert.AreEqual(271.21039153, geoMeasurement.Azimuth.Degrees, StandardTolerance);
            Assert.AreEqual(80.38029386, geoMeasurement.ReverseAzimuth.Degrees, StandardTolerance);
        }