An Array-Type System Applied to Complex Surfaces in Nuclear Pollution Detection

Radioactive pollution detection plays a key role in nuclear technology application. In this paper, an array-type of nuclear pollution detection system is designed for the detection scenario of complex surfaces. Firstly, to get the three-dimensional point cloud of the surface, a complex surface was modeled based on the geometric ranging model of a two-dimensional laser profilometer and the motion model of a two-degree-of-freedom displacement platform. Secondly, an ‘S’ type scanning scheme of profilometer was developed to overcome the problem of limited scanning area of the profilometer. Thirdly, Euclidean distance weighted median filtering was used to solve the impulsive noise that may occur during the point cloud acquisition process. Finally, the 3D point cloud information of the complex surface was used for controlling the movement of the 6 × 6 array channel pollution detector to complete the α and β particle measurement tasks. A mechanical platform was constructed for experiments, the results are as follows. The working range of this system is from −5 cm to 5 cm in elevation difference of surfaces, and the accuracy is 12 μm in surface height measuring. It takes 26.13 s to perform a detection task including surface scanning and the detector moving, and scanning accuracy is 0.35 × 0.35 mm2. The maximum control error of the surface contamination detector is 0.4 mm. Specifically, the detection area of the system reaches 240 × 240 mm2. The results show that the system acquires three-dimensional terrain information, and realizes control over the movement of the pollution detector accurately and then completes the detection of α and β particles effectively.


Introduction
In recent years, the extensive application of nuclear technology has brought a lot of convenience to human survival, but some potential harm to the natural environment that has also attracted the attention of the international community. Radioactive radiation has been listed as one of the main public hazards by the United Nations Conference on Human Environmental Protection (UNCHEP) [1][2][3][4]. Radioactive pollution monitoring technology is an important part of nuclear technology. In the monitoring process, α and β particles are especially targeted due to such characteristics as weak penetration and strong ionization. These are devastating abilities that may destroy body tissues to seriously threaten human health [5]. This is why nuclear safety supervision should always be strictly implemented so as to minimize this sort of adverse effect. It is of great significance to detect α and β particles quickly and accurately on the surface of an object and the surrounding environment.
The traditional detection for nuclear pollution methods mainly include manual hand-held detection, mobile robot for detection and flying robot for detection. The manual detection [6][7][8][9] method has the characteristics of strong environmental adaptability and high precision, but it has high requirements for the staff, and it is difficult to keep the distance between the detector and each detection point consistent. Mobile robots [10][11][12] equipped with detector for detection of nuclear pollution require lower requirements for the staff. Comparing with manual detection, the efficiency is greatly improved. But the complex surface environment will limit the movement of mobile robots, and it is difficult to complete the task of detecting. Flying robots [13][14][15][16][17][18][19] equipped with detector are mainly used for nuclear pollution detection in the air, with a wide detection range, but there are problems such as low detection accuracy for detection of complex surfaces.
Aiming at the problems of the above detection methods in nuclear pollution detection applied to complex surfaces, an array-type of nuclear pollution detection system is designed. Firstly, the system adopts a two-dimensional laser profilometer to get three-dimensional point cloud of the surface, then controls the movement of the 6 × 6 array channel pollution detector to complete the α and β particle measurement tasks based on three-dimensional point cloud. Finally, a mechanical platform was constructed for experiments, the results show that the system completes the detection of α and β particles effectively.
The rest of the paper is organized as follows: In Section 2, there is an overview presentation on related work. In Section 3, a nuclear contamination detection scheme is introduced for the complex surfaces, including measurement methods for complex topographic surfaces, surface contamination detector with motion control strategies, and array detection units for surface contamination. In Section 4, the implementation of the system is elaborated, and three experiments are performed to verify the effectiveness and feasibility of the system. Experimental results show that the system provides a feasible method for nuclear contamination detection of complex surfaces to obtain three-dimensional information. In Section 5, the performance of the proposed method is evaluated and discussed. In Section 6, some work summary and conclusions for the proposed method are presented.

Related Work
In this section, we briefly summarize different types of nuclear pollution detection methods, including manual hand-held detection, mobile robot for detection and flying robot for detection.

Manual Detection
A PRD (personal radiation detector) can detect the dose equivalent rate at the position of the wearer [6]. An alpha/beta surface contamination monitor are mainly used to measure the surface of goods on site whether there are radioactive substances and intensity. Low background alpha and beta measuring instruments are used for accurately measuring activity in a sample in laboratory tests usually. This handheld method is flexible in scheduling, but the action is slow and range is small. To do this, personnel need to participate in a close position. The biggest trouble is the difficulty of the control over the distance between the detector and the target surface [7][8][9]. Firstly, the geometric distance measurement model and the motion model of the profilometer were placed to sample the surface of the point cloud by utilizing the 2D laser profilometer. Through coordinate transformation, the point cloud data was then registered to obtain a complete 3D Point cloud. Subsequently, the point cloud data was smoothened based on weighted median filtering of Euclidean distance. Then building 6 × 6 channel detectors on the measured surface for projection, and rasterize the point cloud into 6 × 6 small equal area blocks to map each detection channel separately, followed by calculating the maximum value of the point cloud corresponding to each measurement channel. The horizontal initial position information of the detector was integrated to calculate the input of the motor. An encoder was set up to record the motor movement count to achieve closedloop control of the pollution detector movement. Finally, after the surface contamination detector moved to the desired position, a special unit started to detect the activity of α and β particles, and the data were transmitted to the host computer for data visualization and storage.

Measurement of Complex Surface Topography
Limited by factors such as light, material, equipment installation control and cost, the commonly used 3D information measurement method for large area is not suitable for this application. Thanks to its high precision and less sensitiveness to light, a two-dimensional laser profilometer was chosen for surface topography measurement. Firstly, the geometric distance measurement model and the motion model of the profilometer were placed to sample the surface of the point cloud by utilizing the 2D laser profilometer. Through coordinate transformation, the point cloud data was then registered to obtain a complete 3D Point cloud. Subsequently, the point cloud data was smoothened based on weighted median filtering of Euclidean distance. Then building 6 × 6 channel detectors on the measured surface for projection, and rasterize the point cloud into 6 × 6 small equal area blocks to map each detection channel separately, followed by calculating the maximum value of the point cloud corresponding to each measurement channel. The horizontal initial position information of the detector was integrated to calculate the input of the motor. An encoder was set up to record the motor movement count to achieve closed-loop control of the pollution detector movement. Finally, after the surface contamination detector moved to the desired position, a special unit started to detect the activity of α and β particles, and the data were transmitted to the host computer for data visualization and storage.

