# Measuring the Angular Velocity of a Propeller with Video Camera Using Electronic Rolling Shutter

Authors: Yipeng Zhao; Yongbin Li; Shijie Guo; Tiejun Li
Journal: Journal of Sensors
 (2018)

---

## Abstract

Noncontact measurement for rotational motion has advantages over the traditional method which measures rotational motion by means of installing some devices on the object, such as a rotary encoder. Cameras can be employed as remote monitoring or inspecting sensors to measure the angular velocity of a propeller because of their commonplace availability, simplicity, and potentially low cost. A defect of the measurement with cameras is to process the massive data generated by cameras. In order to reduce the collected data from the camera, a camera using ERS (electronic rolling shutter) is applied to measure angular velocities which are higher than the speed of the camera. The effect of rolling shutter can induce geometric distortion in the image, when the propeller rotates during capturing an image. In order to reveal the relationship between the angular velocity and the image distortion, a rotation model has been established. The proposed method was applied to measure the angular velocities of the two-blade propeller and the multiblade propeller. The experimental results showed that this method could detect the angular velocities which were higher than the camera speed, and the accuracy was acceptable.


---

## Body

## 1. Introduction

Rotation is one of basic motions, which is common in machines like motors, gears, and other wheels. The rotation should be under control to keep machines in good operation, and a lot of mechanical failures are caused by the rotary movement. So it is important to measure the angular velocity. To measure angular velocities, contact-type sensors are widely used, such as mechanical tachometers, optical tachometers, photoelectric encoders, and optical encoders [1]. These methods are usually based on mechanical contact, and as a result, they are easily affected by the rotation of the target or the small target inertia. In the past twenty years, noncontact methods have been developed like tomography, ultrasound, laser, and computer vision [2]. The advanced sensors can overcome the defects of contact-type sensors, and computer vision could be more widely used compared with the other noncontact sensors.Over the past decade, some researchers have focused on angular measurements based on computer vision. Wang et al. [3] measured motor angular velocities with blur images which had motion information. Angular velocities could be extracted in those motion-blur images in polar coordinates. Ait-Aider et al. [4, 5] obtained the object pose and velocity by an electronic rolling shutter (ERS) camera. His method was based on the assumption that all the lines in the real world were straight and those lines warped in the image captured by the ERS camera. Magerand et al. [6] measured the object pose and motion from a single ERS image with automatic 2D-3D matching. He and Wei [7] measured a shaft’s velocity by using the ERS camera. Zhu and Yu [2] measured angular velocities of the object by Hough transform.Measuring angular velocity with computer vision also has some defects, such as the environmental air quality and time-consuming process. There are two ways to improve the processing speed. One way is to reduce the resolution of the camera, but this way is not consistent with the development trend of cameras [8, 9]. The other way is to reduce the speed of the camera [10]. In other words, it is measuring a high-speed rotation with a low-speed camera. But when the angular velocity is much higher than the camera’s speed, it is difficult to find out the angle for the ambiguity or the distortion in the image. There are few studies on measuring the high-speed rotation with a low-speed camera. The row speed of an ERS camera is much higher than the speed of the camera. This special property would help measure the high-speed rotation of some objects, such as wheels or propellers which have symmetrical structures. Those objects are very common in machines and easy to be analyzed to extract angular velocities. The previous researches did not attempt to measure the high-speed rotation by an ERS camera [4–7].In this paper, an ERS camera works as a sensor to measure the angular velocity of the propeller which rotates faster than the speed of the camera. To measure the angular velocity, a simulation is created to demonstrate through images that the propeller rotates at different speeds. According to the geometrical features in the images, an algorithm is proposed to calculate the angular velocity of the propeller. Experiments are also demonstrated to test and verify the proposed method for the measurement of angular velocity. The rest of the paper is organized as follows. In Section2, the working principle of ERS is presented, and an algorism is proposed and verified on a simulation. In Section 3, two experiments are established to verify the proposed method in the actual environment. Finally, Section 4 concludes the paper.


## 2. Method

