WO2001020898A1 - Method and system for altering defects in a digital image - Google Patents

Method and system for altering defects in a digital image Download PDF

Info

Publication number
WO2001020898A1
WO2001020898A1 PCT/US2000/025478 US0025478W WO0120898A1 WO 2001020898 A1 WO2001020898 A1 WO 2001020898A1 US 0025478 W US0025478 W US 0025478W WO 0120898 A1 WO0120898 A1 WO 0120898A1
Authority
WO
WIPO (PCT)
Prior art keywords
pixel
channel
digital image
pixels
filtered
Prior art date
Application number
PCT/US2000/025478
Other languages
French (fr)
Inventor
Albert D. Edgar
Raymond Shaw Lee
Original Assignee
Applied Science Fiction
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Applied Science Fiction filed Critical Applied Science Fiction
Priority to AU75868/00A priority Critical patent/AU7586800A/en
Priority to EP00965090A priority patent/EP1221256A1/en
Publication of WO2001020898A1 publication Critical patent/WO2001020898A1/en

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/40Picture signal circuits
    • H04N1/409Edge or detail enhancement; Noise or error suppression
    • H04N1/4097Removing errors due external factors, e.g. dust, scratches
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/46Colour picture communication systems
    • H04N1/56Processing of colour picture signals
    • H04N1/58Edge or detail enhancement; Noise or error suppression, e.g. colour misregistration correction

