Abstract
In this article, a method based on gray projection is proposed for automatic measurement of surface braiding angle and pitch length of braided composite preforms. The surface braiding angles are measured by rotated gray projection. First, the original image is preprocessed using Lab transform and block-matching and three-dimensional filter. Then, edge map of gray scale is acquired based on phase congruency and non-maximum suppression. Third, edge direction angles are computed by image rotation and gray projection. Finally, the average surface braiding angles are measured from the edge direction angles. The pitch lengths are measured by cropped gray projection. First, the original image is filtered by an artistic edge and corner enhancing filter. Second, block images are cropped based on gray projection. Finally, pitch lengths are measured by block images gray projection. Experimental results show that the proposed method can achieve the automatic measurement of average surface braiding angle and average pitch length of two-dimensional braided composite preform with high accuracy.
Introduction
As representatives of lightweight materials with high stiffness and strength, fiber-braided composite materials have been widely used in fields of aerospace, automotive, military, energy, sports, and so on.1,2 The braid pattern evolves from an interaction of the yarn movement, which is interlacing three or more threads diagonally to the product axis in order to obtain a thicker, wider, or stronger product or to cover some profile,3,4 forming braided preforms, then the preforms are consolidated by resin transfer molding (RTM) 5 or vacuum assisted resin transfer molding (VARTM),6,7 as composites.
Generally, the braided structures include three-dimensional (3D) braid (see Figure 1(a)) and two-dimensional (2D) braid (see Figure 1(b) to (c)) according to how the yarns interlace into 3D space or 2D plane. 8 The 2D braided structures can be classified as biaxial braid (see Figure 1(b)) and triaxial braid (see Figure 1(c)) according to whether there are axial yarns or not. 9 Biaxial braids is very basic and most common form of 2D braids, which are composed of interlacing yarns arranged in two directions. 10 According to the braiding patterns, biaxial braids are classified as diamond, regular, and Hercules 11 (see Figure 2).

Structures of 3D and 2D braids: (a) 3D, five directional braid; (b) 2D biaxial regular braid; and (c) 2D triaxial braid.

Structures of 2D braids: (a) diamond braid, (b) regular braid, and (c) Hercules braid.
The surface braiding angle θ is defined as the angle between the longitudinal direction of the braided preform and the deposited fibers, 12 and pitch length h is the preform take-up length for one machine cycle,13,14 as shown in Figure 3. Surface braiding angle and pitch length are two important parameters to affect the kinematic analysis of the process, geometry of braided composite.9,14 For example, the surface braiding angle θ can be calculated by Heieck et al., 9 Potluri et al. 15
where R is the radius of cylindrical mandrel (cm), v is take-up speed (cm/s), and ω is the average angular velocity of bobbins around the machine center (rad/s). Sun and Sun 16 indicated that braiding angle went down with increased pitch length. The surface braiding angle and pitch length are also important parameter of influencing the mechanical behavior of braided composite,15,17–20 for instance, as the surface braiding angle increase, the mismatch between longitudinal and transverse stiffness become less pronounced, with a tendency toward balanced Poisson’s ratios for the [0°/± 60°] material. 21 Zou et al. 22 indicated that the radial force is higher with shorter pitch length when surface braiding angle is fixed.

Surface braiding angle and pitch length: (a) in an ideal scenario and (b) in a practical scenario.
In industrial, manual measurement of average surface braiding angle is still the main method, which is time- and manpower-consuming.
At present, several automatic measurement methods23,24 have been proposed for 3D braided composites preforms. For the parameter measurement of 2D braided preforms, Hunt and Carey 25 applied a frequency domain method to achieve automatic real-time of 2D braided preform. This method is more suitable for the fabric images with high contrast between fiber bundle and background and is easy to be affected by the contrast of braided axis. Few methods are proposed for pitch length measurement of 2D braided composite preform. Thus, it is significant to achieve the automatic measurement of the two parameters. In this article, we propose a new automatic method based on gray projection to measure surface braiding angle and pitch length of 2D biaxial regular braided composite preforms.
Materials and experimental methods
Materials
The tested 2D biaxial regular braided tubes were supplied by Shuomin Technology Co. Ltd, China. The braided fibers include carbon, glass, and aramid fiber, and the average diameters of the fiber filaments are, respectively, 5.70, 7.91, and 11.02 µm, tested by the School of Textiles of Tianjin Polytechnic University using the fiber fineness analyzer manufactured by Laizhou Electron Instrument Co., Ltd. (Taiyuan, China). Since the braided tubes are easy to stretch, they are pasted on horizontal plates using double-sided tapes. And the double-sided tapes pasted on the back of the fabrics cover all the areas of the fabrics to make sure that the fiber tows of the fabrics can be well fixed on the plane, as illustrated in Figure 4 (acquired by Redmi 4X mobile phone).

