Digitization of Trail Network Using Remotely-Sensed Data in the CFB Suffield National Wildlife Area

Similar documents
Spatial Analyst is an extension in ArcGIS specially designed for working with raster data.

Remote Sensing. The following figure is grey scale display of SPOT Panchromatic without stretching.

Visualizing a Pixel. Simulate a Sensor s View from Space. In this activity, you will:

Land Cover Analysis to Determine Areas of Clear-cut and Forest Cover in Olney, Montana. Geob 373 Remote Sensing. Dr Andreas Varhola, Kathry De Rego

Module 11 Digital image processing

Land Cover Type Changes Related to. Oil and Natural Gas Drill Sites in a. Selected Area of Williams County, ND

APCAS/10/21 April 2010 ASIA AND PACIFIC COMMISSION ON AGRICULTURAL STATISTICS TWENTY-THIRD SESSION. Siem Reap, Cambodia, April 2010

Remote Sensing And Gis Application in Image Classification And Identification Analysis.

Urban Classification of Metro Manila for Seismic Risk Assessment using Satellite Images

GE 113 REMOTE SENSING

Figure 3: Map showing the extension of the six surveyed areas in Indonesia analysed in this study.

SEMI-SUPERVISED CLASSIFICATION OF LAND COVER BASED ON SPECTRAL REFLECTANCE DATA EXTRACTED FROM LISS IV IMAGE

Enhancement of Multispectral Images and Vegetation Indices

Landsat 8 TIR Bands 10 and 11 Temperature Comparisons

Building Damage Mapping of the 2006 Central Java, Indonesia Earthquake Using High-Resolution Satellite Images

Detecting artificial areas inside reference parcels. A technique to assist the evaluation of non-eligibility in agriculture

ArcGIS Pro: What s New in Analysis

Mod. 2 p. 1. Prof. Dr. Christoph Kleinn Institut für Waldinventur und Waldwachstum Arbeitsbereich Fernerkundung und Waldinventur

Using Color-Infrared Imagery for Impervious Surface Analysis. Chris Behee City of Bellingham Planning & Community Development

Aerial photography and Remote Sensing. Bikini Atoll, 2013 (60 years after nuclear bomb testing)

Chapter 8. Using the GLM

Mapping Open Water Bodies with Optical Remote Sensing

Detecting and Mapping Invasive Phragmites australis in the Coastal Great Lakes with ALOS PALSAR Imagery

Project summary. Key findings, Winter: Key findings, Spring:

!!!! Remote Sensing of Roads and Highways in Colorado

GeoBase Raw Imagery Data Product Specifications. Edition

Using QuickBird Imagery in ESRI Software Products

Terrain Modeling with ArcView GIS

GE 113 REMOTE SENSING. Topic 7. Image Enhancement

How to Access Imagery and Carry Out Remote Sensing Analysis Using Landsat Data in a Browser

Managing and Monitoring Intertidal Oyster Reefs with Remote Sensing in Coastal South Carolina

Quantifying Change in. Quality Effects on a. Wetland Extent & Wetland. Western and Clark s Grebe Breeding Population

The effects of uncertainty in forest inventory plot locations. Ronald E. McRoberts, Geoffrey R. Holden, and Greg C. Liknes

An Analysis of Aerial Imagery and Yield Data Collection as Management Tools in Rice Production

GROUND CONTROL SURVEY REPORT

General report format, ref. Article 12 of the Birds Directive, for the report

Digital Image Processing

Preparing Remote Sensing Data for Natural Resources Mapping (image enhancement, rectifications )

Application of Satellite Image Processing to Earth Resistivity Map

An Introduction to Geomatics. Prepared by: Dr. Maher A. El-Hallaq خاص بطلبة مساق مقدمة في علم. Associate Professor of Surveying IUG

Spectral Signatures of Tombs and their Classification*

REMOTE SENSING OF RIVERINE WATER BODIES

TEMPORAL ANALYSIS OF MULTI EPOCH LANDSAT GEOCOVER IMAGES IN ZONGULDAK TESTFIELD

Overview of how remote sensing is used by the wildland fire community.

Introduction. Introduction. Introduction. Introduction. Introduction

