Cost-Effectiveness of Seven Approaches to Map Vegetation Communities — A Case Study from Northern Australia ’ s Tropical Savannas

Vegetation communities are traditionally mapped from aerial photography interpretation. Other semi-automated methods include pixeland object-based image analysis. While these methods have been used for decades, there is a lack of comparative research. We evaluated the cost-effectiveness of seven approaches to map vegetation communities in a northern Australia’s tropical savanna environment. The seven approaches included: (1). aerial photography interpretation, (2). pixel-based image-only classification (Maximum Likelihood Classifier), (3). pixel-based integrated classification (Maximum Likelihood Classifier), (4). object-based image-only classification (nearest neighbor classifier), (5). object-based integrated classification (nearest neighbor classifier), (6). object-based image-only classification (step-wise ruleset), and (7). object-based integrated classification (step-wise ruleset). Approach 1 was applied to 1:50,000 aerial photography and approaches 2–7 were applied to SPOT5 and Landsat5 TM multispectral data. The integrated approaches (3, 5 and 7) included ancillary data (a digital elevation model, slope model, normalized difference vegetation index and hydrology information). The cost-effectiveness was assessed taking into consideration the accuracy and costs associated with each classification approach and image dataset. Accuracy was assessed in OPEN ACCESS Remote Sens. 2013, 5 378 terms of overall accuracy and the costs were evaluated using four main components: field data acquisition and preparation, image data acquisition and preparation, image classification and accuracy assessment. Overall accuracy ranged from 28%, for the image-only pixel-based approach, to 67% for the aerial photography interpretation, while total costs ranged from AU$338,000 to AU$388,180 (Australian dollars), for the pixel-based image-only classification and aerial photography interpretation respectively. The most labor-intensive component was field data acquisition and preparation, followed by image data acquisition and preparation, classification and accuracy assessment.