Fabric sample: (a) glass-fiber fabric, (b) aramid-fiber fabric, and (c) carbon-fiber fabric.
Image acquisition
The original images are acquired using the image acquisition system shown in Figure 5. First, the CCD camera (Vieworks VA-8MC-M16AO; 3296 × 2472 pixels; Area CCD camera) with camera lens (Nikkor 28mm f/2.8D) and circular polarizing filter (CPF) (NiSiCPL52mm) attached is vertically installed above the fabric, and the micrometer is put on one side of the fabric for calibration, keeping the surfaces of the micrometer and the fabric on the same horizontal plane. Then the dome light source (DLS) (OPT-RID180-RGB) is put in front of the fabric, making the camera’s optical center go through the central axis of the DLS light hole and also the center of the fabric surface. Finally, the camera exposure value and focal length, and the rotated angle of CPF are adjusted, and a clear image is selected.

Image acquisition system: (a) system diagram and (b) physical system.
During the image acquisition, to make sure that the CCD camera is vertically installed above the fabric, the gradienter is stuck to the shelf by its own magnet (see Figure 5(b)), then it is adjusted to be horizontal according to the gradienter reading. To test the accuracy of the acquired images, we compare the image corrected by calibration board (GQ-H56–3.5; expansion coefficient: α = 8.6 × 10−6/°C) (see Figure 6) and HALCON software. 26 The process of distortion correction is described as follows: first, obtain the distortion parameters of the camera by finding the marked points centers of the calibration board images (see Figure 7); then correct the image based on the distortion parameters.

The 7 × 7 HALCON calibration board.

Calibration board images: (a) I1, (b) I2, (c) I3, (d) I4, (e) I5, (f) I6, (g) I7, (h) I8, (i) I9, (j) I10, (k) I11, (l) I12, (m) I13, (n) I14, (o) I15, (p) I16, (q) I17, (r) I18, (s) I19, (t) I20, (u) I21, (v) I22, (w) I23, (x) I24, (y) I25, and (z) I26.
Figure 8 is image before and after distortion correction. It can be seen that there is almost no difference between the image before and after distortion correction. So the images acquired by the system have high accuracy.

Before and after correlation: (a) before correlation, (b) after correlation, (c) local amplification of (a), and (d) local amplification of (b).
Figure 9 is the images acquired based on “CCD camera + nature light source” and based on the proposed system. It can be seen that images acquired based on “CCD camera + nature light source” have reflections with varying degrees (see Figure 9(a), (c), and (e)), and the surface reflections on Figure 9(b), (d), and (f) have been weakened efficiently, and also the images (especially carbon-fiber fabrics) have higher contrast in texture structure. The reasons are described as follows: first, the DLS can make the light uniformly illuminate on the fabric surface; second, the CPF has the characteristic of filtering reflection in a certain direction.27,28