In the field of digital signal processing, the Nyquist–Shannon sampling theorem establishes a sufficient condition for a sample rate that permits a discrete sequence of samples to capture all the information from a continuous-time signal of finite bandwidth. A sufficient sample rate is at least double the highest frequency to be sampled; moreover, in practice, the sample rate is usually fourfold. To measure angular velocity by the camera, there is massive data to be processed. In order to decrease the data generated by the camera, a low-speed camera with ERS could be used. The ERS camera has a special feature that records images row by row, and this feature could increase the speed of the camera in some way.


### 2.1. Working Principle of ERS

The ERS camera is a type of CMOS image sensor and it is very common on cellphones. The amount of signal generated by the image sensor depends on the amount of light that falls on the images, in terms of exposure.Therefore, an on-chip electronic shutter is required to control the intensity and duration of exposure. The shutter of CMOS has two types, global shutter and rolling shutter. With global shutter image sensors, every row of the pixels in the image is exposed at the same time. So there are no motion artifacts in the resulting image. With rolling shutter image sensors, the rows in the image are exposed in sequence starting at the top and proceeding row by row to the bottom. For each row in the image, the time of integration and row delay is fixed, leading to a uniform staggered time across the frame. When the motion direction of the object is orthogonal to the row direction of the image, visual artifacts will arise in the image captured. The working principle of ERS is demonstrated in Figure1. Since the image sensor effectively integrates each row of the pixel array at a different point in time, the static lamppost, which is vertical to the expressway, becomes inclined as shown in Figure 1(b). With such characteristics of the ERS camera, the low-speed camera can record more information of the fast-moving object. Although some works applied this property to the measurement of linear motion or vibration [7, 11], most works tried to eliminate this effect on the image [12–14]. In this paper, this property is used for the measurement of the angular velocity of the propeller.Figure 1
Working principle of ERS: (a) time sequence of row delay and integration and (b) lamppost shot by the ERS camera in a car with a speed about 90 km/h.
(a)
(b)


### 2.2. Background Segmentation

Before the calculation of rotating velocity, the propeller should be segmented from the background. Gaussian mixture model (GMM) [15] is used as an unsupervised image change detection to extract the propeller. Expectation maximization (EM) [16] is a popular technique used to determine the parameters of a mixture with an a priori given number of components. The EM algorithm provides a particular way that implements the maximum likelihood estimation for the parameter in GMM. But the M-step of the EM algorithm cannot evaluate the prior distribution in a closed form because of the complexity of the maximum likelihood estimation. Therefore, for each iteration in the EM algorithm, a reparatory projection step [17] is applied to the M-step for the purpose that the prior probabilities are positive and sum to one.The segmentation by GMM has many defects like holes and noises in the image. Before the segmentation, the images are denoised with the method proposed by Xu et al. [18]. The rotating calculation closely depends on the geometric characteristics of the propeller segmented from the background. Those defects on the propeller would cause errors on the result of the calculation. The holes should be filled and the little isolate area would be erased.


### 2.3. Calculation of Angular Velocity

