LAB EXERCISE #4 Modeling Connectivity
|
|
- Georgia Dennis
- 5 years ago
- Views:
Transcription
1 LAB EXERCISE #4 Modeling Connectivity Instructor: K. McGarigal Overview: In this exercise, you will learn to appreciate the challenges of modeling connectivity and gain practical hands-on experience doing so with R. Specifically, you will: 1) choose among two pre-established applications, 2) examine various connectivity modeling approaches on neutral landscapes, 3) apply various connectivity modeling approaches to your chosen real-world application, and 3) discuss the challenges in choosing and interpreting the results of these analyses. Primary objectives To learn how to model connectivity based using a few common approaches as implemented in R. To gain an appreciation for the many challenges (and alternatives) in quantifying connectivity for a specific application. Part 1: Background information The case study landscape is the Massachusetts portion of the Lower Connecticut River watershed, a HUC6 (Hydrologic Unit Code) level watershed representing the lower half of the Connecticut River watershed. The landscape encompasses 704,882 ha and, as noted in lab1, supports a range of ecological systems and human land uses. Step 1. Establish the objective of the analysis Given a user-specified overall goal pertaining to a pattern-process of interest, the first step of any landscape pattern analysis is to establish a clear objective; specifically, a Specific, Measurable, Attainable, Relevant, and Time-based (SMART) statement of the pattern and process to quantify. Unfortunately, because of the data preparation required for landscape pattern analysis, it is not practical in this brief exercise for you to define your own unique goal and objective and associated landscape definition (step 2). Consequently, while there are numerous possibilities (as you hopefully discovered in lab 1), the goal and objective of the analysis for this exercise were pre-established, as follows. Option 1: Potential vernal pools Goal: Evaluate connectivity of the landscape, as defined below, for vernal pool amphibians at multiple scales and identify and prioritize vernal pools and the associated upland habitat and connections between pools for land protection Pattern: distribution of potential vernal pools and the spatial structure of the intervening landscape, as defined below, within the defined landscape extent. Lab 4 Page 1
2 Process: amphibian seasonal migration (to and from breeding pools), dispersal to non-natal pools, and gene flow across multiple generations. Analysis objective: quantify connectivity of vernal pools to their associated uplands and neighboring vernal pools at the local (i.e., dispersal) and regional (i.e., gene flow) scales within the defined landscape extent in order to prioritize locations for land protection. Option 2: Blackburnian warbler habitat Goal: Evaluate connectivity of the landscape, as defined below, for blackburnian warblers and associated species at multiple scales and identify and prioritize corridors for northward movements to support range shifts under a warming climate scenarios, and locations along the major transportation infrastructure for potential wildlife overpass structures. Pattern: distribution of suitable habitat to support blackburnian warbler movement among potential home ranges and during north-south migration within the defined landscape extent. Process: blackburnian warbler (and associated species) within and among home range movement, seasonal migration, and range shift under a warming climate scenario. Analysis objective: quantify connectivity of quality blackburnian warbler (and associated species) habitat within the defined landscape extent in order to identify priorities for land protection and potential road crossing structures. Step 2. Define the landscape The next step is to define the landscape in accordance with the objective. Option 1: Potential vernal pools Conceptual model.-we will use the patch mosaic model to represent landscape conductance/resistance for vernal pool amphibians. Thematic content.- -the landscape was classified into land cover classes based on the dominant ecological setting and human land use (as represented in the DSLland cover map). Thematic resolution.- the original DSLland cover map containing 65 classes was reclassified into "resistance" values based on previous work (although this can be modified if you are so inclined). Spatial grain.- the landscape was represented as a raster at 30 m resolution - the original dslland resolution. Thus, each 30x30 m cell was classified into one of the 8 land cover classes. Note, for this exercise we maintained the 1 cell minimum mapping unit (i.e., minimum patch size) and did not kernel smooth the raster. Lab 4 Page 2
3 However, for some of the connectivity analyses it was necessary to aggregate the raster to a coarser resolution (e.g., 90 m) based on a focal mean of resistance. Spatial extent.-the landscape extent was defined as an arbitrary rectangular area centered on the Pioneer Valley (see GIS layers below), with the extent chosen to be small enough to allow for processing efficiency. Fragmenting features.- all roads and rivers were treated as major impediments (but not complete barriers) to movement in the reclassification of the DSLland cover into resistance values. Landscape boundary and context.-the boundary was arbitrary, as noted above, and thus offers no meaningful ecological unit for this analysis, other than to provide an area small enough to render the processing practical for this exercise. Option 2: Blackburnian warbler habitat Conceptual model.-we will use the landscape gradient model to represent landscape conductance/resistance for blackburnian warbler habitat. Thematic content.- -the landscape was modeled in a previous lab to produce a resource selection function that generated a continuous surface representing the predicted probability of occurrence. Here, we will use the RSF value as an estimate of conductance for movement, even though in this case the RSF was not based on movement data per se (more on this later). Thematic resolution.- the precision of the RSF probability of occurrence. Spatial grain.- the landscape was represented as a raster at 30 m resolution - the original dslland resolution. Thus, each 30x30 m cell was given a predicted probability of occurrence (0-1), which we interpreted as conductance. However, for some of the connectivity analyses it was necessary to aggregate the raster to a coarser resolution (e.g., 90 m) based on a focal mean of conductance. Spatial extent.-the landscape extent was defined as an arbitrary rectangular area centered on the Pioneer Valley (see GIS layers below), with the extent chosen to be small enough to allow for processing efficiency. Fragmenting features.- no fragmenting features per se were identified other than as already incorporated in the RSF function. Landscape boundary and context.-the boundary was arbitrary, as noted above, and thus offers no meaningful ecological unit for this analysis, other than to provide an area small enough to render the processing practical for this exercise. Step 3. Inspect the GIS data Now that we have defined the landscape, take some time to familiarize yourself with the GIS data associated with the two applications. Open up in ArcMap the following project file: Lab 4 Page 3
4 ...\exercises\lab3\lab4.mxd Take some time to review each of the data layers. The ArcMap table of contents is organized into four groups of layers, from bottom to top, as follows: Lower Connecticut watershed - suite of layers encompassing the Lower Connecticut River HUC6 watershed. These layers provide a spatial context for the focal landscape: o hillshade.tif - backdrop for other layers if desired. o dslland_2010.tif - original 30m resolution land cover map from the DSL project containing 65 classes within the project area. o cover7_2010.tif - dslland reclassified into 7 major cover types based on ecological formation and human land use. Blackburnian warblers - couple of layers with blackburnian warbler point observations and resource selection function: o blbwpred.tif - 30 m raster representing the resource selection function generated from a previous lab, in which the cell values represent the probability of blackburnian warbler occurrence. o blbw - point shapefile containing the ebird sample locations, with the symbology present to display present versus absent locations. Potential vernal pools - point shapefile depicting potential vernal pools, given for the Massachusetts portion of the lower Connecticut River watershed (PVP). Vector overlays - a handful of shapefiles that may be useful as overlays, including: o boundary - boundary of the arbitrary project area centered on the Pioneer Valley in the lower Connecticut River watershed within Massachusetts. o NHDhighRes - streams layer derived from the National Hydrography data based on high resolution 1:24 k streams. o DSLroads - roads layer with the symbology set to display road class. Part 2: Data prep In this section, we will prepare the spatial data we will need for the functions below. Note, the detailed process, and thus the R scripting, differs slightly between the two applications owing to differences in the source data and analysis objective. The detailed steps below will be presented for the potential vernal pool application and significant modifications for the blackburnian warbler application will be noted, but separate R scripts are provided for the two applications to simplify your execution. In both cases, you will need to source the landeco.r package and install (if need be) and load the other required libraries, as given in the provided lab3 R scripts. Lab 4 Page 4
5 Step 2a. read shapefiles First, read in the boundary (polygon) shapefile: setwd('c:\\work\\landeco-umass\\exercises\\lab4\\gisdata\\') boundary<-readogr('.','boundary') Next, read in the PVP (points) shapefile and clip to the boundary extent: pvppoints<-readogr('.','pvp') pvppoints<-pvppoints[boundary,] Next, given the large number of pools, for computation efficiency it is necessary to randomly sample a subset of the pools, otherwise you will be processing for weeks: pvppoints<-pvppoints[sample(1:nrow(pvppoints),size=10,replace=false),] Finally, plot the boundary with the randomly sampled pools: plot(boundary) points(pvppoints) Step 2b. create resistance surface First, read in the dslland cover layer and crop to the boundary extent, and optionally write the resulting surface to disk for viewing in ArcMap: temp<-raster('dslland_2010.tif') temp<-crop(temp,boundary) #temp<-mask(temp,boundary) #only need if boundary irregular writeraster(temp,'dslland_2010_boundary.tif',datatype='int2u',overwrite=true) Next, reclass land cover into resistance values based on the pre-established cross-walk provided in the dsllandkey table (feel free to change the cross-walk): reclass<-read.csv('dsllandkey.csv',header=true) rclmat<-as.matrix(subset(reclass,select=c('export_value','reclasspvp'))) temp<-reclassify(temp,rclmat,include.lowest=true) Next, coarsen resistance raster for processing efficiency with the gdistance functions below, and optionally write the resulting raster to disk for viewing in ArcMap. Note, here the fact=3 argument results in a 3x3 focal window (or 90 m) and the fun=mean says to replace the coarse cell with the mean of the 9 subsumed cells: pvpresist<-aggregate(temp,fact=3,fun=mean,na.rm=true) writeraster(pvpresist,'pvpresist.tif',datatype='flt4s',overwrite=true) Finally, plot the resistance surface with the randomly sampled pools: plot(pvpresist) points(pvppoints,pch=19) Lab 4 Page 5
6 Step 2c. select random sample points [blackburnian warbler application] Note, for the potential vernal pool application we randomly sampled a subset of the preexisting PVP points (see above). However, for the blackburnian warbler application we may not want to use the ebird point locations, as they are somewhat arbitrary and their density and distribution don't really reflect meaningful locations for blackburnian warblers; they merely reflect where the ebirders decided to conduct a survey. It may be preferable in this case to randomly sample the RSF surface so that the sample point locations at least reflect the expected density and distribution of blackburnian warblers. To select points randomly based on inclusion probabilities derived from the RSF surface, we can use the probpointselect() function in the landeco.r library. Specifically, designate the surface containing relative inclusion probabilities (although the surface does not have to already have a probability scale) and the number of points to select: blbwpoints<-probpointselect(r=blbwrsf,n=100) Also, because we are interested in potential corridors for south to north movement, either during seasonal migration or range shift under a warming climate scenario, we may also want to create a set of "from" points on the southern boundary of the landscape and a set of "to" points on the northern boundary of the landscape and ask all the southern points to go visit all of the northern points. The following script will create a set of 20 evenly spaced points along each of the southern and northern boundary of the landscape: xcoord<-seq(xmin(blbwrsf),xmax(blbwrsf),length=20) ycoord<-rep(ymin(blbwrsf),length=20) blbwsouth<-spatialpoints(cbind(xcoord,ycoord)) ycoord<-rep(ymax(blbwrsf),length=20) blbwnorth<-spatialpoints(cbind(xcoord,ycoord)) Finally, plot the resistance surface with the randomly sampled points and the southern and northern boundary points plot(blbwrsf) points(blbwpoints,pch=19) points(blbwsouth,pch=19,col='red') points(blbwnorth,pch=19,col='blue') Part 3: gdistance functions In this section, we will derive least cost paths and random low cost paths using functions in the gdistance library. To begin, we need to create a conductance transition matrix from the resistance or conductance raster. Note, conductance and is the reciprocal of resistance. In the case of the potential vernal pools, we created a resistance matrix, so here we have to convert resistance to conductance by taking the reciprocal of resistance in the transition function that follows: Lab 4 Page 6
7 pvpconduct<-transition(pvpresist,function(x) 1/mean(x,na.rm=TRUE), direction=8) Note, because we used an 8-neighbor rule in the transition matrix, we need to make a geographic distance correction for the diagonals. In addition, it is usually helpful to scale the conductance values using the scl=true argument to allow for more efficient computations: pvpconduct<-geocorrection(pvpconduct,scl=true) Step 3a. least cost surface First, create a least cost surface; specifically, a surface that depicts the accumulated least cost from every cell to the least costly point location. Hence, the cell value represents the accumulated cost along the least cost path to the "closest" point location: pvpcostsurf<-acccost(pvpconduct,pvppoints) Next, plot the least cost surface with the randomly sampled point locations, and optionally write the resulting raster to disk for viewing in ArcMap: plot(pvpcostsurf) points(pvppoints,pch=19) writeraster(pvpcostsurf,'pvpcostsurf.tif', datatype='flt4s',overwrite=true) Step 3b. least cost paths Next, compute least cost paths among the point locations. First, compute the least cost patch between two arbitrarily selected points (1 and 10 in the example below) and plot the result along with the points: lcp<-shortestpath(pvpconduct,pvppoints[1,],pvppoints[10,], output='spatiallines') plot(pvpresist) points(pvppoints,pch=19) points(pvppoints[1,],pch=19,col='red') points(pvppoints[10,],pch=19,col='blue') lines(lcp) Next, compute factorial least cost paths by computing and overlaying the least cost path between every pairwise combination of point locations, and optionally write the resulting raster to disk for viewing in ArcMap (CAUTION, this takes a long time to run, so you might want to let it run while you go for a run): pvpfactlcp<-factoriallcp(r=pvpresist,points=pvppoints,conduct=pvpconduct) writeraster(pvpfactlcp,'pvpfactlcp.tif',datatype='int2u',overwrite=true) Step 3c. Random low cost paths Next, compute random low cost paths among the point locations. Note, in the passage() function used below, the theta= argument controls the magnitude of random deviations Lab 4 Page 7
8 from the straight-line patch between two points, with a theoretical range of close to zero (but not zero) to a maximum of 20. Values close to zero, the default, result in a completely random walk (i.e., a "drunkards" walk) and is analogous to flow in an electrical conduction surface (as implemented in the software Circuitscape). Larger values result in walks that are increasingly focused on the straight-line path between two points. However, I have found that with values above 1 to 5, depending on the landscape, that the resulting passage values are all zero and thus provide no useful output. Consequently, it warrants some experimenting beforehand, but likely values between will work best. First, compute random low cost paths between two arbitrarily selected points (1 and 10 in the example below) and plot the result along with the points. pvppassage<-passage(pvpconduct,pvppoints[1,],pvppoints[10,],theta=0.05) plot(pvppassage) points(pvppoints[1,],pch=19,col='red') points(pvppoints[10,],pch=19,col='blue') Next, compute factorial random low cost paths; i.e., between all pairwise combinations of points and summing the result (CAUTION: this takes bloody forever to run, so it is probably best to run this overnight or over the weekend), and optionally write the resulting raster to disk for viewing in ArcMap: pvpfactrlcp<-factorialrlcp(r=pvpresist,points=pvppoints, conduct=pvpconduct,theta=0.05) writeraster(pvpfactrlcp,'pvpfactrlcp.tif',datatype='flt4s',overwrite=true) Part 4: Resistant kernels In this section, we will compute resistant kernels and derive several different metrics. To begin, we must switch to 32-bit R in order to use the gridio library functions. This can be done in RStudio by selecting Tools Global options General R version and change to 32-bit, and then closing and restarting RStudio. Once you are back in RStudio you will need to resource the landeco.r library and load a few other libraries as indicated in the R script provided. Step 4a. Read shapefiles First, read in the boundary shapefile: setwd('c:\\work\\landeco-umass\\exercises\\lab4\\gisdata\\') boundary<-readogr('.','boundary') Next, read in the points shapefile and clip to the project boundary. Note, here we are NOT going to randomly sample a subset of the points as we did before because the resistant kernel function is efficient enough to work rapidly with thousands of points: pvppoints<-readogr('.','pvp') pvppoints<-pvppoints[boundary,] Lab 4 Page 8
9 Finally, plot the boundary with the points to verify: plot(boundary) points(pvppoints,cex=0.5) Step 4b. Create resistance surface First, read in the dlsland raster and reclass it to represent resistance values based on the cross-walk in the dsllandkey table. Note, here we are NOT going to coarsen the grid as we did before because the resistant kernel function is efficient enough to work with the 30 m resolution data: temp<-raster('dslland_2010_boundary.tif') reclass<-read.csv('dsllandkey.csv',header=true) rclmat<-as.matrix(subset(reclass,select=c('export_value','reclasspvp'))) pvpresist<-reclassify(temp,rclmat,include.lowest=true) Finally, plot the resistance surface with the points to verify? plot(pvpresist) points(pvppoints,cex=0.5) Step 4c. Compute resistant kernels First, in order to use the resistant kernel function spread() in the gridio library, we need to initialize the gridio library with any of the rasters having the desired template (i.e., extent and cell size) as follows: gridinit() pvpresist<-as.grid(pvpresist) setwindow(pvpresist) Next, compute a suite of resistant kernel metrics using the rkernel() function in the landeco.r package, which makes use of the spread() function in gridio: result<-rkernel(r=pvpresist,p=pvppoints,sd=1000) Note, in the function above, the resistant kernel is Gaussian and the sd=1000 specifies a standard deviation of 1,000 m. The resistant kernel extends outward from the focal cell to a maximum of 3 times the standard deviation by default, although this can be changed with the sd.threshold= argument. The rkernel() function returns a list object with two components. The first component is named "cumrk" and contains the cumulative resistant kernel surface in grid form. The following lines extract the cumrk grid from the result object, converts the grid to a raster class object, and then reassigns the original coordinate reference system back to the resulting raster object: pvpcumrk1000<-raster(result$cumrk) crs(pvpcumrk1000)<-crs(pvpresist) Lab 4 Page 9
10 Next, plot the cumulative resistant kernel surface, and optionally write the resulting raster to disk for viewing in ArcMap: plot(pvpcumrk1000) writeraster(pvpcumrk1000,'pvpcumrk1000.tif',datatype='flt4s', overwrite=true) Finally, extract the SpatialPointsDataFrame object from the second list component named "points", and optionally write the resulting object as a shapefile to disk for viewing in ArcMap: writeogr(result$points,dsn='c:\\work\\landeco-umass\\exercises\\lab4\\gisdata', layer='pvprkern',driver='esri Shapefile') The SpatialPointsDataFrame object contains the original point attributes plus three additional attributes: rkern = the raw sum of the resistant kernel derived for the focal point traverse = traversability index, defined as rkern divided by rkern for a corresponding standard kernel (i.e., resistant kernel applied to a uniform nonresistant surface); hence, traverse ranges from 0 for a focal point completely disconnected from neighboring cells to a theoretical maximum of 1 for a focal point surrounded by non-resistant surface out to a distance defined by three times (by default) the specified standard deviation (sd) of the Gaussian kernel, which in the example above equates to 3,000 m. isolate = isolation index (or, inversely, connectedness index), defined as the sum of neighboring resistant kernels minus the resistant kernel derived from the kernel built for the focal point; hence, isolate measures the total flow from neighboring points to the focal point. Isolate ranges from 0 for a completely isolated cell to increasing larger values for focal points that have increasing numbers of connected nearby points. The theoretical maximum is achieved when every neighboring cell is a point and the neighborhood is non-resistant. Assignment: Option 1: Potential vernal pools Q1. Based on the multi-metric connectivity assessment, identify and prioritize vernal pools and associated uplands and intervening matrix lands for land protection. Hint, use all forms of connectivity analysis, including least cost paths, random low cost paths and the various resistant kernel products to help answer this question. Be sure to graphically show your results with respect to priority locations. In deriving your solution, be sure to discuss and interpret what each connectivity measure means to demonstrate your understanding. Q2. Discuss the major challenges and limitations in conducting a meaningful connectivity analysis to meet your goal. What could you do differently or additionally to further address your goal? Lab 4 Page 10
11 Option 2: Blackburnian warbler habitat Q1. Based on the multi-metric connectivity assessment, identify and prioritize lands with high connectivity value for blackburnian warblers and associated species. Specifically, identify and prioritize potential corridors for northward movements to support range shifts under a warming climate scenarios, and locations along the major transportation infrastructure for potential wildlife overpass structures. Hint, use all forms of connectivity analysis, including least cost paths, random low cost paths and the various resistant kernel products to help answer this question. Be sure to graphically show your results with respect to priority locations. In deriving your solution, be sure to discuss and interpret what each connectivity measure means to demonstrate your understanding. Q2. Discuss the major challenges and limitations in conducting a meaningful connectivity analysis to meet your goal. What could you do differently or additionally to further address your goal? Lab 4 Page 11
LAB EXERCISE #2 Quantifying Patch Mosaics
LAB EXERCISE #2 Quantifying Patch Mosaics Instructors: K. McGarigal Overview: In this exercise, you will learn to appreciate the challenges of quantifying patch mosaics and gain practical hands-on experience
More informationLAB EXERCISE #3 Quantifying Point and Gradient Patterns
LAB EXERCISE #3 Quantifying Point and Gradient Patterns Instructor: K. McGarigal Overview: In this exercise, you will learn to appreciate the challenges of quantifying point and gradient patterns and gain
More informationGeographical Information Systems
Geographical Information Systems Geographical Information Systems (GIS) is a relatively new technology that is now prominent in the ecological sciences. This tool allows users to map geographic features
More informationTutorial 8 Raster Data Analysis
Objectives Tutorial 8 Raster Data Analysis This tutorial is designed to introduce you to a basic set of raster-based analyses including: 1. Displaying Digital Elevation Model (DEM) 2. Slope calculations
More informationGIS IN ECOLOGY: ANALYZING RASTER DATA
GIS IN ECOLOGY: ANALYZING RASTER DATA Contents Introduction... 2 Raster Tools and Functionality... 2 Data Sources... 3 Tasks... 4 Getting Started... 4 Creating Raster Data... 5 Statistics... 8 Surface
More informationWEBINAR AGENDA. Review of connectivity mapping progress to date. Describe draft methodology for delineating and validating core areas
WEBINAR AGENDA Review of connectivity mapping progress to date Describe draft methodology for delineating and validating core areas Review examples of the methodology implemented in the Highlands Explain
More informationDelineation of Watersheds
Delineation of Watersheds Adirondack Park, New York by Introduction Problem Watershed boundaries are increasingly being used in land and water management, separating the direction of water flow such that
More informationExercise 6: Using Burn Severity Data to Model Erosion Risk
Exercise 6: Using Burn Severity Data to Model Erosion Risk Document Updated: November 2009 Software Versions: ERDAS Imagine 9.3 and ArcGIS 9.3, Microsoft Office 2007 Introduction A common use of burn severity
More informationUsing the Stock Hydrology Tools in ArcGIS
Using the Stock Hydrology Tools in ArcGIS This lab exercise contains a homework assignment, detailed at the bottom, which is due Wednesday, October 6th. Several hydrology tools are part of the basic ArcGIS
More informationLab 7: Cell, Neighborhood, and Zonal Statistics
Lab 7: Cell, Neighborhood, and Zonal Statistics Exercise 1: Use the Cell Statistics function to detect change In this exercise, you will use the Spatial Analyst Cell Statistics function to compare the
More informationDisplaying and Rotating WindNinja-Derived Wind Vectors in ArcMap 10.5
Displaying and Rotating WindNinja-Derived Wind Vectors in ArcMap 10.5 Chuck McHugh RMRS, Fire Sciences Lab, Missoula, MT, 406-829-6953, cmchugh@fs.fed.us 08/01/2018 Displaying WindNinja-generated gridded
More informationCreating Watersheds from a DEM
Creating Watersheds from a DEM These instructions enable you to create watersheds of specified area using a good quality Digital Elevation Model (DEM) in ArcGIS 8.1. The modeling is performed in ArcMap
More informationWatershed Delineation
Watershed Delineation Jessica L. Watkins, University of Georgia 2 April 2009 Updated by KC Love February 25, 2011 PURPOSE For this project, I delineated watersheds for the Coweeta synoptic sampling area
More informationNatalie Cabrera GSP 370 Assignment 5.5 March 1, 2018
Network Analysis: Modeling Overland Paths Using a Least-cost Path Model to Track Migrations of the Wolpertinger of Bavarian Folklore in Redwood National Park, Northern California Natalie Cabrera GSP 370
More informationGIS IN ECOLOGY: ANALYZING RASTER DATA
GIS IN ECOLOGY: ANALYZING RASTER DATA Contents Introduction... 2 Tools and Functionality for Raster Data... 2 Data Sources... 3 Tasks... 4 Getting Started... 4 Creating Raster Data... 5 Summary Statistics...
More informationChapter 6. Fundamentals of GIS-Based Data Analysis for Decision Support. Table 6.1. Spatial Data Transformations by Geospatial Data Types
Chapter 6 Fundamentals of GIS-Based Data Analysis for Decision Support FROM: Points Lines Polygons Fields Table 6.1. Spatial Data Transformations by Geospatial Data Types TO: Points Lines Polygons Fields
More informationIn this exercise we will learn how to use the analysis tools in ArcGIS with vector and raster data to further examine potential building sites.
GIS Level 2 In the Introduction to GIS workshop we filtered data and visually examined it to determine where to potentially build a new mixed use facility. In order to get a low interest loan, the building
More informationVegetation and Wildlife Habitat Mapping Study in the Upper and Middle Susitna Basin Study Plan Section 11.5
(FERC No. 14241) Vegetation and Wildlife Habitat Mapping Study in the Upper and Middle Susitna Basin Study Plan Section 11.5 Initial Study Report Part C: Executive Summary and Section 7 Prepared for Prepared
More informationGIS CONCEPTS ARCGIS METHODS AND. 2 nd Edition, July David M. Theobald, Ph.D. Natural Resource Ecology Laboratory Colorado State University
GIS CONCEPTS AND ARCGIS METHODS 2 nd Edition, July 2005 David M. Theobald, Ph.D. Natural Resource Ecology Laboratory Colorado State University Copyright Copyright 2005 by David M. Theobald. All rights
More informationLab 1: Landuse and Hydrology, learning ArcGIS II. MANIPULATING DATA
Lab 1: Landuse and Hydrology, learning ArcGIS II. MANIPULATING DATA As you experienced in the first lab session when you created a hillshade, high resolution data can be unwieldy if you are trying to perform
More informationSpatial Units, Scaling and Aggregation (Level 1) October 2017
Spatial Units, Scaling and Aggregation (Level 1) October 2017 Overview: Spatial Units 1. Learning objectives 2. Review of Level 0 (5m) 3. Level 1 (Compilers): Presentation & group exercise Spatial units
More informationNative species (Forbes and Graminoids) Less than 5% woody plant species. Inclusions of vernal pools. High plant diversity
WILLAMETTE VALLEY WET-PRAIRIE RESTORATION MODEL WHAT IS A WILLAMETTE VALLEY WET-PRAIRIE Hot Spot s Native species (Forbes and Graminoids) Rare plant species Less than 5% woody plant species Often dominated
More informationLeon Creek Watershed October 17-18, 1998 Rainfall Analysis Examination of USGS Gauge Helotes Creek at Helotes, Texas
Leon Creek Watershed October 17-18, 1998 Rainfall Analysis Examination of USGS Gauge 8181400 Helotes Creek at Helotes, Texas Terrance Jackson MSCE Candidate University of Texas San Antonio Abstract The
More informationWhat are the five components of a GIS? A typically GIS consists of five elements: - Hardware, Software, Data, People and Procedures (Work Flows)
LECTURE 1 - INTRODUCTION TO GIS Section I - GIS versus GPS What is a geographic information system (GIS)? GIS can be defined as a computerized application that combines an interactive map with a database
More informationReminder that we update the website: with new information, project updates, etc.
1 Reminder that we update the website: http://www.njfishandwildlife.com/ensp/connectivity_proj.htm with new information, project updates, etc. regularly. Any suggestions are welcome on how to make the
More information+ = Spatial Analysis of Raster Data. 2 =Fault in shale 3 = Fault in limestone 4 = no Fault, shale 5 = no Fault, limestone. 2 = fault 4 = no fault
Spatial Analysis of Raster Data 0 0 1 1 0 0 1 1 1 0 1 1 1 1 1 1 2 4 4 4 2 4 4 2 4 4 4 2 4 4 2 4 4 3 4 4 4 2 3 + = 0 = shale 1 = limestone 2 = fault 4 = no fault 2 =Fault in shale 3 = Fault in limestone
More informationLand Cover Data Processing Land cover data source Description and documentation Download Use Use
Land Cover Data Processing This document provides a step by step procedure on how to build the land cover data required by EnSim. The steps provided here my be long and there may be short cuts (like using
More informationHow to Create Stream Networks using DEM and TauDEM
How to Create Stream Networks using DEM and TauDEM Take note: These procedures do not describe all steps. Knowledge of ArcGIS, DEMs, and TauDEM is required. TauDEM software ( http://hydrology.neng.usu.edu/taudem/
More informationExercise 2: Working with Vector Data in ArcGIS 9.3
Exercise 2: Working with Vector Data in ArcGIS 9.3 There are several tools in ArcGIS 9.3 used for GIS operations on vector data. In this exercise we will use: Analysis Tools in ArcToolbox Overlay Analysis
More informationCreating Watersheds from a DEM in ArcGIS 9.x
Creating Watersheds from a DEM in ArcGIS 9.x These instructions enable you to create watersheds (a.k.a. catchments or basins) using a good quality Digital Elevation Model (DEM) in ArcGIS 9.1. The modeling
More informationSpatial Analysis of Raster Data
Spatial Analysis of Raster Data 0 0 1 1 0 0 1 1 1 0 1 1 1 1 1 1 2 4 4 4 2 4 5 5 4 2 4 4 4 2 5 5 4 4 2 4 5 4 3 5 4 4 4 2 5 5 5 3 + = 0 = shale 1 = limestone 2 = fault 4 = no fault 2 =Fault in shale 3 =
More informationProtocol for Prioritizing Conservation Opportunity Areas in Centre County and Clinton County
Protocol for Prioritizing Conservation Opportunity Areas in Centre County and Clinton County Chesapeake Conservancy has developed this methodology to prioritize conservation opportunity areas in Centre
More informationGIS Exercise: Analyses of Home Range and Space Use Patterns
WLF 315 Wildlife Ecology Lab I Fall 2012 GIS Exercise: Analyses of Home Range and Space Use Patterns In this exercise, you will use the Geospatial Modeling Environment (http://www.spatialecology.com/gme)
More informationFR Exam 2 Substitute Project!!!!!!! 1
FR5131 - Exam 2 Substitute Project!!!!!!! 1 Goal We seek to identify county lands that both protect vernal pools and riparian corridors, and that provide open space and recreation within easy reach of
More informationRoad-based Landscape Metrics for Quantifying Habitat Fragmentation
University of South Florida Scholar Commons Graduate Theses and Dissertations Graduate School 2011 Road-based Landscape Metrics for Quantifying Habitat Fragmentation Rebecca Whitehead Loraamm University
More informationUsing ArcGIS for Hydrology and Watershed Analysis:
Using ArcGIS 10.2.2 for Hydrology and Watershed Analysis: A guide for running hydrologic analysis using elevation and a suite of ArcGIS tools Anna Nakae Feb. 10, 2015 Introduction Hydrology and watershed
More informationUsing Python to Automate Map Exports
Using Python to Automate Map Exports Don Katnik Maine Department of Inland Fisheries and Wildlife Northeast ARC Users' Conference, 2016 Falmouth, MA DISCLAIMER: I don't work for Python Background Maine
More informationExercise 2: Working with Vector Data in ArcGIS 9.3
Exercise 2: Working with Vector Data in ArcGIS 9.3 There are several tools in ArcGIS 9.3 used for GIS operations on vector data. In this exercise we will use: Analysis Tools in ArcToolbox Overlay Analysis
More informationAbstract: Contents. Literature review. 2 Methodology.. 2 Applications, results and discussion.. 2 Conclusions 12. Introduction
Abstract: Landfill is one of the primary methods for municipal solid waste disposal. In order to reduce the environmental damage and to protect the public health and welfare, choosing the site for landfill
More informationModeling Disease Transmission in Long-tailed Macaques on Bali
Modeling Disease Transmission in Long-tailed Macaques on Bali Kelly Lane Gerhard Niederwieser Ryan Kennedy University of Notre Dame Macaque Background Coexisted in temples across the island for at least
More informationGRAPEVINE LAKE MODELING & WATERSHED CHARACTERISTICS
GRAPEVINE LAKE MODELING & WATERSHED CHARACTERISTICS Photo Credit: Lake Grapevine Boat Ramps Nash Mock GIS in Water Resources Fall 2016 Table of Contents Figures and Tables... 2 Introduction... 3 Objectives...
More informationTHE CONSERVATION LANDSCAPE CONTEXT TOOLBOX: A Custom ESRI ArcToolbox for the Stewardship of Conservation Lands in Rhode Island
THE CONSERVATION LANDSCAPE CONTEXT TOOLBOX: A Custom ESRI ArcToolbox for the Stewardship of Conservation Lands in Rhode Island Ann Borowik M.E.S.M. Research Project University of Rhode Island May 2008
More informationDesigning Sustainable Landscapes: Aquatic Barriers settings variable
Designing Sustainable Landscapes: Aquatic Barriers settings variable A project of the University of Massachusetts Landscape Ecology Lab Principals: Kevin McGarigal, Professor Brad Compton, Research Associate
More informationThe two foundation ideas or concepts that are of vital importance in GIS work (in my humble opinion) are:
HOME RANGE ESTIMATION AND HABITAT SELECTION USING G.I.S. During this workshop, you will be introduced to ArcGIS 9 software and some of the analytical tools available to go with this software. This workshop
More informationPriority areas for grizzly bear conservation in western North America: an analysis of habitat and population viability INTRODUCTION METHODS
Priority areas for grizzly bear conservation in western North America: an analysis of habitat and population viability. Carroll, C. 2005. Klamath Center for Conservation Research, Orleans, CA. Revised
More informationENMTools: a toolbox for comparative studies of environmental niche models
Ecography 33: 607611, 2010 doi: 10.1111/j.1600-0587.2009.06142.x # 2010 The Authors. Journal compilation # 2010 Ecography Subject Editor: Thiago Rangel. Accepted 4 September 2009 ENMTools: a toolbox for
More informationSteve Pye LA /22/16 Final Report: Determining regional locations of reference sites based on slope and soil type. Client: Sonoma Land Trust
Steve Pye LA 221 04/22/16 Final Report: Determining regional locations of reference sites based on slope and soil type. Client: Sonoma Land Trust Deliverables: Results and working model that determine
More informationSpatial Analyst: Multiple Criteria Evaluation Material adapted from FOR 4114 developed by Forestry Associate Professor Steve Prisley
Spatial Analyst: Multiple Criteria Evaluation Material adapted from FOR 4114 developed by Forestry Associate Professor Steve Prisley Section 1: Data In this exercise we will be working with several types
More informationGeography 281 Map Making with GIS Project Four: Comparing Classification Methods
Geography 281 Map Making with GIS Project Four: Comparing Classification Methods Thematic maps commonly deal with either of two kinds of data: Qualitative Data showing differences in kind or type (e.g.,
More information(THIS IS AN OPTIONAL BUT WORTHWHILE EXERCISE)
PART 2: Analysis in ArcGIS (THIS IS AN OPTIONAL BUT WORTHWHILE EXERCISE) Step 1: Start ArcCatalog and open a geodatabase If you have a shortcut icon for ArcCatalog on your desktop, double-click it to start
More informationLab#8: Working With Geodatabases. create a geodatabase with feature datasets, tables, raster datasets, and raster catalogs
Lab#8: Geodatabase Concepts Page#1 of 25 In this lab, you will learn how to: Lab#8: Working With Geodatabases create a geodatabase with feature datasets, tables, raster datasets, and raster catalogs join
More informationTHE DIGITAL TERRAIN MAP LIBRARY: AN EXPLORATIONIST S RESOURCE
THE DIGITAL TERRAIN MAP LIBRARY: AN EXPLORATIONIST S RESOURCE By I.C.L. Webster, P.J. Desjardins and W.E. Kilby KEYWORDS: digital terrain maps, digital terrain stability maps, surficial geology, GIS, raster
More informationAnalysis of Change in Land Use around Future Core Transit Corridors: Austin, TX, Eric Porter May 3, 2012
Analysis of Change in Land Use around Future Core Transit Corridors: Austin, TX, 1990-2006 PROBLEM DEFINITION Eric Porter May 3, 2012 This study examines the change in land use from 1990 to 2006 in the
More informationLAB EXERCISE #3 Neutral Landscape Analysis Summary of Key Results and Conclusions
LAB EXERCISE #3 Neutral Landscape Analysis Summary of Key Results and Conclusions Below is a brief summary of the key results and conclusions of this exercise. Note, this is not an exhaustive summary,
More informationCompilation of GIS data for the Lower Brazos River basin
Compilation of GIS data for the Lower Brazos River basin Francisco Olivera, Ph.D., P.E. Srikanth Koka Lauren Walker Aishwarya Vijaykumar Department of Civil Engineering December 5, 2011 Contents Brief
More informationPython Raster Analysis. Kevin M. Johnston Nawajish Noman
Python Raster Analysis Kevin M. Johnston Nawajish Noman Outline Managing rasters and performing analysis with Map Algebra How to access the analysis capability - Demonstration Complex expressions and optimization
More informationDetermining the Location of the Simav Fault
Lindsey German May 3, 2012 Determining the Location of the Simav Fault 1. Introduction and Problem Formulation: The issue I will be focusing on involves interpreting the location of the Simav fault in
More informationApplication of GIS and remote sensing in conservation of vernal pools
Priyanka Patel NRS 509 12/16/2014 Application of GIS and remote sensing in conservation of vernal pools Vernal pools are small temporary water bodies that contain water for some parts of the year. They
More informationJJ Munoz. Explanation of the Project/Outline
JJ Munoz Helper GEO 386G 1 December, 2016 GIS Project: Using ArcHydro and Surface Roughness tools to Determine Alluvial Fan Catchment Area Morphology Statistics from 1-meter resolution LiDAR DEM s. Explanation
More informationOverview key concepts and terms (based on the textbook Chang 2006 and the practical manual)
Introduction Geo-information Science (GRS-10306) Overview key concepts and terms (based on the textbook 2006 and the practical manual) Introduction Chapter 1 Geographic information system (GIS) Geographically
More informationOutcrop suitability analysis of blueschists within the Dry Lakes region of the Condrey Mountain Window, North-central Klamaths, Northern California
Outcrop suitability analysis of blueschists within the Dry Lakes region of the Condrey Mountain Window, North-central Klamaths, Northern California (1) Introduction: This project proposes to assess the
More informationHydrology and Watershed Analysis
Hydrology and Watershed Analysis Manual By: Elyse Maurer Reference Map Figure 1. This map provides context to the area of Washington State that is being focused on. The red outline indicates the boundary
More information11. Kriging. ACE 492 SA - Spatial Analysis Fall 2003
11. Kriging ACE 492 SA - Spatial Analysis Fall 2003 c 2003 by Luc Anselin, All Rights Reserved 1 Objectives The goal of this lab is to further familiarize yourself with ESRI s Geostatistical Analyst, extending
More informationNeighborhood Locations and Amenities
University of Maryland School of Architecture, Planning and Preservation Fall, 2014 Neighborhood Locations and Amenities Authors: Cole Greene Jacob Johnson Maha Tariq Under the Supervision of: Dr. Chao
More informationGIS APPLICATIONS IN SOIL SURVEY UPDATES
GIS APPLICATIONS IN SOIL SURVEY UPDATES ABSTRACT Recent computer hardware and GIS software developments provide new methods that can be used to update existing digital soil surveys. Multi-perspective visualization
More information4. GIS Implementation of the TxDOT Hydrology Extensions
4. GIS Implementation of the TxDOT Hydrology Extensions A Geographic Information System (GIS) is a computer-assisted system for the capture, storage, retrieval, analysis and display of spatial data. It
More informationUrban Growth and Development Using SLEUTH: Philadelphia Metropolitan Region
Urban Growth and Development Using SLEUTH: Philadelphia Metropolitan Region ALYSSA LYND CODY HITT GUS FREDERICK 1.0 INTRODUCTION The metropolitan and suburban regions of Philadelphia and its surrounding
More informationGIS Topographic Wetness Index (TWI) Exercise Steps
GIS Topographic Wetness Index (TWI) Exercise Steps October 2016 Jeffrey L. Zimmerman, Jr. GIS Analyst James P. Shallenberger Manager, Monitoring & Protection Susquehanna River Basin Commission Table of
More informationThis lab exercise will try to answer these questions using spatial statistics in a geographic information system (GIS) context.
by Introduction Problem Do the patterns of forest fires change over time? Do forest fires occur in clusters, and do the clusters change over time? Is this information useful in fighting forest fires? This
More informationModeling Fish Assemblages in Stream Networks Representation of Stream Network Introduction habitat attributes Criteria for Success
Modeling Fish Assemblages in Stream Networks Joan P. Baker and Denis White Western Ecology Division National Health & Environmental Effects Research Laboratory U.S. Environmental Protection Agency baker.joan@epa.gov
More informationSuitability Analysis on Second Home Areas Selection in Smithers British Columbia
GEOG 613 Term Project Suitability Analysis on Second Home Areas Selection in Smithers British Columbia Zhengzhe He November 2005 Abstract Introduction / background Data Source Data Manipulation Spatial
More informationNR402 GIS Applications in Natural Resources
NR402 GIS Applications in Natural Resources Lesson 1 Introduction to GIS Eva Strand, University of Idaho Map of the Pacific Northwest from http://www.or.blm.gov/gis/ Welcome to NR402 GIS Applications in
More informationThe Road to Data in Baltimore
Creating a parcel level database from high resolution imagery By Austin Troy and Weiqi Zhou University of Vermont, Rubenstein School of Natural Resources State and local planning agencies are increasingly
More informationHandling Raster Data for Hydrologic Applications
Handling Raster Data for Hydrologic Applications Prepared by Venkatesh Merwade Lyles School of Civil Engineering, Purdue University vmerwade@purdue.edu January 2018 Objective The objective of this exercise
More informationBSEN 6220 GIS LAB #5
BSEN 6220 GIS LAB #5 Dustin Till 10/29/2014 INTRODUCTION: ArcGIS is a geographic information system for working with maps and geographic information. It is used for creating and using maps, compiling geographic
More informationGIS CONCEPTS ARCGIS METHODS AND. 3 rd Edition, July David M. Theobald, Ph.D. Warner College of Natural Resources Colorado State University
GIS CONCEPTS AND ARCGIS METHODS 3 rd Edition, July 2007 David M. Theobald, Ph.D. Warner College of Natural Resources Colorado State University Copyright Copyright 2007 by David M. Theobald. All rights
More informationWorking with Digital Elevation Models in ArcGIS 8.3
Working with Digital Elevation Models in ArcGIS 8.3 The homework that you need to turn in is found at the end of this document. This lab continues your introduction to using the Spatial Analyst Extension
More informationMidterm Exam. 2) What are three important things to keep in mind when selecting colors for a map? (5)
Midterm Exam File Name: LASTNAME_Midterm.pdf Make sure to include your Name, UWNetID, course number, and quarter and year, on the document (points will be deducted if you have not identified your work).
More informationSpatial Effects on Current and Future Climate of Ipomopsis aggregata Populations in Colorado Patterns of Precipitation and Maximum Temperature
A. Kenney GIS Project Spring 2010 Amanda Kenney GEO 386 Spring 2010 Spatial Effects on Current and Future Climate of Ipomopsis aggregata Populations in Colorado Patterns of Precipitation and Maximum Temperature
More informationLaboratory Exercise X Most Dangerous Places to Live in the United States Based on Natural Disasters
Brigham Young University BYU ScholarsArchive Engineering Applications of GIS - Laboratory Exercises Civil and Environmental Engineering 2016 Laboratory Exercise X Most Dangerous Places to Live in the United
More informationSpatial Analysis with Web GIS. Rachel Weeden
Spatial Analysis with Web GIS Rachel Weeden Agenda Subhead goes here Introducing ArcGIS Online Spatial Analysis Workflows Scenarios Other Options Resources ArcGIS is a Platform Making mapping and analytics
More informationMassachusetts Institute of Technology Department of Urban Studies and Planning
Massachusetts Institute of Technology Department of Urban Studies and Planning 11.520: A Workshop on Geographic Information Systems 11.188: Urban Planning and Social Science Laboratory GIS Principles &
More informationHydroacoustic survey and bathymetric map creation for Brant Lake, New York
Hydroacoustic survey and bathymetric map creation for Brant Lake, New York Holly A. Waterfield CLM 1 INTRODUCTION Brant Lake is located in Warren County, New York, within the bounds of the Adirondack Park.
More informationLearning ArcGIS: Introduction to ArcCatalog 10.1
Learning ArcGIS: Introduction to ArcCatalog 10.1 Estimated Time: 1 Hour Information systems help us to manage what we know by making it easier to organize, access, manipulate, and apply knowledge to the
More informationLauren Jacob May 6, Tectonics of the Northern Menderes Massif: The Simav Detachment and its relationship to three granite plutons
Lauren Jacob May 6, 2010 Tectonics of the Northern Menderes Massif: The Simav Detachment and its relationship to three granite plutons I. Introduction: Purpose: While reading through the literature regarding
More informationGIS 2010: Coastal Erosion in Mississippi Delta
1) Introduction Problem overview To what extent do large storm events play in coastal erosion rates, and what is the rate at which coastal erosion is occurring in sediment starved portions of the Mississippi
More informationSpatial Units (Level 1)
Spatial Units (Level 1) Project: Advancing the SEEA Experimental Ecosystem Accounting Overview: Spatial Units 1. Learning objectives 2. Review of Level 0 (5m) 3. Level 1 (Compilers): Presentation & group
More informationA GIS-based Approach to Watershed Analysis in Texas Author: Allison Guettner
Texas A&M University Zachry Department of Civil Engineering CVEN 658 Civil Engineering Applications of GIS Instructor: Dr. Francisco Olivera A GIS-based Approach to Watershed Analysis in Texas Author:
More informationRESERVE DESIGN INTRODUCTION. Objectives. In collaboration with Wendy K. Gram. Set up a spreadsheet model of a nature reserve with two different
RESERVE DESIGN In collaboration with Wendy K. Gram Objectives Set up a spreadsheet model of a nature reserve with two different habitats. Calculate and compare abundances of species with different habitat
More informationGeodatabase Management Pathway
Geodatabase Management Pathway Table of Contents ArcGIS Desktop II: Tools and Functionality 3 ArcGIS Desktop III: GIS Workflows and Analysis 6 Building Geodatabases 8 Data Management in the Multiuser Geodatabase
More informationPredictive Modelling of Ag, Au, U, and Hg Ore Deposits in West Texas Carl R. Stockmeyer. December 5, GEO 327G
Predictive Modelling of Ag, Au, U, and Hg Ore Deposits in West Texas Carl R. Stockmeyer December 5, 2013 - GEO 327G Objectives and Motivations The goal of this project is to use ArcGIS to create models
More informationLecture 9: Geocoding & Network Analysis
Massachusetts Institute of Technology - Department of Urban Studies and Planning 11.520: A Workshop on Geographic Information Systems 11.188: Urban Planning and Social Science Laboratory Lecture 9: Geocoding
More informationWithin this document, the term NHDPlus is used when referring to NHDPlus Version 2.1 (unless otherwise noted).
Exercise 7 Watershed Delineation Using ArcGIS Spatial Analyst Last Updated 4/6/2017 Within this document, the term NHDPlus is used when referring to NHDPlus Version 2.1 (unless otherwise noted). There
More informationGeodatabase An Overview
Federal GIS Conference February 9 10, 2015 Washington, DC Geodatabase An Overview Ralph Denkenberger - esri Session Path The Geodatabase - What is it? - Why use it? - What types are there? Inside the Geodatabase
More informationRaster Spatial Analysis Specific Theory
RSATheory.doc 1 Raster Spatial Analysis Specific Theory... 1 Spatial resampling... 1 Mosaic... 3 Reclassification... 4 Slicing... 4 Zonal Operations... 5 References... 5 Raster Spatial Analysis Specific
More informationTeaching GIS for Land Surveying
Teaching GIS for Land Surveying Zhanjing (John) Yu Evergreen Valley College, San Jose, California James Crossfield California State University at Fresno, Fresno California 7/13/2006 1 Outline of the Presentation
More informationGIS Solutions in Natural Resource Management: Balancing the Technical- Political Equation
GIS Solutions in Natural Resource Management: Balancing the Technical- Political Equation Stan Morain, Editor SUB Gdttlngen 7 208 520 309 98 A14447 0NW0RD PRESS V? % \
More informationSpatial units (Levels 0 and 1)
Spatial units (Levels 0 and 1) Project: Advancing the SEEA Experimental Ecosystem Accounting Overview: Spatial units 1. Learning objectives 2. Level 1: Presentation & group exercise Spatial units (10m
More information06/04/2015. Overview: Spatial units. Advancing the SEEA Experimental Ecosystem Accounting Spatial units (Level 1)
Advancing the SEEA Experimental Ecosystem Accounting Spatial units (Level 1) Advancing the SEEA-EEA Project Overview: Spatial units 1. Learning objectives 2. Level 1: Presentation & group exercise Spatial
More information