Measurement of Complex Surface Topography
Limited by factors such as light, material, equipment installation control and cost, the commonly used 3D information measurement method for large area is not suitable for this application. Thanks to its high precision and less sensitiveness to light, a two-dimensional laser profilometer was chosen for surface topography measurement.
There two-dimensional laser profilometer was only to obtain the X-axis and Z-axis elevation information of the terrain, not prepared to determine the Y-axis coordinate information of the terrain. Therefore, it is necessary to obtain the three-dimensional data of the terrain by controlling the movement of the equipment and integrating the motion posture information. The complex surface topography measurement goes through surface topography 3D modeling, point cloud registration and point cloud filtering. The block diagram of the topographic measurement process is shown in Figure 2. There two-dimensional laser profilometer was only to obtain the X-axis and Z-axis elevation information of the terrain, not prepared to determine the Y-axis coordinate information of the terrain. Therefore, it is necessary to obtain the three-dimensional data of the terrain by controlling the movement of the equipment and integrating the motion posture information. The complex surface topography measurement goes through surface topography 3D modeling, point cloud registration and point cloud filtering. The block diagram of the topographic measurement process is shown in Figure 2. Cell Calibration was used to establish the positional relationship between the coordinate system of the profilometer and the measurement cell. As per the calibration parameters of the measurement window, the motion control strategy of the profilometer and the registration of point cloud data were designed. The impulse noise was then eliminated by the Euclidean distance weighted median filter, and finally the complete three-dimensional information of the surface was obtained.

Three-Dimensional Modeling of Complex Terrain Surface
A two-dimensional laser profilometer was used to model the surface terrain with the help of a two-degree-of-freedom displacement platform motion. This makes sure that the three-dimensional information of the surface be recovered quickly and accurately. First, the geometric distance measurement model of the two-dimensional laser profilometer was analyzed. In the limited measurement range, the laser profilometer was unable to cover the entire measured surface at once. Therefore, an 'S' type scanning scheme was developed to make up for the shortcomings of the equipment range, i.e., detecting the target area in blocks. This 'S' type motion parameter was connected with two factors: the geometric model of the single measurement of the profilometer and the position coordinates of the window calibration. The single measurement model of the laser profilometer is shown in Figure 3. Cell Calibration was used to establish the positional relationship between the coordinate system of the profilometer and the measurement cell. As per the calibration parameters of the measurement window, the motion control strategy of the profilometer and the registration of point cloud data were designed. The impulse noise was then eliminated by the Euclidean distance weighted median filter, and finally the complete three-dimensional information of the surface was obtained.

Three-Dimensional Modeling of Complex Terrain Surface
A two-dimensional laser profilometer was used to model the surface terrain with the help of a two-degree-of-freedom displacement platform motion. This makes sure that the three-dimensional information of the surface be recovered quickly and accurately. First, the geometric distance measurement model of the two-dimensional laser profilometer was analyzed. In the limited measurement range, the laser profilometer was unable to cover the entire measured surface at once. Therefore, an 'S' type scanning scheme was developed to make up for the shortcomings of the equipment range, i.e., detecting the target area in blocks. This 'S' type motion parameter was connected with two factors: the geometric model of the single measurement of the profilometer and the position coordinates of the window calibration. The single measurement model of the laser profilometer is shown in Figure 3. The effective measurement area of the profilometer is distributed in the shape of a trapezoid. In the measurement area of the profilometer, the bottom edge of the trapezoid is the X-axis. The X-axis intersects the vertical line Y at O as the origin to establish a two-dimensional plane coordinate system, where Point A is the end point of the upper base of the trapezoid, Point B is the point on the waist of the trapezoid, and Point C is the end point of the lower base of the trapezoid. In the extreme case, the highest point of the target surface does not exceed Point B, to avoid the height of the surface from overflowing the effective measuring range of the profilometer. The single scanning width of the profilometer cannot exceed d. The position of Points A, O and C may be fixed by the way of window calibration, and the effective scanning width of the profilometer is calculated as: To overcome the range defect of a single device, a profilometer motion path planning was designed in line with the above laser profilometer measurement model and window calibration parameters. Thus, the measurement of all three-dimensional information of the regional terrain was completed. The device movement path is shown in Figure 4. The profilometer works along an 'S' type scanning mode. The motion process from the start point to the end point includes four motion modes of acceleration, uniform speed, deceleration, and commutation. During the sampling period, the motion is maintained at a uniform speed by adopted the acceleration/deceleration control algorithm based on trapezoid-curve, so that a regular threedimensional point cloud image can be obtained to reduce the difficulty and runtime of subsequent point cloud registration. A point cloud image is a digital point set measured by a scanning device on a solid model. Each point in the point set contains its physical information. The profilometer only The effective measurement area of the profilometer is distributed in the shape of a trapezoid. In the measurement area of the profilometer, the bottom edge of the trapezoid is the X-axis. The X-axis intersects the vertical line Y at O as the origin to establish a two-dimensional plane coordinate system, where Point A is the end point of the upper base of the trapezoid, Point B is the point on the waist of the trapezoid, and Point C is the end point of the lower base of the trapezoid. In the extreme case, the highest point of the target surface does not exceed Point B, to avoid the height of the surface from overflowing the effective measuring range of the profilometer. The single scanning width of the profilometer cannot exceed d. The position of Points A, O and C may be fixed by the way of window calibration, and the effective scanning width d of the profilometer is calculated as: To overcome the range defect of a single device, a profilometer motion path planning was designed in line with the above laser profilometer measurement model and window calibration parameters. Thus, the measurement of all three-dimensional information of the regional terrain was completed. The device movement path is shown in Figure 4. The effective measurement area of the profilometer is distributed in the shape of a trapezoid. In the measurement area of the profilometer, the bottom edge of the trapezoid is the X-axis. The X-axis intersects the vertical line Y at O as the origin to establish a two-dimensional plane coordinate system, where Point A is the end point of the upper base of the trapezoid, Point B is the point on the waist of the trapezoid, and Point C is the end point of the lower base of the trapezoid. In the extreme case, the highest point of the target surface does not exceed Point B, to avoid the height of the surface from overflowing the effective measuring range of the profilometer. The single scanning width of the profilometer cannot exceed d. The position of Points A, O and C may be fixed by the way of window calibration, and the effective scanning width of the profilometer is calculated as: To overcome the range defect of a single device, a profilometer motion path planning was designed in line with the above laser profilometer measurement model and window calibration parameters. Thus, the measurement of all three-dimensional information of the regional terrain was completed. The device movement path is shown in Figure 4. The profilometer works along an 'S' type scanning mode. The motion process from the start point to the end point includes four motion modes of acceleration, uniform speed, deceleration, and commutation. During the sampling period, the motion is maintained at a uniform speed by adopted the acceleration/deceleration control algorithm based on trapezoid-curve, so that a regular threedimensional point cloud image can be obtained to reduce the difficulty and runtime of subsequent point cloud registration. A point cloud image is a digital point set measured by a scanning device on a solid model. Each point in the point set contains its physical information. The profilometer only The profilometer works along an 'S' type scanning mode. The motion process from the start point to the end point includes four motion modes of acceleration, uniform speed, deceleration, and commutation. During the sampling period, the motion is maintained at a uniform speed by adopted the acceleration/deceleration control algorithm based on trapezoid-curve, so that a regular three-dimensional point cloud image can be obtained to reduce the difficulty and runtime of subsequent Electronics 2020, 9, 1870 6 of 21 point cloud registration. A point cloud image is a digital point set measured by a scanning device on a solid model. Each point in the point set contains its physical information. The profilometer only obtains one-dimensional position and depth information. The 'S' type scanning mode, combined with laser profilometer ranging model, quickly and accurately may obtain a three-dimensional point cloud. A world coordinate system was established with the upper left corner of the measurement window as the coordinate origin. In the trapezoidal area measured by the profilometer, the center line of the bottom edge was taken as the carrier coordinate system, as shown in Figure 5.
Electronics 2020, 9, x FOR PEER REVIEW 6 of 21 obtains one-dimensional position and depth information. The 'S' type scanning mode, combined with laser profilometer ranging model, quickly and accurately may obtain a three-dimensional point cloud. A world coordinate system was established with the upper left corner of the measurement window as the coordinate origin. In the trapezoidal area measured by the profilometer, the center line of the bottom edge was taken as the carrier coordinate system, as shown in Figure 5. In the measurement window, establish a world coordinate system. − , carrier coordinate system − , make the axis in the carrier coordinate system parallel to the axis in the world coordinate system, and control the laser profilometer on the shaft moves at a constant speed. Nevertheless, there is an acceleration and deceleration process in the movement of the laser profilometer. Therefore, before the axis of the carrier coordinate system coincided with the axis of the world coordinate system, the profilometer must move at a constant speed to ensure that the coordinate system of the terrain world on the axis was evenly sampled, and the 3D model recovered later would not be distorted. Set the measured plane length to , the profilometer scan frequency to f , and the speed of the profilometer to move from 1 C point to 2 C point to v. Then in a single scan, the line resolution on the axis R is expressed as: The profile point resolution of the profilometer is determined by the performance of the sensor, then the networked point cloud set : In the figure, Point represents the three-dimensional information of any point in the measurement window, recorded as ( , , ), where the Point in the world coordinate system has the actual three-dimensional information as: where is determined by Formula (1) as the width of a single scan, and ′ is the height information measured by the laser profilometer. The measurement window width is . To avoid the missed inspection, that is, each scan area needs to overlap, then the number of scans needs to meet: In the measurement window, establish a world coordinate system.O w − X w Y w Z w , carrier coordinate system O c − X c Z c , make the X w axis in the carrier coordinate system parallel to the X w axis in the world coordinate system, and control the laser profilometer on Y w the shaft moves at a constant speed. Nevertheless, there is an acceleration and deceleration process in the movement of the laser profilometer. Therefore, before the X c axis of the carrier coordinate system coincided with the X w axis of the world coordinate system, the profilometer must move at a constant speed to ensure that the coordinate system of the terrain world Y w on the axis was evenly sampled, and the 3D model recovered later would not be distorted. Set the measured plane length to s, the profilometer scan frequency to f , and the speed of the profilometer to move from C 1 point to C 2 point to v. Then in a single scan, the line resolution on the Y w axis R is expressed as: The profile point resolution of the profilometer m is determined by the performance of the sensor, then the networked point cloud set A X w Y w Z w : In the figure, Point A represents the three-dimensional information of any point in the measurement window, recorded as (X w , Y w , Z w ), where the Point A in the world coordinate system has the actual three-dimensional information as: where d is determined by Formula (1) as the width of a single scan, and Z w is the height information measured by the laser profilometer. The measurement window width is l. To avoid the missed inspection, that is, each scan area needs to overlap, then the number of scans N T needs to meet: As shown in Figure 5, the profilometer moves from the initial position point C 1 to point C n , resulting in N T networked 3D point cloud images. After the complete scanning, the profilometer is returned to the origin.