A propeller is selected as the measurand as shown in Figure2. This propeller has two isolated blades and its initial angle is dictated by θ. In order to demonstrate the patterns recorded by the ERS camera, a simulation is illustrated in Figure 3 (in this paper, all the units of angular velocity are r/s). The propeller is at the center of the single image whose size is 1000 × 1000. The exposure sequence of the ERS camera is from top to bottom. The time of row delay Td is 996 μs, and the exposure time of every row is 4996 μs. There is no frame delay, so the whole image time is one second. The propeller rotates clockwise at different speeds as shown in Figure 3. The initial angle of the first-row images in Figure 3 is 0.25π; the initial angle of the second-row images in Figure 3 is 0.5π. To extract the angular velocity from a single image, the procedure includes two parts: searching the center of the propeller and extracting the rotating angle. All the operations are based on an assumption that the two-blade propeller has central symmetry.Figure 2
Two-blade propeller.Figure 3
Simulation of a two-blade propeller: (a) 0.5 r/s, (b) 1 r/s, (c) 2 r/s, (i)θ=0.25π, and (ii) θ=0.5π.The points on the different blade edges cluster into different sets:Hi,i=1,2,…,K, where K is the quantity of the isolated blade. Point Pc is picked out at random in the image. The distance from Pc to the point in Hi is
(1)DiPcPJ→Pj∈Hi.Letdi=maxDi, where max() is a mathematical function to find the maximum values of its input. If Pc can satisfy the condition in (2), the center is Pc.
(2)d1=d2=⋯=dK.But this method for center calculation is sensitive to noise and it is ambiguous in some situations as shown in Figure3(iic, iid). So this method could not be used for finding the center, but it can work as a quick search and narrow the area of the rotating center. Let d¯ be the expectation ofdi:
(3)d¯=1K∑i=1Kdi.The rotating center belongs to the area(4)Hc=Pc1/K∑i=1Kdi−d¯2d¯<δ,where δ is a threshold belonging to the range of (0, 0.1). When δ=0.02, the areas of color in Figure 4 are the range of the center. To find the precise center from Hc, the angular velocity should be considered.Figure 4
Searching the center of the propeller: (a)–(c) correspond to Figures3(a)–3(c), and (i) and (ii) correspond to (i) and (ii) in Figure 3.For the purpose of a better solution, a scanning method is proposed. PointPc is picked from Hc, and it works as the center of a circle whose radius is r=jd¯,0.5≤j≤0.9. This circle scans the blade anticlockwise, and its initial position is at its left horizontal axis. There are rising edges and falling edges on the circle. The points on the rising edges are denoted by P=P1,P2,…,Pm, where m is the quantity of the rising edges on the circle. Meanwhile, the points on the falling edges are denoted by P′=P1′,P2′,…,Pm′. The ordinates of these edge points are Y=y1,y2,…,ym and Y′=y1′,y2′,…,ym′, respectively. The angles between the initial position and these edge points are Φ=ϕ1,ϕ2,…,ϕm and Φ′=ϕ1′,ϕ2′,…,ϕm′, respectively. Let ya=maxY and yb=minY, where min() is a mathematical function to find the minimum values of its input. The propeller rotates from yb to ya as
(5)φ=2πa−bn−ϕa−ϕb,s.t.ϕa−ϕba−b≥2π−ϕa+ϕbm−a+b,2πa−bn−2π−ϕa+ϕb,s.t.ϕa−ϕba−b<2π−ϕa+ϕbm−a+b,where n is the quantity of blades and its value is two in the simulation. When the value of φ is from the first class in (5), the rotational direction is clockwise. When the value of φ is from the second class in (5), the rotational direction is anticlockwise. The angular velocity ω could be calculated out by the points on the rising edges.
(6)ω=f1r=φTdya−yb,where Td is the row delay time of the ERS camera. The points on the falling edges can also be used on the calculation of angular velocity:
(7)ω′=f2r=φ′Tdyc−yd,where yc=maxY′, yd=maxY′, and φ′ is similar to φ. Let ec be the accumulative error corresponding to radius r at point Pc.
(8)ec=∫0.5d0.9d¯rf1r−f2rdr.LetEc=ec∣ec=gPc,Pc∈Hc, where g is a function Pc⟶ec. The value of Ec of different points are normalized and indexed by a color map as shown in Figure 4. Let gPc^=minEc. Point Pc^ is the center of the propeller. In order to promote the accuracy of the calculation, the optimized value of angular velocity ω^ is weighted by r.
(9)ω^=∫0.5d¯0.9d¯rf1r+f2rdr2∫0.5d¯0.9d¯rdr.The result of the simulation is shown in Table1. The actual center position is (500,500), and the most deviation of the center of the propeller is within 5 pixels. The relative errors of angular velocity are no more than 2%.Table 1
Result of the simulation.
Figure
Measured center position (pixel)
Actual velocity (r/s)
Measured velocity (r/s)
(ia)
(500, 500)
0.5
0.510
(ib)
(500, 502)
1
1.006
(ic)
(497, 500)
2
2.003
(iia)
(501, 499)
0.5
0.508
(iib)
(497, 502)
1
0.999
(iic)
(503, 502)
2
2.012


## 2.1. Working Principle of ERS

The ERS camera is a type of CMOS image sensor and it is very common on cellphones. The amount of signal generated by the image sensor depends on the amount of light that falls on the images, in terms of exposure.Therefore, an on-chip electronic shutter is required to control the intensity and duration of exposure. The shutter of CMOS has two types, global shutter and rolling shutter. With global shutter image sensors, every row of the pixels in the image is exposed at the same time. So there are no motion artifacts in the resulting image. With rolling shutter image sensors, the rows in the image are exposed in sequence starting at the top and proceeding row by row to the bottom. For each row in the image, the time of integration and row delay is fixed, leading to a uniform staggered time across the frame. When the motion direction of the object is orthogonal to the row direction of the image, visual artifacts will arise in the image captured. The working principle of ERS is demonstrated in Figure1. Since the image sensor effectively integrates each row of the pixel array at a different point in time, the static lamppost, which is vertical to the expressway, becomes inclined as shown in Figure 1(b). With such characteristics of the ERS camera, the low-speed camera can record more information of the fast-moving object. Although some works applied this property to the measurement of linear motion or vibration [7, 11], most works tried to eliminate this effect on the image [12–14]. In this paper, this property is used for the measurement of the angular velocity of the propeller.Figure 1
Working principle of ERS: (a) time sequence of row delay and integration and (b) lamppost shot by the ERS camera in a car with a speed about 90 km/h.
(a)
(b)


## 2.2. Background Segmentation

Before the calculation of rotating velocity, the propeller should be segmented from the background. Gaussian mixture model (GMM) [15] is used as an unsupervised image change detection to extract the propeller. Expectation maximization (EM) [16] is a popular technique used to determine the parameters of a mixture with an a priori given number of components. The EM algorithm provides a particular way that implements the maximum likelihood estimation for the parameter in GMM. But the M-step of the EM algorithm cannot evaluate the prior distribution in a closed form because of the complexity of the maximum likelihood estimation. Therefore, for each iteration in the EM algorithm, a reparatory projection step [17] is applied to the M-step for the purpose that the prior probabilities are positive and sum to one.The segmentation by GMM has many defects like holes and noises in the image. Before the segmentation, the images are denoised with the method proposed by Xu et al. [18]. The rotating calculation closely depends on the geometric characteristics of the propeller segmented from the background. Those defects on the propeller would cause errors on the result of the calculation. The holes should be filled and the little isolate area would be erased.


## 2.3. Calculation of Angular Velocity

A propeller is selected as the measurand as shown in Figure2. This propeller has two isolated blades and its initial angle is dictated by θ. In order to demonstrate the patterns recorded by the ERS camera, a simulation is illustrated in Figure 3 (in this paper, all the units of angular velocity are r/s). The propeller is at the center of the single image whose size is 1000 × 1000. The exposure sequence of the ERS camera is from top to bottom. The time of row delay Td is 996 μs, and the exposure time of every row is 4996 μs. There is no frame delay, so the whole image time is one second. The propeller rotates clockwise at different speeds as shown in Figure 3. The initial angle of the first-row images in Figure 3 is 0.25π; the initial angle of the second-row images in Figure 3 is 0.5π. To extract the angular velocity from a single image, the procedure includes two parts: searching the center of the propeller and extracting the rotating angle. All the operations are based on an assumption that the two-blade propeller has central symmetry.Figure 2
Two-blade propeller.Figure 3
Simulation of a two-blade propeller: (a) 0.5 r/s, (b) 1 r/s, (c) 2 r/s, (i)θ=0.25π, and (ii) θ=0.5π.The points on the different blade edges cluster into different sets:Hi,i=1,2,…,K, where K is the quantity of the isolated blade. Point Pc is picked out at random in the image. The distance from Pc to the point in Hi is
(1)DiPcPJ→Pj∈Hi.Letdi=maxDi, where max() is a mathematical function to find the maximum values of its input. If Pc can satisfy the condition in (2), the center is Pc.
(2)d1=d2=⋯=dK.But this method for center calculation is sensitive to noise and it is ambiguous in some situations as shown in Figure3(iic, iid). So this method could not be used for finding the center, but it can work as a quick search and narrow the area of the rotating center. Let d¯ be the expectation ofdi:
(3)d¯=1K∑i=1Kdi.The rotating center belongs to the area(4)Hc=Pc1/K∑i=1Kdi−d¯2d¯<δ,where δ is a threshold belonging to the range of (0, 0.1). When δ=0.02, the areas of color in Figure 4 are the range of the center. To find the precise center from Hc, the angular velocity should be considered.Figure 4
Searching the center of the propeller: (a)–(c) correspond to Figures3(a)–3(c), and (i) and (ii) correspond to (i) and (ii) in Figure 3.For the purpose of a better solution, a scanning method is proposed. PointPc is picked from Hc, and it works as the center of a circle whose radius is r=jd¯,0.5≤j≤0.9. This circle scans the blade anticlockwise, and its initial position is at its left horizontal axis. There are rising edges and falling edges on the circle. The points on the rising edges are denoted by P=P1,P2,…,Pm, where m is the quantity of the rising edges on the circle. Meanwhile, the points on the falling edges are denoted by P′=P1′,P2′,…,Pm′. The ordinates of these edge points are Y=y1,y2,…,ym and Y′=y1′,y2′,…,ym′, respectively. The angles between the initial position and these edge points are Φ=ϕ1,ϕ2,…,ϕm and Φ′=ϕ1′,ϕ2′,…,ϕm′, respectively. Let ya=maxY and yb=minY, where min() is a mathematical function to find the minimum values of its input. The propeller rotates from yb to ya as
(5)φ=2πa−bn−ϕa−ϕb,s.t.ϕa−ϕba−b≥2π−ϕa+ϕbm−a+b,2πa−bn−2π−ϕa+ϕb,s.t.ϕa−ϕba−b<2π−ϕa+ϕbm−a+b,where n is the quantity of blades and its value is two in the simulation. When the value of φ is from the first class in (5), the rotational direction is clockwise. When the value of φ is from the second class in (5), the rotational direction is anticlockwise. The angular velocity ω could be calculated out by the points on the rising edges.
(6)ω=f1r=φTdya−yb,where Td is the row delay time of the ERS camera. The points on the falling edges can also be used on the calculation of angular velocity:
(7)ω′=f2r=φ′Tdyc−yd,where yc=maxY′, yd=maxY′, and φ′ is similar to φ. Let ec be the accumulative error corresponding to radius r at point Pc.
(8)ec=∫0.5d0.9d¯rf1r−f2rdr.LetEc=ec∣ec=gPc,Pc∈Hc, where g is a function Pc⟶ec. The value of Ec of different points are normalized and indexed by a color map as shown in Figure 4. Let gPc^=minEc. Point Pc^ is the center of the propeller. In order to promote the accuracy of the calculation, the optimized value of angular velocity ω^ is weighted by r.
(9)ω^=∫0.5d¯0.9d¯rf1r+f2rdr2∫0.5d¯0.9d¯rdr.The result of the simulation is shown in Table1. The actual center position is (500,500), and the most deviation of the center of the propeller is within 5 pixels. The relative errors of angular velocity are no more than 2%.Table 1
Result of the simulation.
Figure
Measured center position (pixel)
Actual velocity (r/s)
Measured velocity (r/s)
(ia)
(500, 500)
0.5
0.510
(ib)
(500, 502)
1
1.006
(ic)
(497, 500)
2
2.003
(iia)
(501, 499)
0.5
0.508
(iib)
(497, 502)
1
0.999
(iic)
(503, 502)
2
2.012


## 3. Experiments and Results

In order to verify the measurement of angular velocity with the ERS camera, two experiments were designed and carried out under controlled conditions. The propeller to be measured had two or multiple blades. All the angular velocities to be detected were higher than the speed of the camera.


### 3.1. Experimental Setup for the Two-Blade Propeller

The measurement system included an ERS camera, an extra light, and a propeller, as shown in Figure5(a). The camera was from Basler in Germany with the model of acA3800-14uc. It had an image sensor of MT9J003 whose shutter type was ERS. The row delay of MT9J003 was 23.09 μs, and its exposure time could be selected from 35 μs to 1,599,535 μs by software. The camera selected resolution was 3840 × 2748, and the maximum frame rate was 14 frames per second. The propeller had two centrosymmetric blades and its diameter was 135 mm. Due to the same color of the propeller and the background, a black poster board was typically used as a background for contrast. The propeller was driven by a motor with encoder and the motor was mounted on the top of a tripod. When the propeller axis was not parallel to the camera axis, the image of the propeller in the video would be declining. Hence, the axis of the camera and the axis of the propeller would be set coaxial. The value of row exposure time was 350 μs and the extra light was adjusted appropriately to get good quality pictures by the camera. A view from the camera was shown in Figure 5(b). Before recording the rotation of the propeller, the ERS camera was calibrated with the method proposed by Heikkila and Silven [19].Figure 5
Rotation experiment on the two-blade propeller: (a) experimental layout and (b) view of the propeller from the RES camera.


