Abstract

In the last decade, a significant number of techniques for three-dimensional (3D) shape measurement have been proposed. There are a large number of measurement demands for metallic workpieces with shiny surfaces in industrial applications; however, such shiny surfaces cannot be directly measured using the conventional structured light method. Therefore, various techniques have been investigated to solve this problem over the last few years. Some reviews summarize the different 3D imaging techniques; however, no comprehensive review exists that provides an insight into high-dynamic range (HDR) 3D shape measurement techniques used for shiny surfaces. We present a survey of recent HDR techniques for the digitization of shiny surfaces and classify and discuss the advantages and drawbacks of different techniques with respect to each other.

1. Introduction

Three-dimensional (3D) shape measurement is an important topic that has been investigated to a large extent. A significant number of techniques for acquiring 3D shapes have been proposed, leading to a wide range of applications, such as range sensing, industrial inspection of manufactured parts, reverse engineering (digitization of complex, free-form surfaces), object recognition, entertainment, 3D map building, biometrics, and documentation of cultural artifacts [1, 2]. The recent release of Microsoft Kinect has enabled 3D technology to penetrate into our daily lives [3], and the combination of 3D printing technology with 3D shape measurement techniques could generate additional excitement across disciplines. These techniques cover a wide range from microscale objects, such as interferometry and confocal techniques, to several meter-sized objects, such as time-of-flight technique [4], and generally can be classified into two categories: surface contact and surface noncontact. The surface contact measurement system, such as mechanical probe-based coordinate measuring machine (CMM), has been used in the industry for many years and can measure geometric shape with high accuracy (usually more accurate than current noncontact systems). However, CMM is relatively slow and only measures a limited number of points in an object’s surface. In contrast, the optical methods have an ability to provide a high-speed 3D inspection of objects in a noncontact manner. These methods include time-of-flight, stereovision, photogrammetry, laser range scanning, and structured light methods. Among these, the structured light method is one of the most widely used 3D shape measurement techniques owing to the advantages of full-field inspection, high resolution, and accuracy that it offers [5].

Typically, a structured light measurement system consists of a camera and a projector that projects a sequence of coded patterns onto object surfaces. The camera captures images of the reflected pattern deformed by the object surface from another direction. The reflected pattern has the height information of the object embedded in its phase distribution. Then, from the captured images, a specially designed algorithm, such as a phase-shifting algorithm, is employed to obtain the absolute phase map, which can be used for determining correspondences between the camera and the projector in a structured light system. Finally, the 3D coordinates of the object surface can be calculated from the absolute phase map using a phase-to-height mapping based on a calibration performed earlier over the measurement volume of interest. The measured point cloud of the object can be further compared with the computer aided design (CAD) model to evaluate the manufacturing quality. However, the use of this system still involves difficulties when dealing with nondiffuse (non-Lambertian) surfaces because the optical signal cannot be properly retrieved. Most of the optical methods assume that the surface being measured is a Lambertian surface, which has a diffuse or near-diffuse reflectance and low reflectivity variations from point to point. Therefore, it is usually very difficult for any optical method to accurately measure shiny objects or objects with a large range of reflectivity variation across the surface. For example, Figures 1(a) and 1(b) show a plastic block and the corresponding captured fringe pattern image, respectively. The camera can capture a clear image of the fringe pattern owing to the diffuse surface of the plastic block. However, for a different surface material, such as a metallic workpiece shown in Figure 2(a), we can see in Figure 2(b) that only regions with strong reflection can provide measurable fringe patterns, while other regions cannot been imaged by the camera. As we increase the exposure time, in Figure 2(c), fringe patterns reflected from regions with low reflection become visible while regions with strong reflection become so bright that image saturation occurs. Thus, insufficient information is available to obtain the 3D coordinates of the metallic workpiece. The fundamental reason is that the metallic workpiece has a much higher level of intensity variation in their reflections compared with what can be captured by the 0–255 gray-level intensity range of the camera.

