LAMINA: a tool for rapid quantification of leaf size and shape parameters
© Bylesjö et al; licensee BioMed Central Ltd. 2008
Received: 07 May 2008
Accepted: 22 July 2008
Published: 22 July 2008
An increased understanding of leaf area development is important in a number of fields: in food and non-food crops, for example short rotation forestry as a biofuels feedstock, leaf area is intricately linked to biomass productivity; in paleontology leaf shape characteristics are used to reconstruct paleoclimate history. Such fields require measurement of large collections of leaves, with resulting conclusions being highly influenced by the accuracy of the phenotypic measurement process.
We have developed LAMINA (Leaf shApe deterMINAtion), a new tool for the automated analysis of images of leaves. LAMINA has been designed to provide classical indicators of leaf shape (blade dimensions) and size (area), which are typically required for correlation analysis to biomass productivity, as well as measures that indicate asymmetry in leaf shape, leaf serration traits, and measures of herbivory damage (missing leaf area). In order to allow Principal Component Analysis (PCA) to be performed, the location of a chosen number of equally spaced boundary coordinates can optionally be returned.
We demonstrate the use of the software on a set of 500 scanned images, each containing multiple leaves, collected from a common garden experiment containing 116 clones of Populus tremula (European trembling aspen) that are being used for association mapping, as well as examples of leaves from other species. We show that the software provides an efficient and accurate means of analysing leaf area in large datasets in an automated or semi-automated work flow.
Leaves are of fundamental importance to plants, representing the power generation facility and aerial environmental sensing units of plants, and by extension ultimately provide the energy for sustaining most terrestrial species on earth. A number of genes known to affect meristomatic pattern formation (e.g. AS1 and WUS, KNOX and CLV, see  for a review of leaf development), the rate of leaf primordia initiation  and that contribute to the determination of leaf length (ROT3 , LNG ) and width (AN ) have now been identified: less is known about the determination of leaf size currently. Despite these advances, it remains clear that leaf area development is a highly complex process that is influenced by genetic, hormonal and environmental factors. Quantitative Trait Loci (QTL) mapping of leaf development and leaf size and shape indicators suggests that these traits are under the control of many genes [6–15], with relatively few genes identified to date . To advance the current understanding of leaf area development and final dimension determination requires the ability to phenotype large collections of leaves from QTL mapping populations, natural populations and forward genetic screens to identify and quantify loci/mutations influencing leaf characteristics. As well as being important to the fields of genetics, physiology, plant breeding and developmental biology, leaf shape parameters are also important as a means of reconstructing historical paleoclimate conditions [16, 17], where information on leaf serration (depth and presence/absence) is used to accurately reconstruct past mean annual temperature [18, 19]. Leaf size and shape parameters (physiognomy) were initially quantified using gridded paper, where a count of squares was used to measure leaf size, or through the development of allometric relationships between length, width and area, with length typically being measured and later used to calculate area using a regression model. This approach can work well within a single species but works poorly when applied to mapping populations, where segregation can lead to extensive variation in both leaf area and shape traits. It is equally inappropriate for forward genetic screens to identify leaf phenotypes, where the induced phenotypic changes are unpredictable. For many species, field-portable leaf scanning equipment can be used to measure leaf area and blade dimensions. However, such equipment cannot be used on large leaves and works poorly on species such as Arabidopsis thaliana due to small leaf area and the proximity of leaves to the soil. Such equipment is also often limited in the range of measurements provided and, as no digital image is captured, retrospective re-analysis using, for example, new software tools is not possible. More recently, methods have concentrated on the capturing of digital images of leaves (or fossils) with subsequent analysis using digital image analysis tools. A number of such tools already exist, but none of the currently available software was able to fulfil our needs. ImageJ  is a widely used application for the analysis of biological images and can be used to analyse area and blade dimensions. However, automated analysis is hard to achieve, as is the simultaneous measurement of area and blade dimensions when leaves are not square within the image. ImageJ offers no method to quantify leaf serration. The development of tools for measuring leaf area was reported in  and , however they offer little to extend the capabilities of ImageJ. More recently  reported the development of LeafAnalyser, which is an excellent tool to facilitate PCA analysis of leaf shape parameters. However, this tool does not report the type of dimensions that are typically required by plant breeders, physiologists, geneticists or palaeontologists and the software was not released as open source, negating the possibility of further development by the community. We additionally found that the implemented thresholding frequently required per-image manual adjustment, making the automated, rapid analysis of leaves more time consuming. We were interested in measuring basic leaf dimension parameters (area, length, width) as well as measures of leaf shape, symmetry, serration number and depth and the missing area within a leaf (as a measure of damage by biting herbivores) in a collection of naturally occurring clones of Populus tremula, the Swedish Aspen (SwAsp) collection, that are being grown in common garden experiments in the south and north of Sweden  and that are being used for association mapping . This species has well defined, characteristic leaf serrations that we had visually observed to show variation between clones within the SwAsp collection. We were therefore interested to see to what degree leaf serration was under genetic control. This required a rapid and reproducible method of quantifying leaf size and shape parameter traits as well as serration characteristics. As was reported in , we were also interested to see how well PCA could be used to describe the variation in leaf area characteristics within this collection of trees.
Main computational steps
The computational processes involved can be described in the following sequential steps.
1. Thresholding. As an initial step, global thresholding is performed to find candidate picture elements (pixels) that putatively represent leaves. In the thresholding process, all pixel intensities are reduced from the typical grayscale range of 0–255 to either 0 (off; pixel is background) or 1 (on; pixel potentially belongs to a leaf object). As input, the inverse of the blue channel intensities are used rather than the entire RGB image. The rationale behind this strategy is that while leaves can be green, orange, red or even black, they are very rarely blue. On a white background, non-blue objects can, with high accuracy, be distinguished from the background using global thresholding.
The process of identifying a suitable global threshold can be performed either manually or automatically. In manual thresholding, the user specifies an arbitrary value t in the range 0–255, where pixels with intensities less than t will be set to 0 (background) and pixels with intensities equal to or greater than t will be set to 1 (putative leaves). The automatic thresholding procedure, on the other hand, tries to automatically determine a value of t that minimizes the variance of the thresholded image . This procedure is generally suitable for images where objects have fairly well-behaved shapes, which is true for most leaf objects (see Artemisia annua section for an exception). The automated search procedure can be greedy, in which a local minima is found based on a greedy search starting from the mean value of the starting image. Alternatively, the procedure can be exhaustive, in which the entire range 0–255 is searched for the value t that minimizes the variance of the thresholded image. The latter procedure is generally more accurate but also considerably slower.
2. Segmentation. Posterior to thresholding, the input image has been reduced to a binary image containing pixels that are either background (0 = off) or potential leaf objects (1 = on). The task of segmentation is to group nearby pixels into segments (objects) that may potentially represent leaves. The segmentation starts by assigning an arbitrary on pixel as the current segment. The segment is then iteratively extended with neighbouring, unassigned on pixels (including diagonal pixels) until no more neighbouring on pixels exist. This procedure is repeated until all on pixels have been assigned to an object.
3. Filtering. Due to measurement noise and presence of contaminants in the image, some objects will not represent actual leaves. To remove dubious objects, filtering can be performed based on both the area of each object (to remove objects that are too small) and based on the density of each object (to remove e.g. black frames surrounding the image). The default filtering is non-stringent and will only remove the smallest objects, likely to represent contaminations in the image.
4. Object boundaries. The boundary of an object is defined by the set of on pixels where at least one neighbour of each on pixel is an off pixel (i.e. pixels on the surface of the object). The identification of the boundary pixels is a straightforward computational process. However, in order to simplify the subsequent steps, the adjacent boundary pixels are internally arranged sequentially (sorted) within each object. This procedure requires that distances are calculated between all boundary pixels and can be time-consuming for highly irregular surfaces, e.g. Artemisia annua images.
5. Cavities. Cavities in the leaf objects can be present due to e.g. biting herbivore damage, which implies that identification and measurement are of interest. A cavity is by definition surrounded by a boundary region that is unconnected to the outer boundary of the object. This distinctive characteristic is used to identify the cavities, seen as 'kinks' in the distances between neighboring boundary pixels. The off pixels that can only be connected to the cavity (inner) boundary define the cavity area. In this sense, cavities are defined as missing leaf area (holes) within the leaf lamina and do not account for herbivory starting at the edge of the leaf, which is computationally more difficult to quantify as it would require retrospective calculation of where the leaf boundary was previously. It is equally hard to distinguish herbivory or wounding at the leaf boundary from serrations. This represents an obvious area of future extension of LAMINA, but is not a trivial task.
6. Serrations and indents. Starting from a boundary pixel, the longest straight line that can be formed without crossing the object formed by non-boundary pixels is sought. The intermediate region between two serrations defines an indent. This is implemented in practice by connecting the starting boundary pixel with boundary pixels of increasing distances until a non-connectable boundary pixel is found. The last connectable pixel, i.e. one that can be connected by a straight line without crossing the object, is the next serration point. The process starts again using the latest serration point as the starting pixel. To allow for small variations in the boundary shape, a consecutive sequence of k non-connectable pixels are allowed before stopping. The parameter k can be adjusted by the user and determines the overall sensitivity of the serration identification algorithm.
7. Indent depths. Each indent is surrounded by two serrations that can be connected by a straight line. The indent depth is measured as the longest line to the base of the indent while being perpendicular to the straight line connecting the surrounding serrations. Due to the discrete nature of images, it is not always possible to achieve perfect perpendicularity, and hence a slight discrepancy in this angle is allowed.
8. Boundary coordinates. From the boundary pixels of each object, a fixed number of boundary coordinate points can optionally be identified. These are defined as equally spaced points around the surface of the object. The boundary coordinates are normalised against the centre coordinate of the object to make the measurements independent of the position of the object in the image.
Output from LAMINA
After processing, LAMINA outputs cropped image files representing the identified objects after thresholding and segmentation. This allows the user to have a record of the results of the image analysis process (Figure 1B shows example cropped images within the LAMINA user iterface and Figure 1C an example generated cropped image). Furthermore, a number of quantitative measurements of the leaves are generated. This includes the leaf area, height, width, circularity, number of serrations, indent widths and depths as well as the boundary coordinates (normalised against leaf centre). For parameters that summarise several measurements, the output includes the mean, median and standard deviation.
Image measurements do not generally contain any information regarding the actual size of the image. In order to convert the pixel-based distances and areas in the leaf image into real quantitative measures, scale calibration has to be performed. The aim of the calibration is to determine the actual size of one pixel in millimeters (mm) and is optimally run once, to find the conversion ratio. LAMINA requires a calibration image to perform this calculation, containing one coloured object (not black) of known size on a white background. Ideally this object should fill the majority of the image area to maximize the accuracy of the calibration. After determining the measured pixel size of the image, and by manual input of the actual size in mm, the pixel-to-mm ratio can be determined and employed for all subsequent image calculations.
Example applications of LAMINA
Exploring leaf physiognomy in the SwAsp Populus tremula collection
Full details of the common garden experiment can be found in . Five leaves per replicate of each clone were sampled on a single day in early August 2007 into paper bags and later scanned using a Canon CanoScan 4400 F A4 flatbed scanner at a resolution of 300 dpi. A 40 × 50 mm yellow rectangle of card was scanned and used for scale calibration. Images were saved as jpeg files. The majority of genets (clones) were represented by four clonal replicates. Our sampling strategy was to select five random leaves from different heights on each replicate as we wanted to know how plastic (i.e. variable) leaf area was within and between both genets and intra-genet clonal replicates. The only criterion applied was that leaves should be mature and should not be from the terminal stem, as these leaves are of a fundamentally different nature in aspens. Images were analysed in LAMINA in a semi-automated work flow to allow for corrections to the orientation of leaves within the image. Default settings were used for all parameters except for the serration detection pixel threshold, where 22 was used. The centre line of each leaf was adjusted where required before proceeding to the next image. In total, 412 images containing 1879 leaves were analysed, with the LAMINA analysis taking 1 working day (8 h).
A random set of 50 leaves were scanned and analysed using ImageJ  and LAMINA. For the ImageJ analysis, images were imported as an image stack, at which point they were transformed to 8 bit (greyscale) and then thresholded using a value of 150 to produce a binary image. The trace tool was then used to select each leaf and the Measure tool used to record the selected area. The scale was set using the line tool to define a known distance using the same calibration image used for the LAMINA analysis. Data generated were analysed and visualised in R . ANOVA tests were performed using the nlme package to test for clone within population and population effects. Principal Components Anlaysis (PCA) was performed in SIMCA P (v11.3, Umetrics, Sweden).
Benchmarking LAMINA using the complex leaves of Artemisia annua
A. annua leaves are highly complex and we deemed them to serve as a comprehensive test of the ability of LAMINA to extract and reliably quantify leaf area and dimension traits. We therefore undertook a more detailed method comparison using either glasshouse-grown or field-grown genotypes of A. annua. One mature leaf from six genotypes grown in a glasshouse was used to compare leaf area meter data to LAMINA. The area of each leaf was measured using a LI-COR LI-3100 Area Meter (LI-COR Environmental, Nebraska, USA) and the same leaves were scanned using an HP Scanjet 3570 c A4 flatbed scanner at 300 dpi. A 100 × 1 mm bar was scanned and used for scale calibration. Three mature leaves from 29 genotypes grown outside at Stockbridge Technology Centre, Cawood, North Yorkshire, U.K. were scanned using the same scanner and used for area analysis in LAMINA and ImageJ. Leaves 20, 21 and 22 (counting down from the top of the plant) were sampled in October 2007.
LAMINA analysis was performed using the following settings: Manual threshold value of 150, no serration detection.
ImageJ analysis was performed by transforming images to 8 bit (greyscale) and then thresholded using a value of 150 to produce a binary image. A polygon was then drawn around a leaf and the Analyze Particles tool used to calculate the area represented by leaf pixels. The scale was set by scanning a standard ruler and using the line tool to define a known distance. The use of the pixel analysis method, rather than the more automated method used for the aspen leaves, was required due to the complex shape of the A. annua leaves. However, this method increases the chance of any noise artifacts in the scanned image being included in the measurement calculations.
Testing LAMINA using species with diverse leaf shapes
In order to ensure that LAMINA functioned for a diverse range of species, we sampled leaves of a number of common European tree species as well as various poplar species and A. thaliana. One to three leaves per species were analysed to ensure that leaves were reliably extracted from the scanned images. All images were scanned as for the SwAsp trees. Additionally, the jpeg format images used as example applications in  and  were downloaded and analysed using LAMINA in order to benchmark our software against these other packages.
Results and discussion
Using LAMINA to explore leaf traits in the SwAsp collection
Overview of leaf size and shape traits in the SwAsp trees
Number of serrations
Benchmarking LAMINA against ImageJ and a leaf area meter using Artemisia annua
Leaves from A. annua plants were by far the most complex in structure of those that we used for testing and developing LAMINA. We therefore examined the results generated for these leaves in more detail as a means of benchmarking LAMINA. Leaf area is an important trait in A. annua as this medicinal crop produces artemisinin, used in anti-malaria drugs, in glandular trichomes found predominantly on the leaf surface. Natural variation in A. annua leaves is being studied using QTL mapping and association studies, while induced mutations and phenotypes are being identified using forward and reverse genetic screens with all of these approaches requiring rapid and reliable quantification of leaf area.
LAMINA is suitable for use in a diverse range of species
Current limitations and future development
We have shown that LAMINA is able to accurately extract and quantify leaf area from scanned images of a diverse range of plant species. However, there are limitations to the use of the provided dimension measurements currently provided by LAMINA, and these limitations represent the most immediate targets for future development and expansion of LAMINA.
Although LAMINA is able to quantify leaf area of the A. annua leaves accurately, there are currently limitations to the use of additional measurements returned by LAMINA, with these being true in a range of leaf forms. Examining the A. annua leaves shown in Figure 3 shows that LAMINA currently returns blade dimensions using only straight lines, which is clearly far from ideal in these leaves. It is also clear that serration and cavity analysis will not return meaningful values from these leaves. There is therefore clear caution and consideration required by end users when making use of values returned by LAMINA. In the case of A. annua, we would suggest that calculated area is certainly reliable and that circularity may also be a useful indicator of how that leaf area is distributed. The use of any other returned values would require careful consideration by the end user.
As is the case for the A. annua leaves, many leaves do not have perfectly straight or symmetrical shapes and as such, the central line deviates from a straight line. Currently LAMINA only returns a measure of the maximum (or user set) straight line distance between the leaf base and tip. The inclusion of a tool to additionally allow manual placement of a non-straight line tracing the centre line (often the central vein) of a leaf is an obvious first target for extension of the current measures provided by LAMINA. At present, the difference between the values of the returned 25 % and 75 % vertical lengths can be used to indicate leaf asymmetry, which will often reflect the degree of leaf curvature and therefore the likely inaccuracy of the returned straight line centre measurement.
The A. annua leaves and the included example A. thaliana leaves shown in Figure 3 identify another important issue to consider when sampling leaves to be analysed using LAMINA – that of petioles: If petioles are sampled as well as actual leaf area, LAMINA will include the petiole as part of the leaf and this will affect generated measurements. In many species, removal of petioles is simple as there is a clearly identifiable boundary between leaf and petiole. If petioles are being removed, it is essential that this is done accurately as any remaining petiole will lead to the mis-identification of a serration either side of the remaining petiole. Figure 3G represents a more complex example, but one that is typical for many A. thaliana plants, where there is no clear boundary between the leaf lamina and the petiole. In such cases, it is often very hard to define where leaf becomes petiole and the sampling strategy must take this into consideration: LAMINA will include any scanned leaf area when generating dimension measurements and the end user must therefore decide what they wish to be included at the point of sampling (or by later manipulation of the generated scanned images to remove e.g. remaining, unwanted petiole area). Although it is not inconceivable that an algorithm could be developed to differentiate between leaf and petiole, this is certainly far from a trivial task, especially if such an algorithm should be generally applicable across species.
The example A. thaliana leaves in Figure 3G highlight another point that users must be aware of: currently, the software will not distinguish between wounding at the leaf boundary and serrations. In the cases shown in Figure 3, boundary damage most likely resulted from flattening the leaves at the point of image collection, as A. thaliana leaves are frequently curved and can not be flattened without tearing the leaf lamina. The use of median rather than mean serration values will limit the influence of such outliers but if serration quantification is being used, users should visually screen through the cropped images produced by LAMINA to identify problem leaves. It is possible that the algorithm for detecting leaf serrations could be extended to differentiate between boundary wounding or grazing herbivore damage (that typically extends from the boundary edge into the leaf) and actual serrations. However, as with other similar problems such as petiole detection, this will not be simple if the algorithm is to be applicable across a wide range of species (for example, many species contain serrations in combination with deeper, more infrequent, leaf lobes). Such algorithmic development would require extensive testing and confirmation across a broad range of species that have been exposed to a range of herbivore damage and wounding. There are also a number of potential extensions to LAMINA that we feel would have broad appeal to leaf researchers, including colour quantification (for example to track senescence), detection of necrotic lesions or flecks, measurement of leaf rust urediospore number and dimensions, and quantification of veinal pattern. As LAMINA has been released as an open source project using the well-supported Java language, it represents an ideal framework for the future integration of such extensions by the community and we hope that the instigation of such an open source project can serve as a means of concentrating development of a powerful phenotyping tool, as has been the case for the analysis of microscopy images since the initial release of ImageJ .
We have developed a new software tool for the automated or semi-automated analysis of leaf morphological traits and have shown that the method is able to extract biologically meaningful data from a range of species with contrasting leaf shapes. The developed software performs equally well as existing software while also providing an extended range of measures of leaf size and shape indicators. We show that the software performs as well as commonly used leaf area meters, even when measuring highly complex leaf forms. Application of this software tool will significantly aid the rapid screening of large-scale collections of genotypes for forward or reverse genetics as well as equally serving plant breeders. This is the first open source tool available for the quantification of leaf serration.
Availability and requirements
Project name: LAMINA: Leaf shApe deterMINAtion
Project home page: http://sourceforge.net/projects/lamina
Operating system(s): Platform independent
Programming language: Java
Other requirements: Java 1.4.x or higher. LAMINA uses the Java Advanced Imaging (JAI) package http://java.sun.com/javase/technologies/desktop/media/jai/ to support common image file formats, which is bundled with the installation and hence no additional installation should be required.
License: GNU GPL2
Leaf shApe deterMINAtion
Quantitative Trait Loci
Principal Component Analysis.
Work on A. annua carried out by VS and AMR is part of the CNAP A. annua Research Project, funded by the Bill and Melinda Gates Foundation.
- Fleming AJ: The control of leaf development. New Phytologist. 2005, 166: 9-20. 10.1111/j.1469-8137.2004.01292.x.PubMedView ArticleGoogle Scholar
- Pien S, Wyrzykowska J, Mason S, Smart C, Fleming A: Local expression of expansin induces the entire process of leaf development and modifies leaf shape. Proceedings of the National Academy of Sciences, USA. 2001, 98: 11812-11817. 10.1073/pnas.191380498.View ArticleGoogle Scholar
- Tsuge T, Tsukaya H, Uchimiya H: Two independent and polarized processes of cell elongation regulate leaf blade expansion in Arabidopsis thaliana (L.) Heynh. Development. 1996, 122: 1589-1600.PubMedGoogle Scholar
- Lee Y, Kim GT, Kim IJ, Park J, Kwak SS, Choi G, Chung WI: LONGIFOLIA1 and LONGIFOLIA2, two homologous genes, regulate longitudinal cell elongation in Arabidopsis. Development. 2006, 133: 4305-4314. 10.1242/dev.02604.PubMedView ArticleGoogle Scholar
- Tsukaya H, Tsuge T, Uchimiya H: The cotyledon: A superior system for studies of leaf development. Planta. 1994, 195: 309-312. 10.1007/BF00199692.View ArticleGoogle Scholar
- Rae AM, Ferris R, Tallis MJ, Taylor G: Elucidating genomic regions determining enhanced leaf growth and delayed senescence in elevated CO2. Plant Cell and Environment. 2006, 29: 1730-1741. 10.1111/j.1365-3040.2006.01545.x.View ArticleGoogle Scholar
- Street NR, Skogstrom O, Sjodin A, Tucker J, Rodriguez-Acosta M, Nilsson P, Jansson S, Taylor G: The genetics and genomics of the drought response in Populus. Plant Journal. 2006, 48: 321-341. 10.1111/j.1365-313X.2006.02864.x.PubMedView ArticleGoogle Scholar
- Juenger T, Perez-Perez J, Bernal S, Micol J: Quantitative trait loci mapping of floral and leaf morphology traits in Arabidopsis thaliana: evidence for modular genetic architecture. Evolution and Development. 2005, 7: 259-271. 10.1111/j.1525-142X.2005.05028.x.PubMedView ArticleGoogle Scholar
- Wullschleger S, Yin T, Difazio S, Tschaplinski T, Gunter L, Davis M, Tuskan G: Phenotypic variation in growth and biomass distribution for two advanced-generation pedigrees of hybrid poplar. Canadian Journal of Forest Research. 2005, 35: 1779-1789. 10.1139/x05-101.View ArticleGoogle Scholar
- El-Lithy M, Clerkx E, Ruys G, Koornneef M, Vreugdenhil D: Quantitative Trait Locus Analysis of Growth-Related Traits in a New Arabidopsis Recombinant Inbred Population. Plant Physiology. 2004, 135: 444-458. 10.1104/pp.103.036822.PubMedPubMed CentralView ArticleGoogle Scholar
- Frary A, Fritz L, Tanksley S: A comparative study of the genetic bases of natural variation in tomato leaf, sepal, and petal morphology. Theoretical and Applied Genetics. 2004, 109: 523-533. 10.1007/s00122-004-1669-x.PubMedView ArticleGoogle Scholar
- Perez-Perez J, Serrano-Cartagena J, Micol J: Genetic Analysis of Natural Variations in the Architecture of Arabidopsis thaliana Vegetative Leaves. Genetics. 2002, 162: 893-915.PubMedPubMed CentralGoogle Scholar
- Jiang C, Wright RJ, Woo SS, Delmonte TA, Paterson AH: QTL analysis of leaf morphology in tetraploid Gossypium (cotton). Theoretical and Applied Genetics. 2000, 100: 409-418. 10.1007/s001220050054.View ArticleGoogle Scholar
- Wu R, Stettler R: Quantitative genetics of growth and development in Populus. III. Phenotypic plasticity of crown structure and function. Heredity. 1998, 81: 299-310. 10.1038/sj.hdy.6883970. [http://www.nature.com/hdy/journal/v81/n3/abs/6883970a.html]View ArticleGoogle Scholar
- Wu R, Bradshaw HD, Stettler RF: Molecular genetics of growth and development in Populus (Salicaceae) .5. Mapping quantitative trait loci affecting leaf variation. American Journal of Botany. 1997, 84: 143-153. 10.2307/2446076.PubMedView ArticleGoogle Scholar
- Wilf P, Wing SL, Greenwood DR, Greenwood CL: Using fossil leaves as paleoprecipitation indicators; an Eocene example. Geology. 1998, 26: 203-206. 10.1130/0091-7613(1998)026<0203:UFLAPI>2.3.CO;2.View ArticleGoogle Scholar
- Krieger JD, Guralnick RP, Smith DM: Generating empirically determined, continuous measures of leaf shape for paleoclimate reconstruction. PALAIOS. 2007, 22: 212-219. 10.2110/palo.2005.p05-079r.View ArticleGoogle Scholar
- Huff PM, Wilf P, Azumah EJ: Digital Future for Paleoclimate Estimation from Fossil Leaves? Preliminary Results. PALAIOS. 2003, 18: 266-274. 10.1669/0883-1351(2003)018<0266:DFFPEF>2.0.CO;2.View ArticleGoogle Scholar
- Royer DL, Wilf P, Janesko DA, Kowalski EA, Dilcher DL: Correlations of climate and plant ecology to leaf size and shape: potential proxies for the fossil record. Am J Bot. 2005, 92: 1141-1151. 10.3732/ajb.92.7.1141.PubMedView ArticleGoogle Scholar
- ImageJ. [http://rsb.info.nih.gov/ij/]
- Igathinathane C, Prakash VSS, Padma U, Babu RG, Womac AR: Interactive computer software development for leaf area measurement. Computers and Electronics in Agriculture. 2006, 51: 1-16. 10.1016/j.compag.2005.10.003.View ArticleGoogle Scholar
- Bakr E: A new software for measuring leaf area, and area damaged by Tetranychus urticae Koch. Journal of Applied Entomology. 2005, 129: 173-175. 10.1111/j.1439-0418.2005.00948.x.View ArticleGoogle Scholar
- Weight , Caroline , Parnham , Daniel , Waites , Richard : LeafAnalyser: a computational method for rapid and large-scale analyses of leaf shape variation. The Plant Journal. 2008, 53: 578-586. 10.1111/j.1365-313X.2007.03330.x.PubMedView ArticleGoogle Scholar
- Luquez V, Hall D, Albrectsen B, Karlsson J, Ingvarsson P, Jansson S: Natural phenological variation in aspen (Populus tremula): the SwAsp collection. Tree Genetics & Genomes. 2007Google Scholar
- Ingvarsson PKK, Garcia V, Luquez V, Hall D, Jansson S: Nucleotide polymorphism and phenotypic associations within and around the phytochrome B2 locus in European aspen (Populus tremula, Salicaceae). Genetics. 2008Google Scholar
- Gonzalez RC, Woods RE: Digital image processing. 2002, Prentice Hall, 2Google Scholar
- Ihaka R, Gentleman R: R: a language for data analysis and graphics. Journal of Computational and Graphical Statistics. 1996, 5: 299-314. 10.2307/1390807.Google Scholar