← Back to branch summary

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

« back to all changes in this revision

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

Tags: 7.0.0~rc1+ds1-1~exp1
http://bugs.debian.org/776349
* 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:
imagery/i.points/description.html
1
 
<H2>DESCRIPTION</H2>
2
 
 
3
 
<EM>i.points</EM> is an
4
 
<A HREF="imagery.html">imagery</A> function that enables
5
 
the user to mark points on a (raster) image to be rectified
6
 
and then input the geographic coordinates of each point for
7
 
calculation of a coordinate transformation matrix.
8
 
<EM>i.points</EM> must be followed by use of the GRASS
9
 
program <EM><A HREF="i.rectify.html">i.rectify</A></EM>, which
10
 
rectifies the image using the transformation matrix
11
 
coefficients calculated by <EM>i.points</EM>.
12
 
 
13
 
<P>
14
 
 
15
 
<EM>Rectification</EM> is the mapping (transformation) of
16
 
an image from one coordinate system to another.  The
17
 
geometry of an image extracted into a GRASS LOCATION having
18
 
an x,y coordinate system is not planimetric.  To create a
19
 
planimetric image, that is, to convert the x,y coordinate
20
 
system into a standard coordinate system (for example, the
21
 
UTM coordinate system or the State Plane coordinate
22
 
system), points from a map having the standard coordinates
23
 
must be associated with the same points on the image to be
24
 
rectified.  <EM>i.points</EM> enables the user to mark
25
 
points on an image and input the standard coordinates for
26
 
that point.  <EM>i.points</EM> then calculates a least
27
 
squares regression using the two coordinate systems (x,y
28
 
and standard) for the marked points.  A matrix containing
29
 
transformation coefficients is the output file for
30
 
<EM>i.points</EM>.
31
 
 
32
 
<P>
33
 
During the process of marking points and entering map
34
 
coordinates, the user can compute the RMS (root mean
35
 
square) error for each point entered.  <EM>i.points</EM>
36
 
does this by calculating the transformation equation (the
37
 
same one that is calculated in the GRASS program <EM>
38
 
<A HREF="i.rectify.html">i.rectify</A></EM>), and then
39
 
plugging these results into an equation for RMS error.
40
 
 
41
 
<P>
42
 
<EM>i.points</EM> offers a zoom option to locate precisely
43
 
the point to be marked on an image.  This program also
44
 
offers the user the option of acquiring standard
45
 
coordinates for a marked point from a map layer in the
46
 
target data base.
47
 
 
48
 
<P>
49
 
<EM><A HREF="i.target.html">i.target</A></EM> must be
50
 
run before running <EM>i.points</EM> to enable the PLOT
51
 
RASTER option to be used and to identify a target data base
52
 
LOCATION_NAME and MAPSET for the rectified image.  To run
53
 
<EM>i.points</EM>, a graphics monitor is required.
54
 
 
55
 
<P>
56
 
The procedure for marking points, entering coordinates, and
57
 
calculating RMS error is described below.
58
 
 
59
 
<P>
60
 
The first prompt in the program asks the user for the 
61
 
<A HREF="i.group.html">imagery group</A> to be
62
 
registered.  Note that if 
63
 
<EM><A HREF="i.target.html">i.target</A></EM> is not run
64
 
before <EM>i.points</EM>, the <EM>i.points</EM> program
65
 
will display the following error message:
66
 
 
67
 
<PRE>
68
 
ERROR: Target information for group [spot] missing
69
 
       Please run i.target for group [spot]
70
 
</PRE>
71
 
 
72
 
After entering the group to be registered the terminal screen
73
 
displays the message: 
74
 
 
75
 
<PRE>
76
 
      Use mouse now...
77
 
</PRE>
78
 
 
79
 
The graphics monitor displays the following screen:
80
 
 
81
 
<PRE>
82
 
 ______________________________________________________________
83
 
| imagery  filename  (mag)        |  target  filename  (mag)  |
84
 
|_________________________________|___________________________|
85
 
|                                 |                           |
86
 
|                                 |                           |
87
 
|                                 |                           |
88
 
|                                 |                           |
89
 
|                                 |                           |
90
 
|                                 |                           |
91
 
|                                 |                           |
92
 
|_________________________________|___________________________|
93
 
|                                 |                           |
94
 
|                                 |                           |
95
 
|                                 |                           |
96
 
|                                 |                           |
97
 