A large number of metallic workpieces are needed to be milled in industrial applications. After milling, the workpiece surface will become very bright that the conventional structured light method cannot directly measure. One method to solve this problem is to use the CMM to measure the workpiece, and the other is to spray a thin layer of powder onto the workpiece to make its surface diffuse prior to measurement by the structured light method. However, this supplementary step is troublesome and time-consuming because the workpiece needs to be cleaned afterward. In addition, the final accuracy is often dependent on the powder thickness and its homogeneity [6]. To avoid the spraying step, various methods have been investigated over the last few years. State-of-the-art methods for 3D shiny surfaces measurement are typically referred to as high-dynamic range (HDR) techniques. HDR 3D shape measurement techniques can be classified into six categories: (1) using multiple exposures; (2) adjusting projected pattern intensities; (3) using polarizing filters; (4) using color invariants; (5) photometric stereo technique; and (6) miscellaneous. Some reviews detail the different techniques for 3D shape measurement [5, 7, 8]. However, no comprehensive review exists that provides an insight into HDR techniques used for shiny surfaces. In this paper, we review various HDR techniques that were recently developed for the measurement of shiny surfaces, discuss the advantages and drawbacks of different techniques with respect to each other, and assess the current state-of-the-art techniques in the field.

The remainder of this paper is organized as follows. Section 2 attempts to classify and discuss the advantages and drawbacks of various HDR techniques. Section 3 compares these HDR techniques based on their main characteristics, and Section 4 draws a conclusion.

2. HDR Techniques

2.1. HDR Techniques Using Multiple Exposures

In 2009, Zhang and Yau [9] proposed a technique called “high-dynamics range scanning” using multiple exposures, where a sequence of images captured at different exposures is combined into a single set of HDR images, that is, phase-shifting images, by selecting the brightest unsaturated intensities at each pixel. Typically, the best method to adjust the exposure to acquire high-quality fringe patterns is to manually adjust the exposure time instead of the lens aperture because manual adjustments of the lens aperture can cause undesired motion between system components that changes the system calibration. The images taken at low exposure times contain useful information from regions on the surface with high reflectivity, whereas the images taken at high exposure times contain useful information from regions on the surface of low reflectivity. Therefore, the specular or dark region can be properly measured without affecting the rest of the regions. Figure 3 shows the synthesis process. This technique is less sensitive to local surface reflectivity variation, benefiting from the phase-shifting algorithm, which is less sensitive to the ambient light. However, for a large range of reflectivity variation, the signal-to-noise ratio (SNR) is very small for low reflectivity regions. Thus, it is necessary to control the ambient light at a very low level in the industrial in situ environment; otherwise, the quality of measurement is hard to be ensured. Since the used exposure time is subjectively selected, it lacks quantitative manner to determine the proper exposure time. Besides, to synthesize the HDR images, a large number of exposure times are required in some cases, which is time-consuming. For example, in Zhang’s experiment, they needed 23 exposures to obtain a high-quality point cloud of a China vase. Even with the three-step phase-shifting algorithm and three-frequency heterodyne phase-unwrapping algorithm, a total of fringe images are required to synthesize the HDR images.

In addition to the use of multiple exposures, Liu et al. [10] used a dual-camera structured light system to acquire the 3D shape of shiny objects and fill in missing data from regions with specular highlights and diffused darkness, that is, two camera-projector monocular systems operated from different viewing angles. An image-segmentation technique was also introduced to distinguish which monocular system is suitable to reconstruct a certainly lost point accurately. This approach worked well because specular highlights are highly dependent on lighting and viewing angles, although additional complexity is required in processing for image masking and registration to merge different surface regions.

For the use of multiple exposures, the initial values of the exposure time are difficult to determine when inspecting an unknown scenario because the initial values cannot be obtained adaptively according to the measured scene. Thus, to adapt different surface reflectivity, it is necessary to use the exposure time as many as possible or rely on experience to select several exposure times. To avoid these, Feng et al. [11] used the histogram to predict one or more optimal exposure times based on the reflectivity distribution of the surface. Then, raw fringe images were captured using these optimal exposure times to synthesize the HDR images.

Due to the influence of camera blurring and random noise, the conventional intensity-based pixel saturation identification method is sometimes susceptible, which determines the brightest intensity reaching the highest quantization level, for example, 255 for an 8-bit camera. Thus, the brightest intensity recorded in a pixel can be unstable and deviates from 255, especially when the saturation of the pixel is slight. Long et al. [12] used the magnitude of a nonprincipal frequency component to identify saturated pixels, leading to more accurate and robust saturation identification and higher measurement accuracy.

