Correlated modulation imaging (CMI)

Abstract

In the present invention of Correlated Modulation Imaging (CMI), the weak optical image signal (and therefore the signal current) is modulated, and the signal integration direction is correlated to the modulation. Therefore, the dark and/or background current, which are not modulated, are cancelled, while the signal current is integrated. As a result, the total integration time of the signal of each pixel is increased, and its signal to noise ratio and dynamic range are improved. Besides, the CMI noise spectrum peaks at the modulation frequency, and therefore, the detector's 1/f and other low frequency noises can be suppressed. In the present invention, the method and theory of CMI, as well as the means and steps for the realization of CMI, are explicitly developed. Two versions of CMOS devices (CMI unit pre-amplifier version 1 and 3), with their circuitry design and testing data are presented as the critical component for correlated modulation imaging. A prototype chip of the CMI circuitries has been fabricated by the MOSIS 0.5-micron foundry. SPICE circuitry simulation has demonstrated the functioning of CMI unit pre-amplifier version 1 and 3 as designed. Experimental work has demonstrated the functioning of CMI unit pre-amplifier version 1 as designed, as well as the key features of CMI—increase of integration time, and improvement of signal to noise ratio and dynamic range.

Description

BACKGROUND OF THE INVENTION

1. Field of the Invention

This invention relates to a new method of electronic imaging and its applications. The invention particularly relates to a new imaging method of focal plane arrays (FPA), in which the image signal is modulated while the background and/or dark current are not modulated, and the image signal integration is correlated to its modulation.

2. Description of Related Art

By using intersubband transition, the quantum well infrared photodetectors (QWIP) operating in photoconductor (PC) mode are made of wide-bandgap materials, such as GaAs/AlGaAs. The maturity of GaAs material growth and its processing technology places the QWIP in a strategically important position. Large-sized (512×512 or greater), multi-color, reliable, reproducible, and low-cost QWIP FPAs working in the important 8-14μ long wavelength infrared (LWIR) atmospheric window seem to be within reach for military, scientific, and commercial applications. However, the effort to develop a competitive LWIR QWIP FPA is still being challenged by QWIP's one impediment: its large dark current. In order to get significant IR response, the PC mode QWIP must be biased at a few volts, which leads to a large dark current due to electron thermal emission, thus making it difficult for its use at T>60 K when the wavelength λ>8μ. Although the 8-12μ QWIP's dark current has been reduced and responsivity increased through innovative designs, it still cannot operate at the convenient liquid nitrogen temperature of 77 K to perform the critically important 8-12 micron thermal imaging.

Note that the condition of

I d or I b >>I s   (1)

exists in other cases, too. In a recent NASA research project, the 15μ QWIP's dark current (working at 55 K) I d is three orders greater than the scene or signal current I s . In solar magnetography, the background photocurrent I b is more than four orders greater than I s . In long wavelength IR spectroscopy and spectroscopic imaging, the narrow bandwidth of signal mandates a low ratio of I s /I d . All these three cases, as well as other scientific and biomedical imaging and spectroscopic applications involving weak signals, have the same feature as the QWIP FPA. The conventional staring array imaging method cannot perform satisfactorily, since its signal integration time is limited by the predominant dark or background current. As the background current I b plays an equivalent role as the dark current I d , only I d is used in the following discussion. Since imaging and spectroscopy share the same signal integration and readout technology, we use the terminology of FPA imaging, which can be applied to both cases.

To deal with the difficulties of imaging a weak signal with presence of a strong background or dark current, the concept of current memory background subtraction (CMBS) was proposed by the Jet Propulsion Laboratory (JPL) Group of Advanced Imager and Focal Plane Technology. The photodetector's large dark current is subtracted by a current sink, which was memorized during the FPA's calibration phase. The result is an increase of effective dynamic range by one order, albeit at the cost of increased shot noise.

Note that in the case of a single photodetector recovering weak optical signal from strong noise, the method of lock-in amplification was invented. In lock-in amplification, a weak signal is modulated together with a strong noise, which mainly consists of the white shot noise and the 1/f flick noise. After going through an amplifier with a bandwidth centered at the modulation frequency, and a band pass filter to maintain the modulated signal while eliminating unwanted noises at frequencies different from the modulation frequency, the modulated signal is fed to a phase sensitive detector (PSD). As a unit square wave multiplier in phase with the modulator, the PSD feeds a DC component, which is proportional to the modulated signal amplitude, to a low pass filter. When the noise is pure white noise, the resulting improvement of the signal to noise ratio due to lock-in amplification is $$ ( S / N ) o ( S / N ) i = B LPF B WN ( 2 )

where (S/N) i is the input signal to noise ratio, (S/N) o the output signal to noise ratio, B LPF the bandwidth of the low pass filter, and B WN the bandwidth of the original white noise.