### 3.2. Results of the Two-Blade Propeller

The photoelectric encoder on the motor recorded the propeller velocities, while the ERS camera shoot the rotation of the propeller. The propeller rotated anticlockwise at different velocities as illustrated in Figure6. The contours of the propeller were extracted by GMM and they were repaired by a technology that included filling holes and filtering noise. The results of the background segmentation are shown in Figure 7. The areas of color in Figure 7 are the range of the center of the two-blade propeller, when the value of δ is 0.05. The different colors in these areas indicate the value of Ec of different points. The camera captured four images at a time and it took four times to capture 16 images. The actual velocities of those images are shown in the third column of Table 2. The second column of Table 2 shows the deviation of the center of the propeller. The maximum deviation of the center is within six pixels. The angular velocities are obtained by the method proposed in this paper, and they are given in the fourth column of Table 2. The related error of velocity is less than 4%.Figure 6
Two-blade propeller rotated at different velocities: (a) 8.51 r/s, (b) 15.91 r/s, (c) 22.00 r/s, and (d) 25.18 r/s.
(a)
(b)
(c)
(d)Figure 7
Searching the center of the two-blade propeller: (a)–(d) correspond to Figures6(a)–6(d).Table 2
Angular velocity of the two-blade propeller.
Frame number
Center deviation (pixel)
Actual velocity (r/s)
Measured velocity (r/s)
Relative error of velocity (%)
1
(5, 2)
8.51
8.27
2.8
2
(0, 1)
8.53
8.36
2.0
3
(3, 5)
8.57
8.44
1.5
4
(5, −2)
8.59
8.51
0.9
5
(6, −1)
15.88
15.37
3.2
6
(3, 5)
15.91
15.72
1.2
7
(−1, 3)
15.93
15.87
0.3
8
(5, −2)
15.95
16.12
1.1
9
(1, −7)
21.87
21.32
2.5
10
(0, 1)
22.00
21.44
2.6
11
(3, 4)
22.11
22.00
0.1
12
(−3, −2)
22.18
22.10
0.3
13
(−8, −2)
25.07
25.31
1.0
14
(0, 3)
25.11
25.32
0.8
15
(−3, 5)
25.18
25.52
1.4
16
(−5, 2)
25.21
25.73
2.1


### 3.3. Experiment on the Multiblade Propeller

The multiblade propeller was also measured by almost the same equipment in the experiment of the two-blade propeller. The ducted propeller was part of the engineering, as shown in Figure8(a), which was a product of QX MOTOR. The propeller had five blades and its diameter was 64 mm. The row exposure time was adjusted to 350 μs and a view from the camera was shown in Figure 8(b). Due to the lack of an encoder, a handheld tachometer was used to measure the angular velocities of the propeller.Figure 8
Rotation experiment on the multiblade propeller: (a) propeller driven by motor and (b) view of the propeller from the RES camera.


### 3.4. Results of the Multiblade Propeller