Images acquired by common acquisition system and by our system: (a) carbon-fiber fabric acquired by “CCD camera + nature light source,” (b) carbon-fiber fabric acquired by our acquisition system, (c) glass-fiber fabric acquired by “CCD camera + nature light source,” (d) glass-fiber fabric acquired by our acquisition system, (e) aramid-fiber fabric acquired by “CCD camera + nature light source,” and (f) aramid-fiber fabric acquired by our acquisition system.
Surface braiding angle measurement
Before surface braiding angle measurement, the original acquired image Fori is first preprocessed by Lab transform 29 and block-matching and 3D filter (BM3D) 30 to enhance contrast and filter noise, respectively.
The CIE-Lab color space introduced by the Commission Internationale de L’Eclairage (CIE) in 1976 has three axes: L, a, and b. It first converts from R, G, and B space into X, Y, and Z space, then converts X, Y, and Z into L, a, and b space, using equations (2) to (4) as follows 29
We use L-channel image FL to achieve image enhancement. Then the L-channel image FL is filtered by BM3D.
BM3D achieves image filter in 3D transform domain by combining sliding-window transform processing with block-matching. First, divide the image into small blocks and choose one block as currently processed one, then search blocks which are similar to the currently processed one. Second, stack the matched blocks together, forming a 3D array, which have high level of correlation. Third, exploit this correlation by applying a 3D decorrelating unitary transform and attenuate the noise by shrinkage of the transform coefficients, and subsequent inverse 3D transform yields estimates of all matched blocks. Finally, repeat this procedure for all image blocks in sliding manner, and compute the final estimate as weighed average of all overlapping block-estimates.
Figure 10 is the process of image preprocessing, where Figure 10(a) is original image acquired by image acquisition system, and Figure 10(b) and (c) are the L-channel image and BM3D filtered image, respectively. It can be seen that original image preprocessed by Lab transform and BM3D filter has higher contrast and weaker fiber texture.

Original image and the preprocessing image: (a) original image Fori, (b) L-channel image FL, and (c) BM3D filtered image FBM3D.
After preprocessing, image FBM3D is then processed for average surface braiding angle measurement based on rotated gray projection. The detailed steps are described as follows:
Step 1. Calculate phase congruency maps 31 of images FBM3D (see Figure 11(a)), getting images FPC (see Figure 11(b)).
Step 2. Calculate non-maximum suppression maps 32 of FPC, getting images Fnom (see Figure 11(c)).

Process of edge direction angles detection: (a) FBM3D, (b) FPC, (c) Fnom, (d) rotated Fnom at −36.1°, (e) gray projection curve of (d), (f) rotated Fnom at 48.8°, and (g) gray projection curve of (f).
From Figure 11(c), it can be seen that the image edges can be classified into two directions, labeled by blue ribbon arrow lines. So the average surface braiding angle θaBM3D can be calculated by the angles θr+ and θr– between the edges of two directions and the horizontal direction (braiding axis direction)
And the angles θr+ and θr– can be calculated based on image rotation, described in the following steps.
Step 3. Calculate edge direction angles θr+ and θr– of image Fnom based on the principle that maximum peak value of gray projection is largest when the edge direction is perpendicular to the projection plane.
Step 3.1. Rotate image Fnom from −89° to 89° with step length 0.1 (the precision of angle can be accurate to 0.1°), getting rotated images Froti (i = 1, 2, . . ., 1781). As shown in Figure 11(d), where the rx–ry are coordinate axes of image rotation, and o located in the center of the image is origin of coordinate zero.
Step 3.2. Calculate gray projection images of Froti projected on the projection plane γ (see Figure 11(d)) that is perpendicular to the image plane, getting image Fgpi.
Step 3.3. Compute the maximum peak values Pi in each Fgpi. Then obtain the angle corresponding to the maximum value Pmax+ of Pi when the rotated angles are positive, as θr+, and the angle corresponding to the maximum value Pmax– of Pi when the rotated angles are negative, as θr–.
As shown in Figure 11(d) to (g), for example, in Figure 11(d) and (e), Pmax– = 24.44, θr– = –36.1°, and the edge direction of direction 1 is perpendicular to the projection plane γ, and in Figure 11(f) and (g), Pmax+ = 25.28, θr+ = 48.8°, and the edge direction of direction 2 is perpendicular to the projection plane γ.
Step 4. Compute average surface braiding angle.
The average surface braiding angle of image FBM3D is calculated by equations (5).
As shown in Figure 11(d) to (g), the average surface braiding angle θaBM3D = (48.8 + 36.1)/2 = 42.45°.
Step 5. Crop block images FCek, FRdk, FLdk, FRuk, and FLuk (k = 1, 2) from image FBM3D, where the subscripts Ce, Rd, Ld, Ru, and Lu for F are used to indicate the block images cropped by starting from center, right down, left down, right up, and left up of FBM3D. FCek, FRdk, FLdk, FRuk, and FLuk can be calculated by equations (6) to (10)
where m and n are, respectively, the row number and column number of image FBM3D. As shown in Figure 12, for intuitive show, the cropped images are superimposed on the original images, from small to large.

