Applications of Spatial Statistics
Hazard Rating Model This image depicts the probability of observing Armilaria root disease on randomly located sample plots in the Black Hills National Forest of South Dakota. Important varibles used in describing the large-scale spatial variability included elevation, site index, and slope. The small-scale spatial variability was modeled using kriging. The final model accounted for 93% of the spatial variability in the sample data.
Forest Inentory Analysis The four images depict some results of a forest inventory of a management unit in the Ejidi el Largo, which is located in the western part of the state of Chihuahua, Mexico. In the summer of 1998, the management unit was sampled to obtain stand level estimates (top figure in each image). The same information was used in developing a spatial model to describe the variable of interest across the management unit at a 60m resolution (bottom figure in each image). Elevation and various Landsat MSS bands (60m resolution) were used in modeling the large-scale spatial variability, and kriging was used to describe the small-scale spatial variability.
Micro-ecological Units
This image demonstates the application of a multivariate spatial clustering
algorithm for identifying micro-ecological units using forest inventory data
from northern Mexico. In this example, areas with the same color have
similar spatial characteristics with respect to habitat suitability for the
wild turkey, standing wood volume, stand density index, and 10 year diameter growth.
Soil Types
This image depicts the results of a model developed to
predict soil types in the Rocky Mountain National Park,
Colorado. Discriminant analysis was used to develop a set
of discriminant functions for 15 soil types identified from
229 soil pits, strategically located throughout the park.
Variables identified as being important in discriminating
among soil types include such variables as, individual
Landsat TM bands, pH, depth to bedrock, percent lamella, and
percent rockiness. A spatial autoregessive model was used
to spatially predict pH, depth to bedrock, percent lamella,
and percent rockiness, as function of percent slope,
elevation, and various Landsat TM bands. The spatially
predicted surfces of soil characteristics, along with the
Landsat TM bands were passed through the discriminant
functions, one pixel at a time, to produce the final map.
Species Richness
These two images depict the spatial distribution of native and non-native
species on 9,500 ha of the Rocky Mountain National Park, Colorado. The
large-scale spatial variability in the number of non-native species was
modeled as a function of various Landsat TM bands, percent sand, percent
silt, carbon, elevation, and the number of native species. The model as fit using a spatial autoregressie model. There was no small-scale spatial
variability associated with the spatial distribution of non-native species.
The large-scale spatial varibility in the number of native species, percent
sand, percent silt, and carbon was modeled using various Landsat TM bands and
elevation. Kriging was used to model the small-scale spatial variability.
When the maps were developed no data was available for the high elevation
tundra and so the models over estimate both the number of native and
non-native species in this area.
Soil Texture A trend surface model was used to spatially predict percent clay as a function of various Landsat TM bands, radar imagery, and elevation. Kriging as used to describe the small-scale spatial varibaility in percent clay. Percent sand was modeled as a function of percent clay, elevation, aspect, and various Landsat TM bands. The prediction of percent clay and sand were constrained to have the same probability density function as the field data. Percent silt was estimated by subtracting percent clay and sand from 100. Finally, the soil texture map was generated by classifying each pixel of the image into one of 10 soil types based on the amount of sand, silt, and clay.
Rock Outcropping
A linear discriminant function was developed to spatially predict
the location of rock outcropping as a function of various Landsat TM bands and
elevation (middle figure). The model falsely predicted the presence
of rock outcroppings on 4% of the area (bottom figure), and falsely predicted the absence
on 10% of the area, for an overall accuracy of 85.5% (middle figure).
Back to Robin's Homepage