Alternatively, Jiang et al. [13] proposed a technique for capturing raw fringe images by simultaneously adjusting the exposure time and the projected light intensity. This technique takes the highest modulation intensity rather than the brightest unsaturated intensity as the criterion for selecting the pixels from the raw fringe images to synthesize the HDR images so that a high SNR is obtained and the ambient light effect is minimized. In addition, a method of autoselection of parameters of the exposure time and the projected light intensity is developed. However, its implementation is relatively complicated, and it achieves the measurement by adding at least 5 times the number of fringe images compared with the conventional phase measurement. To address this time-consuming problem, the same team [14] proposed a very fast HDR fringe pattern projector with a frame rate of 700 Hz, which reduced the routine projecting time by 88%. The dynamic range of the light intensity can vary at least 9 times by using multiple homochromous LED chips and intensity modulation technique.

Similarly, for this time-consuming problem, Ekstrand and Zhang [15] presented an autoexposure technique in which the required exposure time can be predicted automatically according to the reflectivity of the measured object surface. This technique reduces the human intervention and improves the intelligence of the 3D shape measurement system. However, the single predicted exposure time does not always fit the surface with a large range of reflectivity variation, for example, the surface on the metallic workpiece shown in Figure 2. In 2014, a similar approach was presented by Zhong et al. [16], which selects an optimal exposure time as large as possible and meanwhile avoids image saturation in regions of specular reflection during the process of phase shifting. However, these techniques using a single exposure time are limited to improve the SNR in the region of weak reflection.

The advantage of these techniques is that, besides some postprocessing for the captured raw fringe images, there is no need to make too many changes based on the original technique. Adjusting the exposure time would appear to be the best option. However, the exposure time cannot be arbitrarily chosen for a structured light system, where a digital-light-processing (DLP) projector is used to project the fringe patterns. To precisely capture the projected grayscale values, the camera and the projector must be precisely synchronized, and the exposure time must be set as an integer of 1/fp s, where fp is the refresh frequency of the projector (fp is generally 60 Hz).

2.2. HDR Techniques by Adjusting Projected Pattern Intensities

Adjusting the projected fringe pattern intensities is another option. Waddington and Kofman [1719] proposed a technique of projecting sinusoidal fringe patterns with the modified maximum input gray level (MIGL) to accommodate variable ambient illumination that would otherwise cause image saturation and measurement error in 3D shape measurement. However, this technique that only uniformly adjusts the MIGL to the projector to avoid image saturation will reduce the SNR when measuring those surface regions with low reflectivity. Then, the same team [20, 21] combined the MIGL reduction and pixel-by-pixel approaches, that is, by uniformly adjusting the projected fringe pattern intensity and selecting the brightest unsaturated pixels from raw fringe images captured at different illuminations to synthesize the HDR images. This technique can maintain a high SNR for the measurement of the surface with low reflectivity. However, similar to the techniques using multiple exposures, it also could be time-consuming. Recently, Li and Kofman [22] developed an adaptive fringe pattern projection technique by locally modifying the MIGL in projected fringe patterns, according to the local reflectivity of an object surface being measured, which is better than uniform reduction of the MIGL for saturated regions.

Babaie et al. [23] proposed a new technique for measuring 3D shape of objects with a large range of reflectivity variation. Based on the feedback from the reflected images captured by the camera, it recursively adjusts the intensities of the projected fringe pattern at pixel level. Then, the HDR images are acquired, and a four-step phase-shifting algorithm is used to obtain the absolute phase map of the object.

In 2016, Lin et al. [24] developed an adaptive digital fringe projection technique for HDR 3D shape measurement. This technique can handle the shiny surfaces and achieve improved measurement accuracy by adaptively adjusting the pixel-wise intensity of the projected fringe patterns. The adjustment is based on the camera response function according to the reflectivity of the surface and to the illumination of ambient light and surface interreflections. As a novel mathematical model was established to calculate the intensity of the ambient light and surface interreflections, the optimal intensity of each pixel in the fringe pattern to be projected is precisely calculated, which is superior to the previously proposed techniques using raw fringe images captured at different illuminations to synthesize the HDR images. In addition, this technique only needs to capture 31 images to obtain the point cloud of an aluminum alloy workpiece in their experiment. The comparison of results after using this technique is indicated in Figure 4.

