WEEK 10_GIS3015 - FLOW MAPPING

This weeks topic was Flow Mapping.

Process Summary: 

1. Navigate to s:drive.  Create a folder for the week; Week 10_Flow Mapping.
2. Navigate to r:drive and download this weeks shapefiles, 2 excel documents and the pdf file.
3. Paste to Week 10_Flow Mapping and unzip world.zip and states.zip.
4. Open ArcCatalog and Minimize.  Navigate to the Week 10 folder. Bring up the files and view properties of the two shapefiles – No Projection.
5. Open ArcMap Window and minimize.  Drag the two shapefiles into ArcMap from ArcCatolog.
6. Adjust the Map before converting to illustrator. Change to layout view. Change to landscape. Insert Neatline – Inside Margins with a Gap of 10, Border 2.0, Background – Blue, Drop Shadow – Grey 60%.
7. Save:  Week 10 Map_AK.mxd
8. Add Scale Bar: Alternative 2 – miles.  Title “2007 Immigration to the United States”. Author, date and Source (2007 INS Statistical Yearbook).  Legend – allow for 8 different legends for each of the 8 regions to Map (USA, Africa, Asia, Europe, North America, South America, Oceania and Uknown/Australia) – additionally for legend border = 1.0 point, background = grey 40%, drop shadow = grey 60%, Gap = 10, Remove the term “Legend” prior to finishing inserting the legend.
9. Save and then Export File as Week 10 Map_AK.ai under the week 10 folder. Exit out of ArcMap and open iIlustrator from the start menu.  When you receive error message be sure to choose “update”.
10. Insert 4 new layers and title: Map Title, North Arrow, Legend and Author.  Move the respective items to each from the “Other 3” layer.  Rename the Other 3 layer to Scale once have moved all other items out of it.  Add a final layer titled “Flow Lines”.  Lock all the layers that you are done with and do not want to alter.
11. In Excel, calculate the width of your line using the formula; width of line symbol = (maximum line width = 10 here) X (Value/Maximum Value).  Then insert corresponding stroke points to be able to reflect the line widths appropriately on the map Multiply the calculated Flow line width by “3” points.
12. Using the pencil tool draw the lines.  Add arrows to the end pointing at the United states using the Arroheads option from the stroke menu – Scale 30%.  Match arrow color to that of the region.  Add the text to the arrow = population migrating into the united states from the shown region.  Example Details on Arrowheads:
13. Asia – Sroke point 30 – right arrowhead, scale 30, color to match.
14. Europe – Sroke point 17 – right arrowhead, scale 20, color to match.
15. Once done drawing the arrow heads, insert the populations and align on the flow line.
16. Lock all layers.  Save and Export file as a Jpeg.

           

GIS 4043_WEEK 10_DATA EDITING

This week’s goal was to obtain experience in editing feature classes and shapefiles, digitizing features and snapping.

PROCESS SUMMARY:

1.  Log into e_GIS and create folder for the week – Week 10_Data Editing\Work

2.  Copy “Dataediting.zip” on r:\drive to S:\Intro2GIS\Week 10_Data Editing\Work and unzip file.

3.  Open up ArcMap and minimize window.

4.  With ArcCatalog opened and window minimized, navigate to folder with the shape files, GulfCounty_Roads and GulfCounty_Parcels.  Drag both files to ArcMap window to add the two layers.

5.  Maximize ArcMap and click on Editing toolbar, Editor and Start Editing.

6.  With GulfCounty_Parcels layer highlighted, use the Select by Attributes option to select the two recently purchased parcels.  Final “where Clause” should read “PARCELNO”=’03334-000R’Or”PARCELNO”=’03341-000R. Verify the expression then Click Ok to select both parcels.

7.  Zoom to the selected parcels, click on Editor, select Merge, then select GulfCounty_Parcels-RISH WILLIAM J and hit Ok to complete.

8.  Open the GulfCounty_Parcels attribute table, click on Show Selected Records, and make changes as shown: OWN_NAME = Gulf County, OWN_ADDRESS1 = 1000 5^th St., OWN_CITY = Port St. Joe, and OWN_STATE stayed FL.

9.  Close attribute table, click on Editor, Save Edits, and Stop Editing.  Save work.

10.  To find suitable parcels for the extension office site, run new select by attributes query  “OWN_NAME”=’GULF COUNTY’.  74 records are selected.  Right-click GulfCounty_Parcels, select Data, then Export Data as a new shapefile layer titled,  County_Owned.shp

