Recursive Bessel filter for stray light corrections in an imager

1. A system for correcting infrared (IR) radiances comprising:

an imager for receiving IR radiances and outputting intensity data corresponding to pixels in the imager, and

a stray light estimator for receiving the intensity data and estimating stray light in the intensity data,

wherein the stray light estimator includes a Bessel filter for low pass filtering the estimated stray light and a module for providing corrected intensity data as an output to a user, and

the stray light estimator includes an interpolator for receiving multiple bandwidths of the intensity data and interpolating the multiple bandwidths to obtain the estimated stray light in a first bandwidth, and

the Bessel filter is configured to receive the estimated stray light and provide filtered intensity data in the first bandwidth, as an output to the user.

2. The system of claim 1 wherein

the Bessel filter is comprised of a ratio between (a) an (n−1) order polynomial of a complex number and (b) an n order polynomial of a complex number,

where n is a positive integer.

3. The system of claim 1 wherein

the Bessel filter is an n-pole recursive digital filter, where n is a positive integer.

4. The system of claim 3 wherein

n is the number 8.

5. The system of claim 3 wherein

the Bessel filter includes an infinite response curve, and

the response curve monotonically decreases in amplitude as a function of frequency per pixels.

6. The system of claim 1 wherein

the Bessel filter includes a subtraction module for subtracting the estimated stray light from the intensity data and outputting the corrected intensity data in the first bandwidth.

7. The system of claim 1 wherein

the multiple bandwidths include 3.9 um and other bandwidths, and

the first bandwidth is the 3.9 um bandwidth.

8. The system of claim 1 wherein

the imager includes an imager disposed in a satellite,

the IR radiances include radiances reflected from the Earth, and

the Bessel filter is configured to low pass filter a portion of the IR radiances.

9. A filtering system comprising:

a module for outputting pixel data, and

a Bessel filter for low pass filtering the pixel data,

wherein the Bessel filter has an impulse response, H(z), where H(z) is a ratio of polynomials, and z is a variable.

10. The system of claim 9 wherein

the Bessel filter is an n-pole recursive digital filter, where n is a positive integer.

11. The system of claim 10 wherein

n is the number 8.

12. The system of claim 9 wherein

the Bessel filter is comprised of a ratio between (a) an (n−1) order polynomial of a complex number and (b) an n order polynomial of a complex number,

where n is a positive integer.

13. The system of claim 9 wherein

the Bessel filter includes an infinite response curve, and

the response curve monotonically decreases in amplitude as a function of frequency per pixels.

14. The system of claim 9 wherein

the pixel data includes signals in the infrared (IR) region.

15. A method of filtering image data in the infrared region (IR) comprising the steps of:

receiving IR radiances and outputting intensity data;

estimating noisy stray light in the intensity data;

filtering the estimated noisy stray light to obtain smooth stray light;

subtracting the smooth stray light from the received IR radiances;

outputting corrected IR radiances to a user, after the subtracting step.

16. The method of claim 15 wherein

filtering includes using a Bessel filter to low pass filter the estimated noisy stray light.

17. The method of claim 16 wherein

using the Bessel filter includes forming an n-pole recursive digital filter, wherein n is an integer number.

18. The method of claim 16 wherein

using the Bessel filter includes forming a ratio between (a) an (n−1) order polynomial of a complex number and (b) an n order polynomial of a complex number, where n is an integer number.

19. The method of claim 15 wherein

receiving the IR radiances includes receiving multiple radiances having different bandwidths;

estimating the noisy stray light includes interpolating the received multiple radiances to form the noisy stray light in a single bandwidth; and

filtering the estimated noisy stray light includes using a Bessel filter to low pass filter the noisy stray light in the single bandwidth to form the smooth stray light.