Lab 6 Density Analysis

Raster based density calculations represent a type of interpolation.The ArcGIS manual overview of the Density toolset states, "calculating density, you are in a sense spreading the values (of the input) out over a surface". The section of the ArcGIS manual on Understanding Density Analysis states, "[t]he Density tool distributes a measured quantity of an input point layer throughout a landscape to produce a continuous surface." Calculating density is not the same as generating topography from point elevations, but both density and topographic rasters represent an array of discrete points as a field of values. Thus, in a very general sense one can think of Density as a type of spatial interpolation.

ArcGIS provides three Density tools: Kernel, Line, and Point. Kernel densities produce tapered surfaces. Line and Point densities produce magnitude per unit area within the radius of each cell. The two figures below illustrate Point Density and Kernel Density surfaces generated from the same feature classes, that are in this case archaeological sites. Densities can be generated on just about any kind of point or line feature.

 Point Density (90m cell 2.5k radius) Kernel Density (90m cell 2.5k radius)

There are three important variables to consider when generating a density raster.

• Density type: Should the density raster have tapering values (kernel) or constant values (point density)?
• Output raster cell size: What is the size of the cells in the output raster. It is often useful to set the spatial resolution of the output density raster to the same as the digital elevation model (DEM) as this facilitates raster math.
• Search radius: This determines the size of the search neighborhood. In the literature, it is often referred to as "r". To avoid having densities of one entity, consider determining the distance of the minimum nearest feature and setting "r" slightly higher. When considering use areas around a site, a reasonable search radius might be the distance an individual walks in half a day.
• Population field: Should each input point from the feature class be weighted equally or is there a population field that can give greater weight to features that are larger. In the context of archaeological sites, area might make a useful population field. When mapping small objects like flakes, population fields may not be particularly useful. Whether or not a population field is useful depends on the analytical problem.

By default, density will use the extent of the input feature class. This means that points on the edges of the distribution will have truncated densities.

In most cases it is better to use a study area or some other larger bounding polygon to set the limits of the density operation. This is achieved by selecting the appropriate extent under the Environments button and selecting the Processing Extent option. Under Processing Extent, indicate the appropriate feature to employ and click OK.

Beware, that altering the Processing Extent can cause changes in the Output cell size and Search radius. Be sure that these values are correct or the density computation may not be useful. The same cell size as the input DEM is probably a prudent choice. The search radius is best defined by the structure of the data set or research question. For a set of sites on relatively flat terrain, search radius related to walking times can produce interesting results.

1. For the qocha shapefile run the tool Convert feature to point.
1. Be sure that the attributes transfer to the new point shapefile. Think about whether or not to constrain the points inside the feature or not. In this case, it is probably recommendable to place the point inside the feature.
2. Once the features are converted to points, determine the location of the nearest feature.
1. This can be done using the Generate Near Table tool, with ET GeoWizards, and probably some other tools (please report if you find another solution).
2. If using the Generate Near Table tool, it is advisable to use the Find Only Closest Feature option as this will cut down on the creation of unnecessary information.
3. Both the Input and the Near features should be the point shapefile.
4. Open the resulting table of inter-feature distances and sort the entries by the near distance. What is the closest and what is the greatest inter-feature distance? The units are in meters. Is the distribution normal? Is it possible that some outliers could be identified and eliminated?
5. After considering the results from calculating the nearest feature, run the density analysis tool. As we will be running the tool at least twice, it may be worth using the model builder to run the tool.
6. For the first run of this tool, do not set any population field. We will do this in the next version.
7. Set a logical location for the output raster.
8. Set the output cell size to 90m.