Remote Sensing Remote Distinction of a Noxious Weed (musk Thistle: Carduus Nutans) Using Airborne Hyperspectral Imagery and the Support Vector Machine Classifier

Remote detection of non-native invasive plant species using geospatial imagery may significantly improve monitoring, planning and management practices by eliminating shortfalls, such as observer bias and accessibility involved in ground-based surveys. The use of remote sensing for accurate mapping invasion extent and pattern offers several advantages, including repeatability, large area coverage, complete instead of sub-sampled assessments and greater cost-effectiveness over ground-based methods. It is critical for locating, early mapping and controlling small infestations before they reach economically prohibitive or ecologically significant levels over larger land areas. This study was designed to explore the ability of hyperspectral imagery for mapping infestation of musk thistle (Carduus nutans) on a native grassland during the preflowering stage in mid-April and during the peak flowering stage in mid-June using the support vector machine classifier and to assess and compare the resulting mapping accuracy for these two 613 distinctive phenological stages. Accuracy assessment revealed that the overall accuracies were 79% and 91% for the classified images at preflowering and peak flowering stages, respectively. These results indicate that repeated detection of the infestation extent, as well as infestation severity or intensity, of this noxious weed in a spatial and temporal context is possible using hyperspectral remote sensing imagery.


Introduction
Non-native invasive plant species are well known for their successful exploitation of environmental resources (e.g., water, space, light and nutrients), largely due to their aggressive and competitive ability, prolific seed production and seed longevity [1][2][3].Along with having a competitive advantage in resource utilization over native species, many non-native weeds grow in the absence of natural enemies [1,2].This allows non-native invasive weeds to easily establish and compete with native plants in the environment [2].Due to these characteristics, non-native invasive weeds usually move from small, manageable infestations to larger areas, reaching levels where control is economically prohibitive [2].Some non-native invasive species can dominate the vegetative canopy and eventually can form monotypic stands [4,5].Invasion may alter nutrient accumulation and cycling, hydrology, carbon sequestration and grass production [2,4,6,7].Weed infestation has been recognized as a primary reason for loss in global biodiversity and even species extinction [8][9][10][11][12].Therefore, non-native invasive plant species have long been of interest to natural resource managers, ecologists and biological conservationists [13][14][15][16].
One non-native weed of major ecological and economic importance in the US, Canada, Australia, New Zealand and southern Africa is musk thistle or nodding thistle (Carduus nutans L.), an aggressive member of the sunflower (Asteraceae) family that is native to North Africa and Eurasia [17][18][19][20][21][22].Musk thistle infests a considerable amount of land area in the US (found in 45 states) and Canada (found at least in five provinces) and has a long history in these countries as an invasive species [20,21].Musk thistle has been declared a noxious weed in many US states, Canadian provinces and Australia.It infests forest, pasture and rangelands, roadsides, railroad rights-of-way, ditch banks, stream beds, waste or disturbed areas, agricultural fields [17,18,23], as well as residential areas.The presence of this weed can impede grazing by livestock and decreases the forage production due to competition with grasses [17,18].Since livestock do not graze musk thistle, it easily outcompetes grassland forage species and can dominate entire fields with up to 150,000 plants per hectare [1,22].It has been determined that a single musk thistle per 1.5-m 2 reduces pasture yields by 23% [24].
Early determination of the distribution and severity of rapidly spreading invasive populations of musk thistle is needed to implement mitigation treatments, but accurate assessments are often difficult or impossible to obtain with ground surveys, because of the extensive land area involved, time and labor required and inaccessibility of many areas [25][26][27].Therefore, remote sensing has received considerable attention as a rapid, inexpensive and non-destructive method for assessing non-native noxious species invasions [6,27].It is a technology that provides a complete assessment instead of an assumed level based on traditional field-based sub-sampling methods.
A wide range of sensor systems, including aerial photographs, airborne and satellite multispectral and hyperspectral images have been successfully utilized for mapping the distribution of certain species [28][29][30][31][32][33][34][35][36].A large number of investigators [37][38][39][40][41][42] have discussed the superiority of hyperspectral imaging sensors (also known imaging spectrometry), such as Airborne Imaging Spectroradiometer for Applications (AISA; Specim Inc., Oulu, Finland) over multispectral instruments, such as the Landsat Thematic Mapper.The former systems collect spectral information in continuous, narrow spectral channels, while the latter ones record spectral reflectance at a few wide, non-continuous wavelengths separated by spectral segments where no measurements are taken.The hyperspectral system was designed to separate the surface optical properties into tens of bands with spectral resolution less than 20 nm or smaller through the visible, near infrared (NIR) and mid-infrared regions of the electromagnetic spectrum (400-2,500 nm).These narrow and continuous bands provide ample spectral information to identify and distinguish between spectrally similar, but unique, surface materials [39,43,44].Consequently, hyperspectral imagery provides the potential for more detailed and accurate spectral information than multispectral images.Hyperspectral data have been applied to a wide range of spatial modeling and vegetation mapping, including detecting crop stress and classifying non-native invasive plants at the species level across a range of community and ecosystem types [41].Hyperspectral sensors can capture reflectance features of chemical compounds in plants, including lignin, cellulose, nitrogen, chlorophyll, carotenoids and water [6,[45][46][47].Another important characteristic of remotely sensed imagery for mapping plant species is the spatial resolution or pixel size representing the smallest area identifiable in the imagery.Higher spatial resolution imagery often contains detailed intraspecies spectral variability than lower resolution imagery, when target species occur in small patches or cover smaller areas than the pixel size [28,48,49].Timing, location, species composition, spectral and spatial resolutions are all important for vegetation mapping at the species level, because spectral reflectance from vegetation is determined by structural and biochemical properties, including leaf area, leaf angle distribution, architecture, biochemical and pigment (lignin, cellulose, chlorophylls, carotenes, nitrogen and water) contents [46,47] over time and space.
The literature indicates that there is a growing body of work on the hyperspectral mapping of non-native invasive species [50][51][52][53][54][55][56][57].Most of the previous hyperspectral remote sensing studies dealing with mapping invasive plants have commonly utilized Airborne Visible Infrared Imaging Spectrometer, HyMap and Compact Airborne Spectrographic Imager images with spectral resolution ≥ 10 nm and spatial resolution ≥ 1.5 m.In addition to this, spectral angle mapper and mixture tuned matched filtering were commonly employed as the classification algorithms after dimension reduction or transformation of hyperspectral data, such as principal component analysis, minimum noise fraction or vegetation indices.Some of the non-native invasive species frequently distinguished and mapped in the USA include saltcedar (Tamarix spp.) [33,52,58,59], leafy spurge (Euphorbia esula L.) [51,56,57], spotted knapweed (Centaurea maculosa Sam.) [32,48,53] and yellow starthistle (Centaurea solstitialis) [50,54,55].Recently, a mixture of Canada thistle (Cirsium arvense L. Scop.) and musk thistle infestation in a Nebraska floodplain was mapped using AISA hyperspectral imagery [33].
Despite the use of hyperspectral data for identification and mapping of non-native invasive plants and the large distribution of musk thistle in the US and other countries, the application of very narrow banded (very high spectral resolution) and high spatial resolution hyperspectral imagery for mapping musk thistle on grasslands and rangelands is not well studied.Utilizing very narrow banded and high spatial resolution hyperspectral data to discriminate and separate invasive musk thistle from highly variable rangeland vegetation comprised of many species of grasses and broadleaved forbs is needed in order to improve monitoring population dynamics and developing long-term adaptive management practices for musk thistle.Typically, in the southern Great Plains (in Texas and Oklahoma, USA), the diversity of rangeland species at any one site includes species with cool-season (C 3 ) or warm-season (C 4 ) photosynthetic pathways that usually have very different periods of physiological activity during the growing season and would have markedly different spectral properties at any one point in time.
Our primary objective was to evaluate the use of very narrow banded (average about 2.25 nm in the visible, 3.5 nm in the NIR regions and 2.9 nm in the visible and NIR wavelengths; Table 1) and 1-m spatial resolution hyperspectral data and the support vector machine (SVM) classification method for detecting and discriminating musk thistle from the co-occurring vegetation in a rangeland pasture.The secondary objectives of our study were to (1) examine the spectral properties of musk thistle and surrounding vegetative land cover classes for species discrimination, (2) explore the ability of aerial hyperspectral imagery for detecting and mapping musk thistle infestations at the preflowering stage in mid-April and the peak flowering stage in mid-June and (3) compare the mapping accuracy for these two distinctive phenological growth patterns.Our study presents a new application of hyperspectral data in terms of relative value of seasonal images for management applications, spatial and spectral resolutions of hyperspectral data and classification method presented.We initiated this study to fill an important gap in the application of hyperspectral imagery for detecting an invasive plant at the species level.Since the hyperspectral remote detection of musk thistle has not been well documented for preflowering and full flowering phenological stages, our study presents the first application of hyperspectral data for discriminating musk thistle at two distinct growth stages.

Site Description
A musk thistle infestation was identified on a grassland site near the city of Friona in Parmer County located about 114 km south west of Amarillo, Texas, USA (34°38'2.85"N,102°48'4.18"Wand elevation approximately 1,211 m) in early June 2003 (Figure 1) and monitored thereafter until mid-summer 2006.Musk thistle was sprayed in 2004 and mowed by the landowner in 2005 (personal observation); however, the infestation rebounded in 2006.Image collection flights were conducted over the study area when the musk thistle was at preflowering stage on 14 April and at peak flowering stage on 15 June 2006.Herbaceous plants found at the site consisted of C 3 (cool-season) grasses-needle and threads (Stipa spp.) and bromegrasses (Bromus spp.)-C 4 (warm-season) grasses-Johnsongrass (Sorghum halepense L.) and purple three-awns (Aristida spp.)-and broad-leaf forbs-Russian thistle (Salsola iberica L.) and tansy mustard (Descurainia spp.).In addition, a field planted to wheat (Triticum aestivum L.) occupied a small portion of the image in April.Russian thistle, Johnsongrass and brome are non-native invasive species, while tansy mustard is native, but can be invasive.Needle and thread and purple threeawn are native; needle and threads are considered desirable rangeland forage grasses, whereas purple threeawn is not and is an indicator of overgrazing.None of these species were planted by humans, with the exception of wheat.Actively growing and fully green plants in April were tansy mustard, brome, wheat and musk thistle.In June, actively growing species were musk thistle, Johnsongrass and Russian thistle, because cool-season tansy mustard, brome grass and wheat had senesced by mid-June.

Imagery Acquisition
An Airborne Imaging Spectrometer for Applications mounted in a Cessna 172 airplane was used to acquire the images.The AISA is a hyperspectral pushbroom type airborne imaging spectrometer.It had a spectral range from 505 to 900 nm with bands as small as 0.7 nm and as large as 9.4 nm.Under typical conditions in a Cessna 172 airplane, the AISA can acquire between 20 and 60 bands at a spatial resolution between 1-m and 4-m per pixel.It captured 384 pixels across track with a 22.7° field of view.Resolution along the track for a single image was limited by a disk space of two gigabytes.Downwelling irradiance was measured simultaneously via a fiber optic sensor mounted on the roof of the airplane.The downwelling irradiance was used to calibrate the image.Navigation data for georectification were collected by an onboard Boeing C-MIGITS II integrated Inertial Navigation System/Global Positioning System (INS/GPS; Systron Donner Inertial Division, Concord, CA, USA).Radiometric and geometric corrections were done using Caligeo image preprocessing software version 1.7 (Specim Inc., Oulu, Finland) using sensor specific calibration data.The images were atmospherically calibrated using Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) in Environment for Visualizing Images software (ENVI; Exelis Visual Information Solutions, Boulder, CO, USA).Using this atmospheric calibration method, radiance-at-sensor was converted to apparent surface reflectance.Spatial resolution of the image collected over the research site was 1-m by 1-m with a flight height of about 1,000-m, and there were 50 bands ranging from 509 nm to 885 nm (Table 1).

Imagery Classification
The 1-m image allowed clear and visual identification of musk thistle and other vegetation based on the spectral contrast among vegetation and senescent grasses or bare ground to extract training samples for classification.Therefore, training samples manually extracted from the imagery, consisted of 5-10 polygons, each of which had 5-20 pixels from the canopy or each land cover type with 75% canopy cover at locations easily identified on the ground and on the image (Table 2).Training samples were also used to obtain the spectral profile of plant species and bare ground (Figure 2).Before performing classification, the spectral profile of image components for plant species was examined (Figure 2).The spectral profile of each species was the mean values of training samples.Spectral reflectance for plant species in the visible and NIR regions were tested for statistical significance using paired t-tests assuming unequal variance at α = 0.05 for all data pairs using Microsoft Excel.The largest significant difference appeared in the NIR; therefore, the last 15 bands (bands 36-50; 764-885 nm) in this region formed a spectral subset for image classification.
Image classification was performed using the SVM in ENVI.The SVM is a supervised machine learning method that performs classification based on the statistical learning theory.The SVM classifies data by separating a hyperplane that provides the best separation between classes in a multidimensional feature space.This hyperplane is the decision surface on which the optimal class separation takes place.The optimal hyperplane is the one that maximizes the distance between the hyperplane and the nearest positive and negative training example called the margin.From a given set of training samples, the optimization problem is solved to find the hyperplane that leads to a sparse solution.Although the SVM is a binary classifier in its simplest form, implementation of the SVM classifier in ENVI was extended to more than two classes by splitting the problem into a series of binary class separations (ENVI User's Guide).
In order to represent more complex shapes than linear hyperplanes, a variety of kernels, including the polynomial, the radial basis function and the sigmoid, can be used for performing SVM classification in ENVI.The SVM was employed using the radial basis function kernel for performing the pairwise classification.There is also a penalty parameter that can be introduced to the SVM classifier to allow for misclassification during the training process.The penalty parameter was set to its maximum value, whereas a classification probability threshold of zero was used in order to classify all pixels (ENVI User's Guide).
The image taken in April was classified for bare ground, senescent grass, musk thistle, tansy mustard, brome grass and wheat.Since our target species was musk thistle, bare ground and senescent grass classes were combined into a group and reported as non-vegetation, while tansy mustard, brome grass and wheat classes were combined into another group and reported as other vegetation using post classification procedure in ENVI for April image classification.Combining classes resulted in a three-category (musk thistle, other vegetation and non-vegetation) final map of the study site.The same classification and procedure for combining classes were applied to the June image, except for the other vegetation category that included Johnsongrass and Russian thistle (Figure 3).

Accuracy Assessment
Classification of remotely sensed imagery is the process of assigning variables into discrete categories of useful information.However, errors of classification occur because of mixed pixels, similar spectra of different materials, positional errors and human errors [60].The most effective way to represent classification accuracy is via an error matrix [61][62][63][64].Therefore, accuracy assessment for classification was made by constructing an error matrix for each classified image.An error matrix is a square, with the number of columns and rows being equal to the numbers of categories whose classification accuracy is being evaluated [61][62][63][64].Error matrices for each classification map were generated by comparing the classified classes with the ground verification data.Error matrices to evaluate the classification accuracy were calculated, including overall, user's and producer's accuracies.
There is no single set standard for selection of the image and ground areas for comparison [61][62][63]65].Because a pixel in an image has an arbitrary location on the ground and because positional errors of maps and global positioning system receivers become significant with small pixel sizes, areas based on geographic information system polygons are used frequently [63].However, using individual pixels is appropriate if a per-pixel classification is assessed for accuracy, which avoids problems caused by generating -homogeneous‖ polygons on a landscape [66,67].It has been shown that pixel positional error results in a conservative bias of the accuracy assessment [68]; therefore, the unavoidable positional error introduced into this assessment would result in lower or conservative estimates of mapping accuracy [60,69,70].

Field Validation
Field validation (accuracy assessment) of the classification performance was carried out using verification data (ground control points) at the site.Verification data containing 330 locations (points) with an equal spacing were created using the Farm Works software package (CTN Data Service, Hamilton, IN, USA).The verification points were exported to a field computer equipped with Farm Works software and a global positioning system receiver and navigated on the ground on the same dates when images were acquired at the site.At each verification point within a 1-m radius area, presence (at least one musk thistle plant found) and absence (no musk thistle plant found) of musk thistle were recorded.This ground data collection prior to image processing created an unbiased field validation.Of the 330 points, only 208 and 270 points fell within April and June images, respectively.Therefore, only those verification points were used for accuracy assessment.Subsequent to image classification, those points were overlaid on a land cover map, where one-to-one matching was performed to contract an error matrix for each date.Training samples were not included in the validation dataset.

Results
The mean reflectance spectra of training samples extracted from the hyperspectral imagery acquired for musk thistle and associated vegetation in April and June are shown in Figure 2. Paired t-test comparisons of mean reflectance spectra of training samples for the same land cover types are presented in Table 3.In April, musk thistle reflected significantly less light in the green region (510-592 nm) than tansy mustard, brome grass, senescent grass and bare ground (Figure 2(a)).In the same region, musk thistle had significantly more reflectance than wheat.Brome grass had significantly higher reflectance than musk thistle, tansy mustard and wheat and significantly less reflectance than senescent grass and bare ground in this wavelength region.The reflectance spectra of the training samples in the red region of the spectrum exhibited similar reflectance patterns to the green portion of the spectrum.In the 600-691 nm wavelength range, all compared pairs were significantly different, with one exception: that the comparison of musk thistle and wheat was insignificant (Table 3).Around 675 nm, reflectance from senescent grass, brome grass, tansy mustard, musk thistle and wheat started to increase.In the 701-758 nm wavelength range, all pairs of training samples for land cover classes were insignificant, with the exceptions of tansy mustard vs. senescent grass, brome grass vs. senescent grass and senescent grass vs. bare ground that were significant.In the NIR region (763-885 nm), musk thistle had significantly less reflectance than both wheat and brome grass, while it had significantly more reflectance than senescent grass, bare ground and tansy mustard (Figure 2(a) and Table 3).There was no significant difference between brome grass and wheat in the 763-885 nm wavelengths, whereas all other pairs were significant.Reflectance from bare ground gradually increased in the visible and NIR regions examined.
In June, musk thistle reflected significantly more light than Russian thistle and Johnsongrass, whereas it had significantly less reflectance than both senescent grass and bare ground in the 510-592 nm regions (Figure 2(b) and Table 3).There was no difference between Johnsongrass and Russian thistle, while all other pairs were significant in the green portion of the spectrum.In the red region (600-691 nm), reflectance from musk thistle was insignificant when compared to those of Johnsongrass and Russian thistle.In this part of the spectrum, the remaining pairs showed significant difference.Around 690 nm, reflectance from senescent grass, musk thistle, Russian thistle and Johnsongrass started to increase.From 701 nm to 758 nm, reflectance between musk thistle and senescent grass, Johnsongrass and senescent grass and senescent grass and bare ground was significant, while all remaining pairs were insignificant.Musk thistle had significantly lower reflectance than Johnsongrass and Russian thistle, whereas it had significantly higher reflectance than bare ground and senescent grass in the 763-885 nm range.Reflectance from bare ground steadily increased in the full spectrum studied.In both dates, reflectance of senescent grass was significantly lower than bare ground in the visible and NIR spectrums.The classification map in June for musk thistle, other vegetation and non-vegetation is illustrated in Figure 3.The error matrices that compared ground verification data with classified data in April and in June 2006 are contained in Table 4.An overall mapping accuracy of 79% in April and 91% in June was achieved for musk thistle detection.About 28% of the other land cover types were included into musk thistle category, whereas about 21% of the musk thistle were misclassified as other classes by the classification method performed in April.In June, approximately 12% of the musk thistle were excluded from the musk thistle class, while about 10% of other cover types were included in the musk thistle category.

Discussion
All physiologically active vegetation had lower reflectance in the visible region until around 706-730 nm and higher reflectance in the NIR beyond 730 nm when compared to senescent grass and bare ground.The visible region is defined as pigment-related absorption wavelengths that are primarily governed by the presence of chlorophylls a and b and carotenoids [28,29,71].This indicates that actively growing musk thistle, tansy mustard, brome grass, wheat, Johnsongrass and Russian thistle had higher chlorophylls and carotenoids than senescent grass.Musk thistle reflected more light than wheat and less light than tansy mustard and brome grass in the visible spectrum in April.This indicates that musk thistle had a higher amount of photosynthetic pigments than both tansy mustard and brome grass, but it had a lower amount of pigments than wheat during this early spring date, when wheat in this region was physiologically very active.
In June, musk thistle had a lower amount of photosynthetic pigments than Johnsongrass and Russian thistle.The NIR reflectance from vegetation increases with the increasing green leaf or canopy density and with decreasing gaps between them due to the internal scattering at the cell wall-air interfaces within the leaf [72].In both April and June, musk thistle had lower reflectance than wheat, brome grass, Johnsongrass and Russian thistle and higher reflectance than tansy mustard and senescent grass in the 763-885 nm wavelengths.This implies that musk thistle had lower leaf or canopy density when compared to wheat, brome grass, Johnsongrass and Russian thistle and higher leaf or canopy density than senescent grass and tansy mustard.
A similar observation to our study was made that flowering shoots of leafy spurge had higher reflectance in the visible region of the spectrum compared to co-occurring vegetation [60].However, these authors also found higher reflectance from leafy spurge in the NIR spectrum than that from other vegetation.In our study, musk thistle had lower reflectance, starting approximately around 720-740 nm, and stayed lower in the remaining NIR portion of the spectrum examined.Higher NIR reflectance from flowering leafy spurge may be an effect of sample collection from areas with high density of leafy spurge [60].Spectral reflectance of musk thistle in our study also closely agreed with the results previously reported for flowering yellow starthistle [30,55] and perennial pepperweed (Lepidium latifolium) [73].The peak flowering stage of these species was the best time to spectrally identify and differentiate them from background vegetation.Thus, acquiring very narrow banded hyperspectral imagery during the peak flowering of musk thistle could allow in detail identification of this noxious weed from the background vegetation.
There is no set standard for classification accuracy as to what level of accuracy is adequate.The acceptable level of mapping accuracy depends on the project-based choice that differs based on the needs of the end users and utilization of the classification results [60,62].An accuracy target of 85% was recommended by [74], and an overall accuracy target of 85% with no individual class accuracy < 70% was suggested by [75].Although there is no information related to the acceptable classification accuracies specified for remote detection of musk thistle, a threshold of the user's accuracy > 70% is required for mapping leafy spurge [51].Based on these recommended accuracy guidelines, our classification of the April image did not meet the overall map accuracy target of 85%, but provided the recommended levels of the user's and producer's accuracies.Results from the June classification indicated that the overall, producer's and user's accuracies were above recommended guidelines.We achieved these high accuracies because we used a suitable classification algorithm, a very narrow banded hyperspectral imagery with 1-m spatial resolution, a spectral range where separability largely occurred and the ground verification method extended to the entire study area for musk thistle infestation.One of the most likely utilizations of these accuracy maps would be to track changes in musk thistle population density over time and space.With consistent and repeatable classification methodology, the level of accuracies we accomplished in this study would present very accurate musk thistle distribution over the study site.
Image classification methods for weed identification and mapping have a tendency to overestimate the target species [48].In our study, the producer's accuracy of 88% was slightly lesser than the user's accuracy of 90% for the June image classification.Producer's accuracy is more important than user's accuracy for weed detection, because undetected or omission of target weeds (false negative) is the factor that most likely hampers control efforts [73].There is greater potential for misclassification of musk thistle when it grows in close proximity to dense patches of background live vegetation, in locations where musk thistle with canopy coverage is < 15% or when individual musk thistle plants are small in size.All of these factors affected classification accuracy in the present study.Our classification results were consistent with the findings of other studies.For example, hoary cress (Cardaria draba) with 30% cover [76], leafy spurge with cover as low as 10% [51] and spotted knapweed (Centaurea maculosa) with only 1% cover [53] were successfully detected using hyperspectral images with pixel sizes of 3 by 3-m, 3.5 by 3.5-m and 5 by 5-m, respectively.
Accuracy assessment performed for April and June image classifications resulted in a higher mapping accuracy in June than April.This can be attributed to one or all of the following reasons: The C 3 grasses and broad leaf species completed their life cycle by June, which resulted in fewer species present for classification.Phenological differences between preflowering and peak flowering may result in spectral differences.The peak flowering may cause changes in internal leaf structure and morphology, resulting in different NIR reflectance.Differences in water content of the leaves between April and June may also affect reflectance spectra.The collection of training samples for the background vegetation may also affect classification results between April and June.Regardless of the reasons, our classification results were consistent with the literature.For example, classification accuracies were better for yellow hawkweed and oxeye daisy (Chrysanthhemum leucanthemum) in full bloom than at early and post bloom stages [77].In addition, high mapping accuracies were achieved for the peak flowering stage of leafy spurge [60], which indicated that image classification at preflowering and post-flowering stages would have missed some areas occupied by this species because of its similar spectral characteristic with background vegetation.
With the history of remote sensing, researchers have sought the optimal strategy to extract as much information as possible from remotely sensed data by using the best possible physiological, compositional and structural distinctions in target species, classification algorithm, data transformation and band selection.The accuracy of a classification task can be improved by focusing on most informative bands or spectral regions.Just a few species-specific spectral bands may be sufficient to detect target species.Since both image classifications in April and June provided sufficient accuracies to be used for management decisions regarding musk thistle treatment, our study provides an application of band selection from a very narrow banded hyperspectral imagery and classification method for discriminating and mapping musk thistle.
It is more beneficial for the resource manager to have information regarding weed infestations as early as possible.The image classification in April was not as accurate as the June classification, but was of sufficient accuracy to make management decisions for treating this weed.Data from this study indicate that classification and mapping of musk thistle on an early spring image (April) before it reaches peak flowering, seed production and dispersal stage is sufficient to implement treatments when control measures are needed for the current growing season.When control measures are not needed or planned for the current growing season, mapping of musk thistle infestation at the flowering stage is the better option for more accurate mapping and planning for future control measures.

Conclusion
Remote sensing techniques offer important opportunities for detecting plant species and mapping the extent of invasions in a timely and spatial manner, especially when the target species have unique characteristics from the co-occurring vegetation mosaic.This study investigated the use of 1-m spatial resolution hyperspectral imagery and the support vector machine classifier for detecting musk thistle infestation in a rangeland pasture.Given the economic and ecological outcomes of invasion by this weed, our results indicate that hyperspectral imagery may be a useful tool for mapping musk thistle infestation.The accuracy assessment revealed that the overall mapping success were 79% and 91% for the classified images at preflowering and peak flowering stages, with the user's and producer's accuracy > 71% and > 78%, respectively.This indicates that detection of musk thistle during the flowering period is higher than the preflowering stage.Our results have implications for rangeland management and productivity.Maps generated at a 1-m scale provide the spatial extent and dynamics of musk thistle infestation over time.These maps can be used for monitoring, planning and control measures.Such maps have also been shown to provide animal distribution, grazing effects and the effects of natural and anthropogenic disturbances in rangelands.Using maps at a 1-m scale, the reduction in range productivity can be estimated, because a single musk thistle per 1.5-m 2 decreases pasture yields by 23% [24].The use of a 1-m spatial resolution hyperspectral image provided valuable estimates of musk thistle infestation.Therefore, we recommended that this methodology and technology should be considered when high scale maps are needed for rangeland research and management.Currently, the cost of hyperspectral imagery over large areas may be economically prohibitive in a commercial operation on native rangeland, but these costs will likely decrease with improvements in technology and may allow hyperspectral detection of weed invasion to become an economically useful tool for rangeland weed management.

Figure 1 .
Figure 1.Location of study site in Parmer County (a) in TX (b), USA.

Figure 3 .
Figure 3. Color infrared hyperspectral imagery acquired in June 15, 2006 (a), classified image (b) and views of musk thistle infestation (c,d) at the study site.The entire area is native grassland, with the exception of a wheat field on the far right, right of the main road.

Table 2 .
Summary of training samples used to extract spectral signatures and classify land cover classes in images taken over the study site near the city of Friona, TX, USA in April and June 2006.

Table 4 .
Confusion matrices (top for April and bottom for June image classifications) generated from the ground verification and classified data for musk thistle, other vegetation and non-vegetation components.