Range Intensity Proﬁles of Multi-Slice Integration for Pulsed Laser Range-Gated Imaging System

: Pulsed laser range-gated imaging (PLRGI) is one of the most effective methods to achieve underwater high-resolution imaging. When searching target, there are two methods can be used: ﬁxed gate and sliding gate. In practice, ﬁxed gate has a small depth of ﬁeld of view and sliding gate cannot meet the need of real time. In order to overcome these problems, multi-slice integration (MSI) method is proposed in this paper. First, the laser energy received by the PLRGI system is derived from radiative transfer theory. In addition, range intensity proﬁles (RIP) of MSI method is established. Experiments are carried out in lab to validate the RIP model and results show that theoretical data and experimental data are in good coincidence. Then the gate width and the number of pulse assigned to each slice are discussed. Finally, to prove the effectiveness of the MSI method, experiments are carried out in a boat tank. Results show that the MSI method is better than the ﬁxed gate method for surveillance, and can acquire a clear image of the target at 18 m (4.5 times of the attenuation length) in the water with attenuation coefﬁcient of 0.25 m − 1 .


Introduction
A rapid increasing attention has been paid to the marine development and there is an an urgent need for underwater high-resolution imaging technology [1][2][3][4][5][6]. At present, there are mainly two methods for underwater imaging [7]. One is acoustic imaging and the other is optical imaging. Acoustic imaging can work far away from the target, but it is difficult to achieve high-resolution image. Optical imaging has high resolution, but it is limited by the underwater environment (rapid decay of light energy, impact of backscattered light, etc.) and works at a closer distance [8]. Thus, special optical imaging techniques are needed to overcome the shortcomings, and the pulsed laser range-gated imaging (PLRGI) is one of the effective methods [9][10][11][12].
By controlling the delay time and the gate width (GW), the PLRGI system enables a slice imaging approach similar to that of computed tomography [13]. The delay time determines position of the slice and the GW determines slice width. Traditional PLRGI system can acquire an image from only one slice in one acquisition. Thus, two different strategies are assumed for target searching: fixed gate and sliding gate. System with a single fixed gate can be treated as a conventional camera. System with sliding gate scans the field of view (FOV) with more than one gate positions. An example is shown in Figure 1, T1 and T2 are two targets in the FOV, and T2 is father away from the PLRGI system than T1. Figure 1a shows the example of fixed gate. In order to obtain all the targets in one image frame, the value of GW must be big enough. However, the output image of T2 is much darker than that of T1 because light energy transmitted underwater decays rapidly with increasing distance. Example for sliding gate is shown in Figure 1b. Sliding gate can remove the problem brought by the fixed gate, but take more time to scan the FOV and cannot meet the need of real time performance.
Photonics 2022, 9, x FOR PEER REVIEW 2 of 12 much darker than that of T1 because light energy transmitted underwater decays rapidly with increasing distance. Example for sliding gate is shown in Figure 1b. Sliding gate can remove the problem brought by the fixed gate, but take more time to scan the FOV and cannot meet the need of real time performance.
(a) (b) To overcome the problem for searching target, there are two methods can be used. One is multi-channel receivers, which can acquire images from multi-slice at the same time [2,14]. However, this will make the system too complex and expensive. The other one is multi-pulse integration (MPI) method, which split the detection area into multi-areas and assigned pulses to detect each area alone and combined them to a single image [15,16]. However, the model for MPI method is not complete and the number of pulses allocated to each slice is the same in the literature. In addition, this makes it similar to a single slice and suffers from the problem mentioned above.
Inspired by the MPI method, a multi-slice integration (MSI) method is proposed in this work. Different from the MPI method, number of the pulses assigned to each slice and GW of each slice can be different from each other in the proposed MSI method. The paper is organized as follows: following this Introduction section, principles of the traditional PLRGI and MSI method are described and the range intensity profiles (RIP) of the MSI method is established by radiative transfer theory in Section 2. The number of pulse assigned to each slice is determined by the rules that peak energy of each slice are kept as the same. In addition, the best value of GW for every slice is discussed as well. In Section 3, experiments are designed and conducted, and the experimental results are presented and discussed. Conclusions are drawn in the last section.