|                                 |                           |
98
 
|                                 |                           |
99
 
|                                 |                           |
100
 
|_________________________________|___________________________|
101
 
| QUIT  ZOOM  PLOT RASTER  ANALYZE|                           |
102
 
|_________________________________|___________________________|
103
 
</PRE>
104
 
 
105
 
A pop-down menu like that shown below will be superimposed on the
106
 
left half of the screen:
107
 
 
108
 
 
109
 
<PRE>
110
 
        __________________________________
111
 
       | Double click on raster map layer|
112
 
       |          to be plotted          |
113
 
       |   Double click here to cancel   |
114
 
       |_________________________________|
115
 
 
116
 
               ____________________
117
 
              |    Mapset demo    |
118
 
              |___________________|
119
 
              | spotclass|  spot.1|
120
 
              |__________|________|
121
 
              | composite|  spot.2|
122
 
              |__________|________|
123
 
              | spot.3   |        |
124
 
              |__________|________|
125
 
</PRE>
126
 
 
127
 
Any single raster map layer in the 
128
 
 
129
 
<A HREF="i.group.html">imagery group</A> may be used on
130
 
 
131
 
which to mark points, and the user can mark points on more
132
 
than one raster map layer in the 
133
 
 
134
 
<A HREF="i.group.html">imagery group</A> to accumulate the
135
 
 
136
 
suggested minimum number of 12 points.  Any raster map
137
 
layer in the <A HREF="i.group.html">imagery group</A>
138
 
can be rectified (using 
139
 
 
140
 
<EM><A HREF="i.rectify.html">i.rectify</A></EM>) based on the
141
 
transformation matrix computed from these points.
142
 
 
143
 
 
144
 
<P>
145
 
 
146
 
The imagery file chosen by
147
 
the user is displayed in the upper left quadrant of the
148
 
screen.
149
 
 
150
 
<H3>ZOOM</H3>
151
 
 
152
 
To magnify the displayed file, the user must
153
 
place the mouse cross hairs on the word ZOOM.  The following menu
154
 
will then be displayed at the bottom of the screen:
155
 
 
156
 
<PRE>
157
 
   ____________________________________________
158
 
  | Cancel|  Box|  Point|  Select type of ZOOM|
159
 
  |_______|_____|_______|_____________________|
160
 
</PRE>
161
 
 
162
 
 
163
 
The user has the option of identifying the zoom region
164
 
either by using the mouse to make a box, or by using the
165
 
mouse to mark the two diagonal points of the desired
166
 
region.  The terminal screen will display a mouse button
167
 
menu to guide the user in identifying the corner points of
168
 
the region.
169
 
 
170
 
<H3>MARKING POINTS</H3>
171
 
 
172
 
To mark the points on the image that correspond to the
173
 
points on a standard coordinate system map, the user must
174
 
place the mouse cross hairs on the corresponding location
175
 
on the image to be marked and press the left hand button on
176
 
the mouse.  A diamond shaped symbol will be marked on the
177
 
image.  The user's terminal will display the following
178
 
menu:
179
 
 
180
 
 
181
 
<PRE>
182
 
 ___________________________________________________________
183
 
| Point 1 marked on the image at|                          |
184
 
| East: 1023.77                 |                          |
185
 
| North: -164.41                |                          |
186
 
|                               |                          |
187
 
|                               |                          |
188
 
|                               |                          |
189
 
|                               |                          |
190
 
|                               |                          |
191
 
|__________________________________________________________|
192
 
| Enter coordinates as east north:                         |
193
 
|__________________________________________________________|
194
 
</PRE>
195
 
 
196
 
The user then enters the easting and northing (separated by
197
 
a space) for the point marked on the image.  If the user
198
 
wishes not to enter a coordinate, he or she may simply hit
199
 
RETURN to return control to the mouse;  the marked point
200
 
then disappears.
201
 
 
202
 
<H3>PLOT RASTER</H3>
203
 
 
204
 
In addition to acquiring reference points from a standard
205
 
map, the user has the option of acquiring the reference
206
 
points from a raster map layer in the target data base
207
 
LOCATION_NAME.  The data base raster map layer is displayed
208
 
by placing the mouse cross hairs on the words PLOT RASTER.
209
 
The following line is then displayed at the bottom of the
210
 
graphics monitor:
211
 
 
212
 
 
213
 
<PRE>
214
 
       ________________________________________________
