Next Article in Journal
Erratum: Kromer, R.A.; et al. A 4D Filtering and Calibration Technique for Small-Scale Point Cloud Change Detection with a Terrestrial Laser Scanner. Remote Sensing 2015, 7, 13029–13052
Next Article in Special Issue
Quantitative Estimation of Carbonate Rock Fraction in Karst Regions Using Field Spectra in 2.0–2.5 μm
Previous Article in Journal
Accuracy and Precision of Habitat Structural Complexity Metrics Derived from Underwater Photogrammetry
Previous Article in Special Issue
Mapping the Spectral Soil Quality Index (SSQI) Using Airborne Imaging Spectroscopy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of Epidermal Structures on Leaf Spectral Signatures of Ice Plants (Aizoaceae)

by
René Hans-Jürgen Heim
1,*,
Norbert Jürgens
1,
André Große-Stoltenberg
2 and
Jens Oldeland
1
1
Biodiversity Ecology and Evolution of Plants, Biocenter Klein Flottbek and Botanical Garden, University of Hamburg, Ohnhorststraße 18, 22609 Hamburg, Germany
2
Institute of Landscape Ecology, Heisenbergstraße 2, 48149 Münster, Germany
*
Author to whom correspondence should be addressed.
Remote Sens. 2015, 7(12), 16901-16914; https://doi.org/10.3390/rs71215862
Submission received: 1 October 2015 / Revised: 26 November 2015 / Accepted: 6 December 2015 / Published: 15 December 2015
(This article belongs to the Special Issue Field Spectroscopy and Radiometry)

Abstract

:
Epidermal structures (ES) of leaves are known to affect the functional properties and spectral responses. Spectral studies focused mostly on the effect of hairs or wax layers only. We studied a wider range of different ES and their impact on spectral properties. Additionally, we identified spectral regions that allow distinguishing different ES. We used a field spectrometer to measure ex situ leaf spectral responses from 350 nm–2500 nm. A spectral library for 25 species of the succulent family Aizoaceae was assembled. Five functional types were defined based on ES: flat epidermal cell surface, convex to papillary epidermal cell surface, bladder cells, hairs and wax cover. We tested the separability of ES using partial least squares discriminant analysis (PLS-DA) based on the spectral data. Subsequently, variable importance (VIP) was calculated to identify spectral regions relevant for discriminating our functional types (classes). Classification performance was high, with a kappa value of 0.9 indicating well-separable spectral classes. VIP calculations identified six spectral regions of increased importance for the classification. We confirmed and extended previous findings regarding the visible-near-infrared spectral region. Our experiments also confirmed that epidermal leaf traits can be classified due to clearly distinguishable spectral signatures across species and genera within the Aizoaceae.

Graphical Abstract

1. Introduction

The leaf epidermal anatomy, together with its biochemical and biophysical properties, strongly controls the functional properties and the adaptive fitness of plants. Leaf optical properties depend on the leaves’ biochemical and biophysical properties. Optically-measured spectral signatures of leaves are commonly used for analyzing these properties. For example, spectral signatures were used to measure biochemical components, such as pigment concentration [1,2] nutrient content [3] or other biochemicals, such as tannins [4,5]. Cellular leaf structure [6] or leaf water content [7] are typical biophysical properties of leaves that can be captured by spectral signatures. Although epidermal structures (ES) are known to affect the spectral signatures of leaves [8,9,10], they were mostly neglected in recent studies of leaf optical properties.
Already in the 1970s and 1980s, spectral properties of leaf ES, in particular hairs and glaucous waxes of desert plant species, were extensively studied [11,12,13]. The ecological function of ES-like hairs relates to the regulation of transpiration and to increasing the reflectance by up to 70% of the incoming radiation between 400 nm and 700 nm; this results in a significant decrease of heat load compared to glabrous leaves under similar conditions [9]. While the effect of pubescence seems reasonably studied, little is known about the effects of other ES, such as bladder cells and papillary epidermal cells. The state of the art in the 1980s was limited to a spectral range of 1100 nm [11,14]. Today, spectral analysis can range from 450 nm–2500 nm. While wavelengths up to 700 nm are typically reflected by plant pigments, spectral regions beyond the NIR, i.e., the short wave infra-red (SWIR), contain information on biochemical compounds, such as lignin, cellulose and leaf water [15,16,17]. Yet, few studies of these spectral regions with regard to ES are available [18].
Epidermal structures of leaves can be highly diverse in arid ecosystems. In these habitats, the epidermis controls the interaction between the living plant and a harsh environment with extreme values of heat, radiation and very low air humidity [19]. In southwestern Africa, one of the few arid hot-spots of biodiversity on this planet [20], the ice plant family (Aizoaceae) developed a remarkable range of epidermal structures [21]. The Aizoaceae contain 2271 accepted species within 146 genera [22]. Members of the family are dominant in the vegetation in parts of South Africa and Namibia, with regard to biomass and the number of taxa [23]. Ihlenfeldt and Hartmann [21] and Jürgens [23] developed a typology, which was extended by Niesler [24] to finally include eight different types of ES occurring within the Aizoaceae (e.g., big hetero-cellular bladder cells, big papillate hetero-cellular bladder cells, homo-cellular types and hybrid forms). This diversity of ES within one plant family from the same bioclimatic region makes the Aizoaceae an interesting object to study the effects of ES on spectral signatures. A choice of species with a close evolutionary relationship helps to ensure that spectral responses of the sample differ not just due to the line of evolutionary development. An approach such as this has not been conducted with a phylogenetically close set of species thus far.
Many recent studies showed that it is possible to distinguish between plant species based on their spectral signatures [25,26,27]; even cultivars of tea [28], wine [29] or other commercially-relevant crops [30] were spectrally distinguished. However, in all of these cases, the species were a priori known and either physiologically or evolutionary different; these conditions strongly facilitate the separation of spectral responses. Ustin and Gamon [31] argued that it would make more sense to treat plant species or, more precisely, spectrally-measured vegetation canopies as “optical types”. This term was coined in relation to the idea of plant functional types [32] and recently revisited in very comprehensive reviews [33,34]. Epidermal structures could play an important role in defining Ustin and Gamon’s optical types.
In this study, we: (i) test whether Aizoaceae with similar epidermal types share similar spectral responses; (ii) apply a spectral classification in order to test whether different epidermal types are indeed spectrally distinguishable; and (iii) test which spectral regions, e.g., full, VIS, NIR, SWIR I and SWIR II, are most suitable for separating and classifying epidermal types (ETs) based on epidermal structures (ES). In a condensed expression, we want to test if ETs are similar in their spectral properties within an ET and if they significantly differ between other ETs.