Principle of PLRGI
Principle of PLRGI is shown in Figure 2. PLRGI system has three main components: a laser, a receiver (camera with a selective gate), and a synchronous controller. When the laser emits a pulse, the controller records the moment and calculates the delay time for the laser to reach the target and return. Before the pulse reaches the camera, the gate is close. In addition, when the pulse reaches the camera, the gate is opened. The opening time for the gate corresponds to the depth of the slice. In this way, the PLRGI system can remove most of the backscattered light. To overcome the problem for searching target, there are two methods can be used. One is multi-channel receivers, which can acquire images from multi-slice at the same time [2,14]. However, this will make the system too complex and expensive. The other one is multi-pulse integration (MPI) method, which split the detection area into multi-areas and assigned pulses to detect each area alone and combined them to a single image [15,16]. However, the model for MPI method is not complete and the number of pulses allocated to each slice is the same in the literature. In addition, this makes it similar to a single slice and suffers from the problem mentioned above.
Inspired by the MPI method, a multi-slice integration (MSI) method is proposed in this work. Different from the MPI method, number of the pulses assigned to each slice and GW of each slice can be different from each other in the proposed MSI method. The paper is organized as follows: following this Introduction section, principles of the traditional PLRGI and MSI method are described and the range intensity profiles (RIP) of the MSI method is established by radiative transfer theory in Section 2. The number of pulse assigned to each slice is determined by the rules that peak energy of each slice are kept as the same. In addition, the best value of GW for every slice is discussed as well. In Section 3, experiments are designed and conducted, and the experimental results are presented and discussed. Conclusions are drawn in the last section.

Principle of PLRGI
Principle of PLRGI is shown in Figure 2. PLRGI system has three main components: a laser, a receiver (camera with a selective gate), and a synchronous controller. When the laser emits a pulse, the controller records the moment and calculates the delay time for the laser to reach the target and return. Before the pulse reaches the camera, the gate is close. In addition, when the pulse reaches the camera, the gate is opened. The opening time for the gate corresponds to the depth of the slice. In this way, the PLRGI system can remove most of the backscattered light.
PLRGI system enables a slice imaging approach similar to that of computed tomography. In addition, information of the slice depends on the delay time, gate width (GW) and the number of pulse assigned to the slice. Position of the slice depends on the delay time and GW. Intensity of the same slice depends on the number of pulse assigned to the slice. Starting position of the slice is calculated by the TOF of the emitted laser pulse, which is written as, where c w is the light speed in water and t 0 is the delay time. Let duration of the gate be τ, then width of the slice w is calculated as, Photonics 2022, 9, x FOR PEER REVIEW 3 of 12 (a) (b) PLRGI system enables a slice imaging approach similar to that of computed tomography. In addition, information of the slice depends on the delay time, gate width (GW) and the number of pulse assigned to the slice. Position of the slice depends on the delay time and GW. Intensity of the same slice depends on the number of pulse assigned to the slice. Starting position of the slice is calculated by the TOF of the emitted laser pulse, which is written as, where cw is the light speed in water and t0 is the delay time.
Let duration of the gate be τ, then width of the slice w is calculated as, (2)

Multi-Slice Integration Method
Three slices are taken into consideration to present the principle of the MSI method as shown in Figure 3. In the system, a laser with a high pulse repetition frequency is utilized. The laser pulses in a video frame can be split into three groups, and each group corresponds to one slice. Each slice has its delay time and GW.
The time sequence of the MSI method is shown in Figure 3b. The number of pulse assigned to target ①, ② and ③ are n1, n2 and n3, respectively. Delay time of three groups are t1, t2 and t3, which is different from each other. Each pulse has a range intensity profiles (RIP), and pulses belonging to the same group have the same RIPs. Therefore, three sub-RIPs can be integrated into a video frame. With the MSI method, all of the targets in the FOV can be obtained in a video image frame in real time. (a)

Multi-Slice Integration Method
Three slices are taken into consideration to present the principle of the MSI method as shown in Figure 3. In the system, a laser with a high pulse repetition frequency is utilized. The laser pulses in a video frame can be split into three groups, and each group corresponds to one slice. Each slice has its delay time and GW.
PLRGI system enables a slice imaging approach similar to that of computed tomography. In addition, information of the slice depends on the delay time, gate width (GW) and the number of pulse assigned to the slice. Position of the slice depends on the delay time and GW. Intensity of the same slice depends on the number of pulse assigned to the slice. Starting position of the slice is calculated by the TOF of the emitted laser pulse, which is written as, where cw is the light speed in water and t0 is the delay time.
Let duration of the gate be τ, then width of the slice w is calculated as, (2)