It is interesting to note that although the lock-in amplification and the CMI presented in this invention both utilize the modulation of signal, the two methods are different. The lock-in amplifier works only for a single detector, while the correlated modulation imaging technique does for an array of detectors (one-dimensional linear array or two-dimensional -area array). The single photodetector with lock-in amplification operates continuously, with a phase sensitive detector and a low pass filter as the key components. In the method of CMI, the improvement of signal to noise ratio is achieved due to the modulation of signal and correlated multi-cycle integration of the signal, with the increase of total integration time as the key factor.

BRIEF DESCRIPTION OF THE SYMBOLS AND DRAWINGS

Definition of Symbols

I s DC image signal photocurrent generated in the detector by the scene of imaging. I s is defined as the average of the real image signal current i s , which is time varying, through its integration cycle.

I d DC dark current thermally generated. I d is defined as the average of the time-varying real dark current i d .

I b DC background current, which is generated by photons of the same wavelengths under detection. As I s and I d , I b is defined as the average of the time varying real background current i b . When I b can be modulated separately from signal current I s , the background current I b plays the same role as the dark current I d . In this invention, we treat the two as equivalent, using I d to denote I b +I d .

T i Integration time. For CMI, T i =mτ, where m is the number of integration cycles of the CMI unit pre-amplifier, and τ the period of each cycle.

R Signal to noise ratio, defined as $$ R = Signal / Noise = N s n 2 _ ,

 where N s is the number of electrons due to the signal current I s , and {overscore (n 2 )} is the root mean square value of the total number of electrons due to random noise.

D Dynamic range in decibels, defined as D=20 log 10 R max , when the maximum integration time is utilized.

f Frequency.

a) Angular frequency. ω=2πf

ω m Modulation frequency. ω m =2πf m =2π/τ.

T(ω) Noise transmission function of a conventional single cycle integrator for a conventional FPA

T(ω)={square root over (V o (ω,φ)V o *(ω,φ))}, where V o (ω,φ) is the output voltage of the integrator with a unit harmonic current i (t)=e j(ωt+φ) as the input.

T m (ω) Noise transmission function of a CMI unit pre-amplifier version 1 for a CMI FPA

LIST OF FIGURES

FIG. 1 is The Principles of the Method of Correlated Modulation Imaging (CMI)

a. CMI with passive modulator

b. CMI with active modulator

c. Unit cell readout circuit

FIG. 2 is Relative Noise Power Spectra of Integrators

Solid Curve: T m 2 (ω) of a CMI unit pre-amplifier used in CMI FPA

Dashed Curve: T 2 (ω) of a single cycle integrator used in a conventional FPA (Total integration time for both cases: T i =10 ms; modulation frequency f m =1 KHz)

FIG. 3 is Schematic of the CMI Circuitry with CMI Unit Pre-amplifier Version 1

FIG. 4 is Schematic of the CMI Circuitry with CMI Unit Pre-amplifier Version 3

FIG. 5 is The Layout of VLSI Design of the Silicon Chip with CMI Circuitries

FIG. 6 is Computer Simulation Results of the CMI Unit Pre-amplifier Version 3

FIG. 7 is Computer Simulation Results of the CMI Unit Pre-amplifier Version 1

FIG. 8 is Experimental Results of the CMI Unit Pre-amplifier Version 1

FIG. 9 is Improvement of CMI in Terms of Minimum Detectable Signal to Background Ratio

OBJECTS OF THE INVENTION

An object of the invention is the provision of a new method—correlated modulation imaging (CMI)—and the steps and means to realize the method for the detection of weak image signals buried in strong background.

A further object of the invention is the provision of two versions of circuit designs for the performance of the CMI unit pre-amplifier.

DETAILED DESCRIPTION OF THE INVENTION

During the course of this description, like numbers will be used to identify like elements according to different figures which illustrate the invention.

1. The Principles of the Method of CMI

The schematic of the CMI method is shown in FIGS. 1 ( a, b , and c ). Object 10 , which can be a person, a scene, a device, or any other object for imaging, reflects or emits light. The light coming from object 10 is focused by using a lens system 12 to form an image on the focal plane array (FPA) 13 , which is composed of a one-dimensional or two-dimensional array of pixels or unit cells 32 . For a two-dimensional CMI focal plane array, lens system 12 is stationary. For a one-dimensional CMI focal plane array, lens system 12 with a rotating mirror can be used.

Each unit cell is a photodetector (also named as detector) 16 , which is connected to its CMI unit pre-amplifier 31 composed of demodulator 29 and integrator 30 . The signal photocurrent I s is generated in detector 16 by the light from a corresponding area of object 10 . In addition to signal current I s , there are also background current I b generated by background illumination 14 and/or dark current I d due to detector leakage. I b and I d are independent from signal current I s , and play the same role. In the following discussion I b is used to represent the effects of both I b and I d .

In the method of correlated modulation imaging (CMI), the signal current is modulated, while the background and/or dark current are not modulated. Two methods can be used to modulate the signal current—by using a passive modulator ( FIG. 1 a ) or an active modulator ( FIG. 1 b ). The most widely used passive modulator 11 is a mechanical chopper or electro-optic switch, which is placed in the pathway between object 10 and FPA 13 . The most widely used active modulator is a pulsed laser 15 that illuminates object 10 . Passive modulator 11 or active modulator 15 of the FPA 13 , as well as all the unit pre-amplifiers ( 31 ) of the FPA 13 , are controlled by the same correlated controller 28 .

The sum of the modulated signal current I s and the unmodulated background and/or dark current I d is 17 , or I in , which corresponds to a unit cell 32 of FPA 13 , is fed to its corresponding demodulator 29 . The function of demodulator 29 is to offer correlated switch 22 that controls the direction of 17 , or I in . After passing switch 22 , by using integrator 30 , I in is integrated on 24 , an integration or charge storage capacitor C. The voltage accumulated on 24 is 27 , or the output voltage V out , which is fed to a standard amplifier through a standard bus controlled by a standard shift register used in conventional FPA's.

In phase one (φ 1 in FIG. 1 c ) input current I in ( 17 ) is the signal photocurrent I s plus the background (due to unmodulated background illumination 14 ) and/or dark current I d (generated by detector 16 without any illumination). In phase two (φ 2 ), there is no signal current I s . In the case of passive modulation, object 10 is blocked by modulator 11 . In the case of active modulation, object 10 does not emit any light since active modulator or pulsed laser 15 does not illuminate object 10 . Therefore, I in =I d . In CMI, the modulation frequency of correlated controller 28 is selected so that in each cycle of CMI the charge storage capacitor 24 of the integrator 30 is charged only to a portion of the saturation level. Discharged by the dark current I d repeatedly, capacitor 24 is saved for the integration of signal current I s . Thus, C or the charge storage capacitor 24 is not saturated quickly, but reserved for the integration of signal current I s .

2. The Theory or Theoretical Advantages of CMI

a. Improvement of FPA Sensitivity and Dynamic Range (Subtraction of DC Dark Current and Increase of Integration Time)

The sensitivity of a-semiconductor photodetector is characterized by its signal to noise ratio R defined as $$ R = Signal / Noise = N s n ro 2 _ + n 1 / f 2 _ + n th 2 _ + n s 2 _ + n d 2 _ ( 3 )

where N s is the number of integrated photoelectrons generated by a steady signal photon flux, and the denominator is the root mean square (rms) value of the total number of noise electrons. The noise electrons are from various sources. {overscore (n ro 2 )} is the noise due to readout electronics, {overscore (n 1/f 2 )} the 1/f or low frequency noise mostly associated with the fabrication process of the photodetector, {overscore (n th 2 )} the detector thermal noise, and n s 2 and {overscore (n d 2 )} the shot noises associated with the steady signal and dark currents, respectively. When the dark current (DC) is several orders higher than the signal current (DC), the shot noise n d 2 associated with the DC dark current may be predominant among all the components of the noise electrons. Therefore, Eq (3) can be simplified as $$ R = Signal / Noise ≈ N s n d 2 _ = N s N d = I s e ⁢   ⁢ I d ⁢ T i ( 4 )

where N d is the number of integrated electrons due to the DC dark current I d , e the electron charge, T i the signal integration time (approximately equal to the frame period for a staring FPA). Here we use the equation $$ n d 2 _ = N d = I d ⁢ T i e ( 5 )

which is valid for photovoltaic detectors. If the detectors are photoconductors, such as QWIP's, the rms value

$$ n d 2 _

of dark current generated shot noise electrons as expressed by Eq. (5) will be multiplied by a factor of square root of 2. Other equations in the following discussion will be modified with a similar factor of correction, which does not affect substantially our conclusions concerning the method of CMI.

It is readily seen that the signal to noise ratio is proportional to the square root of the total integration time T i of each pixel. However, T i is limited by N sat , the saturation or maximum number of electrons that an integration capacitor can handle. Note that N sat is limited by two factors. First, the capacitance of the signal integration capacitor can only be made in the range of pF due to real estate limit in the FPA; secondly, the increase of the capacitance will induce a higher kTC noise. For a conventional FPA of a single-cycle signal integration, where the capacitance trans-impedance amplifier (CTIA) is the most widely used readout interface, $$ T i = e ⁢   ⁢ N sat I s + I d ≈ e ⁢   ⁢ N sat I d ( 6 )