Registration of 3D Point Cloud Mosaic
Multiple 3D point clouds were acquired under the 3D modeling method of complex terrain surface. Multiple point clouds were matched into a unified coordinate system through rotation dislocation and translation dislocation. The scanning viewpoints in this system were all in the same coordinate system, only related to the X-axis and Y-axis movement, i.e., being translated on the X-axis and reversed on the Y. The schematic diagram of point cloud registration is shown in Figure 6.
Electronics 2020, 9, x FOR PEER REVIEW 7 of 21 As shown in Figure 5, the profilometer moves from the initial position point to point , resulting in networked 3D point cloud images. After the complete scanning, the profilometer is returned to the origin.

Registration of 3D Point Cloud Mosaic
Multiple 3D point clouds were acquired under the 3D modeling method of complex terrain surface. Multiple point clouds were matched into a unified coordinate system through rotation dislocation and translation dislocation. The scanning viewpoints in this system were all in the same coordinate system, only related to the X-axis and Y-axis movement, i.e., being translated on the X-axis and reversed on the Y. The schematic diagram of point cloud registration is shown in Figure 6. Looking at the trajectory of the profilometer in Figure 5, there is a reverse relationship between the two point cloud images of and . The storage is in a sequential order, and they are opposite in space position. The point clouds need to rotate 180° for misalignment registration. Considering avoiding loss of sampled data points, there were resampling units in the two point cloud images, and the stitched boundary point cloud data was obtained by means of averaging.
where ( , ) is the point cloud data after registration, ( 1, 1) is the data obtained by the first scanning, and ( 2, 2) is the data obtained by second. After averaging, the original overlapping redundant data was substituted to get a complete 3D point cloud image.

Three-Dimensional Point Cloud Filtering
Due to errors caused by the measured terrain surface and the measurement system itself, the measured 3D point cloud may contain noises. These noises would cause the quality decreasing of the point cloud data, posing a great impact on the movement control of the probe behind. When it comes to denoising the point cloud, the proposed algorithms for ordered or partially ordered denoising means include Wiener filtering, Kalman filtering, smoothing filtering, median filtering, etc. [20]. For the impulse noise that occurs during the scanning process, the commonly used standard median filter algorithm can effectively filter out regardless of the filter output estimate that may differ from the real information in the center of the filter window. As the filter window moves, this deviation will also participate in the calculation of the center point of other filtering windows, leading to the accumulation of errors. In this work, a center-weighted median filtering algorithm was used to filter out the impulse noise. The similarity between the point clouds in the filter window was taken in to assign the weight of each point to participate in the calculation. In this way, the deviation of the filtered output value from the real value was reduced, and the true depth information of the sampling Looking at the trajectory of the profilometer in Figure 5, there is a reverse relationship between the two point cloud images of N 1 and N 2 . The storage is in a sequential order, and they are opposite in space position. The point clouds need to rotate 180 • for misalignment registration. Considering avoiding loss of sampled data points, there were k resampling units in the two point cloud images, and the stitched boundary point cloud data was obtained by means of averaging.
where A(m, n) is the point cloud data after registration, A(m1, n1) is the data obtained by the first scanning, and A(m2, n2) is the data obtained by second. After averaging, the original overlapping redundant data was substituted to get a complete 3D point cloud image.

