The field of the invention is vision systems that acquire and display a series of image frames, and particularly, vision systems such as those used to enhance night vision.
Night vision systems employ cameras that are sensitive to wavelengths such as infrared (“IR”). The image frames that are produced may be displayed alone in sequence to provide a moving picture of a night-time scene captured by the camera, or the IR image frames may be combined with image frames acquired at optical wavelengths to enhance the optical image frames.
Night vision systems typically employ heavy image intensifier systems. Although image intensifiers can increase brightness and provide modest improvements in detected signal-to-noise ratio (“SNR”) through amplification or spectral changes in the light incident on a detector, they cannot improve on the SNR provided by the light emanating from the low-light-level scene. Thus, such vision systems characteristically produce images having a high noise content. These vision systems are also characterized by difficulties in depicting moving subjects and a changing scene that results from panning of the vision system, for example, from head motion. Such motion makes it difficult to filter noise from the image frames without losing significant information.
It would therefore be desirable to provide a lightweight, optionally non-intensified, system for night vision having an improved signal-to-noise ratio (“SNR”) without the loss of significant information in the image. Additionally, it would be desirable to provide a method for enhancing the visualization of static and moving objects in low light conditions.
The aforementioned drawbacks are overcome by the present invention, which provides a method for displaying images to a user using motion adaptive frame integration with real time digital processing and display. This allows for a lightweight, usually non-intensified, imaging system for night vision.
The present invention is a method for filtering a series of image frames that depict a moving subject, which thereby improves the signal-to-noise ratio (“SNR”) of each image frame. More specifically, a filtered, or “composite,” image is formed by combining, or “integrating,” pixel values in a current image frame with pixel values in previously processed image frames (that is, previously produced composite images) and pixels values in image frames acquired subsequent to the acquisition of the current image frame. This is achieved using a so-called bi-directional recursive filter, in which the contributions from the previous composite images and subsequent image frames to the filtering of the current image frame are weighted with a decay constant. Motion of the imaging system is compensated for by establishing a moving reference frame indicative of the moving imaging system and shifting the image frames to account for this bulk motion, or “panning.” Thus, the pixels in the image frames are registered to the same pixel locations before filtering the current image frame.
It is an aspect of the invention to provide a method for further processing a series of image frames that depict a moving subject such that motion that occurs within the image scene is compensated. For example, when the image scene corresponds to one that would be seen from a moving vehicle, substantially every portion of each image frame is moving at different rates in the image scene. To this end, each image frame is first transformed into the moving reference frame of the imaging system and then motion correction parameters associated with motion occurring within the image scene are determined. In this manner, each image frame is compensated for parallax and the variable speeds with which objects move through the scene.
It is another aspect of the invention to provide a method for filtering a series of image frames that depict a moving subject, which thereby improves the signal-to-noise ratio (“SNR”) of each image frame. More specifically, a filtered, or “composite,” image is formed by combining, or “integrating,” pixel values in a current image frame with pixel values in previously processed image frames (that is, previously produced composite images) and pixels values in image frames acquired subsequent to the acquisition of the current image frame. This is achieved using a so-called bi-directional recursive filter, in which the contributions from the previous composite images and subsequent image frames to the filtering of the current image frame are weighted with a decay constant. Moreover, at least some of the previous composite image frames and subsequent image frames are masked such that only the substantially static portions of these image frames are included in the composite image. The non-masked image frames are weighted with a first decay constant before combining them in the filtering operation and the masked image frames are weighted with a second decay constant before combining them in the filtering operation. This further incorporation of information into the current composite image frame being produced provides a corresponding further increase in SNR in the composite image.
It is yet another aspect of the invention to provide a method for filtering and displaying an image to a user with an imaging device, in which the filtered image has a significantly higher signal-to-noise ratio (“SNR”) than any single image frame acquired with the imaging system. This filtering is achieved by combining, or integrating, pixel values with a method that accounts for complications arising from situations where the imaging system is panning a scene and objects within the scene are moving.
It is yet another aspect of the invention to provide a method for enhancing and displaying an image to a user with an imaging device, in which the enhanced image identifies pixels in the image that are indicative of objects moving within the scene. A motion mask is produced by first establishing a moving reference frame indicative of the motion of the imaging system, and then transforming the acquired image frames into the moving reference frame. In effect, this registers the series of image frames with each other. The motion mask is produced by comparing temporally adjacent image frames and thresholding the comparison to identify those pixels that differ from one image frame to the next. A binary mask is produced using these identified pixels and is then temporally and spatially dilated to expand the region of pixels included in the mask. This mask is then convolved and normalized to blur the edges of the identified pixel regions. The acquired image frames are then enhanced using this produced motion mask. Enhancement occurs by producing a motion image for each image frame by multiplying each image frame by the corresponding motion mask. A static mask is then produced by subtracting the motion mask from a matrix of ones. This static mask is applied to the corresponding composite image frame to produce a “static image frame.” The static and motion image frames are then added together to produce a motion enhanced image frame.
The motion and static image frames are essentially different channels of the enhanced motion image frame, and these channels can be altered to produce different contrast therebetween. For example, the colormap of the motion image can be altered such that the motion image is highlighted with red, while the static image can retain a grayscale colormap. Thus, the motion enhanced images presented to the user will indicate those regions indicative of motion occurring within the scene as highlighted in red.
The foregoing and other aspects and advantages of the invention will appear from the following description. In the description, reference is made to the accompanying drawings which form a part hereof, and in which there is shown by way of illustration a preferred embodiment of the invention. Such embodiment does not necessarily represent the full scope of the invention, however, and reference is made therefore to the claims and herein for interpreting the scope of the invention.
Referring particularly to
Referring now particularly to
The video signal produced by the image array 200 is transmitted to the image processor 204, which receives and processes the video signal in accordance with the present invention. The processed video signal is then communicated to the display device 206, which receives the processed video signal and generates a two dimensional image that can be viewed by the user's eye 208. Exemplary display devices include flat panel microdisplays, which typically measure less than 1.5 inches diagonally. Such microdisplays include emissive, reflective, and transmissive displays. Emissive displays are preferred for the present configuration since they offer the greatest power efficiency; however, reflective and transmissive display types can similarly be employed. Exemplary emissive displays include light emitting diode (“LED”) and organic LED (“OLED”) display panels, while exemplary reflective and transmissive displays include a reflective surface and backlit liquid crystal display (“LCD”) panels, respectively.
A system controller 210 is communicatively connected to the image array 200. The system controller 210 includes, for example, field-programmable gate array (“FPGA”) logic gates that direct the image array 200 to operate in a number of different modes.
The image array 200, system controller 210, and video display device 206 receive electrical energy from a power source 212. Exemplary power sources 212 include batteries, such as one or more “AA” type batteries, or a rechargeable battery system. Exemplary rechargeable battery systems include thin film rechargeable lithium batteries, lithium ion batteries, and high efficiency nano-structured solid-state composite batteries, such as those with a lithium-based cathode and a silicon nanowire coated carbon-black or stainless steel anode.
The image processor 204 provides the capability of performing real-time image processing of the video signal for the purpose of enhancing the picture presented to the user in accordance with the present invention. The image processor 204 is communicatively connected to a system memory 214, such as, for example, random access memory (“RAM”). The memory 214 provides a buffer in which image data is stored as the image processor 204 performs its functions. The image processor also optionally receives motion and orientation information from an integrated motion detector 216 and orientation detector 218. Exemplary motion detectors include three-axis accelerometers, while exemplary orientation detectors include three-axis gyroscopes and three-axis compasses. The information provided to the image processor 204 by these detectors significantly improves the computational efficiency of the image processing methods described herein. However, the information provided by these detectors is not necessary when practicing the present invention. The image processor 204 also carries out conventional image processing functions, such as, for example: contrast stretching, edge detection and enhancement, aperture correction, and mixing of video channel information.
Exemplary image processors include: one or more low power microcontroller units (“MCU”), such as ARM MCUs (ARM Limited, Cambridge, England, UK) and Atmel AVR32 media processors (Atmel Corporation, San Jose, Calif.); digital signal processing (“DSP”) chips, such as TI DaVinci series DSPs (Texas Instrument, Dallas, Tex.) and Analog Devices Shark or Blackfin DSPs (Analog Devices, Cambridge, Mass.); Freescale i.MX multimedia applications processors (Freescale Semiconductors, Inc., Austin, Tex.); field-programmable gate arrays (“FPGA”); and application specific integrated circuits (“ASIC”).
The foregoing components are held in relative position to one another, and are protected from the outside environment, by the display housing 104. In an alternate configuration, and referring particularly now to
While an exemplary imaging system has been presented for utilization when practicing the present invention, other alternative imaging systems can also be readily adapted to employ the present invention. For example, a vision system such as an AN/PVS-14 monocular night vision device can be utilized when practicing the present invention with the appropriate modifications to the image processing capabilities thereof. It should also be appreciated by those skilled in the art that the aforementioned imaging system can be operated to acquire images sensitive to two separate spectral ranges of light. For example, the imaging system can be adapted and configured to acquire images over mid to long infrared wavelengths for thermal imaging. In such an application, a second dedicated imaging array 200 for imaging the second spectral range can be employed. In this manner, images over the separate spectral ranges can be acquired, processed, and presented to the user, or, in the alternative, the images can be combined to provide a fusion image.
Referring particularly to
Referring still to
After the image frames have been filtered to enhance their edge detail, a moving reference frame is established from the image frames, as indicated at step 304. The moving reference frame corresponds to motion of the imaging system FOV, such as, for example, head motion or vehicle motion. For example, if a user is wearing a “goggle system” such as the one described above with respect to
Subsequently, a motion mask is produced using the image frames, as indicated at step 306. The motion mask, Mk, includes pixel locations in the corresponding image frame, Ik, that correspond to motion that is occurring within the image scene. For example, if the FOV corresponds to an image scene in which an animal is running across a street at night time, the motion of the animal is identified across image frames and the pixels in each image frame that correspond to the animal are placed in the motion mask. An exemplary method for producing a motion mask is detailed below with respect to
After the moving reference frame has been established and the motion mask produced, a composite image is produced for presentation to the user, as indicated at step 308. A bi-directional truncated recursive filter is employed to produce the composite image. Referring now to
The selected image frames are co-registered and weighted using a decay constant to produce the composite image, Ck, corresponding to the current image frame, Ik. This recursive filter is performed on a pixel-by-pixel basis as follows:
and when previous composite image frames are utilized in lieu of previous image frames:
where Ck(x,y) is the value of the current composite image, Ck, at the pixel location (x,y); Ik(x,y) is the value of the current image frame, Ik, at the pixel location (x,y); λ is a decay constant; Ĩk−n(x,y) is the value of the nth previous image at the pixel location (x,y) that has been registered with the current image frame, Ik; {tilde over (C)}k−n(x,y) is the value of the nth previous composite image at the pixel location (x,y) that has been registered with the current image frame, Ik; and Ĩk+n(x,y) is the value of the nth succeeding image frame at the pixel location (x,y) that has been registered with the current image frame, Ik. The registration of the previous image frames and the proceeding image frames is performed using the shift values calculated when establishing the moving reference frame in step 304 above and using shift values corresponding to motion that occurs within the image scene, as described below with reference to
Ĩk+n(xk,yk)=Ik+n(xk+n+Δxk+n+δxk+n,yk+n+Δyk+n+δyk+n) Eqn. (3);
Where Δxk+n and Δyk+n are the global shift values for the nth succeeding image frame determined above in step 304 and δxk+n and δyk+n are the local shift values indicative of motion occurring within the nth succeeding image frame, as determined below with reference to
Ĩk−n(xk,yk)=Ik−n(xk−n+Δxk−n+δxk−n,yk−n+Δyk−n+δyk−n) Eqn. (4);
where Δxk−n and Δyk−n are the global shift values for the nth previous image frame determined above in step 304 and δxk−n and δyk−n are the local shift values indicative of motion occurring within the nth previous image frame, as determined below with reference to
The motion mask, Mk, produced above in step 306 can be employed during the production of the composite image, Ck, to allow more information of the static portion of the image scene to be included in the composite image to thereby provide a further increase in SNR. By way of example, and referring now to
MIk−n=(1−Mk−n)·Ik−n Eqn. (5);
where Mk−n is the motion mask produced for the nth previous time frame and 1 is a matrix of ones having the same size as the motion mask, Mk−n. The selected image frames and masked image frames are co-registered and weighted using two decay constants to produce the composite image, Ck, corresponding to the current image frame, Ik. For example, a first decay constant, λ1, is applied to the non-masked image frames 602 and 604, whereas a second decay constant, λ2, is applied to the masked image frames 652 and 654. This recursive filter is performed on a pixel-by-pixel basis as follows:
where, λ1 is a first decay constant; λ2 is a second decay constant; and, as above, {tilde over (C)}k−n(x,y) is the value of the nth previous composite image at the pixel location (x,y) that has been registered with the current image frame, Ik, and Ĩk+n(x,y) is the value of the nth succeeding image frame at the pixel location (x,y) that has been registered with the current image frame, Ik.
It should be appreciated by those skilled in the art that the same decay constant can be employed for both λ1 and λ2. In addition, it should be appreciated that instead of weighting the masked imaged in Eqn. (6) by nλ2, a uniform weighting can be applied as simply λ2. Instead of masking image frames to isolate those regions that are substantially static, only those images that contain substantially no dynamic regions can also be utilized in the aforementioned process. For example, those image frames whose absolute value of displacement is under a threshold value can be included in the bi-directional recursive filter. Likewise, an ongoing sum of the shift values can be formed and the average of all previous frames whose absolute value of the displacement sum is under a threshold value can alternatively be utilized in the filter.
Referring again to
It should be apparent that the magnitudes of the pixel values in this weighted composite image do not accurately indicate the actual pixel values in the scene. If the actual pixel values are to be preserved, an additional normalization is added to the processing. This normalization step is performed by dividing the weighted composite image, Ck, by a blurred version of the composite image, Ck. The blurring of the composite image, Ck, is done in the same manner as that described above for blurring the image frame, Ik. While the normalization step preserves the actual pixel values in the weighted composite image, it adds considerable processing time which can unduly reduce the frame rate at which image frames can be produced.
Exemplary further processing of the composite image, Ck, also includes applying the motion mask produced above in step 306 to isolate and identify those pixels in the composite image that correspond to motion occurring within the image scene. The masking of the composite image, Ck, in this manner is performed in accordance with the following:
Ĉk=(1−Mk)·Ck+Mk·Ik Eqn. (7);
where Ĉk is the masked composite image, Mk is the motion mask, and 1 is a matrix of ones having the same size as the motion mask, Mk. Thus, the pixels corresponding to motion within the image scene are removed from the current composite image, Ck, through multiplication by the static binary mask, (1−Mk), and the pixels not corresponding to motion within the image scene are removed from the current image frame, Ik, through multiplication by the motion binary mask, Mk. These two intermediate images are then combined to produce the masked composite image, Ĉk.
In this manner, the masked composite image, Ĉk, can be viewed as having two channels: a motion channel corresponding to Mk·Ik and a static channel corresponding to (1−Mk)·Ck. These two channels can be modulated to provide a unique contrast therebetween. For example, the static channel can be presented to the user using a gray scale color map while the motion channel can be presented to the user using a “red scale” color map so that those pixels corresponding to motion that has occurred within the image scene are highlighted as red in the masked composite image, Ĉk.
When displaying the composite image, Ck, or masked composite image, Ĉk, to the user, the influence of the image array 200 detector noise can be removed. A calibration step is performed, in which a plurality of image frames are acquired with the imaging system 100 while a lens cap covers the optical assembly, 102. For example, 50 image frames can be acquired in this manner. This arrangement is sufficient to characterize the inherent dark current and bias noise response of the image array 200. These “calibration” image frames are subsequently averaged and the result is stored in the system memory 214 as a “detector noise” map. This detector noise map is then subtracted from each image frame, Ik, prior to producing the composite image, Ck. It should be appreciated by those skilled in the art that a separate detector noise map can be produced for each unique set of settings utilized when operating the imaging system 100. In this manner, the detector noise map corresponding to a particular set of image system settings is subtracted from those image frames produced using the same settings.
The inherent dark current and bias noise of the image array 200 can additionally be characterized by exposing the image array 200 to an optical device such as a uniform intensity light source. An exemplary uniform intensity light source is a so-called integrating, or Ulbricht, sphere. This process results in the acquisition of a so-called uniform intensity image (“UII”). The UII is employed to produce a flat field image, for example, by subtracting the detector noise map from the acquired UII in order to compensate the UII for dark current and bias noise. When utilizing the additional characterization provided by the flat field image, after the detector noise map is subtracted from each image frame, Ik, the result is divided by the flat field image. In this manner, significant pixel structure noise is removed from the displayed composite image, Ck, or masked composite image, Ĉk. As noted above, it should be appreciated by those skilled in the art that a unique detector noise map may be produced for a particular set of imaging system settings. To this end, a particular flat field image can also be produced for the same particular set of imaging system settings when the detector noise map corresponding to those settings is employed during the generation of the flat field image.
Referring now particularly to
After the region that is to be registered has been selected, a roll-off filter is applied to the image frames, as indicated at step 402. This filter is applied so that edge effects that occur when information enters or leaves the FOV in sequential image frames are mitigated, thereby maintaining accuracy in the registration process. An exemplary roll-off filter, Froll-off, has the form:
where α is constant defining the roll-off filter size, Nx is the number of pixel columns in the image frame, Ny is the number of pixel rows in the image frame, and (x,y) indicates the pixel location in the image frame that is being filtered.
After the image frames have been filtered as described above, a gross panning correction is performed to determine initial motion shift values, as indicated at step 404. This is achieved by searching over a coarse shift grid. These motion shift values are indicative of the direction and amount of bulk motion in the FOV, thereby defining an initial “global shift vector”. Such bulk motion is typically caused by movement of the image system, such as panning, from one image frame acquisition to the next.
This is done by determining the direction and distance that the current image frame, Ik, has shifted in position from the position depicted in the previously acquired image frame, Ik−1. A methodical search is made over a coarse grid in the previous image frame, Ik−1, for a region of pixel values that correspond to a central region of pixel values in the current image frame, Ik. The region of pixels is smaller than the image frames, but large enough to detect bulk motion of the subject. During the search the sum of the absolute differences between corresponding pixel values in the two regions is calculated for each candidate region location in the previous image frame, Ik−1. The region location which produces the lowest difference value is chosen and this yields the initial global shift vector, which indicates the amount and direction the subject has shifted within the field of view of the imaging system. Using this initial global shift vector, all of the pixels in the previous image frame, Ik−1, can be shifted in the direction needed to offset the bulk motion. In doing so, the previous image frame, Ik−1, is then said to be “registered” with the current image frame, Ik. However, a refined search is performed around this initial global shift vector to more accurately establish the moving reference frame.
As indicated above, a fine panning correction is performed next to determine more accurate motion shift values, as indicated at step 406. This is achieved by performing shifts on a fine shift grid around the initial shift values calculated above in step 404. Utilizing the initial shift values, a methodical search is performed over a fine grid in the previous image frame, Ik−1, for a region of pixel values that correspond to a central region of pixel values in the current image frame, Ik. As its name implies, the fine grid is finer than the coarse grid employed when determining the initial shift values; thus, a search is performed to find smaller regions of pixels in the current image frame, Ik, that correspond to pixels in the previous image frame, Ik−1. The result of this search is a set of shift values defining a global shift vector that describes how the overall scene has shifted from a previous image frame, Ik−1, and the current image frame, Ik. It is not important that these shifts be exact, since motion that occurs within the image scene, and therefore within in each frame, is subsequently determined. By establishing the moving reference frame, however, this subsequent search is made more efficient by providing an initial estimate of gross motion.
Referring particularly now to
In carrying this out, a block array in the sequential image frames is selected, as indicated at step 702. Exemplary block sizes within this block array are 8×8 or 10×10 pixels. A search is then performed in the registered previous image frame, Ik−1, within the selected block array for the location of the corresponding pixel in the current image frame, Ik, that is being processed. This is done by comparing the values of the pixels in the block surrounding the current pixel being processed in a manner similar to the comparison performed in step 404 described above.
By way of example, and referring particularly now to
Referring now to
The size of the search pattern 806 will depend on the amount of anticipated subject motion from frame-to-frame, which in turn will depend on the particular application. Similar to the global motion correction previously performed, the foregoing local motion detection is repeated using a finer displacement grid on which the search pattern is performed, as indicated at step 710. As before, this comparison is performed by summing the absolute differences between corresponding pixel values in the two block 808 and 810 and storing the resulting “regions difference value” in an updated data table 752. This updated data table 752 includes a shift matrix, which relates the block shift values to the true pixel location shifts in (x,y)-coordinates. This finer local motion correction is subsequently repeated for each pixel location, as indicated at step 712 until the search pattern is completed, as indicated at decision block 714.
Because the image frames are inherently noisy, this noise can affect the calculated shift values, suggesting that there is motion where, in fact, no motion occurred. To mitigate such extraneous shift values, the shift matrix is filtered, as indicated at step 716. An exemplary filtering operation is achieved using a median filter.
Referring particularly now to
Because the central portions of moving objects may be fairly uniform from image frame to image frame, the initial motion mask will often contain pixels only from the leading and trailing edge of the moving object. Thus, in order to produce a motion mask that includes the entire moving object, the initial motion mask is expanded by temporal summation over a specified number of image frames, as indicated at step 506. It will be appreciated by those skilled in the art that other methods for expanding the initial motion mask can be employed, including a spatial dilation binary operation. The expanded binary mask is subsequently convolved to fill in potential undetected motion pixels, as indicated at step 508. The convolved motion mask includes non-unity values, and the pixel values in the mask fall off toward the edge of the mask as a result of the convolution. Therefore, the pixels in the convolved mask are normalized by dividing each pixel in the convolved mask by the maximum value in the convolved mask, as indicated at step 510. This provides a rolled off version of motion information when the current time frame information is multiplied by the motion mask and combined with information from the non-motion regions, as described above.
The present invention has been described in terms of one or more preferred embodiments, and it should be appreciated that many equivalents, alternatives, variations, and modifications, aside from those expressly stated, are possible and within the scope of the invention. While the present invention has been described in the context of a head-mounted camera application, such as a night vision device, it should be apparent that other applications are also possible. More particularly, the invention is applicable when a series of image frames are acquired and processed with a composite image that is produced using a bi-directional recursive filter during the acquisition of image frames, in which bulk motion of the scene within the FOV of each image frame is present along with detailed subject motions. Other such applications include video and medical imaging systems in which such bulk motion within the FOV is a concern. It should be appreciated by those skilled in the art that the present invention can be readily adapted to such exemplary alternate applications.
This application is a continuation of U.S. application Ser. No. 12/538,433, filed Aug. 10, 2009 now U.S. Pat. No. 8,497,914, which issued on Jul. 30, 2013.
This invention was made with government support under EB006393 awarded by the National Institutes of Health. The government has certain rights in the invention.
Number | Name | Date | Kind |
---|---|---|---|
5442407 | Iu | Aug 1995 | A |
5798948 | Rashkovskiy et al. | Aug 1998 | A |
6373055 | Kerr | Apr 2002 | B1 |
6442203 | Demos | Aug 2002 | B1 |
6707044 | Lannestedt et al. | Mar 2004 | B2 |
7103235 | Tener et al. | Sep 2006 | B2 |
7110035 | Hellstrand | Sep 2006 | B1 |
7128705 | Brendley et al. | Oct 2006 | B2 |
7336823 | Guldevall | Feb 2008 | B2 |
7340162 | Terre et al. | Mar 2008 | B2 |
7391450 | Pinto et al. | Jun 2008 | B2 |
20020101510 | Basson | Aug 2002 | A1 |
20020159651 | Tener et al. | Oct 2002 | A1 |
20040252199 | Cheung et al. | Dec 2004 | A1 |
20090022410 | Haskell | Jan 2009 | A1 |
Entry |
---|
McIvor et al., The Background Subtraction Problem for Video Surveillance Systems, pp. 176-183, 2001. |
Adirondack Video Astronomy, Adirondack Flat Field Users Manual [online], 2005. Retrieved from the Internet: <URL: http://www.astrovid.com/technical—documents/AVAflatfielder.pdf>. |
Bolte, Signal-to-Noise in Optical Astronomy [online], 2006. Retrieved from the Internet: <URL: http://www.ucolick.org/˜bolte/AY257/s—n.pdf>. |
Gardner, Characterizing Digital Cameras with the Photon Transfer Curve [online], 2008. Retrieved from the Internet: <URL: http://www.couriertronics.com/docs/notes/cameras—application—notes/Photon—Transfer—Curve—Charactrization—Method.pdf>. |
Heap and Brown, Correction for the STIS Echelle Blaze Function, pp. 114-119, 1997. |
Tulloch, Design and Use of a Novel Flat Field Illumination Light Source, 1996. |
Robinson, Determining the Vignetting and Echelle Blaze Function for the GHRS, pp. 291-303, 1993. |
Kirschman, et al., Flat-field calibration of CCD detector for Long Trace Profiler, Proc. SPIE, vol. 6704, 2007. |
Number | Date | Country | |
---|---|---|---|
20130300883 A1 | Nov 2013 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 12538433 | Aug 2009 | US |
Child | 13938521 | US |