1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
|
<HTML><HEAD><TITLE>Manpage of GEOD</TITLE>
</HEAD><BODY>
<H1>GEOD</H1>
Section: User Commands (1)<BR>Updated: 2000/03/21 Rel. 4.4<BR><A HREF="#index">Index</A>
<A HREF="../index.html">Return to Main Contents</A><HR>
<A NAME="lbAB"> </A>
<H2>NAME</H2>
geod - direct geodesic computations
<BR>
invgeod - inverse geodesic computations
<A NAME="lbAC"> </A>
<H2>SYNOPSIS</H2>
<B>geod</B>
<B>+ellps=<ellipse></B>
[
<B>-afFIlptwW</B>
[
<I>args</I>
] ] [
<B>+args</B>
]
file[s]
<BR>
<B>invgeod</B>
<B>+ellps=<ellipse></B>
[
<B>-afFIlptwW</B>
[
<I>args</I>
] ] [
<B>+args</B>
]
file[s]
<A NAME="lbAD"> </A>
<H2>DESCRIPTION</H2>
<I>Geod</I>
(direct) and
<I>invgeod</I>
(inverse)
perform geodesic (``Great Circle'') computations for determining
latitude, longitude and back azimuth of a terminus point
given a initial point latitude, longitude, azimuth and distance (direct) or
the forward and back azimuths and distance between an initial and
terminus point latitudes and longitudes (inverse).
<P>
The following runline control parameters can appear in any order:
<DL COMPACT>
<DT><B>-I</B>
<DD>
Specifies that the inverse geodesic computation is to be performed.
May be used with execution of
<B>goed</B>
as an alternative to
<B>invgeod</B>
execution.
<DT><B>-a</B>
<DD>
Latitude and longitudes of the initial and terminal points,
forward and back azimuths and distance are output.
<DT><B>-t</B><I>a</I>
<DD>
<I>A</I>
specifies a character employed as the first character to denote
a control line to be passed through without processing.
<DT><B>-le</B>
<DD>
Gives a listing of all the ellipsoids that may be selected with the
<B>+ellps=</B>
option.
<DT><B>-lu</B>
<DD>
Gives a listing of all the units that may be selected with the
<B>+units=</B>
option.
<DT><B>-[f|F]</B><I> format</I>
<DD>
<I>Format</I>
is a
<I>printf</I>
format string to control the output form of the geographic coordinate values
(<B>f</B>) or distance value (<B>F</B>).
The default mode is DMS for geographic coordinates and "%.3f" for distance.
<DT><B>-[w|W]</B><I>n</I>
<DD>
<I>N</I>
is the number of significant fractional digits to employ for
seconds output (when the option is not specified,
<B>-w3</B>
is assumed).
When
<B>-W</B>
is employed the fields will be constant width with leading zeroes.
<DT><B>-p</B>
<DD>
This option causes the azimuthal values to be output as unsigned
DMS numbers between 0 and 360 degrees. Also note -f.
</DL>
<P>
The
<B>+args</B>
run-line arguments are associated with geodetic parameters
for specifying the ellipsoidal or sphere to use.
See
<B>proj</B>
documentation for full list of these parameters and contrl.
The options are processed in left to right order
from the run line.
Reentry of an option is ignored with the first occurrence assumed to
be the desired value.
<P>
One or more
<I>files</I>
(processed in left to right order)
specify the source of data to be transformed.
A - will specify the location of processing standard input.
If no files are specified, the input is assumed to be from
<I>stdin.</I>
<P>
For direct determinations input data must be in latitude,
longitude, azimuth and distance order and output will be
latitude, longitude and back azimuth of the terminus point.
Latitude, longitude of the initial and terminus point are
input for the inverse mode and respective forward and back
azimuth from the initial and terminus points are output along
with the distance between the points.
<P>
Input geographic coordinates
(latitude and longitude) and azimuthal data must be in DMS format and input
distance data must be in units consistent with the ellipsoid
major axis or sphere radius units.
Output geographic coordinates will be in DMS
(if the
<B>-f</B>
switch is not employed) to 0.001"
with trailing, zero-valued minute-second fields deleted.
Output distance data will be in the same units as the ellipsoid or
sphere radius.
<P>
The Earth's ellipsoidal figure may be selected in the same
manner as program
<B>proj</B>
by using
<B>+ellps=, +a=, +es=,</B>
etc.
<P>
<I>Geod</I>
may also be used to determine intermediate points along either
a geodesic line between two points or along an arc of specified distance
from a geographic point.
In both cases an initial point must be specified with
<B>+lat_1=</B><I>lat</I>
and
<B>+lon_1=</B><I>lon</I>
parameters and either a terminus point
<B>+lat_2=</B><I>lat</I>
and
<B>+lon_2=</B><I>lon</I>
or a distance and azimuth from the initial point with
<B>+S=</B><I>distance</I>
and
<B>+A=</B><I>azimuth</I>
must be specified.
<P>
If points along a geodesic are to be determined then either
<B>+n_S=</B><I>integer</I>
specifying the number of intermediate points and/or
<B>+del_S=</B><I>distance</I>
specifying the incremental distance between points must be specified.
<P>
To determine points along an arc equidistant from the initial point both
<B>+del_A=</B><I>angle</I>
and
<B>+n_A=</B><I>integer</I>
must be specified which determine the respective angular increments
and number of points to be determined.
<A NAME="lbAE"> </A>
<H2>EXAMPLE</H2>
The following script determines the geodesic azimuths and distance in
U.S. stature miles from Boston, MA, to Portland, OR:
<DL COMPACT><DT><DD>
<BR> <TT>geod +ellps=clrk66 <<EOF -I +units=us-mi
<BR> 42d15'N 71d07'W 45d31'N 123d41'W
<BR> EOF</TT>
</DL>
which gives the results:
<DL COMPACT><DT><DD>
<BR> <TT>-66d31'50.141" 75d39'13.083" 2587.504
</DL>
where the first two values are the
azimuth from Boston to Portland, the back azimuth from Portland to
Boston followed by the distance.
<P>
An example of forward geodesic use is to use the Boston location and determine
Portland's location by azimuth and distance:
<DL COMPACT><DT><DD>
<BR> geod +ellps=clrk66 <<EOF +units=us-mi
<BR> 42d15'N 71d07'W -66d31'50.141" 2587.504
<BR> EOF</TT>
</DL>
which gives:
<DL COMPACT><DT><DD>
<BR> <TT>45d31'0.003"N 123d40'59.985"W 75d39'13.094"</TT>
</DL>
Note: lack of precision in the distance value compromises
the precision of the Portland location.
<A NAME="lbAF"> </A>
<H2>SEE ALSO</H2>
Thomas, P.D., 1970,
<I>Spheroidal Geodesics, Reference Systems & Local Geometry:</I>
U.S. Naval Oceanographic Office, S-138.
<A NAME="lbAG"> </A>
<H2>HOME PAGE</H2>
<A HREF="http://www.remotesensing.org/proj">http://www.remotesensing.org/proj</A>
<P>
<HR>
<A NAME="index"> </A><H2>Index</H2>
<DL>
<DT><A HREF="#lbAB">NAME</A><DD>
<DT><A HREF="#lbAC">SYNOPSIS</A><DD>
<DT><A HREF="#lbAD">DESCRIPTION</A><DD>
<DT><A HREF="#lbAE">EXAMPLE</A><DD>
<DT><A HREF="#lbAF">SEE ALSO</A><DD>
<DT><A HREF="#lbAG">HOME PAGE</A><DD>
</DL>
<HR>
This document was created by
<A HREF="http://localhost/cgi-bin/man/man2html">man2html</A>,
using the manual pages.<BR>
Time: 13:14:56 GMT, October 14, 2005
</BODY>
</HTML>
|