Multi-Slice Integration Method
Three slices are taken into consideration to present the principle of the MSI method as shown in Figure 3. In the system, a laser with a high pulse repetition frequency is utilized. The laser pulses in a video frame can be split into three groups, and each group corresponds to one slice. Each slice has its delay time and GW.
The time sequence of the MSI method is shown in Figure 3b. The number of pulse assigned to target ①, ② and ③ are n1, n2 and n3, respectively. Delay time of three groups are t1, t2 and t3, which is different from each other. Each pulse has a range intensity profiles (RIP), and pulses belonging to the same group have the same RIPs. Therefore, three sub-RIPs can be integrated into a video frame. With the MSI method, all of the targets in the FOV can be obtained in a video image frame in real time. Take the slice which contains target ① as an example to calculate the RIP of a single pulse. The coordinate system established is shown in Figure 4. Here α is the half divergence angle of the laser beam. The laser is assumed to emit a pulse of energy E0 within a short period τP. Speed of the laser in water is cw. Target is placed at position z away from the imaging system. In this work, multiple scattering and the influence of ambient light are not considered. The time sequence of the MSI method is shown in Figure 3b. The number of pulse assigned to target 1 , 2 and 3 are n 1 , n 2 and n 3 , respectively. Delay time of three groups are t 1 , t 2 and t 3 , which is different from each other. Each pulse has a range intensity profiles (RIP), and pulses belonging to the same group have the same RIPs. Therefore, three sub-RIPs can be integrated into a video frame. With the MSI method, all of the targets in the FOV can be obtained in a video image frame in real time.
Take the slice which contains target 1 as an example to calculate the RIP of a single pulse. The coordinate system established is shown in Figure 4. Here α is the half divergence angle of the laser beam. The laser is assumed to emit a pulse of energy E 0 within a short period τ P . Speed of the laser in water is c w . Target is placed at position z away from the imaging system. In this work, multiple scattering and the influence of ambient light are not considered. Take the slice which contains target ① as an example to calculate the RIP of a single pulse. The coordinate system established is shown in Figure 4. Here α is the half divergence angle of the laser beam. The laser is assumed to emit a pulse of energy E0 within a short period τP. Speed of the laser in water is cw. Target is placed at position z away from the imaging system. In this work, multiple scattering and the influence of ambient light are not considered. Light intensity within the divergence angle is assumed to be uniformly distributed, then the illumination on the target is: where c is the water attenuation coefficient. Let S be a pixel in the FOV at range z. Area of the element S is: where f is the focal length and S0 is the area of a single pixel of the camera. Radiant flux received by element S is, Light intensity within the divergence angle is assumed to be uniformly distributed, then the illumination on the target is: where c is the water attenuation coefficient. Let S be a pixel in the FOV at range z. Area of the element S is: where f is the focal length and S 0 is the area of a single pixel of the camera. Radiant flux received by element S is, Assuming that the target is a Lambert reflector, ρ is the reflectivity of the target. According to the Lambert's law, the intensity of reflection in the normal direction of element S is, Under the assumption of paraxial small angle transmission, the radiant flux reaching the receiver is: where ∆Ω ≈ πD 2 4z 2 and D is the diameter of the receiver. By substituting Equation (6) into Equation (7), we can acquire, where H = ρD 2 E 0 S 0 4πτ p f 2 (tan α) 2 . According to Jens Busck [17], The amount of light energy received by camera can be expressed as, where Φ is the laser pulse radiant flux that reach the camera, G the gate function, z distance of the target, c w the light speed in water and t 1 the delay time of camera gate. Let τ be the opening time of the camera gate, then the gate function can be treated as rectangular function, Let z 1 = c w t 1 /2 be the start position of camera gate, then, where z' = z − z 1 . By integrating from 0 to τ, a single pulse energy received by camera is, Thus, the RIP data of a single pulse can be obtained from Equation (12). In order to verify Equation (12), experiment is carried out in lab. Water attenuation coefficient in the water tank is 0.32/m. The range of interest is from 3 m to 4 m. Delay time t 1 and gate time τ are set as 26.7 ns (z 1 = 3 m) and 8.9 ns, respectively. RIP data comparison between theory and experiment is shown in Figure 5a. From Figure 5a, we can see that data in theory and experiment are in good coincidence. In addition, slice width, influenced by the laser pulse width, is expanded from 1m to 1.5 m. In the MSI method, the laser energy received by camera is, where k means kth slice and Ek can be acquired from Equation (12). According to Equations (12) and (13), three parameters, delay time, number of pulse and GW for each slice, need to be set to acquire the RIP data with the MSI method.

