The Cold Lake oil sands area, Township 56 to 69, Range 1 to 11, west of the 4th Meridian, falls within the Lower Athabasca Regional Plan (LARP). As part of Alberta’s Land-use Framework, LARP was developed in 2012 to set the stage for robust growth, vibrant communities and a healthy environment within the region. One of its implementation objectives is to balance the economic development of oil sands and impacts on ecosystem and environment. This is to be achieved through enhanced science-based monitoring for improved characterization of the environment and collection of the information necessary to understand cumulative effects.
This land classification raster dataset is derived from 2008 Landsat imagery. It contains 13 classes: 0 - unclassified, 1 - exposed land/cut blocks/harvested areas, 2 - water bodies, 3 - transitional bare surfaces, 5 - mixed developed areas, 6 - developed areas, 7- shoal, 8 - shrub land, 9 - grassland, 10 - agriculture areas, 11 - coniferous forest, 12 - broad leaf forest, 13 - mixed forest, and 14- fire scar.
The overall accuracy, defined as the total correct pixels (sum of the major diagonal) divided by the total number of pixels in the error matrix N, is 417436/523627 or 79.72% for this dataset.
The Producer accuracy (Table 1), or omission error, is an indication of the probability of a reference pixel being correctly classified. It is calculated as the total number of correct classified pixels in a class divided by total number of pixels in that class.
Table 1 - Producer accuracy results
Class | Producer Accuracy (%) | Pixel Ratio
----- --------------------- -----------
Developed | 100.00 | 27113/27113
Shrub/grass | 43.98 | 19392/44096
Agriculture | 57.05 | 41935/73505
Coniferous forest | 88.91 | 147073/165411
Broadleaf forest | 85.21 | 181923/213502
The User accuracy (Table 2), or commission error, is defined as the total number of correct classified pixels in a class divided by total number of pixels that were classified in that class. It indicates the probability that a pixel classified on the map represents that class on the ground.
Table 2 - User accuracy results
Class | User Accuracy (%) | Pixel Ratio
----- --------------------- -----------
Developed | 89.66 | 27113/30240
Shrub/grass | 30.12 | 19392/64374
Agriculture | 81.84 | 41935/51238
Coniferous forest | 89.52 | 147073/164290
Broadleaf forest | 85.22 | 181923/213485
Accuracy assessment for the fused classification result was performed with validation data (derived from the ground-reference data in step 2) that include developed areas, mixed developed areas, coniferous forest, broadleaf forest, shrub/grass, and agriculture classes. Transient LULC types with mixed forest, exposed land/cutblocks/harvested areas, transitional bare surfaces, and shoal were not used for validation since ground reference data for these classes were not realistic.
developed - mixed developed area (class 5) and developed area ( class 6) were combined to one class
shrub/grass - shrubland (class 8) and grassland (class 9) were combined to one class
agriculture - agricultural area (class 10)
coniferous - coniferous forest (class 11)
broadleaf - broadleaf forest (class 12)
A major component of this dataset is derived from Ducks Unlimited Canada's drained and boreal Enhanced Wetland Classification inventories.
Http://www.capf.ca/pdfs/AWCS%20Draft%20requesting%20Public%20Feedback.pdf
Process steps performed in ENVI 5.1 to produce classification map:
1. Pre-release versions of annual Landsat Best Available Pixel Composite (LBAPC) datasets (1984 to 2012) were obtained from the Natural Resources Canada/Pacific Forestry Centre for testing purpose. Among these datasets, 2008 LBAPC was used to produce this Land Use Land Cover (LULC) classification result.
2. The Alberta Biodiversity Monitoring Institute (ABMI) LULC data for 2010 and the Alberta Department of Energy (ADoE) oil and gas infrastructure data for 2012 were used to produce the 2008 ground-reference data. Normalized Difference Built-up Index (NDBI) was applied to 2008 and 2010 LBAPC data to produce LULC changes from 2008 to 2010. Ground reference data for 2008 was selected by subtracting these changes from the 2010 ABMI data. Ground-reference vegetation types for 2008 include coniferous forest, broadleaf forest, mixed forest, agriculture, grassland, and shrubland. Ground-reference developed areas for 2008 were queried from the ADoE oil and gas infrastructure data. Training data for vegetation types and developed areas were selected from the ground reference data. Training data for exposed land/cutblocks/harvested land and water were selected manually from the 2008 LBAPC data. Areas with cloud, cloud shadow, and noise were discarded from training and validation data. Maximum Likelihood (ML) classification algorithm, Spectral Angle Mapper (SAM) partial unmixing algorithm, and Constraint Energy Minimization (CEM) partial unmixing algorithm utilized the training data with 2008 LBAPC data to produce LULC classification results.
3. ML classification algorithm was applied to the 2008 LBAPC data to produce LULC classification result with coniferous forest, broadleaf forest, mixed forest, agriculture, grassland, shrubland, exposed land/cutblock/harvested land, and water classes.
4. LULC classification result with developed areas, mixed developed areas, and transitional bare surface classes were produced from the 2008 LBAPC data by dividing the CEM partial unmixing result by the SAM partial unmixing result, followed by the K-Means clustering. Post classification techniques (i.e., majority/minority filter, clump, and sieve) were applied on each of these classes to refine the result and to minimize false detections.
5. A fused classification result was produced by combining LULC classification results from step 3 and 4. A buffer with 90 m was applied on the stream channel datasets from Alberta Energy Regulator (AER) and Alberta Merged Wetland Inventory (AMWI) to track yearly changes associated with landslide activities and the appearance or disappearance of sandbars due to water level fluctuation. This buffer zone is assigned to a separate class called shoal.
6. Accuracy assessment for the fused classification result was performed with validation data (derived from the ground-reference data in step 2) that include developed areas, mixed developed areas, coniferous forest, broadleaf forest, shrub/grass, and agriculture classes. Transient LULC types with mixed forest, exposed land/cutblocks/harvested areas, transitional bare surfaces, and shoal were not used for validation since ground reference data for these classes were not realistic.