Abstract
Many previous studies on cotton maturity used a sole parameter to rank the maturity of a cotton sample containing a large number of fibers. In light of the complexity of maturity distributions, the sole-parameter approach does not appear to be reliable and rational for cotton maturity evaluation. More distributional parameters should be examined and included in the new classification methods. This paper (1) introduces important changes in the image-analysis algorithms for cotton cross-section measurements to enhance the consistency of fiber detections in order to reduce the bias on immature fibers, (2) investigates the characteristics and patterns of cotton maturity distributions, and (3) presents the experimental results on the cross-section images selected from seven cotton varieties that have a wide range of maturities. It is found that the skewness of a maturity distribution is an essential parameter for classifying the distribution pattern and that the dead fiber content and the mature fiber content are the important distributional parameters for assessing cotton maturity.
There is a relatively long history of efforts to establish a set of standard cottons for calibrating fiber maturity and fineness measurements. In 1956, Lord 1 established 100 cottons to calibrate micronaire. In the 1980s, the International Textile Manufacturers Federation (ITMF), in cooperation with other organizations, developed a set of nine calibration cottons for use in large-scale maturity-round trials. 2 In 1999, Thibodeaux and Rajasekaran 3 used a commercial image-analysis system to measure the cross-sectional images of approximately 50 cotton varieties representing a wide range of genetic finenesses that were grown and hand-harvested in the USA. In the 2000s, Hequet's and Thibodeaux's groups performed a multi-year project to create a larger-scale cotton-maturity reference by analyzing images of the 104 cotton varieties 4 collected worldwide using a new fiber cross-sectioning protocol 5 and customized fiber image-analysis software (FIAS). 6 Now, there is a consensus in the cotton research community that fiber cross-sectional analysis provides fundamental measurements that are directly related to cotton maturity and that this is the most reasonable approach for developing the cotton-reference database. The FIAS helped expedite the analysis on a large quantity of fiber cross-section images and improved the data reliability for calibrating other testing methods. 4 However, a recent independent study revealed that there is an overestimation of the maturity level of bale cotton by around 9% using the laboratory protocol and the FIAS, because 10–40% of immature fibers could not be detected correctly in the image processing. 7
Immature fibers, especially dead fibers, have thinner walls, and thus they are more easily scratched or shredded by the cutting blade. They are also likely to be folded up transversely when the lumens collapse, which increases the difficulty of edge detection. Based on the experience of developing and using the FIAS, we know that the FIAS lacks functions that can automatically amend broken edges of cross-sections and that unclosed cross-sections become oblivious in a segmentation routine called “background flooding.” Most of the ignored fibers are immature fibers, which cause a systematic bias in fiber detections, and when a large number of fibers cannot be detected correctly, the cumulative measurements will not realistically reflect maturity distributions of the sample. In addition, the mean value of the cross-section data is often used to indicate the maturity level of a sample.4,6 In light of the complexity of maturity distributions, the sole-parameter approach does not appear to be reliable and rational for ranking the maturity among different samples.
In this paper, we first present the modified algorithms made in the current version of FIAS to improve its robustness and accuracy of detecting fiber cross-sections, including those with obscure and incomplete boundaries, in order to increase the total number of valid measurements and reduce the bias on immature fibers. We then examine the basic patterns of maturity distributions based on the skewness of the distribution and propose a way to calculate various fiber contents (such as dead fiber and mature fiber contents (MFCs)) in the distribution. The paper also presents comparisons between the previous version and the current version of FIAS in fiber detections and in maturity statistics.
Method
In the cross-sectional analysis, the maturity of individual fibers is often measured by a shape factor called circularity or the theta (θ) value.3,4,6,8 The θ is determined by the perimeter (P) and the net area (A) of a fiber cross-section (θ = 4πA/P2), 9 which indicates the degree of wall thickening. Thus, in this paper θ is referred to as cotton maturity. Increasing the accuracy of locating fiber boundaries and lumens is paramount for obtaining reliable maturity measurements.
The images of fiber cross-sections used in this study were provided by the Fiber and Biopolymer Research Institute (FBRI) of Texas Tech University. A total of 15,473 images from the seven cotton varieties were tested with the previous and current FIAS software to examine the differences in the number of detected fibers, the maturity correlation with the data obtained from the Advanced Fiber Information System (AFIS), and the High Volume Instrument (HVI), as well as the distributional parameters. Each of the seven varieties was assigned a unique four-digit symbol (e.g., 2996, 2999…).
Modification of fiber-detection algorithms
To segment individual fiber cross-sections in the image, the FIAS has four major image-processing routines: (1) dynamic thresholding; (2) background flooding; (3) skeletonizing fibers; and (4) identifying lumens. 7 These routines appear to be effective when image illumination is uniform and image contrast is sufficient, but they are not robust enough for fibers that have low contrast, broken edges, or distorted cross-sections due to self-folding/wrapping of fiber ribbons. Routines (2) and (3) are more standardized procedures with fewer adjustable parameters, and therefore we focused on improving routines (1) and (4), which are responsible for locating accurate fiber boundaries and lumens. An edge-amending routine was also added to connect narrowly broken boundaries to prevent them from disappearing in the background flooding.
Adaptive thresholding
Illumination intensity in a captured eight-bit gray-scale image can vary drastically. There is no universal threshold that can generate consistent results across the entire image. The dynamic threshold method with adaptive parameters is a reasonable way to convert an eight-bit gray-scaled image into a binary one. In the previous FIAS,
6
a sub-window of a given size was used to calculate the threshold (Ti) to separate object pixels in the window from the background. Ti was determined by both the mean Mi and the standard deviation SDi of pixel intensities in the window, as follows:
Adaptive thresholding: original image (a); previous binary image (b); and improved binary image (c).
The binary image of Figure 1(a), with the aid of dynamic coefficient ci, is displayed in Figure 1(c), in which the large proportion of the noisy image (Figure 1(b)) was reduced. Fiber boundaries in this image are also more complete than those in Figure 1(b).
Amending broken edges
As seen in Figure 1(c), the binary image still contains much background noise and small objects such as broken edges of fibers, which need to be removed so that only enclosed fiber boundaries are present. This can be done efficiently with a procedure called “background flooding,” which fills the white background with black pixels. Broken fibers are also absorbed by the “flooding” when their inner white regions (fiber walls) are occupied by black pixels. In order to keep those narrowly broken fibers from being eliminated, a step to amend almost-connected edges needs to be added.
There is no generic algorithm for detecting the ends of broken boundaries. The edge-amending routine had to be developed on a case-by-case basis after fiber outer boundaries were traced, and it mainly involved three steps, as follows. Step 1 is to check whether a fiber boundary has a “dead-end” pixel by counting the number of its neighboring pixels. A dead-end pixel is the one that has one neighbor in the boundary chain. Step 2 involves searching for another “dead end” in the same chain. Finally, step 3 consists of checking that the open distance between the two is within an allowable limit and that the boundary length between the two exceeded the threshold on the fiber perimeter. The last step is meant to prevent a connection between two dead ends that are not on the same boundary or are too far apart. The connection distance was limited to five pixels in the new program, and it was simply made with a straight line drawn automatically to join the two identified dead ends.
Figure 2(a) shows a few examples of typical broken boundaries, and Figure 2(b) shows the connections made. The connection lines do not alter the fibers' original shapes but prevent them from being immersed with the background flooding.
Broken boundaries (a) and amended connections (b).
Identifying lumens
In a cross-section image, often some fibers do not possess visible lumens. As seen in Figure 3(a), fibers 1, 2, and 3 are immature fibers with lumens that are totally collapsed; fibers 4, 5, and 6 do not exhibit clear lumens because of low contrasts; and fibers 7, 8, and 9 show partial lumens. Lumen areas in these fibers can be readily overestimated or underestimated.
Original image (a) and detected cross-sections (b).
Once a fiber boundary is located, another threshold is calculated by using only the pixels within the boundary so that the lumen can be more precisely segmented with localized parameters (the mean and standard deviation of the pixel intensities). Since a lumen should be situated in the fiber center, the skeleton of the fiber—that is, the middle axis of the fiber—can be used to locate the lumen if multiple black areas are present inside the boundary. The skeletons of fibers with lumens that are not detected after the thresholding (one-pixel-thick line segments) are placed as lumens. This localized thresholding and skeletonizing enhances the accuracy of lumen size and location (Figure 3(b)).
Comparison of fiber detection results of two versions of fiber image-analysis software
Characteristics of maturity distributions
A regular image of 640 × 480 pixels normally contains a number of fiber cross-sections that varies between 10 and 100. The θ values of cross-sections in multiple images are measured individually, and the frequencies are counted on a scale of [0, 1]. After a sufficient number of fiber cross-sections are measured from multiple images, the descriptive statistics of fiber maturity data (θ), including the mean Mθ, the standard deviation (SDθ), the skewness (Sθ), and the kurtosis (Kθ), can be used to examine the features of maturity distribution. For a unimodal distribution, the mean and standard deviation describe the central tendency and the variability of the distribution, while skewness and kurtosis indicate the symmetry and the “peakedness” of the distribution. According to the imbalance situation in a frequency distribution, skewness can be negative, positive, or zero. A negative skew indicates that the distribution is left skewed (i.e., the data concentrate more on the right-hand side of the mean), a positive skew corresponds to the opposite case, and a zero skew indicates that the tails on both sides of the mean balance out, which is due to the symmetry or the even-out asymmetries. A high kurtosis corresponds to a distribution with a sharp peak and long tails, while a low kurtosis indicates a distribution with a round peak and short, thinner tails. When either skewness or kurtosis significantly deviates from the value corresponding to the normal distribution, the normality of the distribution diminishes. For a non-normal distribution, many classical statistical tests, such as t-tests, F-tests, and chi-squared, tests may not be applicable.
Figure 4 displays the maturity (θ) distributions of two cotton samples (the sample codes are 3044 and 3055), in which the x-axis represents the fiber maturity and the y-axis represents the percentage of totally detected fibers. The θs of cotton 3044 and cotton 3055 have almost the same mean (Mθ = 0.475 and 0.479), standard deviation (SDθ = 0.180 and 0.172), and kurtosis (Kθ = −0.535 and −0.532), but different skewness values (Sθ = 0.367 and 0.232). Obviously, these two cottons have very distinctive θ distributions, as shown in Figure 4, and their maturity levels should not be rated only by the mean θ values. Cotton 3044 contains more immature fibers than cotton 3055, although they have the same Mθ. Thus, distributional information must be taken into account when comparing or ranking maturities of different cottons.
θ distributions of cottons with the same Mθ, SDθ, and Kθ, but different Sθ.
Patterns of maturity distributions
As indicated in Figure 4, the maturity distributions of different cotton samples can have substantially different shapes. Among the parameters describing the shape of a distribution, skewness is the first principal component that distinguishes distribution patterns notably. Five patterns can be defined to categorize maturity distributions based on skewness values (Sθ), as listed in Table 2. From the data of the available samples, it can be observed that when Sθ is in the range of (−0.1 to 0.1), the distribution takes a shape approximate to the normal distribution (pattern III). When Sθ < −0.3 or Sθ > 0.3, the distribution is severely either right-skewed (pattern I) or left-skewed (pattern V). There is one intermediate pattern between patterns I and III, and one between patterns III and V, and each pattern indicates concentrated fiber contents in the distribution. For example, pattern I reveals that the main fiber content is “dead fiber” (extremely immature), while pattern V reflects a large number of extremely mature fibers. Figure 5 displays five maturity distributions to give one sample for each of these patterns. Only the distribution in pattern III is deemed to be approximately normal. These five patterns represent five classes of cotton maturity, which are differentiated by not only the mean values but also the other distribution parameters.
Patterns of maturity distributions. Maturity distribution patterns
Fiber types and contents
Fiber types based on θ values
Compared with mature fibers, dead fibers are more easily entangled to form neps,
10
which decrease the efficiency of yarn spinning and the uniformity of yarn structure. Neps also have an accelerated rate of desorption due to a lack of cellulosic content, and thus they are unable to retain colorants after dyeing.
11
Dead fiber is the main cause of white speck problems in a dyed fabric.12,13 White specks usually result in a major discount in price for the manufacturer. On the other hand, mature fibers offer better strength and dyeability, and should be at a premium in terms of the ultimate assessment of cotton maturity. Therefore, the dead fiber content (DFC), immature fiber content (IFC), and MFC in cotton are the important indicators of fiber quality. The three fiber contents can be obtained by adding the relative probabilities within the three separate θ ranges, as illustrated in Figure 6.
Contents of dead, immature, and mature fibers in a maturity distribution.
Results and discussion
In this paper, the tests performed with the previous and current FIAS analyses were denoted as PA and CA, respectively. Because the previous FIAS was not able to identify correctly many of the fiber cross-sections that had unclear or incomplete boundaries, a manual editing function was added to allow the user to remove mistaken cross-sections from the processed image on the screen. When the mouse was double-clicked on a fiber, the cross-section could be marked, and the corresponding data were deleted from the output. The PA test after the manual removal was symbolized as PAM data, which were provided by FBRI. The current FIAS improves the fiber-detection functions and can automatically delete wrong cross-sections even if they are initially detected. It also has a set of manual editing tools to assist the user in amending unclear boundaries, detaching connected fibers, and deleting anomalous objects using a mouse or a tablet pen. The manual editing can be repeated until the visual examination on all the fibers is satisfied, and the manual editing process after the CA test is called CAM, the results of which can be regarded as the most trustful data for assessing the accuracy of the CA data.
The 2201 images of cotton variety 2996 were selected to perform the CA and CAM tests, and these images were first batch-processed by the current FIAS and then manually edited image by image. Figure 7 displays the θ histograms of the CA and CAM tests and the descriptive statistics (Mθ,, SDθ, Sθ, and Kθ). Although slight discrepancies in fiber frequency can be observed near θ = 0.3, the agreement of the two histograms is extremely high (R2 = 0.999), and the two sets of the descriptive statistics are approximately identical. This proves that the CA test is able to produce reliable fiber detections, except for irreparable immature cross-sections. Because of the high agreement between the CA and CAM data, there is no need to perform manual editing after the CA test.
θ histograms and the descriptive statistics of cotton 2996 in the CA and CAM tests.
Comparisons in number of the identified fibers
Comparisons in mean maturity (Mθ)
From a linear regression.
From a quadratic regression.
AFIS: Advanced Fiber Information System.
Distribution parameters of maturity
Let us take the θ histograms of cotton 2999 in the PA, PAM, and CA tests as an example (see Figure 8). Although the Mθ values in these tests are similar (0.48, 0.45, and 0.45, respectively), their histograms are distinctively different. The CA histogram is severely right-skewed with a sharp peak in the dead fiber range (θ < 0.3), while both the PA and PAM histograms are only slightly right-skewed. Compared with the other two histograms, the PAM histogram has much lower fiber counts because of the removal of wrongly detected fibers.
θ histograms of cotton 2999 in the PA, PAM, and CA tests.
The maturity distributions of the seven cottons in the CA test are displayed in Figure 9, showing diversified distributions from an approximately normal curve (cotton 3074) to a heavily right-skewed curve (cotton 3075). The shape difference of these two cottons are reflected by their Sθ values, which are −0.03 for cotton 3074 and 0.61 for cotton 3075. Of the seven cotton samples, there are two pairs that have similar distributions: cotton 2996/cotton 3009 and cotton 2999/cotton 3008. In each pair, the distribution descriptors Mθ, SDθ, Sθ, and Kθ of the two samples are greatly comparable. The distribution of cotton 3016 has the flattest top among this group, corresponding to its large negative Kθ (−0.60). When assessing the maturity levels of different samples, one must include all of the four descriptors.
θ distributions of the seven cottons in the CA test.
Contents of dead, immature, and mature fibers
R is the correlation coefficient between the fiber contents in this table and the Mθ values of the CA test in Table 5.
Conclusions
The algorithm changes made in the current FIAS improved the consistency in detecting fibers of various maturities and effectively reduced the bias on immature fibers in previous maturity distributions. The maturity distribution generated in the CA test was correlated highly with that of the manually traced fibers (R2 = 0.999). The dramatic reduction in the number of wrong fiber detections eliminates the need for manual editing in a routine test.
The maturity level of a large quantity of cotton fibers should be evaluated with multiple distributional parameters, because any single parameter is insufficient in characterizing the maturity and is often invalid when the distribution is not normal. The basic shapes of cotton maturity distributions can be categorized into five patterns, each representing a major class of cotton maturity. In the same class, the maturities of different samples can be ranked. The skewness of a maturity distribution is an essential parameter to classify the distribution pattern. Most cotton varieties exhibit positively skewed distributions, indicating more IFCs than mature ones.
From a maturity distribution, the contents of dead (DFC), immature (IFC), and mature (MFC) fibers can be calculated by adding the relative probabilities in the specific θ ranges. It seems that IFC is the major body of the maturity distributions for most cotton samples, but only the DFC and the MFC highly correlate with the Mθ—the overall maturity level. Since dead fibers lack strength and dyeability, the DFC should be a discount factor in maturity evaluation. On the other hand, the MFC should be used as a premium in the ultimate rating of fiber quality.
Footnotes
Funding
This work was supported by Cotton Inc. (grant number 13-585).
Acknowledgment
The authors would like to thank Dr Eric Hequet at the FBRI of Texas Tech University for the images and the tested data of the seven cotton samples used in this research.
