Growth analysis and selection of sea cucumbers (
Stichopus japonicas
) is typically performed through length or weight measurements. However, because sea cucumbers continuously change shape depending on the external environment, weight measurement has been the preferred approach. Weight measurements require extensive time and labor, moreover it is often difficult to accurately weigh sea cucumbers because of their wet surface. The present study measured sea cucumber features, including the body length, width, and thickness, by using a vision system and regression analysis to generate
R
^{2}
values that were used to develop a weight estimation algorithm. The
R
^{2}
value between the actual volume and weight of the sea cucumbers was 0.999, which was relatively high. Evaluation of the performance of this algorithm using crossvalidation showed that the root mean square error and worstcase prediction error were 1.434 g and ±5.879 g, respectively. In addition, the present study confirmed that the proposed weight estimation algorithm and single slide rail device for weight measurement can measure weights at approximately 4,500 sea cucumbers per hour.
1. Introduction
The lengths and weights of marine animals are typically measured to determine growth and to select individuals. Weight measurements by scale are the most popular approach
[1]
. The weight measurement method requires frequent calibrations — i.e., taring the scale — and often hinders accurate and rapid measurement because of the webbed surface of the measured object.
To improve these methods, previous studies measured the lengths
[2

5]
and areas
[6

7]
of fish in twodimensional images generated by a computer vision system; fish selection and weight estimation were then conducted. Lee et al. (2001) performed threedimensional measurements and restructured a twodimensional crosssection photography scan into a threedimensional surface using charge coupled device (CCD) cameras and a multiline laser
[8

9]
. In addition, it is possible to calculate the area of an oyster body using the cubic splines method
[10]
to estimate weight. However, unlike other marine animals, the sea cucumber (
Stichopus japonicas
) changes its shape into a footballlike
[11]
or globular structure when it shrinks
[12]
; moreover, its shape is irregular when it is fully extended. Thus, it is difficult to apply the existing method of weight estimation by measuring the body length (
BL
), body width (
BW
), thickness (
T
), and view area
[13]
. In Japan, the body size of the sea cucumber is measured after the administration of
l
menthol anesthesia. However, because this method cannot be easily applied at practical sites, Yamana and Hamano derived a regression equation that uses a digital still camera to measure the
BL
and
BW
of the sea cucumber at any given time; it then converts the
BL
and
BW
to the image that is equivalent to its anesthetized condition
[14]
. Yingst formulated the size index (SI) =
BL
×
BW
to determine body size of the sea cucumber
[15]
. Moreover, Mitsukuri proposed an equation that facilitates calculation of the sea cucumber’s body size using its square length and width
[16]
. However, the above calculation methods are inaccurate because they include no information on the thickness of the sea cucumber; rather, they are simply based on observations of shape changes in which the animal gradually becomes slender, during its extension. Body volume of the animal seems to be relatively stable during body extension; however, the influence of water storage in its respiratory tree remains to be clarified
[14]
.
Existing measurement methods using sea cucumber volume or weight are impractical when a large quantity of farmraised sea cucumbers are to be measured, particularly in a short time period; a substantial amount of time and labor are required for accurate measurements. Therefore, to generalize their everchanging shapes, quickly and accurately analyze their growth or select them at any given time, a technology that automatically estimates sea cucumber weight is required.
The purpose of this study was to develop an algorithm that estimates sea cucumber weight by measuring its features, which include
BL
,
BW
,
T
, and volume, using the image processing technology of vision system. To simply perform the measurement of the sea cucumber — instead of using the existing method of measuring its
BL
and
BW
— it was assumed that the early shape of the sea cucumber in air is a half oblate ellipsoid. Therefore, the volume of the sea cucumber was measured using the
BL
,
BW
, and
T
, which were measured at the time of sample, and by using the volume equation for a modified ellipsoid. Accordingly, regression analyses of the calculated and actual volumes of sea cucumbers were performed for optimization. For optimal design of the weight estimation algorithm, regression analysis was conducted on each feature of the sea cucumbers against the weight, and the
R
^{2}
value and regression equation were examined. To confirm the validity and performance of the final weight estimation algorithm, replicates were performed using the single slider rail, which was produced in a laboratorysized system for the present study.
2. Material and Methods
 2.1 Experimental animals: Sea cucumeber