Definitions

  • This invention relates generally to image processing and more particularly to a method and system for altering defects in a digital image.
  • Tangible images such as photographic images, may have surface defects such as scratches, fingerprints, or dust particles. Such defects may occur, in the case of photographic images, in a transparency or negative as well as in a photographic print of a transparency or negative. Such defects often undesirably degrade a photographic image.
  • a defect channel comprising a digital signal proportional to the defects in a photographic image may be created by scanning the photographic image using an infrared light source and an infrared light sensor. Infrared light will tend to pass through developed photographic film with nearly complete transmission because the dye in various layers of the photographic film does not fully absorb infrared light. On the other hand, where defects are present, a portion of the infrared light will tend to be refracted from the optical path before passing through the film.
  • a defect channel may be obtained by examining the difference between images obtained when the image being scanned is illuminated by light sources at different angles. The challenge is to use the defect channel to automatically alter defects in a digital image, while making as few undesirable changes to the digital image as possible.
  • One aspect of the invention is a method for altering defects in a digital image. At least a first pixel of a first channel of a digital image is filtered using digital circuitry to produce a filtered pixel by averaging the intensity of the first pixel and a plurality of additional pixels in the neighborhood of the first pixel . At least one of the pixels is weighted in response to the intensity value of at least one pixel in a defect channel associated with the digital image. A corrected digital image is produced using the digital circuitry in response to the filtered pixel and the first channel of the digital image.
  • the invention has several important technical advantages. Various embodiments of the invention may have none, one, some, or all of these advantages without departing from the scope of the invention.
  • the invention allows automatic alteration of defects in a digital image based upon a defect channel having a signal proportional to defects in the digital image.
  • the invention allows such alteration in the image substantially with little or no damage to the underlying image. Because the invention filters pixels of the digital image with pixels weighted based upon the expected reliability of each pixel, areas of strong defect are more easily excised without causing significant damage to the area surrounding the defect. Thus, for reasonably sized defects, the invention may automatically alter defects in a way that will most often produce a more pleasing image than the original image that included the defect. The desirability of the resulting image may depend upon the size of the defect and the amount of valid image information remaining in the digital image.
  • FIGURE 1 illustrates a block diagram of a general purpose computer that may be used in accordance with the present invention
  • FIGURE 2 illustrates an example of a scanner that comprises an embodiment of the present invention
  • FIGURE 3 illustrates a flow chart describing the alteration of a defect in a digital image in accordance with one method of the present invention
  • FIGURE 4 illustrates a flow chart describing adjustment of an image in response to weighted averages produced by filtering an image in accordance with the present invention
  • FIGURE 5 illustrates a flow chart describing a second method of altering defects in a digital image in accordance with the present invention
  • FIGURE 6 illustrates a graph of an example weighting function that may be used with the present invention.
  • FIGURES 1-5 of the drawings like numerals being used for like and corresponding parts of the various drawings.
  • FIGURE 1 illustrates a general purpose computer 10 that may be used for image enhancement in accordance with the present invention.
  • general purpose computer 10 may comprise a portion of a digital image processing system and may be used to execute applications comprising image enhancement software.
  • General purpose computer 10 may be adapted to execute any of the well known MS-DOS, PC-DOS, OS2, UNIX, MAC-OS and Windows operating systems or other operating system.
  • General purpose computer 10 comprises processor 12, random access memory (RAM) 14, read only memory (ROM) 16, mouse 18, keyboard 20, and input/output devices such as printer 24, disk drives 22, display 26 and communications link 28.
  • the present invention includes programs that may be stored in RAM 14, ROM 16, or disk drives 22 and may be executed by processor 12.
  • Disk drive 22 may include a variety of types of storage media such as, for example, floppy disk drives, hard disk drives, CD ROM drives, or magnetic tape drives. Although this embodiment employs a plurality of disk drives 22, a single disk drive 22 could be used without departing from the scope of the invention.
  • FIGURE 1 only provides one example of a computer that may be used with the invention. The invention could be used with computers other than general purpose computers as well as general purpose computers without conventional operating systems.
  • General purpose computer 10 further comprises scanner 30 that may be used to scan images that are to be enhanced in accordance with the teachings of the invention.
  • enhancement may be performed by software stored and executed by scanner 30 with the results stored in a storage medium comprising a part of scanner 30 and/or in any of the storage devices of general purpose computer 10.
  • software for image enhancement may be stored in any of the storage media associated with general purpose computer 10 and may be executed by processor 12 to enhance images scanned by scanner 30.
  • image enhancement could occur both internally within scanner 30 and in general purpose computer 10 without departing from the scope of the invention.
  • Scanner 30 may comprise a film scanner or a flatbed scanner of any type without departing from the scope of the invention.
  • Image enhancement may also be performed using special purpose digital circuitry contained either in scanner 30, general purpose computer 10, or in a separate device. Such dedicated digital circuitry which may include, for example, state machines, fuzzy logic, etc .
  • FIGURE 2 illustrates an exemplary scanner 34 constructed in accordance with the invention.
  • Scanner 34 comprises processor 36, storage medium 38 and scanning hardware 40.
  • Processor 36 controls the operation of scanning hardware 40 by executing control software 44 stored in storage medium 38.
  • storage medium 38 may comprise multiple storage mediums as well as comprising storage mediums of different types.
  • control software 44 may be stored in ROM memory, RAM memory, or on a disk drive.
  • Scanning hardware 40 is used to convert an analog image into a digital image utilizing some type of optical circuitry.
  • optical circuitry may also be used to produce a defect channel proportional to defects in the analog image. Any type of optical circuitry could be used for scanning hardware 40 without departing from the scope of the invention.
  • a light source comprised mostly of energy outside the visible spectrum and a matching sensor may be used to create the defect channel.
  • a light source comprised mostly of energy outside the visible spectrum and a matching sensor may be used to create the defect channel.
  • an infrared light source and sensor such as those typically used in image processing applications involving photographic images may be used for this aspect of the scanning hardware.
  • a defect channel can be derived from a plurality of scanned versions of the image in the visible spectrum. Such a defect channel may be derived by illuminating the image being scanned at two or more angles and calculating changes in the scanned image at the plurality of angles. Defects will tend to affect the light differently when illumination is made from different angles. Other types of scanning hardware may be used to create a defect channel without departing from the scope of the invention.
  • scanning hardware 40 After scanning hardware 40 has scanned an image, that image may be enhanced (by altering defects within it) in accordance with the invention using image processing software 42, which is stored in storage medium 38.
  • image processing software 42 which is stored in storage medium 38.
  • the invention may employ other types of digital circuitry comprising any type of dedicated digital hardware to alter defects in the digital image.
  • This hardware may be a part of scanner 34 or general purpose computer 10 as discussed above.
  • the scanned image may be stored in storage medium 38 as may the enhanced image.
  • scanner 34 may not have any image processing software 42.
  • Such software instead may be provided on general purpose computer 10 for enhancement of an image received from scanner 34. Enhancement may occur both in scanner 34 and general purpose computer 10 as well.
  • a scanned image and/or an enhanced scanned image may be provided by scanner 34 to the general purpose computer 10 through a communications port (not explicitly shown) .
  • the defect channel may be similarly provided.
  • FIGURE 3 illustrates a flow chart describing a method employed by one embodiment of the present invention to enhance a digital image.
  • the image enhancement described herein may be carried out using computer software, as can any of the processes described below. That software, as discussed above, may be executed by scanner 34, by general purpose computer 10, or a combination thereof. A digital image received from other than scanner 30 may be enhanced in accordance with the invention.
  • the method described in FIGURE 3 may be used to alter defects in many types of images such as color photographic images (either negative print or transparency) , black and white images (either negative print or transparency and including black and white images derived from photographic film with multiple layers), other monochromatic images, x- rays or any other type of image stored on film.
  • the invention may also be used to alter defects in any image on a tangible medium that may be scanned using a scanner.
  • an image is scanned to create a digital image and defect channel.
  • this step could be omitted and the invention carried out on a image that has been previously scanned and has a defect channel associated with it.
  • the digital image will typically be comprised of three channels: a red channel , a green channel , and a blue channel .
  • Each channel is comprised of a series of pixels with each pixel having an intensity value associated with it corresponding to the intensity of the particular color of light at that spatial location in the original image.
  • Other types of color images can be used without departing from the scope of the invention.
  • the enhanced image with the altered defects could then be converted back into a color image.
  • the methods of the invention could also be applied to a single color channel of a digital image and the same correction could be applied to all channels of the digital image.
  • a black and white image may comprise one or more channels similarly made up of pixels.
  • Other types of images may also comprise one or more similar channels. The invention can be used for any of these types of images.
  • the defect channel comprises a series of pixels, each having an intensity value either directly proportional or inversely proportional to defects (or the lack thereof) in the original image.
  • a defect channel may be created, for example, using an infrared light source and infrared sensor of the kind commonly used in image processing applications involving photographic film. Other types of light sources and sensors may be used to create the defect channel, such as described above in connection with reflective scanners, for example. Any other method may be used to generate a suitable defect channel without departing from the scope of the invention.
  • Such a defect channel will ordinarily produce a signal having a stronger correlation to defects in the original image and a weaker correlation to the visible image itself.
  • Steps 52-58 comprise a method for altering defects in a digital image to produce an enhanced digital image.
  • the invention corrects defects in at least two frequency bands.
  • the term "correction” or “correct” when used in this application refers broadly to altering defects in a digital image. An image defect does not need to be completely removed or completely corrected to fall within the scope of the invention. Accordingly, defect correction includes, but is not limited to, reduction or other alteration of a defect in a digital image.
  • the frequency bands for the image are created using one or more filters which average pixels in the neighborhood of a pixel in question with each pixel weighted according to its expected reliability.
  • the expected reliability is determined by the intensity value in the defect channel of the pixel in question or the pixel in question plus a series of pixels in the neighborhood of the pixel in question. Pixels that appear to be more reliable based upon the defect channel are weighted more heavily than those with a low expected degree of reliability.
  • the defect channel may be similarly divided into frequency bands.
  • each weighted averaging operation is to create a low pass filtered version of the original image.
  • multiple low pass filtered versions of the original image are derived, most often with different bandwidths .
  • These multiple lowpass filtered versions, along with the original image may be used to derive a representation of the original image in multiple contiguous frequency bands.
  • a bandpass representation for one band may be derived where common frequencies are eliminated.
  • a series of bands can be created similarly, as further described below.
  • the filtered representations of the digital image may then be processed such that the original image is divided into a plurality of frequency bands which collectively make up all or substantially all frequencies present in the original image but overlap little, if at all. These bands may then be recombined to produce an image with defects removed .
  • the same filtering operation may optionally be performed on the defect channel .
  • Defects may then be further removed, frequency band by frequency band, by subtracting the defect channel bands from the image frequency bands.
  • An enhanced image is then constructed by combining the corrected individual frequency bands back together again.
  • step 52 data is optionally converted to log space for the defect channel and for each channel of the digital image where defect correction is to be performed.
  • the conversion of data to logarithmic space may enable the remaining steps of the method to be performed more easily than if the conversion were not made.
  • several of the operations described below involve additions and subtractions in log space, but could involve division by zero in non-log space. Because division by zero can lead to erroneous results, it may be more convenient to carry out the process of the invention in log space.
  • the performance of mathematically equivalent functions outside of log space could also be used in connection with the invention.
  • any type of logarithmic calculation can be used, one embodiment of the invention uses the conversion to log space described by Formula 1.
  • x represents the intensity value of the pixel to be converted to log space.
  • Each such pixel in this embodiment comprises a twelve-bit intensity value. Other numbers of bits can be used without departing from the scope of the invention.
  • convenient conversions may be made to take into account the capabilities of digital hardware used to perform calculations in accordance with the invention.
  • a first type of defect correction may be carried out by subtracting the defect channel from each of the visible image channels. Such a subtraction may be a direct subtraction or a bounded subtraction, analogous to the bounded subtraction described below in connection with Figure 5.
  • This subtraction may also take into account red leakage and clear film values where the defect channel is derived using an infrared source and sensor in a film scanner. Any method of taking into account these effects (which are described below) is within the scope of the invention. Defect channels obtained using infrared light in film scanners will often be suitable for this optional enhancement of the image. Even in such a case, however, this step is optional.
  • An alternative to this optional subtraction is to subtract the defect channel from the image channel, frequency band by frequency band, after enhancement by weighted filtering as described below. This alternative is also optional and may be done where a suitable defect channel is available.
  • the defect channel may be suitable where pixel intensity values in the defect channel are proportional, either directly or inversely, to defects and vary approximately linearly with the amount of light blocked by the defect.
  • the image and defect channel are filtered by taking weighted averages of pixels within varying distances from specific pixels.
  • each channel of the image may be filtered, the channels may be combined into a single channel and the combined channel filtered or a subset of the channels can be combined or filtered individually and the defect correction applied based upon results obtained from those channels without departing from the scope of the invention.
  • the filtering of the defect channel is also optional.
  • the weighting applied to a specific pixel in calculating the weighted averages may be based upon the expected reliability of that pixel.
  • the expected reliability of a pixel may be determined by using the defect channel.
  • a weighted average or “averaging” or any similar term refers to any type of average such as a median average, mode average, mean average, arithmetic average, or geometric average.
  • the calculation of weighted averages of pixels around the pixels in each channel of the digital image has the effect of filtering each channel with different strengths of low pass filters.
  • each pixel or a subset of the pixels involved in the calculation of the averages based upon the expected reliability of each pixel dampens the effect of the defect in each filtered version of the original channel of the digital image.
  • An example may illustrate the calculation of the weighted averages.
  • four filtering operations are performed: a 3x3 weighted average, a 5x5 weighted average, a 9x9 weighted average, and a 17x17 weighted average.
  • each of these four weighted averages is computed for each pixel in the digital image. Any pixel outside the boundaries of the image is set to zero for purposes of these calculations.
  • the 3x3 weighted average is computed using a 3x3 matrix of pixels with the pixel in question at the center of the matrix.
  • the weight applied to each pixel corresponds to its expected reliability as determined using a corresponding pixel or pixels in the defect channel.
  • the weighted average is computed by summing the product of the intensity of each pixel times the weight and dividing the total sum by the sum of the weight values that were applied to each pixel .
  • the 5x5, 9x9, and 17x17 filters are calculated similarly. The calculation of the four discussed weighted averages at each pixel results in four low pass filtered versions of the original channel of the digital image.
  • the 3x3 weighted average applied to each pixel of the image channel produces a low pass filtered version of the image having spatial frequencies between essentially zero and one-third of the maximum spatial frequency possible in the image.
  • the 5x5 weighted average produces a low pass filtered version of the original channel of the digital image having frequencies predominantly between zero and one- fifth of the maximum spatial frequency possible in the image.
  • the 9x9 and 17x17 weighted averages produce filtered versions with frequencies predominantly between zero and one-ninth the maximum spatial frequency possible in the image and predominantly between zero and one-seventeenth the maximum spatial frequency possible in the image. Due to the weighting applied to each pixel during this filtering operation, the filtered versions in each frequency band have had the effects of defects reduced.
  • Low pass filtering alone tends to dampen the effects of defects, but may also blur image detail .
  • the weighting may dampen the effects of a defect more than low pass filtering alone would with proportionally less reduction of image detail.
  • This example will be used below to further illustrate adjustment of the image and defect channel .
  • the same filtering operation may also be performed on the defect channel, if desirable for use in further reduction of defects.
  • the filter can have a square shape (as in the example above), a circle shape, or any other type of shape without departing from the scope of the invention.
  • the filter may be symmetric or asymmetric and may be centered or not centered on the pixel in question.
  • the filter may also be a window with feathered edges.
  • any number of filters can be used. The use of more filters will allow division of the image into a greater number of frequency bands. The use of additional filters, however, requires additional calculation. Less filters than those used in the above example can also be used without departing from the scope of the invention.
  • any pixel outside the boundaries of the image was set to an intensity value of zero for purposes of the weighted average calculation.
  • Other boundary conditions could be used without departing from the scope of the invention; for example, soft edge boundary conditions such as a triangle edge or gaussian edge could be used.
  • the weighted average was computed for each pixel in each channel of the digital image and in the defect channel.
  • the weighted average could also be performed on a subset of the pixels in a channel or on all of the pixels in a subset of the channels without departing from the scope of the invention.
  • different weights could be used for each frequency band and could be omitted in some frequency bands.
  • Different weights could also be used for different channels.
  • the weighting used for a particular pixel can be determined in several ways. The weighting could depend upon the intensity of a spatially corresponding pixel in the defect channel. Where the weighting is applied to the defect itself, the weighting could be determined based upon the intensity of each pixel involved in the calculation.
  • the same or different weightings could be used for the defect channel and one or more image channels.
  • the defect channel may be blurred compared to the visible channel due to focal shifts, the nature of the defects, and inconsistent registration between the visible channels of the digital image and the defect channel, for example. Accordingly, it may be desirable to estimate the reliability of a pixel and determine its weight by taking into consideration a plurality of pixels surrounding a pixel in the defect channel corresponding to the pixel in a visible image channel that is the subject of the weighted average calculation. Alternatively, if blurring effects are insignificant and registration error is fairly constant, then the pixel used in the defect channel could be chosen while compensating for the constant registration error.
  • the particular weighting function chosen may depend upon the characteristics of the particular scanner used to create the defect channel.
  • the weighting applied may be a function of the intensity value of a pixel or pixels in the defect channel . Where multiple pixels are used to determine a weight, an average or weighted average of the pixels may be used to come up with an average intensity used to determine a weight.
  • a weight may be assigned for all pixel values greater than a particular high threshold where the threshold provides high predictability that no defect is present.
  • a weight of zero may be assigned when the intensity value is below a certain low threshold indicating a high probability that a defect is present with little or no image detail remaining.
  • a straight line or other curve can be used to connect the two thresholds to establish a continuous function for weights between zero and one. Pixels with intensity in this middle region tend to be ones that are defective but some image detail remains to allow enhancement of the digital image to remove the defect but maintain some image detail .
  • the intensity of the pixels in the defect channel may lack a linear relationship to the amount of light that passes through a defect. With such a defect channel, a discrete set of weights such as 0, 1/2 and 1, may be used. An example of a weight function is illustrated in FIGURE 6.
  • infrared intensity in the defect channel will be high when no defect is present as most or all of the infrared energy is allowed to pass through the film. Where a defect is present, however, infrared intensity will be low.
  • all likely defective pixels (those with an infrared intensity in the defect channel less than 20 percent) are set to a weight of zero.
  • All pixels with a high probability of being nondefective pixels (those with an infrared intensity in the defect channel greater than 64 percent of the maximum intensity are assigned a weight of one) . Pixels in the region in between are assigned a weight based upon the illustrated curve. Again, more complicated functions may be used without departing from the scope of the invention.
  • the weighting applied may also compensate for effects such as leakage from the red channel of a color image into the defect channel .
  • This leakage may result because infrared light sources and/or sensors are residually sensitive to the cyan dye in an image used to modulate the red region of the visible spectrum. This effect manifests itself as an apparent leakage of the red image into the infrared image.
  • the effects of red leakage can be taken into account when establishing a weight value. For example, an overall red leakage value for an image may be calculated. This red leakage value can then be used to establish a constant to be multiplied times a pixel intensity value in the red channel . This product may represent an estimate of the amount of the red channel present in the defect channel for that particular pixel. Thus, in calculating a weight, this product may be subtracted from the intensity of the pixel in the defect channel .
  • the weight may also be adjusted by subtracting the average clear film value for a particular digital image.
  • Step 56 the channels of the digital image and the defect channel are adjusted in response to the weighted averages to lessen the effects of defects in the image.
  • the invention includes any use of the weighted averages computed in Step 54 to adjust an image to lessen the effects of defects in the image.
  • a specific method for adjusting the digital image in response to the weighted averages will be discussed below in connection with FIGURE 4. Such adjustment could occur in the time domain or the frequency domain and, when in the spatial domain, in log space or any other space.
  • another possible method of enhancement is to create a series of essentially contiguous frequency bands based upon the plurality of weighted averages of the original image and the original image itself.
  • the process of creating these bands may be as described below in connection with step 60 of FIGURE 4.
  • This method may be used, for example, where the defect channel is unsuitable for subtraction from the image. This method could, however, be used even where a defect channel is suitable for correction.
  • the enhanced image may also be further enhanced by other operations without departing from the scope of the invention. If the pixel data was converted to log space in step 52 , the enhanced pixel data may be converted back to the space of the original image using a suitable inverse formula such as an inverse to Formula 1.
  • FIGURE 4 illustrates one example of a method for adjusting images in response to weighted averages to lessen the effects of defects in the digital image.
  • the filtered digital image (which comprises the original digital image and the filtered versions of the digital image calculated in Step 54) and the filtered defect channel (which consists of the original defect channel and the filtered versions of the defect channel obtained in step 54) are separated into frequency bands in response to the calculated weighted averages.
  • the filtered images and the original image are used to create a series of frequency bands with little or no overlap representative of the original digital image and original defect channel with defective pixels suppressed due to the weighting that took place during the filtering operation. This method assumes that the defect channel was filtered in FIGURE 3.
  • five frequency bands may be created using the four weighted averages calculated.
  • F m represents the maximum spatial frequency that can exist in the digital image.
  • the five frequency bands corresponding to the channel in question may be created using the weighted averages computed for that channel.
  • the frequency band from approximately 1/3 F m to F m comprises the difference between the original channel and the 3x3 weighted average filtered version of that channel.
  • the frequency band from approximately 1/5 F m through 1/3 F m may be calculated by subtracting the 5x5 weighted average version of the channel from the 3x3 weighted average version of the original channel .
  • the frequency bands from approximately 1/5 F m to 1/9 F m and approximately 1/17 F m to 1/9 F m may be determined by subtracting the 9x9 weighted average version of the original channel from the 5x5 version of the original image channel and by subtracting the 17x17 weighted average version of the original channel from the 9x9 weighted average version of the original channel, respectively.
  • the frequency band from approximately zero to 1/17 F m is represented by the 17x17 weighted average version of the original image channel.
  • the low pass filtered versions of the image and defect channel created using the weighted averages can be used to divide each channel of the image as well as the defect channel into contiguous frequency bands where some or all of the frequency bands have had the effect of defects suppressed using the weighted average calculations.
  • the red residue in each frequency band of the defect channel may be removed.
  • such residue may be removed using a bounded calculation to allow some variance.
  • infrared sources and sensors or other sources and sensors
  • the step may be omitted without departing from the scope of the invention.
  • One option for removing the red leakage from the defect channel is to subtract from the intensity value of each pixel in the defect channel, the product of the intensity of the corresponding pixel in the red channel multiplied by an average red leakage constant representing the average red leakage for the entire digital image. This difference may be divided by the difference between one and the red leakage value to properly normalize the result.
  • red leakage may vary within regions of particular images, it may be useful to use a bounded subtraction to allow for some variance in localized red leakage values.
  • Step 64 defects are further removed from the image by subtracting the relevant frequency band of the defect channel obtained in Step 62 from the relevant frequency band in each image channel obtained in Step 60.
  • a bounded calculation such as that described in connection with Step 62 may be used to allow for some variance caused by registration error, blurring, etc. in the defect channel.
  • Step 66 the frequency bands are recombined through a summation to form an enhanced image with the original defect removed.
  • Step 67 the enhanced image is converted back from log space to the original space from which it was derived using a suitable inverse formula such as an inverse to Formula 1. If the enhanced image was created in a space other than log space, then Step 67 may be omitted without departing from the scope of the invention. If Step 60 through 66 were carried out in the frequency domain, then the enhanced image may be reconverted back to the time domain in Step 67.
  • FIGURE 5 illustrates a flow chart describing a second method of enhancing an image by altering defects in the image in accordance with the invention. In this embodiment, an image may be divided into segments and defects processed within each individual segment .
  • the image is divided into segments and a process similar to that discussed above in connection with FIGURES 3 and 4 is applied to each segment individually as if that segment comprised the entire image.
  • step 68 an image is scanned to create a digital image comprising one or more channels and a defect channel. All of the options discussed above in connection with step 50 are available for this embodiment of the invention.
  • step 70 the data from the defect channel and the channels of the digital image is converted to log space.
  • any of the options discussed above in connection with step 52 may be employed in step 70, including optional subtraction of the defect channel from the visible channel.
  • each channel of the digital image and the defect channel may be divided into segments. Any size or shape of segments may be used without departing from the scope of the invention. In one embodiment, each channel of the digital image and the defect channel are divided into 8x8 segments. This example will be used to illustrate the remainder of the steps of the method illustrated in FIGURE 5.
  • step 74 weighted averages of pixels within varying distances from specific pixels in a segment of the defect channel and a segment of each channel of the image are computed. These weighted averages are computed similarly to the weighted averages that were computed in step 54. However, the weighted averages computed in step 74 are computed for an individual segment of the image, assuming that the segment comprises the entire image. Accordingly, even if the 8x8 segment is surrounded by other 8x8 segments, the pixels beyond the boundaries of 8x8 segment in question are treated as beyond the image boundary and any pixel beyond those boundaries are treated as having an intensity of zero (or other boundary condition) .
  • step 74 All of the options for filtering and weighting discussed above in connection with step 54 may be employed in step 74 for the embodiment disclosed in FIGURE 5
  • an 8x8 filter is the maximum sized filter that will be used in this example.
  • the weight function takes into account red leakage and clear film effects as discussed above in connection with step 54.
  • the weight for a particular pixel for the weighted average calculations may be determined, for example, using Formula 2.
  • x and y represent the coordinates of a particular pixel within an 8x8 segment.
  • Di n represents the defect channel received from the scanner in step 68.
  • R L comprises a constant value representative of average red leakage in either the segment in question or in the entire digital image.
  • R ⁇ N represents the red channel of the digital image.
  • C F is a constant representing the average clear film value for either the segment in question or the entire digital image.
  • the weight function of Formula 2 is constrained such that the weight varies between zero and one. A plot of the weighting function versus the intensity in the defect channel would look similar to the curve illustrated in FIGURE 6. Any other weight function could be used without departing from the scope of the invention.
  • weighted averages are computed for each channel .
  • a 3x3 weighted average is computed.
  • an 8x8 weighted average is computed.
  • the 8x8 average covers the entire segment, and as a result, a single scalar value may be used to represent the result of this weighted average.
  • the 3x3 weighted average for the red, green, blue, and defect channels may be computed using Formulas 3 through 6. ⁇ ⁇ W(x + a,y+b)R m (x + a,y+b) , ,
  • R 3L (x,y) a - l b °-' , j where ⁇ ⁇ W(x + a, y+b) ⁇ 0
  • R ⁇ n , Gi n , B in , D ⁇ n represent the red, green, blue, and defect channels of the digital image that was scanned in step 68.
  • the 8x8 weighted average for the red, green, blue, and defect channels may be calculated using Formulas 7 through 10, respectively. As noted, each of these Formulas produces a scalar result.
  • R 8L (x,y) a ° 7 where ⁇ ⁇ W(a,b) ⁇ 0
  • R 3L represents a weighted low-pass filtered version of the red channel of the digital image having a frequency band between approximately zero and 1/3 F m , where F m is the maximum spatial frequency possible for the segment of the digital image.
  • G 3L , B 3L , and D 3L cover a similar frequency band for the green, blue, and defect channels respectively.
  • R 8L represents a weighted low-pass filtered version of the red channel of the digital image having a frequency band between approximately zero and 1/8 F m .
  • G 8L , B 8L , and D 8L represent weighted low-pass filtered versions of the green, blue, and defect channels, respectively in the same frequency band.
  • step 76 the segment of the image is adjusted to lessen the effects of defects in the image in response to the weighted averages of each image channel and the defect channel obtained in step 74 and the original image channels and defect channels. Any of the options discussed above in connection with step 56 and steps 60 through 67 may be used in connection such adjustment.
  • a version of the original image and original defect channel with effects of the defect partially suppressed due the weighted averages computed in step 74, is created in three separate frequency bands.
  • a first frequency band is obtained having frequencies between approximately zero and 1/8 F m .
  • a second frequency band extends from approximately 1/8 F m to 1/3 F m .
  • a third frequency band extends from approximately 1/3 F m to F m .
  • the first frequency band may be obtained using Formulas 11 through 14.
  • this frequency band was already calculated in step 74, it can be used without any further computation.
  • the third frequency may be computed using Formulas 15 through 18.
  • the weighted low-pass filtered version of each channel in the frequency range of approximately zero to 1/3 F m is subtracted from the unfiltered image channel to produce the frequency band between approximately 1/3 F m and F ra .
  • the second frequency band the one between approximately 1/8 F m and 1/3 F m may be computed using Formulas 19 through 22.
  • the 8x8 weighted filtered version of each image channel is subtracted from the 3x3 weighted filtered version of the original image channel.
  • the process just described for creating a plurality of frequency bands may employ any of the options discussed above in connection with step 60.
  • R 3B (x,y) R 3L (x,y)-R 8L (19)
  • G 3B (x,y) G 3L (x,y)-G 8L (20)
  • the red residue may be removed from the defect channel in each frequency band.
  • any of the options described above in connection with step 62 may be used without departing from the scope of the invention.
  • a bounded subtraction is used for two of the frequency band and a simple subtraction is used for the remaining frequency band.
  • Formula 23 may be used to subtract the red residue in the first frequency band of the defect channel .
  • D 8B D8B ⁇ R8BDK8RL (23) ⁇ -D K8 R L
  • D K8 the constant, D K8 , will normally be set to one, but may vary depending upon the scanner in question.
  • This constant represents the amount of the red channel that should be subtracted from the defect relative to the measured average red leakage for the entire segment and/or the entire image.
  • Formulas 24 through 26 may be used to subtract the red residue from the second frequency band of the defect channel .
  • D 3B ( x > y) T 1 D3B ( ⁇ .y) where
  • D 3B ( ⁇ 5 y) T2 D3B ( ⁇ ; y) where
  • the bounded subtraction provided for in Formulas 24 through 26 allows for local variance in the red leakage value, as well as for some registration error and/or blurring in the defect channel as compared to the visible channel .
  • Formula 24 multiplies a high value constant, D K3H against the red leakage constant, R L .
  • Formula 25 multiplies a low-value constant, D 3L , by the red leakage constant, R L .
  • D K3H and D K3L are the high and low ranges determining how much of the image content in the red channel should be subtracted from the defect channel . These constants may be determined experimentally. In this embodiment, D 3H s chosen to be 1.3 while D K3L is chosen to be 0.6. These constants will tend to average about one.
  • the average is a function of the resolution of the system for the infrared channel versus the red channel .
  • the spread between these two constants is a function of the accuracy of the scanner that produces the defect channel . If D K3H is chosen too large, then small image patterns may cause large matching defects to be erased, rendering defects uncorrected in middle frequencies when they are next to image detail. If D K3H is chosen too small, then image residue may remain in the defect record, causing middle frequency image detail to be erased. If D K3L is chosen too large, then too much image detail may be removed from the defect causing a negative residue that interferes with the ability to distinguish a defect. If D K3L is chosen too small, then visible detail may not be removed in the presence of defect detail of opposite polarity.
  • the revised version of the frequency band is calculated using Formula 26. If the results of Formula 24 and Formula 25 have different signs, then the intensity value is set to zero. If these signs are the same, then the lesser value produced by either Formula 24 or Formula 25 is chosen. Thus, this bounded subtraction tends to drive the revised value closer to zero.
  • Formulas 27 through 32 may be used to remove red residue in the defect channel for the third frequency band created above .
  • T3 D1B (x,y) Tl D1B (x,y) where
  • T3 D1B (x,y) T2 D1B (x,y) where
  • T3 D1B (x,y) 0 where T1 D1B (x,y)andT2 D1B (x,y)havedifferentsigns
  • T4 D1B (x,y) T3 D1B (x,y)-D K1A (30)
  • Dl B ( ,y) T4 D1B (x,y) where
  • D' 1B (x, y) T5 D1B (x, y) where
  • Formulas 27 through 29 are similar to Formulas 24 through 26.
  • the high- constant D K i H is chosen to be 1.0 while the low- frequency constant, D K ⁇ L is chosen to be 0.3. Similar considerations apply to the choice of these constants, as applied to the choice of D K3H and D K3L above. In this case, the effects of these constants will be in the high frequency band.
  • the average of the two constants is less than one and may be chosen as such to the extent that heightened frequency defects are blurred relative to the high frequency image. The range is also wider to accommodate a greater variance that often occurs at higher frequencies.
  • the difference between Formulas 27 through 32 and Formulas 23 through 26 involves the addition of Formulas 30 through 31.
  • the constant for the first frequency band, R K8 is set to one, but could be adjusted.
  • a direct subtraction is used, rather than a bounded subtraction. Blurring effects tend not to affect this frequency band or affect it insignificantly, but a bounded subtraction could be employed if desired.
  • the constants R K3H and R K3L are set to 1.5 and 0.5 respectively.
  • the constants RKIH and R KIL are set to 1.7 and 0.5, respectively. These values are also chosen for the corresponding constants for the green and blue channels. These constants establish the bounds within which correction is made. Again, one selects whichever number between the bounds drives the corrected result closet to zero.
  • the upper and lower bounds will normally average about one, except that one may set the lower bound a bit lower without damaging the image.
  • the defects in the original image may be suppressed even further by multiplying the pixels resulting from a bounded subtraction of the defect channel from the image channel by the weight for a particular pixel calculated using Formula 2.
  • the defects may be removed from each frequency band of the visible image channels, using Formulas 33 through 59.
  • Formulas 33 through 35 may be used to remove defect information from the first frequency band.
  • Formulas 36 through 47 may be used to remove the defect from the middle frequency band of each channel.
  • Formulas 48 through 59 may be used to remove the defect from the third frequency band of each channel .
  • R' 8B R 8B -D' 8B R K8 (33)
  • G 8B G 8B -D 8B G K8 (34)
  • TIMB (X, y) R 3 B ( ⁇ > y)-D 3B ( ⁇ > y) K3 H (36)
  • T R3B ( > y) Tl R3B (x,y) where
  • ⁇ 3 R3B ( ⁇ ,y) T2 R3B (x,y) where
  • T3 R3B (x,y) 0 where T2 R3B (x,y)andTl R3B (x,y) have opposite signs
  • R' 3B (x,y) T3 R3B (x,y)w(x,y) (39)
  • T3 G3B (x,y) Tl G3B (x,y) where
  • T3 G3B (x,y) T2 G3B (x,y) where
  • T3G G 3 3 B B (x,y) 0 where T2 G3B (x,y)andTl G3B (x,y)haveoppositesigns
  • G 3B (x,y) T3 G3B (x,y)w(x,y) (43)
  • Tl B3B (x,y) B 3B (x,y)-D 3B (x ⁇ B ⁇ (44)
  • T3 B3B (x > y) Tl B3B (x,y) where
  • T3 B3B (x,y) T2 B3B (x,y) where
  • T3 B3B (x,y) 0 where T2 B3B (x,y)andTl B3B (x,y)haveoppositesigns
  • T3 R1B (x,y) T2 R1B (x,y) where
  • T3 R1B (x,y) 0 where T2 R1B (x,y)andTl R1B (x,y) have opposite signs
  • R; B T3 R1B (x,y)w(x,y) (51)
  • T1 G1B (x,y) G 1B (x,y)-D; B (x,y) G K1H (52)
  • T3 G i B ( x > y) where
  • T3 GIB ( ⁇ ,y) T G i B ( x >y) where
  • T3 G i B ( x > y) °
  • T2 G1B (x,y)andTl G1B (x,y) have opposite sign
  • G; B T3 G1B (x,y)w(x,y) (55)
  • T3 B1B (x, y) T2 B1B (x, y) where
  • T3 B1B (x,y) 0 where T2 B1B (x,y)andTl B1B (x,y) have opposite signs
  • the frequency bands may be recombined for each channel using Formulas 60 through 62.
  • G out (x,y) G; B (x, y)+G' 3B (x, y)+G' 8B (x, y) (61)
  • step 78 it is determined whether there are any more segments of the digital image to process. If so, then steps
  • step 80 the method continues in step 80.
  • step 80 it is determined if there are more overlaps to perform. If not, then the procedure terminates in step
  • each channel of the enhanced digital image obtained by using steps 74 and 76 for each segment is divided into segments spatially different from the earlier segments. Steps 74 and 76 are then repeated for each of these segments.
  • This aspect of this embodiment may compensate for results that may be obtained for pixels near the boundary of segments that were used when the image was first divided into segments in step 72.
  • the defect may be further suppressed.
  • a second overlap is performed using 8x8 segments comprising a 2x2 corner of four adjoining segments that touch one another at a common point . These 8x8 segments are thus made up of one- fourth of each of four adjoining segments.
  • Any of the options described above with respect to step 72 may be used in step 82 in performing division of the enhanced image into different segments. Different shapes can be used and different size segments can be used without departing from the scope of the invention.

Abstract

One aspect of the invention is a method for altering defects in a digital image. At least a first pixel of a first channel of a digital image is filtered using digital circuitry to produce a filtered pixel by averaging the intensity of the first pixel and a plurality of additional pixels in the neighborhood of the first pixel. At least one of the pixels is weighted in response to the intensity value of at least one pixel in a defect channel associated with the digital image. A corrected digital image is produced using the digital circuitry in response to the filtered pixel and the first channel of the digital image.

Description

METHOD AND SYSTEM FOR ALTERING DEFECTS IN A DIGITAL IMAGE
TECHNICAL FIELD OF THE INVENTION
This invention relates generally to image processing and more particularly to a method and system for altering defects in a digital image.
BACKGROUND OF THE INVENTION
Tangible images, such as photographic images, may have surface defects such as scratches, fingerprints, or dust particles. Such defects may occur, in the case of photographic images, in a transparency or negative as well as in a photographic print of a transparency or negative. Such defects often undesirably degrade a photographic image.
In the field of image processing, digital images derived from photographic images using a scanner most often include the defects present in the underlying photographic image. Because digital images are subject to mathematical manipulation, if image defects may be identified and distinguished from image detail, then those defects can be removed, either partially or completely. A defect channel comprising a digital signal proportional to the defects in a photographic image may be created by scanning the photographic image using an infrared light source and an infrared light sensor. Infrared light will tend to pass through developed photographic film with nearly complete transmission because the dye in various layers of the photographic film does not fully absorb infrared light. On the other hand, where defects are present, a portion of the infrared light will tend to be refracted from the optical path before passing through the film. Thus, defects in the photographic image will tend to show up in a defect channel produced using an infrared light source and infrared sensor. In reflective scanners, a defect channel may be obtained by examining the difference between images obtained when the image being scanned is illuminated by light sources at different angles. The challenge is to use the defect channel to automatically alter defects in a digital image, while making as few undesirable changes to the digital image as possible.
SUMMARY OF THE INVENTION
One aspect of the invention is a method for altering defects in a digital image. At least a first pixel of a first channel of a digital image is filtered using digital circuitry to produce a filtered pixel by averaging the intensity of the first pixel and a plurality of additional pixels in the neighborhood of the first pixel . At least one of the pixels is weighted in response to the intensity value of at least one pixel in a defect channel associated with the digital image. A corrected digital image is produced using the digital circuitry in response to the filtered pixel and the first channel of the digital image.
The invention has several important technical advantages. Various embodiments of the invention may have none, one, some, or all of these advantages without departing from the scope of the invention. The invention allows automatic alteration of defects in a digital image based upon a defect channel having a signal proportional to defects in the digital image. The invention allows such alteration in the image substantially with little or no damage to the underlying image. Because the invention filters pixels of the digital image with pixels weighted based upon the expected reliability of each pixel, areas of strong defect are more easily excised without causing significant damage to the area surrounding the defect. Thus, for reasonably sized defects, the invention may automatically alter defects in a way that will most often produce a more pleasing image than the original image that included the defect. The desirability of the resulting image may depend upon the size of the defect and the amount of valid image information remaining in the digital image.
BRIEF DESCRIPTION OF THE DRAWINGS
For a more complete understanding of the present invention and the advantages thereof, reference is now made to the following descriptions taken in conjunction with the accompanying drawings in which:
FIGURE 1 illustrates a block diagram of a general purpose computer that may be used in accordance with the present invention;
FIGURE 2 illustrates an example of a scanner that comprises an embodiment of the present invention;
FIGURE 3 illustrates a flow chart describing the alteration of a defect in a digital image in accordance with one method of the present invention;
FIGURE 4 illustrates a flow chart describing adjustment of an image in response to weighted averages produced by filtering an image in accordance with the present invention;
FIGURE 5 illustrates a flow chart describing a second method of altering defects in a digital image in accordance with the present invention; and FIGURE 6 illustrates a graph of an example weighting function that may be used with the present invention. DETAILED DESCRIPTION OF THE INVENTION
The preferred embodiment of the present invention and its advantages are best understood by referring to FIGURES 1-5 of the drawings, like numerals being used for like and corresponding parts of the various drawings.
FIGURE 1 illustrates a general purpose computer 10 that may be used for image enhancement in accordance with the present invention. Specifically, general purpose computer 10 may comprise a portion of a digital image processing system and may be used to execute applications comprising image enhancement software. General purpose computer 10 may be adapted to execute any of the well known MS-DOS, PC-DOS, OS2, UNIX, MAC-OS and Windows operating systems or other operating system. General purpose computer 10 comprises processor 12, random access memory (RAM) 14, read only memory (ROM) 16, mouse 18, keyboard 20, and input/output devices such as printer 24, disk drives 22, display 26 and communications link 28. The present invention includes programs that may be stored in RAM 14, ROM 16, or disk drives 22 and may be executed by processor 12. Communications link 28 is connected to a computer network but could be connected to a telephone line, an antenna, a gateway, or any other type of communication link. Disk drive 22 may include a variety of types of storage media such as, for example, floppy disk drives, hard disk drives, CD ROM drives, or magnetic tape drives. Although this embodiment employs a plurality of disk drives 22, a single disk drive 22 could be used without departing from the scope of the invention. FIGURE 1 only provides one example of a computer that may be used with the invention. The invention could be used with computers other than general purpose computers as well as general purpose computers without conventional operating systems. General purpose computer 10 further comprises scanner 30 that may be used to scan images that are to be enhanced in accordance with the teachings of the invention. In this embodiment, enhancement may be performed by software stored and executed by scanner 30 with the results stored in a storage medium comprising a part of scanner 30 and/or in any of the storage devices of general purpose computer 10. Alternatively, software for image enhancement may be stored in any of the storage media associated with general purpose computer 10 and may be executed by processor 12 to enhance images scanned by scanner 30. In addition, image enhancement could occur both internally within scanner 30 and in general purpose computer 10 without departing from the scope of the invention. Scanner 30 may comprise a film scanner or a flatbed scanner of any type without departing from the scope of the invention. Image enhancement may also be performed using special purpose digital circuitry contained either in scanner 30, general purpose computer 10, or in a separate device. Such dedicated digital circuitry which may include, for example, state machines, fuzzy logic, etc .
FIGURE 2 illustrates an exemplary scanner 34 constructed in accordance with the invention. Scanner 34 comprises processor 36, storage medium 38 and scanning hardware 40. Processor 36 controls the operation of scanning hardware 40 by executing control software 44 stored in storage medium 38. Although a single storage medium has been illustrated for simplicity, storage medium 38 may comprise multiple storage mediums as well as comprising storage mediums of different types. Thus, for example, control software 44 may be stored in ROM memory, RAM memory, or on a disk drive. Scanning hardware 40 is used to convert an analog image into a digital image utilizing some type of optical circuitry. In addition, optical circuitry may also be used to produce a defect channel proportional to defects in the analog image. Any type of optical circuitry could be used for scanning hardware 40 without departing from the scope of the invention. For the defect channel, a light source comprised mostly of energy outside the visible spectrum and a matching sensor may be used to create the defect channel. For example, an infrared light source and sensor such as those typically used in image processing applications involving photographic images may be used for this aspect of the scanning hardware.
If scanner 34 comprises a reflective scanner, then a defect channel can be derived from a plurality of scanned versions of the image in the visible spectrum. Such a defect channel may be derived by illuminating the image being scanned at two or more angles and calculating changes in the scanned image at the plurality of angles. Defects will tend to affect the light differently when illumination is made from different angles. Other types of scanning hardware may be used to create a defect channel without departing from the scope of the invention.
After scanning hardware 40 has scanned an image, that image may be enhanced (by altering defects within it) in accordance with the invention using image processing software 42, which is stored in storage medium 38. Alternatively, rather than using software running on a processor (one type of digital circuitry) , the invention may employ other types of digital circuitry comprising any type of dedicated digital hardware to alter defects in the digital image. This hardware may be a part of scanner 34 or general purpose computer 10 as discussed above. Similarly, the scanned image may be stored in storage medium 38 as may the enhanced image. Alternatively, scanner 34 may not have any image processing software 42. Such software instead may be provided on general purpose computer 10 for enhancement of an image received from scanner 34. Enhancement may occur both in scanner 34 and general purpose computer 10 as well. Accordingly, a scanned image and/or an enhanced scanned image may be provided by scanner 34 to the general purpose computer 10 through a communications port (not explicitly shown) . The defect channel may be similarly provided. Although one embodiment of an exemplary scanner 34 that may be used for image enhancement in connection with the invention has been illustrated, other scanners may be used without departing from the scope of the invention.
FIGURE 3 illustrates a flow chart describing a method employed by one embodiment of the present invention to enhance a digital image. The image enhancement described herein may be carried out using computer software, as can any of the processes described below. That software, as discussed above, may be executed by scanner 34, by general purpose computer 10, or a combination thereof. A digital image received from other than scanner 30 may be enhanced in accordance with the invention.
The method described in FIGURE 3 may be used to alter defects in many types of images such as color photographic images (either negative print or transparency) , black and white images (either negative print or transparency and including black and white images derived from photographic film with multiple layers), other monochromatic images, x- rays or any other type of image stored on film. The invention may also be used to alter defects in any image on a tangible medium that may be scanned using a scanner.
In step 50, an image is scanned to create a digital image and defect channel. As noted, however, this step could be omitted and the invention carried out on a image that has been previously scanned and has a defect channel associated with it. In the case of a color image, the digital image will typically be comprised of three channels: a red channel , a green channel , and a blue channel . Each channel is comprised of a series of pixels with each pixel having an intensity value associated with it corresponding to the intensity of the particular color of light at that spatial location in the original image. Other types of color images can be used without departing from the scope of the invention. In addition, it is within the scope of the invention to convert a color image into a black and white image for alteration of defects in the image. The enhanced image with the altered defects could then be converted back into a color image. The methods of the invention could also be applied to a single color channel of a digital image and the same correction could be applied to all channels of the digital image. A black and white image may comprise one or more channels similarly made up of pixels. Other types of images may also comprise one or more similar channels. The invention can be used for any of these types of images.
The defect channel comprises a series of pixels, each having an intensity value either directly proportional or inversely proportional to defects (or the lack thereof) in the original image. Such a defect channel may be created, for example, using an infrared light source and infrared sensor of the kind commonly used in image processing applications involving photographic film. Other types of light sources and sensors may be used to create the defect channel, such as described above in connection with reflective scanners, for example. Any other method may be used to generate a suitable defect channel without departing from the scope of the invention. Such a defect channel will ordinarily produce a signal having a stronger correlation to defects in the original image and a weaker correlation to the visible image itself.
Steps 52-58 comprise a method for altering defects in a digital image to produce an enhanced digital image. Before describing the process in detail, it may be helpful to describe it generally. To enhance a digital image by altering defects in the digital image, the invention corrects defects in at least two frequency bands. The term "correction" or "correct" when used in this application, refers broadly to altering defects in a digital image. An image defect does not need to be completely removed or completely corrected to fall within the scope of the invention. Accordingly, defect correction includes, but is not limited to, reduction or other alteration of a defect in a digital image. The frequency bands for the image are created using one or more filters which average pixels in the neighborhood of a pixel in question with each pixel weighted according to its expected reliability. The expected reliability is determined by the intensity value in the defect channel of the pixel in question or the pixel in question plus a series of pixels in the neighborhood of the pixel in question. Pixels that appear to be more reliable based upon the defect channel are weighted more heavily than those with a low expected degree of reliability. The defect channel may be similarly divided into frequency bands.
The effect of each weighted averaging operation is to create a low pass filtered version of the original image. Where multiple averaging operations are performed, multiple low pass filtered versions of the original image are derived, most often with different bandwidths . These multiple lowpass filtered versions, along with the original image may be used to derive a representation of the original image in multiple contiguous frequency bands. By subtracting one low pass filtered version from another, a bandpass representation for one band may be derived where common frequencies are eliminated. A series of bands can be created similarly, as further described below. The filtered representations of the digital image may then be processed such that the original image is divided into a plurality of frequency bands which collectively make up all or substantially all frequencies present in the original image but overlap little, if at all. These bands may then be recombined to produce an image with defects removed .
The same filtering operation may optionally be performed on the defect channel . Defects may then be further removed, frequency band by frequency band, by subtracting the defect channel bands from the image frequency bands. An enhanced image is then constructed by combining the corrected individual frequency bands back together again.
In step 52, data is optionally converted to log space for the defect channel and for each channel of the digital image where defect correction is to be performed. The conversion of data to logarithmic space may enable the remaining steps of the method to be performed more easily than if the conversion were not made. For example, several of the operations described below involve additions and subtractions in log space, but could involve division by zero in non-log space. Because division by zero can lead to erroneous results, it may be more convenient to carry out the process of the invention in log space. However, the performance of mathematically equivalent functions outside of log space could also be used in connection with the invention. Although any type of logarithmic calculation can be used, one embodiment of the invention uses the conversion to log space described by Formula 1.
γ ( ) = Logιo ( + 1) Logιo 4Q95 )
Log104096
In Formula 1, "x" represents the intensity value of the pixel to be converted to log space. Each such pixel in this embodiment comprises a twelve-bit intensity value. Other numbers of bits can be used without departing from the scope of the invention. Depending upon how the invention is carried out, convenient conversions may be made to take into account the capabilities of digital hardware used to perform calculations in accordance with the invention. If data is converted to log space, and the defect channel is suitable, then a first type of defect correction may be carried out by subtracting the defect channel from each of the visible image channels. Such a subtraction may be a direct subtraction or a bounded subtraction, analogous to the bounded subtraction described below in connection with Figure 5. This subtraction may also take into account red leakage and clear film values where the defect channel is derived using an infrared source and sensor in a film scanner. Any method of taking into account these effects (which are described below) is within the scope of the invention. Defect channels obtained using infrared light in film scanners will often be suitable for this optional enhancement of the image. Even in such a case, however, this step is optional. An alternative to this optional subtraction is to subtract the defect channel from the image channel, frequency band by frequency band, after enhancement by weighted filtering as described below. This alternative is also optional and may be done where a suitable defect channel is available. The defect channel may be suitable where pixel intensity values in the defect channel are proportional, either directly or inversely, to defects and vary approximately linearly with the amount of light blocked by the defect.
In step 54, the image and defect channel are filtered by taking weighted averages of pixels within varying distances from specific pixels. As noted above, each channel of the image may be filtered, the channels may be combined into a single channel and the combined channel filtered or a subset of the channels can be combined or filtered individually and the defect correction applied based upon results obtained from those channels without departing from the scope of the invention. Depending upon the type of enhancement used, the filtering of the defect channel is also optional.
The weighting applied to a specific pixel in calculating the weighted averages may be based upon the expected reliability of that pixel. The expected reliability of a pixel may be determined by using the defect channel. As used herein, "a weighted average" or "averaging" or any similar term refers to any type of average such as a median average, mode average, mean average, arithmetic average, or geometric average. The calculation of weighted averages of pixels around the pixels in each channel of the digital image has the effect of filtering each channel with different strengths of low pass filters. The effect of weighting each pixel (or a subset of the pixels) involved in the calculation of the averages based upon the expected reliability of each pixel dampens the effect of the defect in each filtered version of the original channel of the digital image. An example may illustrate the calculation of the weighted averages. In this example, four filtering operations are performed: a 3x3 weighted average, a 5x5 weighted average, a 9x9 weighted average, and a 17x17 weighted average. In this example, each of these four weighted averages is computed for each pixel in the digital image. Any pixel outside the boundaries of the image is set to zero for purposes of these calculations. For a particular pixel, the 3x3 weighted average is computed using a 3x3 matrix of pixels with the pixel in question at the center of the matrix. Thus, pixels in the neighborhood of the pixel in question are used to calculate the weighted average. The weight applied to each pixel corresponds to its expected reliability as determined using a corresponding pixel or pixels in the defect channel. The weighted average is computed by summing the product of the intensity of each pixel times the weight and dividing the total sum by the sum of the weight values that were applied to each pixel . The 5x5, 9x9, and 17x17 filters are calculated similarly. The calculation of the four discussed weighted averages at each pixel results in four low pass filtered versions of the original channel of the digital image. The 3x3 weighted average applied to each pixel of the image channel produces a low pass filtered version of the image having spatial frequencies between essentially zero and one-third of the maximum spatial frequency possible in the image. Similarly, the 5x5 weighted average produces a low pass filtered version of the original channel of the digital image having frequencies predominantly between zero and one- fifth of the maximum spatial frequency possible in the image. The 9x9 and 17x17 weighted averages produce filtered versions with frequencies predominantly between zero and one-ninth the maximum spatial frequency possible in the image and predominantly between zero and one-seventeenth the maximum spatial frequency possible in the image. Due to the weighting applied to each pixel during this filtering operation, the filtered versions in each frequency band have had the effects of defects reduced. Low pass filtering alone tends to dampen the effects of defects, but may also blur image detail . The weighting may dampen the effects of a defect more than low pass filtering alone would with proportionally less reduction of image detail. This example will be used below to further illustrate adjustment of the image and defect channel . The same filtering operation may also be performed on the defect channel, if desirable for use in further reduction of defects.
Numerous options may be used for calculating the weighted average of a pixel and other pixels in the neighborhood of the pixel . The filter can have a square shape (as in the example above), a circle shape, or any other type of shape without departing from the scope of the invention. The filter may be symmetric or asymmetric and may be centered or not centered on the pixel in question. The filter may also be a window with feathered edges. In addition, any number of filters can be used. The use of more filters will allow division of the image into a greater number of frequency bands. The use of additional filters, however, requires additional calculation. Less filters than those used in the above example can also be used without departing from the scope of the invention.
The filters in the example above had particular dimensions each of which are approximately double the dimensions of the previous filter. Any size filters can be used, however, without departing from the scope of the invention. In the example above, any pixel outside the boundaries of the image was set to an intensity value of zero for purposes of the weighted average calculation. Other boundary conditions could be used without departing from the scope of the invention; for example, soft edge boundary conditions such as a triangle edge or gaussian edge could be used.
In the example above, the weighted average was computed for each pixel in each channel of the digital image and in the defect channel. The weighted average could also be performed on a subset of the pixels in a channel or on all of the pixels in a subset of the channels without departing from the scope of the invention. In addition, different weights could be used for each frequency band and could be omitted in some frequency bands. Different weights could also be used for different channels. The weighting used for a particular pixel can be determined in several ways. The weighting could depend upon the intensity of a spatially corresponding pixel in the defect channel. Where the weighting is applied to the defect itself, the weighting could be determined based upon the intensity of each pixel involved in the calculation. The same or different weightings could be used for the defect channel and one or more image channels. The defect channel, however, may be blurred compared to the visible channel due to focal shifts, the nature of the defects, and inconsistent registration between the visible channels of the digital image and the defect channel, for example. Accordingly, it may be desirable to estimate the reliability of a pixel and determine its weight by taking into consideration a plurality of pixels surrounding a pixel in the defect channel corresponding to the pixel in a visible image channel that is the subject of the weighted average calculation. Alternatively, if blurring effects are insignificant and registration error is fairly constant, then the pixel used in the defect channel could be chosen while compensating for the constant registration error.
The particular weighting function chosen may depend upon the characteristics of the particular scanner used to create the defect channel. The weighting applied may be a function of the intensity value of a pixel or pixels in the defect channel . Where multiple pixels are used to determine a weight, an average or weighted average of the pixels may be used to come up with an average intensity used to determine a weight.
Any type of function may used to relate a weight to the intensity of a pixel in the defect channel or average intensity of a series of pixels in the defect channel. A straight line may be used to establish this function or any other type of curve may be used. In this example, a high threshold and a low threshold may be determined based upon the characteristics of an infrared channel produced by a scanner. For example, in one embodiment, a weight of one may be assigned for all pixel values greater than a particular high threshold where the threshold provides high predictability that no defect is present. A weight of zero may be assigned when the intensity value is below a certain low threshold indicating a high probability that a defect is present with little or no image detail remaining. For points with an intensity in between these threshold values, a straight line or other curve can be used to connect the two thresholds to establish a continuous function for weights between zero and one. Pixels with intensity in this middle region tend to be ones that are defective but some image detail remains to allow enhancement of the digital image to remove the defect but maintain some image detail . In the case of a defect channel derived from a reflective scanner using illumination from multiple angles, the intensity of the pixels in the defect channel may lack a linear relationship to the amount of light that passes through a defect. With such a defect channel, a discrete set of weights such as 0, 1/2 and 1, may be used. An example of a weight function is illustrated in FIGURE 6. In FIGURE 6, infrared intensity in the defect channel will be high when no defect is present as most or all of the infrared energy is allowed to pass through the film. Where a defect is present, however, infrared intensity will be low. In this example, all likely defective pixels (those with an infrared intensity in the defect channel less than 20 percent) are set to a weight of zero. All pixels with a high probability of being nondefective pixels (those with an infrared intensity in the defect channel greater than 64 percent of the maximum intensity are assigned a weight of one) . Pixels in the region in between are assigned a weight based upon the illustrated curve. Again, more complicated functions may be used without departing from the scope of the invention. The weighting applied may also compensate for effects such as leakage from the red channel of a color image into the defect channel . This leakage may result because infrared light sources and/or sensors are residually sensitive to the cyan dye in an image used to modulate the red region of the visible spectrum. This effect manifests itself as an apparent leakage of the red image into the infrared image. The effects of red leakage can be taken into account when establishing a weight value. For example, an overall red leakage value for an image may be calculated. This red leakage value can then be used to establish a constant to be multiplied times a pixel intensity value in the red channel . This product may represent an estimate of the amount of the red channel present in the defect channel for that particular pixel. Thus, in calculating a weight, this product may be subtracted from the intensity of the pixel in the defect channel .
Similarly, a portion of the intensity in the defect channel is proportional to the intensity that would result if infrared light was passed through clear film. Different types of film produce different intensity values when infrared is passed through the clear film. Accordingly, the weight may also be adjusted by subtracting the average clear film value for a particular digital image.
Depending upon the type of filtering employed, one could filter the image once with one or more filters and then filter the corrected image again with different filters. A decision could be made as to whether to apply the second filtering step based upon the estimated size of defects as determined by the defect channel.
In Step 56, the channels of the digital image and the defect channel are adjusted in response to the weighted averages to lessen the effects of defects in the image. The invention includes any use of the weighted averages computed in Step 54 to adjust an image to lessen the effects of defects in the image. A specific method for adjusting the digital image in response to the weighted averages will be discussed below in connection with FIGURE 4. Such adjustment could occur in the time domain or the frequency domain and, when in the spatial domain, in log space or any other space. After the image has been adjusted to lessen the effects of defects, the process terminates in Step 58.
Besides the method described below in connection with Figure 4, another possible method of enhancement is to create a series of essentially contiguous frequency bands based upon the plurality of weighted averages of the original image and the original image itself. The process of creating these bands may be as described below in connection with step 60 of FIGURE 4. After these bands have been created, they may be added together to form an enhanced version of the original image with the defects reduced. This method may be used, for example, where the defect channel is unsuitable for subtraction from the image. This method could, however, be used even where a defect channel is suitable for correction. The enhanced image may also be further enhanced by other operations without departing from the scope of the invention. If the pixel data was converted to log space in step 52 , the enhanced pixel data may be converted back to the space of the original image using a suitable inverse formula such as an inverse to Formula 1.
FIGURE 4 illustrates one example of a method for adjusting images in response to weighted averages to lessen the effects of defects in the digital image. In step 60, the filtered digital image (which comprises the original digital image and the filtered versions of the digital image calculated in Step 54) and the filtered defect channel (which consists of the original defect channel and the filtered versions of the defect channel obtained in step 54) are separated into frequency bands in response to the calculated weighted averages. In general, the filtered images and the original image are used to create a series of frequency bands with little or no overlap representative of the original digital image and original defect channel with defective pixels suppressed due to the weighting that took place during the filtering operation. This method assumes that the defect channel was filtered in FIGURE 3. Using the example discussed above in connection with FIGURE 3 , five frequency bands may be created using the four weighted averages calculated. For purposes of the following, Fm represents the maximum spatial frequency that can exist in the digital image. For a particular channel of the digital image or for the defect channel, the five frequency bands corresponding to the channel in question may be created using the weighted averages computed for that channel. The frequency band from approximately 1/3 Fm to Fm comprises the difference between the original channel and the 3x3 weighted average filtered version of that channel. The frequency band from approximately 1/5 Fm through 1/3 Fm may be calculated by subtracting the 5x5 weighted average version of the channel from the 3x3 weighted average version of the original channel . The frequency bands from approximately 1/5 Fm to 1/9 Fm and approximately 1/17 Fm to 1/9 Fm may be determined by subtracting the 9x9 weighted average version of the original channel from the 5x5 version of the original image channel and by subtracting the 17x17 weighted average version of the original channel from the 9x9 weighted average version of the original channel, respectively. Finally, the frequency band from approximately zero to 1/17 Fm is represented by the 17x17 weighted average version of the original image channel.
In general, the low pass filtered versions of the image and defect channel created using the weighted averages can be used to divide each channel of the image as well as the defect channel into contiguous frequency bands where some or all of the frequency bands have had the effect of defects suppressed using the weighted average calculations.
In Step 62, the red residue in each frequency band of the defect channel may be removed. Optionally, such residue may be removed using a bounded calculation to allow some variance. If infrared sources and sensors (or other sources and sensors) are used that do not leave red residue in the defect channel, then the step may be omitted without departing from the scope of the invention. One option for removing the red leakage from the defect channel is to subtract from the intensity value of each pixel in the defect channel, the product of the intensity of the corresponding pixel in the red channel multiplied by an average red leakage constant representing the average red leakage for the entire digital image. This difference may be divided by the difference between one and the red leakage value to properly normalize the result. Alternatively, because red leakage may vary within regions of particular images, it may be useful to use a bounded subtraction to allow for some variance in localized red leakage values.
An example of a bounded subtraction for red residue will be provided in connection with FIGURE 5 below. In general, two or more subtractions are performed from the intensity value in the defect channel. The product of the intensity value for the pixel in question in the red channel is multiplied by the red leakage constant adjusted upward for one subtraction and adjusted downward for a second subtraction from the intensity value in the defect channel. These results are divided by the difference between one and the appropriate adjusted red leakage constant. If both calculations have results having the same sign (i.e., both are positive or both are negative) then the smaller result is chosen. If the results have opposite signs, then the pixel in question is set to zero. Besides allowing variance in the red leakage value, such a bounded subtraction may compensate for registration error between visible channels of the digital image and the defect channel as well as for blurring that may occur in the defect channel . In Step 64, defects are further removed from the image by subtracting the relevant frequency band of the defect channel obtained in Step 62 from the relevant frequency band in each image channel obtained in Step 60. Optionally, a bounded calculation such as that described in connection with Step 62 may be used to allow for some variance caused by registration error, blurring, etc. in the defect channel. Next, in Step 66, the frequency bands are recombined through a summation to form an enhanced image with the original defect removed. In Step 67, the enhanced image is converted back from log space to the original space from which it was derived using a suitable inverse formula such as an inverse to Formula 1. If the enhanced image was created in a space other than log space, then Step 67 may be omitted without departing from the scope of the invention. If Step 60 through 66 were carried out in the frequency domain, then the enhanced image may be reconverted back to the time domain in Step 67. FIGURE 5 illustrates a flow chart describing a second method of enhancing an image by altering defects in the image in accordance with the invention. In this embodiment, an image may be divided into segments and defects processed within each individual segment . By dividing the image into segments, the amount of storage space used at any one time to enhance the image may be reduced, and the amount of computation required may be reduced. Thus, in this embodiment, the image is divided into segments and a process similar to that discussed above in connection with FIGURES 3 and 4 is applied to each segment individually as if that segment comprised the entire image.
In step 68 an image is scanned to create a digital image comprising one or more channels and a defect channel. All of the options discussed above in connection with step 50 are available for this embodiment of the invention.
Next, in step 70 the data from the defect channel and the channels of the digital image is converted to log space. Again, any of the options discussed above in connection with step 52 may be employed in step 70, including optional subtraction of the defect channel from the visible channel.
In step 72 , each channel of the digital image and the defect channel may be divided into segments. Any size or shape of segments may be used without departing from the scope of the invention. In one embodiment, each channel of the digital image and the defect channel are divided into 8x8 segments. This example will be used to illustrate the remainder of the steps of the method illustrated in FIGURE 5.
In step 74, weighted averages of pixels within varying distances from specific pixels in a segment of the defect channel and a segment of each channel of the image are computed. These weighted averages are computed similarly to the weighted averages that were computed in step 54. However, the weighted averages computed in step 74 are computed for an individual segment of the image, assuming that the segment comprises the entire image. Accordingly, even if the 8x8 segment is surrounded by other 8x8 segments, the pixels beyond the boundaries of 8x8 segment in question are treated as beyond the image boundary and any pixel beyond those boundaries are treated as having an intensity of zero (or other boundary condition) .
All of the options for filtering and weighting discussed above in connection with step 54 may be employed in step 74 for the embodiment disclosed in FIGURE 5
(including the option of not filtering the defect channel) .
However, because the segment in this example is 8x8, an 8x8 filter is the maximum sized filter that will be used in this example. The weight function, in this embodiment, takes into account red leakage and clear film effects as discussed above in connection with step 54. The weight for a particular pixel for the weighted average calculations may be determined, for example, using Formula 2.
W(x,y)=2(Dιn(x,y)-RLRιn(x,y)-CF)+1.4 (2)
Limit 0<W(x,y)<l
In Formula 2, x and y represent the coordinates of a particular pixel within an 8x8 segment. Din represents the defect channel received from the scanner in step 68. RL comprises a constant value representative of average red leakage in either the segment in question or in the entire digital image. RΪN represents the red channel of the digital image. CF is a constant representing the average clear film value for either the segment in question or the entire digital image. The weight function of Formula 2 is constrained such that the weight varies between zero and one. A plot of the weighting function versus the intensity in the defect channel would look similar to the curve illustrated in FIGURE 6. Any other weight function could be used without departing from the scope of the invention.
In this example, two different weighted averages are computed for each channel . A 3x3 weighted average is computed. In addition, an 8x8 weighted average is computed. The 8x8 average covers the entire segment, and as a result, a single scalar value may be used to represent the result of this weighted average. The 3x3 weighted average for the red, green, blue, and defect channels may be computed using Formulas 3 through 6. ∑ ∑ W(x + a,y+b)Rm(x + a,y+b) , ,
R3L(x,y)=a-l b°-' , j where ∑ ∑ W(x + a, y+b)≠0
∑ ∑ W(x + a,y+b) -' "-' a=-l b=-l
(3)
R3L(x>y) = ° where ∑ ∑ W(x + a, y4-b)=0 a=-l b=-l
∑ ∑ W(x + a,y+b)Gιn(x + a,y+b)
G3L(x>y)=,-l b^1 , ; where∑ ∑ W(x + a, y+b)≠0
∑ ∑ W vv(x Λ + T a,,y y+ Tb υ) a=-l b=-l i a=-l b=-l
(4)
1 1
G3L(x,y) = 0 where∑ ∑ W(x + a, y+b)=0
-1 b=-l
∑ ∑ W(x + a,y+b)Bin(x + a, y+b)
B3L(x,y)='=-l b-1 , j where∑ ∑ W(x + a, y+b)≠0
∑ ∑ W vi( iΛxτ + aα,,y y+ -rb u) a=-l b=-l ;
=-1 b=-l
(5) B3L(x,y) = 0 where ∑ ∑ W(x + a, y+b)=0 a=-l b=-l
∑ ∑ W(x + a,y+b)Dιn(x- a,y+b) , ,
D3L(x,y) = '-l b-1 , ; where ∑ ∑ W(x + a, y+b)≠0
∑ ∑ W rv(x Λ + -r a a,,y y+ ~rb v) a=-l b=-l l l b=-ι
(6)
1 1
O3h(x,y)= where∑ ∑ W(x + a, y+b)=0 a=-l b=-l
In these Formulas, Rιn, Gin, Bin, Dιn, respectively, represent the red, green, blue, and defect channels of the digital image that was scanned in step 68. The 8x8 weighted average for the red, green, blue, and defect channels may be calculated using Formulas 7 through 10, respectively. As noted, each of these Formulas produces a scalar result.
∑ ∑ W(a,b)Rm(a,b) 7 7
R8L(x,y)=a° 7 where ∑ ∑ W(a,b)≠0
∑ ∑ W(a,b) a=0 b=0 a=0 b=0
(7)
7 7
R8L (x, y)= 0 where ∑ ∑ W (a, b)= 0 a=0 b=0
∑ ∑ W(a,b)Gm(a,b) 7 7
G8L(x,y)=a=0 b 7 =° 7 where ∑ ∑ w(a,b)≠0
∑ ∑ W(a,b) a=0 b=0 a=0 b=0 (8)
7 7
G8L(x,y)=0 where ∑ ∑ W(a,b)=0 a=0 b=0
7 7
∑ ∑ W(a,b)Bιn(a,b) 7 7
B8L(x,y)=a° b° 7 where ∑ ∑ W(a,b)≠0
∑ ∑ W(a,b) a=0 b=0 a=0 b=0
(9)
B8L (x,y)=0 where ∑ ∑ w(a,b)=0 a=0 b=0
Figure imgf000028_0001
D8L (x, y)=0 where ∑ ∑ W(a,b)=0 a=0 b=0
R3L represents a weighted low-pass filtered version of the red channel of the digital image having a frequency band between approximately zero and 1/3 Fm, where Fm is the maximum spatial frequency possible for the segment of the digital image. G3L, B3L, and D3L cover a similar frequency band for the green, blue, and defect channels respectively. R8L represents a weighted low-pass filtered version of the red channel of the digital image having a frequency band between approximately zero and 1/8 Fm. G8L, B8L, and D8L represent weighted low-pass filtered versions of the green, blue, and defect channels, respectively in the same frequency band.
In step 76, the segment of the image is adjusted to lessen the effects of defects in the image in response to the weighted averages of each image channel and the defect channel obtained in step 74 and the original image channels and defect channels. Any of the options discussed above in connection with step 56 and steps 60 through 67 may be used in connection such adjustment.
In this example, a version of the original image and original defect channel, with effects of the defect partially suppressed due the weighted averages computed in step 74, is created in three separate frequency bands. For each image channel and the defect channel, a first frequency band is obtained having frequencies between approximately zero and 1/8 Fm. A second frequency band extends from approximately 1/8 Fm to 1/3 Fm. A third frequency band extends from approximately 1/3 Fm to Fm.
The first frequency band may be obtained using Formulas 11 through 14.
^ 8B ^ 8L (11)
^ δB ^δL (12)
Figure imgf000029_0001
D8B =D8L (14)
Because this frequency band was already calculated in step 74, it can be used without any further computation.
The third frequency may be computed using Formulas 15 through 18. In these Formulas, the weighted low-pass filtered version of each channel in the frequency range of approximately zero to 1/3 Fm is subtracted from the unfiltered image channel to produce the frequency band between approximately 1/3 Fm and Fra.
RIB (χ,y)=Rιn (χ,y)-R3L (χ,y) (15)
G]B(x,y)=Gm(x,y)-G3L(x,y) (16)
B1B(x,y)=Bm(x,y)-B3L(x,y) (17)
D1B(x,y)=Dιn(x,y)-D3L(x,y) (18) Finally, the second frequency band, the one between approximately 1/8 Fm and 1/3 Fm may be computed using Formulas 19 through 22. In Formulas 19 through 22, the 8x8 weighted filtered version of each image channel is subtracted from the 3x3 weighted filtered version of the original image channel. The process just described for creating a plurality of frequency bands may employ any of the options discussed above in connection with step 60.
R3B(x,y)=R3L(x,y)-R8L (19)
G3B(x,y)=G3L(x,y)-G8L (20)
B3B(x,y)=B3L(x,y)-B8L (21)
D3B(x,y)=D3L(x,y)-D8L (22)
Next, the red residue may be removed from the defect channel in each frequency band. Again, any of the options described above in connection with step 62 may be used without departing from the scope of the invention. In this embodiment, a bounded subtraction is used for two of the frequency band and a simple subtraction is used for the remaining frequency band. Formula 23 may be used to subtract the red residue in the first frequency band of the defect channel .
D8B = D8B~R8BDK8RL (23) ι-DK8RL
In Formula 23, the constant, DK8, will normally be set to one, but may vary depending upon the scanner in question.
This constant represents the amount of the red channel that should be subtracted from the defect relative to the measured average red leakage for the entire segment and/or the entire image.
Formulas 24 through 26 may be used to subtract the red residue from the second frequency band of the defect channel .
Figure imgf000031_0001
T2 D3B (χ,y)-R3B (χ,y)Dκ3L R L
D3B (25)
I-D^RL
D 3B (x >y) = T1D3B (χ.y) where |T1D3B (x,y) < |T2D3B (x,y) andTlD3B (x,y)andT2D3B (x,y)have the same sign
D 3B(χ 5y) = T2 D3B(χ ;y) where |T2D3B (x,y) < |T1D3B (x,y) (26) andTlD3B (x,y)andT2D3B (x,y)have the same sign
D3B(x,y) = 0 where T1D3B (x,y)andT2D3B (x,y)have different signs
The bounded subtraction provided for in Formulas 24 through 26 allows for local variance in the red leakage value, as well as for some registration error and/or blurring in the defect channel as compared to the visible channel . Formula 24 multiplies a high value constant, DK3H against the red leakage constant, RL. Formula 25 multiplies a low-value constant, D 3L, by the red leakage constant, RL. DK3H and DK3L are the high and low ranges determining how much of the image content in the red channel should be subtracted from the defect channel . These constants may be determined experimentally. In this embodiment, D 3H s chosen to be 1.3 while DK3L is chosen to be 0.6. These constants will tend to average about one. The average is a function of the resolution of the system for the infrared channel versus the red channel . The spread between these two constants is a function of the accuracy of the scanner that produces the defect channel . If DK3H is chosen too large, then small image patterns may cause large matching defects to be erased, rendering defects uncorrected in middle frequencies when they are next to image detail. If DK3H is chosen too small, then image residue may remain in the defect record, causing middle frequency image detail to be erased. If DK3L is chosen too large, then too much image detail may be removed from the defect causing a negative residue that interferes with the ability to distinguish a defect. If DK3L is chosen too small, then visible detail may not be removed in the presence of defect detail of opposite polarity.
After the results of Formula 24 and 25 are calculated, the revised version of the frequency band is calculated using Formula 26. If the results of Formula 24 and Formula 25 have different signs, then the intensity value is set to zero. If these signs are the same, then the lesser value produced by either Formula 24 or Formula 25 is chosen. Thus, this bounded subtraction tends to drive the revised value closer to zero.
Formulas 27 through 32 may be used to remove red residue in the defect channel for the third frequency band created above .
Figure imgf000032_0001
Figure imgf000032_0002
T3D1B(x,y) = TlD1B(x,y) where |T1D1B (x,y) ≤ |T2D1B (x,y) andTlD1B (x,y)andT2D]B (x,y)have the same sign
T3D1B(x,y) = T2D1B(x,y) where |T2D1B (x,y) < |T1D1B (x,y) (29) andT2D1B (x,y)andTlD1B (x,y)have the same sign
T3D1B(x,y) = 0 where T1D1B (x,y)andT2D1B (x,y)havedifferentsigns
T4D1B(x,y) = T3D1B(x,y)-DK1A (30)
τ5DiB( >y) = τ3D1B(χ 5y)-DκlA (31)
DlB( ,y) = T4D1B(x,y) where |T4D1B (x,y) < |T5D1B [x,y] andT4D1B (x,y)andT5D1B (x,y) have the same sign
D'1B (x, y) = T5D1B (x, y) where |T5D1B (x, yj < |T4D1B (x, y) (32) andT5D1B (x,y)andT4D1B (x,y)have the same sign
D'1B(x,y) = 0 where T5D1B (x,y)andT4D1B (x,y) have opposite signs
Formulas 27 through 29 are similar to Formulas 24 through 26. For this frequency band, in this example, the high- constant DKiH is chosen to be 1.0 while the low- frequency constant, DKιL is chosen to be 0.3. Similar considerations apply to the choice of these constants, as applied to the choice of DK3H and DK3L above. In this case, the effects of these constants will be in the high frequency band. Here, the average of the two constants is less than one and may be chosen as such to the extent that heightened frequency defects are blurred relative to the high frequency image. The range is also wider to accommodate a greater variance that often occurs at higher frequencies. The difference between Formulas 27 through 32 and Formulas 23 through 26 involves the addition of Formulas 30 through 31. These Formulas make a minor adjustment in the value computed by Formula 29 to adjust for residual noise in the high-frequency portion of the defect channel. In this embodiment, constant DK1A has a value of 0.01. If the residual noise constant, OK1 is set too high, then detail in small defects in the image may not be removed. If the value is set too low, then too little residual electronic noise will be removed from this frequency band of the defect channel and such noise in the defect channel may appear in the visible channel (s) as a negative of this noise, causing the electronic noise in the defect channel to contaminate the processed visible image. Next, the defect is removed from the visible image frequency band by frequency band. Similar bounded subtractions are performed as were performed when the red residue was removed from the defect channel above. In this case, the constant for the first frequency band, RK8 is set to one, but could be adjusted. For this frequency band, a direct subtraction is used, rather than a bounded subtraction. Blurring effects tend not to affect this frequency band or affect it insignificantly, but a bounded subtraction could be employed if desired. The constants RK3H and RK3L are set to 1.5 and 0.5 respectively. The constants RKIH and RKIL are set to 1.7 and 0.5, respectively. These values are also chosen for the corresponding constants for the green and blue channels. These constants establish the bounds within which correction is made. Again, one selects whichever number between the bounds drives the corrected result closet to zero. The upper and lower bounds will normally average about one, except that one may set the lower bound a bit lower without damaging the image. Finally, the defects in the original image may be suppressed even further by multiplying the pixels resulting from a bounded subtraction of the defect channel from the image channel by the weight for a particular pixel calculated using Formula 2.
The defects may be removed from each frequency band of the visible image channels, using Formulas 33 through 59. Formulas 33 through 35 may be used to remove defect information from the first frequency band. Formulas 36 through 47 may be used to remove the defect from the middle frequency band of each channel. Formulas 48 through 59 may be used to remove the defect from the third frequency band of each channel .
R'8B =R8B -D'8BRK8 (33)
G8B =G8B -D8B GK8 (34)
BSB =B8B-D'8BBK8 (35)
TIMB (X, y) = R3B (χ> y)-D3B (χ> y) K3H (36)
T2R3B(χ,y)=R3B (χ,y)-D3B (χ,y)RK3L (37)
T R3B( >y) = TlR3B(x,y) where |T1R3B (x,y) < |T2R3B (x,y) andTlR3B (x,y)andT2R3B (x,y)have the same sign
τ3 R3B(χ,y) = T2R3B(x,y) where |T2R3B (x,y) < |T1R3B (x,y) (38) andT2R3B (x,y)andTlR3B (x, y) have the same sign
T3R3B(x,y) = 0 where T2R3B (x,y)andTlR3B (x,y) have opposite signs
R'3B(x,y)=T3R3B(x,y)w(x,y) (39)
TlG3B(x,y)=G3B(x,y)-D;B(x,y)GK3H (40) T2G3B(x,y)=G3B (x,y)-D'3B (x.yjG^ (41)
T3G3B(x,y) = TlG3B(x,y) where |T1G3B (x,yj < |T2G3B (x,y) andTlG3B (x,y)andT2G3B (x,y) have the same sign
T3G3B(x,y) = T2G3B(x,y) where |T2G3B (x,y) < |T1G3B (x,y) (42) andT2G3B (x,y)andTlG3B (x,y) have the same sign
T3GG33BB(x,y) = 0 where T2G3B(x,y)andTlG3B(x,y)haveoppositesigns
G3B(x,y)=T3G3B(x,y)w(x,y) (43)
TlB3B(x,y) = B3B(x,y)-D3B (x^B^ (44)
T2B3B(x,y)=B3B (x,y)-D'3B (x,y)BK3L (45)
T3 B3B(x>y) = TlB3B(x,y) where |T1B3B (x,y) < |T2B3B (x,yj and T 1 B3B (x, y) and T2 B3B (x, y) have the same sign
T3B3B(x,y) = T2B3B(x,y) where |T2B3B (x,yj < |T1B3B (x,y) (46) andT2B3B (x,y)andTlB3B (x,y) have the same sign
T3B3B(x,y) = 0 where T2B3B (x,y)andTlB3B(x,y)haveoppositesigns
B3B(x,y)=T3B3B(x,y)w(x,y) (47)
T1R1B (x,y)=R1B (x,y)-D'1B (x,y) RK1H (48)
T2R1B (x,y)=R1B (x,y)-D;B (x,y) RK1L (49) T3R1B (x,y) = T1R1B (x,y) where |T1R1B (x,y) ≤ |T2R1B (x,yj andTlR1B (x,y)andT2R1B (x,y)have the same sign
T3R1B (x,y) = T2R1B (x,y) where |T2R1B (x,y) < |T1R1B (x,y) (50) andT2R1B (x,y)andTlR1B (x,y) have the same sign
T3R1B(x,y) = 0 where T2R1B (x,y)andTlR1B (x,y) have opposite signs
R;B=T3R1B(x,y)w(x,y) (51)
T1G1B (x,y)=G1B (x,y)-D;B (x,y) GK1H (52)
T2G1B (χ,y)=G1B (χ,y)-D'1B (x,y) GK1L (53)
T3 GiB(x >y) =
Figure imgf000037_0001
where |T1G1B (x,y)| < |T2G1B (x,yj andTlG1B (x,y)andT2G1B (x,y)have the same sign
T3GIB (χ,y) = T GiB(x>y) where |T2G1B (x,yj < |T1G1B (x,yj (54) andT2G1B (x,y)andTlG1B (x, y) have the same sign
T3 GiB(x >y) = ° where T2G1B (x,y)andTlG1B (x,y) have opposite sign
G;B=T3G1B(x,y)w(x,y) (55)
TIEIE (x >y)=B1B (x,y)-D;B (x,y) BK1H (56)
T2B1B (x,y)=B1B (x,y)-D;B (x,y) BK1L (57) T3B1B (x,y) = TlB1B (x,y) where |T1B1B (x, y) ≤ |T2B1B (x,y) andTlB1B (x,y)andT2B1B (x, y) have the same sign
T3B1B (x, y) = T2B1B (x, y) where |T2B1B (x, yj < |T1B1B (x,y) (58) andT2B1B (x, y)andTlB1B (x,y)have the same sign
T3B1B (x,y) = 0 where T2B1B (x,y)andTlB1B (x,y) have opposite signs
B'1B =T3B1B (x, y) w(x,y) (59)
To obtain an enhanced image with the defect removed, the frequency bands may be recombined for each channel using Formulas 60 through 62.
Rout (χ,y)=R'1B (χ,y)+R3B (x,y)+R'8B ( >y) (60)
Gout (x,y)=G;B (x, y)+G'3B (x, y)+G'8B (x, y) (61)
Bout (x,y)=B'1B (x, y)+B'3B (x, y)+B'8B (x, y) (62)
These enhanced images may be converted back to the original image space using an appropriate inverse logarithmic function. Again, all of the options discussed above in connection with step 67 are applicable to this embodiment as well .
In step 78, it is determined whether there are any more segments of the digital image to process. If so, then steps
74 and 76 are repeated for each remaining segment of the digital image. If no more segments are to be processed, then the method continues in step 80.
In step 80, it is determined if there are more overlaps to perform. If not, then the procedure terminates in step
84. If so, then each channel of the enhanced digital image obtained by using steps 74 and 76 for each segment is divided into segments spatially different from the earlier segments. Steps 74 and 76 are then repeated for each of these segments. This aspect of this embodiment may compensate for results that may be obtained for pixels near the boundary of segments that were used when the image was first divided into segments in step 72. By performing defect removal a second time, the defect may be further suppressed. In this example, a second overlap is performed using 8x8 segments comprising a 2x2 corner of four adjoining segments that touch one another at a common point . These 8x8 segments are thus made up of one- fourth of each of four adjoining segments. Any of the options described above with respect to step 72 may be used in step 82 in performing division of the enhanced image into different segments. Different shapes can be used and different size segments can be used without departing from the scope of the invention.
Although the inventions described herein involve calculations in the spatial domain, analogous calculations in the frequency domain could equivalently be used without departing from the scope of the invention.
Although the present invention has been described in detail, it should be understood that various changes, substitutions and alterations can be made hereto without departing from the sphere and scope of the invention as defined by the appended claims.
To aid the Patent Office, and any readers of any patent issued on this application in interpreting the claims appended hereto, applicants wish to note that they do not intend any of the appended claims to invoke f 6 of 35 U.S.C.
§ 112 as it exists on the date of filing hereof unless "means for" or "step for" are used in the particular claim.

Claims

WHAT IS CLAIMED IS:
1. A method for altering defects in a digital image, comprising: storing on a computer readable storage medium a digital image comprising at least one channel, the at least one channel comprising a plurality of pixels each having an intensity value; storing on the computer readable storage medium a defect channel comprising a plurality of pixels each having an intensity value, at least some of the plurality of pixels in the defect channel having an intensity value proportional to defects in the digital image; filtering a first pixel of the at least one channel of the digital image to produce a filtered pixel the filtering of the first pixel responsive to the intensity of the first pixel and a plurality of additional pixels in the neighborhood of the first pixel, wherein at least one of the pixels used in filtering the first pixel is weighted in response to the intensity value of at least one pixel in the defect channel; and producing a corrected digital image in response to the filtered pixel and the at least one channel.
2. The method of Claim 1 : wherein the digital image comprises three channels, each channel comprising a plurality of pixels; wherein the digital image represents a color image; and wherein the method further comprises filtering at least a second pixel at the same location in each of the three channels to produce a second filtered pixel for each of the three channels, the filtering of the second pixel in each channel, for a particular channel, responsive to the intensity of the second pixel in the particular channel and a plurality of additional pixels in the neighborhood of the second pixel in the particular channel, wherein at least one of the pixels used in filtering the second pixel in the particular channel is weighted in response to the intensity value of at least one pixel in the defect channel; producing a corrected digital image in response to the filtered pixel, the second filtered pixel for each channel, and the three channels.
3. The method of Claim 1, wherein the digital image comprises a single channel and represents a monochrome image .
4. The method of Claim 1, wherein at least one of the pixels used in filtering the first pixel is weighted in response to the intensity value of a plurality of pixels in the defect channel .
5. The method of Claim 1, wherein at least one of the pixels used in filtering the first pixel is weighted in response to the intensity value of a pixel in the defect channel which spatially corresponds to the first pixel.
6. The method of Claim 1, wherein the weighting of the at least one of the pixels used in filtering the first pixel is further responsive to red leakage in the defect channel .
7. The method of Claim 1, wherein the weighting of the at least one of the pixels used in filtering the first pixel is further responsive to a clear film intensity measure .
8. The method of Claim 1, wherein the filtering is performed in log space .
9. The method of Claim 8, wherein the corrected digital image is further responsive to the difference between the digital image and a filtered version of the digital image containing the filtered pixel.
10. The method of Claim 8, wherein the corrected digital image is further responsive to the difference between (a) the difference between the digital image and a filtered version of the digital image containing the filtered pixel and (b) the difference between the original defect channel and a filtered version of the defect channel .
11. The method of Claim 1, wherein the intensity of pixels in the defect channel is responsive to the difference between at least two reflected scans of an image from which the digital image was derived, wherein the image is illuminated by at least one light source at a different angle relative to the surface of the image during each scan.
12. The method of Claim 1, wherein filtering the first pixel further comprises: filtering the first pixel a plurality of times, wherein ones of the plurality of filtering operations cover a different bandwidth, are responsive to the intensity of the first pixel and a plurality of additional pixels in the neighborhood of the first pixel, and wherein at least one of the pixels use in filtering the first pixel is weighted in response to the intensity value of at least one pixel in the defect channel, each of the plurality of filtering operations producing a filtered pixel, collectively comprising a plurality of filtered pixels; and wherein the corrected digital image is further responsive to the plurality of filtered pixels.
13. The method of Claim 1, wherein the digital image comprises three channels, each channel comprising a plurality of pixels; wherein the digital image represents a color image; wherein the at least one channel comprises one of the three channels; and wherein the corrected digital image is produced in response to the filtered pixel by applying a correction in all three channels in response to the filtered pixel determined for the at least one channel .
14. The method of Claim 1, wherein the digital image comprises three channels, each channel comprising a plurality of pixels; wherein the digital image represents a color image; wherein the method further comprises forming a black and white channel based upon the three channels, the black and white channel comprising the at least one channel; and wherein the corrected digital image is produced in response to the filtered pixel by applying a correction in all three channels in response to the filtered pixel determined for the black and white channel .
15. The method of Claim 1, wherein the intensity of pixels in the defect channel is proportional to infrared transmission through film from which the digital image was derived.
16. The method of Claim 15, further comprising: assigning a low weight value to a pixel for the filtering operation if a spatially corresponding pixel in the defect channel has an intensity less than a first threshold.
17. The method of Claim 15, further comprising: assigning a low weight value to a pixel for the filtering operation if the average intensity of a spatially corresponding pixel in the defect channel and other pixels in the defect channel in the neighborhood of the spatially corresponding pixel is less than a first threshold.
18. The method of Claim 1, further comprising: dividing the digital image into a plurality of segments ; wherein the first pixel is in a first segment comprising one of the plurality of segments of the at least one channel of the digital image and wherein the plurality of additional pixels in the neighborhood of the first pixel are in the first segment; filtering a second pixel of the at least one channel of the digital image to produce a second filtered pixel the filtering of the second pixel responsive to the intensity of the second pixel and a plurality of additional pixels in the neighborhood of the second pixel, wherein at least one of the pixels used in filtering the second pixel is weighted in response to the intensity value of at least one pixel in the defect channel; wherein the second pixel is in a second segment comprising one of the plurality of segments of the at least one channel of the digital image and wherein the plurality of additional pixels in the neighborhood of the second pixel are in the second segment ; and wherein the corrected digital image is produced in response the filtered pixel, the second filtered pixel, and the at least one channel .
19. The method of Claim 1, wherein the first pixel is filtered using a square filter.
20. The method of Claim 1, wherein the first pixel is filtered using a circular filter.
21. A method for altering defects in a digital image, comprising : filtering, using digital circuitry, a first pixel of a first channel of the digital image to produce a filtered pixel by averaging the intensity of the first pixel and a plurality of additional pixels in the neighborhood of the first pixel, wherein at least one of the pixels used for averaging is weighted in response to the intensity value of at least one pixel in a defect channel associated with the digital image; and producing a corrected digital image, using digital circuitry, in response to the filtered pixel and the first channel .
22. The method of Claim 21, wherein the digital circuitry comprises dedicated digital hardware.
23. The method of Claim 21, wherein the digital circuitry comprises a microprocessor executing computer software.
24. The method of Claim 21, wherein the digital image comprises three channels, each channel comprising a plurality of pixels; wherein the digital image represents a color image; and wherein the method further comprises filtering at least a second pixel at the same location in each of the three channels to produce a second filtered pixel for each of the three channels, the filtering of the second pixel in each channel, for a particular channel, responsive to the intensity of the second pixel in the particular channel and a plurality of additional pixels in the neighborhood of the second pixel in the particular channel, wherein at least one of the pixels used in filtering the second pixel in the particular channel is weighted in response to the intensity value of at least one pixel in the defect channel; producing a corrected digital image in response to the filtered pixel, the second filtered pixel for each channel, and the three channels.
25. The method of Claim 21, wherein at least one of the pixels used in filtering the first pixel is weighted in response to the intensity value of a plurality of pixels in the defect channel.
26. The method of Claim 21, wherein at least one of the pixels used in filtering the first pixel is weighted in response to the intensity value of a pixel in the defect channel which spatially corresponds to the first pixel.
27. The method of Claim 21, further comprising: filtering the first pixel a plurality of times, wherein ones of the plurality of filtering operations cover a different bandwidth, are responsive to the intensity of the first pixel and a plurality of additional pixels in the neighborhood of the first pixel, and wherein at least one of the pixels use in filtering the first pixel is weighted in response to the intensity value of at least one pixel in the defect channel, each of the plurality of filtering operations producing a filtered pixel, collectively comprising a plurality of filtered pixels; and wherein the corrected digital image is further responsive to the plurality of filtered pixels.
28. The method of Claim 21, further comprising: assigning a low weight value to a pixel for the filtering operation if a spatially corresponding pixel in the defect channel has an intensity indicating a high probability that the pixel is defective.
29. The method of Claim 21, further comprising: assigning a low weight value to a pixel for the filtering operation if the average intensity of a spatially corresponding pixel in the defect channel and other pixels in the defect channel in the neighborhood of the spatially corresponding pixel indicates a high probability that the pixel is defective.
30. The method of Claim 21, further comprising: performing a filtering operation on each pixel in the at least one channel .
31. The method of Claim 21, further comprising: dividing the digital image into a plurality of segments ; wherein the first pixel is in a first segment comprising one of the plurality of segments of the at least one channel of the digital image and wherein the plurality of additional pixels in the neighborhood of the first pixel are in the first segment; filtering a second pixel of the at least one channel of the digital image to produce a second filtered pixel the filtering of the second pixel responsive to the intensity of the second pixel and a plurality of additional pixels in the neighborhood of the second pixel, wherein at least one of the pixels used in filtering the second pixel is weighted in response to the intensity value of at least one pixel in the defect channel; wherein the second pixel is in a second segment comprising one of the plurality of segments of the at least one channel of the digital image and wherein the plurality of additional pixels in the neighborhood of the second pixel are in the second segment ; and wherein the corrected digital image is produced in response the filtered pixel, the second filtered pixel, and the at least one channel .
32. A digital image scanning system comprising: scanning hardware operable to scan an image and convert the image into a digital image having at least one channel and to produce a defect channel responsive to defects in the image ; and computer software associated with the scanning hardware and operable to: filter a first pixel of the at least one channel of the digital image to produce a filtered pixel the filtering of the first pixel responsive to the intensity of the first pixel and a plurality of additional pixels in the neighborhood of the first pixel, wherein at least one of the pixels used in filtering the first pixel is weighted in response to the intensity value of at least one pixel in the defect channel; and produce a corrected digital image in response to the filtered pixel and the at least one channel.
33. The digital image scanning system of Claim 32, wherein the digital image comprises three channels, each channel comprising a plurality of pixels; wherein the digital image represents a color image; and wherein the computer software is further operable to filter at least a second pixel at the same location in each of the three channels to produce a second filtered pixel for each of the three channels, the filtering of the second pixel in each channel, for a particular channel, responsive to the intensity of the second pixel in the particular channel and a plurality of additional pixels in the neighborhood of the second pixel in the particular channel, wherein at least one of the pixels used in filtering the second pixel in the particular channel is weighted in response to the intensity value of at least one pixel in the defect channel; produce a corrected digital image in response to the filtered pixel, the second filtered pixel for each channel, and the three channels.
34. The digital image scanning system of Claim 32, wherein at least one of the pixels used in filtering the first pixel is weighted in response to the intensity value of a plurality of pixels in the defect channel.
35. The digital image scanning system of Claim 32, wherein at least one of the pixels used in filtering the first pixel is weighted in response to the intensity value of a pixel in the defect channel which spatially corresponds to the first pixel .
36. The digital image scanning system of Claim 32, wherein the computer software is further operable to: filter the first pixel a plurality of times, wherein ones of the plurality of filtering operations cover a different bandwidth, are responsive to the intensity of the first pixel and a plurality of additional pixels in the neighborhood of the first pixel, and wherein at least one of the pixels use in filtering the first pixel is weighted in response to the intensity value of at least one pixel in the defect channel, each of the plurality of filtering operations producing a filtered pixel, collectively comprising a plurality of filtered pixels; and wherein the corrected digital image is further responsive to the plurality of filtered pixels.
37. The digital image scanning system of Claim 32, wherein the computer software is further operable to: assign a low weight value to a pixel for the filtering operation if a spatially corresponding pixel in the defect channel has an intensity indicating a high probability that the pixel is defective.
38. The digital image scanning system of Claim 32, wherein the computer software is further operable to: assign a low weight value to a pixel for the filtering operation if the average intensity of a spatially corresponding pixel in the defect channel and other pixels in the defect channel in the neighborhood of the spatially corresponding pixel indicates a high probability that the pixel is defective.
39. The digital image scanning system of Claim 32, wherein the computer software is further operable to: divide the digital image into a plurality of segments; wherein the first pixel is in a first segment comprising one of the plurality of segments of the at least one channel of the digital image and wherein the plurality of additional pixels in the neighborhood of the first pixel are in the first segment; filter a second pixel of the at least one channel of the digital image to produce a second filtered pixel the filtering of the second pixel responsive to the intensity of the second pixel and a plurality of additional pixels in the neighborhood of the second pixel, wherein at least one of the pixels used in filtering the second pixel is weighted in response to the intensity value of at least one pixel in the defect channel; wherein the second pixel is in a second segment comprising one of the plurality of segments of the at least one channel of the digital image and wherein the plurality of additional pixels in the neighborhood of the second pixel are in the second segment ; and wherein the corrected digital image is produced in response the filtered pixel, the second filtered pixel, and the at least one channel .
40. A digital image processing system comprising: a computer readable storage medium; computer software stored on the computer readable storage medium and operable to: filter a first pixel of the at least one channel of the digital image to produce a filtered pixel the filtering of the first pixel responsive to the intensity of the first pixel and a plurality of additional pixels in the neighborhood of the first pixel, wherein at least one of the pixels used in filtering the first pixel is weighted in response to the intensity value of at least one pixel in the defect channel ; and produce a corrected digital image in response to the filtered pixel and the at least one channel.
41. The digital image processing system of Claim 40, wherein the digital image comprises three channels, each channel comprising a plurality of pixels; wherein the digital image represents a color image; and wherein the computer software is further operable to filter at least a second pixel at the same location in each of the three channels to produce a second filtered pixel for each of the three channels, the filtering of the second pixel in each channel, for a particular channel, responsive to the intensity of the second pixel in the particular channel and a plurality of additional pixels in the neighborhood of the second pixel in the particular channel, wherein at least one of the pixels used in filtering the second pixel in the particular channel is weighted in response to the intensity value of at least one pixel in the defect channel; produce a corrected digital image in response to the filtered pixel, the second filtered pixel for each channel, and the three channels.
42. The digital image processing system of Claim 40, wherein at least one of the pixels used in filtering the first pixel is weighted in response to the intensity value of a plurality of pixels in the defect channel.
43. The digital image processing system of Claim 40, wherein at least one of the pixels used in filtering the first pixel is weighted in response to the intensity value of a pixel in the defect channel which spatially corresponds to the first pixel .
44. The digital image scanning system of Claim 40, wherein the computer software is further operable to: filter the first pixel a plurality of times, wherein ones of the plurality of filtering operations cover a different bandwidth, are responsive to the intensity of the first pixel and a plurality of additional pixels in the neighborhood of the first pixel, and wherein at least one of the pixels use in filtering the first pixel is weighted in response to the intensity value of at least one pixel in the defect channel, each of the plurality of filtering operations producing a filtered pixel, collectively comprising a plurality of filtered pixels; and wherein the corrected digital image is further responsive to the plurality of filtered pixels.
45. The digital image scanning system of Claim 40, wherein the computer software is further operable to: assign a low weight value to a pixel for the filtering operation if a spatially corresponding pixel in the defect channel has an intensity indicating a high probability that the pixel is defective.
46. The digital image scanning system of Claim 40, wherein the computer software is further operable to: assign a low weight value to a pixel for the filtering operation if the average intensity of a spatially corresponding pixel in the defect channel and other pixels in the defect channel in the neighborhood of the spatially corresponding pixel indicates a high probability that the pixel is defective.
47. The digital image scanning system of Claim 40, wherein the computer software is further operable to: divide the digital image into a plurality of segments; wherein the first pixel is in a first segment comprising one of the plurality of segments of the at least one channel of the digital image and wherein the plurality of additional pixels in the neighborhood of the first pixel are in the first segment; filter a second pixel of the at least one channel of the digital image to produce a second filtered pixel the filtering of the second pixel responsive to the intensity of the second pixel and a plurality of additional pixels in the neighborhood of the second pixel, wherein at least one of the pixels used in filtering the second pixel is weighted in response to the intensity value of at least one pixel in the defect channel; wherein the second pixel is in a second segment comprising one of the plurality of segments of the at least one channel of the digital image and wherein the plurality of additional pixels in the neighborhood of the second pixel are in the second segment ; and wherein the corrected digital image is produced in response the filtered pixel, the second filtered pixel, and the at least one channel .
48. An altered digital image derived from a digital image having at least one channel, comprising: a computer readable storage medium; an altered digital image stored on the computer readable storage medium wherein the altered digital image was created by: filtering a first pixel of the at least one channel of the digital image to produce a filtered pixel the filtering of the first pixel responsive to the intensity of the first pixel and a plurality of additional pixels in the neighborhood of the first pixel, wherein at least one of the pixels used in filtering the first pixel is weighted in response to the intensity value of at least one pixel in the defect channel; and producing a corrected digital image in response to the filtered pixel and the at least one channel.
49. The altered digital image of Claim 48, wherein the digital image comprises three channels, each channel comprising a plurality of pixels; and wherein the digital image represents a color image and wherein the altered digital image was further created by: filtering at least a second pixel at the same location in each of the three channels to produce a second filtered pixel for each of the three channels, the filtering of the second pixel in each channel, for a particular channel, responsive to the intensity of the second pixel in the particular channel and a plurality of additional pixels in the neighborhood of the second pixel in the particular channel, wherein at least one of the pixels used in filtering the second pixel in the particular channel is weighted in response to the intensity value of at least one pixel in the defect channel; producing a corrected digital image in response to the filtered pixel, the second filtered pixel for each channel, and the three channels.
50. The altered digital image of Claim 48, wherein the altered digital image was further created by: filtering the first pixel a plurality of times, wherein ones of the plurality of filtering operations cover a different bandwidth, are responsive to the intensity of the first pixel and a plurality of additional pixels in the neighborhood of the first pixel, and wherein at least one of the pixels use in filtering the first pixel is weighted in response to the intensity value of at least one pixel in the defect channel, each of the plurality of filtering operations producing a filtered pixel, collectively comprising a plurality of filtered pixels; and wherein the corrected digital image is further responsive to the plurality of filtered pixels.
51. A method for altering defects in a digital image, comprising : storing on a computer readable storage medium a digital image comprising at least one channel, the at least one channel comprising a plurality of pixels each having an intensity value; storing on the computer readable storage medium a defect channel comprising a plurality of pixels each having an intensity value, at least some of the plurality of pixels in the defect channel having an intensity value proportional to defects in the digital image; creating, using digital circuitry, a plurality of filtered versions of the at least one channel, wherein pixels of the filtered versions of the at least one channel are determined by, for a particular pixel in a particular filtered version, computing a weighted average of the intensity values of the particular pixel and pixels within a region adjacent to the particular pixel in the at least one channel, wherein at least some of the intensity values used in the computation of the weighted average are weighted in response to the intensity value of at least one pixel in the defect channel; and producing a corrected digital image in response to the plurality of filtered versions of the at least one channel, the original at least one channel, and the defect channel.
52. The method of Claim 51, further comprising: creating, using digital circuitry, a plurality of filtered versions of the defect channel, wherein pixels of the filtered versions of the defect channel are determined by, for a particular pixel in a particular filtered version, computing a weighted average of the intensity values of the particular pixel and pixels within a region adjacent to the particular pixel, wherein at least some of the intensity values used in the computation of the weighted average are weighted in response to the intensity value of at least one pixel in the defect channel; and wherein the corrected digital image is further produced in response to the plurality of filtered versions of the defect channel .
53. The method of Claim 52, wherein the filtered versions of the at least one channel are created in log space .
54. The method of Claim 51, further comprising: dividing the at least one channel into segments; creating, using digital circuitry, a plurality of filtered versions of each segment of the at least one channel wherein pixels of the filtered versions of each segment of the at least one channel are determined by, for a particular pixel in a particular filtered version, computing a weighted average of the intensity values of the particular pixel and pixels within a region adjacent to the particular pixel and within the same segment as the particular pixel, and wherein at least some of the intensity values used in the computation of the weighted average are weighted in response to the intensity value of at least one pixel in the defect channel; and wherein the the corrected digital image is further produced in response to the plurality of filtered versions of each segment of the at least one channel .
55. The method of Claim 51, further comprising: prior to creating the plurality of filtered versions of the at least one channel, altering a plurality of pixels of the at least one channel by subtracting, for a particular pixel, a first value responsive to the intensity value of a pixel in the defect channel spatially corresponding to the particular pixel's position in the at least one channel from a second value responsive to the intensity value of the particular pixel .
56. The method of Claim 51, wherein the corrected digital image is produced by: creating a series of frequency bands, the frequency bands essentially contiguous and responsive to the at least one channel and the plurality of filtered versions of the at least one channel; and adding the frequency bands together to form the corrected digital image.
57. The method of Claim 55, wherein the corrected digital image is produced by: creating a series of frequency bands, the frequency bands essentially contiguous and responsive to the altered version of the at least one channel and the plurality of filtered versions of the altered version of the at least one channel ; and adding the frequency bands together to form the corrected digital image.
58. The method of Claim 52, wherein the corrected digital image is produced by: creating a first series of frequency bands, the first series of frequency bands essentially contiguous and responsive to the at least one channel and the plurality of filtered versions of the at least one channel; creating a second series of frequency bands, the second series of frequency bands essentially contiguous and responsive to the defect channel and the plurality of filtered versions of the defect channel and wherein each of the second series of frequency bands corresponds to one of the frequency bands of the first series of frequency bands; subtracting each corresponding frequency band in the second series of frequency bands from a frequency band in the first series of frequency bands to form a third series of frequency bands; adding the third series of frequency bands together to form the corrected digital image .
59. The method of Claim 58, wherein the subtracting comprises a bounded subtraction wherein the bounded subtraction, for a particular pixel in a particular frequency band, comprises the selection of a result of a plurality of subtractions including a first subtraction proportional to the difference between a first value proportional to the intensity of the particular pixel in the particular frequency band in the first series of frequency bands and a second value proportional to the intensity of the particular pixel in the particular frequency band in the second series of frequency bands and a second subtraction proportional to the difference between the first value and a third value proportional to the intensity of the particular pixel in the particular frequency band in the second series of frequency bands.
60. The method of Claim 58, wherein the second series of frequency bands is further responsive to red residue in the defect channel such that the second series of frequency bands is adjusted to reduce red residue in the defect channel .
61. The method of Claim 58, wherein the weighted average comprise an average selected from the group consisting of a median average, an arithmetic average, a geometric average, a mean average, and a mode average.
62. The method of Claim 1, wherein the first pixel is filtered using a window filter with feathered edges.
63. The method of Claim 1, further comprising: filtering a second pixel comprising a part of the defect channel to produce a second filtered pixel the filtering of the second pixel responsive to the intensity of the second pixel and a plurality of additional pixels in the neighborhood of the second pixel, wherein at least one of the pixels used in filtering the second pixel is weighted in response to the intensity value of at least one pixel in the defect channel .
64. The method of Claim 63 , wherein the corrected digital image is produced in response to the filtered pixel, the second filtered pixel, the at least one channel, and the defect channel .
65. The method of Claim 1, further comprising: prior to filtering the first pixel, altering at least a second pixel of the at least one channel by subtracting a first value responsive to the intensity value of a pixel in the defect channel spatially corresponding to the second pixel's position in the at least one channel from a second value responsive to the intensity value of the second pixel.
PCT/US2000/025478 1999-09-16 2000-09-15 Method and system for altering defects in a digital image WO2001020898A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
AU75868/00A AU7586800A (en) 1999-09-16 2000-09-15 Method and system for altering defects in a digital image
EP00965090A EP1221256A1 (en) 1999-09-16 2000-09-15 Method and system for altering defects in a digital image

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US15425599P 1999-09-16 1999-09-16
US60/154,255 1999-09-16

Publications (1)

Publication Number Publication Date
WO2001020898A1 true WO2001020898A1 (en) 2001-03-22

Family

ID=22550627

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2000/025478 WO2001020898A1 (en) 1999-09-16 2000-09-15 Method and system for altering defects in a digital image

Country Status (5)

Country Link
US (2) US6487321B1 (en)
EP (1) EP1221256A1 (en)
CN (1) CN1375161A (en)
AU (1) AU7586800A (en)
WO (1) WO2001020898A1 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1343306A2 (en) * 2002-03-04 2003-09-10 Noritsu Koki Co., Ltd. Method, program and apparatus for image processing
EP1511295A2 (en) * 2003-08-29 2005-03-02 Noritsu Koki Co., Ltd. Image processing apparatus and image processing method
CN106570838A (en) * 2016-10-28 2017-04-19 成都西纬科技有限公司 Image brightness optimization method and device

Families Citing this family (50)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TW376627B (en) * 1998-08-05 1999-12-11 Transpacific Optics Llc Dark line modifying device and method for image scanning
US6771833B1 (en) * 1999-08-20 2004-08-03 Eastman Kodak Company Method and system for enhancing digital images
US6711302B1 (en) * 1999-10-20 2004-03-23 Eastman Kodak Company Method and system for altering defects in digital image
JP3620391B2 (en) * 2000-02-23 2005-02-16 日本電気株式会社 Fingerprint input device, image determination method used therefor, and recording medium recording control program thereof
US7113627B1 (en) * 2000-08-09 2006-09-26 Eastman Kodak Company Location of extended linear defects
JPWO2003010525A1 (en) * 2001-07-27 2004-11-18 日本板硝子株式会社 Method for evaluating contamination on the surface of an object and a photographing box used for the method
US7986820B2 (en) * 2001-10-19 2011-07-26 Mitsubishi Electric Research Laboratories, Inc. Method for comparing features extracted from images of fingerprints
US6956976B2 (en) * 2002-01-04 2005-10-18 Warner Bros. Enterianment Inc. Reduction of differential resolution of separations
US6947607B2 (en) * 2002-01-04 2005-09-20 Warner Bros. Entertainment Inc. Reduction of differential resolution of separations
US7031549B2 (en) * 2002-02-22 2006-04-18 Hewlett-Packard Development Company, L.P. Systems and methods for enhancing tone reproduction
CN100399792C (en) * 2002-12-05 2008-07-02 瑞昱半导体股份有限公司 Method and apparatus for removing slide fastener shaped faintness of boundary image
US7639741B1 (en) * 2002-12-06 2009-12-29 Altera Corporation Temporal filtering using object motion estimation
US20040161153A1 (en) * 2003-02-18 2004-08-19 Michael Lindenbaum Context-based detection of structured defects in an image
US7702236B2 (en) 2006-02-14 2010-04-20 Fotonation Vision Limited Digital image acquisition device with built in dust and sensor mapping capability
JP4639037B2 (en) * 2003-07-18 2011-02-23 キヤノン株式会社 Image processing method and apparatus
US7295233B2 (en) 2003-09-30 2007-11-13 Fotonation Vision Limited Detection and removal of blemishes in digital images utilizing original images of defocused scenes
US7315658B2 (en) * 2003-09-30 2008-01-01 Fotonation Vision Limited Digital camera
US7310450B2 (en) * 2003-09-30 2007-12-18 Fotonation Vision Limited Method of detecting and correcting dust in digital images based on aura and shadow region analysis
US7424170B2 (en) 2003-09-30 2008-09-09 Fotonation Vision Limited Automated statistical self-calibrating detection and removal of blemishes in digital images based on determining probabilities based on image analysis of single images
US7369712B2 (en) * 2003-09-30 2008-05-06 Fotonation Vision Limited Automated statistical self-calibrating detection and removal of blemishes in digital images based on multiple occurrences of dust in images
US7676110B2 (en) * 2003-09-30 2010-03-09 Fotonation Vision Limited Determination of need to service a camera based on detection of blemishes in digital images
US7590305B2 (en) * 2003-09-30 2009-09-15 Fotonation Vision Limited Digital camera with built-in lens calibration table
US7206461B2 (en) * 2003-09-30 2007-04-17 Fotonation Vision Limited Digital image acquisition and processing system
US8369650B2 (en) 2003-09-30 2013-02-05 DigitalOptics Corporation Europe Limited Image defect map creation using batches of digital images
US7340109B2 (en) 2003-09-30 2008-03-04 Fotonation Vision Limited Automated statistical self-calibrating detection and removal of blemishes in digital images dependent upon changes in extracted parameter values
US7308156B2 (en) * 2003-09-30 2007-12-11 Fotonation Vision Limited Automated statistical self-calibrating detection and removal of blemishes in digital images based on a dust map developed from actual image data
US7822285B2 (en) * 2004-05-20 2010-10-26 Omnivision Technologies, Inc. Methods and systems for locally adaptive image processing filters
US20050276512A1 (en) * 2004-06-01 2005-12-15 3M Innovative Properties Company Selective deconvolution of an image
US7471820B2 (en) * 2004-08-31 2008-12-30 Aptina Imaging Corporation Correction method for defects in imagers
US7450778B2 (en) * 2005-03-17 2008-11-11 Hewlett-Packard Development Company, L.P. Artifact reduction in a digital video
US7969488B2 (en) * 2005-08-03 2011-06-28 Micron Technologies, Inc. Correction of cluster defects in imagers
US7590303B2 (en) * 2005-09-29 2009-09-15 Samsung Electronics Co., Ltd. Image enhancement method using local illumination correction
US8630498B2 (en) * 2006-03-02 2014-01-14 Sharp Laboratories Of America, Inc. Methods and systems for detecting pictorial regions in digital images
US7864365B2 (en) * 2006-06-15 2011-01-04 Sharp Laboratories Of America, Inc. Methods and systems for segmenting a digital image into regions
US8437054B2 (en) * 2006-06-15 2013-05-07 Sharp Laboratories Of America, Inc. Methods and systems for identifying regions of substantially uniform color in a digital image
US7932938B2 (en) * 2006-08-25 2011-04-26 Micron Technology, Inc. Method, apparatus and system providing adjustment of pixel defect map
US7876959B2 (en) * 2006-09-06 2011-01-25 Sharp Laboratories Of America, Inc. Methods and systems for identifying text in digital images
JP4913075B2 (en) * 2008-01-16 2012-04-11 京セラミタ株式会社 Image reading device
US7961224B2 (en) 2008-01-25 2011-06-14 Peter N. Cheimets Photon counting imaging system
WO2009120830A1 (en) * 2008-03-27 2009-10-01 Tufts University Methods and apparatus for visual sub-band decomposition of signals
WO2009126621A2 (en) * 2008-04-07 2009-10-15 Tufts University Methods and apparatus for image restoration
JP5568884B2 (en) * 2009-04-02 2014-08-13 セイコーエプソン株式会社 Video processing apparatus and video processing method
JP5481908B2 (en) * 2009-04-02 2014-04-23 セイコーエプソン株式会社 Image reading apparatus, control apparatus therefor, program, and control method
GB2479414B (en) * 2010-04-09 2015-11-11 Snell Ltd Repairing scratch impairments to an image
US8351696B2 (en) * 2010-09-17 2013-01-08 Foveon, Inc. Correcting defective pixels in digital color images
WO2013025204A1 (en) * 2011-08-16 2013-02-21 Nikon Corporation Method for detecting existence of dust spots in digital images based on locally adaptive thresholding
US8654398B2 (en) 2012-03-19 2014-02-18 Seiko Epson Corporation Method for simulating impact printer output, evaluating print quality, and creating teaching print samples
US10275971B2 (en) * 2016-04-22 2019-04-30 Ncr Corporation Image correction
US11036978B2 (en) * 2018-05-29 2021-06-15 University Of Electronic Science And Technology Of China Method for separating out a defect image from a thermogram sequence based on weighted naive bayesian classifier and dynamic multi-objective optimization
JP7087792B2 (en) * 2018-07-31 2022-06-21 株式会社リコー Image data generator, image data generation method and program

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1990001240A1 (en) * 1988-07-25 1990-02-08 Eastman Kodak Company Film noise reduction by application of bayes theorem to positive/negative film
WO1991009493A1 (en) * 1989-12-18 1991-06-27 Eastman Kodak Company Method for deriving noise-reduced estimates of colour signal parameters from multiple colour/luminance image sensor outputs
EP0569142A1 (en) * 1992-05-05 1993-11-10 International Business Machines Corporation System and method for image recovery
EP0816833A2 (en) * 1996-06-28 1998-01-07 Fujifilm Electronic Imaging Limited An illumination unit and illumination method

Family Cites Families (93)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2821868A1 (en) 1978-05-19 1979-11-22 Karl Sirowatka Recording damage on moving film or bands - using infrared scanning and pen recorders or warning devices
US4302108A (en) 1979-01-29 1981-11-24 Polaroid Corporation Detection of subsurface defects by reflection interference
US4260899A (en) 1979-06-14 1981-04-07 Intec Corporation Wide web laser scanner flaw detection method and apparatus
US4301469A (en) 1980-04-30 1981-11-17 United Technologies Corporation Run length encoder for color raster scanner
US4462860A (en) 1982-05-24 1984-07-31 At&T Bell Laboratories End point detection
GB2124449B (en) 1982-07-16 1986-11-19 British Broadcasting Corp Concealment of defects in a video signal
US4442454A (en) 1982-11-15 1984-04-10 Eastman Kodak Company Image processing method using a block overlap transformation procedure
JPS60146567A (en) 1984-01-10 1985-08-02 Sharp Corp Color picture reader
JPS61275625A (en) 1985-05-31 1986-12-05 Fuji Photo Film Co Ltd Calibrating method for color photographic image information
NL8501956A (en) 1985-07-09 1987-02-02 Philips Nv IMAGE RECOVERY CIRCUIT.
DE3534019A1 (en) 1985-09-24 1987-04-02 Sick Optik Elektronik Erwin OPTICAL RAILWAY MONITORING DEVICE
US4677465A (en) 1985-11-01 1987-06-30 Eastman Kodak Company Digital color image processing method with shape correction of histograms used to produce color reproduction functions
JPS62116937A (en) 1985-11-16 1987-05-28 Dainippon Screen Mfg Co Ltd Film attaching and detaching device for drum type image scanning and recording device
JPH01143945A (en) 1987-11-30 1989-06-06 Fuji Photo Film Co Ltd Detecting method for defect in tape
US4858003A (en) 1988-01-12 1989-08-15 Eastman Kodak Company Mechanism for handling slides and film strips
US5194950A (en) 1988-02-29 1993-03-16 Mitsubishi Denki Kabushiki Kaisha Vector quantizer
US5047968A (en) 1988-03-04 1991-09-10 University Of Massachusetts Medical Center Iterative image restoration device
CA1309166C (en) 1988-05-20 1992-10-20 Toshinobu Haruki Image sensing apparatus having automatic iris function of automatically adjusting exposure in response to video signal
US5010401A (en) 1988-08-11 1991-04-23 Mitsubishi Denki Kabushiki Kaisha Picture coding and decoding apparatus using vector quantization
US4989973A (en) 1988-10-03 1991-02-05 Nissan Motor Co., Ltd. Surface condition estimating apparatus
WO1990011050A1 (en) 1989-03-28 1990-10-04 Yokogawa Medical Systems, Ltd. Image processing apparatus
US4994918A (en) 1989-04-28 1991-02-19 Bts Broadcast Television Systems Gmbh Method and circuit for the automatic correction of errors in image steadiness during film scanning
US4972091A (en) 1989-05-16 1990-11-20 Canadian Patents And Development Limited/Societe Canadienne Des Brevets Et D'exploitation Limitee Method and apparatus for detecting the presence of flaws in a moving sheet of material
US5058982A (en) 1989-06-21 1991-10-22 Orbot Systems Ltd. Illumination system and inspection apparatus including same
US4937720A (en) 1989-10-13 1990-06-26 Sundstrand Corporation PWM inverter circuit analytically compensating for DC link distortion
JPH07115534A (en) 1993-10-15 1995-05-02 Minolta Co Ltd Image reader
US5003379A (en) 1989-10-16 1991-03-26 Eastman Kodak Company Telecine scanning apparatus with spectrally-shifted sensitivities responsive to negative or print film dyes
DE3942273A1 (en) 1989-12-21 1991-06-27 Broadcast Television Syst METHOD FOR HIDDEN ERRORS IN A VIDEO SIGNAL
US5267030A (en) 1989-12-22 1993-11-30 Eastman Kodak Company Method and associated apparatus for forming image data metrics which achieve media compatibility for subsequent imaging application
US5091972A (en) 1990-09-17 1992-02-25 Eastman Kodak Company System and method for reducing digital image noise
GB9023013D0 (en) 1990-10-23 1990-12-05 Crosfield Electronics Ltd Method and apparatus for generating representation of an image
JP2528383B2 (en) 1990-11-22 1996-08-28 大日本スクリーン製造株式会社 Pinhole erasing method
US5155596A (en) 1990-12-03 1992-10-13 Eastman Kodak Company Film scanner illumination system having an automatic light control
US5465163A (en) 1991-03-18 1995-11-07 Canon Kabushiki Kaisha Image processing method and apparatus for processing oversized original images and for synthesizing multiple images
JPH04291139A (en) 1991-03-20 1992-10-15 Nippon Steel Corp Method for reporting portion with flaw on strip
US5452018A (en) 1991-04-19 1995-09-19 Sony Electronics Inc. Digital color correction system having gross and fine adjustment modes
US5149960B1 (en) 1991-07-03 1994-08-30 Donnelly R R & Sons Method of converting scanner signals into colorimetric signals
EP0527097A3 (en) 1991-08-06 1995-03-01 Eastman Kodak Co Apparatus and method for collectively performing tile-based image rotation, scaling and digital halftone screening
US5200817A (en) 1991-08-29 1993-04-06 Xerox Corporation Conversion of an RGB color scanner into a colorimetric scanner
JP2549479B2 (en) 1991-12-06 1996-10-30 日本電信電話株式会社 Motion compensation inter-frame band division coding processing method
US5371542A (en) 1992-06-23 1994-12-06 The United States Of America As Represented By The Secretary Of The Navy Dual waveband signal processing system
CA2093840C (en) 1992-07-17 1999-08-10 Albert D. Edgar Duplex film scanning
US5583950A (en) 1992-09-16 1996-12-10 Mikos, Ltd. Method and apparatus for flash correlation
US5300381A (en) 1992-09-24 1994-04-05 Eastman Kodak Company Color image reproduction of scenes with preferential tone mapping
US5568270A (en) 1992-12-09 1996-10-22 Fuji Photo Film Co., Ltd. Image reading apparatus which varies reading time according to image density
JPH06311307A (en) 1993-04-22 1994-11-04 Minolta Camera Co Ltd Image forming device
EP0624848A3 (en) 1993-05-04 1994-11-30 Eastman Kodak Company A technique for the detection and removal of local defects in digital continuous-tone images
KR950004881A (en) 1993-07-31 1995-02-18 김광호 Color image processing method and device
JPH0766977A (en) 1993-08-24 1995-03-10 Minolta Co Ltd Picture processing unit
US6128416A (en) 1993-09-10 2000-10-03 Olympus Optical Co., Ltd. Image composing technique for optimally composing a single image from a plurality of digital images
GB2283633B (en) 1993-11-05 1997-10-29 Sony Uk Ltd Anti-alias filter control for a split picture
US5729631A (en) 1993-11-30 1998-03-17 Polaroid Corporation Image noise reduction system using a wiener variant filter in a pyramid image representation
US5477345A (en) 1993-12-15 1995-12-19 Xerox Corporation Apparatus for subsampling chrominance
US5509086A (en) 1993-12-23 1996-04-16 International Business Machines Corporation Automatic cross color elimination
KR100300950B1 (en) 1994-01-31 2001-10-22 윤종용 Method and apparatus for correcting color
US5516608A (en) 1994-02-28 1996-05-14 International Business Machines Corporation Method for controlling a line dimension arising in photolithographic processes
EP0669753A3 (en) 1994-02-28 1995-12-20 Minolta Co Ltd An apparatus for reproducing images.
DE4424577A1 (en) 1994-07-13 1996-01-18 Hoechst Ag Transporter protein for cationic xenobiotic(s) and pharmaceuticals, and related DNA and transformed cells
JPH0877341A (en) 1994-08-29 1996-03-22 Xerox Corp Equipment and method for color image processing
DE4432787A1 (en) 1994-09-15 1996-03-21 Philips Patentverwaltung Method and circuit for detecting and masking errors in a video signal
US5561611A (en) 1994-10-04 1996-10-01 Noran Instruments, Inc. Method and apparatus for signal restoration without knowledge of the impulse response function of the signal acquisition system
EP0707408A1 (en) 1994-10-11 1996-04-17 International Business Machines Corporation Optical scanner device for transparent media
US5565931A (en) 1994-10-31 1996-10-15 Vivo Software. Inc. Method and apparatus for applying gamma predistortion to a color image signal
US5649032A (en) 1994-11-14 1997-07-15 David Sarnoff Research Center, Inc. System for automatically aligning images to form a mosaic image
CH690639A5 (en) 1994-11-29 2000-11-15 Zeiss Carl Fa Apparatus for scanning digitize image templates and processes for their operation.
US6865000B2 (en) 1994-12-06 2005-03-08 Canon Kabushiki Kaisha Image reading apparatus for grouping sensors according to characteristics
US5771107A (en) 1995-01-11 1998-06-23 Mita Industrial Co., Ltd. Image processor with image edge emphasizing capability
US5979011A (en) 1995-04-07 1999-11-09 Noritsu Koki Co., Ltd Dust removing apparatus
US5582961A (en) 1995-06-06 1996-12-10 Eastman Kodak Company Photographic elements which achieve colorimetrically accurate recording
US5710643A (en) 1995-06-29 1998-01-20 Agfa Divisionn, Bayer Corporation Optical path for a scanning system
US6104839A (en) 1995-10-16 2000-08-15 Eastman Kodak Company Method and apparatus for correcting pixel values in a digital image
JP3669448B2 (en) 1995-10-31 2005-07-06 富士写真フイルム株式会社 Image reproduction method and apparatus
US5641596A (en) 1995-12-05 1997-06-24 Eastman Kodak Company Adjusting film grain properties in digital images
US5892595A (en) 1996-01-26 1999-04-06 Ricoh Company, Ltd. Image reading apparatus for correct positioning of color component values of each picture element
EP0794454B1 (en) 1996-03-04 2005-05-11 Fuji Photo Film Co., Ltd. Film scanner
US5982951A (en) 1996-05-28 1999-11-09 Canon Kabushiki Kaisha Apparatus and method for combining a plurality of images
US5963662A (en) 1996-08-07 1999-10-05 Georgia Tech Research Corporation Inspection system and method for bond detection and validation of surface mount devices
US6075905A (en) 1996-07-17 2000-06-13 Sarnoff Corporation Method and apparatus for mosaic image construction
US5808674A (en) 1996-09-09 1998-09-15 Eastman Kodak Company Producing and improved digital image from digital signals corresponding to pairs of photosites
JPH10178564A (en) 1996-10-17 1998-06-30 Sharp Corp Panorama image generator and recording medium
JP3493104B2 (en) 1996-10-24 2004-02-03 シャープ株式会社 Color image processing equipment
US5845017A (en) 1996-12-17 1998-12-01 Eastman Kodak Company Digital image processing method for degraining of film images using distance weighted averaging of target pixel code values
US5982941A (en) 1997-02-07 1999-11-09 Eastman Kodak Company Method of producing digital image with improved performance characteristic
JPH10291139A (en) 1997-04-18 1998-11-04 Amada Co Ltd Work carrying-in and out device for plate material machining machine
EP0893914A3 (en) 1997-07-24 2002-01-02 Nikon Corporation Image processing method, image processing apparatus, and storage medium for storing control process
US6078701A (en) 1997-08-01 2000-06-20 Sarnoff Corporation Method and apparatus for performing local to global multiframe alignment to construct mosaic images
US5969372A (en) 1997-10-14 1999-10-19 Hewlett-Packard Company Film scanner with dust and scratch correction by use of dark-field illumination
US6239886B1 (en) 1998-01-08 2001-05-29 Xerox Corporation Method and apparatus for correcting luminance and chrominance data in digital color images
US6078051A (en) 1998-01-08 2000-06-20 Xerox Corporation Image input device and method for providing scanning artifact detection
US6057040A (en) 1998-01-22 2000-05-02 Vision--Ease Lens, Inc. Aminosilane coating composition and process for producing coated articles
JP4096407B2 (en) 1998-04-22 2008-06-04 株式会社ニコン Image processing apparatus and computer-readable storage medium
JP2000196813A (en) 1998-12-25 2000-07-14 Canon Inc Image reader
ES2234554T3 (en) 1999-07-20 2005-07-01 Tecnoforming S.P.A. LIGHT ALLOY RIM WITH A STAINLESS STEEL FRONT COVERING ELEMENT.

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1990001240A1 (en) * 1988-07-25 1990-02-08 Eastman Kodak Company Film noise reduction by application of bayes theorem to positive/negative film
WO1991009493A1 (en) * 1989-12-18 1991-06-27 Eastman Kodak Company Method for deriving noise-reduced estimates of colour signal parameters from multiple colour/luminance image sensor outputs
EP0569142A1 (en) * 1992-05-05 1993-11-10 International Business Machines Corporation System and method for image recovery
EP0816833A2 (en) * 1996-06-28 1998-01-07 Fujifilm Electronic Imaging Limited An illumination unit and illumination method

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1343306A2 (en) * 2002-03-04 2003-09-10 Noritsu Koki Co., Ltd. Method, program and apparatus for image processing
EP1343306A3 (en) * 2002-03-04 2004-12-15 Noritsu Koki Co., Ltd. Method, program and apparatus for image processing
US7123775B2 (en) 2002-03-04 2006-10-17 Noritsu Koki Co., Ltd. Image processing method, image processing program and image processing apparatus
EP1511295A2 (en) * 2003-08-29 2005-03-02 Noritsu Koki Co., Ltd. Image processing apparatus and image processing method
EP1511295A3 (en) * 2003-08-29 2005-03-23 Noritsu Koki Co., Ltd. Image processing apparatus and image processing method
US7327877B2 (en) 2003-08-29 2008-02-05 Noritsu Koki Co., Ltd. Image processing apparatus and image processing method
CN106570838A (en) * 2016-10-28 2017-04-19 成都西纬科技有限公司 Image brightness optimization method and device
CN106570838B (en) * 2016-10-28 2019-05-03 成都西纬科技有限公司 A kind of brightness of image optimization method and device

Also Published As

Publication number Publication date
CN1375161A (en) 2002-10-16
US20020106133A1 (en) 2002-08-08
EP1221256A1 (en) 2002-07-10
US6650789B2 (en) 2003-11-18
AU7586800A (en) 2001-04-17
US6487321B1 (en) 2002-11-26

Similar Documents

Publication Publication Date Title
US6487321B1 (en) Method and system for altering defects in a digital image
US6711302B1 (en) Method and system for altering defects in digital image
US6285798B1 (en) Automatic tone adjustment by contrast gain-control on edges
US7016080B2 (en) Method and system for improving scanned image detail
US6771833B1 (en) Method and system for enhancing digital images
US8594451B2 (en) Edge mapping incorporating panchromatic pixels
US7876974B2 (en) Method for improving digital images and an image sensor for sensing the same
EP2130176B1 (en) Edge mapping using panchromatic pixels
US8180168B2 (en) One-pass filtering and infrared-visible light decorrelation to reduce noise and distortions
US20040207881A1 (en) Image processing method, image processing apparatus and image processing program
US6075590A (en) Reflection infrared surface defect correction
EP1742178A2 (en) Contrast enhancement of images
US20050249430A1 (en) Image quality improving apparatus and method
JP2004147333A (en) Method of enhancing tonal and spatial characteristic of digital image using selective spatial filter
US7031543B2 (en) Grayscale image de-speckle algorithm
US6792162B1 (en) Method and apparatus to automatically enhance the quality of digital images by measuring grain trace magnitudes
CN101431593A (en) Image processing apparatus, image forming apparatus and image processing method
US20120188418A1 (en) Texture detection in image processing
US8189236B2 (en) Image processing apparatus, image processing method and computer-readable medium
JP3700357B2 (en) Image processing device
US20060056722A1 (en) Edge preserving method and apparatus for image processing
JP2003132352A (en) Image processing method, image processor, image processing program and computer-readable recording medium recording it
US20040169872A1 (en) Blind inverse halftoning
US20040213477A1 (en) Image-processing method and apparatus, computer program for executing image processing and image-recording apparatus
JP2004318425A (en) Image processing method, image processor, and program

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AT AU AZ BA BB BG BR BY BZ CA CH CN CR CU CZ CZ DE DE DK DK DM DZ EE EE ES FI FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NO NZ PL PT RO RU SD SE SG SI SK SK SL TJ TM TR TT TZ UA UG UZ VN YU ZA ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
WWE Wipo information: entry into national phase

Ref document number: 2000965090

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 008129320

Country of ref document: CN

WWP Wipo information: published in national office

Ref document number: 2000965090

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: JP