Next Article in Journal
Holey-Cavity-Based Compressive Sensing for Ultrasound Imaging
Next Article in Special Issue
Maglev Train Signal Processing Architecture Based on Nonlinear Discrete Tracking Differentiator
Previous Article in Journal
Time and Memory Efficient Online Piecewise Linear Approximation of Sensor Signals
Previous Article in Special Issue
Degree-of-Freedom Strengthened Cascade Array for DOD-DOA Estimation in MIMO Array Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Angle Measurement of Objects outside the Linear Field of View of a Strapdown Semi-Active Laser Seeker

College of Intelligence Science and Technology, National university of Defense Technology, Changsha 410073, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(6), 1673; https://doi.org/10.3390/s18061673
Submission received: 25 March 2018 / Revised: 4 May 2018 / Accepted: 18 May 2018 / Published: 23 May 2018
(This article belongs to the Special Issue Sensor Signal and Information Processing)

Abstract

:
The accurate angle measurement of objects outside the linear field of view (FOV) is a challenging task for a strapdown semi-active laser seeker and is not yet well resolved. Considering the fact that the strapdown semi-active laser seeker is equipped with GPS and an inertial navigation system (INS) on a missile, in this work, we present an angle measurement method based on the fusion of the seeker’s data and GPS and INS data for a strapdown semi-active laser seeker. When an object is in the nonlinear FOV or outside the FOV, by solving the problems of space consistency and time consistency, the pitch angle and yaw angle of the object can be calculated via the fusion of the last valid angles measured by the seeker and the corresponding GPS and INS data. The numerical simulation results demonstrate the correctness and effectiveness of the proposed method.

1. Introduction

Semi-active laser guidance, which has high precision and is easy to implement, is widely used in precision-guided weapons and equipment [1,2,3,4]. The core device is the semi-active laser seeker [5,6]. It receives the laser spot reflected by an object and detects the precise coordinates of the laser spot center by using a four-quadrant detector. It then calculates the pitch angle and the yaw angle between the object and the seeker’s optical axis. However, a semi-active laser seeker cannot measure the angles when the object is out of the linear field of view (FOV) [7]. This problem is even worse for a strapdown semi-active laser seeker [8]. It is still one of the bottlenecks that restricts the overall application of strapdown semi-active laser seekers.

1.1. Detection Principle of the Strapdown Semi-Active Laser Seeker

The strapdown semi-active laser seeker is composed of a four-quadrant detector, optical system, circuit system and shell. The four-quadrant detector consists of four 90 photosensitive sectors with the same area and the same photoelectric response [9], which are represented by I, I I , I I I and I V , respectively (as shown in Figure 1). Let R be the radius of the photosensitive surface, r be the radius of the reflected laser spot, ( x , y ) be the coordinates of the laser spot center and I i ( i = I , I I , I I I , I V ) be the output current of the i-th photosensitive sectors, which is proportional to the energy of the received laser spot. There are three cases in which ( x , y ) can be detected using a four-quadrant detector. The first meets the condition R < 2 r < 2 R , that is the laser spot is located in the four-quadrant detector and covers each of the four photosensitive sectors (as shown in Figure 1). Then, ( x , y ) can be calculated by [10,11]:
x = I I + I I V I I I + I I I I I I + I I I + I I I I + I I V y = I I + I I I I I I I + I I V I I + I I I + I I I I + I I V
Obviously, ( x , y ) = ( 0 , 0 ) in the case shown in Figure 1a.
The origin of the seeker’s coordinate system is the origin of the four-quadrant detector, the seeker’s X axis is along the optical axis of the optical system, and the seeker’s Y and Z axes are along the Y and X axes of the four-quadrant detector, respectively. Based on ( x , y ) , the pitch angle ε and the yaw angle θ of the object relative to the seeker’s coordinate system can be calculated using:
ε = arctan y f , θ = arctan x f ,
where f is the focal length of the seeker’s optical system. In this case, the satisfied area on the four-quadrant detector is called the linear area of the detector, and the corresponding FOV of the seeker is called the linear FOV. The second case is when the laser spot is within the photosensitive surface of the four-quadrant detector, but it cannot cover all four quadrants of the photosensitive sector, as shown in Figure 2a. In this case, x , y cannot be precisely determined using Equation (1), and it can only be known in which quadrant the laser spot center is located. This area on the four-quadrant detector is called the nonlinear area, and the corresponding FOV of the seeker is called the nonlinear FOV. The third case is when the laser spot is outside the FOV of the semi-active laser seeker, as shown in Figure 2b. In this case, the four-quadrant detector cannot detect any information regarding the laser spot. Therefore, we must ensure that the object is within the linear FOV of the semi-active laser seeker.

1.2. Related Work