The tachometer recorded the multiblade propeller velocities, when the ERS camera shot the rotation of the propeller. Figure9 shows samples of images which were taken during the rotation with different velocities. When the value of δ is set to 0.06, the color areas in Figure 10 are the range of the center of the multiblade propeller. The ERS camera also captured four images at a time, and it captured 16 images within four times. With the method proposed in this paper, the angular velocities were obtained as shown in Table 3. The maximum deviation of the center is within seven pixels, and the related error of velocity is less than 5%. The maximum angular velocity is six times more than the speed of the camera.Figure 9
Multiblade propeller rotated at different velocities: (a) 29.06 r/s, (b) 57.77 r/s, (c) 75.28 r/s, and (d) 109.50 r/s.
(a)
(b)
(c)
(d)Figure 10
Searching the center of the multiblade propeller: (a)–(d) correspond to Figures9(a)–9(d).Table 3
Angular velocity of the multiblade propeller.
Frame number
Center deviation (pixel)
Actual velocity (r/s)
Measured velocity (r/s)
Relative error of velocity (%)
1
(3, 3)
29.06
28.25
2.8
2
(2, 8)
29.18
29.34
0.5
3
(7, −4)
28.86
29.32
1.6
4
(5, −2)
28.90
28.97
0.2
5
(1, −1)
57.77
56.67
1.9
6
(−2, 6)
58.26
57.93
0.6
7
(−5, 7)
58.10
56.54
2.7
8
(−2, 3)
58.14
56.60
2.6
9
(1, 5)
75.28
74.12
1.5
10
(6, 3)
75.45
73.65
2.4
11
(10, 2)
75.55
73.80
2.3
12
(3, 4)
75.89
73.56
3.1
13
(−2, 3)
109.50
107.10
2.2
14
(−1, 0)
109.45
105.48
3.6
15
(−8, 3)
109.12
104.04
4.7
16
(3, −5)
108.93
103.86
4.7


## 3.1. Experimental Setup for the Two-Blade Propeller

The measurement system included an ERS camera, an extra light, and a propeller, as shown in Figure5(a). The camera was from Basler in Germany with the model of acA3800-14uc. It had an image sensor of MT9J003 whose shutter type was ERS. The row delay of MT9J003 was 23.09 μs, and its exposure time could be selected from 35 μs to 1,599,535 μs by software. The camera selected resolution was 3840 × 2748, and the maximum frame rate was 14 frames per second. The propeller had two centrosymmetric blades and its diameter was 135 mm. Due to the same color of the propeller and the background, a black poster board was typically used as a background for contrast. The propeller was driven by a motor with encoder and the motor was mounted on the top of a tripod. When the propeller axis was not parallel to the camera axis, the image of the propeller in the video would be declining. Hence, the axis of the camera and the axis of the propeller would be set coaxial. The value of row exposure time was 350 μs and the extra light was adjusted appropriately to get good quality pictures by the camera. A view from the camera was shown in Figure 5(b). Before recording the rotation of the propeller, the ERS camera was calibrated with the method proposed by Heikkila and Silven [19].Figure 5
Rotation experiment on the two-blade propeller: (a) experimental layout and (b) view of the propeller from the RES camera.


## 3.2. Results of the Two-Blade Propeller

The photoelectric encoder on the motor recorded the propeller velocities, while the ERS camera shoot the rotation of the propeller. The propeller rotated anticlockwise at different velocities as illustrated in Figure6. The contours of the propeller were extracted by GMM and they were repaired by a technology that included filling holes and filtering noise. The results of the background segmentation are shown in Figure 7. The areas of color in Figure 7 are the range of the center of the two-blade propeller, when the value of δ is 0.05. The different colors in these areas indicate the value of Ec of different points. The camera captured four images at a time and it took four times to capture 16 images. The actual velocities of those images are shown in the third column of Table 2. The second column of Table 2 shows the deviation of the center of the propeller. The maximum deviation of the center is within six pixels. The angular velocities are obtained by the method proposed in this paper, and they are given in the fourth column of Table 2. The related error of velocity is less than 4%.Figure 6
Two-blade propeller rotated at different velocities: (a) 8.51 r/s, (b) 15.91 r/s, (c) 22.00 r/s, and (d) 25.18 r/s.
(a)
(b)
(c)
(d)Figure 7
Searching the center of the two-blade propeller: (a)–(d) correspond to Figures6(a)–6(d).Table 2
Angular velocity of the two-blade propeller.
Frame number
Center deviation (pixel)
Actual velocity (r/s)
Measured velocity (r/s)
Relative error of velocity (%)
1
(5, 2)
8.51
8.27
2.8
2
(0, 1)
8.53
8.36
2.0
3
(3, 5)
8.57
8.44
1.5
4
(5, −2)
8.59
8.51
0.9
5
(6, −1)
15.88
15.37
3.2
6
(3, 5)
15.91
15.72
1.2
7
(−1, 3)
15.93
15.87
0.3
8
(5, −2)
15.95
16.12
1.1
9
(1, −7)
21.87
21.32
2.5
10
(0, 1)
22.00
21.44
2.6
11
(3, 4)
22.11
22.00
0.1
12
(−3, −2)
22.18
22.10
0.3
13
(−8, −2)
25.07
25.31
1.0
14
(0, 3)
25.11
25.32
0.8
15
(−3, 5)
25.18
25.52
1.4
16
(−5, 2)
25.21
25.73
2.1