215
 
      | Cancel|  Indicate which side should be plotted|
216
 
      |_______|_______________________________________|
217
 
</PRE>
218
 
 
219
 
Which side of the graphics monitor is to be plotted
220
 
is indicated by placing the mouse cross hairs on the 
221
 
half of the graphics monitor screen that the user would like to use,
222
 
and pressing the left mouse button.  The following pop-down
223
 
menu will be superimposed on the half of the screen that was chosen:
224
 
 
225
 
 
226
 
<PRE>
227
 
    _________________________________________
228
 
   | Double click on raster (cell) map layer|
229
 
   |              to be plotted             |
230
 
   |       Double click here to cancel      |
231
 
   |________________________________________|
232
 
 
233
 
          ______________________________
234
 
         |         Mapset demo         |
235
 
         |_____________________________|
236
 
         | tm.rectified |              |
237
 
         |______________|______________|
238
 
         | tm.classified|              |
239
 
         |_____________________________|
240
 
         |       Mapset PERMANENT      |
241
 
         |_____________________________|
242
 
         | elevation    |  geology     |
243
 
         |______________|______________|
244
 
         | slope        |  soils       |
245
 
         |______________|______________|
246
 
         | aspect       |              |
247
 
         |______________|______________|
248
 
         | roads        |              |
249
 
         |______________|______________|
250
 
</PRE>
251
 
 
252
 
After the raster map layer is displayed the following message appears
253
 
at the bottom of the graphics monitor:
254
 
 
255
 
 
256
 
<PRE>
257
 
     ______________________________________
258
 
     | input method --&gt;|  keyboard|  screen|
259
 
     |_________________|__________|________|
260
 
</PRE>
261
 
 
262
 
If the user wishes to use the plotted raster map layer only
263
 
as a comparative reference, then the keyboard can be chosen
264
 
as the means to input coordinates corresponding to the
265
 
marked points on the image.  This is done by placing the
266
 
mouse cross hairs on the word KEYBOARD and pressing the
267
 
left button on the mouse.
268
 
 
269
 
 
270
 
<P>
271
 
 
272
 
If the user selects the SCREEN option, then points marked
273
 
on the image will automatically be associated with the
274
 
coordinates from the corresponding points on the target
275
 
data base map layer.  In this option, when the user marks a
276
 
point on the image, the following menu is displayed at the
277
 
terminal:
278
 
 
279
 
 
280
 
<PRE>
281
 
___________________________________________________________
282
 
| Point 5 marked on the image at|                          |
283
 
| East: 1023.77                 |                          |
284
 
| North: -164.41                |                          |
285
 
|                               |                          |
286
 
|                               |                          |
287
 
| Point located at              |                          |
288
 
| East: 679132.57               |                          |
289
 
| North: 4351080.67             |                          |
290
 
|                               |                          |
291
 
|                               |                          |
292
 
|                               |                          |
293
 
|                               |                          |
294
 
|                               |                          |
295
 
|__________________________________________________________|
296
 
| use mouse now...                                         |
297
 
|__________________________________________________________|
298
 
</PRE>
299
 
 
300
 
 
301
 
The user then uses the mouse to mark a corresponding point
302
 
on the displayed image from the target data base.  The
303
 
coordinates for the target data base map layer are
304
 
automatically saved as the coordinates corresponding to the
305
 
marked point on the image.
306
 
 
307
 
<H3>ANALYZE</H3>
308
 
 
309
 
After a number of points have been marked (4 to 7), the
310
 
user can check the RMS error of the points marked on the
311
 
image.  This is done by placing the mouse cross hairs on
312
 
the word ANALYZE at the bottom of the graphics monitor.  An
313
 
error report resembling that shown below is superimposed on
314
 
the graphics monitor:
315
 
 
316
 
 
317
 
<PRE>
318
 
______________________________________________________________________
319
 
|         error                 image                 target         |
320
 
|#     row     col    target    east    north      east       north  |
321
 
|____________________________________________________________________|
322
 
|1     0.0     -0.9    1.0     1048.5   -144.8   679132.5   4351080.6|
323
 
|2     0.4     1.0     1.3     2153.1   -567.2   684314.7   4399001.4|
324
 
|3    -1.2     -0.5     .6     1452.8   -476.5   567841.4   3457682.8|
325
 
|4     1.1     0.5     1.3     1034.0   -109.2   677573.8   4352626.4|
326
 