2. Experimental Section

2.1. Sample Collection

The botanical garden, administrated by the University of Hamburg/Germany, hosts a large living collection of Aizoaceae with several hundred species from 82 genera. We used this living collection as a foundation to sample spectral signatures. Leaves of Aizoaceae species (Table 1) were surveyed ex situ with a 20-fold magnifier to define epidermal types (ET). These types are supposed to be confirmed due to their spectral signature by statistical classification. At the given magnification, we identified five ETs, bladder cells, hairs, flat epidermal cell surface (in the following, known as smooth), convex to papillary epidermal cell surface (in the following, known as rough) and waxes, and assembled them into groups (Figure 1). For each ET, we selected five species, across multiple genera, and ten individuals of each species to measure spectral signatures.
Table 1. Examined species, their assigned epidermal type and their subfamilies.
Table 1. Examined species, their assigned epidermal type and their subfamilies.
Epidermal Type (ET)SpeciesSubfamily
Bladder cellsDrosanthemum curtophyllumRuschioideae
L. Bolus NM
Drosanthemum archeriRuschioideae
L. Bolus NM
Mesembryanthemum crystallinumMesembryanthemoideae
Linnaeus
Delosperma vogtsiiRuschioideae
L. Bolus NM
Trichodiadema rupicolaRuschioideae
L. Bolus NM
Hair coverGibbaeum shandiiRuschioideae
N. E. Brown GC
Cheiridopsis purpureaRuschioideae
L. Bolus NM
Braunsia apiculataRuschioideae
(Kensit)L. Bolus
Juttadineteria albataRuschioideae
(L. Bolus) L. Bolus
Odontophorus marlothiiRuschioideae
N. E. Brown GC
RoughNeohenricia spiculataRuschioideae
S. A. Hammer CSJ
Hereroa puttkamerianaRuschioideae
(Dinter & A. Berger) Dinter & Schwantes
Nananthus aloidesRuschioideae
(Haworth) Schwantes GF
Rhombophyllum dolabriformeRuschioideae
(Linnaeus) Schwantes GF
Stomatium bolusiaeRuschioideae
Schwantes MDK
SmoothGlottiphyllum longumRuschioideae
(Haworth) N. E. Brown
Malotigena frantiskae-niederlovaeRuschioideae
Niederle
Aptenia haeckelianaMesembryanthemoideae
(A. Berger) Bittrich ex Gerbaulet
Schlechteranthus halliiRuschioideae
L. Bolus
Bergeranthus scapigerRuschioideae
(Haworth) Schwantes
Wax coverMalephora purpurea-croceaRuschioideae
(Haworth) Schwantes GF
Oscularia steenbergensisRuschioideae
(L. Bolus) H. E. K. Hartmann Bradleya
Scopelogena verruculataRuschioideae
(Linnaeus) L. Bolus
Amphibolia succulentRuschioideae
(L.Bolus) H. E. K. Hartmann Bradleya
Enarganthe octonariaRuschioideae
(L. Bolus) N. E. Brown GC
Figure 1. Pictures of the five epidermal types; each group is subdivided into two detailed photos and one photo of the plant’s habitus. (a,b,c) Braunsia apiculata, hairs; (d,e,f) Drosanthemum archeri, bladder cells; (g,h,i) Glottiphyllum longum; smooth; (j,k,l) Malephora purpurea-crocea, wax. The detailed photo in (l) is showing the blade of a scalpel used to abrade the wax layer, which Malephora purpurea-crocea deposits on its leaf surface. (m,n,o) Neohenricia spiculata, rough. All photographs were taken by Klaus von Bösekom with an AF-S VR Micro-Nikkor 105 mm 1:2, 8G IF-ED built by Nikon.
Figure 1. Pictures of the five epidermal types; each group is subdivided into two detailed photos and one photo of the plant’s habitus. (a,b,c) Braunsia apiculata, hairs; (d,e,f) Drosanthemum archeri, bladder cells; (g,h,i) Glottiphyllum longum; smooth; (j,k,l) Malephora purpurea-crocea, wax. The detailed photo in (l) is showing the blade of a scalpel used to abrade the wax layer, which Malephora purpurea-crocea deposits on its leaf surface. (m,n,o) Neohenricia spiculata, rough. All photographs were taken by Klaus von Bösekom with an AF-S VR Micro-Nikkor 105 mm 1:2, 8G IF-ED built by Nikon.
Remotesensing 07 15862 g001

2.2. Hyperspectral Data Collection

We recorded spectral signatures using an ASD plant probe attached to an ASD Field-Spec® 4 Hi-Res spectrometer. The spectrometer measures the spectral region from 350 nm–2500 nm with a spectral sampling interval of 1.4 nm in the range of 350–1000 nm and 1.1 nm until 2500 nm, resulting in a total number of 2151 bands. Before each measurement, the spectrometer was calibrated with a Zenith LiteTM Diffuse Reflectance Target, 95%R, from SphereOptics®. The calibration was done initially as soon as the spectrometer reached operating temperature and afterwards every 15 min. Flowering pots with the test object were positioned in a cardboard box (Figure 2) lined with black microcellular rubber. Initial experiments showed that this kind of rubber absorbs radiation very well in the range of the measured wavelengths (see Supplementary Materials Figure S1 for a sample spectrum). Lining the nearby environment with micro-cellular rubber (Figure 2a,d) helped to reduce diffuse light reflection and scattering, thus allowing for clean spectra with low noise.
Figure 2. Overview of the applied measuring arrangement in the greenhouse at the botanical garden of Hamburg. The height-adjustable cover (a) and the cardboard box are lined with black microcellular rubber to reduce the albedo during the measurement. A spectrometer (b) is placed on a trolley, and the flowering pots were positioned on the bottom of the cardboard box. A fiber optic cable is located at the front panel of the spectrometer and leads to the plant probe (c). The plant probe that was used during the experiment is placed between the cardboard box (d) and the spectrometer. The plant probe is illustrated above the cardboard box to emphasize that we measured all samples inside this box and applied the probe from above, in order to have direct physical contact with the sampled leaves. A notebook computer (e) was connected with the spectrometer to visualize and store the results of the measurement.
Figure 2. Overview of the applied measuring arrangement in the greenhouse at the botanical garden of Hamburg. The height-adjustable cover (a) and the cardboard box are lined with black microcellular rubber to reduce the albedo during the measurement. A spectrometer (b) is placed on a trolley, and the flowering pots were positioned on the bottom of the cardboard box. A fiber optic cable is located at the front panel of the spectrometer and leads to the plant probe (c). The plant probe that was used during the experiment is placed between the cardboard box (d) and the spectrometer. The plant probe is illustrated above the cardboard box to emphasize that we measured all samples inside this box and applied the probe from above, in order to have direct physical contact with the sampled leaves. A notebook computer (e) was connected with the spectrometer to visualize and store the results of the measurement.
Remotesensing 07 15862 g002
For data collection, either three different leaves (if the leaves were large enough to cover the light path entrance) or three randomly-selected piles of leaves were measured. For the leaf piles, multiple leaves were pressed against each other to obtain a dense area to cover the light path entrance completely. We always positioned the light path entrance in a way that the sampled leaves had physical contact with the elevated frame of the glass panel. We tried not to maintain the contact for too long because of the raised temperature. However, within one measurement, the temperature difference had no impact for us. Five internal repetitions (IR) per measurement were chosen for later averaging. In total, 750 spectra (10 individuals × 5 species × 3 measures × 5 IR = 750 spectra) per ET were measured.