To ensure that the reflected laser spot lies in the linear area of the four-quadrant detector, the traditional semi-active laser seeker adopts a platform structure, that is the four-quadrant detector is installed on a complicated and high-precision servo control system [12]. The servo control system, which is composed of inertial measurement components and a dynamic follow-up system, can isolate the attitude movements of the seeker and ensure that the object is always located in the linear FOV of the seeker [4,13]. Although the platform-type semi-active laser seeker is a mature product, it has many disadvantages, such as a complex structure, high cost and large volume.
In recent years, the strapdown semi-active laser seeker has become one of the main development directions of the semi-active seeker [14,15]. It removes the high-precision servo control system and places the four-quadrant detector directly onto the longitudinal axis of the seeker. The advantages are its simpler structure, higher reliability, smaller size, lighter weight and lower cost [16]. The main disadvantage is that due to the detector moving with the seeker, the problem of objects going outside the linear FOV is exacerbated. The current approach is to increase the linear FOV of the seeker via a special optics design [17,18,19,20]. However, there are three shortcomings to this approach. First, the amount by which the linear area of the detector can be expanded and the seeker FOV can be increased via optics design is very limited. Second, when the linear FOV increases, the angle measurement accuracy of the seeker will decrease, which will affect the guidance accuracy of the seeker [8]. Third, increasing the FOV of the seeker results in a shortening of the detection range, and the detection range is very important for the terminal guidance of a missile. Therefore, the problem of angle measurement for a strapdown laser seeker when objects are outside the linear FOV still hinders the full application of the strapdown semi-active laser seeker.
Considering the fact that the strapdown semi-active laser seeker is equipped with GPS and inertial attitude measurement equipment on the missile [21,22], in this work, we make full use of GPS and INS data and propose an angle measurement method for the strapdown semi-active laser seeker through data fusion.

2. Proposed Method

When an object is outside the linear FOV, by solving the space consistency problem and the time consistency problem, the pitch angle and the yaw angle can be calculated by fusing the following data: the current GPS and inertial attitude data, the angles measured by the seeker at the last moment when the object is in the linear FOV and the corresponding GPS and inertial attitude data at that moment. The following gives the specific details of this method.

2.1. Definition of Variables

Assume t 0 to be the last moment at which an object is located within the linear FOV of a strapdown semi-active laser seeker, and let the pitch angle and the yaw angle of the object measured by the seeker at t 0 be ε 0 and θ 0 , respectively. Suppose that at time t 1 , the object is outside the linear FOV of the seeker. Then, the pitch angle ε 1 and the yaw angle θ 1 cannot be measured by the seeker and will be calculated using the proposed method. Our method needs the following data: the object position O T (longitude λ T , latitude L T and height h T ), which is given in advance; the position of the seeker at t 0 (longitude λ 0 , latitude L 0 and height h 0 ); the attitude data of the seeker with the yaw-pitch-roll rotation order at time t 0 (yaw angle φ 0 , pitch angle ψ 0 and roll angle γ 0 ) or the quaternions at t 0 ( q 0 _ 0 , q 1 _ 0 , q 2 _ 0 , q 3 _ 0 ); the position of the seeker at t 1 (longitude λ 1 , latitude L 1 and height h 1 ); and theattitude data of the seeker with the yaw-pitch-roll rotation order at time t 1 (yaw angle φ 1 , pitch angle ψ 1 and roll angle γ 1 ) or the quaternions at t 1 ( q 0 _ 1 , q 1 _ 1 , q 2 _ 1 , q 3 _ 1 ). The above positions and attitude data can be obtained via GPS and INS [23,24]. The problem of time consistency can be solved by precisely aligning the above data with the corresponding time.

2.2. Definitions of the Coordinate Systems

To solve the space consistency problem, we define the following coordinate systems.
The Earth-centered frame O e - X e Y e Z e : The origin O e is the center of the Earth. The axis O e Z e is perpendicular to the Earth’s equatorial plane and points toward the North Pole. The axis O e X e lies in the Earth’s equatorial plane and points to the Greenwich meridian. The axis O e Y e is perpendicular to the plane O e X e Z e and forms a right-hand coordinate system with O e X e and O e Z e .
The local navigation frame O n - X n Y n Z n : The local navigation frame is defined as having a north-up-east order. The origin O n is the centroid of the seeker. The axis O n Y n is collinear with the normal of the navigation frame’s reference ellipsoid at the penetration point. The axis O n X n lies in the Meridian plane, is perpendicular to O n Y n and points toward the north. The axis O n Z n is determined according to the right-hand rule.
The body frame O b - X b Y b Z b : The origin O b is the centroid of the seeker. The axis O b X b coincides with the longitudinal axis of the seeker and points toward the forward direction. The axis O b Y b lies in the longitudinal symmetry plane of the seeker, is perpendicular to O b X b and points upward. The axis O b Z b is perpendicular to the O b X b Y b plane and forms a right-hand coordinate system with O b X b and O b Y b .
The on body line-of-sight frame O s - X s Y b Z s : The origin O s is the centroid of the seeker. O s X s points toward the object along the line of sight. The O s Y s axis, which points upward, is on a plane that contains O s X s and is perpendicular to O s X s and the plane O b X b Z b at the same time. O s Z s is determined by the right-hand rule.

2.3. Analysis and Computation of the Proposed Method