Different from the aforementioned techniques using fringe pattern, Zhang et al. [25] utilized monochromatic white and black stripe patterns, leading to increased robustness. An adaptive intensity mask was adopted to dynamically adjust the pattern intensities to prevent overexposure in shiny regions, which is derived from point spread function and camera response function. The point spread function in this method is based on the homography matrix from the camera image plane to the projector image plane, which can be calibrated using the measurement table in advance. However, during the measurement, the position of the object being measured is difficult to ensure that it is the same as that of the checkerboard, so the coordinates after mapping by the homography matrix may not be very accurate.

An advantage of adjusting the projected pattern intensities is that the adjustment can be performed locally pixel by pixel, rather than that of the adjustment of the camera exposure time, which is global. To minimize the influence of the ambient light and surface interreflections, the lens aperture should be minimized and the used projector should provide intense light projection when applying these type of HDR techniques in the measurement.

2.3. HDR Techniques Using Polarizing Filters

To handle objects with shiny surfaces, an alternative method uses polarizing filters to eliminate the effect of highlights by changing the angle between the transmission axes of the polarizers. The basic idea is to limit the reflected light incident on the image sensor at a certain angle. As a result, the intense specular light is effectively removed; however, at the expense of reduction of captured intensity for the whole image, it may result in low SNR. The method introduced in [26, 27] used polarizing filters, which drastically reduce the output light intensity of the projector and the incident light of the camera; therefore, the dark regions of the scene are difficult to measure. On the contrary, a parallel polarizer-analyzer alignment [28, 29] can maintain good fringe quality at the dark regions, but not in the bright places. Additionally, Feng et al. [11] reported a general solution for realizing HDR 3D shape measurement based on digital fringe projection. They combined the techniques using multiple exposures and the techniques using polarizing filters. For the method by polarization, the reduced intensity by polarizers may result in low SNR for captured fringe images. For the method using multiple exposures, they may need to capture many raw fringe images to synthesize the HDR images.

In 2014, Salahieh et al. [30] proposed a multipolarization fringe projection imaging system, shown in Figure 5, in which saturated points are eliminated and the fringe contrast is enhanced by selecting the properly polarized channel measurements. The projected fringes are linearly polarized prior to incidence on the object and are captured after reflection through a camera equipped with a pixelated polarizer array of four states. However, the polarizer must be carefully adjusted for the object without knowing its optical properties, which increases the complexity of the hardware.

2.4. HDR Techniques Using Color Invariants

This method is based on the dichromatic reflection model proposed by Shafer [31], where the color of a point is determined by the color of the object itself and the color of the light source. Diffuse reflection indicates the color information related to the object itself, whereas the specular reflection represents the color information for the light source [32]. Gevers and Smeulders [33] first proposed several color invariants based on the dichromatic reflection model. Then, Benveniste and Ünsalan [3436] used color invariants to solve the problem in scanning the shiny surfaces. In 2011, they [37] introduced a novel color invariant to detect the red stripe, green stripe, and highlights to segment robustly the stripe from the captured image. By using this color invariant in segmentation, the effects of highlights originating from the ambient light were eliminated in scanning shiny surfaces. Recently, with the development of embedded technology, they [38] designed and implemented a novel structured light range scanner system based on the DSP evaluation module, which extracts range data to reduce the computational burden of the PC. Benefiting from color pattern projection and color invariants, their methods based on binary, ternary, and quaternary color-coded patterns can eliminate the effect of highlights and ambient light in the scanning process. An example of such patterns and color invariants is shown in Figure 6. However, this method based on color invariants is inevitably affected by the object itself if the object contains complex textures and multiple colors.

2.5. Photometric Stereo Technique

The photometric stereo technique can recover both the 3D shape and albedo of surfaces from multiple images at a fixed viewpoint under illuminations with different directions. In the earlier work of photometric stereo, the reflectance of surface was required to be Lambertian [39]. However, for non-Lambertian surface, especially for metallic surface, it is very difficult to robustly estimate the reflectance, while the specular region is present. Recently, the conventional photometric stereo techniques have been extended to specular surfaces with varying bidirectional reflectance distribution function (BRDF). Aiming at detecting specular components or separating reflecting components, many methods [4044] are presented in later years.