2.3. Data Preprocessing

The raw data received were converted from radiance to reflectance data with the ViewSpecPro software provided by the manufacturer of the ASD. In order to identify outlier spectra, all spectra of each species were plotted, and obviously, deviating spectra were deleted manually. The spectral regions from 350 nm–400 nm and from 2400 nm–2500 nm were removed because they primarily consisted of stochastic noise. The dataset collected for all species was averaged per individual. The resulting spectral library contained 10 spectra per species and 50 spectra per ET, totaling up to 500 spectra. We deleted one outlier spectrum for the ET smooth and one for wax. Further preprocessing was applied to prepare the spectra for classification.

2.3.1. Jump Correction

We corrected two shifts in each spectral signature. They originated from the construction of the different detector fields of the sensor. All spectra exhibited such shifts at 1000 nm and 1800 nm. These irregularities, also called “jumps”, could be detected by a classification algorithm as a relevant feature, thus negatively influencing classification. Therefore, jumps were corrected using a multiplicative jump correction algorithm [35].

2.3.2. Standard Normal Variate Transformation

Standard normal variate (SNV) transformation is a typical preprocessing step in order to homogenize spectra. Basically, the mathematical transformation, SNV, is used to correct scatter effects in spectra [36,37]. First, each spectrum is centered on zero and afterwards scaled by the standard deviation ±1. However, as the aim of our study was to exploit heterogeneity between the groups, we wanted to test whether SNV transformation could improve our classification. Thus, we produced an additional SNV processed spectral library.

2.4. Classification

To confirm the visually-determined ETs, we applied a partial least squares discriminant analysis (PLS-DA) with the R package “DiscriMiner” [38]. PLS-DA is a commonly-applied algorithm for classifying spectral data [39]. In order to test if the classification quality varied for different regions, we tested subsets between 400 nm and 2400 nm. The PLS-DA was applied on a full spectrum (400 nm–2400 nm), on the VIS region (400 nm–760 nm), on the NIR region (761 nm–1100 nm), on the SWIR I region (1101 nm–1800 nm) and on the SWIR II region (1801 nm–2400 nm). These regions are typical subsets of the whole spectra [40]. To test if distinct bands had a greater impact on the classification than others, we calculated the variable importance in the projection (VIP) score [41,42]. The larger the VIP value, the more important the predictor is for classifying the response [39].

2.4.1. Partial Least Squares Discriminant Analysis

PLS was not originally designed as a tool for statistical discrimination. It is used routinely for regression and is usually applied when there are more predictor variables than observations [43,44]. This makes PLS highly suitable for high-dimensional data, e.g., chemometric or hyperspectral data. For each response variable yk, a regression model on the X-components is expressed as:
y k = h = 1 m ( X w h * ) c h + e = X W * c + e ,
w h * is a p dimension vector (p = number of parameters) containing the weights given to each original variable in the k-th component, and ch is the regression coefficient of yk on the h-th X-component variable [45]. PLS-DA is an extension of the PLS regression by a linear discriminant analysis [46] to be used for classification purposes.
The variable importance in the projection (VIP) scores obtained by the PLS are an important measure of each explanatory predictor [41]. The VIP score for the j-th variable can be calculated by:
V I P j = { p h = 1 m k R 2 ( y k , t h ) w h j 2 h = 1 m k R 2 ( y k , t h ) } 1 2
For each j-th predictor variable j = 1, p, where R 2 (a, b) stands for the squared correlation between items in vector a and vector b; t h = X h 1 w h , where X h 1 is the residual matrix in the regression of X on component   t 1 , …, t h -1; and w h is a vector of norm 1. Note that w h j 2 measures the contribution of each variable j to the h-th PLS component [45]. Thus, V I P j quantifies the influence on the response of each variable summed over all components and categorical responses relative to the total sum of squares of the model. As a rule of thumb, any independent variable with a VIP coefficient smaller than 0.8 had no impact on the classification [47].

2.4.2. Validation

To validate the classification results, we chose two strategies; firstly, a training, test dataset, and, secondly, a cross-validation (CV) approach. The spectral database was split 80 : 20 into a training dataset and a test dataset. The spectral databases contained 50 spectra for each ET, leaving 40 spectra for the training dataset and ten spectra for the test dataset. We deleted one averaged outlier spectrum for the ET smooth and one for wax, resulting in only 49 spectra for each of the two groups. A leave-one-out cross-validation (LOO-CV) was applied on the PLS-DA. LOO-CV is a special case of the k-fold CV where k = N, while N is the number of observations [48]. We used 50 spectra for each group, resulting in k = 50.

2.4.3. Accuracy Assessment

Following classification, performance metrics were chosen. The simplest metric is the overall accuracy (OA) or error rate. OA reflects the agreement between the observed and predicted classes and has the most direct interpretation. However, it does not provide information about the origin of an error [39]. Thus, we also calculated the Kappa statistic [49]. Kappa can take on values between −1 and 1; a value of 0 implies no agreement between the observed and predicted classes, while a Kappa of 1 indicates perfect concordance of the model prediction and the observed classes. A negative Kappa indicates that the prediction is in the opposite direction of the truth, but negative values for Kappa are not very common [39]. Landis and Koch [50] first defined the following standards for the strength of agreement for kappa coefficients: ≤0 = poor; 0.01–0.20 = slight; 0.21–0.40 = fair; 0.41–0.60 = moderate; 0.61–0.80 = substantial; and 0.81–1 = almost perfect. As class-specific metrics, we further calculated producer accuracy (PA) and user accuracy (UA) [51]. PA is the number of correctly-classified spectra for a class divided by the reference number of that class and, thus, represents the accuracy of the classification for a specific class. UA divides the number of correctly-classified spectra in a class by the total number of spectra within that class. A high UA means that spectra within that class can be reliably classified as belonging to that class. User accuracy is often termed to be a measure of reliability, which can be also interpreted as the agreement between repeated measurements within a class [52].