In practice, the seeker moves with the missile all the time; thus, the position and attitude of the seeker at t 1 are different from those at t 0 . As shown in Figure 3, let O b 0 be the seeker position at t 0 and O b 0 - X b 0 Y b 0 Z b 0 be the body frame at O b 0 . Moreover, let O b 1 be the seeker position at t 1 and O b 1 - X b 1 Y b 1 Z b 1 be the body frame at O b 1 . From t 0 to t 1 , the body frame has both attitude movements and position movements simultaneously. Without loss of generality, we assume that the attitude movements occur first; the body frame O b 0 - X b 0 Y b 0 Z b 0 transforms into the intermediate body frame O b 0 - X b 1 Y b 1 Z b 1 ; and the, a translational movement occurs, with O b 0 - X b 1 Y b 1 Z b 1 being translated to O b 1 - X b 1 Y b 1 Z b 1 . Therefore, the analysis and calculation required to solve the space consistency problem can be conducted in two stages.

2.3.1. First Stage of the Proposed Method

In the first stage, we consider only the variation of the pitch angle and the variation of the yaw angle caused by the seeker’s attitude motion from t 0 to t 1 . In this stage, the pitch angle ε 1 and yaw angle θ 1 of the object should be calculated in the intermediate body frame O b 0 - X b 1 Y b 1 Z b 1 , with the calculation involving the following coordinate transformations: O b 0 - X b 1 Y b 1 Z b 1 is transformed into the local navigation frame O n - X n Y n Z n ; then O n - X n Y n Z n is transformed into the body frame O b 0 - X b 0 Y b 0 Z b 0 ; and finally, O b 0 - X b 0 Y b 0 Z b 0 is transformed into the on body line-of-sight frame O s - X s Y s Z s . The specific steps are as follows.
Step 1: O b 0 - X b 1 Y b 1 Z b 1 is transformed into O n - X n Y n Z n by the roll-pitch-yaw rotation order with the rotations of γ 1 , φ 1 and ψ 1 , respectively [25]. The transform matrix C b 1 n 0 is calculated according to Equation (3).
C b 1 n 0 = cos ψ 1 cos φ 1 cos γ 1 sin ψ 1 cos φ 1 + sin γ 1 sin φ 1 sin γ 1 sin ψ 1 cos φ 1 + cos γ 1 sin φ 1 sin ψ 1 cos γ 1 cos φ 1 sin γ 1 cos φ 1 cos ψ 1 sin φ 1 cos γ 1 sin ψ 1 sin φ 1 + sin γ 1 cos φ 1 sin γ 1 sin ψ 1 sin φ 1 + cos γ 1 cos φ 1
To avoid the singularity problem of the Euler angles at about 90 , C b 1 n 0 can be calculated based on the quaternions q 0 _ 1 , q 1 _ 1 , q 2 _ 1 , q 3 _ 1 according to Equation (4).
C b 1 n 0 = q 0 _ 1 q 0 _ 1 + q 1 _ 1 q 1 _ 1 q 2 _ 1 q 2 _ 1 q 3 _ 1 q 3 _ 1 2 q 1 _ 1 q 2 _ 1 + q 0 _ 1 q 3 _ 1 2 q 1 _ 1 q 3 _ 1 q 0 _ 1 q 2 _ 1 2 q 1 _ 1 q 2 _ 1 q 0 _ 1 q 3 _ 1 q 0 _ 1 q 0 _ 1 q 1 _ 1 q 1 _ 1 + q 2 _ 1 q 2 _ 1 q 3 _ 1 q 3 _ 1 2 q 2 _ 1 q 3 _ 1 + q 0 _ 1 q 1 _ 1 2 q 1 _ 1 q 3 _ 1 + q 0 _ 1 q 2 _ 1 2 q 2 _ 1 q 3 _ 1 q 0 _ 1 q 1 _ 1 q 0 _ 1 q 0 _ 1 q 1 _ 1 q 1 _ 1 q 2 _ 1 q 2 _ 1 + q 3 _ 1 q 3 _ 1
Step 2: O n - X n Y n Z n is transformed into O b 0 - X b 0 Y b 0 Z b 0 by the yaw-pitch-roll rotation order with rotations of ψ 0 , φ 0 and γ 0 , respectively. The transform matrix C n 0 b 0 is calculated according to Equation (5).
C n 0 b 0 = cos ψ 0 cos φ 0 sin ψ 0 cos ψ 0 sin φ 0 cos γ 0 sin ψ 0 cos φ 0 + sin γ 0 sin φ 0 cos γ 0 cos φ 0 cos γ 0 sin ψ 0 sin φ 0 + sin γ 0 cos φ 0 sin γ 0 sin ψ 0 cos φ 0 + cos γ 0 sin φ 0 sin γ 0 cos φ 0 sin γ 0 sin ψ 0 sin φ 0 + cos γ 0 cos φ 0
Similar to Step 1, C n 0 b 0 can be calculated based on the quaternions q 0 _ 0 , q 1 _ 0 , q 2 _ 0 , q 3 _ 0 according to Equation (6):
C b 0 n 0 = q 0 _ 0 q 0 _ 0 + q 1 _ 0 q 1 _ 0 q 2 _ 0 q 2 _ 0 q 3 _ 0 q 3 _ 0 2 q 1 _ 0 q 2 _ 0 q 0 _ 0 q 3 _ 0 2 q 1 _ 0 q 3 _ 0 + q 0 _ 0 q 2 _ 0 2 q 1 _ 0 q 2 _ 0 + q 0 _ 0 q 3 _ 0 q 0 _ 0 q 0 _ 0 q 1 _ 0 q 1 _ 0 + q 2 _ 0 q 2 _ 0 q 3 _ 0 q 3 _ 0 2 q 2 _ 0 q 3 _ 0 q 0 _ 0 q 1 _ 0 2 q 1 _ 0 q 3 _ 0 q 0 _ 0 q 2 _ 0 2 q 2 _ 0 q 3 _ 0 + q 0 _ 0 q 1 _ 0 q 0 _ 0 q 0 _ 0 q 1 _ 0 q 1 _ 0 q 2 _ 0 q 2 _ 0 + q 3 _ 0 q 3 _ 0
Step 3: O b 0 - X b 0 Y b 0 Z b 0 is transformed into O s - X s Y s Z s by the pitch-yaw rotation order with rotations of ε 0 and θ 0 . The transform matrix C b 0 s is:
C b 0 s = cos θ 0 cos ε 0 cos θ 0 sin ε 0 sin θ 0 sin ε 0 cos ε 0 0 sin θ 0 cos ε 0 sin θ 0 sin ε 0 cos θ 0
To summarize, we can obtain the transformation matrix from O b 0 - X b 1 Y b 1 Z b 1 to O s - X s Y s Z s a:s
C b 0 s 1 = C b 0 s · C n 0 b 0 · C b 1 n 0
Step 4: According to Equation (8), we can calculate the pitch angle ε 1 using [25]:
ε 1 = arctan C b 0 S 1 [ 1 , 2 ] C b 0 S 1 [ 1 , 1 ] , if C b 0 S 1 [ 1 , 2 ] 0 and C b 0 S 1 [ 1 , 1 ] > C b 0 S 1 [ 1 , 2 ] arctan C b 0 S 1 [ 1 , 1 ] C b 0 S 1 [ 1 , 2 ] + π 2 , if C b 0 S 1 [ 1 , 1 ] > 0 and C b 0 S 1 [ 1 , 1 ] C b 0 S 1 [ 1 , 2 ] arctan C b 0 S 1 [ 1 , 1 ] C b 0 S 1 [ 1 , 2 ] π 2 , if C b 0 S 1 [ 1 , 1 ] 0 and C b 0 S 1 [ 1 , 1 ] > C b 0 S 1 [ 1 , 2 ] arctan C b 0 S 1 [ 1 , 2 ] C b 0 S 1 [ 1 , 1 ] π , if C b 0 S 1 [ 1 , 2 ] < 0 and C b 0 S 1 [ 1 , 1 ] C b 0 S 1 [ 1 , 2 ] arctan C b 0 S 1 [ 1 , 2 ] C b 0 S 1 [ 1 , 1 ] + π , if C b 0 S 1 [ 1 , 2 ] 0 and C b 0 S 1 [ 1 , 1 ] > C b 0 S 1 [ 1 , 2 ] arctan C b 0 S 1 [ 1 , 1 ] C b 0 S 1 [ 1 , 2 ] + π 2 , if C b 0 S 1 [ 1 , 1 ] < 0 and C b 0 S 1 [ 1 , 1 ] C b 0 S 1 [ 1 , 2 ] arctan C b 0 S 1 [ 1 , 1 ] C b 0 S 1 [ 1 , 2 ] + π 2 , if C b 0 S 1 [ 1 , 1 ] 0 and C b 0 S 1 [ 1 , 1 ] < C b 0 S 1 [ 1 , 2 ] arctan C b 0 S 1 [ 1 , 2 ] C b 0 S 1 [ 1 , 1 ] , if C b 0 S 1 [ 1 , 2 ] > 0 and C b 0 S 1 [ 1 , 1 ] C b 0 S 1 [ 1 , 2 ] ,
and the yaw angle θ 1 using:
θ 1 = arc sin C b 0 s 1 [ 1 , 3 ]