In January 2014, 350 sea cucumbers (weight range: 0.23–92.56 g) were collected by divers from the southern coast of Korea. The depth of the sampling location did not exceed 11m and the seawater temperature at the time of sampling was 13.8 ℃. The sea cucumbers were transferred to the nearby laboratory of the Aquaculture Management Division and maintained in a tank, which was equipped with flowing 100 μm filtered seawater at ambient temperature (14 ℃). All animals were starved for 24 h before the experiment to minimize the effect of excretion on weight.
 2.2 Measurement of actual volume and weight
Actual volumes and weights of sea cucumber measurements were obtained using a density meter (Shinko Denshi Co., Ltd., DME220) and digital balance (A&D Co., Ltd., ER180A), respectively, after blotting the water on the external surface of the animals using a paper towel. The average of the actual volumes and weights were obtained after 20 replicate measurements. Using Archimedes’ theory, the actual volumes were measured by applying to Eq. (1) the weights of the sea cucumbers in air and in distilled water, and the water density of the distilled water. The water density in this study was considered as 1, because the density of pure water is 0.997 g/cm
^{3}
at room temperature (25 ℃)
[10]
.
 2.3 Image processing to measureBL,BWandT
Fig. 1
shows an imageprocessing flowchart for measuring the
BL
,
BW
, and
T
of the sea cucumber based on the vision system.
As shown in
Fig. 1(a)
, preprocessing involved capturing images of the sea cucumber, which were projected by the backlights of the CCD cameras installed at the bottom and sides of a tank to measure the sea cucumber features. The brightness information of the blackandwhite images was then converted to an 8bit gray scale. At that point because the
BL
,
BW
, and
T
should be measured after removing the sea cucumber thorns
[14]
, 8bit gray scalebased images were converted to 9×9 matrix structure images by performing a closing calculation using morphological image processing to partially remove the thorns in the images
[21]
. Finally, in the binarization process and based on the 60 level, the pixel brightness was converted to a range between 0 and 255 to improve clarity of the edge lines in the 8bit projected bottom and side images. The base coordinate system was then configured to each output image. As shown in
Fig. 1(b)
, to measure the
BL
,
BW
, and
T
of the sea cucumber after preprocessing, the closed curve was detected from the projected bottom and side images; the center point and angles of the closed curve area were then extracted, whereby a base coordinate system was drawn.
Image processing flowchart: (a) preprocessing; (b) measuring process.
Fig. 2
shows images captured while measuring the
BL
of the sea cucumber. The
BL
of the sea cucumber was measured from the estimated centerline in the projected bottom image. The measurement process was as follows:

1) As shown inFig. 2(a), the circular and linear patterns were extracted from the bottom image of the sea cucumber, which was projected on thexyplane; the bottom image was scanned from the outside to the center point of the area. Then, 100 edge points were extracted.

2) Of the extracted edge points, as shown inFig. 2(b), the left coordinate (x0,y0) and right coordinate (x9,y9), which were the farthest from the area’s center point, were selected.

3) Thexcoordinate area of the selected two points,x0x9, was divided by eight to obtain a straightline equation (x=x1,x2…x8). Then, as shown inFig. 2(c), the center coordinates (x1,y2), … (x8,y8) of the edge points, where brightness was considerably different on the straight line, were detected.

4) TheBLof the sea cucumber was measured by applying each detected coordinate to Eq. (2).
Process of measuring sea cucumber BL: (a) detection of edge points; (b) detection of the left and right coordinates that are the farthest from the center point; (c) connection of the detected coordinates.
Fig. 3
shows the image captured in the sea cucumber
BW
measurement process which was as follows:

1) As shown inFig. 3, in reference to the center point angle in the bottom image, in which the sea cucumber was projected on the xy plane, the straightline equation passing the center of gravity was detected.

2) The two edge points, whose brightness levels were substantially different on the straightline were detected to measure the sea cucumberBW.

