# Ballistic Target PHD Filter Based on Infrared Focal Plane Ambiguous Observation

Authors: Wei Zhao; Shucai Huang
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

Space-based early warning system, the main detection means of which is passive detection based on focal plane, is an important part of ballistic missile defense system. The focal plane is mainly composed of CCD, and its size can reach the micron level, so the pixel is often regarded as point of no area in image postprocessing. The design of traditional tracking methods is based on this, and the observation based on the focal plane is modeled as the azimuth with random noise. However, this modeling is inaccurate. In the context of space-based detection, CCD cannot be simplified as a point, and its size should be considered. And the corresponding observation cannot be treated as azimuth with random noise. In this paper, the observation based on focal plane is modeled as Unambiguously Generated Ambiguous (UGA) measurement. The PHD filter algorithm is redesigned and simplified. The simulation results show that the algorithm based on UGA measurement observation model has better tracking effect compared with that based on traditional observation model. This method provides technical support for more accurate target tracking for space-based early warning system.


---

## Body

## 1. Introduction

The world has always attached great importance to the development of space-based early warning system to deal with the threat of ballistic missile. The United States is the most advanced country that has carried out a lot of research work about MIDAS, DSP, SBIRS, and STSS [1, 2]. DSP has been applied in Gulf War, which proved the strategic significance of space-based early warning system. In recent years, the detection mechanism extends from infrared to visible and ultraviolet, and the detection band also extends from a single band to multiple bands and to multispectral and hyperspectral [3, 4]. Limited by technology, the size of CCD does not meet the actual detection needs [5–7], especially the need for spatial resolution, which will affect the accuracy of target tracking. There are many research works about detection and tracking capability of space-based early warning system and specific target detection and tracking algorithms [8, 9]. After investigating a variety of detection and tracking algorithms based on space-based early warning system, it is found that the infrared passive detection is mostly used in which the direction of arrival (DOA) is detected. The specific data are the azimuth and pitch angle of the target relative to the sensor, and the detection data will also be analyzed in Section 2. Actually, the infrared detection is mostly based on the focal plane [10]; thus the data is recorded in the form of infrared image which describes the detected infrared radiation intensity of the target and background. The pixel is the smallest resolution treated as a point pointing to a precise direction in usual treatment. But this is not realistic and we could not obtain a precise direction from a pixel. In ordinary application the small pointing error may not have a significant impact, but it has a completely different effect in space-based detection. It is for this reason that accurate modeling of space-based infrared detection is proposed, and tracking filter based on this observation model is designed.The observation error is analyzed in Section2, and the observation based on focal plane imaging was modeled. In Section 3, the motion of target image on focal plane is studied, and the constant acceleration model is proved to be reasonable. In Section 4, the PHD filter algorithm based on UGA measurement observation modeling is given. In Section 5, the difficulty of implementing the filter algorithm is analyzed, and a simplified algorithm is given. In Section 6, the PHD filters based on traditional observation modeling and UGA measurement observation modeling are simulated and compared in the same scenario. And the conclusions are given in the last section.


## 2. Traditional Observation and UGA Measurement Observation

Passive observation refers to an observation method that only obtains the target’s direction of arrival rather than its position. It is an important detection method in military application because of its good concealment and high angular resolution. Commonly used passive detection methods include infrared, ultraviolet, visible, electromagnetic wave, and sonar methods which determine the direction of the target relying on receiving the signal emitted or reflected by the target. Passive observation is often modeled as angle containing noise, just as it is shown in (1); this is a generic modeling method that applies to all of the above-mentioned detection methods.(1)zk+1=ηk+1xk+1+Vk+1.In traditional studies on target tracking based on space-based infrared detection, the observation data is concerned as angle containing noise too, which could be calledpoint observation. In fact this is not true, because zk+1 could not be obtained accurately.The orbit of satellite is so high that the detector can get a wide field of view. This is both its advantage and its disadvantage. Due to technical conditions, the infrared CCD cannot be small enough to meet the actual detection needs; what is more, the detector is so high, just as it is shown in Figure1, where the detection area of a CCD can be tens of thousands of square meters. And this is the great difficulty of target positioning and tracking.Figure 1
CCD and its footprint.Therefore the observation based on focal plane is different from traditional passive observation; it belongs to Unambiguously Generated Ambiguous (UGA) measurement. The observer could know the CCD in whose field of view the target appears but could not get the accurate direction. UGA measurement observation has the following characteristics: First, the observation modeling itself has a certain degree of ambiguity; second, the relationship between the observation and the target state can be described by an accurate sensor conversion model; for example,z=η(x).In fact, the observation data obtained by focal plane imaging can be expressed accurately by pixels corresponding to CCDs whose signals exceed the threshold. Such observations can be calledpixel observation, which could be denoted as Ez. The image of the target in the focal plane is in Ez but could not be obtained accurately, limited to resolution. So this kind of observation is UGA measurement observation. The observation equation can be expressed as follows:(2)EZk+1=Eηk+1xk+1+Vk+1.It is assumed that random noise Vk+1 denotes observation noise caused by detector position error and pointing error. Vk+1 is white Gaussian noise, and its covariance is only related to the detector itself; E(Vk+1·Vk+1T)=Rk+1.Generalized observation is defined asΘzk+1≜Ezk+1-Vk+1, and the generalized likelihood function is shown as(3)fΘzk+1∣xk+1≜Pr⁡ηxk+1∈Θzk+1.