2.3.2. Second Stage of the Proposed Method

In the second stage, we analyze the variation of the pitch angle and the variation of the yaw angle caused by translating O b 0 - X b 1 Y b 1 Z b 1 to O b 1 - X b 1 Y b 1 Z b 1 . As shown in Figure 3, to calculate the vector O b 1 O T , we first need to calculate the vector O b 0 O b 1 and the vector O b 0 O T in frame O b 0 - X b 1 Y b 1 Z b 1 . The steps of calculating O b 0 O b 1 in O b 0 - X b 1 Y b 1 Z b 1 are as follows.
Step 1: Calculate the radius R W 0 of curvature in the prime vertical of the Earth and the radius R N 0 of curvature in the meridian of the Earth at O b 0 by:
R W 0 = a e 2 a e 2 cos 2 L 0 + b e 2 sin 2 L 0 1 / 2 R N 0 = R w 0 b e 2 a e 2 ,
and calculate the radius R W 1 of curvature in the prime vertical of the Earth and the radius R N 1 of curvature in the meridian of the Earth at O b 1 by:
R W 1 = a e 2 a e 2 cos 2 L 1 + b e 2 sin 2 L 1 1 / 2 R N 1 = R w 1 b e 2 a e 2 ,
where a e = 6,378,140 m and b e = a e ( 1 1 1 298.257 ) 298.257 ) are the lengths of the Earth’s long and short half-axles, respectively.
Step 2: Calculate the coordinates e x 0 , e y 0 , e z 0 of O b 0 in the Earth-centered frame O e - X e Y e Z e using:
e x 0 = ( R w 0 + h 0 ) cos L 0 cos λ 0 e y 0 = ( R w 0 + h 0 ) cos L 0 sin λ 0 e z 0 = ( R N 0 + h 0 ) sin L 0 ,
and the coordinates e x 1 , e y 1 , e z 1 of O b 1 in O e - X e Y e Z e using:
e x 1 = ( R W 1 + h 1 ) cos L 1 cos λ 1 e y 1 = ( R W 1 + h 1 ) cos L 1 sin λ 1 e z 1 = ( R W 1 + h 1 ) sin L 1
Step 3: Calculate the transformation matrix C e n 0 from the Earth-centered frame O e - X e Y e Z e to the local navigation frame O n - X n Y n Z n at O b 0 by:
C e n 0 = sin L 0 cos λ 0 sin L 0 sin λ 0 cos L 0 cos L 0 cos λ 0 cos L 0 sin λ 0 sin L 0 sin λ 0 cos λ 0 0
Step 4: Calculate the transformation matrix C n 0 b 1 from the local navigation frame at O b 0 to O b 0 - X b 1 Y b 1 Z b 1 by:
C n 0 b 1 = C b 1 n 0
From Step 1 to Step 4, O b 0 O b 1 can be obtained by:
O b 0 O b 1 = C n 0 b 1 · C e n 0 · e x 1 e x 0 e y 1 e y 0 e z 1 e z 0
The vector O b 0 O T in frame O b 0 - X b 0 Y b 0 Z b 0 is calculated as follows.
Step 1: Calculate the radius R W T of curvature in the prime vertical of the Earth and the radius R N T of curvature in the meridian of the Earth at the object position O T by:
R W T = a e 2 a e 2 cos 2 L T + b e 2 sin 2 L T 1 / 2 R N T = R w T b e 2 a e 2
Step 2: Calculate the coordinates e x T , e y T , e z T of O T in the Earth-centered frame by:
e x T = R w T + h T cos L T cos λ T e y T = R w T + h T cos L T sin λ T e z T = R N T + h T sin L T
Step 3: Calculate the distance D i s t 0 between the seeker and the object at t 0 by:
D i s t 0 = e x 0 e x T 2 + e y 0 e y T 2 + e z 0 e z T 2
Step 4: Since the 3-2 order pitch angle and yaw angle of the object in O b 0 - X b 1 Y b 1 Z b 1 are ε 1 and θ 1 , respectively, O b 0 O T is obtained by:
O b 0 O T = D i s _ t 0 · cos θ 1 · cos ε 1 D i s _ t 0 · cos θ 1 · sin ε 1 D i s _ t 0 · sin θ 1
Step 5: The vector O b 1 O T is calculated using:
O b 1 O T = O b 0 O T O b 0 O b 1
Let the three components of O b 1 O T be Δ x _ t 1 , Δ y _ t 1 and Δ z _ t 1 ; then, the pitch angle ε 1 and the yaw angle θ 1 of the object at time t 1 are finally determined as [25]:
ε 1 = arctan Δ y _ t 1 Δ x _ t 1 , if Δ y _ t 1 0 and Δ x _ t 1 > Δ y _ t 1 arctan Δ x _ t 1 Δ y _ t 1 π 2 , if Δ x _ t 1 > 0 and Δ x _ t 1 Δ y _ t 1 arctan Δ x _ t 1 Δ y _ t 1 π 2 , if Δ x _ t 1 0 and Δ x _ t 1 > Δ y _ t 1 arctan Δ y _ t 1 Δ x _ t 1 π , if Δ y _ t 1 < 0 and Δ x _ t 1 Δ y _ t 1 arctan Δ y _ t 1 Δ x _ t 1 + π , if Δ y _ t 1 0 and Δ x _ t 1 > Δ y _ t 1 arctan Δ x _ t 1 Δ y _ t 1 + π 2 , if Δ x _ t 1 < 0 and Δ x _ t 1 Δ y _ t 1 arctan Δ x _ t 1 Δ y _ t 1 + π 2 , if Δ x _ t 1 0 and Δ x _ t 1 < Δ y _ t 1 arctan Δ y _ t 1 Δ x _ t 1 , if Δ y _ t 1 > 0 and Δ x _ t 1 Δ y _ t 1
and:
θ 1 = arc tan Δ z _ t 1 Δ x _ t 1 2 + Δ y _ t 1 2

