5
.TH GEOD 1 "2000/03/21 Rel. 4.4"
7
geod \- direct geodesic computations
9
invgeod \- inverse geodesic computations
37
perform geodesic (\*(lqGreat Circle\*(rq) computations for determining
38
latitude, longitude and back azimuth of a terminus point
39
given a initial point latitude, longitude, azimuth and distance (direct) or
40
the forward and back azimuths and distance between an initial and
41
terminus point latitudes and longitudes (inverse).
43
The following runline control parameters can appear in any order:
46
Specifies that the inverse geodesic computation is to be performed.
47
May be used with execution of
54
Latitude and longitudes of the initial and terminal points,
55
forward and back azimuths and distance are output.
59
specifies a character employed as the first character to denote
60
a control line to be passed through without processing.
63
Gives a listing of all the ellipsoids that may be selected with the
68
Gives a listing of all the units that may be selected with the
76
format string to control the output form of the geographic coordinate values
77
(\fBf\fR) or distance value (\fBF\fR).
78
The default mode is DMS for geographic coordinates and "%.3f" for distance.
82
is the number of significant fractional digits to employ for
83
seconds output (when the option is not specified,
88
is employed the fields will be constant width with leading zeroes.
91
This option causes the azimuthal values to be output as unsigned
92
DMS numbers between 0 and 360 degrees. Also note -f.
96
run-line arguments are associated with geodetic parameters
97
for specifying the ellipsoidal or sphere to use.
100
documentation for full list of these parameters and controls.
101
The options are processed in left to right order
103
Reentry of an option is ignored with the first occurrence assumed to
104
be the desired value.
108
(processed in left to right order)
109
specify the source of data to be transformed.
110
A \- will specify the location of processing standard input.
111
If no files are specified, the input is assumed to be from
114
For direct determinations input data must be in latitude,
115
longitude, azimuth and distance order and output will be
116
latitude, longitude and back azimuth of the terminus point.
117
Latitude, longitude of the initial and terminus point are
118
input for the inverse mode and respective forward and back
119
azimuth from the initial and terminus points are output along
120
with the distance between the points.
122
Input geographic coordinates
123
(latitude and longitude) and azimuthal data must be in DMS format and input
124
distance data must be in units consistent with the ellipsoid
125
major axis or sphere radius units.
126
Output geographic coordinates will be in DMS
129
switch is not employed) to 0.001"
130
with trailing, zero-valued minute-second fields deleted.
131
Output distance data will be in the same units as the ellipsoid or
134
The Earth's ellipsoidal figure may be selected in the same
138
.B "+ellps=, +a=, +es=,"
142
may also be used to determine intermediate points along either
143
a geodesic line between two points or along an arc of specified distance
144
from a geographic point.
145
In both cases an initial point must be specified with
149
parameters and either a terminus point
153
or a distance and azimuth from the initial point with
159
If points along a geodesic are to be determined then either
161
specifying the number of intermediate points and/or
163
specifying the incremental distance between points must be specified.
165
To determine points along an arc equidistant from the initial point both
169
must be specified which determine the respective angular increments
170
and number of points to be determined.
173
The following script determines the geodesic azimuths and distance in
174
U.S. stature miles from Boston, MA, to Portland, OR:
176
\f(CWgeod +ellps=clrk66 <<EOF -I +units=us-mi
177
42d15'N 71d07'W 45d31'N 123d41'W
180
which gives the results:
182
\f(CW-66d31'50.141" 75d39'13.083" 2587.504
184
where the first two values are the
185
azimuth from Boston to Portland, the back azimuth from Portland to
186
Boston followed by the distance.
188
An example of forward geodesic use is to use the Boston location and determine
189
Portland's location by azimuth and distance:
191
\f(CWgeod +ellps=clrk66 <<EOF +units=us-mi
192
42d15'N 71d07'W -66d31'50.141" 2587.504
197
\f(CW45d31'0.003"N 123d40'59.985"W 75d39'13.094"\fR
199
Note: lack of precision in the distance value compromises
200
the precision of the Portland location.
203
.I "Spheroidal Geodesics, Reference Systems & Local Geometry:"
204
U.S. Naval Oceanographic Office, S-138.
206
http://www.remotesensing.org/proj