Exercise 4-1 Image Exploration

Example of Analysis of Yield or Landsat Data Based on Assessing the Consistently Lowest 20 Percent by Using

The New Rig Camera Process in TNTmips Pro 2018

Costal region of northern Peru, the pacific equatorial dry forest there is recognised for its unique endemic biodiversity

In late April of 1986 a nuclear accident damaged a reactor at the Chernobyl nuclear

of Stand Development Classes

High Resolution Multi-spectral Imagery

This week we will work with your Landsat images and classify them using supervised classification.

Combining Images for SNR improvement. Richard Crisp 04 February 2014

Module 3 Introduction to GIS. Lecture 8 GIS data acquisition

DEM GENERATION WITH WORLDVIEW-2 IMAGES

A COMPARISON OF COVERTYPE DELINEATIONS FROM AUTOMATED IMAGE SEGMENTATION OF INDEPENDENT AND MERGED IRS AND LANDSAT TM IMAGE-BASED DATA SETS

Application of Satellite Imagery for Rerouting Electric Power Transmission Lines

GE 113 REMOTE SENSING

DIFFERENTIAL APPROACH FOR MAP REVISION FROM NEW MULTI-RESOLUTION SATELLITE IMAGERY AND EXISTING TOPOGRAPHIC DATA

Image Fusion. Pan Sharpening. Pan Sharpening. Pan Sharpening: ENVI. Multi-spectral and PAN. Magsud Mehdiyev Geoinfomatics Center, AIT

Note: Some squares have continued to be monitored each year since the 2013 survey.

Detecting and Mapping Invasive Phragmites australis in the coastal Great Lakes with ALOS PALSAR imagery

Material analysis by infrared mapping: A case study using a multilayer

Wetlands Investigation Utilizing GIS and Remote Sensing Technology for Lucas County, Ohio: a hybrid analysis.

ISAE - Institute for Studies and Economic Analyses

Saturation And Value Modulation (SVM): A New Method For Integrating Color And Grayscale Imagery

EXAMPLES OF OBJECT-ORIENTED CLASSIFICATION PERFORMED ON HIGH-RESOLUTION SATELLITE IMAGES

A map says to you, 'Read me carefully, follow me closely, doubt me not.' It says, 'I am the Earth in the palm of your hand. Without me, you are alone

Image Extraction using Image Mining Technique

Using Soil Productivity to Assess Agricultural Land Values in North Dakota

Present and future of marine production in Boka Kotorska

On spatial resolution

Documenting Land Cover and Vegetation Productivity Changes in the NWT using the Landsat Satellite Archive

DISTINGUISHING URBAN BUILT-UP AND BARE SOIL FEATURES FROM LANDSAT 8 OLI IMAGERY USING DIFFERENT DEVELOPED BAND INDICES

Ground Truth for Calibrating Optical Imagery to Reflectance

Image interpretation and analysis

University of Technology Building & Construction Department / Remote Sensing & GIS lecture

Unsupervised Classification

Section 6.4. Sampling Distributions and Estimators

Land cover change methods. Ned Horning

GEORGIA WETLANDS TOOL

REMOTE SENSING. Topic 10 Fundamentals of Digital Multispectral Remote Sensing MULTISPECTRAL SCANNERS MULTISPECTRAL SCANNERS

Artificial Neural Network Model for Prediction of Land Surface Temperature from Land Use/Cover Images

Raster is faster but vector is corrector

An investigation of the Eye of Quebec. by means of PCA, NDVI and Tasseled Cap Transformations

GIS Data Collection. Remote Sensing

PROCEEDINGS - AAG MIDDLE STATES DIVISION - VOL. 21, 1988

COMPARISON OF INFORMATION CONTENTS OF HIGH RESOLUTION SPACE IMAGES

GEOGRAPHIC MODELLING AND ANALYSIS

APPENDIX A Vernal Field Office Best Management Practices for Raptors and Associated Habitats

Assessment of Spatiotemporal Changes in Vegetation Cover using NDVI in The Dangs District, Gujarat

Our Quality Promise WHITE PAPER

Remote Sensing Part 3 Examples & Applications