## 3. Target Motion Modeling

To establish a reasonable target motion model, the dynamic model of ballistic target is established firstly. The acceleration, velocity, and trajectory in plane are shown in Figure2 and the sketch map of focal plane imaging is shown in Figure 3. It is assumed that the target is launched from (E 0°, N 0°), the direction is northeast 45°, and the detector is on the GEO satellite above that launch point. After calculation, we can get the image trajectory of the target on focal plane, just as shown in Figure 4.Figure 2
Acceleration velocity and trajectory of a target.Figure 3
Sketch map of focal plane imaging.Figure 4
Trajectory of target image on focal plane.The distance and height of target changing over time are shown in Figure5. As shown in Figure 6, through the quadratic fitting of the distance versus time and height versus time, it is found that the fitting effect is good and the residual is small enough. And the quadratic polynomial corresponds to the constant acceleration (CA) model; hence the CA model is chosen as target motion model in this paper, just as it is shown in(4)xx˙x¨k+1=1ΔΔ2201Δ001xx˙x¨k+Wk+1,where Δ denotes time step, Wk+1 denotes process noise of target motion equation, and E(Wk+1·Wk+1T)=Qk+1.Figure 5
Range and height versus time.Figure 6
Range and height fitting and residual in 10 seconds.


## 4. Filtering Algorithm

In actual combat environment, multitarget tracking is unavoidable for space-based infrared early warning system. Therefore, PHD filter is chosen as the tracking algorithm in this paper. Firstly, PHD filter is introduced briefly.The PHD filter is an advanced approximation technique for multitarget Bayesian filtering and consists of two parts: prediction and update.It is known that probability hypothesis density isDk∣k(x)=abbr.Dk∣k(x∣Z(k)).Prediction Equation. Dk+1∣k(x)=abbr.Dk+1∣k(x∣Z(k))(5)Dk+1∣kx=bk+1∣kx+∫Fk+1∣kx∣x′·Dk∣kx′dx′,where pseudo-Markov transition density is(6)Fk+1∣kx∣x′≜psx′·fk+1∣kx∣x′+bk+1∣kx∣x′.Update Equation. Dk+1∣k+1(x)=LZ(x)·Dk+1∣k(x)(7)LZx≜1-pDx+pDx∑z∈ZLzxλcz+Dk+1∣kpDLz.Under the condition of UGA measurement observation modeling, likelihood function in the update equation should be changed to (8). And the filter schemes for different observation models are shown in Figure 7.Figure 7
The filter schemes for different observation models.In theory,Pr⁡(η(xk+1)∈Θ) in (3) is obtained through integrating η(x) over Θ.(8)Pr⁡ηxk+1∈Θzk+1=∫Θzk+1Φzk+1dz.Φ(z) follows normal distribution in observation space, the mean is η(x), and the covariance is Vk+1.


## 5. Simplification of Generalized Likelihood Function

As described in previous section, in update step, the integration of normal distribution over generalized observation space needs to be calculated. But the integral complexity is not conductive to the implementation of the filter algorithm. Proper simplification is necessary.A physical explanation of the PHD filter is given in literature [11], and another expression of update equation is given as shown in(9)fZk+1∣c∈Ξ,Zk=pDc·fZk+1∣Zk·1μ∑j=1mLzjcfzj∣Zk+1-pDc·fZk+1∣Zk,where zj denotes the jth pixel observation and  c denotes state cell.The box part in (9) corresponds to the summation part of generalized likelihood function in (7). It can be seen that this part describes the ratio of the single-target likelihood function to the spatial density of the observation process, that is, the ratio of probability of zj origins from state cell c to probability of zj origins from others state cells and false alarms. Therefore, it is appropriate to choose a method to replace the ratio of probability here.In this paper, the particle swarm is chosen to replace the pixel observation to describe its relationship with the prediction observation. As shown in Figure8, nine particles are sampled in the pixel, and the sampling criterion is that the particles are evenly distributed in the pixel. The number of particles is denoted by v, and the particle is denoted as pij(i=1⋯v), where i denotes the ith particle. The likelihood function of particle and prediction observation can be expressed as(10)Lpijηxk+1∣k=12πV1/2exp-12pij-ηxk+1∣kτV-1pij-ηxk+1∣k.Finally, the pseudo-likelihood function between the pixel observation and the prediction observation can be obtained by summation.(11)LΦzk+1ηxk+1∣k=∑i=1vLpijηxk+1∣k.Figure 8
Sampling particles in image element.