3) The distance between the coordinates of the two edge points was then designated as theBW.
Measurement of sea cucumber BW.
Fig. 4
shows the image captured in the process of measuring the sea cucumber
T
.
Fig. 5
shows the thickness errors according to measurement points observed while moving a 10 mm regular hexahedron standardized specimen along the
y
axis. The measurement point was at the area center point of the sea cucumber’s bottom image.
Measurement of the sea cucumber T.
As shown in
Fig. 5
, measuring thickness using CCD cameras, in general, led to variations because of the perspective between the point of the measured specimen and the CCD camera installed at the side. Hence, an additional modification process was required for this study. To minimize variations when measuring
T
, the linear regression analysis between the location of the sea cucumber and the measured
T
was performed. The revision equation was derived, as shown in Eq. (3), where
t
is the measured thickness of the sea cucumber;
y
is the location of the sea cucumber; and
α
and
β
are regression coefficients, which were 0.044 and 0.009, respectively.
Variations in thickness according to measurement points.
The process of measuring sea cucumber
T
was as follows:

1) As shown inFig. 4, 100 edge points were detected in the side image. The sea cucumber was projected in thexz plane by scanning the image from white to black in the direction of the center of gravity using the straightline pattern.

2) Of the detected 100 edge points, the farthest edge point from thexyplane and center of gravity — i.e., the location of the sea cucumber from the bottom image of the sea cucumber projected on thexyplane — were simultaneously detected.

3) The sea cucumberTwas measured by the distance between the detected edge point andxyplane, and it was calibrated as shown in Eq. (3).
The sea cucumber shrinks to a half oblate ellipsoid at its maximum as it is handled in the air. Thus, assuming that the early shape of the sea cucumber is a half oblate ellipsoid, as shown in
Figs. 2
,
3
, and
4
, the
BL, BW
, and
T
were simultaneously measured using the vision system. They were then calculated using Eq. (4), which is a modified ellipsoid equation. The actual volume (
A_{v}
) and calculated volume (
C_{v}
) were compared and verified through linear regression analysis using Eq. (5), where
α
and
β
are regression coefficients.
 2.4 Regression analysis for designing a sea cucumber weight estimation algorithm
For design optimization of the sea cucumber weight estimation algorithm, regression analyses on the measured
BL
,
BW
,
T
, and
A_{v}
, and on the
C_{v}
were performed based on the assumption that the sea cucumber form is a half oblate ellipsoid. The
R
^{2}
and regression equation were then obtained. The power curve regression model was used for the regression analyses of the
BL
s,
BW
s, and
T
s against the weights; the linear regression model was used for the actual sea cucumber volumes against the weights.
 2.5 Experimental systems
We produced a vision systembased device (size: 600×500×1300 mm) for measuring sea cucumber features,
BL
s,
BW
s,
T
s, and volumes, as shown in
Fig. 6
. The vision system (PV200, Panasonic) and CCD cameras were installed at the bottom and sides, respectively, to capture twodimensional images of the bottom and sides of the sea cucumber. These images were projected on the twodimensional plane. In the opposite direction of each CCD camera, a 30 W infrared backlight was installed to generate a silhouette of the sea cucumber. A data collection device (National Instruments, cDAQ9174) was installed to collect the visual information required to measure the
BL
,
BW
, and
T
and to calculate the volume, based on LabVIEW software (version: LabVIEW 2013, National Instruments).
Device for measuring sea cucumber features.
Fig. 7
shows the single slide rail system for measuring sea cucumber weights that we produced in a laboratory size to verify the performance of the optimized sea cucumber weight estimation algorithm through repeated testing. The system components were the same as those of the sea cucumber feature measurement device, except for the addition of an installed photosensor (Autonics Co., LTD, BYS500TDT) that generates trigger signals to automatically capture images. The weight of each sea cucumber whose respiratory tree was emptied of water in advance was measured by dropping it in the single side rail linked to the tank. At that point, the optimized algorithm for sea cucumber’s weight estimation was used to continuously measure the weights of sea cucumbers moving in the rail.
Single slide rail device for measuring sea cucumber weights.
3. Results and Discussion
 3.1 Repeatability of measurement