Cross-Talk in the ACS WFC Detectors. II: Using GAIN=2 to Minimize the Effect

GST 101: Introduction to Geospatial Technology Lab Series. Lab 6: Understanding Remote Sensing and Aerial Photography

High Fidelity 3D Reconstruction

Modeling Nightscapes of Designed Spaces Case Studies of the University of Arizona and Virginia Tech Campuses

INFORMATION CONTENT ANALYSIS FROM VERY HIGH RESOLUTION OPTICAL SPACE IMAGERY FOR UPDATING SPATIAL DATABASE

Remote Sensing in an

Transcription:

Digitization of Trail Network Using Remotely-Sensed Data in the CFB Suffield National Wildlife Area Brent Smith DLE 5-5 and Mike Tulis G3 GIS Technician Department of National Defence 27 March 2007

Introduction The Canadian Forces Base Suffield National Wildlife Area (NWA) contains an extensive trail network created naturally by wildlife, grazing cattle, and artificially by oil and gas exploitation. The trail network has never been monitored in the past, and experience shows that increased human activity in NWA has a direct impact on expansion of the trail network. This is a concern from a wildlife conservation perspective because such trails cause habitat fragmentation. Objectives The objectives of this project were twofold: 1. to create a baseline digital database of trails using remotely-sensed imagery in a GIS system, to allow for year-to-year monitoring; and 2. to determine whether trails could be manually classified into distinct categories, based on pixel brightness and trail width characteristics identified from remotely-sensed imagery. Methods GIS Digitizing The NWA trail network database was created by heads-up digitizing 2005 SPOT ortho-rectified panchromatic imagery (2.5 meter resolution) in ESRI ArcGIS 8.3. In addition to digitizing the spatial extent of the network, trails were manually classified into 5 categories in an unsupervised approach; no pre-existing knowledge of ground-based trail characteristics was used to classify them. Classification was carried out using pixel brightness values (based on unprocessed 8-bit digital number/dn) and width of visible trail--the combination of the trail s width and brightness resulted in respective trail classification. After the trails were digitized and classified, ground-truthing was conducted in the Spring of 2006 to verify the classification accuracy. GIS Analysis In order to identify trends between trail categories and SPOT imagery, zonal statistics were created for each trail type using the ArcGIS Spatial Analyst extension. Trail network polylines were buffered by 2.5 m. Resulting buffered polygons were analyzed against 2005 imagery to identify and correlate SPOT DN values with each trail category. Ground-truthing Analysis In order to determine whether heads-up digitizing is a suitable method for digitizing and categorizing trails, a ground-truthing exercise was carried out in May 2006. Ground-based assessment included the field identification and characterization of digitized trail types, and an accuracy assessment. Initially, a field reconnaissance was employed to identify general characteristics of each of the categories identified by heads-up digitizing.