Cropped block images: (a) FCek, (b) FRdk, (c) FLdk, (d) FRuk, and (e) FLuk.
Step 6. Repeat step 1–4 to detect the average surface braiding angles θaCek, θaRdk, θaLdk, θaRuk, and θaLuk of block images FCek, FRdk, FLdk, FRuk, and FLuk (k = 1, 2), respectively, and use the average value θa as the final average surface braiding angle, as computed by
Pitch length measurement
Before pitch length measurement, the original acquired image Fori is first preprocessed by an artistic edge and corner enhancing filter (AECE) 33 to filter the effects of fiber texture and noise.
AECE is a noniterative edge and corner preserving smoother. As shown in Figure 13, the original image Fori (see Figure 13(a)) filtered by AECE has less noise and fiber texture (see Figure 13(b)), which is better for the flowing image processing.

Original image and filtered image: (a) original image Fori and (b) AECE filtered image Fae.
Then the AECE filtered image Fae is processed for average pitch length measurement based on cropped gray projection. The detailed steps are described as follows:
Step 1. Project image Fae on plane γ, getting gray projection curve Fpo. As shown in Figure 14(a), where the plane γ is perpendicular to image Fae, then smooth Fpo, getting Fp (see Figure 14(b)).
Step 2. Crop block images from Fae based on projection image Fp. First, in Fp, calculate the locations dpti (i = 1, 2, . . .) of peak and trough. Then cut image Fae by using dpti as the row border, getting block images Fbi (see Figure 14(c)).
Step 3. Project block images Fbi to the plane β, getting projection curves, and smooth these curves, getting Fbpk (k = 1, 2, . . .). The plane β is shown in Figure 14(a) and is perpendicular to block images Fbi. Then, calculate the peak-to-peak values and trough-to-trough values in Fbpk (k = 1, 2, . . .), getting dplpj (j = 1, 2, . . .). Then obtain the average values

Process of average pitch length measurement: (a) image Fae and projection planes, (b) Fp, (c) block images Fbi, and (d) Fbp1.
Figure 14(d) is projection image of one block image Fb1, for example.
The principle of this step is described as follows: first, in Figure 14(c), the clock image Fb1, for example, the distance between Ai and Aj (labeled by yellow points) or distance between Bi and Bj (j = i + 1) (labeled by blue points) (i = 1, 2, 3, . . .), is equal to the pitch length in pixels. Second, as shown in Figure 14(d), for example, the peak and trough values of Fbp1 approximately reflect points Ai and Bi of Fb1 (see Figure 14(c)), respectively. Thus, the single pitch length in pixels can be obtained by calculating the peak-to-peak values and trough-to-trough values of projection image Fbp1. Similarly, single pitch lengths in pixels in other block images can be calculated, and the average pitch length in pixels can be obtained based on all the single pitch lengths.
Step 4. Calculate average pitch length by
where dc defined as calibration value represents the actual distance per pixel.
Results and discussion
In this section, five types of images are tested using the proposed method, and each type contains 35 images, as shown in Figure 15, where Figure 15(a) to (e) are, respectively, glass-fiber, aramid-fiber, carbon (6k), carbon-fiber (3k), and carbon-fiber (3k) 2D biaxial regular braided images.