## 3.3. Experiment on the Multiblade Propeller

The multiblade propeller was also measured by almost the same equipment in the experiment of the two-blade propeller. The ducted propeller was part of the engineering, as shown in Figure8(a), which was a product of QX MOTOR. The propeller had five blades and its diameter was 64 mm. The row exposure time was adjusted to 350 μs and a view from the camera was shown in Figure 8(b). Due to the lack of an encoder, a handheld tachometer was used to measure the angular velocities of the propeller.Figure 8
Rotation experiment on the multiblade propeller: (a) propeller driven by motor and (b) view of the propeller from the RES camera.


## 3.4. Results of the Multiblade Propeller

The tachometer recorded the multiblade propeller velocities, when the ERS camera shot the rotation of the propeller. Figure9 shows samples of images which were taken during the rotation with different velocities. When the value of δ is set to 0.06, the color areas in Figure 10 are the range of the center of the multiblade propeller. The ERS camera also captured four images at a time, and it captured 16 images within four times. With the method proposed in this paper, the angular velocities were obtained as shown in Table 3. The maximum deviation of the center is within seven pixels, and the related error of velocity is less than 5%. The maximum angular velocity is six times more than the speed of the camera.Figure 9
Multiblade propeller rotated at different velocities: (a) 29.06 r/s, (b) 57.77 r/s, (c) 75.28 r/s, and (d) 109.50 r/s.
(a)
(b)
(c)
(d)Figure 10
Searching the center of the multiblade propeller: (a)–(d) correspond to Figures9(a)–9(d).Table 3
Angular velocity of the multiblade propeller.
Frame number
Center deviation (pixel)
Actual velocity (r/s)
Measured velocity (r/s)
Relative error of velocity (%)
1
(3, 3)
29.06
28.25
2.8
2
(2, 8)
29.18
29.34
0.5
3
(7, −4)
28.86
29.32
1.6
4
(5, −2)
28.90
28.97
0.2
5
(1, −1)
57.77
56.67
1.9
6
(−2, 6)
58.26
57.93
0.6
7
(−5, 7)
58.10
56.54
2.7
8
(−2, 3)
58.14
56.60
2.6
9
(1, 5)
75.28
74.12
1.5
10
(6, 3)
75.45
73.65
2.4
11
(10, 2)
75.55
73.80
2.3
12
(3, 4)
75.89
73.56
3.1
13
(−2, 3)
109.50
107.10
2.2
14
(−1, 0)
109.45
105.48
3.6
15
(−8, 3)
109.12
104.04
4.7
16
(3, −5)
108.93
103.86
4.7


## 4. Conclusions

The computer vision provides an effective way to measure angular velocities. In this paper, an algorithm has been proposed for the measurement of angular velocity of propeller with the low-speed camera. It bases on the principle of rolling shutter. By analyzing the exposure procedure of the camera, a simulation model was established to simulate the motion of propeller captured by the ERS camera. According to the recorded contour shape of the propeller, a method was proposed to search the rotational center and angular velocity. To validate this method, the angular velocities of a two-blade propeller and a three-blade propeller were measured with this technique. The experimental results show that the proposed method can be used to measure the high-speed rotation of the propeller with the low-speed camera. The maximum angular velocity could be six times more than the speed of the camera. Furthermore, the simulation and experiments show that this method is efficient, applicable, noncontact, and less expensive for the measurement of angular velocity of the propeller.The method proposed in this paper still has some limitations and imperfect aspects. The detection accuracy is sensitive to illumination. Improving exposure time can get good quality images, but the captured image will be very blur when the propeller rotates fast. In the future work, automatic exposure hardware and related algorithm should be added to the system. This method bases on the precise contour of the propeller. So, the background segmentation affects the result greatly. To improve the quality of background segmentation, an artificial neural network will be applied in the future.



---

*Source: 1037083-2018-03-21.xml*