Introduction
There is an increased requirement for reliable and up-to-date vegetation community information across the globe [1,2].Due to an increasing population, natural resources are compromised as more land is required to accommodate agricultural, residential, industrial and commercial activities.This results in extensive land clearing of native, often remnant vegetation [1][2][3].Without reliable vegetation spatial information, sound decision making for natural resource management and biodiversity conservation is jeopardized.Furthermore, the vegetation datasets that are available are too coarse to make decisions at a local or property management scale.Vegetation community mapping is being practiced in most parts of the world.The numerous methods used continue to evolve and adapt to the diversity of vegetation types found.
There are a number of techniques commonly used to map vegetation communities.The traditional method is aerial photography interpretation (API) and intensive field sampling.Two semi-automated methods include pixel-based image analysis (PBIA) and geographic object-based image analysis (GEOBIA).API is a method that requires manual stereoscopic interpretation of aerial photograph pairs.The use of aerial photography was anticipated to decrease as technological advances improved in the capabilities of new airborne and satellite sensors, however, it continues to be a universal and acceptable method for numerous mapping applications [4,5].API is documented as an expensive, time consuming and inefficient process even though aerial photography is available digitally [6][7][8][9][10].
Semi-automated techniques are becoming increasingly available as a result of developing image processing software.These software packages now offer a realm of pixel and object-based functions for discriminating vegetation communities from commercially available image datasets.Pixel-based classification assigns individual pixels to a class according to their statistical characteristics.Object-based classification first groups or segments neighboring pixels with similar characteristics and then assigns each segment to a class [11].While these methods have been used for decades, there is a lack of research that compares the accuracy and cost-effectiveness of the common approaches applied to various image datasets and spatial scales.Although GEOBIA emerged as early as the 1970s, GEOBIA was not used extensively within the field of remote sensing until recently [12][13][14].Some comparative studies have highlighted strengths and weaknesses of each approach and further recommend a combination of the two can be beneficial [11,15].
Existing studies have compared aspects of land cover and vegetation community maps, such as spatial scale [2], spatial versus spectral resolution [16][17][18], the influence of ancillary data on classification results [6][7][8][19][20][21][22][23][24], and the accuracy of these types of maps [25].However, few studies have presented an analysis of cost or even generic discussion of cost-effectiveness of remote sensing applications.Cost-effectiveness herein refers to the accuracy of thematic mapping products and the costs and time associated to generate the end mapping product.The evaluation of cost appears to be limited by the large number of variables to consider, including survey purpose, hardware availability, size of study area, technical expertise of staff, cost of data acquisition, intended accuracy, and so on.In tropical savanna environments, there are few publications on the cost of evaluations for vegetation community mapping applications.The majority of cost evaluations for vegetation community mapping programs seem to reside internally as unpublished documents [26].Also, mapping programs in some Australian jurisdictions have assessed the costs in artificial situations but are not related so much to actual operational costs that happen in practice [27].
Research that does exist on cost-effectiveness of remote sensing techniques is focused more on tropical coastal habitat mapping [28][29][30].A toolkit developed by United Nations Educational, Scientific and Cultural Organization (UNESCO) is also available for assessing the cost-effectiveness of remote sensing for coastal management [31].Several coastal habitat mapping projects acknowledge it is difficult to place a monetary value on remotely sensed information despite its great utility in coastal systems [29].The frameworks that have been developed for coastal environments have generally identified key cost categories, including set-up costs, field survey costs, image acquisition costs, and costs associated with field data analysis and image processing.A framework has also been developed, initially for coastal environments, to enable spatial data collection needs to be linked to a suitable form of remotely sensed data [30]; however, this framework does not incorporate cost-effectiveness.
Accuracy assessment is another essential constituent of cost-effectiveness and is regarded as a fundamental component of thematic mapping from remotely sensed data.Many problems are encountered when assessing the accuracy of image classification as there is no single, universal and acceptable measure [32][33][34].The most widely used approach is the derivation of confusion or error matrices [8,33,35], and is applied in this study.After generating an error matrix, accuracy assessment metrics, such as overall accuracy, producer and user accuracies, can be derived [8].It has also been customary to report Kappa coefficient which became a standard component of almost every accuracy assessment, however, recent papers by Pontius and Millones [36] and Foody [37] report the redundancy of Kappa as a result of misleading information for practical decision making, and we reflect this premise accordingly.Accuracy assessment of vegetation community maps, and similarly land cover maps generated from remotely sensed data, provides important data quality information [8,32,33,35,38].Given the complexity of remotely sensed classifications, there is a requirement to assess the reliability.Traditionally, API has been accepted to be correct without any form of accuracy assessment [35].However, for semi-automated approaches such as PBIA and GEOBIA, accuracy assessment is common practice.Thus, there is a glaring inconsistency in the application of accuracy assessment applied to various mapping techniques which presents a need for comparative research.We address this gap in knowledge herein by comparing the accuracy of API, PBIA and GEOBIA generated vegetation community thematic maps.
Existing cost-effective frameworks highlight the urgency of cost-effective evaluations if managers are to make appropriate choices when faced with the growing plethora of remote sensing technologies [28], especially where the potential benefits and costs are very large.Project managers prefer the least expensive remote sensing technique that will achieve a given task with an acceptable degree of accuracy [39].This study assesses the cost-effectiveness of seven vegetation community mapping approaches applied to 1:50,000 aerial photography, SPOT5 and Landsat5 TM multi-spectral data.Some of the approaches integrate ancillary data in the classification processes, including a digital elevation model (DEM), DEM-derived slope model, normalized difference vegetation index (NDVI) and hydrology information.Field data were used for each approach, 50% for the classification processes and the remaining 50% for the accuracy assessment.Two spatial scale maps (1:25,000 and 1:100,000) were generated for each approach and image dataset.The cost-effectiveness comprised four main cost components and was assessed in terms of overall accuracy and the costs associated with each approach and image dataset.We provide the first cost-effectiveness study for a tropical savanna environment to map the floristic and structural components of vegetation communities at a property management scale (1:25,000) and regional scale (1:100,000), whilst also comparing numerous techniques and image datasets.Three major comparisons are included in this study in terms of cost-effectiveness: seven mapping approaches, three image datasets, and two spatial resolution outputs.

Study Area
The study area is located on Bullo River Station in the Victoria River District in north western Northern Territory, Australia (Figure 1).The study area covers 530 km 2 and is situated in the Bullo River catchment, representing three broad landform types: rugged sandstone hills and escarpment; low hills, rises and plains; and alluvial plains towards the intertidal fringes of the Bullo and Victoria rivers.These landform types support a range of habitats typical of northern Australia's tropical savannas including a variety of eucalypt communities, riparian zones, paperbark swamps, mangrove communities and saline coastal flats subject to tidal inundation.

Field Data Sampling
The field sampling was conducted over four years (2006 to 2009) and six sampling efforts, where access was achieved by helicopter and four wheel drive vehicle.A systematic approach was used to pre-select sites covering the geographic and environmental range across the study area using a Geographic Information System (GIS).To represent the various vegetation patterns, site selection was based on tonal variation, color and texture of the area as represented in 1:50,000 aerial photography and SPOT5 imagery.Disturbed areas (i.e., recent fire and grazing) were avoided.
Across the study area two site types were sampled: full floristic sites and less detailed sites, called "road notes" (Figure 1). Figure 1 illustrates field sites across Bullo River Station that were used for the multi-variate analysis of this study, including 392 sites and 957 plant species.At each full floristic site, all plant species present in a 20 × 20 m quadrat were recorded with associated structural information (cover, height, and growth form across three strata).Cover was estimated as canopy cover (crowns treated as opaque) for the upper strata and projective foliage cover (PFC-vertical projection of foliage only) for the mid-and ground strata.Mean height and range were measured for species greater than 2 m tall and visually estimated for those less than 2 m (for species greater than 1% cover).Road notes were qualitative and included a description of the dominant species and an estimate of the structural formation (including cover and height) for homogeneous vegetation patterns.Within the study area, we sampled 137 full floristic sites and 104 road notes, 241 sites in total.Fifty percent of the field dataset were for training the seven mapping approaches and the remaining 50% were reserved for the accuracy assessment.Sampling intensity was dependent on accessibility, funding and resources.Detailed field sampling methods for the full floristic and road notes sites is provided by Lewis et al. [40] and is based on Australia's national standards for vegetation field sampling and vegetation classification [40,41].