Sample images: (a) F1 (glass-fiber), (b) F2 (aramid-fiber), (c) F3 (carber fiber (6k)), (d) F4 (carber fiber (3k)), and (e) F5 (carber fiber (3k)).
Manual measurements of average surface braiding angle and pitch length are used to benchmark our proposed method. Manual measurements of average surface braiding angle is based on electronic protractor, as demonstrated in Figure 16(a), for example, one angle θ is measured to be 62.85°, and thus the single surface braiding angle θM is set to θM = θ/2 = 62.85/2≈31.43°. The average surface braiding angle θaM is the average value of all single surface braiding angles, calculated by
where n is the number of single surface braiding angles measured. Figure 16(b) is the labeled single angles θi (i = 1, 2, 3, . . ., 35) to be measured, and here n = 35. The manual measurement of pitch length is based on clicking corner points of fabric image on computer screen, as displayed in Figure 16(c), for example, the coordinates of clicked points are A(97, 192) and B(99, 248), respectively. So one pitch length is calculated as dplAB = dc × dAB, where dc is calibration value, that is, the true distance per pixel, and
where k is the number of single pitch lengths measured. Figure 16(d) is the labeled corner points to be measured, here we only labeled the first two lines’ corners, for example, pitch lengths dplAiAj (i = 1, 2, 3, . . .; j = i + 1) can be calculated as dplAiAj = dc × dAiAj, where dAiAj is distance between corners Ai and Aj. Similarly, pitch lengths dplBiBj and pitch lengths in other lines can be defined and calculated. So the average pitch length of Figure 16(d) can be computed by equation (14).

Manual measurement of surface parameters: (a) manual measurement of surface braiding angles, (b) labeled angles, (c) manual measurement of pitch lengths, and (d) labeled corners.
The errors of the average surface braiding angle and pitch length measured by our method relative to the manual measurements are illustrated in Figure 17, where the errors are sorted in ascending order. Tables 1 and 2 show the average values, relative errors, standard deviations, and average relative errors of average surface braiding angle and pitch length, respectively, where
where in equations (17), N = 35, ei is each error of average surface braided angles (or average pitch lengths) measured by our method relative to manual measurement, and
Average values, relative errors, standard deviation, and average relative errors of average surface braiding angle.
Average values, relative errors, standard deviation, and average relative errors of average pitch length.

Relative error curves of average surface braiding angles and pitch lengths: (a) average surface braiding angle relative errors of image F1; (b) average pitch length relative errors of image F1; (c) average surface braiding angle relative errors of image F2; (d) average pitch length relative errors of image F2; (e) average surface braiding angle relative errors of image F3; (f) average pitch length relative errors of image F3; (g) average surface braiding angle relative errors of image F4; (h) average pitch length relative errors of image F4; (i) average surface braiding angle relative errors of image F5; and (j) average pitch length relative errors of image F5.
The following observations are made from Figure 17, Table 1, and Table 2:
Standard deviations of relative errors, and average values of relative errors of average surface braiding angle and pitch length are small, which shows that the stability and accuracy of the results measured by our method.
In the measurement of average surface braiding angle, sample image F4 has smaller relative errors and standard deviation than other samples, for the edges in sample F4 are more closed to straight line than other types of samples, which has smaller impact in gray projection.
The relative errors of average pitch length are smaller than the values of surface braiding angles, for the average pitch lengths are measured based on local features of fabric image while average surface braiding angles are measured based on the overall features of fabric image.
Conclusion
This article presented a new method for automatically measuring the average surface braiding angles and pitch lengths for biaxial regular braided glass-, aramid- and carbon-fiber preforms. Some conclusions are drawn as follows:
The proposed method can not only measure surface braiding angles but also measure pitch lengths.
The phase congruency used in initial edge detection improves the accuracy of edge detection, so as to improve the accuracy of average surface braiding angle measurement.
In the calculation of average surface braiding angle, the measurement results depend on the average surface braiding angles of multiple block images, which is helpful to reduce the errors of measuring only one image.
The cropping method based on overall image projection well crop image into block images according image structure texture, which is helpful to extract local pitch length of fabric.
The proposed automatic method can be used as a foundation of computer parameter measurement for 2D braided texture analysis.
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by Applied Basic Research Programs of China National Textile and Apparel Council (No. J201509) and the Program for Innovative Research Team in University of Tianjin (No. TD13-5034).
