GIS Tutorial2
GIS Tutorial2
E00 is one of the most popular formats in GIS worldwide. The majority of GIS
software has a function to input and convert the *.E00 file format. Particularly when
you download free GIS data via the Internet, it would be very likely that you would
obtain an E00 file. In this section, you will become familiar with how to convert a
*.E00 file to a shape file.
This data was freely downloaded from one of the most well-known GIS websites,
GeoComunity. From GeoComunity, you can download country-wide data such as
administrative boundaries, rivers, roads and others, for any countries throughout the
world. The address of the website is as follows:
http://www.geocomm.com/
2. Start ArcMap > open a new map > set the coordinate system of the data frame to
WGS_1984 >click the Add Data button > browse to the same folder you have
decompressed the
E00 file > find the
decompressed E00
file > click the Add
button.
3-1
You will find that compressed E00 files will also appear in the Add Data window
to show that they are recognized by ArcGIS, but if you try to add them before
decompression you will encounter an error.
- poly: polygon feature (in this data set, shows lakes and very wide part of rivers)
You can use any of these together with other data to create your map. For example,
you can add the feature of Eth to create a river distribution map in Ethiopia.
4. Add the feature of Eth which you created in Assignment 3. Set the symbology as no
color for fill and red for outline. Then you will find the river route stretch and lakes
do not match the boundary of Ethiopia well; in the northern part are rivers and lakes
which lie outside of Ethiopia.
5. To make the river route stretch and lakes match the country’s boundary, open
ArcToolbox > Analysis Tools > Extract > Clip. Then, using the dropdown menu,
set the Input Features to the layer DNNET_arc, and Eth for Clip Features. Next,
set the Output location to your Result folder with the Output name as eth_river.
3-2
6. Using the same procedure as in Step 5, you can clip the DNNET_poly by Eth to
obtain a polygon feature of Eth_lake.
Note: the point layer is useless in this exercise and the text layer can not be
converted to a shape file.
7. If you would like to perform any editing, you may, for example, express all river
routes as only one line no matter how big the river width is. To do this, Start
Editing by clicking the Edit Tool button then click Start Editing in the
dropdown menu for Editor.
Notice at this moment that it is not possible to edit all layers in your data frame,
only Eth and those clipped features can be selected from the Target layer list.
Therefore, after you add E00 data into your data frame, it is recommended to export
the necessary layers to your Result folder so that these data will be more
conveniently used and also be editable.
GIS software “Global Mapper” (USGS) has a very powerful format conversion
function. Not only *.E00 files, but many other common GIS data formats can also be
read and then converted to other formats to meet your needs. If you convert a *.E00 file
with ArcGIS, you need to decompress the source data first because, generally, *.E00
data is distributed in a compressed format, especially when distributed over the Internet.
However, if you use Global Mapper to perform the conversion, it’s unnecessary to
decompress the source data, which allows you to keep your database light and simplify
the tasks you must perform.
You can check the File Type list to find how many formats are supported, and
consequently, all supported formats can be converted by Global Mapper.
3-3
2. In the “Open File(s)” window, select the file DNNET.E00.gz from the following
folder:
C:/EthiopiaArcGIS/Data/CountriesData/Ethiopia.
When you use Global Mapper for format conversion you don’t need to decompress
the original file.
3. If the file you select to open does not have a coordinate system indicator, you will
be requested to select projection/datum. At this point the only option is to simply
click OK button, so do so.
- Select Geographic
(Latitude/Longitude) in the
textbox for Projection
- Make sure the box is checked for “Use Selected Projection for All Selected
Files”
- Click OK
5. All the contents within the compressed data file DNNET.E00.gz will appear in the
window of Global Mapper.
3-4
6. In order to convert the *.E00 file to a commonly used format, i.e. shape file, you
can use the export function of Global Mapper. Click the File menu > click Export
Vector Data > choose the export type by clicking Export Shapefile.
7. In the Shapefile Export Options window, check the box for Export Areas and
then in the following Save As window, set the necessary parameters as:
- Location: C:/GIS_Exercise/Database/My_data
3-5
The contents of E00 files are as follows:
RDLINE : roads
RRLINE : railway
2. Create a background of your country by using the country boundary file that you
created in Assignment 3.
You are free to set the color for fill and outline as you like. However, please avoid
the setting of “no color”. Also, it is suggested that a dark color be set for the outline
and a light color for fill.
Set the width of the outline as 2.
3. Add layers above the background and set their symbology as the following
sequence:
- road (top layer, dark brown, width 1)
- lake (forth layer, dark blue for fill, the same dark or no color for outline
and the width of the outline not over 1)
3-6
3-2: Coordinate System Conversion
Exercise 1: Understand the Difference of Coordinate Systems
The survey area of the EWTEC project can be expressed by different coordinate
systems according to different purposes. The two kinds of basic coordinate systems are
Geographic (Longitude/Latitude) and UTM (Universal Transverse Mercator). Within
this exercise, you will learn the basic differences between them.
1. Start ArcMap, open a new map, save it as Prj_check in the Maps folder
When you add features into the map, a dialog box will appear to inform you if the
feature is missing coordinate system (Spatial reference) information. At this point
just click the OK button and continue with the exercise.
5. Select the layer SA_prj1, then zoom in to it (right-click the layer and then click
Zoom to Layer). Check the coordinate values in task bar below. You will find
nothing is strange. The coordinates of the study area appear to be properly
expressed in longitude and latitude.
6. Select the layer SA_UTM and also zoom in to it. This time you will find the
coordinate system shows thoroughly strange values. Obviously the values shown in
the task bar are anything but longitude and latitude.
3-7
Actually the unit in UTM projection is "meter". This unit is necessary in case you
would like to get an accurate calculation in length of lines and areas of polygons.
7. Click the Full Extent tool try to see both layers at the same time. You will find
nothing could be seen in the Data View. That is, you can not express the two layers
properly within the same Data View without processing the necessary coordinate
systems.
2. Zoom To Layer SA_prj1, then add all features named SA_prj# from the GeoDB:
C:/GIS_Exercise/Database/OtherData.mdb
3. Select the fill color for all the layers as No Color, and then set a unique outline
color for each layer so that you can easily distinguish between them. All the layers
within the Data View will be seen as shown in the left side of the figure below.
All layers give a similar expression of the study area of EWTEC, but each layer is a
little different in spatial extent from the other layers.
4. To see clearly the difference between layers, you can zoom in to any part in the
Data View, e.g. upper right corner, as shown in the right side of the figure below.
3-8
6. Open the Data Frame Properties window, if it was closed after the coordinate
system setting >
select the
General tab >
set the Units of
Display as
Meters. Then
click the OK
button.
7. Click the Measure tool > measure the difference between layers. You should
find the difference for the same point in different layers will be in the range of
about 50 to 1000 m.
1. Continue Exercise 1. Click the Insert Menu > click Data Frame. A new data frame
is added into the table of contents with the name New Data Frame. At the same
time, you will find that all layers have disappeared from the Data View.
The default setting dictate that when new data is added, the new data frame will be
the active data frame. Thus, from the Data View, you can only see the layers in the
new data frame.
3. Go back to the old data frame and rename that frame as Geographic.
4. Select the layer SA_UTM from data frame of Geographic > right-click the layer >
click Copy.
5. Right the data frame of UTM > click Paste Layer(s). The layer SA_UTM will be
pasted into the new data frame and is shown in the Data View.
6. Remove the SA_UTM layer from the old data frame (now titled Geographic).
3-9
C:/EthiopiaArcGIS/Data/OtherData/Prj_E
The two layers in the Data View are expected to have exactly the same shapes and
spatial extents.
8. Because the two layers have the same name, it is recommended that you change the
expression of the table of contents to Source View to better be able to distinguish
between them.
You can alternatively change the name of one or both layers to make them easily
distinguishable.
9. Open the attribute table of the original data SA_UTM, which is from the folder
C:/EthiopiaArcGIS/Data/OtherData/Prj_E. You will find there are three (3)
fields: FID, Shape, and NAME.
You will also find an asterisk attached to the Shape field, this asterisk indicates that
the value in this field can not be edited.
10. Open the attribute table of the other layer that is from your GeoDB OtherData.
You will find that five (5) fields are stored there. The first three fields correspond to
the three fields of the original data. The last 2 fields, Shape_length and
Shape_Area, give the length and area in the coordinate system unit and are
automatically added into the attribute table when you import the feature into the
GeoDB.
Therefore, if you find the default method of area calculation is troublesome, you
can simply import the feature into a GeoDB, after you perform the appropriate
coordinate system conversion.
11. You may decide the length of the polygon is not so important, but the area would be
very important data in many cases for analysis, so a reliable method is necessary for
area calculation.
In this exercise you can see from the attribute data that the area of the study area of
EWTEC is 6,759,590,300m2.
12. To calculate area by the default method of ArcGIS, you must follow the steps
below:
3-10
- Open the attribute table of a polygon feature layer for which you would
like to calculate the area. In this case it is SA_UTM from
C:/EthiopiaArcGIS/Data/OtherData/Prj_E.
- Click the Options button > click Add Field.... The Add Field dialog box
will appear.
- Set the Name of the new field as Area, Type as Float, Precision as 20 and
Scale as 2, then click the OK button.
- Right-click the title of the new field Area > click Calculate Values... and
click Yes in the following Field Calculator message box. The Field
Calculator window will appear
dblArea
- The calculation result is shown in the attribute table and the value is also
6,759,590,300m2.
If you find it difficult to remember the above statements for area calculation,
you can access them through the following procedure:
• Click ArcGIS Desktop Help > click the Key Word tab > input area
as the key word > select calculating for polygons > click the Display
button > scroll to the section of Updating area for a Shapefile
3-11
• You will find all the statements there. You can use the help menu as a
reference to input the statements yourself, or you can even select the
statements > right-click the selection area > click Copy > go back to
the Field Calculator window > put your curser in the corresponding
text box > right-click > click Paste.
Now go back to the table of contents and make sure it is in expression of Display.
Then go to the data frame of Geographic.
13. Select one of the layers with the geographic coordinate system e.g. SA_prj1. Open
its attribute table. You can find the two fields Shape_length and Shape_Area have
also been added into the attribute table when you imported them into GeoDB.
However, the values of both length and area are nonsense due to the fact that the
calculation is based on the unit of the coordinate system, that is, the decimal
expression of longitude and latitude.
14. You can try to use the procedure given in Step 12 to calculate the area once again,
and you will find the result is identical to that shown in the current attribute table.
In performing the above exercise, you are requested to keep in mind the following
very essential concepts of coordinate system:
- You can not make GIS data match each other by simply adding them into
the same data frame without performing the necessary coordinate system
processing.
3-12
- Coordinate systems can be generally classified into two basic types:
Geographic and Projected.
- The projected coordinate system generally uses meters as map unit. Also,
the most widely used projected coordinate system is UTM with differing
zones for countries or even areas.
When you obtain GIS data, you have to check the projection (coordinate system) by
which the data has been created. In most cases the projection data is stored in metadata,
that is, the information about the GIS data.
In the above exercises (Exercise 1 and 2), 7 features were used to express the same
location of the EWTEC study area. These 7 features are created in different projections
as follows:
3-13
- SA_UTM.shp : UTM Zone37N (world wide)
3. Right-click the feature class of SA_UTM, and then click Properties.... The
Feature Class Properties window will appear.
4. Click the Fields tab > click the Data Type of Geometry corresponding to the Field
Name of Shape. The Field Properties will appear in the corresponding list.
5. You will find the property is Unknown from the bottom row of Spatial Reference.
6. Because the coordinate system of the feature SA_UTM is UTM WGS1984 zone
37N, you can set the coordinate system by clicking the button next to the
property item. The Spatial Reference Property window will appear. The
Coordinate System tab will be selected as the default setting.
7. Now the Name of the coordinate system is Unknown, so you have to set it by
clicking one of the buttons in the lower part of the window. Because most
3-14
coordinate systems have been saved in the coordinate system list of ArcGIS, the
most convenient method to set the coordinate system is by clicking the Select...
button.
8. The Browse for Coordinate System window will appear, which you have
probably encountered several times in previous exercises. Expand the folder of
Projected Coordinate System > expand the folder of Utm > expand the folder of
Wgs 1984. The UTM WGS1984 zone list will appear.
9. Slide the scroll button to find the destination zone WGS 1984 UTM Zone 37N.prj
and then click the Add button.
10. When you go back to the Fields tab of the Feature Class Properties window you
can check again if you have correctly selected the coordinate system from the Field
Properties partition, and then click the OK button.
11. When you go back to the Coordinate System tab of the Spatial Reference
Property window you still have a chance to check if you have correctly selected
the coordinate system from the Name textbox, and then click the OK button.
12. Open ArcMap and open the map of prj_check that you created in the previous
exercise.
13. If the data frame of Geographic is inactive, make it active (so it is current showing
in Data View) by right-clicking on it and then clicking Activate.
14. Make all layers invisible except the layer SA_prj1, because you already know the
feature is created in the same coordinate system of WGS1984 that you set for the
dataset in Exercise1.
3-15
15. Adjust the symbology of the SA_prj1 layer to no color for fill color and red for
outline color.
17. This time you will find the layer AS_UTM can be shown in the same Data View
and exactly matches the layer of AS_prj1.
In this exercise you will have learned that if your GIS data are created in widely used
coordinate systems such as WGS84 and WGS84 UTM, and you have set the
coordinate system of those GIS data and the data frame, then the data may be
automatically matched to each other in the coordinate system of the data frame.
3. Expand the menu of Data Management Tools > expand Projections And
Transformations > open the Define Projection tool. The Define Projection
window will appear.
4. Select the SA_prj1 layer into the textbox of Input Dataset or Feature Class
5. Click the coordinate system selection button next to the Coordinate System
textbox , and the same Spatial Reference Properties window will appear.
6. Refer to the procedure you performed in the previous exercise, select Geographic
Coordinate System > World > WGS 1984.
3-16
7. Activate the data frame UTM.
9. Adjust the symbology of the layer SA_UTM as no color for Fill color and red for
Outline.
12. This time you will find the layer AS_prj1 can be shown in the same Data View
and exactly matches the layer of AS_UTM in UTM projection and map unit of
meters.
13. Right-click the layer of AS_prj1 > click Data > click Export Data. The Export
Data dialog box will appear.
3-17
15. When the exported feature Temp is added into the same data frame UTM, you can
check its coordinate system by opening the Properties > click the Source tab. You
will find the coordinate system has been changed to
WGS_1984_UTM_Zone_37N.
In this exercise you will have learned that ArcGIS provides several methods for
coordinate systems conversion and one of them is to add a feature to a data frame
with a different data frame, then export the feature with the coordinate system of
the data frame. But you will also learn in later exercises that this method is not
always effective. For many coordinate systems which are not as common as
WGS1984 or UTM, this method is often ineffective.
2. After you select the folder Geographic Coordinate System > double click the
Africa folder. The list of projections generally used in Africa will appear.
3. Since you already know the feature SA_prj2 is created in the coordinate system of
Abidjan 1987, you can find the corresponding prj file from the list, which will be
the first one. Select the Abidjan 1987.prj and then click the Add button.
5. Make all layers invisible except the layers SA_prj1 and SA_prj2 since you have
3-18
finished defining the coordinate systems for these two features. However, you will
find that the two layers still do not match each other even after their coordinate
system has been defined. In order to make SA_prj2 (in the coordinate system of
Abidjan 1987.prj) match SA_prj1 (in the coordinate system of WGS84) you have
to perform a project conversion process.
In this exercise you will practice converting the coordinate system of SA_prj2 from
Abidjan 1987.prj to WGS84.
7. Expand the menu for Data Management Tools > expand Projections And
Transformations > expand Feature > select and open the Project tool. The
Project tool window will appear.
8. Select the layer SA_prj2 in the textbox for Input Dataset or Feature Class
If you would like to practices more to master coordinate systems processing, including
define and conversion, you can practice by converting the other 4 features referring to
the coordinate system information given in the beginning of Exercise 3.
When you perform the coordinate system conversion, you might become confused how
to select the parameter for geographic transformation (optional). Here are some hints
for that step:
3-19
- SA_prj4 : Minna_To_WGS_1984_2
- SA_prj5 : Schwarzeck_To_WGS_1984
- SA_prj6 : Point_Noire_To_WGS_1984_2
What a geodesist can do is to find a substitute surface or a mathematical model that can
approximate the Earth's surface accurately enough to create a map. The simplest model
of the Earth is a sphere, but the accuracy of this model is usually not high enough
except when making a small scale map, such as 1/1,000,000. A more complicated model
is the geoid, an idealization by expanding the oceans to stretch below the landmasses
and making the assumption that the water there can also flow freely.
If we then neglect tidal and current effects on this "global ocean", a geoid was expected
to be a smooth, regular shape, however, it isn't. Water surface is affected by gravity,
while the direction of gravity depends on the mass distribution inside the Earth. Due to
irregularities or mass variation in distribution, the "global ocean" is forced to be an
undulated surface as shown in the following concept diagram. Consequently, the geoid
is not an analytic surface and it is thereby not suitable as a reference surface to
3-20
determine locations or create maps.
On the other hand, it is known that if the Earth was of uniform density and the Earth's
topography didn't exist, the geoid would have the shape of an ellipsoid.
The following figure offers the concept of the real Earth surface, geoid surface and
ellipsoid surface. The geoid surface is far smoother than the real surface; however, the
deviation of the geoid to the ellipsoid surface can be as large as +/- 100 meters. The
deviation between the geoid and an ellipsoid is called the Geoid undulation.
Because the ellipsoid exists only in theory and not in real life, it is completely smooth
and does not take into account any irregularities, such as mountains or valleys.
3-21
The shape of an ellipsoid is usually defined by its semi-major axis a and semi-minor
axis b; or semi-major axis a and flattening f. Flattening f is defined by the following
formula.
f = (a - b) / a
Datums
As mentioned above, the ellipsoid model of the Earth is based on the geoid and would
deviate from the geoid in various locations. Therefore, the reference ellipsoid used for
various regions of the Earth are usually determined by the local geoid as shown in the
following figure.
3-22
b) A defined orientation, position and scale of the geodetic system.
Datums are the basis for all geodetic survey work. They act as reference points in the
same way that starting points do when you give someone directions. For instance, when
you want to tell someone how to get to your house, you give them a starting point that
they know, like a road or a building.
Main ellipsoids
The coordinates of a point in the Earth will differ when different datums are used to
define the location of the point. If we plot these coordinates for the same point, for
example, the location of the Texas Capital Dome, as shown in the following figure, the
position shift from datum differences would be as large as about 1km.
3-23
The following figure gives the difference expressions of the Ethiopian boundary from
different datums. A base map uses the datum of WGS84, while other datums are either
used in Ethiopia or neighboring countries. The location shift between different datums
can be as much as 300m or more.
Geographic Coordinates
3-24
Projection
The term 'projection' comes from the notion of placing a light source inside a
transparent globe and projecting shadows of the meridians (longitude), parallels
(latitude) and other geographic features onto a sheet of paper. Mathematically speaking,
map projections are
transformations of geographic
coordinates (latitude, longitude)
into the Cartesian (x, y)
coordinates of a map.
Light Source
Different projections can be produced by changing the position of the light source.
Gnomonic projections place the light source at the centre of the globe, Stereographic
projections place the light source at the antipode of the point of tangency, and
Orthographic projections place the light source an infinite distance from the point of
tangency, resulting in parallel light rays.
3-25
Projection Surface
Map projections can be thought of as being drawn on one of three developable surfaces
- plane, cylinder and cone. Although a cylinder and a cone are not flat surfaces, either
can be flattened by cutting along its length and then unrolling the surface.
Aspect
Aspect refers to the orientation of the projection surface relative to the generating globe.
Polar Aspect (Normal aspect) aligns the axis of rotation of the projection surface with
the rotation axis of the Earth. Equatorial Aspect (Transverse aspect) rotates the
projection surface 90 degrees relative to the Earth’s rotation axis. The result is a plane
tangent at some point on the Equator, a cylinder tangent along an opposing pair of
meridians of longitude, or a cone whose axis of rotation lies in the plane of the Equator.
Oblique Aspect refers to any other aspects and are frequently used to generate maps
centred on particular locations.
The following figure gives the basic combinations of surfaces and aspects.
3-26
Distortion of Projections
Even though it is not possible to represent all distances or all angles correctly on a
single map, the map can be created to keep one of the selected geometric properties.
Consideration of basic properties results in four classes of map projections: equivalent,
equidistant, azimuthal and conformal.
Equivalent Projections
Maps that maintain constant areal scale over the entire surface.
Equidistant Projections
Equidistant projections preserve correct linear scale over some portion of the map. It is
not possible to represent all distances correctly.
Azimuthal Projections
Conformal Projections
The defining property of conformal projections is that at any point on the map, the scale
is the same in all directions about the point. The implication of this property is that
angles about the point are shown correctly and one would therefore expect shapes of
areas to be correctly represented. However, because map scale varies from point to point,
the property of correct representation of angles applies only to angles having infinitely
short sides. Shapes of infinitely small areas are preserved on the map but shapes of
larger areas are distorted.
3-27
Compromise Projections
While many map projections seek to optimize one of the above four geometric
properties, some projections do not attempt to maintain any of the above properties but
instead attempt to seek a balance between different properties.
Method of Construction
Projections can also be classified based on their method of construction. This involves
consideration of the type of projection surface used as well as the positioning and
orientation of the projection surface relative to the Earth.
The type of projection surface determines the basic pattern of the map. Positioning the
projection surface tangent or secant to the generating globe modifies the pattern of
distortion in the map. Re-orienting the projection surface can be used to achieve
minimal distortion within the region of greatest interest.
With increasing demands for global surveying activities, several worldwide useful
datum and projections have been developed. The World Geodetic System 1984
(WGS84) is the most famous and widely used of the world wide datums.
Basic parameters of WGS84 are shown in the following figure. WGS84 provides the
basic reference frame for GPS (Global Positioning System) measurements and almost
all open source GIS data are created on the basis of the GWS84.
Two kinds of worldwide coordinate systems Geographic and UTM have been
developed on the basis of the WGS84.
3-28
Geographic Coordinates
Geographic coordinates are expressed in terms of latitude and longitude as have been
mentioned in the previous part of this section.
UTM Coordinates
3-29
3-3 Data Integrating (Image Georeference)
Add Raster Dataset and Reference Data
1. Start ArcMap. Open a new map and save it as Georeference into the Maps folder.
3. Add the polygon feature of Eth from your GeoDB of Ethiopia/ Geographic_C.
4. Set the symbology of the layer Eth as no color for Fill and red color for Outline.
5. Add the image feature Ethiopia.jpg into the map from the GeoDB Image.
3-30
You already know that the coordinate system for feature class Eth is defined as
WGS1984, while the image Ethiopia.jpg does not contain a proper coordinate
system. The information of the coordinate system has to be added otherwise the
image layer can not match the feature layer of Eth. The process of defining a
coordinate system for an image requires a reference with a defined coordinate
system. In this exercise, you can take the layer of Eth as a reference. The process
used in this type of matching is usually referred to as georeference.
7. Right-click anywhere on the toolbar and a dropdown menu of ArcMap will appear.
8. Click Georeferencing to open the Georeferencing toolbar. Drag and drop the
Georeferencing toolbar somewhere it will be convenient to use.
9. In the Georeferencing toolbar, click the dropdown arrow next to the Layer textbox
and select the raster layer you want to georeference. Here it is ethiopia
Because you have only one raster datum in your map, the raster data of Ethiopia
will be automatically selected as the target layer of georeferencing.
11. Click dropdown arrow next to Georeferencing > click Fit To Display
The image size and location are changed to the full extent View as shown in the
following image. Since the view has been set to the full extent of the reference layer
Eth, the image layer is actually at an extent similar to the reference layer of Eth.
3-31
Add Control Points
To conduct georeferencing the basic method is to find several points in the target image
and their corresponding points in the reference layer. These points are referred to as
control points.
1. Click the Window menu in the toolbar > click Magnifier. A Magnification
window will appear.
2. Drag the magnification window to a location between a pair of control points in the
image and reference layer.
3. Click the Control Points tool in the Georeferencing tool bar. The mouse
pointer turns into a crosshair.
4. Center the crosshair over one control point in the image and click. A green control
point is added to the image.
5. Move the mouse pointer away from the control point but don’t click again. A line
stretches from the control point as you move the mouse. This is the link with which
you’ll connect the other end of the line to the corresponding point in the reference
layer.
6. Center the mouse pointer over the corresponding point on the layer of Eth - you
can see the link stretch as you do this - and click.
As soon as you click on the control point in the reference layer the map will move
to make the point in the image layer fit the point in the reference layer.
7. Drag and drop the Magnifier window to the next pair of control points.
3-32
If necessary you can use the Zoom In button or Zoom Out button to
make it easier to find the next control point.
8. Center the crosshair over the second control point on the image and click. Then
center the crosshair over the corresponding control point on the reference layer and
click to add the second control point.
9. By repeating the steps 4 to 8, you can add enough control points to make the image
fit the reference layer well.
Perform Georeferencing
If you are satisfied with the result you see in the Data
View window, save the result by clicking the
dropdown arrow of Georeferencing and then clicking
Update Georeferencing.
The information of the coordinate system will be saved into the image. You can check it
> remove it from the data frame and then Add it again. This time you will find the
image is located in the Data View just the same as it was before you removed it.
The above exercise is the simplest case of georeferencing because the image was
created by using the same coordinate system as the reference layer. The reason the
image did not match the reference before georeferencing is because when the image was
saved a projection file to store coordinate system information was not created.
Georeferencing by Transformation
On the other hand, sometimes the image is created in a coordinate system different from
the reference layer. Also, it may be seriously deformed as a consequence of the original
paper map being stored poorly or improperly scanned. In such cases you will find it
impossible to make the image neatly match the reference layer no matter how many
control points you add between the image and the reference layer. For this kind of image
georeference you need to conduct a transform process.
3-33
2. Add another image EthiopiaN. Repeat the procedure in the previous exercise to
make the image fit the display and then add a sufficient number of control points.
3. You will find that the image Ethiopa.jp could not be adjusted to match the
reference simply by adding control points.
4. When you have added enough control points, you must pay special attention to add
control points with homogeneous distribution in the image as much as possible.
Here “enough” is a relative term with the minimum number being 3 points, but
there is no maximum number. How many control points are needed for
georeferencing depends on the condition of the image. The more deformed the
image is, the more control points will be needed. Experience shows a number
between 10-20 is enough to obtain acceptable georeferencing results.
5. Click the dropdown arrow for Georeferencing > click Transformation > select a
method from the three orders of Polynomial.
The minimum number of control points needed for a different polynomials order is:
3-34
Assignment 5: Georeference a Base Map for Your Country
1. Add the base map of your country in ArcMap.
The base map of your country is an image in .jpg format with the name of
OOObase. Find the image from the folder:
C:/EthiopiaArcGIS/Data/CountriesData/OOO
2. Add the country boundary file into the same data set as the reference layer.
4. When you finish georeferencing, use the image as the background and the country
boundary file as the overlay layer. Set the symbology of the country boundary as no
color for Fill, red (or some other color that stands out) for Outline. The width
setting of the outline depends on your preference.
6. Export the map in .jpg format to the folder C:/GIS_Exercise/Result and set the
resolution to make sure pixels for both width and height are greater than 1000.
Name the export file as follows:
7. Upload the resulting image to the folder labeled with you name on the server
computer.
3-35
3-4: Point Feature Creation
When you finish georeferencing an image, the image becomes available in a number of
ways to support map creation or data creation. You can use the image as a background
for your map or create various kinds of vector data based on the image. Extracting point
data from a paper map is one of the essential processes in GIS data preparation (or
creation).
In this exercise you will become familiar with how to extract point data an image (i.e. a
scanned paper map).
1. Start ArcMap > open a new map > set the coordinate system of the data frame to
WGS84 > save the map as Temp (Overwrite the old one).
2. Add the image map_ethiopia_small.tif into the map from your GeoDB Image.
Within this image you will find the main cities of Ethiopia. Your task in this
exercise is to extract these cities from the image to create a point feature of the
main cities. The image has been georeferenced in the WGS84 coordinate system.
4. Expand the menu for Data Management Tools > expand Feature Class > open
the Create Feature Class tool. The Create Feature Class window will appear.
6. After the above settings have been made click the OK button. A New empty layer
Maincities is added into the map.
3-36
Extract Points from an Image
1. If the Editor toolbar is not showing click the Editor Toolbar button .
2. Start Editing.
3. Make sure the Task is Create New Feature and the Target layer is Maincities.
4. Click the Sketch Tool > move your curser into the Data View. Your curser
3-37
5. Put your curser on one point of a city and click. A new point is added into the Data
View and highlighted in blue.
6. Move your curser to a point of a another city, clicking them one by one.
9. Zoom in so that it is possible to check if the points were added to the best location.
10. To adjust a point, change the editing Task to Reshape Features > click the Editor
Tool > select a point you would like to move > move it to the best location by
11. After you finish adjustmenting all the points Stop Editing and save the Edits.
Your task in this exercise is to create a point feature for maincities; at least one field is
necessary to store the name of the cities.
1. Open the attribute table of the layer Maincities > click the Options button > click
Add Field. The Add Field dialog box will appear.
- Name : Name
- Type : Text
- Field Properties/Length: 25
3-38
3. Click the OK button. Now you have to input the name of the cities into the new
field you added.
7. Click the Attributes Tool . The Attribute dialog box will appear. Click the cell
within the Value column that corresponds to the Name field > referring to the
background image input the city name into the cell.
8. Repeat Step 7 to add all the city names into the attribute table.
9. After you finish adding attributes of city names, close the Attributes dialog box.
10. Open the Properties of the layer Maincities > click the Labels tab > check the box
for Label features in this layer > accept the default Method of Label all the
features the same
way > set Name as
Label Field > (if you
like you can also set
the Text symbol)
>click OK to close
the Properties
window.
3-39
11. Right-click the layer of Maincities > click Label
Features. All the data you have added into the Name
field of the attribute table will appear. Check if you have
properly added the names of the cities.
Using a table is the most common method to store point data in our work. You might
have tables for wells, water quality survey results, village information, and so on.
Within some tables the location information is generally stored in two columns, X and Y.
Certainly, for different coordinate systems, X and Y can be expressed in different ways
such as longitude and latitude for the Geographic coordinate system or Easting and
Northing for the UTM coordinate system. Given that you have at least these two kinds
of coordinate information, you can create a point feature using a table.
One of the most popular formats of tabular data is the data sheet using MS Excel,
therefore a *.xls table is used as a sample for this exercise.
C:/EthiopiaArcGIS/Data/OtherData.
WQ_AA.xls is a data table for storing groundwater quality analysis results of many
wells in the Addis Ababa area.
You can find the first 3 columns are initiated by Long followed by 3 columns
initiated by Lati. You might have already made out that the meanings of Long and
Lati are abbreviations of Longitude and Latitude, respectively. These 6 columns are
essential coordinate information to add this table into ArcMap.
However, this coordinate information can not be used to add the table into ArcGIS
because ArcGIS can only recognize decimal format no matter if the coordinate
system is Geographic (longitude/latitude) or UTM (meter). Therefore, you have to
perform data format conversion according to the followed steps.
To the right side of the coordinate information you will find several water quality
3-40
analyses results are stored, such as TDS, EC and so on.
3. Use the following equation to change the coordinate information from separated
degree, minute and second into decimal degree:
- input the above equation into the top cell of the column
- delete the columns in which the dataset for longitude and latitude are
stored.
- set the format of these two column as Value with 6 decimal places
The above format setting is needed only to save the data sheet as a *.dbf file. If you
would like to save the data sheet as a *.txt file, you can set the format for all data as
standard.
4. After you finish the data format conversion in Step 3, you should save the data
sheet into text format or *.dbf (dBase) format since those are the only two formats
in which ArcGIS can read table data.
- click Save As
3-41
- browse to the folder My_data
- click Save
1. Start ArcMap > open a new map > set the coordinate system of the data frame as
WGS1984 > save the map as Temp (Overwrite the old one)
3. Click the Tools menu > click Add XY Data. The Add XY Data window will
appear.
- click Add
You will find that if the columns for coordinate information are titled as longitude
and latitude, they will be automatically selected by the Add XY Data tool.
If the columns for storing coordinate data are titled X and Y, they will also be
recognized automatically. However, if the columns of coordinate data are titled
with other words, you have to select them by clicking the dropdown arrow next to
the textboxes of X Field and Y Field.
5. Click the Edit button to define the coordinate system; in this exercise WGS_1984.
3-42
6. Click OK.
At the present moment you just finished adding the XY data into ArcMap, but the
layer you can see from your map is an in-memory feature. That is, it can not be
used for other maps or programs.
8. Right-click the layer of QW_AA_data Event. The dropdown menu will appear.
- Click Export Data. The Export Data dialog box will appear.
- Click OK
3-43
9. When the exporting is complete you will be asked if you would like to add the
export feature to the map. Click Yes.
1. Start ArcCatalog > browse to the folder My_data > right-click the text file
WQ_AA_data.txt > click Create Feature Class > click From XY Table. The
corresponding window will appear.
2. The same way you added XY data into ArcMap, the two columns of coordinates are
automatically selected into the text boxes of X field and Y Field.
4. Set the output location as your Temp folder, and set the output name to be the same
as the original text file, WQ_AA_data, then click OK.
5. Browse to the Temp folder, although you might not find the feature WQ_AA_data
that you just created. If so, click the View menu > click Refresh. The feature will
appear in your folder.
6. You can drag and drop it into ArcMap to compare the feature WQ_AA_data with
Wq_aa, and you will find them to be identical.
3-44
Assignment 6: Create a Point File from Place Name Data
1. Add your country’s boundary file in ArcMap.
Use the country boundary file that you created in Assignment 3. Set the symbology
as you like, but avoid using “no color” as a color setting.
C:/EthiopiaArcGIS/Data/CountriesData/OOO
http://earth-info.nga.mil/gns/html/cntry_files.html
Please refer to the table attached to this section to get the details of the field in the
place name file.
3. Add your country’s place name file into ArcMap as an event layer.
5. Add the exported feature into the same dataset. Inactivate the events layer (although
in the meantime the events layer should cover the new layer exactly). Set the
symbology of the new layer as you like.
6. If there are points outside the country’s boundary, use the boundary layer to clip the
point layer.
To clip you can use the Clip tool in the ArcToolbox menu for Analysis Tools/
Extract
7. Use the column of FC in the new layer as the target field to create a selection based
on one value in the field, for example, A.
3-45
8. Export the map in .jpg format into the folder C:/GIS_Exercise/Result and set the
resolution to make sure pixels for both width and height are greater than 1000, with
the file name as follows:
Yourlastname_E5. (such as Lei_E5.jpg)
9. Upload your resulting image to folder with your name in the server computer.
3-46
Field Field Default
Field Description
Name Type Selection
negative sign (-) = West.
Latitude of the Geoname feature in ±degrees, minutes, and
DMS_ seconds; DMS; (±ddmmss):
number Y
LAT no sign (+) = North;
negative sign (-) = South.
Longitude of the Geoname feature in ±degrees, minutes, and
DMS_ seconds; DMS; (±dddmmss) :
number Y
LONG no sign (+) = East;
negative sign (-) = West.
Universal Transverse Mercator coordinate grid reference. (This
field is derived from the system generated and stored Military
UTM character Y
Grid Reference System (MGRS) coordinates when this file is
generated.)
JOG Joint Operations Graphic reference. character Y
Feature Classification: Nine (9) major Geoname feature categories
into which similar feature designations are grouped.
A = Administrative region type feature
P = Populated place type feature
V = Vegetation type feature
FC L = Locality or area type feature character Y
U = Undersea type feature
R = Streets, highways, roads, or railroad type feature
T = Hypsographic type feature
H = Hydrographic type feature
S = Spot type feature
Feature Designation Code. A two to five-character code used to
DSG character
identify the type of Geoname feature a name is applied to.
Populated Place Classification. A graduated numerical scale
denoting the relative importance of a populated place. The scale
PC ranges from 1, relatively high, to 5, relatively low. The scale number
could also include NULL (no value) as a value for populated
places with unknown or undetermined classification.
Primary Country Code. A two alphabetic character code from the
CC1 FIPS 10-4 Standard that uniquely identify a geopolitical entity character
(countries, dependencies, and areas of special sovereignty).
First-order administrative division code. A two alpha-numeric
code from the FIPS 10-4 Standard describing a primary
ADM1 character
administrative division of a geopolitical entity, such as a state in
the United States.
Second-order administrative division. The name of a subdivision
variable
ADM2 of a first-order administrative division, such as a county in the
character
United States.
POP Population figures. number
ELEV Elevation in meters. number
Secondary Country Code. A two alphabetic character code from
the FIPS 10-4 Standard uniquely identifying the country code of a
CC2 particular Geoname feature name if different than that of the character
primary country code of that Geoname feature (International
features).
Name Type:
NT C = Conventional name character
N = BGN Standard name
3-47
Field Field Default
Field Description
Name Type Selection
NS = BGN Standard name in non-Roman script
P = Provisional name
PS = Provisional name in non-Roman script
H = Historic name
HS = Historic name in non-Roman script
D = Not verified or daggered name
DS = Not verified name in non-Roman script
V = Variant or alternate name
VS = Variant name in non-Roman script
Language Code. A three alphabetic character code (ISO 639-3)
LC uniquely identifying the language assigned to a Geoname feature character
name.
SHORT_ variable
A part of the full name that could substitute for the full name.
FORM character
GENERI The descriptive part of the full name such as Mount, River, or Sea variable
C (generally does not apply to populated place names). character
A form of the full name that allows for alphabetical sorting of the
file into gazetteer sequence. All character/diacritic combinations
and special characters are substituted with QWERTY (visible U.S.
English keyboard) characters, all characters are upper-cased,
SORT_ variable
numerals are converted to lower-case characters (0-9 = a-j), spaces
NAME character
are stripped out, and commas replaced with space. The
substitution is done according to a character mapping table that is
Region Font Code specific. This field is included for the benefit of
the end user of the data to aid in the sorting of names if required.
The full name is the complete name that identifies a named
FULL_N feature. The full name is output in reading order, "Mount variable
AME Everest", vs. reversed generic, "Everest, Mount", as stored in the character
database.
Same as the full name but the character/diacritic combinations and
special characters are substituted with QWERTY (visible U.S.
FULL_
English keyboard) characters while still maintaining casing and variable
NAME_
spaces. The substitution is done according to a character mapping character
ND
table that is Region Font Code specific ND = No Diacritics /
Stripped Diacritics.
MODIFY The date a new feature was added or any part of an existing
date-time
DATE feature was modified (YYYY-MM-DD).
3-48
3-5: Polyline and Polygon Feature Creation
Exercise :1 Create a Polyline or Polygon Feature by Tracing an Image (optional)
The method of tracing a polyline or polygon from an image is basically the same as
creating a point feature from an image.
1. Start ArcMap, open a new map > set the coordinate system of the data frame, for
example WGS_1984 UTM Zone 37N > save it as Temp (overwrite the old one).
2. Add the reference image Baseline_primary_data_a3 into the map from your
GeoDB Image.
Your task in this exercise is to delineate the main roads within the study area of
EWTEC.
4. After the above settings, click OK. A new empty layer of SA_Mroad is added to
the map.
3-49
8. Edit the traced result (refer to the Exercise 3 in Module 2, Section 5)
The procedure of creating a polygon feature by tracing an image is almost the same as
the above exercise. The difference is:
When your base map image contains a lot of lines or polygons, the tracing method as
shown in the above exercise would be a very tedious and time consuming task.
Therefore, as you might expect, there are some methods to carry out the delineating
process more effectively and more automatically. The following exercise will introduce
some of these methods.
Preparation
1. Scan the paper as an image with the extension *.tif, *.jpg, *.bmp or other formats. In
this exercise you should have already completed scanning and have the file
afoutline.tif.
2. Start ArcMap, open a new map, set the WGS84 as the coordinate system for data
frame, save the map as Temp (overwrite the old one).
You will find the image shows the boundary lines for all African countries. In this
exercise, you will try to delineate these countries’ boundary lines to create a
polygon feature.
4. Start ArcCatalog > browse to and right-click your Temp folder > click New > click
Shape file. The Create New Shape file window will appear.
- Name: templine
- Feature Type : polyline
- Spatial Reference: WGS_1984.
3-50
6. Click OK. A new empty polyline shape templine file will be added into your Temp
folder.
7. Drag and drop the newly created shape file templine into ArcMap.
8. Referring to Steps 4 to 7, create and add a new empty polygon tempgon.shp into
ArcMap.
1. Go to ArcMap > open the Layer properties window of the layer afoutline > click
the Symbology tab > set the Show as Classified and then set the classes as 2. Click
OK.
3. Start Editing. Make sure the Task is Create New Feature and the Target is
templine that you created and added earlier in this exercise.
4. Click the Tools menu > click Extensions.... make sure the extension ArcScan has
been checked.
5. Right-click anywhere on the tool bar or menu bar to bring up the toolbar list, select
ArcScan, if it is unselected.
6. Click the dropdown arrow next to Vectorization in the ArcScan Toolbar > click
Vectorization Setting.... The Vectorization Settings window will appear.
3-51
7. Specify the values for each parameter as shown in the figure, and then click Apply.
Notice: If you don’t set any parameter in the above window the Apply button will
remain unavailable.
8. Don’t close the Vectorization Settings window. Move it to one side of the display
to allow you to see the Data View.
9. Click the dropdown arrow next to Vectorization and then click Show PreView.
10. Generally the result will be acceptable, but if you are not satisfied with the result
you can go back to the Vectorization Settings window and change the parameter
setting, then regenerate the lines for the layer of templine.
In the full extent of Data View, you might find it difficult to confirm if the
vectorization has been performed properly because in some part the red line appears
to deviate from the line in the background image (refer to the left side of the figure
below). Then you can zoom in to a scale large enough to confirm the result of
vectorization (refer to the right side of the figure below)
3-52
11. After you have checked the result and find it acceptable, close the Vectorization
Settings window.
12. Click the dropdown arrow next to Vectorization and then click Generate Features
to bring up the Generate Features window.
13. Make sure the generated line class will be added into the layer of templine, accept
the default setting for “Generate polygon where the maximum line width setting
is exceeded”, and make sure those polygons will be added into the newly created
empty polygon layer tempgon.
14. Click OK and wait until the vectorization result is added into the layer of templine.
Since you have completed the vectorization, it is suggested you hide the ArcScan
3-53
toolbar to keep your window uncluttered and easy to work in.
You will find that when you stopped editing, the ArcScan tools of were dimmed. In
other words, the extension of ArcScan only work in editing mode.
If your task were only to extract lines from an image to create a polyline feature, you
have completed your task with the exercise above. However, the source map in this
exercise is a boundary map of African countries, and generally the country’s boundary is
expressed as a polygon shape. Therefore, the results from the last exercise should be
further processed to convert to a polygon feature.
2. Expand the Data Management Tools menu > expand Features, open the Feature
to Polygon tool. The Feature To Polygon window will appear.
3-54
4. Click OK and wait a moment for the conversion process to complete.
5. The new polygon feature will be added into your map and then you can check and
modify it to complete the vectorization.
You can also use ArcCatalog to perform the feature type conversion from polyline to
polygon with the following steps.
1. Start ArcCatalog.
2. Right-click your GeoDB OtherData > click New > click Feature Dataset > name
the new feature dataset as temp and set the coordinate to WGS_1984.
3. Import the polyline feature templine from your Temp folder to the new dataset
temp.
4. Right-click the new feature dataset of temp > click New > click Polygon Feature
Class From Lines…. The Polygon Feature Class From Line window will appear.
5. In the Enter a name for the feature class textbox input a name such as
Af_cntry_temp1 for the new polygon feature. Then check the only polyline file
3-55
Templine in the Select the features classes that will contribute lines feature list
box. Click OK.
You have completed the creation of a polygon feature for African countries, but the
feature is not particularly useful at this moment because it doesn't contain any of the
necessary attributes. The feature of Af_cntry_temp will not even tell you the name of
the countries. Therefore, the next task should be adding the necessary attributes to the
polygon feature of Af_cntry_temp. There are many methods to add attributes to a
feature; a common method is to add attribute data by obtaining them from another
related feature.
In this exercise, you will learn the process of adding attributes from a point feature to a
polygon.
This point feature gives the location of centroid points of all countries on the
African continent and stores basic attributes like country name, area and so on.
2. Right-click the layer of Af_cntry_temp > click Joins and Relates > click Join.
The Join Data window will appear.
3-56
3. In the top textbox, select “Join data from another layer based on spatial
location”.
5. Check the second radio button to make all the attributes of the point layer cntry_pt
added into the polygon
layer Af_cntry_temp.
- Name: Temp
- Save as type:
Shapefile
8. A new layer Temp is saved in the destination folder and added to your data frame.
This layer is a polygon feature, and has the same spatial extent as the layer of
Af_cntry_temp.
9. Open the attribute table of the layer Temp, you will find the attribute data has been
added from the cntry_pt into Temp.
3-57
The field will go to the leftmost side of the attribute table and when you move the
button of the slide bar in a horizontal direction so that the COUNTRY field will
not disappear. It will now be easy for you to compare the data in field of
COUNTRY with data in any other field.
11. Sort the COUNTRY field in ascending order, make width of the attribute table
window as narrow as can to show only two fields, and make the last field of
Distance shown together with the field of COUNTRY.
- Several countries for which you have created two or more polygons, so
you can to combine (merge) them into one.
- Every country that has a feature with the data value 0 in field of Distance.
This value indicates to you that there is a point from layer cntry_pt
located within the polygon feature.
These values indicate to you that there are no points from layer cntry_pt
located within the corresponding polygon features. You will find that all
polygon features with these values correspond to more than one of the
same country name in the COUNTRY field. Therefore, these features
should be modified by combining them with another polygon feature or
deleting them.
12. Based on the indications in the attribute table, it will be easy for you to find and
reshape the polygon layer Af_cntry_temp.
- Inactivate all layers except the background image afoutline and Temp.
- Put the Attributes table window adjacent to the Data View so that you
can easily change between them.
- Start Editing. Set Task to Reshape Feature and the Target layer to
Temp.
3-58
- Zoom in to the selected features to check if all the features should be
combined into one country. If so, merge them; if not, decide how you will
deal with it. In this exercise you have only one option other than merging,
which is deleting. In other cases you might have another option, leaving it
without reshape. In that case you have to add attribute data by another
method.
13. After you finish editing the polygon feature, you can complete your task by Stop
Editing > Save the editing result > export the edited polygon feature to your data
base.
In the above exercises you practiced how to convert a single color image into a line
feature or even a polygon feature. However, in your stocks of paper maps, which will be
used to create GIS data, you may have some color classified maps, such as geology
maps, landuse maps and so on. To make color classified maps available in our GIS, the
method of using ArcScan is not a suitable selection because the method is based on the
two class classification for symbology of the map.
In this exercise, one nonstandard method will be introduced for your reference.
3-59
Preparation
1. Start ArcMap > open a new map > set coordinate system of the data frame to:
WGS_1984 UTM Zone 37N > save the map as Temp (overwrite the old one).
The image is a geology map which covers the majority of the study area of EWTEC.
Just glancing at the map, you might imagine how hard it would be to manually trace
all polygons from the map.
2. Expand the Conversion Tools menu > expand From Raster > open the Raster to
Polygon tool. The corresponding window will appear.
- Field: Value (the tool will automatically set this factor for you)
The converted result of Temp_geology will be added into the data frame and you
can find the boundary line zigzags and is too thick. If you check the attribute data of
the layer, you will also find the number of polygons included in the layer is as
3-60
much as 450,000 or greater. If you check the size of the shape file from the folder
you saved it to, you will find it is as big as about 88M.
1. Open the attribute table of the layer Temp_geology > click Options > click Add
Field. The Add Field dialog box will appear.
5. Sort the Area field in ascending direction. You should find that most polygons have
a very small area, less than 1000. These features consist of boundary and slivers.
The geological entities are those features with relatively large areas.
3-61
Remove Large Area Features
1. Click the Selection menu, make the Temp_geology layer the only selectable layer.
3. Confirm that the layer Temp_geology is the target layer, set the selection equation
as: "Area" > 1000000. Then click the Apply button.
4. When the selection complete, close the Select By Attributes window. Check the
selected result.
The result of your selection will be highlighted in the Data View. All the main
geological entities have been selected except a very small one on the upper lefthand
corner.
Also, you can check the attribute table to see that your selection contains about 32
features.
5. Start Editing, set the Task to Reshape Feature and the Target layer as
Temp_geology.
6. Click the Delete button in the toolbar. When deletion is complete, the main
geological entities and the boundaries in left, top and right directions disappear. The
Data View will look like the figure shown below.
7. Stop Editing.
3-62
Add a Frame and Export the Data View
1. Set the Symbology of the layer Temp_gfology as black for Fill and no color for
Outline.
2. Display the Draw tool > select the New Rectangle tool > draw a rectangle in
the Data View.
3. Make the Outline of the rectangle black, 1 for outline width, and no color for Fill.
4. Zoom in to each side of the rectangle and adjust it to make it fit the boundary of the
layer Temp_gfology.
5. After you finish adjusting the rectangle, hide the Draw toolbar.
6. Adjust the extent of your Data View to make it fit the layer of Temp_gfology and
the rectangle you drew. The Data View will look like the figure shown below.
7. Click the File menu > click Export Map.... The Export Map window will appear.
Input the settings within the window as follows:
- Name: Temp_geology
- Type: TIFF(*.tif)
3-63
- Resolution: Adjust to make the pixels for both height and width greater than
1500.
- Write World File: checked
In the Format tab:
- Color Mode : 1-bit Monochrome Threshold
- Compression: None
8. Add the exported result of Temp_geology. You will be asked if you would like to
build a Pyramid, click Yes.
Generally, the exported image should have a similar spatial extent as the original
Data View. However, sometimes you might find the spatial extent of the
exported .tif format image is not exactly the same as its original View. If this is the
case, you should conduct georeferencing by the method you practiced in Section 3
of this Module.
The procedure of this part is almost the same as you practiced to vectorize the image of
African countries in Exercise 2 of this section. So only a digest is given in this text
book as a reference for you to practice vectorization.
2. Expand Data Management > expand Feature class > open Create Feature Class
3. Refer to the process above to create a new, empty polygon feature Temp_gon.
Both new, empty polyline and polygon features are added into the same data frame
as the image Temp_geology.tif.
3-64
4. Set the symbology of the layer Temp_line as red. And set the symbology of the
layer Temp_gon as pink for outline; no color for fill.
5. Start Editing. Set the Task to Create New Feature and the Target to the layer of
Temp_line.
6. Set the Symbology of the image Temp_geology.tif as classified and the classes as
2.
7. Display the ArcScan toolbar and make sure the Target Raster is
Temp_geology.tif.
9. Set the Maximum Line Width as 50. (Even though the setting might be 50, you
have to set it once again, otherwise the Apply button will not function). Click the
Apply button.
10. Click the dropdown arrow next to Vectorization > Show Previous.
11. Click the dropdown arrow next to Vectorization > Generate Features. Set the
Temp_line and Temp_gon as the two necessary layers feature generation. Click
OK.
12. Stop Editing with the edited result saved. Hide the ArcScan toolbar.
13. When you check the result from Temp_line and Temp_gon, respectively you will
find that the main geological entities are stored in the layer of Temp_line (refer to
the upper side of the figure below), and the layer of Temp_gon contains mainly
slivers (refer to the lower side of the figure below). Maybe you would agree that the
polyline feature would be better than the polygon feature for use as a vectorization
result, even though you still have to edit over the polyline feature.
3-65
Vectorization Result Editing
The main task of polyline feature editing is to delete the extra line features that
correspond to slivers such as labels or defect segments that occurred in all of the above
processes. You can use the method that you practiced in Exercise 3, Section 5, Module
2. You can also utilize the characteristic that most slivers have a relatively short length
to simplify the editing process by the following procedure.
1. Start by converting the temp_line into a polygon feature. Display the ArcTool
window if it is hidden.
2. Expand Data Management Tools > expand Features > open Feature To Polygon.
In the Feature to Polygon window, input the settings as follows:
3-66
The new polygon will be added into the same data frame as its source data of
Temp_line. You will find many slivers were eliminated in the conversion process.
Only a few extra polygons representing labels and slivers need to be edited further
to delete them.
Although you can use the method polygon editing feature that you practiced in
Exercise 6, Section 5, Module 2, here your main purpose is to attach small-area
polygons to their adjacent large-area polygon so you can use the Elimination tool
of ArcGIS.
3. Add a new field and calculate the area as you have practiced several times already.
4. Sort the new field and define selection criteria to distinguish small sliver polygons
from real geological entities. You might set 7000m2 as a temporary value of the
criteria.
5. Perform Selection By Attributes to select all polygons with an area less than 7000
in the layer of geo_line2gon.
3-67
6. Expand of the Data Management Tools menu > expand Generalization > open
the Eliminate tool.
The eliminate result geo_line2gon_E is added into your data frame. You can check
the result of the elimination and find the there are less than 100 polygons remaining
in the layer geo_line2gon_E.
8. You can eliminate slivers further by selecting the slivers directly from Data View
or setting new selection criteria.
For example, if you set the new selection criteria as area < 1,000,000 and then
perform elimination once again you would get a result as shown in the figure
below.
3-68
The new eliminated result contains 32 features; to check all of them would not be
too difficult a task. After you set the symbology of the new eliminated result as no
color for fill and red for outline, you will be able to compare the eliminated result
with the original geology map. You will find that you still need to do some manual
editing to complete the polygon feature creation from a color classed image.
In this exercise you practiced creating a polygon feature based on a color classified
image. One thing that you should understand is that right now there is no standard
method for color classified image vectorization. That is why high precision GIS data
requires that you perform digitizing using very expensive machinary and consume much
more time. The method introduced here is only to be used in the following conditions:
- You don’t have enough time to manually digitizing the paper map
- The precision of this semi-automatic method for digitizing is sufficient for your
needs of GIS data base creation.
A lot of features such as lakes, ponds, farm land, etc. are base data for water resources
management. For example, farmland distribution is one of the most important factors
for groundwater management because irrigation in many countries is the largest portion
of water use. In many countries excessive groundwater pumpage for irrigation has
3-69
caused serious groundwater head depression, and thereby resulted in a number of
environmental problems such as reduction in natural vegetation or even extinction,
abandonment of existing water use facilities (qanats, wells, etc.), expantion of
contaminated groundwater and so on. Therefore, irrigation water use control has
become a very important issue in groundwater management in these countries. To carry
our irrigation water use control, data concerning the distribution of farmland is
indispensable. But this kind of map is not a common one and might be difficult for you
to find and obtain. Therefore, you have to create such maps by yourselves.
Developments in remote sensing technology enable you to obtain the image data,
including farmland distribution. In this exercise you will become familiar with the basic
concept of satellite photograph processing.
Many types of satellite photography have become available for purchase, but most of
them are too expensive. Fortunately, one of the most famous satellite photograph types,
Landsat images from USNASA, have been released for free. You can download Landsat
imagery from the following homepage:
http://glcfapp.umiacs.umd.edu:8080/esdi/index.jsp
Another very useful source of data you should know is MrSID (Multi Resolution
Seamless Image Database). MrSID (pronounced Mister Sid) is one of the most excellent
formats for image compression available and thereby used by NATO, NASA USGS and
many other well-known organizations and agencies. MrSID has also become one of the
standard formats of GIS data and adopted by ESRI, Autodesk, Mapinfo, ERDAS and
other famous GIS applications.
Some of MrSID images have been released by NASA. For general utility, you can
obtain (download) 2 sets of MrSID data from the following Web site:
https://zulu.ssc.nasa.gov/mrsid/
The two sets of MrSID data are MrSID 1990 and MrSID 2000. Both of the data are
based on US satellite image Landsat and cover the whole world.
One scene from MrSID covers a far more vast extent than one scene from Landsat itself.
MrSID extends north-south over 5 degrees of latitude, and spans east-west for the full
width of the UTM zone (6 degrees of longitude). The pixel size of MrSID is 14.25m for
MrSID 2000 and 28.5m for MrSID 1990.
3-70
The spectral bands for MrSID are red, green and blue, which correspond to the Landsat
bands of 7, 4 and 2.
1. Start ArcCatalog > click the Tools menu > click Options > click the General tab >
uncheck the box for Hide File Extensions > click OK.
This is a MrSID scene which covers Addis Ababa and a vast area within Ethiopia.
3. Right-click the data in the directory tree > open the Properties of the data > move
the slide button to find the item Spatial Reference. You will find it is currently
Undefined.
4. Click the Edit button and the now familiar Spatial Reference Property window
will appear.
5. Set the coordinate to the UTM_WGS_1984, and the UTM zone, which is just the
same as the data name, here is 37N.
6. Start ArcMap > open a new map > set the coordinate system of the new map to
UTM_WGS_1984_Zone_37N > save the map as SatilleteImage in your Maps
folder.
7. Add the MrSID data N-37-05_2000.sid into your map, then add the features
A_Ababa and Eth as a reference.
3-71
8. Even you don’t carry out any more of the process, the map you have added into your
data set can be used as a background.
Tip: when you add the MrSID data into your map, you can add the individual band
form of red, green and blue. When you reach MrSID, double click it and the three
bands will be shown in your Add Data window; then you can select one of them to
add into your data frame.
To extract a feature from a satellite image you have to know the area you're interested in
clearly, otherwise you will not be able to select the color range to extract a special
feature. One water resource survey JICA project was carried out in the northwest part of
China, so the color range of farmland has been made clear there. This data will be used
here for the farmland extraction exercise.
1. Open a new map > set the coordinate system of the data frame to UTM_
WGS_1984_Zone45N > save the map as SatilleteImage (overwrite the old one).
C:/EthiopiaArcGIS/Data/MrSID/N-45-40
This feature you have used in the exercise of conversion polygon to centroid. The
feature gives the shape and extent of a simplified grid for a groundwater simulation
model. Your task here is to extract the farmland distribution in the simulation area
so that the irrigation amount can be calculated for model establishment.
3-72
4. Confirm the structure of MrSID within the Data View. Right-click the layer of
N-45-40.sid > click Data > click Export Data. The Export Raster Data window
will appear. You can confirm the structure of the data is as follows (don't save):
After confirming the structure, click the Cancel button to close the window.
5. Open the Properties of the layer N-45-40.sid > click the Symbology tab > make
sure the Show type is RGB Composite > select Stretch Type to Standard
Deviations > set the number of n to 0.01 > click OK.
You will find that the layer of N-45-40.sid changes its expression from full color to
something like a color classified map.
7. The next step is to convert the MrSID expression to an image using the map export
function with which you have become familiarized. However, here you have to
specify some factors that are different from what you set in your assignment.
3-73
- Set name and format of the output image to Sp_image.jpg
- Click the Format tab and then select the 8-bit grayscale as Color Mode
- Click the General tab > adjust the Resolution to make the width and
height of the pixels near the numbers of column and row that you checked
in Step 4.
- Don’t forget to put a check in the box for Write World File, otherwise the
exported result will not match other data in your data frame.
8. Add the exported image Sp_image.jpg. Check the value of the area that
corresponds to the farmland by using the Identify Tool . You will find the value
10. Expand the Spatial Analyst menu > expand Extraction > open the Extract by
Attribute tool. The corresponding tool window will appear.
11. In the Extract by Attributes window, set the necessary factors as follows:
3-74
- Click the OK button twice to perform the extraction.
12. The extracted result well be added into the data frame.
If you make all layers inactive except the newly created layer of Sp_fmland, you will
find the farmland area has been extracted, but it contains many small slivers. One
method to reduce these slivers is to resample the image.
1. Expand the toolbox menu of Data Management Tool, expand Raster > open the
Resample tool. The corresponding window will appear.
4. Expand the toolbox menu of Conversion Tools > expand From Raster > open the
Raster to Polygon tool.
3-75
- Input raster: Sp_farmland
6. Make the reference layer nss_gridbase active, set the symbology of the feature
layer Sp_farmland as red for fill; no color for outline; set the symbology of
MrSID layer to its default setting. The result of farmland area extraction from
MrSID data will appear as shown in the figure below.
1. Make sure the symbol of source raster data N-45-40_2000.sid is set to default.
(RGB Composite for Show and None for Stretch Type)
2. Make sure the Data View extent is the same as in your last exercise.
3. Check the extent of groundwater simulation model area by opening the Properties of
layer nss_gridbase > click Source tab.
3-76
As shown in the figure below, the limit of the groundwater simulation area is:
759000 for left, 4756000 for top, 769000 for right and 4749000 for bottom
5. Click the dropdown arrow next to Spatial Analyst > click Options > click the
Extent tab > click the dropdown arrow next to the textbox of Analysis extent and
then select As Specified Below > use the results from checking the extent to define
the analysis extent > click the OK button.
- Type: JPEG
3-77
8. Click the Add Data button > browse to the Temp folder > double-click the
newly exported raster data N_45_40_image > expand the image and then select all
three (3) bands > click the Add button. The three bands which represent the colors
red, green and blue are added into your data frame.
- Make layers of Band 2 and Band 3 inactive, then you can confirm that the
relatively small value in Band 1 (red) corresponds to the farmland area.
- Set the color of red for Class 1 (class for small value range), no color for
Class 2.
- Adjust the Break Value until almost all farmland areas are covered by
Class 1 of Band 1.
11. Click the dropdown arrow next to Spatial Analyst > click the Raster Calculator.
The corresponding window will appear. Referring to the equations in the following
figure, create farmland calculation equations. The values for extraction from each
band should be based on notes that you will have taken in Steps 9 and 10.
3-78
The basic rules for raster calculation are as follows:
12. All the equations you build in the Raster Calculator will create a raster data layer.
As the result of the calculation, 4 raster data layers of B1, B2, B3 and Farmland
are added into your data frame. B1 to B3 are interim results, which are not so useful
so you can remove them from your data frame.
13. The symbol of Farmland is automatically set as Unique Value. You can set all
values to no color except value of 3. According to your calculation, the value of the
raster data Farmland is as follows:
You will find this result obtained by raster calculation is a little better than the
result of the raster data simplifying method.
3-79
14. Referring to the previous exercise, complete the farmland polygon feature creation:
- Extract the farmland raster data by using the tool Extract by Attribute in
the ArcToolbox Spatial Analyst/ Extraction
- Convert the extracted farmland raster data to polygon feature by using the
Raster to Polygon tool in the ArcToolbox Conversion Tools/ From
Raster
3-80