3. Results

3.1. Partial Least Squares Discriminant Analysis: Classification

Out of twenty method combinations, we found that a cross-validated PLS-DA without a preceding SNV transformation, applied on a full spectrum, delivered a Kappa value of 0.90 (Table 2), indicating an almost perfect concordance of the model prediction and the observed classes [50]. The second best model had a Kappa value of 0.77; here, we applied a test validation, including an SNV transformation applied on a full spectrum. Both models still indicate the substantial concordance of the model prediction and the observed classes.
Across all approaches, the overall accuracy OA ranged from 0.21–0.92. Within the training validation dataset, it ranged from 0.21–0.81, and within the cross validation dataset, it ranged from 0.34–0.92. The Kappa range across all approaches was 0.01–0.90. Within the training validation dataset, Kappa ranged from 0.01–0.77 and within the cross-validation from 0.17–0.90. By calculating the accuracies across all approaches to compare the tested subsets, it becomes obvious that the regions VIS (OA = 0.52, Kappa = 0.40), NIR (OA = 0.41, Kappa = 0.27) and SWIR I (OA = 0.30, Kappa = 0.13) were outperformed by the accuracies of the full spectrum (OA = 0.72, Kappa = 0.66) and SWIR II (OA = 0.64, Kappa = 0.55). The influence of an SNV transformation on our classification can be judged by comparing, again, the average accuracies for SNV (OA = 0.55, Kappa = 0.43) and non-SNV (OA = 0.49, Kappa = 0.36).
Table 2. Results deriving from a PLS-DA applied on the spectral database. The summary of the accuracy assessment statistics is displayed, including the best classification result written in bold.
Table 2. Results deriving from a PLS-DA applied on the spectral database. The summary of the accuracy assessment statistics is displayed, including the best classification result written in bold.
SNV/no SNVSegmentValidationOA95% CIp-ValueKappa
MinMax
SNVFullTraining0.810.670.91<0.050.77
SNVFullCross0.550.480.61<0.050.44
SNVVISTraining0.540.390.69<0.050.43
SNVVISCross0.520.460.58<0.050.40
SNVNIRTraining0.520.370.67<0.050.40
SNVNIRCross0.500.440.57<0.050.38
SNVSWIR ITraining0.210.100.350.01
SNVSWIR ICross0.370.310.43<0.050.21
SNVSWIR IITraining0.670.520.80<0.050.59
SNVSWIR IICross0.770.710.82<0.050.71
no SNVFullTraining0.600.450.74<0.050.51
no SNVFullCross0.920.880.95<0.050.90
no SNVVISTraining0.460.310.61<0.050.32
no SNVVISCross0.540.480.60<0.050.43
no SNVNIRTraining0.290.170.440.11
no SNVNIRCross0.340.280.40<0.050.17
no SNVSWIR ITraining0.270.150.420.09
no SNVSWIR ICross0.350.290.410.19
no SNVSWIR IITraining0.560.410.71<0.050.45
no SNVSWIR IICross0.550.480.61<0.050.44
OA, overall accuracy; CI, confidence interval; Full, full spectrum; VIS, visible spectral wavelength; NIR, near infrared spectral wavelength; SWIR I, short wave infrared spectral wavelength Section I; SWIR II, short wave infrared spectral wavelength Section II; T, training validation; C, cross validation; SNV, standard normal variate transformation applied; no SNV, no standard normal variate transformation applied.
Averaged spectral signatures per group (Figure 3) show that hairs and waxes differ strongly from the other groups in the VIS region, in particular, they have a higher reflectance around 550 nm. Bladder cells and the rough texture show similar pattern in the spectral region from 850 nm–1100 nm, while the other groups show very similar spectral patterns. Bladder cells also display an interesting feature of increased reflection around 1700 nm. Hairs generally have, over the whole spectrum, a high reflectance, while the smooth texture generally has the lowest reflectance across the whole spectrum (except in the VIS region). A clear ranking of spectra according to reflectance was not possible, because no spectrum is always higher than any other. Confusion matrices for the best two models (Table 3) revealed that the ETs bladder cells and waxes had in both datasets high user accuracies (UA), while their producer accuracies (PA) were moderate. Generally weak were the classes smooth and rough. Both had the lowest UA, depending on the dataset. The lowest absolute UA value had smooth (67%) in the test dataset. PA was lowest for smooth (80%) in the cross-validated dataset, while wax had the lowest PA in the test dataset (75%). The cross-validated dataset generally had higher accuracy values than the test dataset. In terms of PA, bladder cells, hairs and wax cover were consistently, highly accurate for the two different methods (Table 3), while rough and smooth often got mixed up, but also with the class hairs. Across the two methods, consistently high UA values are only found for hairs and bladder cells (Table 3). The other classes differ between the two methods with the lowest values found for smooth in the CV/no SNV and wax for the training/SNV confusion matrix. We checked the robustness of the PLS-DA models by comparing the ranking of the different method combinations presented in Table 3 based on the Kappa values. We inspected the rankings separately for the two validation steps, i.e., CV and training dataset. The order of the method combinations seemed relatively robust with respect to the subset used (full versus different wavelength regions). However, the effect of SNV transformation changed more often, in particular within the first ranks for each method. Hence, we can say that full spectra and SWIR II worked equally well in both methodologies, yet depending on the SNV transformation, the method combinations got ranked differently. With this unsteadiness based on the transformation, we decided to compare the confusion matrices for the best model of each validation technique. For a complete overview of the accuracy assessment, see the Supplementary Materials (Table S1).
Figure 3. Reflectance spectra. The image displays five averaged spectra that representing each functional group. The spectral regions (1–7), revealed by variable importance for projection (VIP) analysis, are highlighted in the white and unshaded areas. The numeration, 1–7, can be found directly above the center of each region.
Figure 3. Reflectance spectra. The image displays five averaged spectra that representing each functional group. The spectral regions (1–7), revealed by variable importance for projection (VIP) analysis, are highlighted in the white and unshaded areas. The numeration, 1–7, can be found directly above the center of each region.
Remotesensing 07 15862 g003
Table 3. Confusion matrices for the best two models, each recorded over the full spectrum. The upper table belongs to the cross-validated and no SNV transformation, the lower table to the test dataset and SNV transformation. The columns denote observed (ground truth) cases and the rows predicted (classification) cases. For each dataset, the user’s and producer’s accuracies are indicated.
Table 3. Confusion matrices for the best two models, each recorded over the full spectrum. The upper table belongs to the cross-validated and no SNV transformation, the lower table to the test dataset and SNV transformation. The columns denote observed (ground truth) cases and the rows predicted (classification) cases. For each dataset, the user’s and producer’s accuracies are indicated.
Ground Truth (GT)
DatasetBladderHairsRoughSmoothWax# of Classified SpectraUA
CV/no SNVClassificationBladder4701205094.0%
Hairs0470014897.9%
Rough0243104693.5%
Smooth3164415580.0%
Wax0002474995.9%
# of GT Spectra5050504949248
PA94.0%94.0%86.0%89.8%95.9%
Ground Truth (GT)
BladderHairsRoughSmoothWax# of Classified SpectraUA
Training/SNVClassificationBladder1001101283.3%
Hairs07100887.5%
Rough01710977.8%
Smooth00160785.7%
Wax020191275.0%
# of GT Spectra1010109948
PA100.0%70.0%70.0%66.7%100.0%
CV, cross-validation; Training, training validation; P, predicted; O, observed; SNV, standard normal variate transformation applied; no SNV, no standard normal variate transformation applied.