3. Numerical Simulation Results

3.1. Numerical Simulation Setups

We carry out two simulation experiments using MATLAB to verify and evaluate the proposed method. The setups of the simulations are as follows: the measurement period of the seeker is 50 ms; the linear FOV of the seeker is ± 10 ; and the FOV of the seeker is ± 20 . In the simulations, the seeker performs a sinusoidal-like motion to place the object at different positions within the seeker’s FOV. As shown in Table 1, the object is in the linear FOV in the first frame, and the pitch angle and yaw angle are 6.752 and 7.187 , respectively. The object exits the linear FOV in the second frame (0.05 s) and stays in the nonlinear FOV from 0.05 s to 0.90 s. Then, the object leaves the FOV of the seeker at 0.95 s, re-enters the nonlinear FOV at 2.80 s and stays in the nonlinear FOV until 3.65 s. Next, it enters the linear FOV again and stays in the linear FOV until 5.05 s. Subsequently, the object again exits the linear FOV at 5.10 s and stays in the nonlinear FOV until 6.25 s. The object then stays outside the FOV between 6.30 s and 7.40 s before entering the nonlinear FOV at 7.45 s. In the simulations, we calculate the pitch angle and the yaw angle via the proposed method when the object exits the linear FOV of the seeker and evaluate the method’s performance by comparing it to the ground truth.

