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 l and width h. These variables provide convenient estimations of the corresponding blade area a, through the leaf length times width proxy [1]. If we used the symbol o as a subscript to represent observed values for the above named variables, then estimations of leaf area obtained through this proxy are given by
(1)ao=lo·ho
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 [11–13]. 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 lo and ho 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 amc denote these estimations, then they are obtained through
(2)amc=LPNUPN2,
where LPN is the number of points placed inside the considered leaf area and UPN2 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 ld and hd, respectively, denote leaf length and width obtained from the associated digital image, then these variables can be estimated through
(3)ld=nplunp,(4)hd=nphunp,
where npl and nph are, respectively, obtained by counting the number of points contained over the length and width dimensions of the leaf, and unp is the number of points contained in the appropriate distance measurement unit. Therefore, denoting by means of ad the associated leaf area, we will correspondingly have
(5)ad=ld·hd.
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 Vp(x,y) of the point P(x,y) is defined through
(6)Vp(x,y)=(x+1,y),(x-1,y),(x,y+1),(x,y-1),(x+1,y+1),(x+1,y-1),(x-1,y+1),(x-1,y-1).
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 Cmax which identifies the colors number to be used through the whole image processing task; for an RGB color space we have Cmax=256. Usually, the intensity of each of the primary colors appearing in a mix is measured on a scale ranging from 0 to Cmax-1. 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 P=(p1,p2,…,pn) and Q=(q1,q2,…,qn) in an RGB color space, we will define the distance dE(P,Q) between them through
(7)dE(P,Q)=∑i=1n(pn-qn)2.
Moreover, given a point P in an RGB color space, a second one Q with the greatest similarity to P is the one placed at the smallest distance dE(P,Q). Furthermore, let ST(x)=[0,x] be a color tonality range, with x being the number of different colors included. Then, we must have 1≤x≤Cmax-1 and we will say that two pixels P and Q are similar to a tolerance limit ST(x) if the inequality
(8)dE(P,Q)≤x
is satisfied. In what follows the range ST(x) will be simply called “tolerance of similarity” and the upper bound x 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 P and Q are connected with tolerance of similarity ST(x) 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 ld, hd, and amc requires two stages. On a first one we create a digital image for each one of the collected leaves. Then, we set Cmax and continue by choosing an interval of tolerance of similarity ST(x) with 0≤x≤Cmax-1; we use this to obtain the peripheral contour of each one of the available leaf images and from them the linked ld, hd, and amc values. Different intervals ST(x) will produce different estimations for ld, hd, and amc, and consequently we must rely on a criterion for the selection of the ST(x) 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 q and use that arrangement to obtain related statistics βa and λa that are used to implement what we call the ISx 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.
Different symbols used in the digital image processing method.
Symbol
Description
l
Leaf length
h
Leaf width
Cmax
Number of colors in a format of a digital image
ST(x)=[0,x]
Interval of Tolerance of Similarity
lmax
Maximum observed leaf length
q=lmaxn
Norm of the partition for the interval [0,lmax]
Ik
Partition interval of the form [q(k-1),qk) for 0≤k≤n
P0lmax=⋃1n(Ik)
Collection of n intervals Ik that cover [0,lmax]
Gk(l)
Group of leaves whose lenghts (l) lie in Ik
nk
Number of leaves in the group Gk(l)
CG=⋃1nGk(l)
Collection of all groups Gk(l) of leaves
Symbols for observed, digitally obtained variables and related averages.
Description
Observed data
Digital data
Monte Carlo data
Leaf length (l)
lo
ld
—
Leaf width (h)
ho
hd
—
Leaf area (a)
ao
ad
amc
Length of the jth leaf in group Gk(l)
lojk
ldjk
—
Width of the jth leaf in group Gk(l)
hojk
hdjk
—
Area of the jth leaf in group Gk(l)
aojk
adjk
amcjk
Average length of the leaves in group Gk(l)
l-ok
l-dk
—
Average width of the leaves in group Gk(l)
h-ok
h-dk
—
Average area of the leaves in group Gk(l)
a-ok
a-dk
a-mck
Approximation errors.
Symbols
Formal expression
Description
eljk
lojk-ldjk
Difference of observed and image obtained leaf lengths in group Gk(l).
ehjk
hojk-hdjk
Difference of observed and image obtained leaf widths in group Gk(l).
eajk
aojk-adjk
Difference of observed and image obtained leaf areas in group Gk(l).
emcjk
aojk-amcjk
Difference of observed and Monte Carlo estimated leaf areas in group Gk(l).
Estimation errors for observed and image obtained variables averages and standard deviations.
Symbol
Formal expression
Description
δlk
∑1nkeljknk
Leaf length average deviation in the group Gk(l).
δhk
∑1nkehjknk
Leaf width average deviation in the group Gk(l).
δak
∑1nkeajknk
Leaf area average deviation in the group Gk(l).
δ-l
∑1nδlkn
Leaf length average deviation in CG.
δ-h
∑1nδhkn
Leaf width average deviation in CG.
σδl
(∑1nk(eljk-δlk)2(nk-1))1/2
Standard deviation of eljk.
σδh
(∑1nk(ehjk-δhk)2(nk-1))1/2
Standard deviation of ehjk.
Auxiliary statistics λl, λh, θl, θh, λa, βa, λamc and βamc used to obtain the set of leaves with estimation errors in range for a reliable estimation.
Symbol
Description
Reference equation
λl
Proportion of leaves in CG for which δlk satisfies:δ-l-σδl≤δlk≤δ-l+σδl
(E.5)
λh
Proportion of leaves in CG for which δhk satisfies:δ-h-σδh≤δhk≤δ-h+σδh
(E.6)
θl
Proportion of leaves in CG for which δlk do not satisfies:δ-l-σδl≤δlk≤δ-l+σδl
(E.7)
θh
Proportion of leaves in CG for which δhk do not satisfies:δ-h-σδh≤δhk≤δ-h+σδh
(E.8)
λa
Proportion of leaves in CG for which δlk and δhk satisfies:δ-l-σδl≤δlk≤δ-l+σδl and δ-h-σδh≤δhk≤δ-h+σδh, and eajk≥0.
(E.9)
βa
Proportion of leaves in CG for which δlk and δhk do not satisfies:δ-l-σδl≤δlk≤δ-l+σδl and δ-h-σδh≤δhk≤δ-h+σδh
(E.10)
λamc
Proportion λamc equivalent to λa respectively but linked to leaf area estimation by Monte Carlo method (cf. (2)).
(E.11)
βamc
Proportion βamc equivalent to βa but linked to leaf area estimation by Monte Carlo method (cf. (2)).
(E.12)
3.1. The Procedure to Obtain ld, hd, and amc Assessments
Choose a color format and set Cmax.
Load the leaf image.
Enter an interval of tolerance of similarity ST(x); 1≤x≤Cmax-1.
Select a starting point inside the loaded leaf image.
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 ST(x)).
Obtain ld, hd, and ad by using (3), (4), and (5), respectively.
Obtain amc by using (2).
Record ST(x) and the associated ld, hd, ad, and amc estimations.
Repeat steps 2–8 for each one of the available leaf images.
Change the ST(x) interval and jump to step (a.3).
Different ST(x) intervals will produce through the above procedure different estimations for ld, hd, ad, and amc. We now outline a procedure for the selection of the image that produces the most accurate estimations ad or amc for the observed leaves area ao. This requires the identification of the interval of tolerance of similarity ST(x) that yields the smallest values of the selection index ISx defined by (9) below.
3.2. The Method for the Selection of an Optimal ST(x) Interval
For the entered ST(x) interval, use (E.9) to calculate λa (this value gives the proportion of leaves for which ad produces consistent estimations of ao).
For the entered ST(x) interval, use (E.10) to calculate βa (this value yields the proportion of leaves for which ad overestimates observed leaf area ao).
For the entered ST(x) interval, calculate the value of the image selection index ISx according to
(9)ISx=βaλa.
Record both ST(x) and ISx.
Change the ST(x) interval and repeat steps (b.1) to (b.3) until all the ST(x) intervals generated in Section 3.1 are exhausted.
Choose the ST(x) interval that produces the smallest value of ISx for image processing and leaf area ad estimations.
The above selection index ISx criterion can be adapted for Monte Carlo method estimations of leaf area. It becomes
(10)ISxmc=βamcλamc,
where λamc and βamc are, respectively, given by (E.11) and (E.12) in Appendix E and are equivalent to λa and βa correspondingly.
4. Results4.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 lo to the nearest 1/10 mm taken as the distance from the top of the sheath to the leaf tip. Meanwhile, ho was measured at a point halfway between the top of the sheath and the tip. Observed leaf area estimations ao were calculated by means of (1).
We obtained lmax=460mm. For data grouping we selected n=46, so we acquired q=10mm and for the interval [0,lmax] we formed a partition P0460 of disjoint intervals Ik of the form Ik={l∣q(k-1)≤l<qk}, with 1≤k≤46. Hence, as described in the appendices for each value of the index k, we formed a group Gk(l) containing leaves with sizes varying in the interval Ik (Table 6). Longer and older leaves displayed darker tonalities than younger and shorter ones. Moreover, leaves with lengths varying on a given partition interval Ik 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).
Numbers nk of whole leaves classified in groups Gk(l) formed by leaf sizes varying in corresponding length intervals Ik.
k
Ik
nk
k
Ik
nk
1
[0, 10)
10
24
[230, 240)
24
2
[10, 20)
43
25
[240, 250)
23
3
[20, 30)
38
26
[250, 260)
15
4
[30, 40)
38
27
[260, 270)
21
5
[40, 50)
32
28
[270, 280)
16
6
[50, 60)
37
29
[280, 290)
12
7
[60, 70)
43
30
[290, 300)
10
8
[70, 80)
32
31
[300, 310)
9
9
[80, 90)
34
32
[310, 320)
9
10
[90, 100)
38
33
[320, 330)
4
11
[100, 110)
28
34
[330, 340)
7
12
[110, 120)
40
35
[340, 350)
3
13
[120, 130)
28
36
[350, 360)
4
14
[130, 140)
29
37
[360, 370)
3
15
[140, 150)
19
38
[370, 380)
3
16
[150, 160)
27
39
[380, 390)
1
17
[160, 170)
19
40
[390, 400)
1
18
[170, 180)
14
41
[400, 410)
1
19
[180, 190)
17
42
[410, 420)
1
20
[190, 200)
21
43
[420, 430)
1
21
[200, 210)
19
44
[430, 440)
0
22
[210, 220)
20
45
[440, 450)
0
23
[220, 230)
14
46
[450, 460)
1
4.2. Image Selection Procedure
For each one of the leaves in the collection CG of groups Gk(l), 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 ld, hd, ad, and amc 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 Cmax=256. Therefore, different tolerances of similarity ranges, ST(x)=[0,x] with 1≤x≤255, were used. Moreover, the procedure was automatically applied up to 256 times on each individual leaf image. For every tolerance of similarity interval ST(x) 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 ldjk, hdjk, and adjk assessments as well as leaf area estimations amcjk were acquired (see B.3 in Appendix B).
For each group Gk(l) of leaves determined by the partition P0460 we calculated deviation values δhk and δlk and their averages δ-h and δ-l taken over CG. We also calculated the associated standard deviation values σδh and σδl (see Appendix D) and the proportions, λh, λl. θh, θl, λaβa, λamc, and βamc (see Appendix E); calculated values are presented in Tables 7 and 8. Values of the image selection index ISx were obtained and compared. For easy of presentation we focus on the results obtained for ST(68), ST(128), and ST(192) which include the smallest obtained value for ISx (Table 8). Figures 1, 2, and 3 display comparisons of the averages l-ok and h-ok of observed leaf lengths and widths in groups Gk(l) versus estimations l-dk and h-dk from images acquired using these tolerance of similarity intervals.
Direct comparison statistics for different ST(x) range values.
ST(x)
δ-h
σδh
δ-l
σδl
θl
θh
λl
λh
ST(68)
0.4493
0.2721
24.0157
23.6548
0.0161
0.1038
0.9839
0.8962
ST(128)
0.2599
0.2576
5.0342
13.7282
0.0049
0.0445
0.9951
0.9555
ST(192)
−0.1291
0.2496
3.8965
12.9700
0.0049
0.1669
0.9951
0.8331
Proportions of overestimation and underestimation of leaf area and selection index values for a given ST(x) range.
ST(x)
λa
βa
ISx
λamc
βamc
ISxmc
ST(68)
0.6820
0.3180
0.4662
0.6666
0.3344
0.5016
ST(128)
0.7005
0.2995
0.4275
0.7197
0.2803
0.3894
ST(192)
0.4982
0.5018
1.0072
0.4917
0.5083
1.0337
(a) Comparison of observed h-ok and image estimated h-dk width averages taken over groups Gk(l). (b) Comparison of observed l-ok and image estimated l-dk length averages taken over groups Gk(l) (see Table 6). The values obtained from digitized leaves were estimated by using ST(68).
(a) Comparison of observed h-ok and image estimated h-dk width averages taken over groups Gk(l). (b) Comparison of observed l-ok and image estimated l-dk length averages taken over groups Gk(l) (see Table 6). The values obtained from digitized leaves were estimated by using ST(68). The values obtained from digitized leaves were estimated by using ST(128).
(a) Comparison of observed h-ok and image estimated h-dk width averages taken over groups Gk(l). (b) Comparison of observed l-ok and image estimated l-dk length averages taken over groups Gk(l) (see Table 6). The values obtained from digitized leaves were estimated by using ST(68). The values obtained from digitized leaves were estimated by using ST(192).
For 1≤x≤67, λl and λh values were greater than those calculated for ST(68). Moreover, λl and λh values obtained using 1≤x≤127 were greater than those obtained using ST(128) but smaller than those obtained for 1≤x≤67. Nevertheless, λl values produced for 1≤x≤196 were smaller than those linked to 1≤x≤128, but generally λh values increased implying a greater concentration of a particular tonality within the range of colors forming the color spectrum of the image. That is, ST(x) broadens or reduces the collection of colors which can be taken into account for pixel identification within the extent of the image. Whenever x 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 δ-h and δ-l values in a direct way. Indeed for ST(192), δ-h was negative (Table 7) which means that most hdjk values were greater than hojk values. Nevertheless, this outcome is limited by the maximum value that x can attain and also by image size. For 1≤x≤255, the value of λl remained unchanged.
By using the proportion values λa, βa, λamc and βamc, we assessed that for ST(192) leaf area was overestimated by the method. And although for ST(128) the method slightly overestimated leaf area, for ST(x) intervals beyond ST(128) 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 ST(x) intervals. We can assess that a ST(128) 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 ISx selection index. Therefore, we suggest that a similarity ST(128) interval must be set to process the present Zostera marina leaf image set through the method proposed.
RMSD calculated by using observed versus image calculated variables.
ST(x)
RMSD(ho,hd)
RMSD(lo,ld)
RMSD(ao,ad)
RMSD(ao,amc)
ST(68)
0.4590
26.4500
360.8746
151.2869
ST(128)
0.4016
12.9587
99.1725
90.6759
ST(192)
0.7303
10.8674
155.3371
160.6715
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 ISxselection index. Moreover, the ST(x) interval selected using the ISx 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 ISx 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 ST(x) 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: lo and ho
Identify the maximum observed leaf length (lmax).
Chose a positive integer n and define a partition of the interval [0,lmax] with norm q=(lmax/n).
Form the collection ⋃1n(Ik) of n disjoint intervals of the form Ik=[q(k-1),qk) with 1≤k≤n. This collection is denoted thought P0lmax.
For each value of the index k identify the group Gk(l) of leaves whose lengths are contained in Ik. Notice that Gk(l) holds the leaves whose size differences are bounded by q.
For each value of the index k obtain and record nk standing for the number of leaves in the group Gk(l).
For each value of the index k introduce an index j such that 1≤j≤nk and label as lojk, hojk and aojk respectively, the straight length, width and area of the jth leaf in Gk(l). The character aojk denotes the associated estimations of leaf area obtained by means of (1).
Form and record the collection ⋃1nGk(l) of all groups of leaves Gk(l). This collection is denoted by means of CG.
Obtain the average length l-ok for each group of leaves Gk(l). That is, calculate and record
(A.1)l-ok=1nk∑1nklojk.
Obtain the average width h-ok for each group of leaves Gk(l). That is, calculate and record
(A.2)h-ok=1nk∑1nkhojk.
Obtain the average area a-ok for each group of leaves Gk(l). That is, calculate and record,
(A.3)a-ok=1nk∑1nkaojk.
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 Cmax of colors.
For processing the digital images of all collected leaves, we choose different intervals of tolerance of similarity ST(x)=[0,x], with the upper bound x satisfying 0≤x≤Cmax-1.
For a picked ST(x) interval, for 1≤j≤nk use the algorithm (a) described in the method section to obtain ldjk and hdjk, which respectively denote the length and width of the image of the jth leaf in Gk(l). Also, obtain adjk and amcjk, 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 l-dk, h-dk, a-dk and a-cmk (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 ST(x) interval and for 1≤j≤nk, calculate the leaf length approximation errors through
(C.1)eljk=lojk-ldjk.
For the picked ST(x) interval and for 1≤j≤nk, calculate the individual leaf width approximation errors through
(C.2)ehjk=hojk-hdjk.
For the picked ST(x) interval and for 1≤j≤nk, calculate the leaf area approximation errors through
(C.3)eajk=aojk-adjk.
For the picked ST(x) interval and for 1≤j≤nk, calculate the leaf area approximation errors linked to the Monte Carlo method through
(C.4)emcjk=aojk-amcjk.
D. Obtain the Average Deviations Produced by the Individual Estimation Errors
For the picked ST(x) interval, obtain the average leaf length deviations δlk, which are calculated by averaging the eljk values. That is,
(D.1)δlk=1nk∑1nkeljk.
We notice that δlk=(l-ok-l-dk) and also that negative values of δlk imply that in a lot Gk(l) most image assessments ldjk overestimate observed lojk values.
For the picked ST(x) interval, obtain the average leaf width deviations δhk, which are calculated by averaging the ehjk values. That is,
(D.2)δhk=1nk∑1nkehjk.
notice that since δhk=(h-ok-h-dk), negative values of δhk imply that in a group Gk(l) most image assessments hdjk overestimate observed hojk values.
For the picked ST(x) interval, calculate the average leaf area deviations δak by averaging the eajk values. That is,
(D.3)δak=1nk∑1nkeajk.
again since δak=a-ok-a-dk negative values of δak imply that in a group Gk(l) most image assessments adjk overestimate observed aojk values.
For the picked ST(x) interval, calculate δ-l, the average value of deviations δlk taken over CG. Calculate also the associated standard deviation σδl.
For the picked ST(x) interval, calculate δ-h, the average value of deviation δhk taken over CG. Calculate also the associated standard deviation σδh.
E. Criteria for Selecting the ST(x) 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 ST(x)=[0,x], identify the leaves satisfying the conditions
(E.1)δ-h≥0,(E.2)δ-l≥0,(E.3)δ-l-σδl≤δlk≤δ-l+σδl,(E.4)δ-h-σδh≤δhk≤δ-h+σδh.
Calculate the proportion λl of leaves in CG that comply with the condition (E.3) through
(E.5)λl=∑k=1n∑j=1nk[ldjk∣leavesinGkthatcomplywithcondition(E.3)ldjk]×(∑k=1n∑j=1nklojk)-1.
Calculate the proportion λh of leaves in CG that comply with the condition (E.4) through
(E.6)λh=∑k=1n∑j=1nk[hdjk∣leavesinGkthatcomplywithcondition(E.4)hdjk]×(∑k=1n∑j=1nkhojk)-1.
Calculate the proportion θl of leaves in CG that do not comply with the condition (E.3) through,
(E.7)θl=1-λl.
Calculate the proportion θh of leaves in CG that do not comply with the condition (E.4) through,
(E.8)θh=1-λh.
Obtain the concomitant proportions of leaves in CG that provide consistent leaf area estimations by the proxy of (1)
(E.9)λa=∑k=1n∑j=1nk[adjk∣leavesinGkthatcomplywithcondition(E.3),(E.4)andeajk≥0]×(∑k=1n∑j=1nkaojk)-1.
Calculate the proportion of leaves in CG for which image estimated blade length and width measurements overestimate leaf area calculated through (1)
(E.10)βa=1-λa.
The proportions λamc and βamc equivalent to λa and βa respectively but linked to leaf area estimation by Monte Carlo method (cf. (2)) that is,
(E.11)λamc=∑k=1n∑j=1nk[amcjk∣leavesinGkthatcomplywithcondition(E.3),(E.4)andemcjk≥0]×(∑k=1n∑j=1nkaojk)-1,(E.12)βamc=1-λamc.
Conditions (E.1) and (E.2) grant bounded estimation errors, for h and l respectively. Moreover, the groups of leaves that also satisfy conditions (E.3) and (E.4) can be identified as those groups for which image ld and hd estimations are closer to directly obtained lo and ho measurements. Therefore, groups in CG, 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.
Echavarria-HerasH.Solana-ArellanoE.Leal-RamirezC.Franco-VizcaínoE.The length-times-width proxy for leaf area of eelgrass: criteria for evaluating the representativeness of leaf-width measurements20112176046132-s2.0-8215519588510.1002/aqc.1219WilliamsT. P.BubbJ. M.LesterJ. N.Metal accumulation within salt marsh environments19942852772892-s2.0-002830767010.1016/0025-326X(94)90152-XShortF. T.ShortC. A.KennedyV. S.The seagrass filter: purification of coastal water1984Academic Press395413WardL. G.Michael KempW.BoyntonW. R.The influence of waves and seagrass communities on suspended particulates in an estuarine embayment1984591–4851032-s2.0-0021642476FonsecaM. S.FisherJ. S.A comparison of canopy friction and sediment movement between four species of seagrass with reference to their ecology and restoration198629152210.3354/meps029015TuyaF.Ribeiro-LeiteL.Arto-CuestaN.CocaJ.HarounR.EspinoF.Decadal changes in the structure of Cymodocea nodosa seagrass meadows: natural vs. human influeces20141374149DennisonW. C.RonaldC. P.McRoyC. P.Leaf production19907779GaeckleJ. L.ShortF. T.A plastochrone method for measuring leaf growth in eelgrass, Zostera marina L2002713123712462-s2.0-0037742213Echavarria-HerasH. A.Solana-ArellanoM. E.LeeK. S.HosokawaS.Franco-VizcainoE.An evaluation of leaf biomass: length ratio as a tool for nondestructive assessment in eelgrass (Zostera marina L.)2012201254373010.1100/2012/543730Echavarría-HerasH.Solana-ArellanoE.Franco-VizcaínoE.An allometric method for the projection of eelgrass leaf biomass production rates2010223158652-s2.0-7264910181510.1016/j.mbs.2009.10.008IgathinathaneC.WomacA. R.SokhansanjS.PordesimoL. O.Mass and moisture distribution in aboveground components of standing corn plants2006491971062-s2.0-33645456964Rico-GarcíaE. .Hernández-HernándezF.Soto-ZarazúaG. M.Herrera-RuizG.Two new methods for the estimation of leaf area using digital photography200911397400Femat-DiazA.Vargas-VazquezD.Huerta-ManzanillaE.Rico-GarciaE.Herrera-RuizG.Scanner image methodology (SIM) to measure dimensions of leaves for agronomical applications20111010184018472-s2.0-79952637803PatilS. B.BodhS. K.Image processing method to measure sugarcane leaf area20113863946400LüC.RenH.ZhangY.ShenY.Leaf area measurement based on image processing2Proceedings of the International Conference on Measuring Technology and Mechatronics Automation (ICMTMA '10)March 2010IEEE Computer Society5805822-s2.0-7795317371510.1109/ICMTMA.2010.141RamfosA.GazisA.KatselisG.Development and evaluation of an automated digital image analysis software for obtaining seagrass leaf metrics2012556601610JuneauK. J. .TarasoffC. S.Leaf area and water content changes after permanent and temporary storage201278e42604EngleJ. M.MillerK. A.GarcelonD. K.ScwemmC. A.Distribution and morphology of eelgrass (Zostera marina L.) at the California channel islandsProceedings of the 6th California Islands Symposium2005405414National Park Service Technical Publication CHIS-05-01OrthR. J.CarruthersT. J. B.DennisonW. C.DuarteC. M.FourqureanJ. W.HeckK. L.Jr.HughesA. R.KendrickG. A.KenworthyW. J.OlyarnikS.ShortF. T.WaycottM.WilliamsS. L.A global crisis for seagrass ecosystems200656129879962-s2.0-3384553318110.1641/0006-3568(2006)56[987:AGCFSE]2.0.CO;2BjörkM.ShortF.McLeodE.BeerS.Managing seagrasses for resilience to climate ChangeIUCN Resilience Science Group Working Paper Series no 3, 2008