The inspection or measurement method of a certain object based on images leads to respective inspection or measurement errors depending on the brightness of lighting and external environments. Therefore, the repeatability of the designed algorithm for sea cucumber features was measured to verify its reliability. After selecting standardized regular hexahedron specimens with a regular surface (dimensions: 10×10×10 mm), their lengths, widths, and heights were simultaneously measured 300 times within an interval of a second under the same measurement conditions.
Fig. 8
shows the results of 300 measurements of the standardized regular hexahedron specimens.
Table 1
indicates the mean, standard deviation, and coefficient of variation of the features. As a result, the means, standard deviations, and coefficients of variation were within 10.004 mm, 0.017mm, and 0.170%, respectively, thereby confirming high repeatability.
Test results of repeated measurements of the standardized regular hexahedron specimens.
Results of measurements: mean, standard deviation, and coefficient of variation
Results of measurements: mean, standard deviation, and coefficient of variation
 3.2 Design of sea cucumber weight estimation algorithm
For optimized design of the sea cucumber weight estimation algorithm, regression analyses were performed on the
BL
s,
BW
s,
T
s, actual volumes, and calculated volumes — where sea cucumber form was assumed to be a half oblate ellipsoid — against the sea cucumber weights.
Figs. 9
through
11
show the correlations of the
BL
,
BW
, and
T
against the weight. Each correlation was nonlinearly distributed. The power curve regression model was used to examine the regression equation and
R
^{2}
value, as shown in
Table 2
. Here,
y
denotes weight,
x
represents the measurement item, and
α
and
β
are regression coefficients. The test results confirmed that the
R
^{2}
values of the
BL
,
BW
, and
T
against the weight of the sea cucumber were all >0.871, where the
R
^{2}
value between the
T
and weight was rather high (0.925).
Equation fit result for sea cucumber weight compared toBL, BW, andT
Equation fit result for sea cucumber weight compared to BL, BW, and T
However,
Figs. 9
,
10
, and
11
suggest that the smaller size of the sea cucumber indicates a smaller weight deviation. Thus, the larger the size of the sea cucumber is, the bigger the weight deviation is, implying that a bigger error would occur in estimating the weight of a large sea cucumber. In addition, the weight against the
BL
,
BW
, and
T
was nonlinear; therefore, a more precise measurement technology would be required when classifying sea cucumbers according to their weights.
Plot of power curve fit of weight compared with BL and weight.
Plot of power curve fit of weight compared with BW and weight.
Plot of power curve fit of weight compared with T and weight.
Figs. 12
and
13
show the correlations between the actual volume and calculated volume where the sea cucumber was assumed to be a half oblate ellipsoid. Each correlation was linearly distributed. The regression equation and
R
^{2}
value were investigated using the linear regression model, as shown in
Table 3
.
Plot of linear fit of weight compared with A_{v} and weight.
Plot of linear fit of weight compared with C_{v} and weight.
Equation fit result for sea cucumber weight compared toAvandCv
Equation fit result for sea cucumber weight compared to A_{v} and C_{v}
The
R
^{2}
value between the actual volume and weight was 0.999, which was relatively very high. The
R
^{2}
value between the calculated volume and weight was 0.995, which was likewise high. The results of comparing each regression equation confirmed that the differences in the slope and yintercept were approximately 0.017 and 0.944, respectively.
As shown in
Fig. 14
, we ran a linear regression analysis using Eq. (5) to compare and verify the correlation between the actual volumes and calculated volumes.
Plot of linear fit of weight compared with C_{v} and A_{v}.
The coefficients of variation
α
and
β
were 0.984 and 0.912, respectively, where the
R
^{2}
value was 0.995. The regression analysis, where
β
was 0, showed
α
and
R
^{2}
were 1.010 and 0.997, respectively. Thus, the regression equation slope and
R
^{2}
between the actual volume and calculated volume were both close to 1. Moreover, the
R
^{2}
value was also high. Based on this finding, we confirmed that there was no significant difference between the actual volume and the calculated volume in experimental range.
The results of performing the regression analyses of the
BL
,
BW
,
T
,
A_{v}
, and
C_{v}
against weight confirmed that the method of measuring the actual sea cucumber volume was effective in estimating sea cucumber weight. Therefore, using the regression equation between the actual volume and weight, as shown in
Table 3
and Eqs. (4) and (5), and the algorithm — which estimates sea cucumber weight by measuring its volume (Eq. (6), where
C
_{1}
and
C
_{2}
denote the constants calculated with the regression coefficients between the actual volume and weight and between the actual and calculated volume) — the
C
_{1}
and
C
_{2}
of our experiments were determined to be 0.534 and 0.958, respectively.
The performance of this algorithm was evaluated through leaveoneout crossvalidation
[22
,
23]
, moreover, as shown in
Fig. 15
, the root mean square error and weight estimation error spread were 1.434 g and ±5.879 g, respectively. These finding prove that the accuracy of this algorithm improved compared to that of existing estimation methods that measure threedimensional projected volume of the fish
[24
,
25]
. It was confirmed that the developed weight estimation algorithm and single slide rail system for weight measurement can weigh approximately 4,500 sea cucumbers per hour.
Sea cucumber weight prediction error distribution, evaluated by leaveon out crossvalidation.
4. Conclusion
Regression analysis between each feature and the weight of the sea cucumber to estimate weight showed that the correlation between the actual volume and weight of the sea cucumber was high. The proposed algorithm was shown to estimate the weight of sea cucumbers by measuring their volumes using the vision system. This method calculated the volume of the sea cucumber by assuming its early form as a half oblate ellipsoid; moreover, it proved to be more convenient than existing measurement methods in which the line laser and vision system were combined. Moreover, optimization was performed through linear regression analysis between the calculated and actual volume of the sea cucumber.
The developed sea cucumber weight estimation algorithm can rapidly and accurately measure the growth of sea cucumbers that are raised within the recently developed integrated multitrophic aquaculture (IMTA) system; facilitates automated management. The proposed algorithm shows potential applications to the selection of sea cucumbers.
Acknowledgements
This work was supported by a grant (RP2014FE013) from the National Fisheries Research and Development Institute (NFRDI), Republic of Korea.
BIO
Donggil Lee received B.S., M.S. and Ph.D. degree in Electrical Engineering from Pusan National University in 2008, 2010 and 2013, respectively. He has been a senior researcher in NFRDI, Busan, Korea. His research interests are in the field of vision algorithm, machine vision, automatic control system design, circuit design, and robotics.
Seonghun Kim received his BE and ME degree in Fisheries Physics in Pukyong National University. in 2000 and 2002, respectively and his Ph.D. degree in Fisheries Engineering from Hokkaido Univ. in 2008. He is a researcher at Fisheries system engineering Division, NFRDI, busan, Korea. His research interests include Fishing gear design, mesh size selectivity, echofriendly fishing gear.
Miseon Park received her Ph.D. in Shellfish Disease from Pukyong National University in 1989. Since April 1982, she has been a junior researcher, a senior researcher, and a director general in NFRDI, Busan, Korea. Her research interests include shellfish aquaculture, shellfish disease, and Integrated MultiTrophic Aquaculture.
Yongsu Yang received his Ph.D. degree in Fishery Engineering from Jeju National University in 2000. He has been a senior researcher in NFRDI (National Fisheries Research & Development Institute), Busan, Korea. His research interest is the fishing gear design, and underwater acoustics.
Fulton T. W.
1904
The Rate of Growth of Fishers
Fisheries Board Scotland
Edinburgh
141 
241
Ibrahim M. Y.
,
Wang J.
2009
“Mechatronics applications to fish sorting  part 1: fish size identification”
IEEE International Symposium on Industrial Electronics
1978 
1983
Jeong S. J.
,
Yang Y. S.
,
Lee K.
,
Kang J. G.
,
Lee D. G.
2013
“Visionbased automatic system for noncontact measurement of morphometric characteristics of flatfish”
Journal of Electrical Engineering & Technology
8
(5)
1194 
1201
DOI : 10.5370/JEET.2013.8.5.1194
White D. J.
,
Svellingen C.
,
Strachan N. J. C.
2006
“Automated measurement of species and length of fish by computer vision,”
Fish. Res.
80
203 
210
DOI : 10.1016/j.fishres.2006.04.009
Balaban M. O.
,
Şengör G. F. Ünal
,
Soriano M. G.
,
Ruiz E. G.
2010
“Using image analysis to predict the weight of Alaskan salmon of different species”
J. Food Sci.
75
E157 
E162
DOI : 10.1111/j.17503841.2010.01522.x
Lee D. J.
,
Lane R. M.
,
Chang G. H.
2001
“Threedimensional reconstruction for highspeed volume measurement,” Intelligent Systems and Smart Manufacturing
International Society for Optics and Photonics
Lee D.J.
,
Eifert J.
,
Zhan P.
,
Westover B.
2003
“Fast surface approximation for volume and surface area measurements using distance transform”
Opt. Eng.
42
2947 
2955
DOI : 10.1117/1.1605737
Damar S.
,
Yagiz Y.
,
Balaban M.O.
,
Ural S.
,
Oliveira A. C. M.
,
Crapo C.A.
2007
“Prediction of oyster volume and weight using machine vision”
J. Aquat. Food Prod. Technol.
15
3 
15
DOI : 10.1300/J030v15n04_02
Feindel S.
,
Bennett T.
,
Kanwit K.
2011
“The Maine sea cucumber (Cucumaria frondosa) fishery”
Submitted to Standing Legislative Committee on Marine Resources. Department of Marine Resources
Maine
http://www.maine.gov/dmr/cukes/feinde/eta/2011.pdf
Choe S.
,
Ohshima Y.
1961
“On the morphological and ecological differences between two commercial forms. ‘green’ and ‘red’, of the Japanese common sea cucumber, Stichopus japonicas Selenka,”
Nippon Suisan Gakkaishi
27
97 
106
DOI : 10.2331/suisan.27.97
Gümüş B.
,
Balaban M. O.
2010
“Prediction of the weight of aquacultured rainbow trout (Oncorhynchus mykiss) by image analysis,”
J. Aquat. Food Prod. Technol.
19
227 
237
DOI : 10.1080/10498850.2010.508869
Yamana Y.
,
Hamano T.
2006
“New size measurement for the Japanese sea cucumber Apostichopus japonicas (Stichopodidae) estimated from the body length and body breadth”
Fish. Sci.
72
585 
589
DOI : 10.1111/j.14442906.2006.01187.x
Yamana Y.
,
Hamano T.
2006
“New size measurement for the Japanese sea cucumber Apostichopus japonicas (Stichopodidae) estimated from the body length and body breadth,”
Fish. Sci.
72
585 
589
DOI : 10.1111/j.14442906.2006.01187.x
Mitsukuri K.
1903
“Notes on the habits and lifehistory of Stichopus japonicus Selenka,”
Dong Y. W.
,
Dong S. L.
,
Tian X. L.
,
Wang F.
,
Zhang M. Z.
2006
“Effects of diel temperature fluctuations on growth, oxygen consumption and proximate body composition in the sea cucumber Apostichopus japonicus Selenka”
Aquaculture
255
514 
521
DOI : 10.1016/j.aquaculture.2005.12.013
Zamora L. N.
,
Jeffs A. G.
2012
“Feeding metabolism and growth in response to temperature in juveniles of the Australasian sea cucumber, Australostichopus mollis”
Aquaculture
15
358 
359
Ji T.
,
Dong Y.
,
Dong S.
2008
“Growth and physiological responses in the sea cucumber, Apostichopus japonicus Selenka: aestivation and temperature”
Aquaculture
283
180 
187
DOI : 10.1016/j.aquaculture.2008.07.006
Paltzat D. L.
,
Pearce C. M.
,
Barnes P. A.
,
McKinley R. S.
2008
“Growth and production of California sea cucumbers (Parastichopus californicus Stimpson) cocultured with suspended Pacific oysters (Crassostrea gigas Thunberg)”
Aquaculture
275
124 
137
DOI : 10.1016/j.aquaculture.2007.12.014
Gonzalez R.
,
Woods R.
2007
Digital Image Processing
Prentice Hall
New Jersey
Theodoridis S.
,
Koutroumbas K.
1999
Pattern recognition
Academic Press
San Diego, California
Mathiassen J. R.
,
Misimi E.
,
Toldnes B.
,
Bondø M.
,
Østvik S. O.
2011
“Highspeed weight estimation of whole herring (Clupea harengus) using 3D machine vision”
J. Food Sci.
76
E458 
E464
DOI : 10.1111/j.17503841.2011.02226.x