GeoDestination
GeoDestination[loc,{d,α}]
gives the end position of the geodesic of length d starting from loc with azimuthal direction α.
Details
- The location loc in GeoDestination[loc,…] can be specified as {lat,lon} coordinates, as a GeoPosition, GeoPositionXYZ, GeoPositionENU, or GeoGridPosition object, or as a geographical entity Entity[…].
- The distance d can be given as a Quantity length or as a number in meters.
- The bearing or azimuthal direction α is an angle measured clockwise from true north. It can be given as a Quantity angle or as a number in degrees.
- GeoDestination[loc,{{d1,d2,…},α}] returns the list of positions at distances di along the geodesic starting from loc with bearing α.
- GeoDestination[loc,{d,{α1,α2,…}}] returns the list of points at distance d along the geodesics starting from loc with bearings αi.
- GeoDestination[loc,{d,α}] can also be written as GeoDestination[loc,GeoDisplacement[{d,α}]].
- Use GeoDestination[loc,GeoDisplacement[{d,α},"Rhumb"]] to find the endpoint along a path of length d and constant bearing α.
- GeoDestination[{lat,lon},…] gives the latitude-longitude destination for the default reference ellipsoid.
- GeoDestination[loc,…] finds the destination on the reference ellipsoid associated with the datum for loc. Heights are ignored.
- GeoDestination solves the geodetic direct or forward problem.
Examples
open all close allBasic Examples (2)
Compute the position reached by moving 100 kilometers north from a starting position:
GeoDestination[GeoPosition[{37., -109.}], {100000, 0}]Move 4000 kilometers along a geodesic with initial bearing 105 degrees:
p = GeoPosition[Entity["City", {"Seattle", "Washington", "UnitedStates"}]]GeoDestination[p, {Quantity[4000, "Kilometers"], Quantity[105, "AngularDegrees"]}]GeoGraphics[Arrow[GeoPath[{p, %}]]]Scope (9)
Specify the initial position as a {lat,lon} pair in degrees:
GeoDestination[{40.385, -30.216}, {10 ^ 6, 45}]GeoDestination preserves the type of geodetic position:
pos = GeoPosition[{40.385, -30.216}]GeoDestination[%, {10 ^ 6, 45}]Transform to GeoPositionXYZ:
GeoPositionXYZ[pos]GeoDestination[%, {10 ^ 6, 45}]Transform to GeoGridPosition. The destination is returned in the same projection:
GeoGridPosition[pos, "Mercator"]GeoDestination[%, {10 ^ 6, 45}]Transform to GeoPositionENU. The destination is given with respect to the same ENU origin:
GeoPositionENU[pos, GeoPosition[{40, -30}]]GeoDestination[%, {10 ^ 6, 45}]Use DMS strings to specify the original position or the path to follow:
GeoDestination[GeoPosition[{"40°23'6''", "30°12'57.6''W"}], GeoDisplacement[{10 ^ 6, "45d 32m 10s"}]]Use a compass point specification of the initial bearing:
GeoDestination[GeoPosition[{"40°23'6''", "30°12'57.6''W"}], GeoDisplacement[{10 ^ 6, "ESE"}]]Find destination along a geodesic:
GeoDestination[GeoPosition[{40.385, -30.216}], {10 ^ 6, 45}]GeoDestination[GeoPosition[{40.385, -30.216}], GeoDisplacement[{10 ^ 6, 45}, "Geodesic"]]Find destination along a line of constant bearing:
GeoDestination[GeoPosition[{40.385, -30.216}], GeoDisplacement[{10 ^ 6, 45}, "Rhumb"]]Height and time information is returned unmodified:
GeoDestination[GeoPosition[{40.385, -30.216, 100.}], {10 ^ 6, 45}]GeoDestination[GeoPosition[{40.385, -30.216, 100., 3.6 10 ^ 9}], {10 ^ 6, 45}]Computations are performed using the specified datum:
GeoDestination[GeoPosition[{40.385, -30.216, 100.}, "WGS72"], {10 ^ 6, 45}]GeoDestination[GeoPosition[{40.385, -30.216, 100.}, "ITRF00"], {10 ^ 6, 45}]Differences are small, in this case under 1 arc second:
DMSString /@ (QuantityMagnitude@LatitudeLongitude[%%] - QuantityMagnitude@LatitudeLongitude[%])Specify a list of distances to obtain a list of points along the same geodesic:
GeoDestination[GeoPosition[{0, 0}], {Range[0, 10 ^ 6, 10 ^ 5], 45}]Specify a list of bearings to obtain a list of points on a geodesic circle:
GeoDestination[GeoPosition[{0, 0}], {10 ^ 6, Range[0, 90, 15]}]Specify both a list of distances and a list of bearings. The result is a matrix of points:
GeoDestination[GeoPosition[{0, 0}], {Range[0, 10 ^ 6, 10 ^ 6 / 3], Range[0, 90, 30]}]Properties & Relations (4)
The inverse of GeoDestination is GeoDisplacement:
origin = GeoPosition[{50, 130}];destination = GeoDestination[origin, {5 10 ^ 6, 60}]GeoDisplacement[origin, destination]It can also be obtained as a combination of GeoDistance and GeoDirection:
GeoDistance[origin, destination]UnitConvert[%, "Meters"]GeoDirection[origin, destination]Compute multiple points along the same geodesic:
GeoDestination[{50, 150}, {Range[0, 10 ^ 7, 10 ^ 6], -75}]GeoGraphics[{Red, Line[%]}]Compute points on a geodesic circle of 2000 kilometers around Stockholm:
stockholm = Entity["City", {"Stockholm", "Stockholm", "Sweden"}];
distance = ["2000 km"];circle = GeoDestination[stockholm, {distance, Range[0, 360, 6]}];Compare with the corresponding GeoCircle:
GeoGraphics[{GeoCircle[stockholm, distance], Blue, Point[stockholm], Red, Point[circle]}]Compute consecutive destination points, all moving 1000 kilometers:
d = Quantity[1000, "Kilometers"];course = Table[{d, α}, {α, 0, 315, 45}]FoldList[GeoDestination, GeoPosition[{50, -100}], course]Due to the Earth surface's curvature, the path does not close:
GeoGraphics[{Red, Line[%]}]Neat Examples (1)
Points 20000 kilometers from geo position {0,0} in all directions:
points = GeoDestination[GeoPosition[{0, 0}], {20000000, Range[0, 360, 4]}];GeoGraphics[Point[points], GeoRange -> {{-0.5, 0.5}, {179.5, 180.5}}, Frame -> True]a = GeodesyData["ITRF00", "SemimajorAxis"];
b = GeodesyData["ITRF00", "SemiminorAxis"];
{a, b, c = a - b}ListPointPlot3D[GeoPositionXYZ[points]["XYZ"], PlotRange -> QuantityMagnitude[{{-a - 2c, -a + 2c}, {-2c, 2c}, {-c, c}}, "Meters"], Axes -> True, BoxRatios -> 1]Text
Wolfram Research (2008), GeoDestination, Wolfram Language function, https://reference.wolfram.com/language/ref/GeoDestination.html (updated 2014).
CMS
Wolfram Language. 2008. "GeoDestination." Wolfram Language & System Documentation Center. Wolfram Research. Last Modified 2014. https://reference.wolfram.com/language/ref/GeoDestination.html.
APA
Wolfram Language. (2008). GeoDestination. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/GeoDestination.html