Close Home
iForest - Biogeosciences and Forestry
vol. 7, pp. 85-91
Copyright © 2014 by the Italian Society of Silviculture and Forest Ecology
doi: 10.3832/ifor1035-006

Technical Advances

Detecting mistletoe infestation on Silver fir using hyperspectral images

Mario Ančić (1)Corresponding author, Renata Pernar (1), Milan Bajić (2), Ante Seletković (1), Jelena Kolić (1)


The narrow ecological niche and the extensive exposure to pollution makes silver fir (Abies alba Mill.) a highly endangered species. It is currently the most damaged tree species in Croatian forests ([30], [29]). In addition to insects and fungi, fir is increasingly attacked by mistletoe (Viscum album L. ssp. abietis [Weisb.] Abrom), which physiologically weakens the tree by its remarkable consumption of water ([13], [21], [11]) and mineral nutrients ([41]), thus predisposing the host to pest and disease attacks. Mistletoe is also considered an indicator of existing stress ([2], [33], [22], [40]). It has been reported that an excessive number of shrubs may even cause tree death ([28], [5], [42], [6], [36], [24]). Morevover, mistletoe-infected trees have lower increment than uninfected trees ([19], [37], [38], [23], [4]) and should therefore have a cutting priority as long as their timber still retains its commercial value.

A high correlation between fir infestation and mistletoe distribution has been reported ([15], [18], [4]) with an almost exponential trend in the mistletoe population increase ([23]). An increase in the intensity of mistletoe infestation was also noted, favoring its spread to healthy trees ([16]). Tree health status can be determined in two ways: by ground and by remote sensing observations, of which color infrared (CIR) images are particularly important ([25], [26]). In Croatia, forest damage over large areas has been assessed by means of CIR aerial imagery since 1988, though mistletoe could not be detected based on these images. Preliminary pilot tests aimed at assessing fir damage using multispectral digital images have been carried out in 2004 ([27]), confirming that mistletoe could be successfully detected using a high-resolution multispectral digital camera. In contrast, standard field methods for mistletoe detection are generally based on visual survey (using binoculars) of edge trees along roads or in clearings, from which only the upper crown is visible ([17]), and therefore do not often reflect the real conditions in the field. For example, Diminić et al. ([10]) suggests that ground analyses of mistletoe presence in the crowns of standing trees are not fully reliable.

The limitations of standard field methods call for new approaches to assess tree health and quantify mistletoe infestation. The importance and efficiency of high-resolution multispectral and hyperspectral imagery for the goals mentioned above have been highlighted by Entcheva-Campbell et al. ([12]), Somers et al. ([34]) and Cho et al. ([7]). By recording a high number of narrow bands and increasing the spectral resolution, hyperspectral images make it possible to discern very narrow spectral features ([1], [32]). Hyperspectral research in Croatia was launched within the Croatian mine action center (CROMAC) programme focused on the reduction of mine-polluted areas in Croatia ([8], [3]). A line hyperspectral scanner provided the first reflectance curves of fir and mistletoe, which were used as entry parameters for further research ([26], [27]).

The aim of this research was to develop an efficient and reliable method of mistletoe detection over large areas using a hyperspectral scanner and supervised classification methods.

Materials and methods 