3.2. Variable Importance for Projection

To test if distinct bands had a greater impact on the classification than others, we calculated the VIP score. We modified the suggestion of Leung [47] to delete VIP values smaller than 0.8; we considered VIP values >1 to be significant. The calculation of the VIPs revealed seven spectral regions of increased importance for the classification (Figure 3). One large region covers the whole VIS region, including a small bit of the NIR region (400 nm–820 nm). Regions 2 (1325 nm–1347 nm), 3 (1394 nm–1460 nm), 4 (1717 nm–1733 nm), 5 (1899 nm–1985 nm), 6 (2288 nm–2320 nm) and 7 (2339 nm–2356 nm) all lie beyond the NIR region, while Regions 3 and 5 are located in the water absorption regions.

4. Discussion

The classification of our spectral signatures generally approved the visually-defined epidermal types (ET) as indicated by the high Kappa values. That implies species grouped by their epidermal structures (ES), like hairs, waxes; bladder cells, etc. are indeed spectrally separable; or in other words, ETs are similar in spectral properties within an ET, but do significantly differ between other ETs. The fact that this pattern also holds across species from different genera within one ET supports the idea of Ustin and Gamon [31] that optical types could be a new form of grouping species with a more functional interpretation.
Our findings carry weight because they confirm and expand the findings of earlier studies regarding the effects of ES on spectral signatures [9,10,19,53]. We identified five different functional groups based on their ES: hairs, bladder cells, rough and smooth textures and epicuticular waxes. Other studies defined similar groups, but there might have been differences in the interpretation of epidermal structures. ES of succulents were frequently described by botanists; e.g., Opel and Hammer [54] characterized three different epidermal types: trichrome/papillae, bladder cells and stomata alone for the genus Conophytum with numerous subdivisions. Ihlenfeldt and Hartmann [21] defined xeromorphic epidermis types for succulents of the Mesembryanthemaceae. This concept was advanced by Jürgens [23] and Niesler [24]. Yet, they did not make any spectral measurements. Tanner and Eller [14] analyzed the spectral signatures of the epidermal structures of the Mesembryanthemaceae, i.e., bladder cells and smooth leaf surfaces. However, they did not find any conspicuous differences in the spectral signatures. Our averaged spectra per ET (Figure 3) show a similar pattern for bladder cells and smooth surfaces in the VIS region that were studied by Tanner and Eller [14]. However, bladder cells clearly differ from apparently smooth surfaces in the spectral region between 1600 nm and 1700 nm, a spectral region that was not accessible in the days of Tanner and Eller. Hence, it can be assumed that our findings provide an extension of previous findings.
Leaf pubescence increases the reflectance of global radiation and is very effective between 400 nm and 700 nm [9]. Besides affecting reflectance, hairs on leaf surfaces help to diminish transpiration and support thermal convection processes [10]. Eller and von Willert [10,19] both stated that hairs have no effect on spectral regions larger than 750 nm, though we found six spectral regions beyond 750 nm where hairs actually had an impact on the spectral classification. In fact, our ET hairs had, except from 750 nm–1200 nm, constantly high reflectance values compared to other ETs. Epicuticular waxes on leaves are, like hairs, able to increase the spectral reflectance between 400 nm and 750 nm [11]. Furthermore, enclosed oxalate crystals inside of wax layers can amplify the reflection of photons up to 33% [10,21]. Another process of leaf surfaces influencing an ecological strategy was shown by Haberland [55]. He found that papillae have a lens effect on incoming photons. Tanner and Eller [14] and von Willert [19] specified this and explained that due to the lens effect, photons are able to reach assimilation tissues located in deeper cell layers of leaves. Lastly, plants may exhibit papillae on their leaf surfaces to form a rough epidermal structure; as a consequence, they are able to modify the convective energy transfer by producing Reynolds’ turbulences [10,19,56]. That bladder cells also interfere with photons was shown by Tanner and Eller [14]. This time, the reflectance between 400 nm and 1200 nm was lowered, and the absorption was increased.
All of these examples emphasize a strong connection between spectral signatures and plants’ epidermal structures, which represent different adaptions of plants to their environment. However, these examples also just focus on the IR (infra-red) region of the electromagnetic spectrum. Today, it is well known that spectral responses of the IR region are mainly caused by plant pigments of the palisade mesophyll and leaf cell structures. Beyond this region (1300 nm–2600 nm), it is known that the plant’s water content has a great effect on the reflected signature [40]. We are happy to confirm reflection patterns in the IR region, but it is even more gratifying that our applied VIP calculation revealed the already mentioned bands, which supports our classification of epidermal types (ETs). This finding shows that the bands beyond the IR region might contain not just information regarding leaf water content, but also regarding leaf surface structures.
Our classification had also several weaknesses. First of all, biochemical properties or varying leaf angles during the measurement could have an impact on the spectral response, as well, which were not taken into account in this study. Our focus was to show that, basically, one can separate ETs spectrally from each other and put them into categories. Still, the biochemical components would be worthwhile to include in a future multivariate regression model combined with functional groupings as our ETs. This would better allow discerning the different effects of biochemical and biophysical properties. It would also be useful to test the regions, detected by the VIP calculation, for a classification to increase accuracies. Again, the spectral differences between the groups cannot be explained by the reflectance of epidermal structures alone. Further studies are required to understand the proportion to which ES have an impact on the detected regions. A new technology called hyperspectral microscopy has a high potential for disentangling the biochemical and biophysical effects on spectral signatures of leaves [57]. It combines high resolution microscopy, even on the nanoscale, with hyperspectral resolution from 400 nm to approximately 1000 nm.
In addition, we had problems classifying the group “smooth”, which had the lowest user accuracy. Initially, we classified the ETs with a 20-fold magnifier, but the spectral classification was able to detect details that we could not detect with the magnifier. Roberts [58] stated that spectral features change with the scale of the measurement. At each larger scale, new features are included in the measurements that were either absent or of minor importance at lower scales. Hartmann [59,60] summarized properties of the family Aizoaceae and described epidermal structures of the plants we investigated in this study. Having this source, we once more compared the plants of the group “smooth”, Glottiphyllum longum, Malotigena frantiskae-niederlovae, Aptenia haeckeliana, Schlechteranthus hallii and Bergeranthus scapiger, with the classification by Hartmann [59,60], who describes Glottiphyllum longum and Aptenia haeckeliana as exhibiting inconspicuous bladder cells. Schlechteranthus hallii and Bergeranthus scapiger are described to have a wax layer besides slightly elevated cells (more or less rough) on the leaf surface. These descriptions conceivably explain why the class smooth could not be perfectly classified as a single ET. However, the species were correctly classified into other respective groups, i.e., bladder cells, rough and smooth (cf. Table 2) groups, in which these plants belong instead. The positive aspect of this mistake is that the spectral classification corrected our initial classification and, therefore, emphasized the usefulness of our classification.