11.  Open attribute table for County_Owned, select Add Field, Name it Acres, select the type Double, click Ok, ANOTHER WAY TO DO THIS:

Close out ArcMap. Go to ArcCatalog and right click on shapefile. Under properties select fields and then add new field as specified above. Once done, reopen arcmap to proceed.

12.  Right-click column title of new field and select calculate geometry, Select Units - Acres US and hit Ok. Right-click again and select Field Calculator, and typed in the expression round([Acres],2) to round to two decimal places.

13.  Open properties of County_Owned, select the Definition Query tab, type in the expression, Acres >=20 to remove parcels less than 20 acres and click Ok.  Only 10 parcels are selected. 

14.  Click Add Data, and add  VICD.dbf. Right-click the County_Owned layer, select Join, then Join Attributes from a Table, with the Parcel_ID.

15.  To get the available vacant parcels, Select by Attributes of the County_Owned layer with the expression, “VICD.VI_CD1”=’V’, this gave the available parcels that were vacant. Export the selected parcels as a new layer called SuitableParcels.shp.

16.  To adjust the roads, Click on Start Editing, select GulfCounty_Road, and select by attribute, roads FID 4168, with expression “FID”=4168.

17.  Double-click on the selected road, turn off vertex and edge snapping and drag the road to the center of it proper location in the middle of the roadway.  May take some playing  around with using the “snapping option” until you get the road in the center.  Save Edits. 

18.  Click the create Features Button in the Edit toolbar, double-click on GulfCounty_Roads, click Ok, then click on the new road to draw the new street with vertices.

19.  Click on the attributes for this new road and fill in the new attributes as shown on point #7 of pg. 9 of the lab summary.  Click Save Edits and stopping editing.

20.  Add all essential Map elements; Change to layout view, landscape, insert neatline, author/date and source, scale bar(be sure to adjust units to feet to match shapefiles), legend.

21.  Click to Add a Basemap, Bing Maps Aerial.

22.  Adjust the suitable parcels layer to have no fill color and red outline with a width of 2. 

23.  Adjust the Gulf County Parcels layer also with no fill color and a black outline with a 1.5 outline.  Make the roads green in color with a width size of 2.

24.  Insert label in the suitable Parcels layer by clicking on the open attribute table and under, “OWN_NAME” type in “03334_000R”.

25.  Save map and export as SuitableParcels_Gulf CountyBOCC_Ext.jpg

WEEK 9: GIS 3015_ISOPLETH MAPPING

GOAL:  Create an isohyet map illustrating mean annual precipitation for the state of Georgia.

PROCESS SUMMARY DESCRIPTION