Method to Determine Parameters in MSI Method
Assuming that the range of interest is from z1 to z0, then delay time t1 is, From Figure 5a, we can see that data in theory and experiment are in good coincidence. In addition, slice width, influenced by the laser pulse width, is expanded from 1m to 1.5 m. In the MSI method, the laser energy received by camera is, where k means kth slice and E k can be acquired from Equation (12). According to Equations (12) and (13), three parameters, delay time, number of pulse and GW for each slice, need to be set to acquire the RIP data with the MSI method.

Method to Determine Parameters in MSI Method
Assuming that the range of interest is from z 1 to z 0 , then delay time t 1 is, Delay time t 2 and t 3 are, For the number of pulse, we follow the rules that peak energy of each slice are kept as the same. From Equation (12), it can be seen that the peak energy of each slice appear at the distance z 1 , z 2 and z 3 , whose delay time corresponding to be t 1 , t 2 and t 3 . Take the 1st and the 2nd slice as an example, peak energy of 1st slice is, Peak energy of 2nd slice is, To keep peak energy of the 1st and the 2nd slice the same, we can acquire, Then the proportion of n 1 , n 2 is, Another experiment is carried out in lab to verify the pulse allocation strategy. The range of interest is from 3 m to 6 m. According to Equations (14)-(16), delay time t 1 , t 2 and t 3 are set as 26.7 ns, 35.6 ns and 44.4 ns, respectively. Gate time τ for all slices is 8.9 ns. Image of the 1st slice is appropriate when the number of pulse is 28. Following the peak energy rule, the number of pulse n 2 and n 3 are set to be 95 and 280. The difference between RIP data in theory and experiment is shown in Figure 5b and they are also in good coincidence. Moreover, comparing to a single slice, the MSI method can extend the maximal range. It can be seen from Figure 5 that the maximal range is extended from 1.5 m to 3.5 m.
Another parameter need to be set is GW. Light energy received by the system at maximal range is chosen to measure the performance of GW. In simulation, delay time t 1 , t 2 and t 3 are set as 26.7 ns, 35.6 ns and 44.4 ns, respectively. Gate time τ is changed from 8.9 ns to 44.4 ns. In addition, intensity of the image varies with the change of τ. Tendency of the intensity is shown in Figure 6.
cidence. Moreover, comparing to a single slice, the MSI method can extend the maximal range. It can be seen from Figure 5 that the maximal range is extended from 1.5 m to 3.5 m.
Another parameter need to be set is GW. Light energy received by the system at maximal range is chosen to measure the performance of GW. In simulation, delay time t1, t2 and t3 are set as 26.7 ns, 35.6 ns and 44.4 ns, respectively. Gate time τ is changed from 8.9 ns to 44.4 ns. In addition, intensity of the image varies with the change of τ. Tendency of the intensity is shown in Figure 6. From Figure 6, we can see that intensity increases with GW and the maximum intensity appears firstly when GW is the sum of the range of interest and the laser pulse width. Thus, ideal GW for the MSI method is the sum of the range of interest and the laser pulse width.

Experiment Setup
As shown in Figure 7, PLRGI system used in the experiment includes four parts: a pulsed laser, an ICMOS camera, a synchronous controller and an image processor. The laser works at a repetition rate of 10 kHz with a pulse length of 5 ns. In addition, the laser divergence angle can be adjusted when needed. The minimum gate width of the ICMOS is 5 ns. Furthermore, the maximum frame rate is 30 Hz. A FPGA is used for synchronous controlling and a NVIDIA TX2 core module for image processing. The system is controlled by a computer to set its working parameters and to show the imaging result. More information about the system can be acquired from literature [18]. From Figure 6, we can see that intensity increases with GW and the maximum intensity appears firstly when GW is the sum of the range of interest and the laser pulse width. Thus, ideal GW for the MSI method is the sum of the range of interest and the laser pulse width.