5. Conclusions

With this study, we aimed at contributing to the concept of plant functional classification, i.e., the “optical types” by Ustin and Gamon [31], which directly links ecological pattern and remotely-sensed spectral data. We demonstrated that five different functional types of epidermal leaf structures sampled across 25 species from a single plant family could be distinguished based on leaf spectral signatures. We found six spectral regions, all beyond 1100–1300 nm, that had a significant impact on the classification of leaf epidermal structures according to the variable of importance (VIP) measure. Epidermal structures, like bladder cells, hairs and waxes, were accurately classified (producer accuracy between 70% and 100%), while rough and smooth structures had lower classification rates (producer accuracy between 66% and 89%). In particular, we could demonstrate that the morphological pattern can translate directly into spectral signals. However, the effects of biochemical constituents, such as pigments, which are also known to influence spectral signatures, were not analyzed in this study. Future studies should take into account biophysical, i.e., leaf epidermal structures and biochemical information for discriminating plant species based on leaf spectral signals.

Supplementary Files

Supplementary File 1

Acknowledgments

We thank Katrin Baumgärtner for taking care of the Aizoaceae living collection and her advice on selecting appropriate plants, Ute Schmiedel and Ingeborg Niesler for their knowledge on epicuticular structures of Aizoaceae and Klaus von Bösekom for taking the pretty photos. Funding was provided by DFG (QUEEN, OL 417/1-1, TH 1295/4-1).

Author Contributions