Copy from R:drive and paste to S:drive this week’s file (Georgia.ai)>GOAL:  Create an isohyet map illustrating mean annual precipitation for the state of Georgia>b.  Include all map elements>c.  Use the Pen and Pencil Tools>d.  Know how to Reshape a Path Outline>e.  isohyet interval = 5 inches, starting value = 0>f.  Data Source = Southeast Regional Climate Center and National Oceanic and Atmospheric Administration> THE LAB: - Open up Adobe Illustrator and opened file (Georgia.ai)>Printed out hard copy of the file and sketched contours by hand as well as determined where labels would go>Back to Adobe Illustrator: Open Layers Window and add a new layer called Contours (NOTE: Location to save all contours and labels)>Begin with the contours and Begin to carefully use the Pencil tool to draw the contours which should be black with a width of 2 points> When drawing the contours with the pencil tool, be sure to break before a turn so that if you have to erase, it is not the entire drawn contour but just until the last lift of the pencil (NOTE:  To delete, Edit>Undo)>As you complete each contour then name the group (double click on the sub-layer and type in name under the name box) so you can come back later and label it>NOTE:  Countours got really close at close to 0 inches>Save work periodically>Next create a new Layer called “Labels”>Begin to label each contour, midway and any other location on the contour that will help the reader know what contour they are looking at>TO LABEL: >a.  Click on Text Tool and click into artboard>Type in the Contour label, ex. 40>Adjust font to myriad pro – size 6>Using the rectangle tool, draw a rectangle around this>Change fill color to same color as the map background>under the layers window – Labs Layer, drag the  Path layer to be under the number layer, so that the number can be seen on the contour map>Using the direct tool, click to  highlight the “Path Layer” and while holding down the shift key, click to highlight the “number” layer> Select the Object Menu>Group (NOTE:  This groups these two layers into one, so that when you move, it moves as one)>Using the direct select tool, move the entire group and position strategically over the contour line to where you are not blocking much any of the other features, precipitation, Rainfall amount>Repeat process for all contours (NOTE:  Due to how close the contour lines are for the 0 to 40, at 5 inches interval, it was decided not to label these on the actual map to avoid crowding)>Once done move the Labels layer directly underneath the Contours Group>Click to Lock the Layer>Click to lock all the other layers too (Layers, Scale, Precipitation, Rain Amount and Contours) as they are not being used>Add a new layer can call this Neatline>Using the Rectangle box tool draw a rectangle around the map and remove the fill>Drag to move the neatline to the most bottom layer in the TOC under layers>This allows all other map features to be shown on the map>Make a new layer for the remaining items as above and again add the respective features; North Arrow, Source/date and author, Title and Legend (NOTE:  For the legend use symbols, text and pencil tools to complete then put a boundary box around it)>Save File>Before Exporting be sure to check that the file is a Black and White file as requested: >Unlock and click the layers layer to highlight>Using drop down menu remove the fill>Save and lock that layer>b. Be sure to remove the fill and stroke from the contour labels too.  However here we are not really removing the fill but rather making it the same color as the background = white>Click on fill dropdown menu>Select Swatch Options Icon at bottom of window>Fill in the information as shown: a. Swatch Name: K=0, >b.  Color Type=Process Color>c.  Color Mode = Grayscale (NOTE:  Drag the scale until it reads 0% to get it to white and where it matches the background>OK>Be sure to check that there is no Color selected for the stroke>Complete for all labels in the same manner>i.e. click on the “path” layer>go to the drop down menu for the fill>select the white fill to match map background>select the no fill on the stroke>Repeat one by one for all labels(NOTE:  unfortunately at this point selecting globally will also eliminate the labels from the contours so you have to be sure to only select the path and the various groups> Export File as a Jpeg file>Open up the Jpeg file and adjust any issues with margins>Save in the project folder ready to post to blog>Archive.

WEEK 9: GIS 4043 VECTOR ANALYSIS II

PROCESS SUMMARY DESCRIPTION

NOTE:  A little more detailed, simply because it helps highlight every single step taken and makes it easier to reference areas that may have been complex when first done.

Open ArcMap and minimize the window>download this week’s file (NOTE: Vector II.zip)>unzip>Open ArcCatalog 10 and minimize the window then drag “Roads.shp” to Blank ArcMap Template>Maximize window to be able to see ArcTool Box>ArcToolbox: >a. Analysis Tools>b.  Proximity>c. Buffer(OBSERVE: Buffer window comes up)>input options as shown: >a. Input Features = Roads>b.  Output Feature Class = S:\Intro2GIS\WEEK 9_VECTOR ANALYSIS II\Roads_Buffer.shp (NOTE:  Record the output under this weeks file as “Roads_Buffer”)>c.  Distance (value or Field) = check the Linear Unit box and enter distance as 300 meters>d.  Side Type = FULL>e.  End Type=ROUND>f.  Dissolve Type = ALL>OK (OBSERVE:  A message comes up at the bottom right of the screen saying “buffer”.  When done a green checkmark is shown at the same location at the bottom of the screen and a new layer [Roads_Buffer] appears on your map.  Also a buffer feature is now around the initial roads on the map.)>Turn off the Roads_Buffer layer under TOC>Save File (Vector Analysis II_AK.mxd under this week’s folder) and move on to step 2: CREATING A VARIABLE DISTANCE BUFFER: >To ArcCatalog Minimized Window: drag into your ArcMap the shapefile “Water.shp”>Save and Close File in ArcMap in order to be able to update it’s attribute table in ArcCatalog10 (NOTE:  Will be unable to change attribute table when file is open in ArcMap10, or any file for that matter is open in ArcMap10.  Completely exit from all instances of ArcMap10>ArcCatalog10>under contents tab, right click on file water.shp>properties>scroll down to empty field and click on the row to highlight Field Name and Data Type (NOTE:  Data Type will be highlighted in Black while the first field, “Field Name” will be highlighted in white)>Start typing (NOTE:  Will begin to type in the first box = Field Name Box>type, “BUFFDIST”>Tab to Data Type Column (NOTE: Drop down menu option pops up) – select “short integer”>(NOTE:  To delete a field name and data type if you make a mistake or do not want it, simply highlight the row and hit delete on keyboard)>Reopen file in ArcMap 10>Open  Attribute table (OBSERVE: the new attribute is automatically added to the properties for the water layer.  Also notice that under Type have both lakes and rivers (13 lakes and 1 river)>TO ASSIGN DIFFERENT VALUES TO LAKES AND RIVERS:>Under attributes table – click select by attributes and input the following information: >a.  Method = Create a new selection (NOTE:  You can also select from current selection but before doing this you will have preselected  ALL rows in the attribute table – select first row>hold shift key down and scroll to last row and click to select – highlights all columns)> b.  Double click on type to drop it in the “select from water where:” box>c.  click on the “=” sign>d. Click on “Get Unique Values” button>e.  Double click on ‘Lake’ to drop it in the “select from water where:” box>f.  verify (NOTE:  A message comes up “the expression was successfully verified”>OK>Apply(NOTE:  All 13 lakes are selected – the number is displayed under the (13 out of 14 selected)>scroll to the buffdist column and right click on column heading>select “Field Calculator Function (OBSERVE:  A Field Calculator window pops up)>Enter the following information in the Field Calculator window: >a.  Parser = VB Script>b.  Fields = Leave as is>c.  Show Codeblock = leave as is>d.  Type = see that Number is checked>e.  Functions = Leave as is>buffdist = box = type in 150>OK(OBSERVE:  Get a Processing window with the message “calculating records”.  Once done should see all the 13 lakes assigned the value 150 under the buffdist column>Repeat for the river (NOTE:  To deselect under the attribute table, click on the far top left box at the start of the column headings)>(NOTE: when the select by attributes window comes up be sure to delete the previous expression from the “select from water where:” box and select river instead of lakes this time. Should see that one selection is made>Assign the value 500 as done above with 150>Now that both have been assigned values, deselect attributes table as shown above (NOTE:  Be sure to do this so that the buffering is applied to all and not just where you left selected)>Close  the Attributes table and open ArcTool Box to the Buffer window as indicated above>Under the Buffer window, enter the shown information:>a.  Input Features = Water> b.  Output Feature Class = S:\Intro2GIS\WEEK 9_VECTOR ANALYSIS II\Water_Buffer.shp (NOTE:  Record the output under this weeks file as “Water_Buffer”)>c.  Distance (value or Field) = check the Field  box and select buffdist>d.  Side Type = FULL>e.  End Type=ROUND (NOTE:  This is greyed out)>f.  Dissolve Type = LIST>g. Dissolve Field(s) = BUFFDIST>OK/Apply (OBSERVE:  A message comes up at the bottom right of the screen saying “buffer”.  When done a green checkmark is shown at the same location at the bottom of the screen and a new layer [Water_Buffer] appears on your map.  Also a buffer feature is now shown on the map.)>Bring the water layer above the buffer layer under TOC (NOTE:  If having issues doing this, through the drag feature, save the initial file and open another, then remove both layers>Then add them back onto the map, but add the Water_Buffer layer first and then the Water layer (OBSERVE:  it now shows up to where you can see both layers and the water layer is not covered up as happened before.)>To make the water layer more visible, change the color to a dark blue including the outline >Turn the Roads and Roads_Buffer Layers on to be able to see all layers at once(NOTE:  In arranging your layers, it may be best to arrange such that the roads is the top layer, followed as shown – Roads>b. Roads Buffer>c. Water>d. Water_Buffer >Save File (Vector Analysis II_AK.mxd under this week’s folder) and move on to step 3.  USING ArcPython TO RUN BUFFER ANALYSIS:  Maximize the ArcMap window to be able to see the ArcPy menu button>before clicking on arcpy menu button, go to the Geoprocessing menu>Geoprocessing Options>Check that the first box is checked “Overwrite the outputs of geoprocessing operations”>OK to close window>click on arcpy menu button (OBSERVE:  Window called Python, comes up)>Type  the first line of code that you need to write when working with ArcPy: “import arcpy” (NOTE:  space between the two but no quotation marks.  If a selection box comes up as you type, double click on the choice of word to add to the type line>Enter key to move to next line>Type “arcpy.(include  the period, but no quotation marks)(NOTE:  drop menu box comes up and this narrows down as you continue to type.  If does not then you may have made a typing error, check your spelling)>Type Buffer_analysis (NOTE:  can  type “buffer” then double click on buffer_analysis to select when it comes up in the drop down menu  box)>Type an Open Parenthesis (NOTE: A window automatically pops up giving you a list of possible “in_features”.   Review notes on commands and using open parenthesis pg. 8 of lab instructions)>Double click on “Roads”>type a comma to move to the out_feature_class>type opening quotation marks>copy and paste in path of Roads_buffer (leave out the ” .shp” part of the path, as it will automatically create a shape file)  file from process summary>type closing quotation marks>type comma >Open quotation marks> type the distance 100>space>type Meters (NOTE:  with initial capital)>closing quotation marks>comma>a selection window pops up, select “FULL” (NOTE:  comes already with the quotation marks so no need to type those>comma (OBSERVE: a second pop up window opens)>select “ROUND”>comma(OBSERVE: another pop up window opens)>select “ALL”>type closing parenthesis (NOTE:  quotation marks and comma, no space [except for the distance] are key in your python commands)>Press enter key(OBSERVE: words Buffer appear on bottom right hand window and message appears to show your results having been successfully added, see message below,

>>>import arcpy

>>> arcpy.Buffer_analysis("Roads",”S:\Intro2GIS\WEEK 9_VECTOR ANALYSIS II\Roads_Buffer”,”100 Meters”,"FULL","ROUND","ALL")

<Result ‘S:\Intro2GIS\WEEK 9_VECTOR ANALYSIS II\Roads_Buffer.shp’>

>>>

PART 2: OVERLAY IN ArcMap

Save and Close File in ArcMap in order to be able to update it’s attribute table in ArcCatalog10 (NOTE:  Will be unable to change attribute table when file is open in ArcMap10, or any file for that matter is open in ArcMap10.  Completely exit from all instances of ArcMap10>ArcCatalog10>under contents tab, right click on file water_buffer.shp>properties>scroll down to empty field and click on the row to highlight Field Name and Data Type (NOTE:  Data Type will be highlighted in Black while the first field, “Field Name” will be highlighted in white)>Start typing (NOTE:  Will begin to type in the first box = Field Name Box>type, “INSD_WBUF” (NOTE: Skip the period as shown in the lab)>Tab to Data Type Column (NOTE: Drop down menu option pops up) – select “short integer”>Reopen file in ArcMap 10>Open  Attribute table (OBSERVE: the new attribute is automatically added to the properties for the water_buffer  layer.  >TO ASSIGN A  VALUE FOR EVERY FEATURE IN THE LAYER:>Skip the select by attribute and scroll to the INSD_WBUF column and right click on column heading>select “Field Calculator Function (OBSERVE:  A Field Calculator window pops up)>Enter the following information in the Field Calculator window: >a.  Parser = VB Script>b.  Fields = Leave as is>c.  Show Codeblock = leave as is>d.  Type = see that Number is checked>e.  Functions = Leave as is>INSD_WBUF = box = type in 1>OK(OBSERVE:  Get a Processing window with the message “calculating records”.   Once done should see all features assigned the value 1 under the INSD_WBUF column>Close Attributes Table and Save file>Repeat for Roads_Buffer layer using 1 for the value and INSD_RBUF for the new attribute>Close  Attributes table and Save file>Open ArcToolbox>Analysis Tools>Overlay>Union (OBSERVE:  A union window pops up.  Follow the instructions to complete):> a. Input Features = One by one, click on Roads_Buffer then Water_Buffer to add to the Input Features Selection>b.  Output Feature Class = S:\Intro2GIS\WEEK 9_VECTOR ANALYSIS II\Buffers_Union.shp (Name the file “Buffers_Union”)>c.  Join Attributes= ALL>d.  Gaps Allowed = Check box>e.  XY Tolerance = leave blank and leave units for meter>OK (OBSERVE:  A message comes up at the bottom right of the screen saying “union”.  When done a green checkmark is shown at the same location at the bottom of the screen and a new layer [Buffers_Union] appears on your map.  Also a buffer feature is now around the initial roads on the map.)>Save File (Vector Analysis II_AK.mxd under this week’s folder) and move on to the next step:  SELECTING RECORDS BY ATTRIBUTE: Under attributes table – click select by attributes and input the following information: >a.  Method = Create a new selection (NOTE:  You can also select from current selection but before doing this you will have preselected  ALL rows in the attribute table – select first row>hold shift key down and scroll to last row and click to select – highlights all columns)> b.  Double click on INSD_WBUF  to drop it in the “select from Buffer_Union where:” box>c.  click on the “=” sign>d. Click on “Get Unique Values” button>e.  Double click on ‘1’ to drop it in the “select from Buffer_Union where:” box>Add the 2^nd query using the “And” button>after clicking the and button repeat the above for INSD_RBUF with a valueof 1 too>f.  verify (NOTE:  A message comes up “the expression was successfully verified”>OK>Apply(NOTE:  only the records that meet the shown query show up)>EXPORT THE SELECTED FEATURES TO A NEW SHAPEFILE USING THE SHOWN INSTRUCTIONS: Right click on Buffers_Union>Data>Export data (OBSERVE: An export data window pops up)>a.  Export = selected features>b. Use the same coordinate system as = check this layer’s source data>c.  Output Feature Class = Path to this week’s folder and save as Buffer_Union_Export>Yes to add new layer to shown map>New layer with the shown name shows up>uncheck all layers but for Roads, Water and Buffer_Union Export>Repeat using the Update Overlay tool (NOTE:  Results are the same)>Move to the final step: SEPARATING AND EXCLUDING AREAS: With ArcMap and ArcCatalog windows both minimized, drag into ArcMap the conservation areas shapefile (NOTE: New layer = Conservation Areas, shows up under TOC and the areas are shown on map)>can modify color – change to green (best for showing conservation areas – personal preference)>Using overlay tool exclude the conservation areas that need protection of plant and animals (NOTE:  Use the Erase tool)>a. Input Feature = Buffer_Union_Export>b.  Erase Features = Conservation_Areas> Output Feature Class = S:\Intro2GIS\WEEK 9_VECTOR ANALYSIS II\Possible_Sites.shp (NOTE:  Record the output under this weeks file as “Possible_Sites”)>OK (OBSERVE: The word Erase scrolls in bottom right of window and a new layer shows up on map called Possible_Sites>NEXT – TURN MULTIPLE PART LAYERS INTO SINGLE PART LAYERS: Right click on the Possible_Sites layer>ArcToolbox>Data Management Tools>Features>Multipart To Singlepart (OBSERVE: A multipart to singlepart window pops up)>Input information as shown: >a.  Input Features = Possible Sites>b.  Output Feature Class = S:\Intro2GIS\WEEK 9_VECTOR ANALYSIS II\Possible_Sites_Single.shp (NOTE:  Record the output under this weeks file as “Possible_Sites_Single”)>OK (OBSERVE: The word Multipart to singlepart scrolls in bottom right of window and a new layer shows up on map called Possible_Sites once the process is done>Add a new attribute to the Possible_Sites_Single table>Close all sessions of ArcMap and be sure you only have on session of ArcCatalog running>From ArcCatolog insert as shown above, a new field “Area” with a type “Double”>Once inserted and back in ArcMap open up attributes table for this layer>select all records>right click on the area column>select calculate geometry (OBSERVE:  A window opens up)>populate window as shown>a.  Property = Area>b.  Coordinate system = use one specified>c.  Units = square Meters (NOTE:  Given that meters has been the linear unit for the shapefiles, for Area we need Square meters>d.  Calculate selected records only = check>OK (OBSERVE:  Window pops up saying calculating records)>once done, the areas of each of the areas is added to the column>To find minimum and maximum use the sort ascending and sort descending features respectively (NOTE:  Be sure to drag the column to allow for the entire width of the column showing all the areas)>deselect the fields>add all the map elements remembering that for the scale and scale text, you should use units>the legend will be made of only 3 items; roads, water and possible_sites_single>Save file>Export as a Jpeg File and lastly archive.

SEE FINAL MAP BELOW

WEEK 8 - GIS 4043 Data Analysis/Vector Analysis 1

INTRODUCTION TO SPATIAL ANALYSIS

The following is this week's project summary:

INSTRUCTIONS ON HOW TO FIND NEARBY FEATURES:

1)      Query for Criteria #1: Land Parcels within 2000 feet of a school - (Begin by setting up the query to chose features from the target layer) Open ArcMap>Open up VectorAnalysis_1.mxd file>Selection>Select by Location>Selection Method>Select Features From (Because you have this "select feature from" it is prompting you to identify a target layer so have to do so in the next step)>  (NEXT - identify your target layer)Target Layer(s) - check box next to Pensacola Parcels>(Select your Source Layer = layer that you use to relate to your target layer )Source Layer - choose Pensacola Schools>(Set up a Spatial Relationship)Spatial Selection Method>Target layer(s) features are within a distance of the source layer feature (avoid the 3D option since we are using 2D for this lab)>Check the "Apply a search distance" box>Type "2000" in the "Apply a search distance" box AND choose "feet" in the dropdown menu>OK

2)      Query for Criteria #2:  Land Parcels within one mile of a hospital - with ArcMap still open>Selection>Select by Location>Selection Method>Select from the currently selected features(NOTE:  Here you are omitting the step to identify your target layer coz you are selecting from the currently selected features (SO SAYING TARGET LAYER = Pensacola Parcels selected in initial Query and just carried on)>Source Layer - choose Escambia_hospitals> Spatial Selection Method>Target layer(s) features are within a distance of the source layer feature (avoid the 3D option since we are using 2D for this lab)>Check the "Apply a search distance" box>Type "1" in the "Apply a search distance" box AND choose "miles" in the dropdown menu>APPLY(not Ok, to apply the query and keep the window open to perform the next query)

3)      Query for Criteria #3:  Remove Features that are within 1 Kilometer of a Railroad - with ArcMap still open>Selection>Select by Location>Selection Method>Remove from the currently selected features in(NOTE:  Here you are omitting the step to identify your target layer coz you are REMOVING from the currently selected features (SO SAYING TARGET LAYER = Pensacola Parcels selected in initial Query and just carried on)>Source Layer - choose Rails >Spatial Selection Method>Target layer(s) features are within a distance of the source layer feature (avoid the 3D option since we are using 2D for this lab)>Check the "Apply a search distance" box>Type "1" in the "Apply a search distance" box AND choose "kilometers" in the dropdown menu> APPLY(not Ok, to apply the query and keep the window open to perform the next query)

4)      Query for Criteria #4:  Land Parcels within 3000 feet of a cellular tower  - with ArcMap still open>Selection>Select by Location>Selection Method>Select from the currently selected features >Source Layer - choose Escambia_Cellular >(Set up a Spatial Relationship)Spatial Selection Method>Target layer(s) features are within a distance of the source layer feature (avoid the 3D option since we are using 2D for this lab)>Check the "Apply a search distance" box>Type "3000" in the "Apply a search distance" box AND choose "feet" in the dropdown menu>OK(Since we are done and not applying any more queries)

Once you have found these features - see how many have been selected - Instructions:

5)      Right click on the Pensacola Parcels layer in the TOC>Open Attribute Table>see text at bottom of attribute table window telling you how many of the 27547 features in this layer have been selected (DO NOT CLICK ON ANYTHING in the Attribute table window - Could change your selection)>close when done

Preserving your work - Instructions:

6)      Rightclick on the Pensacola Parcels layer in the TOC>Selection>Create Layer from selected features (A new layer is automatically created called "Pensacola Parcels Selection")>(Rename this layer to "Parcels Nearby")>Save your map>Uncheck your Parcels Nearby layer in the TOC to hide it from display>Save Map

If have to return to the Parcels Nearby Layer coz something goes wrong along the way - Instructions:

7)      Query to retrieve Nearby Parcels Layer/Query to retrieve Land Parcels within the Nearby Parcels Layer:   with ArcMap still open>Selection>Select by Location>Selection Method>Select Features From >  (NEXT - identify your target layer)Target Layer(s) - check box next to Pensacola Parcels> Source Layer - Parcels Nearby >(Set up a Spatial Relationship)Spatial Selection Method>Target layer(s) features are within the source layer feature (NOTE:  avoid the 3D option since we are using 2D for this lab)>OK

Instructions on how to find features Within other features (find out if a feature is inside or outside a boundary)

8)      Query for finding Land Parcels within a given Zone:  With ArcMap still open>Selection>Select by Location>Selection Method>Select from the currently selected features>(Select your Source Layer = layer that you use to relate to your target layer )Source Layer - choose Enterprise Zones polygon layer >(Set up a Spatial Relationship)Spatial Selection Method: Target layer(s) features are within source layer feature (NOTE:  Also - avoid the 3D option since we are using 2D for this lab)>OK

Instructions on how to find more/different kind features within other features

9)      Query for finding certain features (here = superact_risk_features) within the given land parcels:  (NOTE:  Here then the target (=superact_risk_features)and source (Land Parcels = Pensacola Parcels) layers change. >Expand the "Water" group in the TOC and check to display the "superact_risk_features layer">(To now find features within any of the currently selected parcels) with ArcMap still open>Selection>Select by Location>Selection Method>Select Features From >  (NEXT - identify your target layer)Target Layer(s) - check box under water layer and next to Superact_risk_features (check that this is the only layer selected and all other boxes are unchecked including Pensacola Parcels)>(Select your Source Layer =  layer that you use to relate to your target layer )Source Layer - choose Pensacola Parcels (ONLY Now at this point since finding features within, CHECK the "use selected features" box, because you have changed the feature (ie. the feature is no longer land parcels and in particular Pensacola Parcels))>(Set up a Spatial Relationship)Spatial Selection Method>Target layer(s) features are within the source layer feature (NOTE avoid the 3D option since we are using 2D for this lab)>OK

Finding Intersecting Features (includes the category intersection): 

10)  Query for finding intersecting features (specifically land Parcels that intersect with the Flood Zones):  With ArcMap still open>Selection>Select by Location>Selection Method>Select Features From (Because you have this "select feature from" it is prompting you to identify a target layer so have to do so in the next step)>  (NEXT - identify your target layer)Target Layer(s) - check box next to Pensacola Parcels>(Select your Source Layer = layer that you use to relate to your target layer )Source Layer - choose Flood_Zones layer >(Set up a Spatial Relationship)Spatial Selection Method: Target layer(s) features intersect the source layer feature (NOTE:  Also - avoid the 3D option since we are using 2D for this lab)>OK

Combining Location Queries with Attribute Queries

11)  Query combining an attribute Query and a Location Query: Finding Residential Parcels within Flood Zones - With ArcMap still open>Selection>Select by Attribute>Layers - Pensacola Parcels (leave unchecked box that reads "only show selectable layers in this list">Method - Select from current selection>Double click "USE_GEN" to drop into box> select "=">Click Get Unique Values>double-click "Residential">Verify>OK (when get message "This expression was successfully verified"the currently selected features  (SO SAYING TARGET LAYER = Pensacola Parcels selected in initial Query and just carried on) unlike when you started where you were "selecting features from" so were prompted to identify a target layer)>(Select your Source Layer = layer that you use to relate to your target layer )Source Layer - choose Flood_Zones layer(LEAVE UNCHECKED the "use selected features" box)>(Set up a Spatial Relationship)Spatial Selection Method: Target layer(s) features are within source layer feature (NOTE:  Not typing in "Target layer(s) features are within A DISTANCE of the source layer feature - So the "apply a distance box should be UNCHECKED and skip any previous additional steps after this"  Also - avoid the 3D option since we are using 2D for this lab)>OK (As part of your attribute query, specify a field name, operator and attribute value)>"USE_GEN" = "Residential">OK (Query returns a # = 666)To find the total value of these residential parcels) - Open attribute table>right-click on the Market_Value Column Heading>Statistics

Spatially Joining Tables:

12)  Spatial Join to find the number of risk features within each drainage basin in the county:  Display all 3 layers under the water group in the TOC>rightclick on Superact_risk_features>Joins and Relates - Select Join>Under "what do you want to join to this layer" Select - Join data from another layer based on spatial location>Select the Layer "Drainage Basins" (This will be joined to the superact_risk_feature layer)>It falls inside (this radio button defines the spatial relationship)>Click on browse folder button>Select file path and specify the name of your output file>Save new layer as a feature class in Lab7Data.gdb database in Vector1 folder on s:drive (be sure to check that the "save as file type" is set to File or Personal Geodatabase Feature Class to be able to see the Lab7Data.gdb)>OK (A new layer is added to the map with points located in the same places as the input file superact_risk_feature layer.  Attribute table now also includes information about the drainage basin that each point is located within - find these "Joined fields" at the far right of the attributes table after all the original fields.)>(To find number of risk features within each basin) - scroll through the attributes table columns and find the "Basin" column heading>right click and select "summarize" (DO NOT check any of the boxes in the summarize window)>specify file name and location for output table (same location as previous file.  Name is Basins_sum_Output)>OK>Yes (to add table to map)>(To Find output table in TOC and open it to see results) - right click (do not double click) on table and select Open (Table will have 3 columns, :"Object ID", "Basin" and "Count_Basin" (this field lists the # of superact_risk_features found inside each basin)>(To sort descending) Right click on the Count_Basin column heading and select sort descending>(To find the name and location of the nearest risk features to each of the water wells in the superact_wells_esc layer) Right click on superact_wells_esc layer (check it projection information for its linear units)>Join>Under "what do you want to join to this layer" Select - Join data from another layer based on spatial location>Superact_risk_feature(the join to layer)>select lower radio button to join all attributes from each point layer and a new distance field >specify file name and location for output table>OK(New Layer is added to Map)>(To find the nearest risk feature to the local drinking water water well run by the People's Water Service Company) - open attribute table for join layer just created>Locate "NAME" column>Right click on the NAME column heading and Locate "People's Water Service Company">scroll to the right in the table to find column "Superact_risk_features_NAME" and "Distance" - record this information in Process summary.