Experiment Setup
As shown in Figure 7, PLRGI system used in the experiment includes four parts: a pulsed laser, an ICMOS camera, a synchronous controller and an image processor. The laser works at a repetition rate of 10 kHz with a pulse length of 5 ns. In addition, the laser divergence angle can be adjusted when needed. The minimum gate width of the ICMOS is 5 ns. Furthermore, the maximum frame rate is 30 Hz. A FPGA is used for synchronous controlling and a NVIDIA TX2 core module for image processing. The system is controlled by a computer to set its working parameters and to show the imaging result. More information about the system can be acquired from literature [18]. Experiments are carried out in a boat tank, which size is 175 m × 6 m × 4 m (length × width × depth). The wall of the boat tank is rough and the reflectivity of light is low. The PLRGI system is placed at 1m underwater as shown in Figure 8.  Experiments are carried out in a boat tank, which size is 175 m × 6 m × 4 m (length × width × depth). The wall of the boat tank is rough and the reflectivity of light is low. The PLRGI system is placed at 1m underwater as shown in Figure 8. Experiments are carried out in a boat tank, which size is 175 m × 6 m × 4 m (length × width × depth). The wall of the boat tank is rough and the reflectivity of light is low. The PLRGI system is placed at 1m underwater as shown in Figure 8.

Comparison of Single Slice and Multi-Slice
Two targets, a resolution target and a checkerboard target, are arranged at 10 m and 15 m directly in front of the system. Size of the resolution target and the checkerboard target are 1 m × 1 m, 0.5 m × 0.5 m, respectively. The albedo for the two regions (black and white) on both targets are 10% and 90%. Range from 10 m to 15 m is the region of interest. Delay time of single slice is 88 ns and GW is 50 ns. For MSI method, the region of interest is split into 3 slices. In addition, delay time of the slices are 88 ns, 103 ns and 117 ns, respectively. GW of all the slices is the same. Furthermore, typical value 25 ns, 43.75 ns and 50 ns are adopted to measure performance of the system. Water attenuation coefficient is estimated to 0.25/m by the method we proposed in literature [19]. So the number of pulse n1, n2 and n3 are set to be 36, 108 and 310. Output of the single slice and MSI method are shown in Figure 9.

Comparison of Single Slice and Multi-Slice
Two targets, a resolution target and a checkerboard target, are arranged at 10 m and 15 m directly in front of the system. Size of the resolution target and the checkerboard target are 1 m × 1 m, 0.5 m × 0.5 m, respectively. The albedo for the two regions (black and white) on both targets are 10% and 90%. Range from 10 m to 15 m is the region of interest. Delay time of single slice is 88 ns and GW is 50 ns. For MSI method, the region of interest is split into 3 slices. In addition, delay time of the slices are 88 ns, 103 ns and 117 ns, respectively. GW of all the slices is the same. Furthermore, typical value 25 ns, 43.75 ns and 50 ns are adopted to measure performance of the system. Water attenuation coefficient is estimated to 0.25/m by the method we proposed in literature [19]. So the number of pulse n 1 , n 2 and n 3 are set to be 36, 108 and 310. Output of the single slice and MSI method are shown in Figure 9. Experiments are carried out in a boat tank, which size is 175 m × 6 m × 4 m (length × width × depth). The wall of the boat tank is rough and the reflectivity of light is low. The PLRGI system is placed at 1m underwater as shown in Figure 8.

