Abstract

Despite the ecological importance of eelgrass, nowadays anthropogenic influences have produced deleterious effects in many meadows worldwide. Transplantation plots are commonly used as a feasible remediation scheme. The characterization of eelgrass biomass and its dynamics is an important input for the assessment of the overall status of both natural and transplanted populations. Particularly, in restoration plots it is desirable to obtain nondestructive assessments of these variables. Allometric models allow the expression of above ground biomass and productivity of eelgrass in terms of leaf area, which provides cost effective and nondestructive assessments. Leaf area in eelgrass can be conveniently obtained by the product of associated length and width. Although these variables can be directly measured on most sampled leaves, digital image methods could be adapted in order to simplify measurements. Nonetheless, since width to length ratios in eelgrass leaves could be even negligible, noise induced by leaf humidity content could produce misidentification of pixels along the peripheral contour of leaves images. In this paper, we present a procedure aimed to produce consistent estimations of eelgrass leaf area in the presence of the aforementioned noise effects. Our results show that digital image procedures can provide reliable, nondestructive estimations of eelgrass leaf area.

1. Introduction

Zostera marina also known as eelgrass is a relevant seagrass species, which supplies significant amounts of organic materials to food webs in shallow coastal environments and provides habitat (in bays, lagoons, or estuaries) for many fishes and their larvae [1]. Eelgrass beds can also help remediate contaminated sediments [2], filter and retain nutrients from the water column [3], help in the stabilization of sediments [4], and reduce erosion forces by stumping wave energy, thus promoting the stabilization of adjacent shorelines [5]. However, the permanence of eelgrass beds—as well as those formed by other seagrass species—is currently being threatened by anthropogenic influences to such an extent that special conservation efforts are needed [6]. This requires the development of accurate and cost-effective procedures aimed at obtaining scientific knowledge about the pertinent growth dynamics. This is particularly relevant in assessments of restoration projects, where the use of noninvasive data gathering techniques turns out to be of fundamental importance.

The characterization of eelgrass biomass and its dynamics is an important input for the assessment of the overall status of both natural and transplanted eelgrass populations. In eelgrass the basic unit for studying biomass and its production is the shoot, which includes sheaths, leaves, rhizomes, and roots. Biomass consists of an aboveground component formed by sheaths and leaves and a belowground constituent formed by rhizomes and roots. Root emergence occurs at leaf scars, also known as rhizome nodes. The production of leaves and rhizome nodes is connected such that each leaf produced is linked to a rhizome node. Hence, the overall production of shoots can be estimated by measuring the production of leaves [7]; this makes us know the growth rate of leaves fundamental to the assessment of eelgrass populations [8]. Moreover, estimations of leaf biomass and leaf-growth rates are keys to assessing the reestablishment of ecological functioning in restored areas. Nevertheless, traditional methods for the estimation of eelgrass leaf biomass and the related leaf growth rates are destructive and time consuming. Even though these procedures do not damage natural seagrass populations, they could produce undesirable effects on transplant experiments. Favorably, the conspicuous growth form of eelgrass makes it possible to introduce proxies that allow assessments while avoiding invasive interference. Moreover, estimations of leaf biomass or productivity in eelgrass can be efficiently obtained using allometric alternatives, which state these variables in terms of leaf length or area [9, 10]. But even though leaf architecture in eelgrass makes length a consistent descriptor of area, allometric models that express leaf biomass in terms of linked area perform relatively better than those involving leaf length as an independent variable. Therefore, for consistent allometric estimations of leaf biomass or productivity of eelgrass it is convenient to produce reliable estimations of leaf area. The observed ribbon-like appearance of the leaves in Zostera marina is a feature that permits obtaining direct and fairly accurate estimations of blade length and width . These variables provide convenient estimations of the corresponding blade area , through the leaf length times width proxy [1]. If we used the symbol as a subscript to represent observed values for the above named variables, then estimations of leaf area obtained through this proxy are given by which combined with allometric methods could simplify assessments of eelgrass leaf biomass and productivity [9, 10].

Digital image processing techniques were initially aimed to calculate the area of leaves for terrestrial plants [1113]. These methods provide simplified estimations of biologically relevant variables. For example, Patil and Bodh [14] used area of sugarcane leaves for plant growth monitoring to analyze manure scarcity and environmental stress and to assess disease severity. Lü et al. [15] used leaf area measurement to assess long-term influences on yield and because it is a fundamental index in crop growth and nurturing practice. Although, leaf area in eelgrass can be conveniently obtained by means of (1), and both and can be directly measured on most sampled leaves, methods based on digital image processing could be adapted in order to simplify these tasks. Moreover, eelgrass leaf area can be directly estimated from digital imagery by using the Monte Carlo method [1]: if we let denote these estimations, then they are obtained through where LPN is the number of points placed inside the considered leaf area and stands for the number of points contained in a unit area.

Besides, the Monte Carlo method eelgrass leaf area could be also obtained from digital images by using the length times width proxy of (1). Indeed, if and , respectively, denote leaf length and width obtained from the associated digital image, then these variables can be estimated through where and are, respectively, obtained by counting the number of points contained over the length and width dimensions of the leaf, and is the number of points contained in the appropriate distance measurement unit. Therefore, denoting by means of the associated leaf area, we will correspondingly have Nevertheless, when using either (2) or (5) to produce estimations of leaf area we must be aware that some Zostera marina leaves could be very long or present curvatures, among irregularities caused by environmental factors like grazing or drag forces. The influence of these factors could affect image quality, which could produce biased estimations for leaf area. These effects have been partially addressed by Ramfos et al. [16], who proposed a method based on image processing techniques for measurements of a Zostera marina leaf by taking into account the effects of curvature on accuracy. Yet another important factor which we address here concerns the effects that the humidity contents of a leaf can originate in image processing. In fact, once leaves are removed from a shoot they begin to lose water and degrade. Hence, if leaves cannot be processed immediately after being collected, it is important to keep them in a manner that reduces changes in shape [17]. Therefore, an efficient digitalizing of a Zostera marina blade requires maintenance of an optimal humidity for increased image fidelity. On the other hand, humidity contents in a leaf can induce noise to an image by adding extraneous information, which usually manifests by pixel value misidentification.

Data published by Echavarria-Heras et al. [1], taken over a comprehensive sampling experiment, show that measured maximum width for a Zostera marina leaf is 6 mm. Surely, other authors report a variation range from 1.5 to 12 mm for this estimation [18]. A wide variation range in width in conjunction with noise due to humidity content can increase uncertainty in blade width measurements obtained from digitalized leaf images. This makes it necessary to devise a way that allows discriminating the concomitant error spreading over leaf area assessments. So far, an approach that integrates among others techniques, one aimed to handle noise effects induced by the humidity contents on a Zostera marina leaf, has not been produced. In this study, we conceived a method which using criteria based on statistical analysis techniques reduces the effects that noise linked to the humidity contents of a Zostera marina leaf produces on the accuracy of associated area estimations obtained from a digital image.

2. Conceptual Framework for Image Processing

Our arrangement depends in a fundamental way on the concept of the peripheral or bordering contour of a bidimensional enclosure or domain. Several definitions of peripheral contour exist, being each one appropriate for different settings. Our interpretation is similar to the perimeter definition of a regular pattern in geometry. More accurately, a peripheral contour in the present settings will be defined as the sequence of boundary pixels of a digitalized eelgrass leaf. Moreover, for a reasonable identification of the area of the pertinent blade it is imperative that in the extents of the corresponding image minimal changes of color levels occur, even though around its outer contour abrupt changes of a color levels could be shown. The effect projected by humidity adds to leaf area pixels placed between the pixels captured by the digital image of the leaf itself and others belonging to its background. Hence, a reliable imbedding of the area of a leaf into an image requires the unambiguous identification of the pixels on its surrounding contour.

Our design is aimed to the aforesaid identification in the presence of noise due to the humidity content in the leaf. For the incumbent characterization, our system uses a quantitative setup developed on the basis of the concepts of adjacency, vicinity, connectivity, and tolerance of similarity between pixels. We briefly describe these notions in what follows.

Two pixels are adjacent if and only if they share one of their borders, or at least one of their corners. Two pixels are neighbors if they fulfill the definition of adjacency. Formally, the vicinity of the point is defined through Without loss of generality, we explain the notion of tolerance of similarity by referring to the RGB format description of a color. This allows quantifying tonality in terms of the intensities of the constituting primary colors: red, green, and blue. To indicate at which amount each one of these colors is mixed to produce a given tonality a value is assigned to each prime color; for example, the value 0 means that a given basic color does not appear in the mix, but if a chief color component is nonvanishing it means that it contributes to the mix in a given intensity. We set which identifies the colors number to be used through the whole image processing task; for an RGB color space we have . Usually, the intensity of each of the primary colors appearing in a mix is measured on a scale ranging from 0 to . The set of all color intensities can be represented in the form of a cube in the cartesian coordinate system, where each color is a point on the surface or in its interior. Given points and in an RGB color space, we will define the distance between them through Moreover, given a point in an RGB color space, a second one with the greatest similarity to is the one placed at the smallest distance . Furthermore, let be a color tonality range, with being the number of different colors included. Then, we must have and we will say that two pixels and are similar to a tolerance limit if the inequality is satisfied. In what follows the range will be simply called “tolerance of similarity” and the upper bound can be interpreted as the maximum distance that two points located within the extent of an object can attain in an RGB color space in order to be considered similar. Connectivity between pixels is used to identify the limits in objects and regions in an image. We will say that two pixels and are connected with tolerance of similarity if they fulfill the definition of adjacency and also if inequality (8) holds.