Three-Dimensional Point Cloud Filtering
Due to errors caused by the measured terrain surface and the measurement system itself, the measured 3D point cloud may contain noises. These noises would cause the quality decreasing of the point cloud data, posing a great impact on the movement control of the probe behind. When it comes to denoising the point cloud, the proposed algorithms for ordered or partially ordered denoising means include Wiener filtering, Kalman filtering, smoothing filtering, median filtering, etc. [20]. For the impulse noise that occurs during the scanning process, the commonly used standard median filter algorithm can effectively filter out regardless of the filter output estimate that may differ from the real Electronics 2020, 9, 1870 8 of 21 information in the center of the filter window. As the filter window moves, this deviation will also participate in the calculation of the center point of other filtering windows, leading to the accumulation of errors. In this work, a center-weighted median filtering algorithm was used to filter out the impulse noise. The similarity between the point clouds in the filter window was taken in to assign the weight of each point to participate in the calculation. In this way, the deviation of the filtered output value from the real value was reduced, and the true depth information of the sampling point was better retained.
Using the 3 × 3 filter window, the first step is to use the similarity function to calculate the similarity between each point cloud in the window and the point cloud in the center of the window, and sort their sizes. Each pixel in the filter window was adaptively given a corresponding weight. Finally, perform weighted median filtering on the point cloud of the filter window center to obtain the output response of the window center. The three-dimensional networked point cloud may be projected in the two-dimensional space. In this background, the magnitude of similarity is reflected as the Euclidean distance between the center point cloud of the filter window and the other point clouds in the window. The closer the distance from the center point cloud of the filter window, the greater the similarity and the same to weight. The principle diagram of Euclidean distance weighted filtering is shown in Figure 7. and sort their sizes. Each pixel in the filter window was adaptively given a corresponding weight. Finally, perform weighted median filtering on the point cloud of the filter window center to obtain the output response of the window center. The three-dimensional networked point cloud may be projected in the two-dimensional space. In this background, the magnitude of similarity is reflected as the Euclidean distance between the center point cloud of the filter window and the other point clouds in the window. The closer the distance from the center point cloud of the filter window, the greater the similarity and the same to weight. The principle diagram of Euclidean distance weighted filtering is shown in Figure 7.
In Equation (8), stands for that the Point Cloud has the similarity between a point and We use the similarity function to calculate the similarity of each point cloud in the window. The similarity function uses the monotonically decreasing function ϕ(x) = 1 1+x in [0, +∞], where x is the distance between the point cloud and the center of the window. The similarity between the point A (i−1,j−1) and A (i,j) is defined as S (i−1,j−1) : Then 1+b . Similarly, the similarity between each point cloud in the filter window and the center point of the filter window. The distance a is related to the line resolution of the profilometer m, a = c m , b = S R , R are the scanning resolution of the profilometer in the Y w axis direction, R see Equation (2). Let P A be the sum of the similarities of the point clouds in the filter window except for the center point cloud, and then assign the weight of the point cloud in the window according to the similarity W: Electronics 2020, 9, 1870 9 of 21 In Equation (8), P A stands for that the Point Cloud A has the similarity between a point and Point Cloud A in the spatial neighborhood. u, v is the point coordinate position of the Point Cloud A area, excluding the Point Cloud A. The formula W (u,v) represents the filter weight of the Point Cloud A points in the spatial neighborhood, and finally, the value after A. The X axis and the Y axis of this system were different as far as the resolution is concerned, and the resolutions were higher. Additionally, the elevation information between adjacent point clouds had a high similarity. Each point in the traditional median filter window shared the same weight value. As a result, a large deviation occurred between the output response of the filter window and the true value. A deviation of this kind was easy to cause error transmission with the movement of the filter window. Comparatively, the median weighted filtering based on the Euclidean distance of the point cloud assigned the weight of the point cloud by the distance between the point clouds. The filter output obtained by this method may effectively reduce this deviation.

Motion Control Strategy of Surface Contamination Detector
The motion control of the probe requires terrain data as a control input. There are 36 channels in the array-type pollution detection unit. The 36 probe channels are respectively mapped in the measurement area to calculate the maximum depth value in each small area.
Ideally, the measurement area may be evenly divided into 36 parts, which, respectively, are mapped 36 measurement channels. Due to the mechanical assembly and mechanical processing, the probes were not arranged equidistantly or without gaps, and there was a large positional deviation. As shown in Figure 8. Each point in the traditional median filter window shared the same weight value. As a result, a large deviation occurred between the output response of the filter window and the true value. deviation of this kind was easy to cause error transmission with the movement of the filter window. Comparatively, the median weighted filtering based on the Euclidean distance of the point cloud assigned the weight of the point cloud by the distance between the point clouds. The filter output obtained by this method may effectively reduce this deviation.

Motion Control Strategy of Surface Contamination Detector
The motion control of the probe requires terrain data as a control input. There are 36 channels in the array-type pollution detection unit. The 36 probe channels are respectively mapped in the measurement area to calculate the maximum depth value in each small area.
Ideally, the measurement area may be evenly divided into 36 parts, which, respectively, are mapped 36 measurement channels. Due to the mechanical assembly and mechanical processing, the probes were not arranged equidistantly or without gaps, and there was a large positional deviation. As shown in Figure 8. The surface contamination detector moves under the control of high-precision stepping motor and screw drive. The communication works on the RS485 bus that a full-duplex communication, configured with a number of detector channels. To ensure communication efficiency and response speed, two groups of 485 buses were combined to communicate with 18 channel motors. In addition, a code disk was added to the stepper motor to achieve closed-loop control, and a brake device was engaged on the shaft to complete the motor self-locking function. The terrain information of each measurement was transmitted to the main control system, which distributed the supply of the stepper motor mapped to the terrain with the counter on the motor code wheel as control feedback. The block diagram of the motor control system is shown in Figure 9.  The surface contamination detector moves under the control of high-precision stepping motor and screw drive. The communication works on the RS485 bus that a full-duplex communication, configured with a number of detector channels. To ensure communication efficiency and response speed, two groups of 485 buses were combined to communicate with 18 channel motors. In addition, a code disk was added to the stepper motor to achieve closed-loop control, and a brake device was engaged on the shaft to complete the motor self-locking function. The terrain information of each measurement was transmitted to the main control system, which distributed the supply of the stepper motor mapped to the terrain with the counter on the motor code wheel as control feedback. The block diagram of the motor control system is shown in Figure 9. speed, two groups of 485 buses were combined to communicate with 18 channel motors. In addition, a code disk was added to the stepper motor to achieve closed-loop control, and a brake device was engaged on the shaft to complete the motor self-locking function. The terrain information of each measurement was transmitted to the main control system, which distributed the supply of the stepper motor mapped to the terrain with the counter on the motor code wheel as control feedback. The block diagram of the motor control system is shown in Figure 9. There are two inputs for motor control, one for the measured terrain height and the other for the initial position of the motor. By calculating the two inputs as the motor control feed, the motor code disc count value is read as control feedback to control the motor movement. Due to the mechanical installation errors of the detector, lead screw and motor, the initial position of each detector may be different from each other. The mechanical position of each channel detector should be calibrated and compensated for errors. The actual differences are seen in Figure 10. There are two inputs for motor control, one for the measured terrain height and the other for the initial position of the motor. By calculating the two inputs as the motor control feed, the motor code disc count value is read as control feedback to control the motor movement. Due to the mechanical installation errors of the detector, lead screw and motor, the initial position of each detector may be different from each other. The mechanical position of each channel detector should be calibrated and compensated for errors. The actual differences are seen in Figure 10.

Surface Pollution Array Detection Unit
The surface pollution detection unit is a device used to measure the radioactivity on the target surface. A semiconductor α and β particles detector was arranged to detect the target. The surface contamination detection subsystem consists of a front-end sensor measurement channel, a probe data processing box and an industrial computer. Distributed in 6 × 6 array, a single channel is composed of a surface contamination detector, an adjustable screw and a stepper motor. The single channel is shown in Figure 12.