Hyperspectral scanning was performed in two areas characterized by the presence of beech-fir forests in the Dinaric region of Croatia (Lipovlje: 44° 50’ N, 15° 11’ E; Ramljani: 44° 44’ N, 15° 23’ E). At the time of aerial survey (April), deciduous species had not yet burst into leaves. Data collected with a hyperspectral line scanner ImSpector V9 (spectral range 430-900 nm, spectral resolution 7 nm) mounted on an helicopter were recorded as individual lines, together with spatial data from the Position and Orientation System (POS). These data were used to conduct parametric geocoding and obtain a raw geocoded hyperspectral cube (Fig. 1). After formulating the real geocoded hyperspectral cube, an image similar to a color infrared was obtained by selecting the following bands: 94 (895 nm, the near infrared part of the spectrum); 54 (695 nm, the red part of the visible spectrum); 22 (535 nm, the green part of the visible spectrum). Hyperspectral scanning provided pure spectral samples (endmembers), and showed a difference between the spectral curves of fir and mistletoe. The spectral samples formed a part of the spectral booklet created by means of the TNTmips Free® Hyperspectral module. The booklet served as a basis for interpreting the data collected by aerial surveys over large areas. Samples according to three classes (vegetation, shadow, and rock) were taken from the whole area. These samples differ not only visually but also in terms of spectral curves. Based on the collected samples, the following supervised classifications were conducted: Feature Mapping, Hyperspectral analysis (Spectral Angle Mapper - SAM, Cross-correlation, Linear Unmixing, Fuzzy C means), as well as an unsupervised Fuzzy C means classification. Satisfactory results were obtained using SAM and Cross-correlation, of which the SAM method was the best (Fig. 2). The SAM algorithm is an easy and rapid method for mapping the similarity of image spectra to reference spectra. It is also a very powerful classification method because it holds back the influence of shading effects and enhances the target reflectance characteristics ([9], [14]). The reference spectra can either be taken from laboratory or field measurements or extracted directly from the image. SAM measures the spectral similarity by calculating the angle between the two spectra treating them as vectors in a n-dimensional space ([20], [39], [31], [14]). Small angles between the two spectra indicate high similarity, while large angles indicate low similarity. A supervised classification procedure was conducted for the selected areas in the vegetation class in order to delineate the spectral curves for fir and mistletoe. Pixels differing from their neighbors (surrounding pixels) based on the shape of their spectral curve were selected in the supervised classification process (Fig. 3). The spectral curve of the selected pixels followed the shape of the reference spectral curve for mistletoe obtained by scanning under controlled conditions (Fig. 4). The procedures of selecting the class of mistletoe and testing within the same class included the selection of the samples for SAM classification separately with 4°, 5° and 6° of threshold value (degrees were selected by testing). The arithmetic means were calculated separately for each threshold. As the degrees within the SAM classification increased, the number of the selected samples also increased, but a shift also occurred in the spectral curves of the selected mistletoe samples with 4°, 5° and 6° thresholds. The largest difference occurred in the near infrared part of the spectrum, which showed a downward tendency towards the fir spectrum, whereas the spectral curves in the visible and “Red Edge” part of the spectrum retained their shape and exhibited only very slight fluctuations. The selected curves of the arithmetic means for 4°, 5° and 6° SAM classifications were compared each other and with the spectral curves obtained for fir (Fig. 5). After the mistletoe samples were selected and SAM classifications performed, a SAM classification within the vegetation class was carried out using pure spectral samples (endmembers) for fir and mistletoe as reference. Comparison for all estimations was performed using a repeated measures analysis of variance (ANOVA) with α = 0.05 using the package SUAUJCUJDA 7.1® ([35]).

Fig. 1 - 3D image of the hyperspectral cube in TNTmips Free, wavelength 0.610 μm.
Fig. 2 - Result of the application of the SAM classification method with 7° threshold value for (a) vegetation, (b) rock and (c) shadows.
Fig. 3 - Comparison of mistletoe and fir spectra in the magnified part.
Fig. 4 - Mean values of spectral reflectance for silver fir and mistletoe over the wavelength range considered in this study.
Fig. 5 - Comparison of arithmetic mean values of reflectance for mistletoe (using SAM classification with threshold 4°, 5° and 6°) with mean values obtained for silver fir (green line).

The SAM classification procedure with threshold value 5°, which proved to be optimal, was performed over all the image segments. In order to test for any major changes in solar irradiation during hyperspectral scanning data collection, mean values of insolation were calculated based on data recorded by FODIS. Spectral ranges where large portions of electromagnetic radiation are absorbed by the atmosphere (atmospheric maxima) were excluded from data processing to avoid spurious results in the spectrum values.