Field Data Analysis and Vegetation Community Classification
Multi-variate routines were applied to a subset (upper strata) of the full floristic dataset.Species contributing less than 0.1% cover were of low abundance and not diagnostic of the vegetation communities in the study and were removed for the analysis.A square root transformation was then applied to the upper strata dataset [40].The most commonly used similarity coefficient (Bray-Curtis) was calculated and multi-dimensional scaling plots were used as a visual aid to remove 39 outlier sites.A combination of multi-variate analysis and intuitive classification identified 22 discrete and mappable vegetation communities across the study area (Appendix 1).The Similarity of Percentages (SIMPER) procedure was then used to discern species typical of the vegetation communities.SIMPER was also used to rank species in order of their relative contributions to determine community patterns for each floristic group [40].Vegetation attributes were summarized to construct vegetation community descriptions [40] and were described using the National Vegetation Information System (NVIS) Information Hierarchy Level VI-sub-association, the highest level of detail floristically and structurally [42].Finally, all field sites were assigned a vegetation community number prior to API, PBIA and GEOBIA.

Image Data Acquisition and Pre-Processing
Aerial photography was captured in May 2006 at a scale of 1:50,000 covering the extent of Bullo River Station (Figure 1).Differential Global Positioning System (GPS) centers and exterior orientation points from Applanix, color contact prints and photo scans at 15 micron resolution were also sourced.Every second digital aerial photograph was ortho-rectified using an exterior orientation method.Eight fiducial points were selected with a root mean square error of below 0.10 units.The output coordinate system was Geocentric Datum of Australia (GDA94), Map Grid of Australia (MGA52) eastings/northings at a scale of 1:50,000 (0.00002) and cell size 2 × 2 m pixels using a nearest neighbor algorithm.The ortho-rectified aerial photos were mosaiced and color balanced with a blue haze filter [43].
The SPOT5 (10 m pixels) and Landsat5 TM (25 m pixels) scenes were also captured in May 2006 and supplied by Geoscience Australia and Geoimage respectively.The scenes were geometrically corrected to minimize distortions by using 24 ground control points.The SPOT5 imagery was ortho-rectified by Geoscience Australia and the Landsat5 TM scene was geometrically corrected using a standard technique developed by the Queensland Government Remote Sensing Centre (http://www.derm.qld.gov.au/slats/index.html).

Training Areas and Ancillary Data
Training areas were created from 50% of the field data full floristic and road note sites.A systematic method was used to differentiate which field sites to use for the classification of the API, PBIA and GEOBIA approaches, and which ones for the accuracy assessment.Odd number sites were selected for the classifications and even site numbers for the accuracy assessment.There were two vegetation communities (13 and 30) which did not have field sites for the accuracy assessment, thus additional training areas were delineated based on interpreter knowledge of the study area.These communities were uncommon and small in extent across the study; therefore were under-sampled.
Four ancillary datasets were used in the integrated PBIA, NN-GEOBIA and step-wise GEOBIA approaches, including NDVI, DEM, DEM-derived slope model, and the 1:250,000 hydrology layer.The NDVI was calculated for Landsat5 TM and SPOT5 multi-spectral scenes according to the equation NIR-Red/NIR+Red, where NIR and R are the pixel values for the Near-Infrared and Red bands, respectively.The DEM was acquired through the Australian Defence Force as a 30 × 30 m post spacing, 10 m contour file and currency 1990-1994.These ancillary datasets were subset to the study area.The hydrology polyline shapefile was sourced from the Water Resources Survey of Bullo River Station [44].From the original hydrology dataset, two shapefiles were produced by buffering the major rivers and creeks at 80 m, and minor creeks at 30 m to capture the extent of the two riparian vegetation communities present in the study area (21 and 4).

Classification Methods
The cost-effectiveness of seven vegetation community mapping approaches were assessed, including: (1).aerial photography interpretation, (2).pixel-based image-only classification, (3).pixel-based integrated classification, (4).object-based image-only nearest neighbor classification, (5).object-based integrated nearest neighbor classification, (6).object-based image-only step-wise classification, and (7).object-based integrated step-wise classification.The first approach was applied to 1:50,000 aerial photography and the remaining ones (approaches 2 to 7) to both SPOT5 and Landsat5 TM multi-spectral data.The integrated approaches incorporated ancillary data; whereas the image-only approaches only used the image datasets.Field data were used for each approach, 50% for the classification processes and the remaining 50% for the accuracy assessment.Two spatial scale maps (1:25,000 and 1:100,000) were generated for each approach and image dataset.

Aerial Photography Interpretation
The aerial photography stereo pairs were examined under a stereoscope to delineate vegetation communities.Line-work was digitized as a polyline shapefile (GDA94, decimal degrees) using the mosaic as an interpretive base in a GIS.The spatial scale was set to 1:10,000 for line-work digitizing.The polyline dataset was smoothed using a smooth polylines algorithm and converted to a polygon shapefile [43].Preliminary map attributes were assigned to the polygons of the original 1:10,000 polyline dataset and updated once the final vegetation community groups were determined based on 76 full floristic and 50 road notes and visual interpretation.Topography was also evaluated to define landform and land surface characteristics.Polygons less than 0.25 ha were eliminated to create the 1:25,000 thematic map and 4 ha for 1:100,000 [43].

Pixel-Based Image Analysis
Two pixel-based methods were applied to the Landsat5 TM and SPOT5 image datasets: image-only approach (image and field data), and integrated approach using ENVI® software.The integrated approach combined image data with field data, and various combinations of ancillary data, including: (1).DEM only, (2).DEM and NDVI, (3).DEM and slope, (4).DEM, slope and NDVI, and (5).DEM, slope and hydrology.Combining the datasets involved layer stacking, with the exception of the hydrology dataset.The hydrology dataset was used in the classification process by assigning vegetation communities 4 and 21 to the riparian areas from the hydrology map.A supervised Maximum Likelihood Classifier (MLC) was applied to the Landsat5 TM and SPOT5 for the pixelbased image-only and integrated approaches.MLC is one of the most widely used image classification algorithms [11,15,34,[45][46][47][48] and was therefore selected for the PBIA component of this study.
The classification results were smoothed using a majority filter (3 × 3 kernel).Two vegetation community maps were generated at 1:25,000 and 1:100,000 spatial scales, where polygons less than 0.04 and 0.25 ha respectively were eliminated.Since ENVI MLC classification auto generates its own class numbers, these classes had to be relabeled to reflect the 22 communities on the four resultant maps prior to accuracy assessment [49].

Object-Based Image Analysis
Four object-based approaches were applied to the Landsat5 TM and SPOT5 multi-spectral image datasets to map 22 vegetation communities using eCognition® software, including: (

Segmentation
Multi-resolution segmentations were applied to the Landsat5 TM and SPOT5 image datasets for the four object-based approaches.For the NN classifications we used one level of segmentation, called level one, segmented after incorporating a shapefile including the training areas (named test/training area (TTA) mask).For the step-wise classifications three hierarchical levels of segmentation were used to obtain contextual information and were called levels one, two and three.Thus, level one in the step-wise approaches was used for identifying landform patterns, level two was used to segment homogeneous features across the rises, low hills and hills, and the final classification of vegetation communities was carried out in level three aided by the information contained in levels one and two [50].

Nearest Neighbor Classifications
Training areas were selected as samples once the Landsat5 TM and SPOT5 image datasets were segmented, the.Samples were manually selected based on the discrete objects from the 241 training areas representing the 22 vegetation communities.The 22 vegetation communities were incorporated in the class hierarchy for both NN classifications (image-only and integrated).The object-features used for the image-only NN classification included the mean and standard deviation (SD) of all image layers.For the integrated NN classification, the feature space also included the mean and SD of all the image layers, plus the mean and SD of the NDVI, DEM and slope model.The GEOBIA NN classification algorithm was applied to level 1 and considered all vegetation communities.The minimum membership function was set to 0 [50].

Image-only Step-Wise Classification
The Landsat5 TM and SPOT5 images were segmented at three levels to capture the diversity of vegetation communities across the study area.At level one, the landform patterns were captured, including the riparian communities, using a range of object features, such as band ratios and layer values, and class-related features, including relations to neighboring objects.The landform patterns were incorporated into the class hierarchy as parent classes, including alluvial plains, plains and rises, low hills and hills, slope and riparian [50].
The second step captured the homogeneous vegetation communities which occurred on the plains and rises, and low hills and hills in level two.The objects classified at level two also included the riparian areas that had high vegetative reflectance and objects that were small (linear and compact).Finally, level three, the highest level of detail, was used to classify the vegetation communities that were small in extent and of heterogeneous nature using a variety of object features including layer values, geometry and position, and customized features, and class related features including relations to neighbor objects, relations to super objects and relations to classification [50].

Integrated Step-Wise Classification
The integrated step-wise classifications incorporated the ancillary datasets into the classification processes.Similar to the image-only step-wise classification, objects at level one captured the landform patterns including alluvial plains, plains, rises and low hills, hills, slope and riparian.The two riparian communities (4 and 21) were first classified based on the thematic attributes of the two hydrology shapefiles.The landform patterns were classified based on the mean values of the DEM.The slope model was used to classify the vegetation communities that occur in areas with high slope across the rises and low hills, and hills.At level two, vegetation communities occurring across the plains, rises and low hills, and hills were classified based on a range of object and class-related features, similar to the image-only step-wise classification [50].Vegetation communities across the alluvial plains were classified at level three.A final step in the classification process was the refining of other vegetation communities classified at previous levels, using contextual information.

Cost-Effectiveness
Two components were used for the cost-effectiveness: accuracy assessment which was based on overall accuracy and cost associated with each mapping approach.
Confusion matrices were calculated for the 28 resultant vegetation community maps using ENVI® software for the accuracy assessment component.The training areas reserved for the accuracy assessment were used to match with the classes in the classified images.Overall accuracy was computed and results analyzed for each approach, image dataset and spatial scale map.
The cost incorporated four main components: field data acquisition and preparation; image data acquisition and preparation; image classification; and accuracy assessment.The cost-effectiveness evaluation assumed hardware and software requirements were already in place.The cost and time associated with each component and subcomponent were documented in preparation for the costeffectiveness evaluation.

Vegetation Communities
Twenty two vegetation communities were described for the study area (Appendix 1), based on the field data sampling, multi-variate analysis and vegetation classification.Some vegetation communities were more extensive than others.The most common and widespread vegetation community was Eucalyptus tectifica low woodland (1).This community occurred across a range of landform patterns and substrates, predominantly plains, rises and low hills in association with community 7. Corymbia dichromophloia low open woodland (2) was widespread on a range of sandstone landform patterns, mainly low hills and hills.Acacia spp., Grevillea spp., Gardenia spp.tall sparse shrubland (22) was also extensive and characteristic of broken sandstone plateaus and hills.Eucalyptus phoenicea low open woodland (10) was another widespread community, generally in association with community 22 and 2. Eucalyptus miniata mid open woodland (6) existed as three forms, an influence of substrate and landform.A typical form was on the plateaux, a second form occurred on rugged sandstone hill slopes and the third was on heavier soils adjacent to drainage lines on the alluvial plains.
The major river system (Bullo River) and its perennial tributaries were dominated by Melaleuca leucadendra mid woodland (21), also including significant paperbark swamps not associated with the riparian areas.The second riparian community (4) was dominated by Lophostemon grandiflorus mid woodland on ephemeral stream channels across plains, rises, low hills, hills and plateaux, usually in association with community 21.On the drainage depressions, Corymbia polycarpa mid open woodland (11) and Melaleuca viridiflora low woodland (20) were either discrete or intergraded.Adjacent to drainage depressions and major riparian systems were relict levees dominated by Corymbia bella mid woodland (3) and neighbouring on the plains was Corymbia foelscheana mid woodland (18).
The community dominating the slopes was Buchanania obovata and Terminalia latipes low open woodland (12).Also common on scarps and the heads of gullies hills was Xanthostemon paradoxus, Pouteria sericea, and Ziziphus quadrilocularis mid woodland (28), dry vine thicket.There were a number of less extensive communities, largely including on permanent springs (13) and swamps dominated by either tussock grasses (8) or sedges (30).Eucalyptus brevifolia low open woodland (15) and Melaleuca sericea low open woodland (16) were confined to hills in the south-west corner of the study area on a distinctive geology type [40].

Cost-Effectiveness
The results of the cost-effectiveness assessment are divided into: (1).accuracy assessment which was based on overall accuracy and defined as the 'benefit', and (2). the cost and time associated with the approaches and image datasets.

Accuracy Assessment
The number of sites selected for map attribution (API)/classification training (PBIA and GEOBIA) and accuracy assessment had the same number of sites and were equal in area (km 2 ) for the majority of vegetation communities.The number of sites per vegetation community for the accuracy assessment ranged from one to 27, with an average of six.The spatial distribution of field sites (Figure 1) shows that they were not randomly distributed throughout the study area, but the majority were confined to the main access track and non-gazetted tracks.However, to sample the required number of sites for quantitative accuracy assessment, field sampling costs would double, which was unaffordable.Nonetheless, the total cost of field sampling was AU$120,000, including staff salaries, travel allowance, and vehicle costs.Helicopter hire (wet rate) alone exceeded AU$28,000.
API presented the highest overall accuracy with 66.67% (Figure 2) [43].On average, PBIA presented the lowest overall accuracies which ranged from 30 to 53% [49].For both OBIA and PBIA, the dataset demonstrating the highest overall accuracies was the Landsat5 TM with DEM and slope integrated, and at the 1:25,000 spatial scale.For GEOBIA, the overall accuracy ranged from 32% for the NN image-only Landsat5 TM classification to 53% for the integrated step-wise classification applied to the SPOT5 image [50].
The results of the PBIA image-only and integrated approaches indicate that the incorporation of ancillary data considerably improved overall classification results.The incorporation of a DEM and slope model improved classification results for both image datasets.Overall, any combination of image and ancillary data improved overall accuracy when compared to the image-only approach applied to Landsat5 TM and SPOT5.Other studies report that the incorporation of ancillary data improves classification results for a range of land cover and vegetation community mapping applications, and in urban areas [19][20][21][22]51,52].The incorporation of the DEM and slope model improved the classification of the image-only approach by 10% [50].For both PBIA image-only and integrated approaches, Landsat5 TM gave higher overall accuracy results compared to SPOT5.This was attributed to the additional spectral information of Landsat5 TM, particularly its visible and mid-infrared bands, even though SPOT5 has a higher spatial resolution.This finding is consistent with other studies by Gibbes et al. [16] and Stickler et al. [18].Similarly, with the GEOBIA approaches, the incorporation of ancillary data considerably improved overall classification accuracy for both the step-wise and the NN classifications, and for the two image datasets.The dataset presenting the highest overall accuracy for the NN classifications was SPOT5 with ancillary datasets integrated; a 10% improvement to using the image data in isolation.When compared to Landsat5 TM, incorporating ancillary data increased the overall accuracy result by 3%.For all GEOBIA approaches, SPOT5 gave higher overall accuracy results, on average, compared to Landsat5 TM.This finding was attributed to the higher spatial resolution of SPOT5.The NN classifications applied to Landsat5 TM and SPOT5 failed to depict the floristic detail of vegetation communities.The image-only step-wise classification displayed an improvement to the NN classifications.This was improved by 7% when the ancillary data was included in the integrated step-wise approach.The image dataset and technique illustrating the highest overall accuracy was the SPOT5 integrated step-wise classification.
In summary, it appears that the spectral information of the Landsat sensor compensated for its coarser spatial resolution in the PBIA approaches, which is consistent with findings by Stickler et al. [18].Conversely, the higher spatial resolution of the SPOT5 compensated for the fewer spectral bands in the four GEOBIA approaches.The benefit of GEOBIA is the ability to use a combination of spectral and object-based information to build context into rulesets [53]; therefore, is not solely reliant upon spectral information alone, unlike PBIA.Gibbes et al. [16] also reports the benefits of using high resolution imagery, although the limitation is the cost and is also dependent on the research question.

Overall Accuracy (%) Mapping Approach, Image Dataset and Spatial Scale
In terms of spatial scale, a slight difference was apparent in overall accuracy for all PBIA and GEOBIA techniques and image datasets.Appendix 2 contains 25 thematic maps for the seven mapping approaches, three image datasets and at the two spatial scales.On average, 1:100,000 thematic maps produced slightly higher accuracies.This result can be attributed to the spatial and attribute detail eliminated from the maps.Both API maps generated the same accuracy as none of the eliminate polygons intersected any validation field site.If additional sites were sampled for the accuracy assessment, the result may be different, although not significantly, given that 51 small polygons (less than 4 ha) were eliminated from a total 700.There are no formally published studies that provide a comparison of cost effectiveness using API, PBIA and GEOBIA that capture the floristic and structural components of vegetation communities.The majority of studies that do exist focus on capturing land cover which comprise very broad classes, fewer classes (14 or less), and minimal structural information.These studies show higher overall accuracies for GEOBIA compared to PBIA in a range of environments, including native vegetation to agricultural landscapes.In a northern Australia's tropical savanna, results showed a higher overall accuracy of the GEOBIA over PBIA, 78.5% versus 69.5% respectively [54].The results by Whiteside et al. [54] suggest that GEOBIA has potential as an alternative method (over per-pixel) approaches for extracting land cover information from medium to high resolution satellite imagery.This is consistent with findings within the literature [55][56][57].However, in agricultural areas there was no significant difference in overall accuracy between GEOBIA and PBIA [58].

Cost
Four main cost components were documented and evaluated, including: field data acquisition and preparation; image data acquisition and preparation; image classification; and accuracy assessment.Two of the components were consistent for each approach, including the field data acquisition and preparation and the accuracy assessment.These components do not influence the cost-effectiveness comparison for the seven tested approaches, the three image datasets, or the two spatial scales.However, they are a major component of any vegetation community mapping project, irrespective of the approach, image dataset, and to a degree, spatial scale of the output thematic map.The importance of the field data sampling, plant identification, data basing, multi-variate analysis and vegetation classification is often not reported on in terms of cost and time.Thus, it is presented in this study to highlight the fact that it is one of the most labor-intensive and expensive steps in vegetation community mapping.
The results of the consistent cost components are presented in Table 2.The consistent component displaying the highest cost was field data acquisition and preparation, equating to AU$324,710 (Australian dollars), while the accuracy assessment came to AU$1,890 (100% salaries).Field data acquisition and preparation is made up of three subcomponents including field sampling, plant identification and data basing, and multi-variate analysis and vegetation classification.
Of the three consistent subcomponents, plant identification and data basing was the most expensive and equated to AU$201,700 (100% salaries).This included sorting specimens to family and generic levels, identifying all specimens collected that were unable to be identified in the field, and data basing the plot sheets.The cost of plant identification would fluctuate based on botanical expertise, time of sampling (wet season versus dry season), and full floristic or dominant species inventories.A recent study indicates a full floristic inventory is not necessary for vegetation community mapping [40], thus this cost would be significantly reduced.This cost was followed by the field sampling, which totaled AU$120,000, and included staff salaries (including on-costs and travel allowance) and vehicle costs (including helicopter hire, vehicle lease, maintenance and fuel).Helicopter hire (wet rate) alone was AU$28,700.Increased field sampling is thought to increase overall accuracy, however, few studies have been undertaken to address this.Benson [59] found costs vary depending on the complexity of the vegetation and the aims of the project.For instance, local spatial scales (1:25,000) far exceeded the cost of regional surveys (1:100,000).The less labor-intensive and expensive consistent subcomponents were the multi-variate analysis and vegetation classification, which came to AU$3,010 (100% salaries).The remaining two components including image data acquisition and preparation, and image classification are relevant to the cost-effectiveness comparison of the seven mapping approaches.Figure 3 illustrates the costs associated with the image data acquisition and preparation of the three image datasets.Pre-processing of the 1:50,000 aerial photography equated to AU$611, SPOT5 AU$216 and Landsat5 TM AU$36.The cost of image acquisition was the highest for the aerial photography (AU$24,000), followed by SPOT5 (AU$10,000) and Landsat5 TM (AU$760).
Figure 4 demonstrates the cost associated with the image classification component.The API is the most labor-intensive, thus the most expensive approach equating to AU$11,200 (100% salaries).The least expensive approach was the PBIA image-only approach, with AU$140 (100% salaries).The GEOBIA approaches ranged between AU$420 and AU$3,325.The integrated step-wise classification was the most expensive GEOBIA approach, almost eight times more expensive than the NN classification.This large difference in cost was due to the difference in hours of work required to set up the step-wise ruleset.Although setting up a step-wise rule set from scratch is very demanding in terms of hours, once this has been set up it can be easily adapted to adjacent areas.In that case, the total amount of time invested would be considerably reduced.

Segmentation Classification
Figure 5 summarizes the results obtained in the cost-effectiveness analysis for the seven approaches based on the four cost components.Figure 6 illustrates the highest overall accuracy results for the seven approaches.Figures 5 and 6 are based on the following image datasets and spatial scales for the seven mapping approaches: (1) API: 1:50,000 aerial photography mosaic at 1:100,000, (2) PBIA image-only: Landsat5 TM at 1:100,000, (3) PBIA integrated: Landsat5 TM (DEM and slope) at 1:100,000, (4) GEOBIA NN image-only: SPOT5 at 1:100,000, (5) GEOBIA NN integrated: SPOT5 (DEM and slope) at 1:100,000, (6) GEOBIA step-wise image-only: SPOT5 (contextual information) at 1:100,000, and (7) GEOBIA step-wise integrated: SPOT5 (DEM, slope and contextual information) at 1:100,000.Taking into account the four cost components, API is the most expensive and labor-intensive mapping approach and also has the highest overall accuracy.The total cost was AU$388,180.PBIA was the cheapest mapping approach (AU$338,000); however had the lowest average overall accuracy for the image-only classification (37%).The cost of the GEOBIA approaches ranged from AU$338,000 (NN classifications) to AU$341,000 for the step-wise integrated classification, which also displayed the highest overall accuracy of the GEOBIA approaches.There is a correlation with cost and accuracy in that the lower the cost of classification approach, the lower the overall accuracy.
Training of the analyst was required for the PBIA classification and involved seven hours.More training was required for the GEOBIA segmentation and classification.The total number of hours spent on analyst training was 61 and included a five-day training course, online tutorials and one-on-one training with experienced staff.No training was required for the analyst to undertake the API mapping approach.However, training and experience would be required in such a study if an analyst had no experience with API.Training should be factored into future mapping programs if staff are inexperienced with any of the components.There have been very few cost-effectiveness studies undertaken on vegetation community mapping programs around the globe.There are, however, a number of examples from Australia, where Neldner [60] undertook a study highlighting the value and utility of vegetation survey and mapping programs in Queensland, and suggested costs can be calculated based on time, staff, transport, data and laboratory analysis, and publication costs.The field work component was found to account for 60% to 70% of total costs.This result is consistent with the findings of this study where the field survey component was the most expensive and labor-intensive.Other land resource mapping programs across northern Australia have been calculated for a number of surveys, yet nothing has been formally documented.
In south eastern Australia, cost estimates have been tested in pilot vegetation surveys across various 1:100,000 map sheets in the Murray-Darling Basin [59].Costing was based on API methods and included costs of field survey, botanical support, vehicle hire, aerial photograph digitizing and attributing and the cost of each photo.The study found costs vary depending on the complexity of the vegetation and the aims of the project.For instance, local spatial scales (1:25,000) far exceeded the cost of regional surveys (1:100,000).Studies have acknowledged that the benefit to cost ratio of vegetation surveys is approximately 40:1 [59].The benefits of vegetation surveys and maps are numerous largely providing a basis for land use planning, environmental impact assessment, community participation, research and teaching, conservation, fauna distribution, forestry and military uses, amongst others [60].
The selection of approach to map vegetation communities is dependent on budget, resources, staff expertise, imagery and ancillary dataset availability, time constraints (deadlines) and the required degree of accuracy.These constraints are also influenced by the perceived application of an output vegetation community map.Thus, to generate the most accurate map, based on this study, aerial photography interpretation is the most appropriate method, although sufficient funds and adequate timeframe is required.For projects with limited timeframes, low budgets and a perceived accuracy of less than 55%, PBIA applied to Landsat5 TM would be acceptable.If, however, experienced staff were available, GEOBIA presents an alternative to aerial photography interpretation, using either high or moderate spatial resolution imagery.The results of this study suggest that GEOBIA is a potential technique to capture the floristic component and associated structural elements of vegetation in tropical savanna vegetation communities using the step-wise approach.

Conclusions
This study compares the cost-effectiveness of mapping vegetation communities using: (1).seven mapping approaches, (2).three image datasets at varying degrees of spatial resolution, and (3).two spatial scales including 1:25,000-property management scale and 1:100,000-regional scale.The seven approaches included: (1).aerial photography interpretation, (2).pixel-based image-only classification, (3).pixel-based integrated classification, (4).object-based image-only NN classification, (5).object-based integrated NN classification, (6).object-based image-only step-wise classification, and (7).object-based integrated step-wise classification.The approaches were applied to 1:50,000 aerial photography (approach 1), and SPOT5 and Landsat5 TM multi-spectral data (approaches 2 to 7).Overall accuracy ranged from 28 to 67% for the seven approaches.Aerial photography presented the highest overall accuracy (67%) for the 1:25,000 spatial scale map, and the image-only pixel-based approach applied to Landsat5 TM at 1:25,000 demonstrated the lowest (28%).Of the four main cost components, the most labor-intensive was the field data acquisition and preparation, followed by image data acquisition and preparation.The cost of image acquisition was the highest for the aerial photography, followed by SPOT5 and Landsat5 TM.The accuracy assessment component was the least expensive.Overall, the API was the most labor-intensive, thus most expensive approach.The least expensive approach was the PBIA image-only classification.Of the four GEOBIA approaches, the integrated step-wise classification was the most expensive.When the four components are factored into the total cost for the seven mapping approaches, API remains the most expensive and labor intensive approach, and also has the highest accuracy.PBIA was the cheapest mapping approach, however, had the lowest overall accuracy.Of the four GEOBIA approaches, the NN image-only gave the lowest overall accuracy and was the cheapest approach.The GEOBIA step-wise integrated classification displayed the highest overall accuracy and highest cost.There was a definite correlation with overall accuracy and associated costs for the seven mapping approaches, where an increase in overall accuracy reflected an increase in total costs.Suggestions to improve this work include the use of a fuzzy accuracy assessment to compensate for continuous variables like vegetation communities.The application of fuzzy rules to this study would strengthen the results of vegetation communities with poor separability and similar landscape position.WorldView-2 imagery, comparable in the number of spectral bands to Landsat5 TM and a much higher spatial resolution to SPOT5 (8-band multispectral and a spatial resolution of 46 cm), would be of interest in this study region and expected to yield better cost-effective results.An extension of the cost-effectiveness evaluation could potentially include set-up costs such as hardware and software.
1) an image-only NN classification (image data and training areas), (2) an integrated NN classification (image and ancillary data and training areas) (3) an image-only step-wise classification (image data), (4) an integrated step-wise classification (image and ancillary data).

Figure 2 .
Figure 2.Overall accuracy for the seven vegetation community mapping approaches (aerial photography interpretation (API); pixel-based image analysis (PBIA) (integrated and image-only); geographic object-based image analysis (GEOBIA) (integrated and image only; nearest neighbor and step-wise)); applied to the three image datasets (aerial photography, Landsat5 TM and SPOT5), for two spatial scales (1:25,000 and 1:100,000).

Figure 3 .
Figure 3. Costs in Australian dollars (AU$) associated with the image data acquisition and preparation component, including its two sub-components (image acquisition and image pre-processing) for the three image datasets (aerial photography, SPOT5 and Landsat5 TM).

Figure 5 .
Figure 5. Distribution of the total costs, in Australian dollars (AU$), for the seven mapping approaches.Cost for each component (field data acquisition and preparation; image data acquisition and preparation; image classification, and accuracy assessment) is presented in a different color.Note that, in order to represent different orders of magnitude, the scale used for each cost component is different (color and scale used for each cost component is specified in the legend).

Figure 6 .
Figure 6.Highest overall accuracy results for the seven mapping approaches tested.

Table 1 .
Table 1 outlines how the cost and time were estimated for each component/subcomponent and mapping approach.Components and subcomponents used for estimating the cost and time associated with each mapping approach.
* denotes consistent component for the seven mapping approaches.

Table 2 .
Consistent cost components for vegetation community mapping subcomponents, detailed cost and total cost in Australian dollars (AU$).