Surface Pollution Array Detection Unit
The surface pollution detection unit is a device used to measure the radioactivity on the target surface. A semiconductor α and β particles detector was arranged to detect the target. The surface contamination detection subsystem consists of a front-end sensor measurement channel, a probe data processing box and an industrial computer. Distributed in 6 × 6 array, a single channel is composed of a surface contamination detector, an adjustable screw and a stepper motor. The single channel is shown in Figure 12.

Surface Pollution Array Detection Unit
The surface pollution detection unit is a device used to measure the radioactivity on the target surface. A semiconductor α and β particles detector was arranged to detect the target. The surface contamination detection subsystem consists of a front-end sensor measurement channel, a probe data processing box and an industrial computer. Distributed in 6 × 6 array, a single channel is composed of a surface contamination detector, an adjustable screw and a stepper motor. The single channel is shown in Figure 12. The stepping motor was available to adjust the distance between the probe surface and the measured surface. The set distance was in a range of 1 mm to 5 mm, the counting time may be 1 s, 6 s, 10 s, 30 s and 60 s in several manual setting modes. Each channel and the digital communication controller were connected by a power line and a data line, which may simultaneously measure two kinds of particles: α, β. The digital communication controller sampled the data of the data probe and processes it into a digital signal. Via UART (Universal Asynchronous Receiver/Transmitter), this signal was then transmitted to the digital display platform that presents the data to the user.
The surface pollution array detection unit adopts α and β particles scintillator detectors. The scintillator detector is a kind of widely used radiation detector which uses the flash of radiation produced in certain substances for radiation detection. The incident radiation would get into the scintillator and loses energy in the scintillator, which could cause ionization and excitation of the scintillator atoms and generate visible light. The photon hits the photomultiplier tube through the light transmission process to produce photoelectrons. Finally, the electrons reach the anode and form a voltage pulse. Suppose is the multiplication coefficient of the photomultiplier (PMT) tube, is the number of photoelectrons generated by photons hitting the cathode, then the number of electrons collected at the output terminal is , and the corresponding pulse charge number is can be obtained by = Then the voltage pulse : where is the output capacitance, = / is the number of photoelectrons produced per unit energy of the incident particle. is the amount of charge a unit photoelectron has, and is the incident particle energy. The α and β particle detectors used in the proposed approach select a composite scintillator composed of ZnS (Ag) scintillators and plastic scintillators. The composite scintillator is uniformly coated with ZnS (Ag) on a thin plastic scintillator. The ZnS (Ag) crystal can absorb all incident α particles, and the β particles with lower energy, which can pass through the ZnS (Ag) crystal. Then, the β particles will enter the plastic scintillator and will be recorded. In addition, this composite scintillator can be made into a variety of sizes and shapes, which is suitable for the surface pollution detection unit designed herein. The count rate is obtained by the statistical number of voltage pulses to characterize the radioactive activity of the radioactive source. The detector's count rate for the standard planar source divided by the standard planar source unit area activity is the surface activity response of the pollution instrument (s ). The stepping motor was available to adjust the distance between the probe surface and the measured surface. The set distance was in a range of 1 mm to 5 mm, the counting time may be 1 s, 6 s, 10 s, 30 s and 60 s in several manual setting modes. Each channel and the digital communication controller were connected by a power line and a data line, which may simultaneously measure two kinds of particles: α, β. The digital communication controller sampled the data of the data probe and processes it into a digital signal. Via UART (Universal Asynchronous Receiver/Transmitter), this signal was then transmitted to the digital display platform that presents the data to the user.
The surface pollution array detection unit adopts α and β particles scintillator detectors. The scintillator detector is a kind of widely used radiation detector which uses the flash of radiation produced in certain substances for radiation detection. The incident radiation would get into the scintillator and loses energy in the scintillator, which could cause ionization and excitation of the scintillator atoms and generate visible light. The photon hits the photomultiplier tube through the light transmission process to produce photoelectrons. Finally, the electrons reach the anode and form a voltage pulse. Suppose M is the multiplication coefficient of the photomultiplier (PMT) tube, N is the number of photoelectrons generated by photons hitting the cathode, then the number of electrons collected at the output terminal is MN, and the corresponding pulse charge number is Q can be obtained by Then the voltage pulse V: where C is the output capacitance, n = N/E is the number of photoelectrons produced per unit energy of the incident particle. e is the amount of charge a unit photoelectron has, and E is the incident particle energy. The α and β particle detectors used in the proposed approach select a composite scintillator composed of ZnS (Ag) scintillators and plastic scintillators. The composite scintillator is uniformly coated with ZnS (Ag) on a thin plastic scintillator. The ZnS (Ag) crystal can absorb all incident α particles, and the β particles with lower energy, which can pass through the ZnS (Ag) crystal. Then, the β particles will enter the plastic scintillator and will be recorded. In addition, this composite scintillator can be made into a variety of sizes and shapes, which is suitable for the surface pollution detection unit designed herein. The count rate is obtained by the statistical number of voltage pulses to characterize the radioactive activity of the radioactive source. The detector's count rate for the standard planar source divided by the standard planar source unit area activity is the surface activity response of the pollution instrument (s −1 Bq −1 cm 2 ).
where R a is the surface activity response of the surface contamination detector for α and β. N is the counting rate of the detector to the standard source, N b is the counting rate of the detector to the standard source. A s is the activity per unit area of the standard plane source. In practical applications, the reference value of the α and β standard planar sources is the surface emissivity, and the surface source emissivity expressed as q is related to the radioactivity through the planar source efficiency ε.
where A is the radioactivity of the source, q is the surface emissivity of the planar source. Following this, we divide the radioactivity of the planar source by the surface area of the planar source to get A s . Normally for α-plane sources, ε(α) = 0.51, ε(β) = 0.62. Relatively inherent errors are usually used to describe the accuracy of system measurements.

System Implementation and Experimental Verification
This section describes the specific implementation of the proposed method and discusses the experiment of related indicators. As the control foundation of the system, the experimental platform must be reasonable and reliable. The experimental platform mainly includes the overall structure, the motion platform of laser profilometer, the control platform for detector motion and the computer processing platform. We designed three experiments to verify the effectiveness of the system. The first experiment tested the time consumption and scanning accuracy of the system, and analyzed the effectiveness of the surface profile measurement subsystem; the second experiment separately tested the control accuracy of 6 × 6 channels, plus the reliability of our method in different scenarios; the third experiment aimed at the activity response of the surface contamination detection system to a standard planar source, testing the effect of different measurement heights on the surface contamination detection error, qualitatively analyzing the feasibility of the system.