Moreover, in order to verify the impact of the relief on vegetation classification results, a digital elevation model was constructed. Three directions (north-south - NS, the flight direction - DF, and orthogonal to the flight direction - ODF) were considered (Fig. 6), each including two sub-categories (before hill and after hill). Results were verified by aerial recording with an unmanned aircraft (hexacopter), for which all flight parameters were recorded (position, height, direction, speed, etc). Survey was performed at several flight altitudes above the crowns (~30 m, ~50 m, ~70 m).

Fig. 6 - Flight directions for which vegetation spectrum values were measured.

Results and discussion 

Of the applied methods of supervised (Cross correlation, SAM, Fuzzy C means, Feature Mapping) and unsupervised (Auto-classify) classifications, satisfactory results were obtained using both Cross correlation and SAM methods, though the latter proved to be the best. Based on the mistletoe samples identified using SAM classification, it was possible to detect and select mistletoe bushes within the fir crown. However, a high variability was observed in the data collected, as reflected in the misclassification of a part of the overlaying selected pixels as fir classes obtained by automatic classification methods. This could be due to the partially overlapping spectral curves obtained for fir, not clearly differentiated as the spectral curve for mistletoe.

Using the SAM method with 4°, 5° and 6° threshold values in the selected samples within the mistletoe class, mean values of spectral curves for the mistletoe were obtained, clearly distinct from the fir spectrum for all the tested degrees. By increasing the degree, a partial overlap of the near infrared part of the spectrum was observed in the reflectance curves for mistletoe and silver fir. This implies that for threshold degrees higher than 6° fir would be discerned instead of mistletoe (Fig. 5).

The final outcome of the investigation carried out was the identification of the spatial distribution of mistletoe-infected silver fir trees over the whole studied area (Fig. 7). The results of testing the effects of insolation on classification outcomes are presented graphically in Fig. 8, after subtraction of maximum absorbance values. As expected, insolation varied during scanning, though no abrupt leaps or large fluctuations were observed around the central values. Therefore, we concluded that insolation did not affect the classification obtained.

Fig. 7 - Results of the SAM classification using endmembers (pure spectral samples) as reference for mistletoe (white pixels).
Fig. 8 - Graphic representation of measured insolation (solar irradiation) values and its arithmetic mean (red line) across the analyzed spectral range (430-809 nm).

The impact of the relief on the classification obtained was assessed by comparing reflectance values measured in the different directions (NS, DF, ODF - see Fig. 6). The largest deviation was found for the NS direction in the visible part of the spectrum, up to the “Red Edge” area (430-700 nm), and again from the center of the near infrared area to the end of the spectrum (835-900 nm). In the visible part of the spectrum, differentiation in the blue (455-475 nm), cyan (476-495 nm) and green (496-535 nm) are particularly conspicuous (Fig. 9a). Comparison of values for the three measuring directions showed that NS- and DF-curves did fairly overlap, while the ODF spectral curve diverged in the visible part of the spectrum over the range 455-650 nm (Fig. 9b).

Fig. 9 - Comparison among arithmetic means and their mean values in NS direction (a) and comparison of arithmetic means by measuring directions (b).

Finally, data collected from the flight with the unmanned aerial vehicle at ~30 m fully supported the spatial distribution of mistletoe-infected silver fir trees drawn by the SAM classification with 5° threshold value, therefore confirming that a clear identification of mistletoe bushes in fir crowns using hyperspectral images is precisely feasible (Fig. 10). Accordingly, we recommend the use of hyperspectral images and the SAM classification (threshold: 5°) as a reliable method for mistletoe detection on fir.

Fig. 10 - Aerial photograph of the recording site taken at an altitude of approx. 30 m. Red circles indicate the position of mistletoe bushes on the silver fir crowns.

Future research are needed to corroborate the approach used here and to extend its validity to other site/ecosystems. However, the method applied has several advantages over the standard field methods, that imply felling a large numbers of trees and expensive field work ([10]). Hyperspectral image analysis allows a rapid and economic detection of mistletoe over large areas. Moreover, early detection of infestation hotbeds is fundamental in both commercial forest harvesting (reducing losses in roundwood value caused by mistletoe infestation) and forest protection.