The BRDF detailed in [45] is a function that defines how light is reflected at an opaque surface. For practical surfaces, the general BRDFs are complicated nonlinear equations involving seven or more parameters [46, 47]. Moreover, the reflectance and geometrical parameters vary spatially across the surfaces. Conventional techniques normally use dense images or 3D model of measured object; however, these will limit their application. In 2010, Li et al. [48] presented an efficient photometric stereo method to estimate robustly BRDF parameters without the 3D geometric model of object. Given the small number of input images, they estimate the BRDF parameters from the best window of one image by using the recovered surface normal. In [49], Goldman et al. proposed a photometric stereo method designed for surfaces with spatially varying BRDFs, including surfaces with both varying diffuse and specular properties. They modeled the surface reflectance as a linear combination of fundamental materials and recovered the shape, material BRDFs, and weight maps using the well-known Ward model [47].

In order to tackle rough specular surfaces, on which the microfacets reflect light in different directions, producing wide specular lobes under different lightings, Hin-Shun and Jiaya [50] presented a new photometric stereo method aiming to efficiently estimate BRDF parameters and reconstruct glossy surfaces. In their approach, they did not apply unreliable highlight separation and specularity estimation. Instead, an important visual cue, that is, the cast shadow silhouette of the object, was employed to optimally recover global BRDF parameters and thereby was capable of reconstructing convex object surfaces with different level of roughness. However, their method requires high computing. For example, in their experiments, it took 50 minutes in computation to reconstruct a surface of the mouse, whereas the RMS error is 0.0326 mm. Even, in [51], it used hours to produce a surface.

In 2015, Meng et al. [52] developed a gonio-plenoptic imaging system that realizes a single-shot shape measurement for specular surfaces. The system comprised of a collimated illumination source and a plenoptic camera captures the BRDF variation of the object surface in a single image in addition to the light field information from the scene. In their experiment, they oriented the illumination in 30 different directions. However, due to the fact that their system only captures about ±7° cone of light, it constrains the range of surfaces that can be reconstructed. Furthermore, their system only works for relief surfaces and will have problems to handle the surfaces with steep slopes.

2.6. Miscellaneous Techniques

To handle specular reflections that appear in isolated regions of an image, other techniques, such as multiple cameras, color light sources, color filters [53], and multiple light projection directions [54], have been used to construct the HDR image without saturation from the different captured images. Kowarschik et al. [54] used up to 15 light projection directions to compensate the influence of specular reflections or shadowed areas. Moreover, the object can be rotated, yielding other views of the object. The camera captured different patches of the object that converted into local coordinates and then synthesized to the global coordinates. However, the numerous parameters of the system make it time-consuming to perform the calibration. In addition, additional complexity in hardware, system setup, and processing of image mask and registration is required to synthesize the different surface patches.

Alternatively, Jeong et al. [55] used the spatial light modulator to control the camera exposure time at individual pixels. The final exposure setting at saturated pixels is determined from multiple iterations. In this manner of adjusting the transmittance level of the liquid crystal spatial light modulator located at the conjugate position of the CCD plane, the projected fringes are well imaged and synthesized to the HDR images. However, additional optical and control hardware is required.

Ri et al. [56] proposed an intensity extension method using a digital micromirror device camera, where the optimum exposure time for each CCD pixel can be adjusted; thus the HDR images can be acquired.

Compared with the image intensity-based sinusoidal pattern that is vulnerable to the reflective nature of shiny surfaces, the edges of binary stripe pattern are more localized and better preserved despite the specular nature of the object surfaces; therefore, it makes the decoding procedure more robust. Song et al. [57] presented a novel structured light coding strategy for shiny objects measurement, where stripe edges instead of raw image intensities are encoded in the projected patterns. An improved zero-crossing feature detector using both positive and negative stripe patterns was proposed for stripe edge localization, leading to subpixel accuracy. However, if the saturation regions in the captured image are too large to be included in a single stripe, no constructible information can be found in these regions.

More recently, Jiang et al. proposed a technique [58] that does not require the change of exposures. The fundamental idea is that, besides capturing regular fringe patterns, inverted fringe patterns are captured to complement regular fringe patterns for phase retrieval. If not all of the regular fringe patterns are saturated, inverted fringe patterns are used instead of original saturated patterns for phase retrieval. If all of the regular fringe patterns are saturated, both the original and inverted fringe patterns are all used for phase computation to reduce phase error. Though not as robust as the previously proposed HDR technique using multiple exposures, the proposed technique can substantially increase the measurement quality for high-contract surfaces in real time.