3. The Image Selection Method

The procedure to obtain efficient estimations of , , and requires two stages. On a first one we create a digital image for each one of the collected leaves. Then, we set and continue by choosing an interval of tolerance of similarity with ; we use this to obtain the peripheral contour of each one of the available leaf images and from them the linked , , and values. Different intervals will produce different estimations for , , and , and consequently we must rely on a criterion for the selection of the range that produces the most accurate estimations. To carry out this task in a second stage of the method we arrange leaf length data into groups of leaves whose size differences are bounded by a preferred tolerance and use that arrangement to obtain related statistics and that are used to implement what we call the selection index. In what follows we describe pseudo-codes for the above referred stages. Detailed formulae are presented in the appendices. Tables 1, 2, 3, 4, and 5 summarize the involved notation.

3.1. The Procedure to Obtain , , and Assessments

(a.1)Choose a color format and set .(a.2)Load the leaf image.(a.3)Enter an interval of tolerance of similarity ; .(a.4)Select a starting point inside the loaded leaf image.(a.5)Find the contour of the leaf image through (6), (7), and (8) (these equations identify all adjacent pixels falling within the selected interval of tolerance of similarity .(a.6)Obtain , , and by using (3), (4), and (5), respectively.(a.7)Obtain by using (2).(a.8)Record and the associated , , , and estimations.(a.9)Repeat steps 2–8 for each one of the available leaf images.(a.10)Change the interval and jump to step (a.3).

Different intervals will produce through the above procedure different estimations for , , , and . We now outline a procedure for the selection of the image that produces the most accurate estimations or for the observed leaves area . This requires the identification of the interval of tolerance of similarity that yields the smallest values of the selection index defined by (9) below.

3.2. The Method for the Selection of an Optimal Interval

(b.1)For the entered interval, use to calculate (this value gives the proportion of leaves for which produces consistent estimations of ).(b.2)For the entered interval, use to calculate (this value yields the proportion of leaves for which overestimates observed leaf area ).(b.3)For the entered interval, calculate the value of the image selection index according to (b.4)Record both and .(b.5)Change the interval and repeat steps (b.1) to (b.3) until all the intervals generated in Section 3.1 are exhausted.(b.6)Choose the interval that produces the smallest value of for image processing and leaf area estimations.

The above selection index criterion can be adapted for Monte Carlo method estimations of leaf area. It becomes where and are, respectively, given by and in Appendix E and are equivalent to and correspondingly.

4. Results

4.1. Leaf Data Grouping

The present data set was obtained by randomly sampling 5 shoots biweekly from January through December 2009 in a Zostera marina field at Punta Banda estuary, a shallow coastal lagoon located near Ensenada, Baja California, Mexico (31° 43–46 N and 116° 37–40 W). For each sampled leaf, a millimeter ruler was used to obtain to the nearest 1/10 mm taken as the distance from the top of the sheath to the leaf tip. Meanwhile, was measured at a point halfway between the top of the sheath and the tip. Observed leaf area estimations were calculated by means of (1).

We obtained . For data grouping we selected , so we acquired and for the interval we formed a partition of disjoint intervals of the form , with . Hence, as described in the appendices for each value of the index , we formed a group containing leaves with sizes varying in the interval (Table 6). Longer and older leaves displayed darker tonalities than younger and shorter ones. Moreover, leaves with lengths varying on a given partition interval displayed a similar color distribution. For some of the partition intervals there was at most one leaf with length placed in the linked variation range. Therefore, these groups are not taken into account because they do not provide information for the statistical analysis (see bold in Table 6).

4.2. Image Selection Procedure

For each one of the leaves in the collection of groups , we applied the procedure described in the pseudo-code (a) aimed to detect the points on the associated peripheral contour and to get the concomitant , , , and estimations. For that purpose we used a variety of equivalence of tones, which permitted an unambiguous framing of the extent of the leaf. A RGB 256 color format was used for all leaves images. Hence, we set . Therefore, different tolerances of similarity ranges, with , were used. Moreover, the procedure was automatically applied up to 256 times on each individual leaf image. For every tolerance of similarity interval we selected a starting point inside a chosen leaf image and we identified all adjacent pixels falling within the named similarity range. This identified the peripheral contour of the leaf image so the linked , , and assessments as well as leaf area estimations were acquired (see B.3 in Appendix B).

For each group of leaves determined by the partition we calculated deviation values and and their averages and taken over . We also calculated the associated standard deviation values and (see Appendix D) and the proportions, , . , ,   , , and (see Appendix E); calculated values are presented in Tables 7 and 8. Values of the image selection index were obtained and compared. For easy of presentation we focus on the results obtained for , , and which include the smallest obtained value for (Table 8). Figures 1, 2, and 3 display comparisons of the averages and of observed leaf lengths and widths in groups versus estimations and from images acquired using these tolerance of similarity intervals.

For , and values were greater than those calculated for . Moreover, and values obtained using were greater than those obtained using but smaller than those obtained for . Nevertheless, values produced for were smaller than those linked to , but generally values increased implying a greater concentration of a particular tonality within the range of colors forming the color spectrum of the image. That is, broadens or reduces the collection of colors which can be taken into account for pixel identification within the extent of the image. Whenever stands for a greater amount of colors relative to the assortment defining the image, the pixel selection procedure will lead to subjective identification. This influences and values in a direct way. Indeed for , was negative (Table 7) which means that most values were greater than values. Nevertheless, this outcome is limited by the maximum value that can attain and also by image size. For , the value of remained unchanged.

By using the proportion values , , and , we assessed that for leaf area was overestimated by the method. And although for the method slightly overestimated leaf area, for intervals beyond leaf area overestimation always increased. Moreover, in Table 9 we provide root-mean-square deviation (RMSD) values for comparisons of directly versus image obtained estimations when using intervals. We can assess that a interval produced the highest RMSD values for the comparison of directly versus image obtained estimations of leaf width, length or area. This is consistent with the smallest value obtained for the selection index. Therefore, we suggest that a similarity interval must be set to process the present Zostera marina leaf image set through the method proposed.

5. Discussion

Deleterious effects derived from anthropogenic influences are currently increasing worldwide, threatening the health of many eelgrass meadows [19]. Remediation efforts have included transplant projects [20] and the valuation of their status, requires the estimations of key variables such as standing stock or productivity. Although shoot removal for the measurement of these variables does not damage natural seagrass populations, these procedures could produce undesirable effects on transplants. Therefore, when standing stock or productivity assessment are performed over the early stages of an eelgrass restoration experiment, data gathering approaches that avoid disruptive interference are essential. Allometric methods can provide convenient proxies, which reduce leaf biomass and growth assessments to simple blade length or area measurements. What is more, if the estimation of these leaf attributes can be done without removing the blades allometric approaches could furnish truly non-destructive assessments. Modern electronic scanning technologies could be used under water in order to produce reliable images of marine plants leaves, which guarantee non-destructive sampling of leaves length, width or area. However, insitu scanning of eelgrass leaves could add extraneous information mainly due to the inherent humidity content or to materials attached to blades like mud related particles. Hence, for estimating leaf attributes such as length or area, we must take into account that image noise effects could render biased assessments. Moreover, if we strive to use (6), the digital image counterpart of (1), we must take into account that a Zostera marina leaf does not show a perfect rectangular shape. We should also notice that since width attains small values, noise produced by humidity could increase the width of the digital leaf in a noticeable way which could certainly heighten uncertainty in leaf area estimations produced through (5). What is more, even when leaf area is estimated from images using Monte Carlo methods, noise effects could produce ambiguity in peripheral contour identification rendering biased results.

In the other hand, we need to be aware that since the power functions involved in allometric approaches are highly sensitive to parameter uncertainty we require consistent estimations of blade length or area [10]. Therefore, when we estimate eelgrass leaf area, using digital imagery in combination with allometric proxies, we must be aware that image noise could certainly reduce the accuracy of estimates. This makes it necessary to rely on efficient image selection methodologies for uncertainty reduction. Our results show that the present methods produce reliable results. This conclusion is mainly substantiated by the obtained values for the RMSD. We used these statistics to determine consistency between directly obtained measurements and image estimated assessments. Table 9 shows that the highest RMSD values corresponded with the smallest values obtained for the selection index. Moreover, the interval selected using the criteria produced also the highest RMSD values for the comparison of observed values and those obtained by means of Monte Carlo method. This justifies our claim that the proposed procedure abridged by the selection index can be expected to produce consistent estimations of the leaf attributes necessary for allometric estimation of relevant variables required to assess the status of an eelgrass population. Moreover, the presented procedure could be straightforwardly applied to other eelgrass populations or seagrass species that exhibit similar leaf architectures making leaf length times width a reliable proxy for the pertinent area.

Appendices

We now explain how to identify the interval of tolerance of similarity that yields accurate estimations for the observed leaf area values. The task, is achieved through statistical methodologies, which requires the completion of the following steps.

A. Grouping Leaf Data: and

Identify the maximum observed leaf length .Chose a positive integer and define a partition of the interval with norm .Form the collection of disjoint intervals of the form with . This collection is denoted thought .For each value of the index identify the group of leaves whose lengths are contained in . Notice that holds the leaves whose size differences are bounded by .For each value of the index obtain and record standing for the number of leaves in the group .For each value of the index introduce an index such that and label as , and respectively, the straight length, width and area of the th leaf in . The character denotes the associated estimations of leaf area obtained by means of (1).Form and record the collection of all groups of leaves . This collection is denoted by means of .Obtain the average length for each group of leaves . That is, calculate and record Obtain the average width for each group of leaves . That is, calculate and record Obtain the average area for each group of leaves . That is, calculate and record,

B. Obtaining Length, Width and Area from the Image of each Leaf

For processing all digital images, we chose a specified color format with a number of colors.For processing the digital images of all collected leaves, we choose different intervals of tolerance of similarity , with the upper bound satisfying .For a picked interval, for use the algorithm (a) described in the method section to obtain and , which respectively denote the length and width of the image of the th leaf in . Also, obtain and , which respectively stand for leaf area obtained from the image and calculated by means of (1) and (2) respectively. Record these values.interval obtain and record the concomitant averages , , and (cf. (A.1) through (A.3)).

C. Obtain Estimation Errors between the Observed and Image Obtained Values in Steps A and B

For the picked interval and for , calculate the leaf length approximation errors through For the picked interval and for , calculate the individual leaf width approximation errors through For the picked interval and for , calculate the leaf area approximation errors through For the picked interval and for , calculate the leaf area approximation errors linked to the Monte Carlo method through

D. Obtain the Average Deviations Produced by the Individual Estimation Errors

For the picked interval, obtain the average leaf length deviations , which are calculated by averaging the values. That is, We notice that and also that negative values of imply that in a lot most image assessments overestimate observed values.For the picked interval, obtain the average leaf width deviations , which are calculated by averaging the values. That is, notice that since , negative values of imply that in a group most image assessments overestimate observed values.For the picked interval, calculate the average leaf area deviations by averaging the values. That is, again since negative values of imply that in a group most image assessments overestimate observed values.For the picked interval, calculate , the average value of deviations taken over . Calculate also the associated standard deviation .For the picked interval, calculate , the average value of deviation taken over . Calculate also the associated standard deviation .

E. Criteria for Selecting the Interval That Produces the Highest Correspondence Level between Image Obtained Measurements and Those Obtained Directly from Collected Leaves

For a given range of similarity values , identify the leaves satisfying the conditions Calculate the proportion of leaves in that comply with the condition (E.3) through Calculate the proportion of leaves in that comply with the condition (E.4) through Calculate the proportion of leaves in that do not comply with the condition (E.3) through, Calculate the proportion of leaves in that do not comply with the condition (E.4) through, Obtain the concomitant proportions of leaves in that provide consistent leaf area estimations by the proxy of (1) Calculate the proportion of leaves in for which image estimated blade length and width measurements overestimate leaf area calculated through (1) The proportions and equivalent to and respectively but linked to leaf area estimation by Monte Carlo method (cf. (2)) that is, Conditions (E.1) and (E.2) grant bounded estimation errors, for and respectively. Moreover, the groups of leaves that also satisfy conditions (E.3) and (E.4) can be identified as those groups for which image and estimations are closer to directly obtained and measurements. Therefore, groups in , which do not comply with conditions (E.1)–(E.4), denote the set of leaves with estimation errors out of range for a reliable estimation.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

We thank Jose Maria Dominguez and Francisco Ponce for the art work.