André Große-Stoltenberg assisted in setting up the experiments regarding the handheld spectrometer and in processing the spectral data. Norbert Jürgens contributed to the manuscript and granted access to the living collection of Aizoaceae. Jens Oldeland conceived of the idea of the study, edited the manuscript and designed the statistical analysis. Rene Hans-Jürgen Heim conducted the experiments and the statistical analysis, performed the literature research and wrote the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sims, D.A.; Gamon, J.A. Relationships between leaf pigment content and spectral reflectance across a wide range of species, leaf structures and developmental stages. Remote Sens. Environ. 2002, 81, 337–354. [Google Scholar] [CrossRef]
  2. Blackburn, G.A. Hyperspectral remote sensing of plant pigments. J. Exp. Bot. 2007, 58, 855–867. [Google Scholar] [CrossRef] [PubMed]
  3. Hellmann, C.; Grosse-Stoltenberg, A.; Laustro, V.; Oldeland, J.; Werner, C. Retrieving nitrogen isotopic signatures from fresh leaf reflectance spectra: Disentangling delta δ15N from biochemical and structural leaf properties. Front. Plant Sci. 2015, 6, 307. [Google Scholar] [CrossRef] [PubMed]
  4. Asner, G.P.; Martin, R.E.; Tupayachi, R.; Anderson, C.B.; Sinca, F.; Carranza-Jimenez, L.; Martinez, P. Amazonian functional diversity from forest canopy chemical assembly. Proc. Natl. Acad. Sci. USA 2014, 111, 5604–5609. [Google Scholar] [CrossRef] [PubMed]
  5. Lehmann, J.R.K.; Grosse-Stoltenberg, A.; Romer, M.; Oldeland, J. Field spectroscopy in the vnir-swir region to discriminate between mediterranean native plants and exotic-invasive shrubs based on leaf tannin content. Remote Sens. 2015, 7, 1225–1241. [Google Scholar] [CrossRef] [Green Version]
  6. Kalacska, M.; Bohman, S.; Sanchez-Azofeifa, G.A.; Castro-Esau, K.; Caelli, T. Hyperspectral discrimination of tropical dry forest lianas and trees: Comparative data reduction approaches at the leaf and canopy levels. Remote Sens. Environ. 2007, 109, 406–415. [Google Scholar] [CrossRef]
  7. Sims, D.A.; Gamon, J.A. Estimation of vegetation water content and photosynthetic tissue area from spectral reflectance: A comparison of indices based on liquid water and chlorophyll absorption features. Remote Sens. Environ. 2003, 84, 526–537. [Google Scholar] [CrossRef]
  8. Da Luz, B.R. Attenuated total reflectance spectroscopy of plant leaves: A tool for ecological and botanical studies. New Phytol. 2006, 172, 305–318. [Google Scholar] [CrossRef] [PubMed]
  9. Ehleringer, J.R.; Bjorkman, O. Pubescence and leaf spectral characteristics in a desert shrub, encelia farinosa. Oecologia 1978, 36, 151–162. [Google Scholar] [CrossRef]
  10. Eller, B.M. Epidermis and the spectral properties of plant-surfaces. Ber. Deut. Bot. Ges. 1985, 98, 465–475. [Google Scholar]
  11. Eller, B.M.; Willi, P. The significance of leaf pubescence for the absorption of global radiation by tussilago farfara l. Oecologia 1977, 29, 179–187. [Google Scholar] [CrossRef]
  12. Ehleringer, J.R.; Mooney, H.A. Leaf hairs: Effects on physiological activity and adaptive value to a desert shrub. Oecologia 1978, 37, 183–200. [Google Scholar] [CrossRef]
  13. Mulroy, T.W. Spectral properties of heavily glaucous and non-glaucous leaves of a succulent rosette-plant. Oecologia 1979, 38, 349–357. [Google Scholar] [CrossRef]
  14. Tanner, V.; Eller, B.M. Epidermis structure and its significance for the optical-properties of leaves of the mesembryanthemaceae. J. Plant Physiol. 1986, 125, 285–294. [Google Scholar] [CrossRef]
  15. Huntjr, E.; Rock, B. Detection of changes in leaf water content using near- and middle-infrared reflectances. Remote Sens. Environ. 1989, 30, 43–54. [Google Scholar] [CrossRef]
  16. Martin, M.E.; Aber, J.D. High spectral resolution remote sensing of forest canopy lignin, nitrogen, and ecosystem processes. Ecol. Appl. 1997, 7, 431–443. [Google Scholar] [CrossRef]
  17. Martin, M.E.; Newman, S.D.; Aber, J.D.; Congalton, R.G. Determining forest species composition using high spectral resolution remote sensing data. Remote Sens. Environ. 1998, 65, 249–254. [Google Scholar] [CrossRef]
  18. Carter, G.A.; Knapp, A.K. Leaf optical properties in higher plants: Linking spectral characteristics to stress and chlorophyll concentration. Am. J. Bot. 2001, 88, 677–684. [Google Scholar] [CrossRef] [PubMed]
  19. Von Willert, D.J. Life Strategies of Succulents in Deserts: With Special Reference to the Namib Desert; CUP Archive: New York, NY, USA, 1992; p. 372. [Google Scholar]
  20. Myers, N.; Mittermeier, R.A.; Mittermeier, C.G.; da Fonseca, G.A.; Kent, J. Biodiversity hotspots for conservation priorities. Nature 2000, 403, 853–858. [Google Scholar] [CrossRef] [PubMed]
  21. Ihlenfeldt, H.D.; Hartmann, H.E.K. Leaf surfaces in mesembryanthemaceae. In The Plant Cuticle; Cutler, D.F., Alvin, K.L., Price, C.E., Eds.; Academic Press: London, UK, 1982; pp. 397–423. [Google Scholar]
  22. List, T.P. The Plant List (2013). Version 1.1. Available online: http://www.theplantlist.org/ (accessed on 1 May 2014).
  23. Jürgens, N. Untersuchungen zur Ökologie sukkulenter Pflanzen des südlichen Afrika. Mitt. Inst. Allg. Bot. Hamburg 1986, 21, 139–365. [Google Scholar]
  24. Niesler, I.M. Untersuchungen zur Ontogenie und zum Bau der Blätter bei Mesembryanthemaceae Fenzl unter besonderer Berücksichtigung der Epidermis. Ph.D. Thesis, Universität Hamburg, Hamburg, Germany, 1998. [Google Scholar]
  25. Clark, M.L.; Roberts, D.A.; Clark, D.B. Hyperspectral discrimination of tropical rain forest tree species at leaf to crown scales. Remote Sens. Environ. 2005, 96, 375–398. [Google Scholar] [CrossRef]
  26. Ghiyamat, A.; Shafri, H.Z.M.; Mandiraji, G.A.; Shariff, A.R.M.; Mansor, S. Hyperspectral discrimination of tree species with different classifications using single- and multiple-endmember. Int. J. Appl. Earth Observ. Geoinf. 2013, 23, 177–191. [Google Scholar] [CrossRef] [Green Version]
  27. Prospere, K.; McLaren, K.; Wilson, B. Plant species discrimination in a tropical wetland using in situ hyperspectral data. Remote Sens. 2014, 6, 8494–8523. [Google Scholar] [CrossRef]
  28. He, Y.; Li, X.; Deng, X. Discrimination of varieties of tea using near infrared spectroscopy by principal component analysis and bp model. J. Food Eng. 2007, 79, 1238–1242. [Google Scholar] [CrossRef]
  29. Edelmann, A.; Diewok, J.; Schuster, K.C.; Lendl, B. Rapid method for the discrimination of red wine cultivars based on mid-infrared spectroscopy of phenolic wine extracts. J. Agric. Food Chem. 2001, 49, 1139–1145. [Google Scholar] [CrossRef] [PubMed]
  30. Rao, N.R.; Garg, P.K.; Ghosh, S.K. Development of an agricultural crops spectral library and classification of crops at cultivar level using hyperspectral data. Precis. Agric. 2007, 8, 173–185. [Google Scholar] [CrossRef]
  31. Ustin, S.L.; Gamon, J.A. Remote sensing of plant functional types. New Phytol. 2010, 186, 795–816. [Google Scholar] [CrossRef] [PubMed]
  32. Duckworth, J.C.; Kent, M.; Ramsay, P.M. Plant functional types: An alternative to taxonomic plant community description in biogeography? Prog. Phys. Geogr. 2000, 24, 515–542. [Google Scholar] [CrossRef]
  33. Homolova, L.; Maenovsky, Z.; Clevers, J.G.P.W.; Garcia-Santos, G.; Schaeprnan, M.E. Review of optical-based remote sensing for plant trait mapping. Ecol. Complex. 2013, 15, 1–16. [Google Scholar] [CrossRef]
  34. Houborg, R.; Fisher, J.B.; Skidmore, A.K. Advances in remote sensing of vegetation function and traits. Int. J. Appl. Earth Observ. Geoinf. 2015, 43, 1–6. [Google Scholar] [CrossRef]
  35. Dorigo, W.; Bachmann, M.; Heldens, W. AS Toolbox & Processing of Field Spectra-user’s manual; German Aerospace Center, DLR-DFD. Imaging Spectroscopy Group: Wessling, Germany, 2006. [Google Scholar]
  36. Fearn, T.; Riccioli, C.; Garrido-Varo, A.; Guerrero-Ginel, J.E. On the geometry of snv and msc. Chemom. Intell. Lab. Syst. 2009, 96, 22–26. [Google Scholar] [CrossRef]
  37. Barnes, R.J.; Dhanoa, M.S.; Lister, S.J. Standard normal variate transformation and de-trending of near-infrared diffuse reflectance spectra. Appl. Spectrosc. 1989, 43, 772–777. [Google Scholar] [CrossRef]
  38. Sanchez, G.; Determan, C.; Sanchez, M.G. R Package DiscriMiner (Version 0.1–29). Available online: https://cran.r-project.org/web/packages/DiscriMiner/ (accessed on 13 May 2015).
  39. Kuhn, M.; Johnson, K. Applied Predictive Modeling; Springer: Berlin, Germany, 2013. [Google Scholar]
  40. Thenkabail, P.S.; Lyon, J.G.; Huete, A. Hyperspectral Remote Sensing of Vegetation; CRC Press: Boca Raton, FL, USA, 2011; p. 766. [Google Scholar]
  41. Chong, I.G.; Jun, C.H. Performance of some variable selection methods when multicollinearity is present. Chemom. Intell. Lab. Syst. 2005, 78, 103–112. [Google Scholar] [CrossRef]
  42. Wold, S.; Kettaneh, N.; Tjessem, K. Hierarchical multiblock pls and pc models for easier model interpretation and as an alternative to variable selection. J. Chemom. 1996, 10, 463–482. [Google Scholar] [CrossRef]
  43. Wold, S.; Martens, H.; Wold, H. The multivariate calibration problem in chemistry solved by the PLS method. In Matrix Pencils; Kågström, B., Ruhe, A., Eds.; Springer: Berlin, Germany, 1983; pp. 286–293. [Google Scholar]
  44. Barker, M.; Rayens, W. Partial least squares for discrimination. J. Chemom. 2003, 17, 166–173. [Google Scholar] [CrossRef]
  45. Perez-Enciso, M.; Tenenhaus, M. Prediction of clinical outcome with microarray data: A partial least squares discriminant analysis (PLS-DA) approach. Hum. Genet 2003, 112, 581–592. [Google Scholar] [PubMed]
  46. Fisher, R.A. The use of multiple measurements in taxonomic problems. Ann. Eugen. 1936, 7, 179–188. [Google Scholar] [CrossRef]
  47. Leung, K.; Cheong, F.; Cheong, C.; O’Farrell, S.; Tissington, R.; Ou, C.M. A Comparison of Variable Selection Techniques for Credit Scoring. In Proceedings of the 7th International Conference on Computational Intelligence in Economics and Finance, Taiwan, 4 Augest 2008.
  48. Efron, B.; Tibshirani, R. Improvements on cross-validation: The 632+ bootstrap method. J. Am. Stat. Assoc. 1997, 92, 548–560. [Google Scholar]
  49. Cohen, J. A coefficient of agreement for nominal scales. Educ. Psychol. Meas. 1960, 20, 37–46. [Google Scholar] [CrossRef]
  50. Landis, J.R.; Koch, G.G. The measurement of observer agreement for categorical data. Biometrics 1977, 33, 159–174. [Google Scholar] [CrossRef] [PubMed]
  51. Story, M.; Congalton, R.G. Accuracy assessment—A users perspective. Photogramm Eng. Remote Sens. 1986, 52, 397–399. [Google Scholar]
  52. Jones, H.G.; Vaughan, R.A. Remote Sensing of Vegetation: Principles, Techniques, and Applications; Oxford University Press: New York, NY, USA, 2010; p. 381. [Google Scholar]
  53. Opel, M.R. Leaf anatomy of conophytum n.E.Br. (aizoaceae). Haseltonia 2005, 11, 27–52. [Google Scholar] [CrossRef]
  54. Opel, M.R.; Hammer, S.A. Conophytum leaf structures. In SA Hammer, Dumpling and His Wife: New Views of the Genus Conophytum; EAE Creative Color Ltd.: Norwich, UK, 2002; pp. 300–321. [Google Scholar]
  55. Haberland, G. Physiologische pflanzenanatomie, Leipzig 1909, IV. Aufl. S 1909, 557, 671. [Google Scholar]
  56. Barthlott, W.; Wollenweber, E. Zur Feinstruktur, Chemie und taxonomischen Signifikanz epicuticularer Wachse und ähnlicher Sekrete. Trop. Subtrop. Pflanzenwelt. 1981, 32, 67. [Google Scholar]
  57. Roth, G.A.; Tahiliani, S.; Neu-Baker, N.M.; Brenner, S.A. Hyperspectral microscopy as an analytical tool for nanomaterials. Wiley Interdiscip. Rev. Nanomed. Nanobiotechnol. 2015, 7, 565–579. [Google Scholar] [CrossRef] [PubMed]
  58. Roberts, D.A.; Ustin, S.L.; Ogunjemiyo, S.; Greenberg, J.; Dobrowski, S.Z.; Chen, J.Q.; Hinckley, T.M. Spectral and structural measures of northwest forest vegetation at leaf to landscape scales. Ecosystems 2004, 7, 545–562. [Google Scholar] [CrossRef]
  59. Hartmann, H.E.K. Illustrated Handbook of Succulent Plants: Aizoaceae A–E; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2002; p. 456. [Google Scholar]
  60. Hartmann, H.E.K. Illustrated Handbook of Succulent Plants: Aizoaceae F–Z; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2002; p. 422. [Google Scholar]