3. Discussion

Table 1 summarizes the pros and cons of the representative HDR techniques. The adaptive digital fringe projection technique developed by Lin et al. [24] provides a good and accurate 3D shape measurement of shiny surfaces and seems to have high potential for industrial application. The HDR digital fringe projection technique proposed by Jiang et al. [13, 14] is also promising, which can obtain a fast and effective measurement of objects with high-reflective surfaces under ambient light. Feng et al. [11] presented a general solution involving three concrete techniques for realizing HDR 3D shape measurement of three types of objects, which is effective. The HDR techniques using polarizing filters are accurate and fast but relatively complex in hardware compared to other techniques. Good results were also obtained by Benveniste and Ünsalan [38], but this technique is easily affected by the colors and complex textures of the object. The gonio-plenoptic imaging system developed by Meng et al. [52] can measure the specular surfaces in high accuracy. However, due to the limitation of their system’s structure, only a small range of surface can be measured at a time; thus, it will take a long time to measure the whole surface. Therefore, it is not suitable for online applications. In addition, their system will have measuring problems to handle the surfaces with steep slopes.

Table 2 summarizes the complexity of the representative HDR techniques. The technique proposed by Jiang et al. [13] can provide an accurate 3D shape measurement of the standard porcelain balls with high-reflective surfaces. This technique takes the highest modulation intensity rather than the brightest unsaturated intensity as the criterion for selecting the pixels from the raw fringe images to synthesize the HDR images. As a result, it requires a large number of square root operations to calculate the modulation intensity at pixel level. Therefore, their algorithm is complicated. The adaptive digital fringe projection technique developed by Lin et al. [24] can measure the workpiece made of aluminum alloy and have a good performance with high accuracy. As their algorithm is simple, it seems to have high potential for industrial application. Additionally, Feng et al. [11] reported a general solution for realizing HDR 3D shape measurement based on digital fringe projection. For the proposed method by polarization, in their experiment of measuring a planar black/white calibration board, the used polarizers reduced intensity of the whole fringe images resulting in low SNR; thus, it reduces the measurement accuracy. Benefiting from using quaternary color-coded patterns, the algorithm proposed by Benveniste and Ünsalan [38] is the simplest; however, the resolution is limited, so the measurement accuracy is not high. In contrast, the gonio-plenoptic imaging system developed by Meng et al. [52] can achieve very high accuracyl however, their system can only measure a small range of surface at a time; thus, it will take a long time to measure the whole surface. Therefore, it is not suitable for online applications.

4. Conclusion

This paper presents a survey of recent HDR techniques for the digitization of shiny surfaces. All techniques were briefly recalled and compared with each other, with some highlights on their respective advantages and drawbacks.

The HDR techniques using multiple exposures are the first type of digital fringe projection technique used to obtain 3D data from shiny surfaces. The advantage of these techniques is that there is no need to make too many changes based on the original technique. They only require some postprocessing for the captured raw fringe images, and then a high precision 3D point cloud can be obtained. However, they may be time-consuming due to a large number of images needed to be captured in some cases. Fortunately, Jiang et al. [13, 14] improved this technique by developing a very fast HDR fringe pattern projector, thus reducing the time consumed. The HDR techniques using polarizing filters were also interesting, but the only drawback is the complexity of the hardware due to the use of the polarizer. The HDR techniques using color invariants provide good results but suffer from the colors and complex textures of the object being measured. Adjusting the projected fringe pattern intensities is another option to handle shiny surfaces. The adaptive digital fringe projection technique recently proposed by Lin et al. [24] has achieved high accuracy, and its novel mathematical model for shiny surface measurement shows some potential for industrial applications. A scanning device, which will merge these techniques, would also be interesting to investigate in the near future.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China and the Guangdong Provincial Natural Science Foundation under Grants no. 51675106, no. U1601202, and no. 2015A030312008, the Guangdong Provincial R&D Key Projects (no. 2015B010104008, no. 2016A030308016, and no. 2015A010104009), and the Young Talent Innovation Projects from Guangdong Educational Department (no. 2015KQNCX147).