1GMTDIGITIZE(1) Generic Mapping Tools GMTDIGITIZE(1)
2
6 gmtdigitize - Digitizing and Inverse map transformation of map x/y
7 coordinates
8
10 gmtdigitize -Jparameters -Rwest/east/south/north[r] [ -A ] [ -Cdevice ]
11 [ -Dlimit ] [ -F ] [ -H[i][nrec] ] [ -Llpi ] [ -Nnamestem ] [ -S ] [ -V
12 ] [ -Zk|v ] [ -bo[s|S|d|D[ncol]|c[var1/...]] ] [ -m[flag] ] [ > out‐
13 put.d ]
14
16 gmtdigitize digitizes points from a digitizer via a serial line connec‐
17 tion and computes map coordinates using the specified map projection.
18 The program is interactive and will take you through the setup proce‐
19 dure and how you will digitize points. The program will determine the
20 actual map scale as well as rotation of the paper that is taped to the
21 digitizer table. By default the output will go to stdout.
22 No space between the option flag and the associated arguments.
23 Use upper case for the option flags and lower case for modifiers.
24 -J Selects the map projection. Scale is UNIT/degree, 1:xxxxx, or
26 width in UNIT (upper case modifier). UNIT is cm, inch, or m,
27 depending on the MEASURE_UNIT setting in .gmtdefaults4, but this
28 can be overridden on the command line by appending c, i, or m to
29 the scale/width value. When central meridian is optional,
30 default is center of longitude range on -R option. Default
31 standard parallel is the equator. For map height, max dimen‐
32 sion, or min dimension, append h, +, or - to the width, respec‐
33 tively.
34 More details can be found in the psbasemap man pages.
35 CYLINDRICAL PROJECTIONS:
37 -Jclon0/lat0/scale (Cassini)
39 -Jcyl_stere/[lon0/[lat0/]]scale (Cylindrical Stereographic)
40 -Jj[lon0/]scale (Miller)
41 -Jm[lon0/[lat0/]]scale (Mercator)
42 -Jmlon0/lat0/scale (Mercator - Give meridian and standard paral‐
43 lel)
44 -Jo[a]lon0/lat0/azimuth/scale (Oblique Mercator - point and
45 azimuth)
46 -Jo[b]lon0/lat0/lon1/lat1/scale (Oblique Mercator - two points)
47 -Joclon0/lat0/lonp/latp/scale (Oblique Mercator - point and
48 pole)
49 -Jq[lon0/[lat0/]]scale (Cylindrical Equidistant)
50 -Jtlon0/[lat0/]scale (TM - Transverse Mercator)
51 -Juzone/scale (UTM - Universal Transverse Mercator)
52 -Jy[lon0/[lat0/]]scale (Cylindrical Equal-Area)
53 CONIC PROJECTIONS:
55 -Jblon0/lat0/lat1/lat2/scale (Albers)
57 -Jdlon0/lat0/lat1/lat2/scale (Conic Equidistant)
58 -Jllon0/lat0/lat1/lat2/scale (Lambert Conic Conformal)
59 -Jpoly/[lon0/[lat0/]]scale ((American) Polyconic)
60 AZIMUTHAL PROJECTIONS:
62 -Jalon0/lat0[/horizon]/scale (Lambert Azimuthal Equal-Area)
64 -Jelon0/lat0[/horizon]/scale (Azimuthal Equidistant)
65 -Jflon0/lat0[/horizon]/scale (Gnomonic)
66 -Jglon0/lat0[/horizon]/scale (Orthographic)
67 -Jglon0/lat0/altitude/azimuth/tilt/twist/Width/Height/scale
68 (General Perspective).
69 -Jslon0/lat0[/horizon]/scale (General Stereographic)
70 MISCELLANEOUS PROJECTIONS:
72 -Jh[lon0/]scale (Hammer)
74 -Ji[lon0/]scale (Sinusoidal)
75 -Jkf[lon0/]scale (Eckert IV)
76 -Jk[s][lon0/]scale (Eckert VI)
77 -Jn[lon0/]scale (Robinson)
78 -Jr[lon0/]scale (Winkel Tripel)
79 -Jv[lon0/]scale (Van der Grinten)
80 -Jw[lon0/]scale (Mollweide)
81 NON-GEOGRAPHICAL PROJECTIONS:
83 -Jp[a]scale[/origin][r|z] (Polar coordinates (theta,r))
85 -Jxx-scale[d|l|ppow|t|T][/y-scale[d|l|ppow|t|T]] (Linear, log,
86 and power scaling)
87 For geographic projections you can give 1 as the scale will be
88 solved for anyway.
89 -R xmin, xmax, ymin, and ymax specify the Region of interest. For
91 geographic regions, these limits correspond to west, east,
92 south, and north and you may specify them in decimal degrees or
93 in [+-]dd:mm[:ss.xxx][W|E|S|N] format. Append r if lower left
94 and upper right map coordinates are given instead of w/e/s/n.
95 The two shorthands -Rg and -Rd stand for global domain (0/360
96 and -180/+180 in longitude respectively, with -90/+90 in lati‐
97 tude). Alternatively, specify the name of an existing grid file
98 and the -R settings (and grid spacing, if applicable) are copied
99 from the grid. For calendar time coordinates you may either
100 give (a) relative time (relative to the selected TIME_EPOCH and
101 in the selected TIME_UNIT; append t to -JX|x), or (b) absolute
102 time of the form [date]T[clock] (append T to -JX|x). At least
103 one of date and clock must be present; the T is always required.
104 The date string must be of the form [-]yyyy[-mm[-dd]] (Gregorian
105 calendar) or yyyy[-Www[-d]] (ISO week calendar), while the clock
106 string must be of the form hh:mm:ss[.xxx]. The use of delim‐
107 iters and their type and positions must be exactly as indicated
108 (however, input, output and plot formats are customizable; see
109 gmtdefaults).
110
112 -A Give an audible signal each time the digitizer mouse/puck is
113 clicked [Default is silent].
114 -C Specify the device (port) to read from [Default is /dev/ttyS0].
116 -D Only output a point if it is further than limit units from the
118 previous point. Append c, i, m, p for cm, inch, meter, or
119 point, respectively [Default is no limit].
120 -F Force the program to ask for 4 arbitrary calibration points
122 [Default is to use the 4 corners of the map, if possible].
123 -H This option allows you to write out any number of header records
125 to the beginning of the output file. Each record will automati‐
126 cally start with a #-character to indicate comment. Headers are
127 not written if multiple output files are selected with -N -m.
128 -L Set the digitizer table resolution in lines per inch [2540].
130 -N Set name for output file(s). If a regular filename is given,
132 then all digitized data will be written to that file. If the
133 file contains a C-format for an integer (i.e., %d) then the file
134 is used as a format statement to create unique filenames based
135 on the current segment number (e.g., line_%d.d will yield files
136 line_0.d, line_1.d, etc). By default, all output is written to
137 stdout. Multiple segment files requires specifying the -m
138 option.
139 -S Suppress points that fall outside the specified map region
141 [Default outputs all points].
142 -V Selects verbose mode, which will send progress reports to stderr
144 [Default runs "silently"]. The program will also duplicate data
145 output to stderr for monitoring.
146 -Z Append v to prompt for a z-value and output it as a third data
148 column. Append k to output the button key as the final data
149 column. Both -Zk and -Zv can be specified. [Default is just 2
150 column x,y output].
151 -bo Selects binary output. Append s for single precision [Default
153 is d (double)]. Uppercase S or D will force byte-swapping.
154 Optionally, append ncol, the number of desired columns in your
155 binary output file.
156 -m Multiple segment file(s). Segments are separated by a special
158 record. For ASCII files the first character must be flag
159 [Default is '>']. For binary files all fields must be NaN and
160 -b must set the number of output columns explicitly. By default
161 the -m setting applies to both input and output. Use -mi and
162 -mo to give separate settings to input and output.
163
165 To digitize lines from a mercator map made for a given region, and save
166 each line segment in individual files called segment_000.xy, seg‐
167 ment_001.xy etc, try
168 gmtdigitize -R20/50/12/25 -Jm1:1 -m -Nsegment_%3.3d.xy
170 To digitize seismically defined interfaces from a multichannel seismic
172 section, with horizontal distances from 130 to 970, and vertical times
173 from 0 to 10 seconds, write out the button code, and save all line seg‐
174 ment to a single multisegment file, and beep at each click, try
175 gmtdigitize -R130/970/0/10 -Jx1/-1 -m -A -Z > interfaces.d
177
179 This applies to the Calcomp DrawingBoard III hooked up to a RedHat
180 Linux workstation. We use /dev/ttyS0 as the serial port and change
181 permissions so that it is world read/write-able. Then, stty -F
182 /dev/ttyS0 evenp will set the terminal settings, which can be checked
183 with stty -F /dev/ttyS0 -a. Setup of digitizer: We use the CalComp
184 2000 ASCII (Save 3) setup, which has:
185 Mode: Point
186 Baud Rate: 9600
187 Data Bits: 7
188 Parity: Even
189 Data Rate: 125 pps
190 Resolution: 200 lpi
191 Output Format: Format 0
192 Emulation: CalComp 2000 ASCII
193 (A)We need to make a slight modification to the Preset No 3 settings:
194 (1) 2450 LPI instead of 200, and (2) None instead of yes for added CR.
195 These modifications can be changed and saved to Preset 3 on the digi‐
196 tizer but a power outage may reset in back to the factory defaults,
197 necessitating a manual reset of those two settings. (B) Setup tty
198 port. stty -F /dev/ttyS0 evenp (C) Run gmtdigitize. Map scale does
199 not matter; it is computed from the region and plot size.
200
202 gmtdefaults(l), GMT(l), gmtstitch(l), mapproject(l), project(l)
203GMT 4.5.6 10 Mar 2011 GMTDIGITIZE(1)