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: emp).
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. A ll 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 / . 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. 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. 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