## 6. Simulation and Comparison

To verify the proposed PHD filter based on pixel observation or UGA measurement observation, the following simulation scenario is designed. The PHD filter based on UGA measurement observation is compared with traditional image target tracking filter. The OSPA is chosen as the index to evaluate the performance of multitarget tracking, and the GM-PHD filter is used.Scenario. It is assumed that the earth geosynchronous orbit satellite is equipped with infrared detector. The subastral point is located at (E 0°, N 0°), and the focal plane of the detector is area array of 1024 × 1024 infrared CCD. The size of every CCD is 30 μm × 30 μm, and the frame rate is 1 Hz. The LOS random error follows Gaussian distribution, the mean is zero, and derivation is 20 μrad. From time t0, three ballistic targets appear in the field of view and are recorded by focal plane. The launch parameters of these three targets are shown in Table 1.Table 1
Launch point, direction, and time of three targets.
Target 1
Target 2
Target 3
Launch point
(W 5°, N 5°)
(E 5°, N 5°)
(E 5°, S 5°)
Direction
135°
225°
315°
Discovery time
t
0
t
0
+
10
s
t
0
+
20
sAs shown in Figure9, three targets are denoted by T1, T2, and T3, respectively, and their starting points are indicated by “O.” The trajectories in Figure 9 are projection of the three ballistic targets on the focal plane. The movements of the three targets are affected by the thrust of the engine, the gravity of the earth, and the aerodynamics. So their movement on the focal plane is a complex movement. And the transition from 3-dimensional space to focal plane is shown in Figure 3. Through the study of the target image trajectory, taking into account the computational complexity, the motion of the target image can be modeled as the constant acceleration model with Gaussian noise.Figure 9
Target image trajectories on focal plane.
(a)
(b)
(c)
(d)In fact, the sensor is often affected by clutter and noise, resulting in some false alarms. Therefore, the false alarms are considered in the simulation scenario. The number of false alarms obeys Poisson distribution with a parameter of 5, and they are evenly distributed throughout the field of view.Figure10 shows the false alarms on the focal plane in this scenario. The logos are enlarged for easy display.Figure 10
False alarms on focal plane.The tracking effects of PHD filters based on two kinds of observation models are shown in Figure11. It is not easy to see difference between these two figures, so OSPA of these two multitarget tracking filters are given in Figure 12, and the statistical results of OSPA and runtime of two filters are shown in Table 2.Table 2
OSPA and runtime.
Point observation
Pixel observation
OSPA
30.5
22.5
Runtime
1.048 s
3.282 sFigure 11
Estimated and true trajectories of three targets: (a) is target 1, (b) is target 2, and (c) is target 3.
(a)
(b)
(c)Figure 12
OSPA based on point observation and pixel observation.Running Environment. The running environment is as follows: Intel® Core™ i3-2120 CPU @ 3.30 GHz 3.29 GHz, 3.24 GB RAM.Running Platform. The running platform is MATLAB R2014a.It is easy to see that the PHD filter based on pixel observation can obtain a more accurate target trajectory. And this is of great importance for ballistic target tracking based on space-based infrared focal plane detection and a breakthrough in traditional deterministic observation.Obviously, the running time is closely related to the number of particles sampled from the pixel observation, and the specific statistical results are shown in Figure13.Figure 13
Running time and OSPA with different number of particles.It is easy to see that the running time increases linearly with the increase in the number of particles. OSPA continues to decline as the number of particles increases; however, when the number of particles increases to 5, OSPA’s decline has become very slow. Therefore, the number of particles used to replace the pixel observation needs to be chosen reasonably for different situations to meet certain accuracy.


## 7. Conclusions

The UGA measurement observation model based on focal plane imaging is studied. A more accurate observation model was established, and PHD filter algorithm based on pixel observation is given. Simulation results verify the accuracy of the observation modeling, and more accurate tracking is achieved. Unfortunately, the computation spent in the implementation of PHD filter based on pixel observation is several times more than that based on point observation. And the computation is closely related to the number of particles used in the simplification of the generalized likelihood function. To obtain better computation efficiency, we will continue to study the simplification of UGA measurement observation integral method.



---

*Source: 5702592-2017-11-06.xml*