An assessment of mistletoe infection on silver fir was conducted in two beech-fir forests in the Dinaric region of Croatia by the use of hyperspectral images and supervised classification methods. The following conclusions were made:

  1. Mistletoe attacks can be successfully detected by hyperspectral images. The best classification method for detection was Spectral Angle Mapper (SAM) classification with 5° threshold value.
  2. Limitations of standard field methods currently in use call for new approaches based on remote-sensed hyperspectral images to assess fir damage and detect mistletoe.
  3. The methods used here may allow a cost-effective mapping of the distribution and intensity of mistletoe infection in silver fir, and can be considered an helpful tool for forest management and planning, forest protection and harvesting.


Aspinall RJ, Marcus WA, Boardman JW (2002). Considerations in collecting, processing and analysing high spatial resolution hyperspectral data for environmental investigations. Journal of Geographical Systems 4: 15-29.
::CrossRef::Google Scholar::
Atzberger C, Werner W (1998). Needle reflectance of healthy and diseased Spruce stands. In: Proceedings of the “1 EARSeL workshop on imaging spectroscopy”. Remote Sensing Laboratories, University of Zurich, Switzerland, pp. 271-283.
::Google Scholar::
Bajić M, Gold H, Božičković K (2008). Application of the airborne remote sensing methods in the humanitarian countermine action (Croatian orig.). In: “Marked the 10 Anniversary of HCR-CTRO Scientific Council” (Knapp V ed). Croatian Mine Action Centre, Centre for Testing, Development and Training, Zagreb, Croatia, pp. 41-88.
::Google Scholar::
Barbu C (2009). Impact of mistletoe attack (Viscum album ssp. abietis) on radial growth of silver fir. A case study in the north of eastern Carpathians. Annals of Forest Research 52: 89-96.
::Online::Google Scholar::
Cervera JM, Villaescusa R (1977). Inventario de los abetares del Valle de Arán afectados por el muérdago. Boletín de la Estación Central de Ecología 6: 41-48.
::Google Scholar::
Chéret V, Dagnac J, Fromard F (1987). Le dépérissement du sapin dans les Pyrénées luchonnaises. Revue Forestière Française 39: 12-24.
::CrossRef::Google Scholar::
Cho MA, Debba P, Mutanga O, Dudeni-Tlhone N, Magadla T, Khuluse SA (2012). Potential utility of the spectral red-edge region of Sumbandila sat imagery for assessing indigenous forest structure and health. International Journal of Applied Earth Observation and Geoinformation 16: 85-93.
::CrossRef::Google Scholar::
CROMAC ARC team (2002). Trial II Evaluation Report. Deliverable D15, Report Version: 2.0.0, Airborne Minefeld Area Reduction - ARC, European Commission project 2001-2003, IST-2000-25300, European Commission, Brussels, Belgium, pp. 83.
::Google Scholar::
De Carvalho OA, Meneses PR (2000). Spectral correlation mapper (SCM); an improvement on the spectral angle mapper (SAM). In: Summaries of the “9 JPL Airborne Earth Science Workshop”. JPL Publication 00-18, Pasadena, CA, USA, pp. 9.
::Google Scholar::
Diminić D, Potočić N, Jazbec A, Zupanić M (2011). Infestation of common mistletoe and nutrition status of silver fir in Gorski Kotar (Croatia). Croatian Journal of Forest Engineering 32 (1): 223-237.
::Google Scholar::
Ehleringer JR, Cook S, Tieszen LL (1986). Comparative water use and nitrogen relationships in the mistletoe and its hosts. Oecologia 68: 279-284.
::CrossRef::Google Scholar::
Entcheva-Campbell PK, Rock BN, Martin ME, Neefus CD, Irons JR, Middletin EM, Albrechtova J (2004). Detection of initial damage in Norway spruce canopies using hyperspectral airborne data. International Journal of Remote Sensing 24: 5557-5583.
::CrossRef::Google Scholar::
Fisher JT (1983). Water relations of mistletoe and their hosts. In: “The biology of mistletoe” (Calder M, Bernhardt P eds). Academic Press, Sidney, Australia, pp. 161-184.
::Google Scholar::
Girouard G, Bannari A, El Harti A, Desrochers A (2004). Validated spectral angle mapper algorithm for geological mapping: comparative study between Quickbird and Landsat-TM. In: Proceedings of the “XX ISPRS Conference” (Altan O ed). Istanbul (Turkey) 12-23 July 2004, International Society for Photogrammetry and Remote Sensing, pp. 599-604. [online]
::Online::Google Scholar::
Hofstetter M (1988). Über die Verbreitung der Mistel in der Schweiz. Schweiz Z Forstwes 139 (2): 97-127.
::Google Scholar::
Idžojtić M, Pernar R, Kauzlarić Z, Abramović M, Janković D, Pleše M (2003). Intensity of infestation of silver fir (Abies alba Mill.) with mistletoe (Viscum album L. ssp. abietis (Wiesb.) Abrom.) on the area of the forest administration Delnice (Croatian orig.). Forestry Journal 127 (11-12): 545-559.
::Google Scholar::
Idžojtić M, Glavaš M, Zebec M, Pernar M, Dasović M, Pavlus M (2005). Infestation of silver fir (Abies alba Mill.) with mistletoe (Viscum album L. ssp. abietis (Wiesb.) Abrom.) in Croatia (Croatian orig.). Forestry Journal 129 (11-12): 559-573.
::Google Scholar::
Idžojtić M, Pernar R, Glavaš M, Zebec M, Diminić D (2008). The incidence of mistletoe ( L. ssp. (Wiesb.) Abrom.) on silver fir ( Mill.) in Croatia. Biologia 63 (1): 81-85.
::CrossRef::Google Scholar::
Klepac D (1955). Influence du gui sur l’accroissement des sapinières (Croatian orig.). Forestry Journal 79 (11-12): 231-244.
::Google Scholar::
Kruse FA, Boardman JW, Lefkoff AB, Heidebrecht KB, Shapiro AT, Barloon PJ, Goetz AFH (1993). The spectral image processing system (SIPS) - Interactive visualization and analysis of imaging spectrometer data. Remote Sensing of Environment 44: 145-163.
::CrossRef::Google Scholar::
Lamont B (1983). Mineral nutrition of the mistletoes. In: “The Biology of Mistletoe” (Calder M, Bernhardt P eds). Academic Press, Sidney, Australia, pp. 185-204.
::Google Scholar::
Moorthy I, Miller J, Noland TL (2008). Estimating chlorophyll concentration in conifer needles with hyperspectral data: an assessment at the needle and canopy level. Remote Sensing of Environment 112: 2824-2838.
::CrossRef::Google Scholar::
Noetzli KP, Müller B, Sieber TN (2003). Impact of population dynamics of white mistletoe ( ssp. ) on European silver fir (). Annals of Forest Science 60(8): 773-779.
::CrossRef::Google Scholar::
Olivia J, Colinas C (2007). Decline of silver fir ( Mill.) stands in Spanish Pyrenees: role of management, hystoric dynamics and pathogens. Forest Ecology and Management 252: 84-97.
::CrossRef::Google Scholar::
Pernar R (1997). Application of results of aerial photograph interpretation and geographical information system for planning in forestry. Annales experimentatis silvarum culturae provehendis, Universitas studiorum Zagrebiensis 34: 141-189.
::Google Scholar::
Pernar R, Ančić M, Seletković A (2007a). Application of colour infrared aerial photographs for the assessment of forest damage in the Gospić forest administration (Croatian orig.). Forestry Journal 131 (11-12): 507-522.
::Google Scholar::
Pernar R, Bajić M, Ančić M, Seletković A, Idžojtić M (2007b). Detection of misteletoe in digital colour infrared images of infested fir trees. Periodicum biologorum 109 (1): 67-75.
::Google Scholar::
Plagnat F, Brossier J (1969). Silver fir stands with mistletoe. Silviculture of infested fir stands. Tentative conclusions. Thoughts on fir mistletoe. Revue Forestiere Française 26: 553-561.
::CrossRef::Google Scholar::
Potočić N, Seletković I, Ugarković D, Jazbec A, Mikac S (2008). The influence of climate properties on crown condition of Common beech ( L.) and Silver fir ( Mill.) on Velebit. Periodicum Biologorum 110 (2): 145-150.
::Google Scholar::
Prpić B, Seletković Z (2001). Ecological constitution of the Silver fir (Croatian orig.). In: “Silver fir ( Mill.) in Croatia” (Vukelić J, Gračan J, Dundović J eds). Academy of Forestry Sciences, Zagreb, Croatia, pp. 255-276.
::Google Scholar::
Rowan LC, Mars JC (2003). Lithologic mapping in the Mountain Pass, California area using Advanced Spaceborn Thermal Emission and Reflection Radiometer (ASTER) data. Remote Sensing of Environment 84: 250-266.
::CrossRef::Google Scholar::
Shaw G, Manolakis D (2002). Signal processing for hyperspectral image exploitation. IEEE Signal Processing Magazine 19: 12-16.
::CrossRef::Google Scholar::
Smith KL, Steven MD, Colls JJ (2004). Use of hyperspectral derivative ratios in the red-edge region to identify plant stress responses to gas leaks. Remote Sensing of Environment 92: 207-217.
::CrossRef::Google Scholar::
Somers B, Verbesselt J, Ampea E M, Sims N, Verstraetena W W, Coppina P (2010). Spectral mixture analysis to monitor defoliation in mixed-aged Labill plantations in southern Australia using Landsat 5-TM and EO-1 Hyperion data. International Journal of Applied Earth Observation and Geoinformation 12: 270-277.
::CrossRef::Google Scholar::
StatSoft (2006). Statistica (v. 7.1). Data analysis software system, Tulsa, OK, USA.
::Online::Google Scholar::
Tsopelas P, Angelopoulos A, Economou A, Souilioti N (2004). Mistletoe () in the fir forest of Mount Parnis, Greece. Forest Ecology and Management 2002: 59-65.
::CrossRef::Google Scholar::
Uščuplić M (1992). Influence of forest management on the occurrence of silver fir mistletoe ( L.). University of Belgrade, Bulletin of the Faculty of forestry 74: 7-18.
::Google Scholar::
Uščuplić M, Dautbašić M (1998). Diseases and pests that threats forest ecosystems in Bosnia and Herzegovina (Bosnian orig.). Works of the Faculty of Forestry University of Sarajevo 1: 19-26.
::Google Scholar::
Van der Meer F, Vasquez-Torres M, Van Dijk PM (1997). Spectral characterization of ophiolite lithologies in the troodos ophiolite complex of Cyprus and its potential in prospecting for massive sulphide deposits. International Journal of Remote Sensing 18 (6): 1245-1257.
::CrossRef::Google Scholar::
Wang Q, Li P (2012). Hyperspectral indices for estimating leaf biochemical properties in temperate deciduous forests: Comparison of simulated and measured reflectance data sets. Ecological Indicators 14: 56-65.
::CrossRef::Google Scholar::
Watson DM (2001). Mistletoe - a keystone resource in forests and woodlands worldwide. Annual Review of Ecology and Systematics 32: 219-249.
::CrossRef::Google Scholar::
Zuber R (1983). Forstschutz. In: “Wald- und Forstwirtschaft in Graubünden”. Ründerwald, Beiheft 12, pp. 115-119.
::Google Scholar::


Paper Contents

Paper Sections

Paper Figures



Ančić M, Pernar R, Bajić M, Seletković A, Kolić J (2014).
Detecting mistletoe infestation on Silver fir using hyperspectral images
iForest - Biogeosciences and Forestry 7: 85-91. - doi: 10.3832/ifor1035-006
First Previous Next Last
© iForest

Download Reference

Paper ID# ifor1035-006
Title Detecting mistletoe infestation on Silver fir using hyperspectral images
Authors Ančić M, Pernar R, Bajić M, Seletković A, Kolić J
Close Download