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
|
proj for DEBIAN
---------------
This is unusual, but let me provide some QuickStart on using this stuff.
proj
~~~~
I use `proj' all the time to draw maps in various projections using
graphing tools that cannot themselves handle fancy map projections
by piping all my data (coastline, etc) through proj.
For example, drawing Eastern Canada in Lambert Conformal Conic, I'd use:
proj +proj=lcc +ellps=clrk66 +lat_o=48.5 +lon_0="-62.5" -m 1:111120 -f '%.5f' -t9
invproj +proj=lcc +ellps=clrk66 +lat_o=48.5 +lon_0="-62.5" -m 1:111120 -f '%.5f'
-r : reverses the order of the expected input from longitude-latitude or
x-y to latitude-longitude or y-x.
-f : Format is a printf format string to control the form of the output
values.
-m : scale
-t : specifies a character employed as the first character to denote a
control line to be passed through without processing. This option
applicable to ascii input only. (# is the default value).
I use this to skip over the missing data parameter.
I pass Lat&Lon data segment to proj like so:
-65 49
-65.1 49.1
99 99
-66 49
-66 49.1
using 99 in this example to separate segments (missing data). After piping
the above through proj, I get :
$ proj +proj=lcc +lat_o=48.5 +lon_0="-62.5" -m 1:111120 +ellps=clrk66 -f '%.5f' l.tmp
-1.66314 51.62596
-1.72676 51.72892
* *
-2.32812 51.64792
-2.32423 51.74899
So, I might tell my graph-making package to use 99 as missing data
parameter, and tell proj "-t9" (the first character of my missing value
parameter) and get:
$ proj +proj=lcc +lat_o=48.5 +lon_0="-62.5" -m 1:111120 -t9 ellps=clrk66 -f '%.5f' l.tmp
-1.66314 51.62596
-1.72676 51.72892
99 99
-2.32812 51.64792
-2.32423 51.74899
geod
~~~~
geod is useful for calculating a latitude & longitude given another
position, along with distance and bearing. The man page has a good
example. It's also useful to calculate the distance and bearing between two
locations. Let's verify that one nautical mile is (approx) equal to one
minute in latitude (between 42 degN and 42 degN 1 minute):
$ geod -I +ellps=clrk66 +units=kmi <<EOF
42d00 72d00 42d01 72d00
42 72 42.016666 72
EOF
It outputs two lines, showing that it's pretty liberal about input data
format:
0d 180d 1.000
0d 180d 1.000
So 1 nm at 180deg bearing.
nad2nad and cs2cs:
~~~~~~~~~~~~~~~~~
The Debian proj package includes built-in datum conversion. (It was
compiled with proj-datumgrid-1.4.zip unpacked in the nad source directory.)
nad2nad is a filter to convert data between North America Datum 1927
(NAD27) and North American Datum 1983. nad2nad can optionally process both
State Plane Coordinate System (SPCS) and Universal Transverse Mercator
(UTM) grid data as well as geographic data for both input and output.
cs2cs is more general. It performs transformation between the source and
destination cartographic coordinate system on a set of input points. The
coordinate system transformation can include translation between projected
and geographic coordinates as well as the application of datum shifts.
Test nad2nad with the following example:
$ nad2nad -i 83 -o 27 -r conus << EOF
71d14'58.27"W 44d20'15.227"N
EOF
It should produce: 71d15'W 44d20'15"N
While the nad2nad program can be used in some cases, the cs2cs is now the
preferred mechanism. The following example demonstrates using the default
shift parameters for NAD27 to NAD83:
$ cs2cs +proj=latlong +datum=NAD27 +to +proj=latlong +datum=NAD83
-117 30
producing:
117d0'2.901"W 30d0'0.407"N 0.00
See http://trac.osgeo.org/proj/wiki/FAQ
Peter S Galbraith <psg@debian.org>
Francesco P Lovergine <frankie@debian.org>
Updated May 2009
|