Comparison of Single Slice and Multi-Slice
Two targets, a resolution target and a checkerboard target, are arranged at 10 m and 15 m directly in front of the system. Size of the resolution target and the checkerboard target are 1 m × 1 m, 0.5 m × 0.5 m, respectively. The albedo for the two regions (black and white) on both targets are 10% and 90%. Range from 10 m to 15 m is the region of interest. Delay time of single slice is 88 ns and GW is 50 ns. For MSI method, the region of interest is split into 3 slices. In addition, delay time of the slices are 88 ns, 103 ns and 117 ns, respectively. GW of all the slices is the same. Furthermore, typical value 25 ns, 43.75 ns and 50 ns are adopted to measure performance of the system. Water attenuation coefficient is estimated to 0.25/m by the method we proposed in literature [19]. So the number of pulse n1, n2 and n3 are set to be 36, 108 and 310. Output of the single slice and MSI method are shown in Figure 9. (c) Output of multi-slice, whose GW is 25 ns; (d) Output of multi-slice, whose GW is 43.75 ns; (e) Output of multi-slice, whose GW is 50 ns; (f) Output of multi-slice, whose GW is 55 ns.
Output of single slice shows that target 2 is much darker than that of target 1 because light decays rapidly underwater with increasing distance. Output of MSI method are much better than single slice. Furthermore, in all three MSI method, GW of 50 ns is the best.
To evaluate the quality of the image, intensity and gradient are used as metrics in this paper. The gradient is a metric used to indicate sharpness of an image, and the larger its value, the higher sharpness of the image obtained. We can acquire gradient G from Equation (21), Figure 9. Performance of single slice and multi-slice. (a) Two targets in experiment; (b) Output of single slice; (c) Output of multi-slice, whose GW is 25 ns; (d) Output of multi-slice, whose GW is 43.75 ns; (e) Output of multi-slice, whose GW is 50 ns; (f) Output of multi-slice, whose GW is 55 ns.
Output of single slice shows that target 2 is much darker than that of target 1 because light decays rapidly underwater with increasing distance. Output of MSI method are much better than single slice. Furthermore, in all three MSI method, GW of 50 ns is the best.
To evaluate the quality of the image, intensity and gradient are used as metrics in this paper. The gradient is a metric used to indicate sharpness of an image, and the larger its value, the higher sharpness of the image obtained. We can acquire gradient G from Equation (21), where M, N represent the number of columns and rows of the image. Metric of intensity is the mean intensity of the area within red dotted lines. It can be expressed as, Quantified comparisons between Figure 9c-f are shown in Figure 10. best.
To evaluate the quality of the image, intensity and gradient are used as metric paper. The gradient is a metric used to indicate sharpness of an image, and the larger the higher sharpness of the image obtained. We can acquire gradient G from Equation where M, N represent the number of columns and rows of the image.
Metric of intensity is the mean intensity of the area within red dotted lines. expressed as,  It can be seen from Figure 10 that performance of intensity and gradient incre GW. When GW larger than 50 ns, which is almost the sum of range of interest a pulse width, there are not much change in intensity and gradient. This is in goo ment with conclusion in Section 2.

Performance of MSI Method
To test the performance of MSI method, another experiment is carried out. A in Figure 11, we placed 6 targets within the FOV, whose distance are 4 m, 6 m, 9 It can be seen from Figure 10 that performance of intensity and gradient increase with GW. When GW larger than 50 ns, which is almost the sum of range of interest and laser pulse width, there are not much change in intensity and gradient. This is in good agreement with conclusion in Section 2.

Performance of MSI Method
To test the performance of MSI method, another experiment is carried out. As shown in Figure 11, we placed 6 targets within the FOV, whose distance are 4 m, 6 m, 9 m, 12 m, 15 m and 18 m, respectively. There are three different sizes of targets: size of 1 is 0.1 m × 0.8 m, size of 2 is 0.08 m × 0.8 m and size of 3 is 0.05 m × 0.8 m. Two targets of 3 are placed at 4 m and 6 m. Two targets of 2 are placed at 9 m and 12 m. Two targets of 1 are placed at 15 m and 18 m. The albedo for the two regions (black and white) on all targets are 10% and 90%.
In the experiment, 6 slices are used in MSI method. Delay time of 6 slices are 26.7 ns, 53.3 ns, 80.1 ns, 106.6 ns, 133.3 ns and 160 ns. GW of all the slices are 160 ns. According to Section 2, the number of pulse assigned to each slice are 1, 4, 36, 288, 2000 and 13,018. However, total number of pulse can be assigned is 2500 at video frame of 4 fps. So the number of the pulse assigned to each slice are 1, 4, 36, 288, 2000 and 171.
As shown in Figure 11, all the targets can be seen in the output image and the farthest target at 18 m is distinguishable. Thus, the maximal range of identification for MSI method can reach at least 4.5 attenuation length while that is 2-3 attenuation length for conventional camera system [8,20,21]. are placed at 15 m and 18 m. The albedo for the two regions (black and white) on all targets are 10% and 90%.
In the experiment, 6 slices are used in MSI method. Delay time of 6 slices are 26.7 ns, 53.3 ns, 80.1 ns, 106.6 ns, 133.3 ns and 160 ns. GW of all the slices are 160 ns. According to Section 2, the number of pulse assigned to each slice are 1, 4, 36, 288, 2000 and 13,018. However, total number of pulse can be assigned is 2500 at video frame of 4 fps. So the number of the pulse assigned to each slice are 1, 4, 36, 288, 2000 and 171.