GRIDVECTOR (version 1) : AN ARC/INFO AML PROGRAM TO EXTRACT
LINEAR FEATURES FROM A GRAY-SCALE IMAGE OF A PAPER
GEOLOGICAL MAP

by

Feliks Persits

 

Open-File Report 97-713

This report is preliminary and has not been reviewed for conformity with U. S. Geological Survey editorial standards.Any use of trade, product, or firm names is for descriptive purposes only and doesnot imply endorsement by the U.S. Government.


Although this program has been used by the U. S. Geological Survey, no warranty,expressed or implied, is made by the USGS as to the accuracy and functioning of the program and related program material, nor shall the fact of distribution constitute any such warranty, and no responsibility is assumed by the USGS in connection therewith.

Figures
Source Code

Previous Next

GRIDVECTOR, an ARC/INFO macro language program, provides a more effective method of accomplishing the first step of on-screen digitization of large and complicated geologic maps. The traditional method used is first, hands-on tracing of selected lines on a transparent mylar sheet. The mylar is then scanned in black and white mode and the scanned image converted to lines by the ARC/INFO GRIDLINE or ARCSCAN commands.

Disadvantages of the traditional method are:

  1. Much more time is needed to prepare drawings on mylar when features on the original geologic map are very tiny and the number of polygons is large (e.g., hundreds or thousands).
  2. Poor accuracy is common because of accidental shifts of the mylar while tracing lines and because lines may be inadvertently omitted.
  3. Adding these omitted features requires rescanning the mylar and consequently introducing new spatial errors when matching the new features to the existing data.
  4. Identifying the same line on the screen and on the original paper map is time consuming and tedious work.

 

GRIDVECTOR is designed to trace boundaries of geologic features from a gray-scale scanned image of a paper geologic map. The advantages of this method are:

  1. A gray-scale scanned image is far less expensive than a color image and the image itself needs much less disk space.
  2. Even at relatively low resolution (200 dots/inch), most geological features are retained in the gray-scale image.
  3. The color separates of the original color map, if available, can be used as input to GRIDVECTOR to trace colored boundaries.
  4. Adding newly requested features or correcting existing ones is simple because the same gray scale image is used.

Scanned images are transformed into a grid, using the ARC/INFO IMAGEGRID command and georeferenced using the ARC/INFO REGISTER and RECTIFY commands. The georeferencing is done mainly by using latitude/longitude crosses or other features with known coordinates as references.

The output of GRIDVECTOR is a line and polygon ARC/INFO coverage. Here lines are supposed to be one of the following features on the gray-scale image of the original paper map:

a) the centerlines of the dark, relatively thin, strips (polygon outlines, faults, etc. on the paper map);

b) boundaries between areas with different gray levels (different color or pattern on the paper map).

The method presented in this report extracts lines from a gray-scale image. These lines (strips of pixels) can be easily seen on a computer screen in ARCEDIT, when the image has been used as the background. If the original map is of good quality (i.e., good line contrast and a small number of different patterns and text) the GRIDVECTOR program can extract 50-60% of the lines, saving a great amount of handwork and time. The rest of the digitization, and correction, can be done effectively along with the geologic attributing in ARCEDIT. This helps avoid errors caused by separation of these processes in time.

The program is based on the assumption that a dark (mainly black) line outlines every polygon (geologic feature) on a paper map. Such a line in gray-scale is represented by a thin strip of pixels with minimum gray values. Features that can be treated as noise are other dark lines, such as rivers, and text annotations. To reduce the number of lines and polygons created because of such noise, the average slope of gray-scale pseudo relief was calculated. The average gray levels are about the same on both sides of the dark line representing a river, text, or coordinate grid. By contrast, a geologic polygon gives a jump in gray level because of the different color or pattern on either side of the line. The user has a choice of selecting lines by the minimum of gray values, by average slope, or by both methods. Before calculation of an average slope, several consecutive smoothing filters are applied.

The program is based entirely on standard ARC/INFO, ARCSCAN and ARCEDIT commands. ARC/INFO’s GRID program is limited to a maximum filter size of 7 cells (pixels). Therefore, the optimal image resolution is about 200 dots/inch, giving a width of 4-6 pixels for dark lines.

The program uses ARC/INFO’s GRIDLINE command to create a line coverage; the final results are very dependent on the effectiveness of GRIDLINE.

An easy way to apply the GRIDVECTOR AML is to select parameters on a relatively small portion of the input grid. The user can then experiment as many times as necessary, because the program deletes all previously made files. The goal is to select parameters that allow an optimum output that satisfactorily filters out the noise without removing too many geological boundaries.


Previous Next

U.S. Geological Survey Open-File Report 97-713


USA.gov logo  Take Pride in America button