Lecture ERS 120: Principles of GIS - Practice 5: Advanced Analysis I - N.D.Bình
It is possible to select all the cells with certain characteristics - in this example where the rating
equals '6'. With 'Map Calculation l ' active go to the 'Analysis' theme and select the 'Map Query'
function. Enter the following query, by double clicking on 'Map Calculation L" (top layer), then
single clicking on the '=' operator and finally double clicking on the unique value of '6': -
[Map Calculation] = 6.AsGrid
Click on the 'evaluate' button and then close the 'map query' window. A new Boolean theme
will be produced - indicating those cells that fulfill the map query criteria (grid cell will become
'true' = 1) and those that do not ('false' = 0). Your final map of the most suitable locations should
resemble the figure illustrated below:
You can now close ArcView. There is no need to save any files - they are no longer required.
All files should have been saved to the temp directory and as such will automatically be deleted
when the computer is re-started.
9 trang |
Chia sẻ: thucuc2301 | Lượt xem: 657 | Lượt tải: 0
Bạn đang xem nội dung tài liệu Lecture ERS 120: Principles of GIS - Practice 5: Advanced Analysis I - N.D.Bình, để tải tài liệu về máy bạn click vào nút DOWNLOAD ở trên
Practice 5: Advanced Analysis I
ArcView Steps ................................................................................................................ 26
Step 1 GeoProcessing - Dissolve .............................................................................. 26
Step 2 GeoProcessing - Merge ................................................................................. 27
Step 3 GeoProcessing - Clip ..................................................................................... 27
Step 4 GeoProcessing - Intersect .............................................................................. 28
Step 5 GeoProcessing - Union .................................................................................. 29
Step 6 GeoProcessing - Spatial Join ......................................................................... 30
Step 7 Spatial Analyst - Find Distance ...................................................................... 30
Step 8 Spatial Analyst - Assign Proximity .................................................................. 31
Step 9 Spatial Analyst - Density ................................................................................ 31
Step 10 Spatial Analyst - Summarise Zones ............................................................. 32
Step 11 Spatial Analyst - Reclassify .......................................................................... 32
Step 12 Spatial Analyst - Map Calculator .................................................................. 33
Step 13 Spatial Analyst - Map Query ........................................................................ 33
In this practical you will discover how to manipulate your spatial data through the use of the
GeoProcessing wizard; as well as exploring the main functionality of the 'Spatial Analyst'
extension. Firstly, you will need to start ArcView, start a new project, add a view and set the
working directory (to 'c:\temp').
ArcView Steps
Step 1 GeoProcessing - Dissolve
We need to load the GeoProcessing extension - go to the 'File' menu and then 'Extensions'.
Put a tick in the box corresponding to 'GeoProcessing' and then click on 'OK'. We mw require
some sample dala to experiment with. Add the following theme from the directory
'c:\arcv32\arcview\esridata\usa\ ': 'counties.shp '. You may wish to hide the legend within the view
(with the relevant theme active, go to the 'Theme' menu and select 'Hide/show Legend'.
We are now going to aggregate the counties boundaries into states. This can be achieved
because each county has a corresponding state name within the attribute table. Go to the
'GeoProcessing Wizard' within the 'View' menu. Because we've presently only got one theme
open all but the dissolve GeoProcessing option are currently unavailable. Click on the 'Next'
burton and then dissolve using the attribute of 'State name'. Save the theme to the following
location: 'c:\temp\us state.;' (as seen below), and then click on the 'Next' button. We can now
choose the attribute data we wish to keep -highlight 'Area by sum' and finally click on the 'Finish'
button. You should mw have a theme that contains the US state boundaries.
Practice 5: Advanced Analysis I ERS 120: Principles of GIS
N.D. Bình 27/59
Step 2 GeoProcessing - Merge
We are now going to combine the Canadian Provinces with our US state data. With the
'merge' function only attributes with identical theme names are kept (based only on one of the
themes attribute titles). Therefore, we need to rename the attribute titles of our US state theme.
Go to the US states table and then select 'Properties' from the 'Table' menu. Change the
field aliases to the following:
Once you are finished dick on 'OK'. This will only save the attribute tables as we see them;
the flies themselves are still really entitled using the original tags. To overcome this we need to
dose the table and then select 'Convert to Shapefile' under the 'Theme' menu. Save your theme
as c:\Itemp\us_states2'. We can now discard the old flies, to do this make both the 'us_states'
and 'counties.shp' active (you can select multiple themes by holding down the shift button when
activating your themes), go to the 'Edit' menu and select 'Delete Themes'.
Open the Canadian Provinces theme ‘c:\arcv32\arcview\esridata\canada\province.shp', go to
the Geoprocessing wizard (under the 'View' menu) and select 'merge'. Highlight the required
themes and select the options as illustrated below:
You should now have both the US states and Canadian Provinces in the same theme.
Examine the new theme using the 'identify' tool .Remember that the identity
tool only works for the active theme.
Step 3 GeoProcessing - Clip
The 'clip' function allows us to cut a theme based on another themes extent - so that we keep
only the required study area from a larger theme. Open the river theme for the US
('c:\arcv32\arcview\esridata\usa\rivers.shp j. Now select any four US adjoining states that all
contain a river using the 'Select feature' tool . Do not unselect the selected features until
Practice 5: Advanced Analysis I ERS 120: Principles of GIS
N.D. Bình 28/59
informed to do so. For an example, see below:
Go to the GeoProcessing Wizard and select' Clip one theme based on another' and clip the
river theme based on the extent of the selected features within the US states theme - as
demonstrated below:
A new theme will be created containing only the rivers from the selected four states. We can
now delete the un-required themes (River_Sample.shp, N_America.shp & Province.shp). To do
this, make the themes active and select 'Delete Themes' from the 'Edit' menu.
Step 4 GeoProcessing - Intersect
The 'intersect' function allows the user to divide a theme up into new segments based on
features from the overlay theme. The output theme will contain attributes from both themes. Open
up the 'c:\arcv32\arcview\esridata\usaladi.shp’ which represents the overlay theme. Now select
the 'intersect two themes' option from the GeoProcessing wizard and select the following:
Practice 5: Advanced Analysis I ERS 120: Principles of GIS
N.D. Bình 29/59
Notice how we now have both boundaries combined for the selected states as illustrated in
the figure below. Once finished you can delete the following themes: 'adi.shp' and
'Adi_States.shp'.
Step 5 GeoProcessing - Union
The 'union' function is different to those previously used for it produces an output theme that
contains the attributes and full extent of both themes. Open up the theme containing data on the
major US lakes ('c:\arcv32\arcview\esridata\usa\lakes.shp') and clear the selected features for the
US states theme - when the theme is active). Now go to the GeoProcessing wizard, select
the 'union two themes' function and fill out the criteria as depicted below: -
Practice 5: Advanced Analysis I ERS 120: Principles of GIS
N.D. Bình 30/59
You can now delete both the 'lakes. slip' and ' lake union. slip, themes.
Step 6 GeoProcessing - Spatial Join
Last session you saw how to perform manual geographic spatial joins, the option 'Assign data
by location' is a more automated approach. Open up the following theme:
't:\arcv32\arcview\esridata\usa\cities.shp' containing information pertaining to US cities. We will
row join both themes so that the relevant state information is available for each city.
Sellect the 'Assign data by location' from the GeoProcessing wizard and assign data to the
'cities. slip' ,from 'us states2.slIp'. "When you enter the cities attribute table some new fields will
be added to the right-hand side of the table.
Step 7 Spatial Analyst - Find Distance
Within this section the main functiouality of the spatial analysis extension will be explored
Firstly, we need to load the relevant extension ('Spatial Analyst'). You should see that a new
menu is now visible entitled' Analysis' .We will row calculate the distance from the ten largest
cities. We need to enter the table for the 'cities.shp' ( and sort the field entitled 'Pop1990'
descending ( . Now select the first ten records by clicking on them while holding the shift key -
remember they should turn yellow to indicate they are selected. Close the table to return to the
view.
Go to the 'Analysis' menu and select 'Find Distance' - enter the grid extent to be the same as
that of the 'US-states2.shp' theme. Now we are required to inform Arcview of the grid resolution
we require - obviously a greater resolution will be more accurate but will incur longer calculation
times. Enter the number of rows to be 200 and press return - the grid cell size and the number of
columns should be automatically re-calculated. The resulting raster grid should resemble the
following:
Practice 5: Advanced Analysis I ERS 120: Principles of GIS
N.D. Bình 31/59
Each cell now represents the distance from the nearest of the ten selected cities, in miles.
Now repeat the process for the 'rivers.shp' theme, using a grid extent to be the same as that ofthe
'US-states2.shp' theme and the number of rows to equal 200. Your results should be as follows:
Step 8 Spatial Analyst - Assign Proximity
We can also determine the catchment area for each of the selected cities using the 'assign
proximity' function. This assigns the identifier field of the nearest feature to each cell. With the
'cities/shp' theme active, go to 'Analysis' and then 'Assign Proximity'. Again enter the grid extent
to be the same as that of the 'US-states2.shp' theme and the number of rows to equal 200. Now
we need to select the field that will be used to differentiate between the selected features -
usually the feature ill, in this case we need to highlight 'City_fips'. The result should resemble:
It is often useful to convert the proximity from its original grid format to vector polygons. If you
require this you should convert the theme into a shapefile ('Convert to shapefile' under the
'Theme' menu).
Step 9 Spatial Analyst - Density
Densities can be calculated using the 'Spatial Analyst' extension. Clear the selected features
of the 'cities.shp' theme ( ) and then select 'Calculate Density' from the 'Analysis' menu. Again
enter the grid extent to be the same as that of the 'us_states2.shp' theme and the number of rows
to equal 200. Now select the following options:
Practice 5: Advanced Analysis I ERS 120: Principles of GIS
N.D. Bình 32/59
Your result should resemble that illustrated those over the page:
Step 10 Spatial Analyst - Summarise Zones
This function calculates summary attributes for features using a grid theme. Of the most useful
applications is to calculate summary statistics of a grid that relates topolygon boundaries
contained within another theme.
Make sure the boundary layer is active (US_states2.shp) and then select 'Summarise Zones'
from the 'Analysis' menu. We now have to identify which field will be used to uniquely identify the
state polygons - only 'name' can be used so click on OK.
Then we must choose the required grid - select 'densities from cities' and click OK. Seeming
as there are more than 25 states ArcView will not draw a chart of selected statistics. A new table
should appear illustrating summary information for each state.
We can now examine the number of grid cells within each state ('count') but also the
minimum, maximum, range, mean, standard deviation and sum of the grid cells, broken down by
the state name. You can now delete the following themes: 'density from cities.shp': 'proximity
form cities.shp"; 'rivers.shp'; and 'cities.shp':
Step 11 Spatial Analyst - Reclassify
Reclassify can be used to classify grid cells into explicit
groupings. For the remainder Of this practical we will use the
'distance from cities' and 'distance from rivers.shp' for site
analysis. We need to find a suitable location that is both close to
one of the ten largest US cities but also away from major rivers.
We will use 'reclassify' to assign a suitability rating based on
the distance grids. Make the 'Distance to cities' grid active and
then go to 'reclassify' within the 'Analysis' menu. Values close to
the cities will be given a low score to represent that these
locations are most suitable. Fill in the 'reclassify values' table to
hold the following information:
The result should resemble the following:
Practice 5: Advanced Analysis I ERS 120: Principles of GIS
N.D. Bình 33/59
Now repeat the process for the 'Distance from rivers' grid, but reverse the rating given for the
new values (i.e. 0 - 0.5 = 9 and 10 - 100 = 1). By assigning distances close to rivers with a high
rating we are indicating that these sites are undesirable locations. Your resulting grid should look
like the one below:
Step 12 Spatial Analyst - Map Calculator
We can manipulate raster grids using the 'map calculator' - here we want to add to grids
together. Go to the 'Analysis' menu and then the 'Map Calculator', now double click on 'Reclass of
Distance to Rivers.shp' (top layer), single click on the '+' operator and then double click on
'Reclass of Distance to Cities.shp' (third layer down). Your calculation should read: -
[Reclass of Distance to Cities.shpJ + [Reclass of Distance to Cities.shpJ
Click on the 'evaluate' button and then close the map calculator window. The resulting grid is a
combination of the values from both grids. The lower the value the more acceptable the site will
be, the lowest rating appears to be 6 - representing the most optimal locations.
Step 13 Spatial Analyst - Map Query
It is possible to select all the cells with certain characteristics - in this example where the rating
equals '6'. With 'Map Calculation l ' active go to the 'Analysis' theme and select the 'Map Query'
function. Enter the following query, by double clicking on 'Map Calculation L" (top layer), then
single clicking on the '=' operator and finally double clicking on the unique value of '6': -
[Map Calculation] = 6.AsGrid
Click on the 'evaluate' button and then close the 'map query' window. A new Boolean theme
will be produced - indicating those cells that fulfill the map query criteria (grid cell will become
'true' = 1) and those that do not ('false' = 0). Your final map of the most suitable locations should
resemble the figure illustrated below:
Practice 5: Advanced Analysis I ERS 120: Principles of GIS
N.D. Bình 34/59
You can now close ArcView. There is no need to save any files - they are no longer required.
All files should have been saved to the temp directory and as such will automatically be deleted
when the computer is re-started.
Last modified: Oct 25, 2009
ERS 120: Introduction to Geographic Information Systems /
Các file đính kèm theo tài liệu này:
- prac05_analysisi_8895_4575_2053323.pdf