System Implementation
The overall mechanical structure of the system is like a hand-push box, equipped with four universal wheels as the center of gravity support, featuring high maneuverability. The mechanical structure of the system includes a laser profilometer motion platform, which is a two-dimensional LLT-2900-100 high-speed profilometer produced by Germany-based Micro-Epsilon. It operates with a semiconductor laser having a wavelength of 658 nm (visible/red). When the measured object is a complex surface with high transmittance, such as glass, the depth/geometric profiling will be affected. The measurement scenes involved in this paper are mainly metal surface and concrete surface, in generally the system will not be used in the case that reflected light is blocked obviously, so the profilometer is chosen to meet the requirements of the system. The motion platform consists of an X-axis motion platform, a Y-axis rectangular coordinate motion platform and a cable protection towline. The combination of two displacement platforms enforces the complete movement of the profilometer in the horizontal plane of the detection area. The control platform for 6 × 6-channel detectors motion contains a total of 36 detection channels that can be moved simultaneously or separately. Each of these channel packages incorporates MOONS's STM11-3RE series integrated stepper motors and the custom couplings and radiation detectors, aiming at high-precision control. The user operation platform is composed of industrial computers and located in the most mechanical structure of the system upper layer. The mechanical structure of the system is shown in Figure 13. According to the mechanical structure of the system, the corresponding physical platform was built, as shown in Figure 14. _ is the Z-axis ready to control the up and down movement of the 6 × 6 array detection channel platform, _ is the X-axis and _ is the Y-axis, and _ is responsible for the single-channel detector movement. The overall support frame was built with anodized aluminum profiles in conformity with the European Commission's standard; _ and _ worked with the TL170 single-arm motion platform. This is a semi-closed screw series with the motor power up to 400 W and the repeat accuracy ±0.02 mm; _ used the M2DV-4D52R servo motor as the drive, and the mechanical layout of the detector was organized in a 6 × 6 array. As to the experimental platform, the overall dimensions are L900 × W800 × 1300 mm 3 , the X-axis travel is 400 mm, the Y-axis travel is up to 350 mm, the Z-axis travel is 450 mm, and the surface contamination detector travels up and down approximately 110 mm. The software implementation part operated on Windows 10 as its system environment, being equipped with Visual Studio 2015 and Qt5.8. The software design part mainly includes acquisition and processing of point cloud data, motion control of the 6 × 6 detection channel, data reading and storage of surface pollution detector as well as terrain surface data and visualization of pollution detector data. The software flow chart is shown in Figure 15.
After the program starts, the initial position of the horizontal level of the detection channel will be self-checked. If the self-checking fails, it will prompt the coordinates of the abnormal channel position, and then perform the reset operation. After the self-test passes, the profilometer is initialized, and the point cloud data is collected and processed. The maximum height of the terrain corresponding to the 6 × 6 channels is calculated, and then converted to the motor feed. Subsequently, According to the mechanical structure of the system, the corresponding physical platform was built, as shown in Figure 14. M_Z is the Z-axis ready to control the up and down movement of the 6 × 6 array detection channel platform, M_X is the X-axis and M_Y is the Y-axis, and M_S is responsible for the single-channel detector movement. The overall support frame was built with anodized aluminum profiles in conformity with the European Commission's standard; M_X and M_Y worked with the TL170 single-arm motion platform. This is a semi-closed screw series with the motor power up to 400 W and the repeat accuracy ±0.02 mm; M_Z used the M2DV-4D52R servo motor as the drive, and the mechanical layout of the detector was organized in a 6 × 6 array. As to the experimental platform, the overall dimensions are L900 × W800 × 1300 mm 3 , the X-axis travel is 400 mm, the Y-axis travel is up to 350 mm, the Z-axis travel is 450 mm, and the surface contamination detector travels up and down approximately 110 mm. According to the mechanical structure of the system, the corresponding physical platform was built, as shown in Figure 14. _ is the Z-axis ready to control the up and down movement of the 6 × 6 array detection channel platform, _ is the X-axis and _ is the Y-axis, and _ is responsible for the single-channel detector movement. The overall support frame was built with anodized aluminum profiles in conformity with the European Commission's standard; _ and _ worked with the TL170 single-arm motion platform. This is a semi-closed screw series with the motor power up to 400 W and the repeat accuracy ±0.02 mm; _ used the M2DV-4D52R servo motor as the drive, and the mechanical layout of the detector was organized in a 6 × 6 array. As to the experimental platform, the overall dimensions are L900 × W800 × 1300 mm 3 , the X-axis travel is 400 mm, the Y-axis travel is up to 350 mm, the Z-axis travel is 450 mm, and the surface contamination detector travels up and down approximately 110 mm. The software implementation part operated on Windows 10 as its system environment, being equipped with Visual Studio 2015 and Qt5.8. The software design part mainly includes acquisition and processing of point cloud data, motion control of the 6 × 6 detection channel, data reading and storage of surface pollution detector as well as terrain surface data and visualization of pollution detector data. The software flow chart is shown in Figure 15.
After the program starts, the initial position of the horizontal level of the detection channel will be self-checked. If the self-checking fails, it will prompt the coordinates of the abnormal channel position, and then perform the reset operation. After the self-test passes, the profilometer is initialized, and the point cloud data is collected and processed. The maximum height of the terrain The software implementation part operated on Windows 10 as its system environment, being equipped with Visual Studio 2015 and Qt5.8. The software design part mainly includes acquisition and processing of point cloud data, motion control of the 6 × 6 detection channel, data reading and storage Electronics 2020, 9, 1870 14 of 21 of surface pollution detector as well as terrain surface data and visualization of pollution detector data. The software flow chart is shown in Figure 15.
α and β particles are collected, the detector data is read and stored in the database, and finally the terrain surface data and the detector measurement data are visualized in 3D and the data report is printed out to complete the above process prior to resetting the system. Figure 16 is the effect diagram of the software interface, where (a) is the main interface of the system, including six buttons for such functions as detection, reset, and stop, etc. The stop button has the highest priority and it interrupts the system process after the system is found to be abnormal, protecting the system from being damaged; (b) is a 3D histogram of the location of each detection channel and its corresponding surface terrain height and pollution meter detection value; (c) shows the interface as a visual interactive interface for historical data query, querying historical data at each time period, available for the printout of data reports.  After the program starts, the initial position of the horizontal level of the detection channel will be self-checked. If the self-checking fails, it will prompt the coordinates of the abnormal channel position, and then perform the reset operation. After the self-test passes, the profilometer is initialized, and the point cloud data is collected and processed. The maximum height of the terrain corresponding to the 6 × 6 channels is calculated, and then converted to the motor feed. Subsequently, the pollution detector is controlled to move downwards. After moving to a given position, the surface α and β particles are collected, the detector data is read and stored in the database, and finally the terrain surface data and the detector measurement data are visualized in 3D and the data report is printed out to complete the above process prior to resetting the system. Figure 16 is the effect diagram of the software interface, where (a) is the main interface of the system, including six buttons for such functions as detection, reset, and stop, etc. The stop button has the highest priority and it interrupts the system process after the system is found to be abnormal, protecting the system from being damaged; (b) is a 3D histogram of the location of each detection channel and its corresponding surface terrain height and pollution meter detection value; (c) shows the interface as a visual interactive interface for historical data query, querying historical data at each time period, available for the printout of data reports.

Evaluation of Runtime Overhead Test of Array Detection Method
The authors tested the system running time overhead and surface topography scanning accuracy. According to the requirements of the system time index, the movement speed of each degree of freedom needs to be specially set. For too low a speed, the overhead may increase. For a higher one, there may be a step loss under load, and so is the control accuracy. In the worst case, the inertial motion carrier may exceed the limit switch. _ axis and _ axis are driven by belt, _ axis and _ axis by screw. See the transmission device parameters of the experimental platform are shown in Table 1. When testing the time overhead of the analysis system, we divided the movement of the profilometer into 7 segments ( 1 to 7), and the downward movement of the 6 × 6 array channel was divided into 8 to 9, as shown in Figure 17, where O is the initial position.