Share and Cite

MDPI and ACS Style

Heim, R.H.-J.; Jürgens, N.; Große-Stoltenberg, A.; Oldeland, J. The Effect of Epidermal Structures on Leaf Spectral Signatures of Ice Plants (Aizoaceae). Remote Sens. 2015, 7, 16901-16914. https://doi.org/10.3390/rs71215862

AMA Style

Heim RH-J, Jürgens N, Große-Stoltenberg A, Oldeland J. The Effect of Epidermal Structures on Leaf Spectral Signatures of Ice Plants (Aizoaceae). Remote Sensing. 2015; 7(12):16901-16914. https://doi.org/10.3390/rs71215862

Chicago/Turabian Style

Heim, René Hans-Jürgen, Norbert Jürgens, André Große-Stoltenberg, and Jens Oldeland. 2015. "The Effect of Epidermal Structures on Leaf Spectral Signatures of Ice Plants (Aizoaceae)" Remote Sensing 7, no. 12: 16901-16914. https://doi.org/10.3390/rs71215862

APA Style

Heim, R. H. -J., Jürgens, N., Große-Stoltenberg, A., & Oldeland, J. (2015). The Effect of Epidermal Structures on Leaf Spectral Signatures of Ice Plants (Aizoaceae). Remote Sensing, 7(12), 16901-16914. https://doi.org/10.3390/rs71215862

Article Metrics

Back to TopTop