|5    -2.7     14.0    14.2    1048.6   -144.9   679132.6   4351080.7|
327
 
|                                                                    |
328
 
|____________________________________________________________________|
329
 
|    overall   rms    error:   4.46                                  |
330
 
|____________________________________________________________________|
331
 
</PRE>
332
 
 
333
 
The following menu then appears at the bottom of the graphics monitor:
334
 
 
335
 
<PRE>
336
 
 ____________________________________________________________________
337
 
| DONE|  PRINT FILE|   Double click on point to be included/excluded|
338
 
|_____|____________|________________________________________________|
339
 
</PRE>
340
 
 
341
 
The RMS error for the image is given under the column
342
 
TITLEd "error" and subTITLEd "row" and "col".  In the above
343
 
report, point number 1 is 0.0 rows and -0.9 columns from
344
 
the predicted location calculated from the transformation
345
 
equation.  The RMS error for the target raster map layer is
346
 
listed under the heading "target".  This is the RMS error
347
 
for the east and north coordinates of the target map layer,
348
 
but it is presented in the table using one general value.
349
 
The overall RMS error is displayed at the bottom of the
350
 
screen in meters.  Points that create high RMS error are
351
 
displayed in red on the graphics monitor (represented here
352
 
in italics).
353
 
 
354
 
<P>
355
 
The location of the point marked on the 
356
 
<A HREF="i.group.html">imagery group</A> file is given
357
 
under the heading "image" and the subheadings "east" and
358
 
"north".  The location of the point in the target data base
359
 
is given under the heading "target" and the subheadings
360
 
"east" and "north".  If the user would like to exclude or
361
 
include a point, this can be accomplished by placing the
362
 
mouse cross hairs on the point number to be included (if
363
 
the point is absent) or excluded (if the point is
364
 
displayed) and pressing the left button on the mouse
365
 
twice.  When a point is excluded, it is not afterwards
366
 
included in the calculation of the RMS error, or included
367
 
in the final transformation matrix.  However, it can be
368
 
retrieved within <EM>i.points</EM> at any time by double
369
 
clicking with the mouse as described above.
370
 
 
371
 
<H3>QUIT</H3>
372
 
 
373
 
To end the <EM>i.points</EM> program place the mouse cross
374
 
hairs on the word QUIT;  the marked points (including
375
 
coordinates) will be saved.
376
 
 
377
 
 
378
 
<H2>NOTES</H2>
379
 
 
380
 
A good rule of thumb is to mark at least 12 to 15 points
381
 
which are evenly distributed over the entire 
382
 
<A HREF="i.group.html">imagery group</A> file in order to
383
 
obtain an accurate transformation equation for the
384
 
rectification process.  The RMS error may increase with
385
 
more points added, but the transformation equation will be
386
 
more accurate.
387
 
 
388
 
<P>
389
 
An RMS error of less than or equal to approximately one
390
 
resolution unit (pixel or cell) for the image being
391
 
rectified is generally considered acceptable.
392
 
 
393
 
<P>
394
 
In order to use a digitizer with <EM>i.points</EM>, at
395
 
least one digitizer driver besides "none" (the on-screen
396
 
digitizer) must be available in the digitcap file.
397
 
 
398
 
<P>
399
 
This program is interactive.
400
 
 
401
 
 
402
 
<H2>SEE ALSO</H2>
403
 
 
404
 
The GRASS 4 <em>
405
 
<A HREF="http://grass.itc.it/gdp/imagery/grass4_image_processing.pdf">Image
406
 
Processing manual</A></EM>
407
 
 
408
 
<P>
409
 
<EM>
410
 
<A HREF="g.mapsets.html">g.mapsets</A><br>
411
 
<A HREF="i.group.html">i.group</A><br>
412
 
<A HREF="i.rectify.html">i.rectify</A><br>
413
 
<A HREF="i.target.html">i.target</A><br>
414
 
<A HREF="i.vpoints.html">i.vpoints</A>
415
 
</EM><br>
416
 
<em><a href="gm_georect.html">gis.m: GEORECTIFY TOOL</a></em>
417
 
 
418
 
 
419
 
<H2>AUTHOR</H2>
420
 
 
421
 
Michael Shapiro,
422
 
U.S.Army Construction Engineering Research Laboratory
423
 
<p>
424
 
<i>Last changed: $Date: 2007-10-29 12:12:37 +0100 (Mon, 29 Oct 2007) $</i></p>