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 2013 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 397917/541861 or 73.44% 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 | 21387/21387
Shrub/grass | 53.87 | 18557/34445
Agriculture | 18.87 | 6711/35571
Coniferous forest | 84.08 | 165908/197320
Broadleaf forest | 73.16 | 184787/252571
Fire scar | 100.00 | 567/567
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 | 93.92 | 21387/22772
Shrub/grass | 15.49 | 18557/119803
Agriculture | 73.73 | 6711/9102
Coniferous forest | 89.81 | 165908/184731
Broadleaf forest | 92.24 | 184787/200322
Fire scar | 11.05 | 567/5131
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)
Process steps performed in ENVI 5.1 to produce classification map:
1. Publicly available annual Landsat Best Available Pixel Composite (LBAPC) datasets for 2012 and 2013 were obtained from the the Department of Geographical Sciences, University of Maryland. Of these, 2012 LBAPC data was used to produce this Land Use and 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 2013 ground-reference data. Normalized Difference Built-up Index (NDBI) was applied to 2013 and 2010 LBAPC data to produce LULC changes from 2013 to 2010. Ground reference data for 2013 was selected by subtracting these changes from the 2010 ABMI data. Ground-reference vegetation types for 2013 include coniferous forest, broadleaf forest, mixed forest, agriculture, grassland, and shrubland. Ground-reference developed areas for 2013 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 2013 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 2013 LBAPC data to produce LULC classification results.
3. ML classification algorithm was applied to the 2013 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 2013 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.