~ubuntu-branches/ubuntu/wily/grass/wily

« back to all changes in this revision

Viewing changes to imagery/i.vpoints/description.html

Tags: 7.0.0~rc1+ds1-1~exp1
* New upstream release candidate.
* Repack upstream tarball, remove precompiled Python objects.
* Add upstream metadata.
* Update gbp.conf and Vcs-Git URL to use the experimental branch.
* Update watch file for GRASS 7.0.
* Drop build dependencies for Tcl/Tk, add build dependencies:
  python-numpy, libnetcdf-dev, netcdf-bin, libblas-dev, liblapack-dev
* Update Vcs-Browser URL to use cgit instead of gitweb.
* Update paths to use grass70.
* Add configure options: --with-netcdf, --with-blas, --with-lapack,
  remove --with-tcltk-includes.
* Update patches for GRASS 7.
* Update copyright file, changes:
  - Update copyright years
  - Group files by license
  - Remove unused license sections
* Add patches for various typos.
* Fix desktop file with patch instead of d/rules.
* Use minimal dh rules.
* Bump Standards-Version to 3.9.6, no changes.
* Use dpkg-maintscript-helper to replace directories with symlinks.
  (closes: #776349)
* Update my email to use @debian.org address.

Show diffs side-by-side

added added

removed removed

Lines of Context:
1
 
<H2>DESCRIPTION</H2>
2
 
 
3
 
This program enables the user to identify coordinate pairs of points from a
4
 
vector map or keyboard entry and corresponding points in an image to be
5
 
rectified. The map coordinate values of each point are used to calculate a
6
 
transformation matrix.  The operator may then use the 
7
 
<EM><a href="i.rectify.html">i.rectify</A></EM> program to rectify the image using
8
 
the transformation matrix coefficients calculated from the control point
9
 
file created in <EM>i.vpoints</EM>.  The
10
 
<EM> <a href="i.rectify.html">i.rectify</A></EM> program performs a first, second or
11
 
third order transformation of the image.
12
 
 
13
 
<P>
14
 
The first step is to display the unrectified image and
15
 
corresponding vector map data.  The operator would then
16
 
mark corresponding control point locations on the image and
17
 
map. To identify the precise location of a point to be
18
 
marked, <EM>i.vpoints</EM> has a zoom option.  In addition to
19
 
marking control points on an image to be rectified and
20
 
inputting their world coordinate values using the
21
 
keyboard,  <EM>i.vpoints</EM> has the option to simultaneously
22
 
display vector map data available in the targeted database,
23
 
and identify on the vector map the location of the
24
 
corresponding marked points.  When this option is chosen,
25
 
the coordinate values are input automatically.  Any GRASS
26
 
map layer or vector map in the targeted database LOCATION
27
 
can be displayed using <EM>i.vpoints</EM>.  The <EM>i.vpoints</EM> program
28
 
also has the capability of overlaying (i.e., warping) the
29
 
vector data onto the raster image to visually check the
30
 
accuracy of the registration based on the current set of
31
 
active control points.  During the process of marking
32
 
points and entering map coordinates, the user can compute
33
 
the RMS (root mean square) error for each point entered.
34
 
The <EM>i.vpoints</EM>  program does this by calculating a
35
 
transformation equation (the same one that is calculated in
36
 
the GRASS program
37
 
 
38
 
<A HREF="i.rectify.html">i.rectify</A>
39
 
2).  Coefficients are
40
 
computed for the equation.  The coefficients are then used in the
41
 
equation along with the x,y coordinates of the marked points.  The
42
 
results are plugged into an equation for RMS error.  The
43
 
interpretation of RMS error is described in the ANALYZE subsection.
44
 
 
45
 
 
46
 
<P>
47
 
The procedures for marking control points (registration
48
 
points), displaying vector map layers,  overlaying vector maps onto
49
 
the raster image, and calculating RMS error are described in the
50
 
following sections .
51
 
 
52
 
 
53
 
<P>
54
 
To enter the program (the <EM>i.vpoints</EM> program requires the use of a
55
 
graphics monitor) type <EM>i.vpoints</EM>:
56
 
 
57
 
 
58
 
<DL>
59
 
     <DD><B>i.vpoints</B>
60
 
</DL>
61
 
 
62
 
The first prompt in the program asks for the imagery group to
63
 
be registered . . .
64
 
 
65
 
 
66
 
<PRE>
67
 
     Enter imagery group to be registered
68
 
     Enter 'list' for a list of existing imagery groups
69
 
     Enter 'list -f' for a verbose listing
70
 
     Hit RETURN to cancel request
71
 
     &gt;
72
 
</PRE>
73
 
 
74
 
For example,   
75
 
 
76
 
<DL>
77
 
<DD>&gt;<B>list</B>     
78
 
</DL>
79
 
 
80
 
might produce the following response:
81
 
 
82
 
<PRE>
83
 
     &lt;list&gt;
84
 
     Available groups
85
 
     - - - - - - - - - - - - - - - - - - - - - - - - - -
86
 
     test
87
 
     - - - - - - - - - - - - - - - - - - - - - - - - - -
88
 
     hit RETURN to continue --&gt;
89
 
</PRE>
90
 
whereas,  
91
 
<DL>
92
 
<DD>&gt;<B>list -f</B>       
93
 
</DL>
94
 
 
95
 
might produce the following response
96
 
 
97
 
<PRE>
98
 
     Available groups
99
 
     - - - - - - - - - - - - - - - - - - - - - - - - - -
100
 
     test
101
 
         photo in PERMANENT
102
 
         res2 in tifftest
103
 
     - - - - - - - - - - - - - - - - - - - - - - - - - -
104
 
     hit RETURN to continue --&gt;
105
 
</PRE>
106
 
 
107
 
The imagery group entered above should contain the files
108
 
that you wish to rectify.  After entering the group to be
109
 
registered, the terminal screen displays the message:
110
 
<PRE>
111
 
     &gt;test
112
 
     &lt;test&gt;
113
 
     Use mouse now . . .
114
 
</PRE>
115
 
And the color graphics monitor displays the following
116
 
screen:
117
 
 
118
 
<P>
119
 
     Any single file in the imagery group may be used to mark
120
 
points, and points can be marked on more than one file in the
121
 
imagery group to accumulate the suggested minimum number of points
122
 
(3 for a 1st order transformation, 6 for a 2nd order, and 10 for a
123
 
3rd order).  Any file in the imagery group can be subsequently
124
 
rectified (using <EM><A HREF="i.rectify.html">i.rectify</A></EM>)
125
 
based on the transformation matrix computed from these points.
126
 
 
127
 
<P>
128
 
     The chosen file is displayed in the upper left quadrant of the
129
 
monitor at a default magnification based on the extent of the
130
 
current active window.
131
 
 
132
 
<P>
133
 
 
134
 
<B>RASTER IMAGE</B>
135
 
 
136
 
<P>
137
 
     The raster image option on the menu at the bottom of the
138
 
window allows the user to display any single file in the imagery
139
 
group in the upper right quadrant of the window screen.  The option
140
 
provides the same file selection pick list as is presented when you
141
 
first enter the <EM>i.vpoints</EM> program.  When you select this option,
142
 
the program will erase the data contained in all of the four
143
 
quadrant windows and will reinitialize all program values.
144
 
 
145
 
 
146
 
<P>
147
 
 
148
 
<B>VECTOR MAPS</B>
149
 
<P>
150
 
     The vector maps option on the menu at the bottom of the screen
151
 
allows the user to display vector map data in the upper right
152
 
quadrant of the screen.  After selecting the vector map layer to
153
 
display, a menu selection bar appears along the bottom on the
154
 
screen. This pick list is used to select the line color (blue,
155
 
gray, green, red, white, or yellow) for the selected vector data
156
 
layer.
157
 
<P>
158
 
<B>Refresh</B>
159
 
<P>
160
 
     The refresh option on the main menu allows the user to
161
 
"refresh" or re-draw the displayed vector data. This function will
162
 
erase all outlines showing the limits of previously zoomed areas. 
163
 
A "yes/no" roompt will appear:
164
 
<PRE>
165
 
     Refresh Map ?       NO   YES
166
 
</PRE>
167
 
<P>
168
 
<B>Zoom</B>
169
 
<P>
170
 
     To enlarge a raster or vector image, place the mouse cross
171
 
hairs on the word zoom on the main menu and press the left button. 
172
 
The following menu will be displayed at the bottom. of the screen:
173
 
<PRE>
174
 
     CANCEL   BOX   POINT      Select type of zoom
175
 
</PRE>
176
 
You have the option to identify the map extent of the zoom window
177
 
using either the mouse to define a box, or the mouse to mark a
178
 
center point from which to enlarge the image.  The box option first
179
 
prompts you to identify a starting corner for the zoom region and
180
 
then allows you to define the area to be zoomed using a rubber band
181
 
box.  The prompts appear as follows...
182
 
<PRE>
183
 
     CANCEL    Mark the first corner of region
184
 
 
185
 
     CANCEL     Define the region
186
 
</PRE>
187
 
After marking the first corner of the region to be enlarged, hold
188
 
down the left button and move the mouse to change the size and
189
 
shape of the rubber band box.  After defining the area to be
190
 
enlarged, press the right button  to accept it.
191
 
 
192
 
 
193
 
<P>
194
 
 
195
 
The point method for enlarging an image will display a mouse menu
196
 
to guide you in selecting the appropriate enlargement.  To enlarge
197
 
or reduce the magnification factor, place the cursor on the "+" or
198
 
"-" box and press the left button on the mouse.
199
 
<P>
200
 
<P>
201
 
You may zoom either the raster or the vector display.
202
 
<P>
203
 
Upon accepting the new region limits, the raster or vector data are
204
 
redisplayed in either the lower left (raster) or lower right
205
 
(vector) windows.
206
 
<P>
207
 
The extent of the zoomed area is outlined on the unzoomed image in
208
 
the main window area.
209
 
<P>
210
 
     While the main menu is displayed,  you can mark
211
 
corrosponding control points on the raster and vector
212
 
images or enter map coordinates from the keyboard.  If you
213
 
are using coordinates taken from a reference map, circle
214
 
these points and then use whatever means you have available
215
 
to identify as precisely as possible the coordinate values
216
 
for these points.  Digitizing software is recommended,
217
 
especially GRASS v.digit.  Once you have
218
 
determined the standard coordinates (for example, UTM's) of
219
 
each circled point, you are ready to mark the points on the
220
 
displayed image.
221
 
<P>
222
 
     To mark the points on the image, that correspond to the points
223
 
on the standard coordinate map, place the mouse cross hairs on the
224
 
point on the image to be marked (you will probably have to ZOOM to
225
 
find the exact spot) and press the left hand button on the mouse. 
226
 
A diamond shaped symbol will be marked on the image.  The text
227
 
monitor will display the following screen:
228
 
<PRE>
229
 
          Point 1 marked on the image at
230
 
          East:  1023.77
231
 
          North:  -164.41
232
 
 
233
 
          Enter coordinates as east north:
234
 
</PRE>
235
 
<B>Analyze</B>
236
 
<P>
237
 
 
238
 
     After a number of points have been marked (a minimum of 4 for
239
 
a 1st order transformation, 7 for a 2nd order, and 11 for a 3rd
240
 
order), the RMS error of the points marked on the image can be
241
 
checked.  This is done by placing the cross hairs on the word
242
 
ANALYZE on the main menu at the bottom of the monitor.  The
243
 
following error report is superimposed on the monitor:
244
 
<PRE>
245
 
 
246
 
   error                      image                target
247
 
  #    col  row target   east      north       east      north 
248
 
     
249
 
  1   -0.9  0.0  1.0   1048.5     -144.8   679132.5  4351080.6
250
 
  2    1.0  0.4  1.3   2153.1     -567.2   684314.7  4399001.4
251
 
  .
252
 
  .
253
 
  .
254
 
          Overall rms error  76.85
255
 
</PRE>
256
 
 
257
 
     The RMS error for the image being rectified is recorded under
258
 
the column "error" and subTITLEd "row' and "col".  In the above
259
 
report, the marked point number 1 is 0.0 rows and -0.9 columns from
260
 
the predicted location calculated by the transformation equation.
261
 
     The RMS error for the target database map is recorded under
262
 
the heading "error"  and the subheading "target".  This is the RMS
263
 
error for the east and the north coordinate values of the target
264
 
map, but it is represented in the table using one general value.
265
 
     The overall RMS error for the image is displayed at the bottom
266
 
of the screen in meters.  Points that generate a high RMS error are
267
 
displayed in red on the monitor.
268
 
     The x,y coordinate values of the point marked on the image
269
 
being rectified are recorded under the heading "image" and the
270
 
subheadings "east" and "north".
271
 
     The standard coordinate values of the point in the target
272
 
database are recorded under the heading "target" and the
273
 
subheadings "east" and "north".
274
 
     If the user would like to exclude or include a point, this can
275
 
be accomplished by placing the mouse cross hairs on the point
276
 
number to be included (if the point is absent) or excluded (if the
277
 
point is displayed) and then pressing the left button on the mouse
278
 
twice.  When a point is excluded, it is not included in the
279
 
calculation of the RMS error, or included in the final
280
 
transformation matrix.  However, it can be retrieved within
281
 
<EM>i.vpoints</EM> at any time by double clicking with the mouse as
282
 
described above.
283
 
<P>
284
 
The following menu appears at the bottom of the monitor:
285
 
<PRE>
286
 
   DONE   PRINT   FILE   OVERLAY   DELETE ON       Transformation -
287
 
-&gt;   1st ORDER   Double click on point to be DELETED
288
 
</PRE>
289
 
Selecting DELETE ON will toggle the option to DELETE OFF, the
290
 
toggle option is used to allow the user to physically remove a
291
 
control point from the POINTS file instead of just flagging it as
292
 
an non-active reference point.
293
 
<P>
294
 
<B>Overlay</B>
295
 
<P>
296
 
Overlay allows the user to overlay the vector map(s) onto the
297
 
raster image.  Overlay can be used to warp (register) and display
298
 
the selected vector map data on top of the raster image contained
299
 
in the upper left window of the color screen.  An inverse
300
 
coordinate transformation is performed using the currently active
301
 
order of transformation (i.e., first, second, or third).
302
 
<PRE>
303
 
     Overlay vectors on raster image   NO   YES
304
 
</PRE>
305
 
By selecting the 1st ORDER option, the user may select the order of
306
 
transformation desired:
307
 
<PRE>
308
 
     Select order of transformation --&gt; 1st Order 2nd Order 3rd Order
309
 
</PRE>
310
 
The program will immediately recalculate the RMSE and the number of
311
 
points required.
312
 
<P>
313
 
<B>Quit</B>
314
 
<P>
315
 
     To exit the <EM>i.vpoints</EM> program, place the mouse cross hairs on
316
 
the word QUIT at the bottom of the monitor and all of the marked
317
 
points (including coordinates) will be saved.
318
 
 
319
 
<H2>SEE ALSO</H2>
320
 
 
321
 
The GRASS 4 <em>
322
 
<A HREF="http://grass.itc.it/gdp/imagery/grass4_image_processing.pdf">Image
323
 
Processing manual</A></EM>
324
 
 
325
 
<P>
326
 
<EM>
327
 
<A HREF="g.mapsets.html">g.mapsets</A><BR>
328
 
<A HREF="i.group.html">i.group</A><BR>
329
 
<A HREF="i.points.html">i.points</A><BR>
330
 
<A HREF="i.rectify.html">i.rectify</A><BR>
331
 
<A HREF="i.target.html">i.target</A>
332
 
</EM><br>
333
 
<em><a href="gm_georect.html">gis.m: GEORECTIFY TOOL</a></em>
334
 
 
335
 
 
336
 
<H2>AUTHOR</H2>
337
 
 
338
 
William R. Enslin, Michigan State University Center for Remote Sensing<br>
339
 
Radim Blazek (update to GRASS 5.7)<BR>
340
 
Hamish Bowman (finish GRASS 6 porting)
341
 
<p>
342
 
<i>Last changed: $Date: 2007-10-29 12:12:37 +0100 (Mon, 29 Oct 2007) $</i></p>