Using signal integration time defined by Eq. (6), we can obtain the optimized signal to noise ratio for a conventional FPA using CTIA $$ R Conv , max = ( Signal / Noise ) Conv , max ≈ I s I d ⁢ N sat ( 7 )

In CMI, the integration of the DC dark current is self-cancelled. Therefore, the integration time T i is maximized, and so is the signal to noise ratio: $$ T i ≈ e ⁢   ⁢ N sat I s ( 8 )

The maximum achievable signal to noise ratio is $$ R CMI , max = ( Signal / Noise ) CMI , max ≈ I s ⁢ N sat I d ( 9 )

Note that in deriving equation (9) it is assumed that the signal integration continues during the second half of integration period when the DC dark current of the previous half period is being cancelled out. This ideal assumption cannot be realized in the case discussed here: the signal current is modulated while the dark current is not. Thus, the maximum signal to noise ratio is reduced to $$ R CMI , max = 1 2 ⁢ I s ⁢ N sat I d ( 10 )

As the sensitivity of image detection of an FPA is defined by

Signal/Noise=1  (11)

the weakest detectable image signal current of a FPA using the CMI method is

( I s /I d ) min, CMI =4 /N sat   (12)

The typical value of N sat being in the order of 10 6 to 10 7 , therefore, the weakest detectable signal photocurrent I s is in the order of 10 −6 I d . Compared.with the weakest detectable image signal photocurrent of a FPA using the conventional single-cycle integration method $$ ( I s / I d ) min , Conv = 1 N sat ( 13 )

the weakest detectable signal of CMI has an improvement of more than 3 orders.

In addition to signal to noise ratio and weakest detectable signal, another important figure of merit of a focal plane array is its dynamic range. In terms of decibels, the dynamic range D of a conventional staring FPA is $$ D Conv = 20 ⁢   ⁢ log 10 ⁢   ⁢ R Conv , max = 20 ⁢   ⁢ log 10 ⁡ ( I s I d ⁢ N sat ) ( 14 )

whereas for a FPA using CMI method it is $$ D CMI = 20 ⁢   ⁢ log 10 ⁡ ( 1 2 ⁢ I s ⁢ N sat I d ) ( 15 )

The improvement of dynamic range is $$ D CMI - D Conv = 20 ⁢   ⁢ log 10 ⁡ ( 1 2 ⁢ I d I s ) ( 16 )

b. Reduction of Low Frequency Noise (Spectrum of the Transmission Function of CMI Unit Pre-amplifier Peaking at ω m )

As shown in FIG. 1 , in the method of CMI both the signal current (modulated) and the dark current (not modulated) are fed to the CMI unit pre-amplifier 31 . The noise current generated in the detector also goes through the CMI unit pre-amplifier. To study the output noise spectrum, we consider the output voltage V o (ω, φ) with a unit current component i (t)=e j(ωt+ω) as the input. The correlated switching functions as a multiplier m (t), which is a harmonic function or other more complicated functions. The most widely used one is a series of square waves depicting the effect of switching: $$ m ⁡ ( t ) = {   ⁢ 1 k ⁢   ⁢ τ < t < ( k + 1 ⁢ / ⁢ 2 ) ⁢ τ   - 1 ( k + 1 ⁢ / ⁢ 2 ) ⁢ τ < t < ( k + 1 ) ⁢ τ k = 0 , 1 , … ⁢   , m ⁢ - ⁢ 1   ⁢ 0 otherwise   ( 17 )

where τ is the period of each integration period, f m =1/τ the frequency of modulation, m the number of cycles of modulation as well as integration.

T i =m τ  (18)

as the total integration time of the CMI unit pre-amplifier version 1. The Fourier transform of m (t) in frequency domain is given by $$ M ⁡ ( ω ) = 1 2 ⁢ π ⁢ ∫ - ∞ ∞ ⁢ m ⁡ ( t ) ⁢ ⅇ - jω ⁢   ⁢ t ⁢   ⁢ ⅆ t = j πω ⁢ tan ⁡ ( ω ⁢   ⁢ τ 4 ) ⁢ sin ⁡ ( m ⁢   ⁢ ω ⁢   ⁢ τ 2 ) ⁢ ⅇ - j ⁢   ⁢ m ⁢   ⁢ ωτ / 2 ( 19 )

Note that m (t) can also be expressed as the inverse Fourier transform of M (ω) $$ m ⁡ ( t ) = ∫ - ∞ ∞ ⁢ M ⁡ ( ω ) ⁢ ⅇ jω ⁢   ⁢ t ⁢   ⁢ ⅆ ω = j π ⁢ ∫ - ∞ ∞ ⁢ tan ⁡ ( ωτ 4 ) ⁢ sin ⁡ ( m ⁢   ⁢ ωτ 2 ) ⁢ ⅇ jω ⁡ ( t - m ⁢   ⁢ τ / 2 ) ω ⁢   ⁢ ⅆ ω ( 20 )

Since m (t)=0 for 0>t≧T i , we can extend integration time of the input current from (0, T i ) to (−∞, ∞) so that Fourier transform can be used $$ V 0 = - 1 C ⁢ ∫ 0 Ti ⁢ ⅈ ⁡ ( t ) ⁢ m ⁡ ( t ) ⁢   ⁢ ⅆ t = - 1 C ⁢ ∫ - ∞ ∞ ⁢ ⅈ ⁡ ( t ) ⁢ m ⁡ ( t ) ⁢   ⁢ ⅆ t ( 21 )

where C is the capacitance of the integrator's charge storage capacitance. The output voltage for an input current e j(ωt+φ) is $$ V o ⁡ ( ω , φ ) =   ⁢ - j π ⁢   ⁢ C ⁢ ∫ - ∞ ∞ ⁢   ⁢ ⅆ ϖ ⁢ tan ⁡ ( ϖτ 4 ) ⁢ sin ⁡ ( m ⁢   ⁢ ϖτ 2 ) ϖ ⁢ ⅇ j ⁡ ( φ - m ⁢   ⁢ ϖτ / 2 ) ⁢ ∫ - ∞ ∞ ⁢ ⅇ j ⁡ ( ω + ϖ ) ⁢ t ⁢   ⁢ ⅆ t =   ⁢ 2 ⁢ j C ⁢ tan ⁡ ( ω ⁢   ⁢ T i 4 ⁢ m ) ⁢ sin ⁡ ( ω ⁢   ⁢ T i 2 ) ω ⁢ ⅇ j ⁡ ( ω ⁢   ⁢ T i / 2 + φ ) ( 22 )

The phase φ of the input noise current being random, we can obtain the rms value of the output voltage by averaging V o (ω, φ) through the phase φ. Thus, the noise transmission function of the CMI unit pre-amplifier version 1 is $$ T m ⁡ ( ω ) = V 0 ⁡ ( ω , φ ) ⁢ V 0 ′ ⁡ ( ω , φ ) = T i C ⁢ &LeftBracketingBar; tan ⁡ ( ω ⁢   ⁢ T i 4 ⁢   ⁢ m ) ⁢ sin ⁡ ( ω ⁢   ⁢ T i 2 ) ω ⁢   ⁢ T i 2 &RightBracketingBar; ( 23 )

FIG. 2 shows the distribution of the spectrum of relative noise power, which is the square of the noise transmission function (23). For comparison, also shown in FIG. 2 is the square of the noise transmission function (24) of a single-cycle integrator used in a conventional FPA $$ T ⁡ ( ω ) = V 0 ⁡ ( ω , φ ) ⁢ V 0 * ⁡ ( ω , φ ) = &LeftBracketingBar; 2 ⁢ sin ⁡ ( ω ⁢   ⁢ T i 2 ) &RightBracketingBar; ω ⁢   ⁢ C ( 24 )

To verify the results of Eq. (23) and (24), we calculate the total number of electrons due to shot noise. As a white noise, the shot noise associated with dark current I d in a photovoltaic diode has the power spectrum of

w ( f )=2 eI d   (25)

In CMI, through the total integration time Ti, the rms value of the number of shot noise electrons accumulated in the storage capacitor of the CMI unit pre-amplifier is $$ n d 2 _ = V n 2 _ ⁢ C e = C e ⁡ [ ∫ 0 ∞ ⁢ T m 2 ⁡ ( ω ) ⁢ w ⁡ ( f ) ⁢   ⁢ ⅆ f ] 1 / 2 = C e ⁡ [ 2 ⁢ w π ⁢   ⁢ C 2 ⁢ ∫ 0 ∞ ⁢ tan 2 ⁡ ( ωτ 4 ) ⁢ sin 2 ⁡ ( m ⁢   ⁢ ωτ 2 ) ω 2 ⁢   ⁢ ⅆ ω ] 1 / 2 = C e ⁢ ( wT i 2 ⁢ C 2 ) 1 / 2 = I d ⁢ T i / e = N d ( 26 )

In the case of a single cycle integrator, for the same integration time T i , the rms value of the number of shot noise electrons accumulated in the storage capacitor is also $$ n d 2 _ = V n 2 _ ⁢ C e = C e ⁡ [ ∫ 0 ∞ ⁢ T 2 ⁡ ( ω ) ⁢ w ⁡ ( f ) ⁢   ⁢ ⅆ f ] 1 / 2 = I d ⁢ T i / e = N d ( 27 )

Both Eq. (26) and (27) are identical to the widely used Eq. (5). As expected, the CMI method does not reduce the total output rms value of the number of shot noise electrons, since shot noise is completely random, and not correlated with the modulation and correlated integration utilized in the CMI.

Comparing the two spectra of FIG. 2 , however, we note that the transmission function T (ω) of a single cycle integrator is analogous to a low pass filter with bandwidth equal to ½T i , while that of the CMI unit pre-amplifier is analogous to a band pass filter peaking at the center of modulation frequency ω m , with satellite windows centered at the odd harmonics of ω m . An apparent advantage of CMI then lies in its capability of suppressing the flicker or 1/f noise, which, among all the noise sources, is usually important, and even dominant for some types of widely used detectors, such as HgCdTe IR photodetectors.

Besides 1/f noise, other low frequency noises may also play important roles. For example, in solar magnetography or magnetic field imaging, the atmospheric turbulence induced transmission fluctuation has a typical cut-off frequency of 1 kHz. If the image signal modulation frequency is higher than 1 kHz, the noise due to atmospheric turbulence may be eliminated by the method of CMI.

3. The Circuit Designs of CMI

The critical component for the functioning of the method of CMI is the CMI unit pre-amplifier 31 ( FIG. 1. c ). Two versions of CMI unit pre-amplifier are presented in this patent filing.

FIG. 3 is the first type of circuitry designed for CMI technology—CMI unit pre-amplifier version 1, which comprises two stages. The second stage S 2 is integrator 30 (the same numbering as in FIG. 1. c ), the capacitive transimpedance amplifier (CTIA) widely used in the readout circuit of infrared (IR) FPA's. In this stage, 41 or feedback capacitor C 2 of operational amplifier 43 is the main integrator, which accumulates signal throughout multi-cycle integration. 42 is the reset switch M 4 enabled at beginning of integration with a narrow pulse, and reset after readout of the integrated signal. The first stage S 1 is a special demodulator 29 (the same numbering as in FIG. 1. c ). In this stage, 35 or capacitor C 1 is a temporal integrator. M 1 , M 2 , and M 3 are the three MOSFET switches or gates of integrator 35 . The function of this stage is to integrate input 34 or input current I in from detector 16 (the same numbering as FIG. 1. c ), and to transfer the accumulated charges across 35 a and 35 b to C 2 with opposite polarities depending on the gate timings. When M 1 , M 3 are off and M 2 is on (Φ 2 in FIG. 1. c ), the input current I in from detector 16 is integrated on C 1 first. At the end of this half period of integration, by turning off M 2 followed by turning on M 1 and M 3 , the charges stored on 35 b , the right plate of C 1 is transferred to 41 a , the left plate of C 2 . At the same time, C 1 is reset. Assuming that the charges stored on 35 b are Q 1 , the change of output voltage V o at node 41 after the transfer is ΔV o1 =−Q 1 /Q 2 . At the next phase or half period of integration, with M 1 and M 2 off and M 3 on (Φ 1 in FIG. 1. c ), input current I in from detector 16 charges C 1 and C 2 . Assuming at the end of this period of integration the charges accumulated on 35 b or the right plate of C 1 is Q 2 , the change of output voltage V o at node 41 after the transfer is ΔV o2 =Q 2 /C 2 . By modulating the signal current I s and not modulating the background and/or dark current I d , we have Q i =Q d and Q 2 =Q d +Q s , where Q s is the integrated charge due to signal current I s , and Q d the integrated charge due to background and/or dark current I d . Therefore, the total change of output voltage V o at node 41 after the two consecutive integrations is ΔV o =Q s /C 2 . Note that for some type of signal modulation ΔV o =2Q s /C 2 . The function of the CMI method—integration of signal and cancellation of background—is thus fulfilled.

FIG. 4 is the second type of circuitry designed for CMI technology—CMI unit pre-amplifier version 3, which comprises four stages. The first stage S 1 converts the input current I in of detector 50 (the same as 16 in FIG. 1. c .) to voltage. S 1 is composed of 54 (an operational amplifier OP 1 ) and 53 (a feedback resistor with resistance of R). Feedback resistor 53 is connected to output node 55 of S 1 and input node 51 , the negative terminal of OP 1 . The positive terminal of OP 1 is connected to grounded node 52 . The output voltage of S 1 is therefore equal to −I in R. The second stage S 2 inverts the output voltage at node 55 of S 1 to current. It is composed of 59 (a second operational amplifier OP 2 ) and 56 and 60 (two resistors with equal resistance R). Since the negative input node 57 of OP 2 is virtually grounded, the current passing through 56 and 60 is −I in R/R=−I in . Consequently, the output voltage of S 2 at node 60 is I in R. Controlled by the phase φ of correlated controller 28 of FIG. 1. c , the third stage S 3 either passes the voltage I in R or the inverted voltage −I in R to the input of the fourth stage S 4 . S 3 is a 2-1 multiplexer composed of 63 (an n-channel MOSFET switch M 1 ) and 64 (a p-channel MOSFET switch M 2 ). When Φ is high (Φ 1 in FIG. 1. c ), with M 1 turned on and M 2 turned off, voltage I in R passes through this stage. When Φ is low (Φ 2 in FIG. 1. c ), with M 2 on and M 1 off, voltage −I in R passes through this stage. The fourth stage S 4 is a CTIA integrator with 71 , the third operational amplifier OP 3 , as its main component. The input resistor 66 of S 4 first turns the output voltage of S 3 back to the original input current I in , the direction of which is controlled by Φ. When Φ is high, the current I in flows in from 65 to 67 . When Φ is low, the current I in flows out from 67 to 65 . The reconverted current I in is integrated on 70 , the feedback or integration capacitor C. The change of the output voltage of S 4 is $$ - 1 2 ⁢ I i ⁢   ⁢ n ⁢ τ

when Φ is high, and $$ 1 2 ⁢ I i ⁢   ⁢ n ⁢ τ

when Φ is low, where τ is the period of modulation and integration. Assuming that during the first half period of τ, Φ is low, we have I in =I d . During the next half period of τ, with Φ turned high, I in =I d +I s , where I d is the unmodulated background and/or dark current, I s the modulated signal current. The net change of output voltage ΔV o at node 72 is $$ 1 2 ⁢ I d ⁢ τ - 1 2 ⁢ ( I d + I s ) ⁢ τ = - 1 2 ⁢ I s ⁢ τ

after one period of modulation. Note that for some type of signal modulation $$ Δ ⁢   ⁢ V 0 = 1 2 ⁢ ( I d - I s ) ⁢ τ - 1 2 ⁢ ( I d + I s ) ⁢ τ = - I s ⁢ τ .

The function of the CMI method—integration of signal and cancellation of background—is thus fulfilled.

FIG. 5 is part of the layout of VLSI design of the silicon CMOS chip, which has various CMI circuitries, including single pixels, and linear and area arrays with multiplexers, as well as the necessary device parameter test fixtures. The figure shows one pixel of the correlated modulation imager, including both CMI unit pre-amplifier version 1 and 3. The chip was designed using HP's 0.5μ CMOS design rules, and fabricated at HP's 0.5μ CMOS facility, which was managed by the Foundry of MOSIS. A total of 20 chips sized 1 mm×1 mm were delivered and tested.

4. The Circuit Simulation of CMI

FIG. 6 is a computer simulation result of the CMI unit pre-amplifier version 3, by using the circuit simulation package HSPICE. 25 cycles of integration, with a total integration time of 5 ms, are shown in the figure. With the level of dark current as shown in FIG. 6 , the maximum integration time will be less than 0.5 ms if a single cycle integrator is used.

FIG. 7 is a computer simulation result of the CMI unit pre-amplifier version 1 by using the same circuit simulation package HSPICE. 15 cycles of integration, with a total integration time of 30 ms, are shown in the figure. During the odd number of milliseconds (1st, 3rd, 5th, . . . ), or the first phase of integration, the charge accumulation on C 1 due to integration of the signal current plus the dark current, as well as the charge transfer from C 1 to C 2 , are shown in the plotting. During the even number of milliseconds (2nd, 4th, 6th, . . . ), or the second phase of integration, only the dark current is integrated, directly on C 2 and in the opposite direction. The function of correlated integration of modulated signal current and cancellation of background and/or dark current is achieved as designed.

5. The Test Results of CMI

FIG. 8 is experimental results of the CMI unit pre-amplifier version 1. The Pulse Instrument Focal Plane Array Test Equipment Setup is used to control the input and output of the Integrator. The modulation and correlated integration frequency of {fraction (1/32)}×10 6 Hz is generated by a pulse generator. With 30 cycles, a total time of 960μs is used for integration. The data of test results were fed to a tracer. FIG. 8 is a copy of the photograph taken from the screen of the tracer. The output of the integrator of the CMI unit pre-amplifier has the same features as the simulation results shown in FIG. 7 . It is therefore demonstrated experimentally that the CMI unit pre-amplifier version 1 functions as designed.

In the experimental test of the operation of correlated modulation imaging (CMI), a near infrared (NIR) light emitting diode (LED) powered by a DC source generated the background photocurrent I b of 2.4×10 −8 A, which was equivalent to a dark current I d in the NIR photodiode of the same magnitude. A second NIR LED, which was connected with a programmable square wave power source, generated modulated signal photocurrent I s several orders lower than I b . FIG. 9 shows the CMI unit pre-amplifier version 1 output voltage V o of the NIR photodiode as a function of the input I s /I d . Note that with the measured output error bar due to random fluctuation being estimated as 0.035 V, it is obvious that the measured weakest detectable signal I s /I d

( I s /I d ) min,CMI,Measured <5×10 −6   (28)

The integration capacitor of each pixel is 2 pF, and $$ N sat = CV sat e = 3.75 × 10 7 ( 29 )

in which we assume V sat =3 V. From Equations (12) and (13), we can obtain theoretical limits of the weakest detectable signal currents for the CMI unit pre-amplifier version 1 of a correlated modulation imager (CMI) and for the single-cycle integrator of a conventional FPA, respectively.

( I s /I d ) min, CMI, Theoretical =4 /N sat =1.07×10 −7   (30)

$$ ( I s / I d ) min , Conv , Theoretical = 1 N sat = 1.63 × 10 - 4 ( 31 )

A comparison of (28) and (31) shows that under the same background (or dark current), the weakest measured signal using CMI method is almost two orders smaller than the theoretical limit of the measurable signal using the conventional single-cycle integration FPA. Since signal to noise ratio and dynamic range of a FPA are directly related to its weakest detectable signal, therefore we conclude that we have experimentally demonstrated that compared to a conventional FPA, a FPA using CMI technology can improve its signal to noise ratio and dynamic range substantially. The theoretical limit of Eq. (30) was not achieved, since we used only a portion of the available saturation output voltage of 3 V to avoid the saturation, which is caused by slow drifting of CMOS device parameters and other instabilities of the electronics involved. We expect that the theoretical limit as depicted by Equation (30) will be approached if we can improve the stability and uniformity of our device and electronics, and as a result can use more cycles of integration and longer total integration time.

Claims

1. A circuit design of CMI unit pre-amplifier using one operational amplifier to perform the function of correlated modulation imaging (CMI), which is composed ofa. a first stage S2 including an integrator (30), said integrator (30) comprising an operational amplifier (30), a feedback capacitor C (24), and a reset switch M5 (42), wherein the capacitor C (24) accumulates signal throughout multi-cycle integration, while the reset switch M5 enabled at beginning of integration with a narrow pulse, and reset after readout of the integrated signal; and b. a second stage S1 including a demodulator (29), said demodulator comprising a temporal integrator C1 (35) and three MOSFET switches M2, M3, and M4 being gates of integrator (35); when M2 and M3 are off and M4 is on, an input current Iin (15) is integrated on C1 first; at the end of the half period of integration, by turning off M4 followed by turning on M2 and M3, the charges stored on the right plate (35b) of C1, will be transferred to the left plate (24a) of C, while at the same time, C1 is reset; assuming that the charges stored on said right plate (35b) are Q1, the change of output voltage Vo at node (27) after the transfer is ΔVo1=−Q1/C; at the next phase or half period of integration, with said M3 and M4 off and M2 on, said input current Iin charges said C1 and C; assuming at the end of this period of integration the charges accumulated on the right plate (35b) of C1 is Q2, the change of output voltage Vout at node (27) after the transfer is ΔVo2=Q2/C; thus, the total change of output voltage Vout at node (27) after the two consecutive integrations is ΔVout=Qs/C and the function of the CMI method, integration of signal and cancellation of background, is fulfilled.

2. A circuit design of CMI unit pre-amplifier using 3 operational amplifiers to perform the function of correlated modulation imaging (CMI), which is composed of:a. a first stage S1, which is composed of an operational amplifier OP1 (54) and a feedback resistor with resistance of R (53); said first stage converts an input current Iin to voltage, and with said resistor connected to the output node and the node of the negative terminal of said OP1, the output voltage of S1 is therefore equal to −IinR; b. a second stage S2, which is composed of a second operational amplifier OP2 (59) and two resistors with the same resistance R (56, 60), said second stage inverts the output voltage at node (55) of said first stage S1, and renders an output voltage of S2 at node (61) IinR; c. a third stage S3 controlled by the phase φm of a correlated controller (28), said third stage S3 including a 2-1 multiplexer composed of an n-channel MOSFET switch M1 and a p-channel MOSFET switch M2, either passes the voltage IinR when said φm is high, or the inverted voltage −IinR when said φm is low; d. a fourth stage S4, which is a CTIA integrator, gives an output voltage of -12⁢Ii⁢ ⁢n⁢τ when φm is high, and 12⁢Ii⁢ ⁢n⁢τ when φm is low, where τ is the period of modulation and integration, and therefore a net change of output voltage Δ⁢ ⁢V0=12⁢Id⁢τ-12⁢(Id+Is)⁢τ=-12⁢Is⁢τ after one period of modulation, the function of the CMI method, integration of signal and cancellation of background, is thus fulfilled.