Evaluation of Runtime Overhead Test of Array Detection Method
The authors tested the system running time overhead and surface topography scanning accuracy. According to the requirements of the system time index, the movement speed of each degree of freedom needs to be specially set. For too low a speed, the overhead may increase. For a higher one, there may be a step loss under load, and so is the control accuracy. In the worst case, the inertial motion carrier may exceed the limit switch. M_X axis and M_Y axis are driven by belt, M_Z axis and M_S axis by screw. See the transmission device parameters of the experimental platform are shown in Table 1. When testing the time overhead of the analysis system, we divided the movement of the profilometer into 7 segments (t1 to t7), and the downward movement of the 6 × 6 array channel was divided into t8 to t9, as shown in Figure 17, where O is the initial position.

Evaluation of Runtime Overhead Test of Array Detection Method
The authors tested the system running time overhead and surface topography scanning accuracy. According to the requirements of the system time index, the movement speed of each degree of freedom needs to be specially set. For too low a speed, the overhead may increase. For a higher one, there may be a step loss under load, and so is the control accuracy. In the worst case, the inertial motion carrier may exceed the limit switch. _ axis and _ axis are driven by belt, _ axis and _ axis by screw. See the transmission device parameters of the experimental platform are shown in Table 1. When testing the time overhead of the analysis system, we divided the movement of the profilometer into 7 segments ( 1 to 7), and the downward movement of the 6 × 6 array channel was divided into 8 to 9, as shown in Figure 17, where O is the initial position.  The movement time of each process was tested separately and compared with the theoretical calculation value, see Table 2. t8 is the descent time of the Z-axis with the 6 × 6 array platform, the descent distance is 15 cm; t9 is the descent movement time of a single measurement channel. If 36 channels are moving at the same time, the maximum descent range is 10 cm according to the index requirements, so the theoretical value of the maximum movement time is 5 s. From Table 2, the actual test time consumption increased compared to the theoretical value. The theoretical total motion time is 24.13 s and the actual measured time is 26.13 s. Both the theoretical and actual measured values meet the system time requirements. The reasons why the test result was higher than the theoretical runtime are as follows: 1.
Two groups of 485 communication method was adopted, each way mounted with 18 devices, the communication latency was 200 ms.

2.
Unlike the value calculated theoretically for no-load, the motor was loaded during the actual test.
The speed was lower than the theory under the load with the same parameters. 3.
In the comparison table, it is particularly found that the difference between the test value and the theoretical value in t8 is the largest, and the storage and calculation of the three-dimensional point cloud is included in this movement time period.
The motion parameters of the profilometer may also affect the sampling rate of the topography of the measured surface. The resolution of a single contour of the profilometer line is 320, the length of the line is 112 mm, and the operating frequency is 270 Hz. The sampling platforms are t2, t4 and t6. Table 3 is the resolution we actually measured. From Table 3, the actual value of the X-axis resolution is consistent with the theoretical value, which is determined by the built-in parameters of the profilometer. The theoretical resolution of the Y-axis can be calculated by Equation (2), which is not much different from the measured value. The slight difference is due to a slight fluctuation in the operating frequency of the profilometer, and there is a small error between the measured value and the theoretical value due to the accumulation of time. To facilitate the registration of the point cloud, we unified the resolution on the Y-axis. It would be more appropriate to measure the Y-axis resolution to 689 through 100 experiments, so the surface point cloud resolution is 689 × 686 (L × W), and the scanning accuracy is 0.35 × 0.35 mm 2 .

Evaluation of Motion Control Effect of Detector
A standard horizontal platform is selected as the measurement object, the system is placed on the horizontal platform in parallel. The system scans the platform, then controls the detector to descend, and sets the distance between the detector horizontal plane and the horizontal platform at 1 mm. Some standard measurement tools are adopted to measure the distance between the detector and the ground, as shown in Figure 18.  Figure 18a is a standard thickness circular gasket verified by a vernier caliper, ranging from 0.2 mm to 1.5 mm. After the end of detector descending, the standard tool above was used to measure the detector motion control error. The error test of motion control was carried out on 36 channels, each for 50 times. The test statistical results are shown in Figure 19. The test results are displayed in the paper in the form of a box chart, describing the distribution of error data for each channel. In Figure 19, the error box chart for channels 0 to 35. The motion control errors of all detection channels were basically distributed from −0.2 mm to +0.2 mm, with the maximum error of 0.4 mm.

Evaluation of Detection Efficiency for Surface Pollution
The surface activity response is the most important parameter to measure the performance of the surface pollution meters with α and β particles. During the measurement process, accurate control of distance is critical to the accuracy of surface activity response magnitude. This experiment was designed to verify the performance of the surface contamination detection subsystem. The standard plane sources of α particle and β particle were used to give a qualitative and quantitative analysis of the influence of the distance on the surface activity response of the surface contamination detectors for α and β particles. Figure 20 is the physical map of the surface pollution detection system. In addition to a 36-channel pollution detector, the system is also equipped with a digital communication controller dedicated to storing and uploading the measurement data of the detector.  Figure 18a is a standard thickness circular gasket verified by a vernier caliper, ranging from 0.2 mm to 1.5 mm. After the end of detector descending, the standard tool above was used to measure the detector motion control error. The error test of motion control was carried out on 36 channels, each for 50 times. The test statistical results are shown in Figure 19.  Figure 18a is a standard thickness circular gasket verified by a vernier caliper, ranging from 0.2 mm to 1.5 mm. After the end of detector descending, the standard tool above was used to measure the detector motion control error. The error test of motion control was carried out on 36 channels, each for 50 times. The test statistical results are shown in Figure 19. The test results are displayed in the paper in the form of a box chart, describing the distribution of error data for each channel. In Figure 19, the error box chart for channels 0 to 35. The motion control errors of all detection channels were basically distributed from −0.2 mm to +0.2 mm, with the maximum error of 0.4 mm.

Evaluation of Detection Efficiency for Surface Pollution
The surface activity response is the most important parameter to measure the performance of the surface pollution meters with α and β particles. During the measurement process, accurate control of distance is critical to the accuracy of surface activity response magnitude. This experiment was designed to verify the performance of the surface contamination detection subsystem. The standard plane sources of α particle and β particle were used to give a qualitative and quantitative analysis of the influence of the distance on the surface activity response of the surface contamination detectors for α and β particles. Figure 20 is the physical map of the surface pollution detection system. In addition to a 36-channel pollution detector, the system is also equipped with a digital communication controller dedicated to storing and uploading the measurement data of the detector. The test results are displayed in the paper in the form of a box chart, describing the distribution of error data for each channel. In Figure 19, the error box chart for channels 0 to 35. The motion control errors of all detection channels were basically distributed from −0.2 mm to +0.2 mm, with the maximum error of 0.4 mm.