3.2. Numerical Simulation Results

The purpose of the first numerical simulation is to verify the correctness of the proposed method. In the numerical simulation, the GPS and INS data do not contain errors. The results are shown in Figure 4. Figure 4a shows the pitch angle results, while Figure 4b shows the yaw angle results. The blue `∘’ represents the ground truth, the red `+’ the result of the proposed method and the green `*’ the error between the proposed algorithm and the ground truth. It can be seen that both the pitch angle errors and the yaw angle errors are very close to zero throughout the simulation. To be more precise, the absolute values of both the pitch angle errors and the yaw angle errors are less then 1 × 10 6 , which are caused by the numerical truncation of the simulation software. Therefore, this simulation result proves the correctness of the proposed method.
The purpose of the second numerical simulation is to evaluate the angle measurement accuracy of the proposed method when the object is outside the linear FOV. In the numerical simulation, both the GPS data and INS data contain errors. The error in terms of the GPS latitude, longitude and height is 10 m. To make this numerical simulation more challenging, we assume that the INS uses a low-precision MEMS gyroscope [26,27] and that the angle drift ratio is 20 / h. In addition, we assume that the INS has been working for 60 s after the initial alignment. Thus, the initial attitude error of the INS is 0.333 . The attitude error of the INS during the numerical simulation is represented by the green curve in Figure 5c. The simulation results are shown in Figure 5: Figure 5a shows the pitch angle results; Figure 5b shows the yaw angle results; and Figure 5c shows the error between the proposed method and the ground truth. In Figure 5a,b, the blue curve represents the ground truth, while the red curve represents the result of the proposed method. In Figure 5c, the red `+’ represents the pitch angle error, and the blue `∘’ represents the yaw angle error. It can be seen that as time progresses, the object enters the nonlinear FOV or comes out of the FOV, and the angular measurement error of this method increases with the increase of the attitude error of the INS. Specifically, within the 3.6 s when the object leaves the linear FOV for the first time, as the attitude error of the INS increases to 0.354 , the absolute value of the pitch angle error of the proposed method increases from 0 to 0.14 , and the absolute value of the yaw angle error increases from 0 to 0.068 . Furthermore, we can conclude that under the above GPS and INS error conditions, this method can ensure that the angular measurement error is less than 0.2 in the 6.5 s when the object is outside the linear FOV.
From the theoretical derivation and simulation process, it can be seen that when the object is outside the linear FOV, the angle measurement accuracy of the method increases as the GPS accuracy and INS accuracy increase and as time reduces. In practice, the accuracy of the INS is higher than 20 / h, and the amount of time the object spends outside the linear FOV does not exceed 3.6 s. Thus, the proposed method can achieve better angle measurement performance than exhibited during the numerical simulation.

4. Conclusions

To solve the problem in which a strapdown semi-active laser seeker cannot measure the angles of objects outside the linear FOV, we make full use of GPS and INS data and propose an angle measurement method based on information fusion. When an object is within the nonlinear FOV or outside the FOV, the pitch angle and the yaw angle of the object can be calculated via a fusion of the last valid angles measured by the seeker and the corresponding GPS and INS data. The numerical simulation results show that the proposed method can tolerate a certain amount of GPS and INS errors and ensure the angular measurement error is less than 0.2 in the 6.5 s when the object is outside the linear FOV. In general, the proposed method is simple, accurate and effective for angle measurement of objects outside the linear FOV of a strapdown semi-active laser seeker.

Author Contributions

Y.Z. proposed the method. Y.Z. and Z.Z. conceived of and designed the experiments. H.C. performed the experiments and helped perform the data analysis. Y.Z. wrote the paper. All the authors reviewed and approved the entire manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (No. 61403412).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
FOVlinear field of view
INSinertial navigation system
MEMSmicro electromechanical systems
References

References

  1. Robert, A.N.; Harold, L. Copperhead Semiactive Laser Guidance System Development. J. Guid. Control Dynam. 1979, 2, 374–381. [Google Scholar] [CrossRef]
  2. Hubbard, K.; Katulka, G.; Lyon, D.; Petrick, D.; Fresconi, F.; Horwath, T. Low-Cost Semi-Active Laser Seekers for US Army Application. In Proceedings of the International Telemetering Conference, San Diego, CA, USA, 27–30 October 2008. [Google Scholar]
  3. Allen, R.; Volpi, J.; Roemerman, S. COTS lens and detector characterization for low cost, miniature SAL seekers. In Proceedings of the Aiaa Guidance, Navigation, and Control Conference, Boston, MA, USA, 19–22 August 2013. [Google Scholar]
  4. Titterton, D.H. Military Laser Technology and Systems; Artech House: Norwood, MA, USA, 2015. [Google Scholar]
  5. English, J.E.; Pittman, W.C. Polarized Semi-Active Laser Last Pulse Logic Seeker Using a Staring Focal Plane Array. U.S. Patent 6,987,256, 17 January 2006. [Google Scholar]
  6. Taylor, B.; Schaub, M.; Jenkins, D. Projectile Guidance System Including a Compact Semi-Active Laser Seeker. U.S. Patent 8,207,481, 26 June 2012. [Google Scholar]
  7. Zhang, L.; Yang, Y.; Xia, W.; Zhu, X.; Chen, W.; Lu, Y. Linearity of quadrant avalanche photodiode in laser tracking system. Chin. Opt. Lett. 2009, 7, 728–731. [Google Scholar] [CrossRef]
  8. Zhang, X.; Du, Z.Y.; Qiao, Y.F. Study on linear field of strapdown semi-active laser seeker. Chin. Opti. 2015, 8, 415–421. [Google Scholar] [CrossRef]
  9. Zhang, H.; Chen, Y.S.; Geng, T.W.; Tao, C. Study on Main Factors Affecting Position Detection Accuracy of Four-Quadrant Detector. Chin. J. Lasers 2015, 42, 1217002. [Google Scholar] [CrossRef]
  10. Tang, Y.; Gu, G.; Qian, W.; Chen, Q.; Jun, Z. Laser spot center location algorithm of four-quadrant detector based on Gaussian distribution. Infrared Laser Eng. 2017, 46, 206003. [Google Scholar] [CrossRef]
  11. Cheng, X.; Yang, Y.; Hao, Q. Analysis of the Effects of Thermal Environment on Optical Systems for Navigation Guidance and Control in Supersonic Aircraft Based on Empirical Equations. Sensors 2016, 16, 1717. [Google Scholar] [CrossRef] [PubMed]
  12. Waldmann, J. Line-of-sight rate estimation and linearizing control of an imaging seeker in a tactical missile guided by proportional navigation. IEEE Trans. Control Syst. Technol. 2002, 10, 556–567. [Google Scholar] [CrossRef]
  13. Gapiński, D.; Krzysztofik, I.; Koruba, Z. Analysis of the dynamics and control of the modified optical target seeker used in anti-aircraft rocket missiles. J. Theor. Applied Mech. 2014, 52, 629–639. [Google Scholar]
  14. Lee, C.H.; Hyun, C.; Lee, J.G.; Choi, J.Y.; Sung, S. A hybrid guidance law for a strapdown seeker to maintain lock-on conditions against high speed targets. J. Electr. Eng. Technol. 2013, 8, 190–196. [Google Scholar] [CrossRef]
  15. Fan, J.; An, X.; Lin, P.; Dong, S. A modified biased pursuit guidance for missile with strapdown seeker. In Proceedings of the 2015 Chinese Automation Congress, Wuhan, China, 27–29 November 2015; pp. 1815–1818. [Google Scholar] [CrossRef]
  16. Wei, C.Z.; Han, Y.P.; Cui, N.G.; Xu, H.C. Fifth-Degree Cubature Kalman Filter Estimation of Seeker Line-of-Sight Rate Using Augmented-Dimensional Model. J. Guid. Control Dynam. 2017, 40, 1–8. [Google Scholar] [CrossRef]
  17. Hartman, R.L.; Griffin, R.W. Faceted Ball Lens for Semi-Active Laser Seeker. U.S. Patent 7,185,845, 6 March 2007. [Google Scholar]
  18. Zhang, J.Q.; Xie, F.; Xue, Q.S.; Xin, J.X. Laser guided lens based on four-quadrant detector. Chin. Opt. 2015, 8, 471–479. [Google Scholar] [CrossRef]
  19. Chen, M.; Yang, Y.; Jia, X.; Gao, H. Investigation of positioning algorithm and method for increasing the linear measurement range for four-quadrant detector. Optik Int. J. Light Electr. Opt. 2013, 124, 6806–6809. [Google Scholar] [CrossRef]
  20. Wei, Z.; He, Z.; Zhang, X. Optimization design of small optical detection system with large instantaneous field of view. Infrared Laser Eng. 2016, 45, 518002. [Google Scholar]
  21. Vergez, P.L.; Mcclendon, J.R. Optimal control and estimation for strapdown seeker guidance of tactical missiles. J. Guid. Control Dynam. 1982, 5, 225–226. [Google Scholar] [CrossRef]
  22. Hyslop, G.; Gerth, D.; Kraemer, J. GPS/INS integration on the standoff land attack missile (SLAM). In Proceedings of the IEEE Symposium on Position Location and Navigation. A Decade of Excellence in the Navigation Sciences, Las Vegas, NV, USA, 20–20 March 1990; pp. 407–412. [Google Scholar]
  23. Chu, H.J.; Tsai, G.J.; Chiang, K.W.; Duong, T.T. GPS/MEMS INS data fusion and map matching in urban areas. Sensors 2013, 13, 11280–11288. [Google Scholar] [CrossRef] [PubMed]
  24. Wu, Y.X.; Wang, J.L.; Hu, D.W. A New Technique for INS/GNSS Attitude and Parameter Estimation Using Online Optimization. IEEE Trans. Signal Process. 2014, 62, 2642–2655. [Google Scholar] [CrossRef]
  25. Groves, P.D. Principles of GNSS, Inertial, and Multisensor Integrated Navigation Systems; Artech House: Norwood, MA, USA, 2013. [Google Scholar]
  26. Piyabongkarn, D.; Rajamani, R.; Greminger, M. The development of a MEMS gyroscope for absolute angle measurement. IEEE Trans. Control Syst. Technol. 2005, 13, 185–195. [Google Scholar] [CrossRef]
  27. Du, S.; Sun, W.; Yang, G. MEMS IMU Error Mitigation Using Rotation Modulation Technique. Sensors 2016, 16. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Detection principle of a four-quadrant detector. (a) The spot lies in the center; (b) The spot lies within the linear area.
Figure 1. Detection principle of a four-quadrant detector. (a) The spot lies in the center; (b) The spot lies within the linear area.
Sensors 18 01673 g001
Figure 2. Examples of laser spot distribution on the four-quadrant detector. (a) The spot lies in the nonlinear area; (b) The spot is located outside the FOV.
Figure 2. Examples of laser spot distribution on the four-quadrant detector. (a) The spot lies in the nonlinear area; (b) The spot is located outside the FOV.
Sensors 18 01673 g002
Figure 3. Equivalent decomposition of the seeker’s movement from t 0 to t 1 .
Figure 3. Equivalent decomposition of the seeker’s movement from t 0 to t 1 .
Sensors 18 01673 g003
Figure 4. Results of the first numerical simulation. (a) The pitch angle results; (b) The yaw angle results.
Figure 4. Results of the first numerical simulation. (a) The pitch angle results; (b) The yaw angle results.
Sensors 18 01673 g004
Figure 5. Results of the second numerical simulation. (a) The pitch angle results; (b) The yaw angle results; (c) The error between the proposed method and the ground truth.
Figure 5. Results of the second numerical simulation. (a) The pitch angle results; (b) The yaw angle results; (c) The error between the proposed method and the ground truth.
Sensors 18 01673 g005
Table 1. Relationship between the object position and the seeker FOV.
Table 1. Relationship between the object position and the seeker FOV.
Time t (s)00.05–0.900.95–2.752.80–3.653.7–5.055.10–6.256.30–7.407.45–8.5
Position relative to FOVLinearNonlinearOutsideNonlinearLinearNonlinearOutsideNonlinear

Share and Cite

MDPI and ACS Style

Zheng, Y.; Chen, H.; Zhou, Z. Angle Measurement of Objects outside the Linear Field of View of a Strapdown Semi-Active Laser Seeker. Sensors 2018, 18, 1673. https://doi.org/10.3390/s18061673

AMA Style

Zheng Y, Chen H, Zhou Z. Angle Measurement of Objects outside the Linear Field of View of a Strapdown Semi-Active Laser Seeker. Sensors. 2018; 18(6):1673. https://doi.org/10.3390/s18061673

Chicago/Turabian Style

Zheng, Yongbin, Huimin Chen, and Zongtan Zhou. 2018. "Angle Measurement of Objects outside the Linear Field of View of a Strapdown Semi-Active Laser Seeker" Sensors 18, no. 6: 1673. https://doi.org/10.3390/s18061673

APA Style

Zheng, Y., Chen, H., & Zhou, Z. (2018). Angle Measurement of Objects outside the Linear Field of View of a Strapdown Semi-Active Laser Seeker. Sensors, 18(6), 1673. https://doi.org/10.3390/s18061673

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop