Performance Assessment of Pre-processing Filters for Infrared Search and Track Applications

To enhance detection probability and to reduce false alarms, infrared imagery is pre-processed before subjecting it to detection algorithms in infrared search and track systems. Pre-processing algorithms are used to predict the complex background and then to subtract the predicted background from the original image. The difference image is passed to the detection algorithm to further distinguish between the target and the background and/ or noise more accurately. A number of pre-processing algorithms have been reported in literature, with their relative advantages and disadvantages. This paper brings out the computational complexities and simulation results of various algorithms for assessing their relative performances. Based on these parameters, statistical algorithms in general and max-min algorithms in particular, are recommended to be used for infrared search and track systems.

The detection of small moving targets in clutter is an important subject in the area of signal/image processing. The algorithms for the detection of small and point targets are fundamental part of infra-red search and track (IRST) system and play vital role in the success of such systems. Typically, the spatial pre-processing step is performed on the input image to predict the background, and consequently enhance the signal-to-clutter ratio (SCR). The detection algorithm may result in many false targets and this requires using the post-processing algorithms to reduce the false alarms. Pre-processing algorithms are used to predict the complex background and then to subtract the predicted background from the original image. The difference image is passed to the detection algorithm to further distinguish between the target and the background and/or noise more accurately. The aim is to fit the background as closely as possible in the original image without diminishing the target signal. Many pre-processing algorithms have been reported in the literature1. It has experimentally been verified that the detection of dim point size targets in cluttered background is not possible without increasing the signal-to-clutter plus noise ratio (SCNR) by pre-processing of IR data. The output of pre-processing algorithm is passed to the detection algorithm which actually detects the targets and it may result in many false targets and this may require using the post-processing algorithms to reduce false alarms and to generate candidate target list. The block diagram of the image processing algorithms for detection of point and/or small targets is given in Fig. 1.

The pre-processing algorithms can be broadly classified as morphological and statistical filters1.

2.1 Morphological Filters

Morphological filters are based on successive use of erosion and dilation operations using a structuring element. Mathematical morphology provides an approach to the processing of images which is based on shape. The morphological operations tend to simplify image data, preserving their essential shape characteristics and eliminating irrelevances. Opening process (erosion of image by structuring element, followed by dilation of the result by the same structuring element) eliminates bridges (narrow edge like set of pixels) connecting two regions of the image. It rounds outward-pointing corners in the image while leaving the inward-pointing corners unaffected. The closing process (dilation of image by structuring element, followed by erosion of the result by the same structuring element) rounds the inward-pointing corners in the image while leaving the outward corners unaffected. The problem with morphological filters is that the result is highly dependent on size and shape of the structuring element. The adaptive selection of size and shape of structuring element is a challenging task.

The gray-scale morphological operations can be defined as follows: Let f(x) and g(x) be 1-D gray-scale functions of coordinate x, where f(x) is the image and g(x) is the structuring element. Let E represents Euclidean space, and, F, GE2. Then, f: FE and g : G → E , and the basic gray-scale dilation, erosion, closing and opening operations are defined in Eqns (1), (2), (3), and (4), respectively.

${D}_{g}\left(f\right)\text{\hspace{0.17em}}=\text{\hspace{0.17em}}\mathrm{max}\left[f\left(x\text{\hspace{0.17em}}-z\right)\text{\hspace{0.17em}}+\text{\hspace{0.17em}}g\left(z\right)\right],\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall z\text{\hspace{0.17em}}\in \text{\hspace{0.17em}}G,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\text{\hspace{0.17em}}-\text{\hspace{0.17em}}z\text{\hspace{0.17em}}\in \text{\hspace{0.17em}}F$ (1)

${E}_{g}\left(f\right)\text{\hspace{0.17em}}=\text{\hspace{0.17em}}\mathrm{min}\left[f\left(x\text{\hspace{0.17em}}+\text{\hspace{0.17em}}z\right)\text{\hspace{0.17em}}-\text{\hspace{0.17em}}g\left(z\right)\right],\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall z\in \text{\hspace{0.17em}}G,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\text{\hspace{0.17em}}+\text{\hspace{0.17em}}z\text{\hspace{0.17em}}\in \text{\hspace{0.17em}}F$ (2)

${O}_{g}\left(f\right)\text{\hspace{0.17em}}=\text{\hspace{0.17em}}\left(f\text{\hspace{0.17em}}\circ \text{\hspace{0.17em}}g\right)\text{\hspace{0.17em}}·\text{\hspace{0.17em}}g\left(x\right)$ (3)

${C}_{g}\left(f\right)=\left(f•g\right)\circ g\left(x\right)$ (4)

The computational load for some of the common composite morphological filters is shown in Table 1. The composite morphological pre-processing filters (in reference to Sr. Nos. 1, 2, 3 and 4 of the Table 1) used for simulation are defined by the Eqns (5), (6), (7), and (8), respectively. The Selective-Morpho CO_OC is computed similar to Morpho CO_OC filter with the difference that only alternate elements of the structuring elements are considered for processing.

$Morpho\text{\hspace{0.17em}}C\text{\hspace{0.17em}}_\text{\hspace{0.17em}}O\text{\hspace{0.17em}}=\text{\hspace{0.17em}}Mean\text{\hspace{0.17em}}\left({C}_{g}\left(f\right),\text{\hspace{0.17em}}{O}_{g}\left(f\right)\right)$ (5)

$Morpho\text{\hspace{0.17em}}CO\text{\hspace{0.17em}}_\text{\hspace{0.17em}}OC\text{\hspace{0.17em}}=\text{\hspace{0.17em}}Mean\left({C}_{g}\left({O}_{g}\left(f\right)\right),\text{\hspace{0.17em}}{O}_{g}\left({C}_{g}\left(f\right)\right)\right)$ (6)

(7)

$Morpho\text{\hspace{0.17em}}ECOC\text{\hspace{0.17em}}_\text{\hspace{0.17em}}DOCO\text{\hspace{0.17em}}=\text{\hspace{0.17em}}Mean\text{\hspace{0.17em}}\left(\begin{array}{l}{E}_{g}\left[{C}_{g}\left({O}_{g}\left({C}_{g}\left(f\right)\right)\right)\right],\\ {D}_{g}\left[{O}_{g}\left({C}_{g}\left({O}_{g}\left(f\right)\right)\right)\right]\end{array}\right)$ (8)

2.2 Statistical Filters

Statistical filters are based on the fact that the statistical behaviour of the target and the background are not identical to predict the background. Filters like mean, median, maxmedian, max-min, mean median, selective median, 2-D Gaussian etc.2-7. The methodology here is to replace the background pixel by mean/median, etc. of the neighborhood. This assumes that one odd pixel occupied by the target (with higher intensity than the background) will not affect the statistics drastically. The computational loads of some of the common statistical pre-processing filters are shown in Table 2.

This study aims to find out a filter that is consistent and gives good performance, which can be used in the detection process. Towards this end, clouds were simulated using standard software program. Targets with known intensity and statistics were embedded in these video sequences.

The challenge is to decide a parameter for comparative analysis. The generally accepted parameter for this is the change in signal-to-noise ratio (SCNR). However, as brought out later in the paper, no clear-cut direction could be found "based on this parameter. Next, we tried contrast sensitivity analysis. The ratio of "the output contrast sensitivity to the input contrast sensitivity provided a good direction for selection of the pre-processing filter.

SCNR8 is defined as given in the following equation:

$SCNR\left(dB\right)\text{\hspace{0.17em}}=\text{\hspace{0.17em}}10\mathrm{log}\left(\frac{\left({S}_{\mathrm{max}}\text{\hspace{0.17em}}-\text{\hspace{0.17em}}\text{μ}}{{\text{σ}}^{2}}\right)$ (9)

where, Smax is the signal peak value, ì is the mean, and σ is the standard deviation of the background. The contrast sensitivity, also known as Weber’s fraction, is proposed to be the parameter for deciding the efficacy of pre-processing filters. Contrast sensitivity9 is defined as given in Eqn (10).

${C}_{s}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}\frac{\Delta I}{I}$ (10)

where, ΔΙ is the difference value between the object and the background (mean value), I is the mean background intensity.

To compare various pre-processing methods, a number of synthetic clouds were generated by means of infrared scene simulator7,10. The infrared scene simulator (IRSS) uses the modified Gardener’s Fourier series method, selfsimilar method, and Perlin’s noise method to generate the synthetic infra-red clouds. The IRSS typically does the object creation and its modelling, background scene creation and it’s modeling, and the integration of created IR object into simulated background scene and image rendering process. The snapshot of the IRSS with generated cloud is shown in Fig. 2.

The real infra-red video data containing the long-range small/point air targets is not available for study and analysis. Therefore the small/point targets were synthetically embedded in these clouds through a target embedding routine in MATLAB11 environment with known target and background profile, as shown in Fig. 3.

Finally, these images were subjected to various filtering processes. Appropriate pre-processing filters were selected for simulation process. The parameters for comparative study (SCNR, differential intensity and input-output contrast values) were automatically generated and logged in a file for analysis. The run time snapshot of the pre-processing simulation program is shown in Fig. 4.

A total number of 50 types of background cloud sequences (each sequence with 25 image frames) were generated and targets with known target characteristics (target position, shape and gray-level difference between the target and the neighbouring background) were embedded in each of the image frame.

The clouds generated as above, were embedded with a point target with input SCNRs of 5 dB, 10 dB and 20dB (wrt the surrounding window of 5x5 pixels. The pre-processing filters discussed1 were simulated on infra-red image frames as above. SCNR of the input and output images was calculated and the improvement in SCNR was found. The SCNR of the filtered output image was calculated in the same window of 5x5 pixels. The average performance of these filters over the full sample of clouds is presented in Fig. 5.

As seen in Fig. 5, no clear cut pattern is emerging with respect to the performance of any of the filters. Hence, there is a need to find an alternate parameter for comparison of algorithms. Contrast sensitivity, as defined above, was chosen as a parameter for comparing various filters, specifically, the ratio between the output contrast sensitivity and input contrast sensitivity was plotted.

The same target-background scenarios that were used for SCNR calculations were used for calculating contrast sensitivity. Encouraging results were seen. The five test cases were analysed on each of 50 image sequences with gray-level difference (GLD) of 2, 5, 10, and 20 between the target and the background on 8-bit gray-scale images. The results are similar for other test cases also. It is seen that max-min filter is almost always outscoring all other filters in terms of contrast sensitivity. It was also observed that there is a constant improvement in output contrast by various filters. Whenever, the input contrast is high, the output contrast saturates at 8-bit (255) value.

The max-min filter is defined as given in the following equation:

$y\left(m,\text{\hspace{0.17em}}n\right)\text{\hspace{0.17em}}=\text{\hspace{0.17em}}\mathrm{max}\left[{z}_{1},\text{\hspace{0.17em}}{z}_{2},\text{\hspace{0.17em}}{z}_{3},\text{\hspace{0.17em}}{z}_{4}\right]$ (11)

where y(m, n) is the output of the (2N+1) th order max- min filter and z1, z2, z3 and z4 are minimum values of the middle row vector, middle column vector and the two diagonal vectors of the convolution kernel centered around the pixel of interest, respectively. The larger the value of N, the more was the spike suppression. This property was used to detect small targets. The maximum operation was chosen with criteria to preserve edges and discontinuities in the filtered signals. If other ranks are chosen rather than the maximum, it can be easily shown that the filter will tend to blur edges of the original signal. In many applications, all the zi, s estimates may not be needed. The minimum operations on middle row vector, middle column vector, and the two diagonal vectors initially predict the background as the minimum intensity value over the horizontal, vertical or diagonal direction. This operation preserves the low-intensity targets, and hence helps to enhance the SCNR in the absolute difference image.

The graphs showing the performance of the simulated filters for GLD of 2, 5, 10, and 20 are shown in Figs 6, 7, 8, and 9, respectively.

Pre-processing filters are very useful to suppress the background and enhance the contrast sensitivity of the difference image frame for detection of small targets. From the simulation, it is seen that the statistical filters have superior performance in the present application. It was also found that max-min filter outperforms the other statistical filters. Moreover, the computational requirements are also O(n), making it an attractive choice for real-time application. Selective-median, improved hybrid morphological filter, median filter and max-median filter also perform satisfactorily if the input contrast of the target is good. However, their performance degrades for very low contrast targets.

1. Srivastava, H.B.; Kumar, V.; Verma, H.K. & Sundaram, S.S. Image pre-processing algorithms for detection of small/point airborne targets. Def. Sci. J., 2009, 59(2), 166-74.

2. Boyle, R. & Thomas, R. Computer vision: A first course. Blackwell Scientific Publications, 1988. pp. 42-44.

3. Gonzalez, R. & Woods, R. Digital image processing, Addison-Wesley Publishing Company, 1992. 191 p.

4. Helstrom, C. Probability and stochastic process for engineers. Macmillan, New York, 1991.

5. Boyle, R. & Thomas, R. Computer vision: A first course. Blackwell Scientific Publications. 1988. pp. 32-34.

6. Deshpande, Suyog D.; Er, M.H.; Ronda, V. & Phillip, Chan. Max-mean and max-median filters for detection of small-targets. In Proceedings of SPIE: Signal and Data Processing of Small Targets, 3809, 1999. pp. 74-83.

7. Desai, Uday B. & Merchant, S.N. Development of signal processing/image processing techniques and IR scene simulator for infra-red search and track (IRST) systems- Phase II. Final Completion Report, SPANNN Lab, Department of Electrical Engineering, IIT Bombay, Mumbai. October 2006. pp. 23-24.

8. Diani, M.; Baldacci, A. & Corcini, G. Novel background removal algorithm for navy infrared search and track systems. Optical Engineering, 2001, 40(8), 1729-734.

9. Pratt, William K. Digital image processing. Ed. 3. Wiley- Interscience Publication, 2001.

10. Desai, Uday B. & Merchant, S.N. Infrared search and tracking. SPANNN Lab, Department of Electrical Engineering, IIT Bombay, Mumbai. Report-2, June 2002.

11. Image processing toolbox: For use with MATLAB-User’s guide. Ver. 5, October 2005.

 Mr Hari Babu Srivastava obtained his BE (Electrical Engg) from University of Roorkee and MTech (Electrical Engg) from Indian Institute of Technology, Kanpur, in 1983 and 1989, respectively. Currently as a Scientist ‘G’ he is heading the Fire Control System Design Centre, IRDE, Dehradun. He is life member of Instrument Society of India and Optical Society of India and member of Computer Society of India and IEEE, USA. He has presented 20 research papers in various national/international conferences and published six research papers in refereed international/national journals. Prof Vinod Kumar did his MTech and PhD in Electrical Engg. from University of Roorkee in 1975 and 1984, respectively. He is Head, Department of Electrical Engineering and Continuing Education Centre and the Coordinator of the Quality Improvement Centre at IIT, Roorkee. His research interests include ECG signal processing and analysis, digital signal processing and transducer instrumentation. He is the Life Fellow of the Institution of Engineers (I), Institution of Electronics and Telecommunication Engineers and Systems Society of India, Life Member of Indian Society for Continuing Engineering Education, Ultrasonic Society of India, Indian Society for Technical Education, Instrument Society of India, Biomedical Engineering Society of India, IEE UK and Senior Member of IEEE, USA. Prof H.K. Verma is Head of the Department of Electrical Engineering, Dean of Research and Industry Liaison, Sponsored Research and Industrial Consultancy and Faculty Affairs at IIT, Roorkee. His research interests include measurement and instrumentation, digital protective relays, microprocessor and computer applications in instrumentation, power system protection, monitoring and control, medical instrumentation, telemedicine and medical informatics, telemetry and remote / supervisory control, automated test equipment and systems. He is a member of National Working Group on Philosophy and Application of Backup Protection for Transmission Lines, Reactors and Busbars and CIGRE National Study Committee on Information System and Telecommunication for Power Systems. Mr S.S. Sundaram is presently the Director of IRDE, Dehradun. He is engaged in the development of optical instrumentation, night vision instrumentation and laser instrumentation for the Services in addition to realisation of fire control systems for tanks and ships. He has vast experience in the field of communication electronic warfare and radar electronic warfare and has fielded integrated systems of high complexity, operating on wide area network, as leader of major programmes/projects. As Chief Designer, he has led technology teams in the development of communication electronic countermeasure systems including the design of fast settling frequency synthesizers, surveillance and monitor receivers and jammer exciters at DLRL, Hyderabad. He is a fellow of IETE, Optical Society of India, and Instrument Society of India.