Evaluation of Detection Efficiency for Surface Pollution
The surface activity response is the most important parameter to measure the performance of the surface pollution meters with α and β particles. During the measurement process, accurate control of distance is critical to the accuracy of surface activity response magnitude. This experiment was designed to verify the performance of the surface contamination detection subsystem. The standard plane sources of α particle and β particle were used to give a qualitative and quantitative analysis of the influence of the distance on the surface activity response of the surface contamination detectors for α and β particles. Figure 20 is the physical map of the surface pollution detection system. In addition to a 36-channel pollution detector, the system is also equipped with a digital communication controller dedicated to storing and uploading the measurement data of the detector. The standard plane source of α particle used in the experiment was Am-241. The rest data included an active area 150 cm 2 , surface activity 14.58 Bq cm 2 , Source of β particle was T1-204, including an active area 150 cm 2 , surface activity 11.69 Bq cm 2 . The reading value of the detector is the count rate ( ), the area of the detection surface is 4 × 4 cm. It is believed that each point within plane sources of α particle and β particle shared the same surface activity, and 36 channels may be detected at the same height at the same time. Using the upper left corner of the detector as the starting point, the resolution number was 0 to 35, and the surface activity response of each channel was calculated separately. Then adjust the distance between the detection plane of the detector and the plane source, continue to test the surface activity response of each channel, adjust the distance to 3-7 mm. Meanwhile, the test value at 3 mm was agreed as baseline, with the detection efficiency here set to 100%. Next, the two sources (α, β) were separately detected by the surface pollution detection system. The experiments were performed at different distances, and the count values of 36 channels were averaged.
In this paper, when the detector measures α particle, the distance adjustment range is 3-5 mm. For β particle, the distance adjustment is 8-12 mm. The activity response of the radioactive source is measured at different distances, and the measurement efficiency of the instrument at different distances is obtained. According to Tables 4 and 5, it can be concluded that the smaller the distance between the detector and the source, the higher the surface activity response value and the higher the detection efficiency. The distance has a great influence on the detection result. In actual use, the accuracy of the distance needs to be strictly maintained. The system can meet the requirements of precise distance control  The standard plane source of α particle used in the experiment was Am-241. The rest data included an active area 150 cm 2 , surface activity 14.58 Bq cm 2 , Source of β particle was T1-204, including an active area 150 cm 2 , surface activity 11.69 Bq cm 2 . The reading value of the detector is the count rate (s −1 ), the area of the detection surface is 4 × 4 cm. It is believed that each point within plane sources of α particle and β particle shared the same surface activity, and 36 channels may be detected at the same height at the same time. Using the upper left corner of the detector as the starting point, the resolution number was 0 to 35, and the surface activity response of each channel was calculated separately. Then adjust the distance between the detection plane of the detector and the plane source, continue to test the surface activity response of each channel, adjust the distance to 3-7 mm. Meanwhile, the test value at 3 mm was agreed as baseline, with the detection efficiency here set to 100%. Next, the two sources (α, β) were separately detected by the surface pollution detection system. The experiments were performed at different distances, and the count values of 36 channels were averaged.
In this paper, when the detector measures α particle, the distance adjustment range is 3-5 mm. For β particle, the distance adjustment is 8-12 mm. The activity response of the radioactive source is measured at different distances, and the measurement efficiency of the instrument at different distances is obtained. According to Tables 4 and 5, it can be concluded that the smaller the distance between the detector and the source, the higher the surface activity response value and the higher the detection efficiency. The distance has a great influence on the detection result. In actual use, the accuracy of the distance needs to be strictly maintained. The system can meet the requirements of precise distance control.  Table 5. The source of β particle surface activity response test at difference distance.

Performance of System and Discussion
Experiments on the experimental platform showed that it could automatically measure complex surfaces with a height difference of ±5 cm, being able to better restore the three-dimensional information of the terrain surface. The scanning accuracy was as good as 0.35 mm, and the surface height measurement accuracy was 12 µm, along with the single measurement area up to 240 × 240 mm 2 . The runtime of completing a surface 3D information measurement and motion control was 26.13 s, and the maximum control error of the detector movement was 0.4 mm. Tests were also separately performed on different lighting, surface materials and venues. Experiments have shown that the system was able to respond to some basic environmental changes for motion control of the pollution detector. To verify the feasibility of the detection of surface nuclear pollution intensity, tests were also done on condition of standard planar sources (α, β). The effects of different detection distances on the detection efficiency and maximum relative error of the system were analyzed accordingly.
In complex surface environments, our method has obvious advantages. Compared with manual detection [6], our method has the characteristics of large detection range and high efficiency. For a mobile robot [10], our method is less affected by the ground environment. For a flying robot [17], our method has higher efficiency and precision.
In the motion control of the profilometer, the acceleration/deceleration control algorithm based on trapezoid-curve is adopted to ensure that the profilometer is in a uniform motion state during the scanning process. Further, the motion scanning scheme of an 'S' shaped curve is applied to ensure the integrity of the scanning area, which can reduce the operating time. Besides, two groups of 485 communication method and multi-threading are utilized to improve the real-time response time of the system.
The advantages of resilient detection are manifested in the following two aspects. Firstly, by controlling 36 stepping motors descending to different heights, the proposed array-type system can detect the pollutions of complex nuclear surfaces, which is difficult for conventional devices to achieve. Secondly, the XY two-axis motion platform and 'S' type scanning scheme are utilized in the horizontal direction, which can detect a larger range.
However, there are still some places that need to be improved. When testing the runtime overhead, the system was found to run over more time for processing the point cloud, and needs to be further optimized in this regard. Limited by the range of the profilometer, the effective range is trapezoidal. The higher the height, the smaller the horizontal range. If the height difference of the measured surface exceeds ±5 cm, it may cause a measurement dead zone within the set scanning area range. The restoration of the three-dimensional information on the surface could therefore be harmed, resulting in abnormal detection.

Conclusions
In this paper, an array system applied to complex surfaces was constructed for nuclear pollution detection. The measurement model of the two-dimensional laser profilometer was also analyzed. An 'S' scanning mode of profilometer, combined with the laser profilometer and the motion parameters, was used to obtain the three-dimensional terrain information. Due to the limited working area of a single laser profilometer, the system measured the target area in blocks. By the way of point cloud registration, complete point cloud data was obtained. Besides, Euclidean distance median filtering was introduced to solve the impulse noise problem of point cloud. On the surface contour, the rasterized 3D point cloud data was mapped to the lower plane of the pollution detector. Due to the position error caused by the machinery and installation, the edges must be blurred when rasterized. After the rasterization, the maximum value of the point cloud and the initial position of the probe were calculated as the motor closed-loop control input, with the code wheel count of the stepper motor as feedback. Likewise, due to mechanical and installation errors, the initial position of the probe was compensated to guarantee more robustness of the system control. Finally, considering the friendly human-computer interaction, a visual interface with data visual output was also included on the experimental platform. The experiment proved that the proposed methodology may effectively obtain three-dimensional terrain information. The scanning accuracy was 0.35 mm in X-axis and Y-axis, and the measuring accuracy was 12 µm in Z-axis. It also realized control over the movement of the pollution detector accurately, regardless of the maximum control error 0.4 mm. The single measurement area even got up to 240 × 240 cm 2 , and the time for single complete detection was 26.13 s. As a bonus, this new system completed the detection of α and β particles. For future work, we plan to optimize method of point cloud processing for improving the detecting efficiency.