Several days of fieldwork were then required to collect information from 150 randomly-selected trail sites, stratified across all trail types. To minimize potential bias, GPS locations for each site did not include the predicted trail category; this was confirmed only after all of the field data had been collected. Field data included the collection of digital photographs, trail measurements (including rut widths, trail topography--the depth of the rut relative to surrounding areas--percentage of rut vegetative cover, and rut soil exposure. Results GIS Digitizing Based on the completed database, the following was digitized for each trail category: GIS Analysis a. category 1: 113.6 km; b. category 2: 271.8 km; c. category 3: 884.7 km; d. category 4: 370.8 km; and e. category 5: 49.0 km. Based on zonal statistics of raw DN values of 2005 SPOT panchromatic imagery underlying each trail type, there are insignificant differences between pixel brightness of each trail type (Table 1). Table 1 Trail Categories and DN attributes CATEGORY MIN DN MAX DN RANGE MEAN DN STANDARD DEVIATION MEDIAN 1 50 147 97 90.0 7.9 89 2 0 160 160 89.7 8.5 89 3 0 171 171 90.7 10.3 89 4 51 191 140 94.9 12.1 93 5 0 181 181 96.1 31.7 99 Ground-truthing Based on pixel brightness and trail width characteristics, five major categories were identified and described as follows: a. Category #1 represents barely visible trails, minor continuous differences in pixel brightness in comparison to surrounding area, dark pixels with low pixel value (for example animal/cattle trails; b. Category #2 represents minor trails including heavily vegetated trails or pipeline in final stages of recovery, or multiple animal trails;

c. Category #3 represents moderate size trails with brighter pixels, interspersed with a few very bright pixels (visible, distinct trails); d. Category #4 represents wider trails with high pixel brightness values, including heavily-used trails with high soil exposure, poorly vegetated pipelines, or pipeline and trail combinations; and e. Category #5 represents trails with the highest pixel brightness values, and are not distinguishable from major roads except that trail widths are narrower. Accuracy Assessment Based on the accuracy assessment (Table 2), overall classification accuracy is 71%, with user accuracy ranging from 57 to 80%. Table 2 Accuracy Assessment of Digitized Trails Producer's Accuracy - - - - - > Ground Condition (Field Observation) CATEGORY 1 2 3 4 5 Total 1 24 4 2 30 User's Accuracy 2 7 18 5 30 Digitized 3 2 2 25 1 30 4 2 3 23 2 30 V 5 13 17 30 Total 35 24 35 37 19 150 Producer's Accuracy User's Accuracy Overall Accuracy = 71% Categ 1 = 24 / 35 = 69% Categ 1 = 24 / 30 = 80% Categ 2 = 18 / 24 = 75% Categ 2 = 18 / 30 = 60% Statistic Value Categ 3 = 25 / 35 = 71% Categ 3 = 25 / 30 = 83% N = 150 Categ 4 = 23 / 37 = 62% Categ 4 = 23 / 30 = 77% Part A = 107 Categ 5 = 17 / 19 = 89% Categ 5 = 17 / 30 = 57% Part B = 4500 Khat = 64.17% OMISSION - - - - - > Ground Condition (Field Observation) CATEGORY 1 2 3 4 5 Total 1 24 4 2 30 2 7 18 5 30 COMMISSION Digitized 3 2 2 25 1 30 4 2 3 23 2 30 5 13 17 30 V Total 35 24 35 37 19 150 Omission Commission Overall Accuracy = 71% Categ 1 = 4 + 2 / 30 = 20% Categ 1 = 7+2+2/30 = 37% Categ 2 = 7 + 5 / 30 = 40% Categ 2 = 4+2/30 = 20% Statistic Value Categ 3 = 2+2+1/30 = 17% Categ 3 = 2+5+3/30 = 34% N = 150 Categ 4 = 2+3+2/30 = 23% Categ 4 = 1+13/30 = 47% Part A = 107 Categ 5 = 13 / 30 = 43% Categ 5 = 2 / 30 = 7% Part B = 4500 Khat = 64.17%

Discussion Based on the combined attributes of field observations and GIS analysis, each trail category is described below, using trail width, soil micro-topography, vegetation cover, soil exposure of a single rut, and raw DN characteristics, including mean, and standard deviation (Std). Category 1: Width 20-40 cm Topography 0-5 cm Vegetation Cover 80-100 % (with exception of animal trails 5-15) Soil Exposure 0-20 % (with exception of animal trails 85+ %) Mean DN +/- Std 90.0 +/- 7.9 Note: Cattle trails have significantly greater soil exposure than access trails, but significantly narrower disturbance. Nonetheless, pixel brightness values are indistinguishable between animal trails and access trails using SPOT imagery. Category 2: Width 30-60 cm Topography 2-20 cm Vegetation Cover 70-100 % Soil Exposure 0-20 % Mean DN +/- Std 89.7 +/- 8.5 Note: Multiple parallel cattle trails could be digitized as category 2 trail, because of the increased soil exposure and higher pixel brightness values. Category 3: Width 35-75 cm Topography 5-30 cm Vegetation Cover 10-60 % Soil Exposure 10-100 % Mean DN+/- Std 90.7 +/- 10.3 Note: Average size and the most common trail. Category 4: Width 90-200 cm Topography 5-35 cm Vegetation Cover 0-50 % Soil Exposure 50-100 % Mean DN+/- Std 94.9 +/- 12.1 Note: Wide, eroded trail with little vegetation.

Category 5: Width 80-200 cm Topography 0-60 cm Vegetation Cover 0-15 % Soil Exposure 80-100 % Mean DN +/- Std 96.0 +/- 31.7 Note: Typical category 5 trail has very wide ruts with no vegetation. Conclusions The field component and accuracy assessment of this study positively supported the objectives of this project. We conclude that manually digitizing and classifying trails based on remote-sensing imagery is effective and sufficiently accurate for creating and monitoring the trail network within the NWA. Because of the subtle and insignificant differences in pixel brightness between each of the trail categories, we conclude computer-based classification is not possible. The separation between category 4 and category 5, and separation between category 1 and category 2 trails during heads up digitizing is difficult to the inexperienced eye. Because of subtle differences in pixel brightness values, sufficient training is required to ensure that these categories can be accurately classified during heads-up digitizing.

Field Examples The following provide examples of field-based descriptions of each category. Category 1: Cattle Trail Grid 514649 5576018 Category 1 Pictures No. 134-3446 to 47 Width (cm) 30-40 Topography (cm) 0-5 Vegetation (%) 0-15 Soil Exp. (%) 80-100 Heavily vegetated trail Grid 515832 5574770 Category 1 Pictures No. 134-3454 to 55 Width (cm) 220-280 Topography (cm) 0-5 Vegetation (%) 80-100 Soil Exp. (%) 10-20 Completely recovered trail Grid 516376 5574735 Category 1 Pictures No. 134-3457 Width (cm) 20-25 Topography (cm) 0-5 Vegetation (%) 95+ Soil Exp. (%) 5

Category 2: Heavily vegetated trail Grid 514658 5576301 Category 2 Pictures No. 134-3441 to 43 Width (cm) 35-60 Topography 0-10 (cm) Vegetation (%) 75-95 Soil Exp. (%) 0-20 Heavily vegetated trail Grid 514578 557630 Category 2 Pictures No. 134-3444 to 45 Width (cm) 30-40 Topography (cm) 0-5 Vegetation (%) 90-100 Soil Exp. (%) 0-10 Heavily vegetated trail Grid 516793 557375 Category 2 Pictures No. 134-3458 Width (cm) 45-55 Topography (cm) 5-20 Vegetation (%) 70-90 Soil Exp. (%) 10-20

Category 3: Moderately-eroded trail Grid 514689 5576368 Category 3 Pictures No. 134-3422 to 25 Width (cm) 60+ Topography (cm) 10-30 Vegetation (%) 10-60 Soil Exp. (%) 40-100 Moderately-eroded trail Grid 515177 557594 Category 3 Pictures No. 134-3450 to 52 Width (cm) 45-65 Topography (cm) 5-20 Vegetation (%) 10-50 Soil Exp. (%) 20-95 Eroded trail and pipeline Grid 515849 5575092 Category 3 Pictures No. 134-3453 Width (cm) 40-60 Topography (cm) 10-20 Vegetation (%) 30-60 Soil Exp. (%) 10-35

Category 3 to 4 transition: Heavily eroded trail on side slope Grid 520908 5570518 Category bottom 3 / top 4 Pictures No. 134-3472 Width (cm) 35 / 90 Topography (cm) 5-10 / 10-20 Vegetation (%) 30-60 / 0-20 Soil Exp. (%) 30-55 / 80-100 Category 4: Heavily eroded trail Grid 514708 5576423 Category 4 Pictures No. 134-3426 to 35 Width (cm) 100-210 Topography (cm) 15-35 Vegetation (%) 0-30 Soil Exp. (%) 70-100 Eroded pipeline right-away Grid 521755 5572425 Category 4 Pictures No. 134-3469 Width (cm) 50-200 Topography (cm) 5-25 Vegetation (%) 0-50 Soil Exp. (%) 50-100

Category 5: High bare ground exposure Grid 514739 5576644 Category 5 Pictures No. 134-3438 Width (cm) 80+ Topography (cm) 20-60 Vegetation (%) 0-15 Soil Exp. (%) 80-100 High bare ground exposure Grid 521131 5569947 Category 5 Pictures No. 134-34 Width (cm) 125 Topography (cm) 5 Vegetation (%) 0 Soil Exp. (%) 100