WO2015103953A1 - Method and device for extracting iris image image under condition of non-uniform illumination - Google Patents

Method and device for extracting iris image image under condition of non-uniform illumination Download PDF

Info

Publication number
WO2015103953A1
WO2015103953A1 PCT/CN2015/070058 CN2015070058W WO2015103953A1 WO 2015103953 A1 WO2015103953 A1 WO 2015103953A1 CN 2015070058 W CN2015070058 W CN 2015070058W WO 2015103953 A1 WO2015103953 A1 WO 2015103953A1
Authority
WO
WIPO (PCT)
Prior art keywords
iris
image
channel
component
extraction
Prior art date
Application number
PCT/CN2015/070058
Other languages
French (fr)
Chinese (zh)
Inventor
孙贤军
董记平
夏鹏
Original Assignee
上海帝仪科技有限公司
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 上海帝仪科技有限公司 filed Critical 上海帝仪科技有限公司
Publication of WO2015103953A1 publication Critical patent/WO2015103953A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/136Segmentation; Edge detection involving thresholding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/18Eye characteristics, e.g. of the iris
    • G06V40/193Preprocessing; Feature extraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10024Color image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20056Discrete and fast Fourier transform, [DFT, FFT]

Definitions

  • Embodiments of the present invention generally relate to the field of image processing, and more particularly to an iris extraction method and apparatus under uneven illumination conditions.
  • Biometric identification technology refers to the use of the physiological characteristics or behavioral characteristics inherent in the human body for personal identification.
  • the iris identification technology is an important branch of biometric identification technology.
  • the center of his eye is a black pupil
  • the iris is the annular tissue between the outer edges of the pupil, which presents interlaced texture features similar to spots, filaments, stripes, and crypts.
  • the iris is hardly changed in one's life, and the iris of different people is completely different.
  • Iris identification technology is the application of computer image processing technology and pattern recognition technology in the field of personal identification. Because of its high stability and high accuracy, it can also make people get rid of the cumbersome memory credit card number, bank account number, ID number, and network access number. Therefore, it is widely used in banking, public security, airport, network security and other industries. .
  • the typical iris identification is mainly composed of image acquisition, image preprocessing, feature coding and feature matching.
  • image preprocessing iris extraction is especially critical, and its execution time and accuracy will directly affect the recognition speed and accuracy of the entire iris identification process.
  • Traditional iris extraction methods include Hough transform circle detection method, Daugman's method based on differential integral operator, and boundary detection method adopted by Wildes et al. They both model the iris boundary as a circular ring with two inner and outer boundaries, and extract the iris by extracting the inner and outer circles of the ring.
  • the iris is usually not a complete circular ring since the upper eyelid always blocks the upper half of the iris.
  • the outer boundary of the iris is often It is more ambiguous, that is, the transition between the iris and the sclera (the white part of the sclera, that is, the outer part of the eyeball, which is the outermost layer of the eye) is not obvious, so it is difficult to detect the outer boundary by the above-mentioned conventional boundary detection method.
  • the traditional method is time-consuming and laborious, and it is difficult to meet the requirements of real-time fast and accurate detection.
  • iris extraction In the invention patent application CN103246871A, the inventor first uses the existing Daugman method to extract the inner boundary of the image, and then selects the rectangular regions on the left and right sides of the inner boundary such that the rectangular region covers the left and right sides outside the inner boundary as much as possible. The iris is then gray-scale transformed into the rectangular region, and the Canny operator is used to detect the boundary of the rectangular region. Specifically, when performing boundary detection on a rectangular area, a pixel value lower than the threshold is set to 0 according to a predetermined threshold, and otherwise set to 1, thereby obtaining a boundary point set. It can be seen that in the invention patent application CN103246871A, a common image processing basic operation means called threshold segmentation is used to extract the target of interest.
  • threshold segmentation a common image processing basic operation means called threshold segmentation is used to extract the target of interest.
  • embodiments of the present invention provide an iris extraction method and apparatus under uneven illumination conditions.
  • an iris extraction method under uneven illumination conditions comprising: filtering out an uneven illumination intensity component from an image containing an iris; calculating a color difference in the filtered image; Iris extraction is performed according to the color difference.
  • the filtering out the uneven illumination intensity component from the image containing the iris comprises: converting a color mode of the image such that the uneven illumination intensity component in the converted image and The color components are separated.
  • the image is converted from an RGB color mode to a Lab color mode to obtain corresponding image components in the L channel, the a channel, and the b channel such that The uneven illumination intensity component and the color component in the converted image are separated, wherein the image component in the L channel includes the uneven illumination intensity component, an image component in the a channel, and the The image component in the b channel includes the color component.
  • the calculating the color difference in the filtered image comprises calculating the color difference using an image component in the a channel and an image component in the b channel.
  • the color difference is calculated using image components in the a channel and image components in the b channel according to the following formula: Where ⁇ E ab represents the chromatic aberration, ⁇ (a) represents the difference between the two values in the a channel, and ⁇ (b) represents the difference between the two values in the b channel.
  • the performing iris extraction according to the color difference comprises: comparing the color difference with a predetermined color difference threshold to obtain a first comparison result; and performing iris extraction according to the first comparison result.
  • the filtering out the uneven illumination intensity component from the image comprising the iris further comprises filtering the non-uniform illumination intensity component from the image components in the L channel.
  • filtering the uneven illumination intensity component from the image components in the L channel comprises: performing a logarithmic operation on the image components in the L channel; the log transformed image Performing a fast Fourier transform on the component; filtering the low frequency portion of the image component after the fast Fourier transform by high-pass filtering; performing inverse transform of the high-pass filtered image component on the fast Fourier transform; The inversely transformed image component is subjected to an exponential operation.
  • the performing iris extraction according to the color difference further comprises: comparing the color difference with a predetermined color difference threshold to obtain a first comparison result; and calculating the filtering after filtering the uneven light intensity component a gray value corresponding to the image component in the L channel; comparing the gray value with a predetermined gray threshold to obtain a second comparison result; calculating the first comparison result and the first of the second comparison result Intersection; and performing iris extraction based on the first intersection.
  • the performing iris extraction according to the color difference further comprises: calculating a brightness of any point in the image component in the L channel after filtering the uneven light intensity component minus a point around the point a difference between the average values of the brightness of the predetermined number of points; comparing the difference with a predetermined brightness difference threshold to obtain a third comparison result; calculating the first intersection and the second of the third comparison result Intersection; and performing iris extraction based on the second intersection.
  • the difference is calculated according to the following formula: Where (x, y) represents the coordinates of any point in the image component in the L channel after filtering out the uneven illumination intensity component, and f(x, y) represents the point of the coordinate (x, y) Brightness value, Represents a non-negative integer set and n represents a non-negative integer.
  • performing iris extraction according to the second intersection includes: determining, respectively, whether each foreground pixel in the binary image corresponding to the second intersection has an adjacent foreground pixel, and if so, Deriving a foreground pixel in the same target area as the adjacent foreground pixel to obtain a plurality of target areas; respectively calculating a number of foreground pixel points owned by the plurality of target areas; and the target area having the largest number Determined to be an iris region; and perform iris extraction in the iris region.
  • performing iris extraction in the iris region includes: eliminating holes in the iris region, non-adjacent boundaries, resulting in a complete iris region; and performing iris extraction in the intact iris region.
  • the holes in the iris region, non-adjacent boundaries are eliminated according to the following formula:
  • A represents the iris region
  • represents any point in A
  • B represents a collection of structural elements
  • b represents any point in B. Represents a two-dimensional integer grid.
  • the holes in the iris region are also eliminated according to the following formula, Adjacent borders are based on:
  • A represents the iris region
  • represents any point in A
  • B represents a collection of structural elements. Represents a two-dimensional integer grid.
  • an iris extraction apparatus under uneven illumination conditions, comprising: a filtering device for filtering out uneven illumination intensity components from an image containing an iris; and a computing device for Calculating the color difference in the filtered image; and an iris extracting means for performing iris extraction based on the color difference.
  • the filtering device includes a conversion unit for converting a color mode of the image such that the uneven illumination intensity component and the color component in the converted image are separated.
  • the conversion unit converts the image from an RGB color mode to a Lab color mode to obtain corresponding image components in the L channel, the a channel, and the b channel, such that in the converted image
  • the uneven illumination intensity component and the color component are separated, wherein the image component in the L channel includes the uneven illumination intensity component, and the image component in the a channel and the image component in the b channel include The color component.
  • the computing device includes a first computing unit for calculating the color difference using image components in the a channel and image components in the b channel.
  • the first calculating unit calculates the color difference by using an image component in the a channel and an image component in the b channel according to the following formula: Where ⁇ E ab represents the chromatic aberration, ⁇ (a) represents the difference between the two values in the a channel, and ⁇ (b) represents the difference between the two values in the b channel.
  • the iris extraction device includes: a first comparison unit configured to compare the color difference with a predetermined color difference threshold to obtain a first comparison result; and an iris extraction unit configured to The results were compared for iris extraction.
  • the filtering device further comprises: a filtering unit for filtering out the uneven illumination intensity component from image components in the L channel.
  • the filtering unit includes: a first operation module, configured to perform a logarithm operation on image components in the L channel; and a second operation module, configured to convert the logarithmically transformed image
  • the component performs a fast Fourier transform;
  • a third operation module is configured to filter the low frequency portion of the image component after the fast Fourier transform by high-pass filtering; and
  • a fourth operation module configured to filter the high-pass
  • the image component performs an inverse transform of the fast Fourier transform; and a fifth operation module is configured to perform an exponential operation on the inverse transformed image component.
  • the iris extraction device includes: a first comparison unit, configured to compare the color difference with a predetermined color difference threshold to obtain a first comparison result; and a second calculation unit, configured to calculate the filtering a gray value corresponding to the image component in the L channel after the uneven light intensity component; a second comparing unit, configured to compare the gray value with a predetermined gray threshold to obtain a second comparison result; a third calculating unit, configured to calculate a first intersection of the first comparison result and the second comparison result; and the iris extraction unit is configured to perform iris extraction according to the first intersection.
  • the iris extraction device further includes: a fourth calculation unit, configured to calculate a brightness of any one of the image components in the L channel after filtering the uneven illumination intensity component, minus the a difference between the average values of the brightness of the predetermined number of points around the point; a third comparing unit for comparing the difference with a predetermined brightness difference threshold to obtain a third comparison result; the fourth calculating unit, Computing a second intersection of the first intersection and the third comparison result; and the iris extraction unit is configured to perform iris extraction according to the second intersection.
  • a fourth calculation unit configured to calculate a brightness of any one of the image components in the L channel after filtering the uneven illumination intensity component, minus the a difference between the average values of the brightness of the predetermined number of points around the point
  • a third comparing unit for comparing the difference with a predetermined brightness difference threshold to obtain a third comparison result
  • the fourth calculating unit Computing a second intersection of the first intersection and the third comparison result
  • the iris extraction unit is configured to perform iris extraction according to the second
  • the fourth calculating unit calculates the difference according to the following formula:
  • (x, y) represents the coordinates of any point in the image component in the L channel after filtering out the uneven illumination intensity component
  • f(x, y) represents the point of the coordinate (x, y) Brightness value
  • the iris extraction device further includes: a determining unit, configured to divide Determining whether each foreground pixel in the binary image corresponding to the second intersection has an adjacent foreground pixel, and if present, dividing the foreground pixel and the adjacent foreground pixel into the same target area a plurality of target regions are obtained thereby; a fifth calculating unit, configured to separately calculate a number of foreground pixel points owned by the plurality of target regions; and a determining unit configured to determine the target region with the largest number as an iris region; And the iris extraction unit performs iris extraction in the iris region.
  • the iris extraction device further includes: an elimination unit for eliminating holes in the iris region, non-adjacent boundaries, to obtain a complete iris region; and the iris extraction unit at the complete Iris extraction is performed in the iris region.
  • the elimination unit eliminates holes, non-adjacent boundaries in the iris region according to the following formula:
  • A represents the iris region
  • represents any point in A
  • B represents a collection of structural elements
  • b represents any point in B. Represents a two-dimensional integer grid.
  • the elimination unit also eliminates holes, non-adjacent boundaries in the iris region according to the following formula:
  • A represents the iris region
  • represents any point in A
  • B represents a collection of structural elements. Represents a two-dimensional integer grid.
  • the extraction is gradually performed to avoid the influence of the uneven illumination condition on the iris extraction.
  • FIG. 1 illustrates a flow chart of an iris extraction method under uneven illumination conditions, in accordance with an embodiment of the present invention
  • FIG. 2 illustrates a schematic diagram of raw RGB before color mode conversion, in accordance with an embodiment of the present invention
  • FIG. 3 illustrates a schematic diagram of a color mode converted L channel in accordance with an embodiment of the present invention
  • FIG. 4 illustrates a schematic diagram of a channel after color mode conversion, in accordance with an embodiment of the present invention
  • FIG. 5 illustrates a schematic diagram of a b-channel after color mode conversion in accordance with an embodiment of the present invention
  • FIG. 6 illustrates a schematic diagram of chromatic aberrations of a channel and b channel in accordance with an embodiment of the present invention
  • FIG. 7 illustrates a schematic diagram of a left eye region color difference effect according to an embodiment of the present invention
  • FIG. 8 illustrates a schematic diagram of a first comparison result obtained by comparing the color difference in FIG. 7 with a predetermined color difference threshold according to an embodiment of the present invention
  • FIG. 9 illustrates a schematic diagram of filtering image components in an L channel after uneven illumination, in accordance with an embodiment of the present invention.
  • FIG. 10 illustrates a schematic diagram of a second comparison result obtained by comparing the grayscale value corresponding to FIG. 9 with a predetermined grayscale threshold value according to an embodiment of the present invention
  • FIG. 11 illustrates a schematic diagram of performing an operation of removing a local mean for the gray value corresponding to FIG. 9 according to an embodiment of the present invention
  • FIG. 12 is a diagram illustrating a third comparison result obtained by comparing a luminance difference value after the partial mean subtraction operation in FIG. 11 with a predetermined luminance difference threshold value according to an embodiment of the present invention
  • Figure 13 illustrates a schematic diagram of a first intersection of the first comparison result in Figure 8 and the second comparison result in Figure 10, in accordance with an embodiment of the present invention
  • FIG. 14 illustrates a schematic diagram of a second intersection of the first intersection of FIG. 13 and the third comparison result of FIG. 12, in accordance with an embodiment of the present invention
  • FIG. 15 illustrates a schematic diagram of a plurality of target regions obtained by dividing each foreground pixel into different target regions, respectively, according to an embodiment of the present invention
  • FIG. 16 illustrates a schematic diagram of the iris region having the largest number of foreground pixel points in FIG. 15 in accordance with an embodiment of the present invention
  • FIG. 17 illustrates a schematic view of a complete iris region obtained by eliminating holes, non-adjacent boundaries in an iris region, in accordance with an embodiment of the present invention
  • FIG. 18 illustrates a structural block diagram of an iris extraction apparatus under uneven illumination conditions according to an embodiment of the present invention.
  • FIG. 1 illustrates a flow chart of an iris extraction method under uneven illumination conditions, including steps S102 through S106 as follows, in accordance with an embodiment of the present invention.
  • Step S102 filtering out the uneven light intensity component from the image containing the iris.
  • step S104 the color difference is calculated in the filtered image.
  • step S106 iris extraction is performed according to the color difference.
  • each point in the image has a non-uniform illumination intensity component.
  • the influence of the uneven illumination condition on the iris extraction can be avoided.
  • embodiments of the present invention are capable of improving non-uniform illumination between regions within an image.
  • the luminance difference between the image frames is often used, for example, the luminances in the respective image frames are each extracted and the average value is calculated, and then the average value is used as the average luminance of each image frame.
  • the difference in luminance between image frames can be improved, it is not possible to improve non-uniform illumination between regions within the image.
  • filtering the uneven light intensity component in step S102 can be This is achieved by converting the color mode of the image, for example, by converting the uneven illumination intensity component and the color component in the converted image, thereby facilitating filtering out the uneven illumination intensity component.
  • the image can be converted from RGB color mode to Lab color mode (Lab is often used as an informal abbreviation for CIE 1976 (L*, a*, b*) color mode, ISO 11664-4: 2008 (E)/CIE S 014-4/E: 2007). Since the Lab color mode only has the uneven illumination intensity component in the image component of its L channel, the separation of the uneven illumination intensity component and the color component in the converted image is achieved.
  • iris extraction can be performed using the image component in the a channel of the Lab color mode and the color difference of the image component in the b channel. That is, the color difference is compared with a predetermined color difference threshold to obtain a first comparison result; and iris extraction is performed according to the first comparison result.
  • the color difference can be calculated according to the following formula and the extraction effect thereof can be illustrated in FIG. 6.
  • ⁇ E ab represents the color difference
  • ⁇ (a) represents the difference between the two values in the a channel
  • ⁇ (b) represents the difference between the two values in the b channel.
  • the unevenness can also be filtered out from the image components in the L channel.
  • the uniform light intensity component is calculated, and the gray value corresponding to the image component in the L channel after filtering the uneven light intensity component is calculated, and the iris extraction is performed more accurately according to the gray value and the color difference.
  • performing more precise iris extraction based on the gray value and the color difference may be achieved by comparing the gray value with a predetermined gray threshold to obtain a second comparison result; And comparing the result of the comparison with the first comparison result; and extracting the iris according to the intersection.
  • the joint extraction effect according to this embodiment can be referred to FIG. As is apparent from the image illustrated in Fig. 13, most of the iris regions have been whitened and can be clearly distinguished except for the eyeglass frame.
  • filtering the uneven illumination intensity component from the image components in the L channel can perform homomorphic filtering on the image components in the L channel, that is, performing logarithmic operations, fast Fourier transforms, High-pass filtering, inverse transform of fast Fourier transform, and exponential operation. Since the relative variation of the uneven illumination intensity in the image is small, it can be regarded as a low frequency component of the image, so that the uneven illumination intensity component can be filtered out by high-pass filtering.
  • the operation of removing the local mean may be performed, that is, calculating the brightness of any point in the image component in the L channel after filtering out the uneven illumination minus the reservation around the point.
  • the difference between the average values of the brightness of the number of points is calculated according to the following formula: Where (x, y) represents the coordinates of any point in the image component in the L channel after filtering out the uneven illumination intensity component, and f(x, y) represents the point of the coordinate (x, y) Brightness value, Represents a non-negative integer set and n represents a non-negative integer. Comparing the difference with a predetermined brightness difference threshold to obtain a third comparison result; calculating a second intersection of the first intersection and the third comparison result; and performing iris extraction according to the second intersection. In this way, the influence of various objects including the eyeglass frame, such as eyelashes, on iris extraction can be removed.
  • FIGS. 7 to 14 only take the right eye as an example, and the case of the left eye is similar to the right eye.
  • FIG. 7 illustrates a schematic diagram of chromatic aberration between image components in a channel and image components in a b channel, which has been described in detail above, and thus will not be described again herein, in accordance with an embodiment of the present invention.
  • FIG. 8 illustrates a schematic diagram of a first comparison result obtained by comparing the color difference in FIG. 7 with a predetermined color difference threshold, in which there is a spectacle frame and an enlarged iris edge region, in accordance with an embodiment of the present invention.
  • FIG. 9 illustrates a schematic diagram of filtering image components in an L channel after uneven illumination according to an embodiment of the present invention, and how to acquire image components in the L channel has been described in detail above, and thus will not be described herein.
  • FIG. 10 illustrates a schematic diagram of a second comparison result obtained by comparing the grayscale value corresponding to FIG. 9 with a predetermined grayscale threshold value, in accordance with an embodiment of the present invention.
  • FIG. 11 is a schematic diagram showing an operation of removing a local mean for the gray value corresponding to FIG. 9 according to an embodiment of the present invention, and how to remove the local portion has been described in detail above. Mean, so I won't go into details here.
  • FIG. 12 is a diagram illustrating a third comparison result obtained by comparing the luminance difference value after the partial mean subtraction operation in FIG. 11 with a predetermined luminance difference threshold value, in which the clearing can be clearly performed, according to an embodiment of the present invention. See the position of the eyeglass frame.
  • FIG. 13 illustrates a schematic diagram of a first intersection of the first comparison result in FIG. 8 and the second comparison result in FIG. 10, wherein the second comparison result is as shown in FIG. 9 according to an embodiment of the present invention.
  • the image component in the L channel is obtained by filtering the grayscale value corresponding to the uneven illumination intensity component and comparing it with a predetermined grayscale threshold. In Fig. 13, except for the eyeglass frame, there are almost no other unrelated areas.
  • FIG. 14 illustrates a second intersection of the first intersection in FIG. 13 and the third comparison result in FIG. 12 (ie, the intersection of the first comparison result, the second comparison result, and the third comparison result, according to an embodiment of the present invention).
  • the schematic diagram of the image mainly contains the iris area and few remaining fragments of the spectacle frame.
  • performing iris extraction according to the second intersection may be implemented as follows: respectively determining whether each foreground pixel in the binary image corresponding to the second intersection has an adjacent foreground pixel, and if present, the foreground pixel a point is divided into the adjacent target pixel in the same target area to obtain a plurality of target areas; respectively calculating a number of foreground pixel points owned by the plurality of target areas; determining the most-numbered target area as an iris area; Iris extraction is performed in the iris region.
  • FIG. 15 illustrates a schematic diagram of a plurality of target regions obtained by dividing respective foreground pixel points into different target regions, respectively, according to an embodiment of the present invention.
  • the determining step may be implemented as follows: whether each foreground pixel in the binary image is determined one above the other, below, below, left, right, upper left, upper right, lower left, and/or lower right. There is an adjacent foreground pixel, if present, dividing the foreground pixel into the same target area in the same target area, for example, using a uniform connectivity flag for each target area (eg, as an example in FIG. 15 The numbers 1, 2, 3) are used; if not present, the foreground pixel is the separated foreground pixel and is typically the residual background area in the non-iris region.
  • FIG. 16 illustrates the owned foreground pixels of FIG. 15 in accordance with an embodiment of the present invention.
  • holes in the iris region, non-adjacent boundaries can also be eliminated to obtain a complete iris region; and iris extraction is performed in the intact iris region.
  • the holes in the iris area and the non-adjacent boundaries are eliminated according to the following formula:
  • A represents the iris region
  • represents any point in A
  • B represents a collection of structural elements
  • b represents any point in B. Represents a two-dimensional integer grid.
  • A represents the iris region and B represents a set of structural elements with a radius R (for example, a circle with a radius of 3)
  • the above formula can be visually understood as sliding the center of B along the edge point of A to obtain an overlapping region;
  • the region performs a union operation with B to obtain an iris region (also referred to as an expanded iris region) whose edges are enlarged, holes, and non-adjacent boundaries are eliminated.
  • the center of B is slid along the edge point of the expanded iris region to obtain another overlapping region; the other overlapping region is intersected with B to obtain the above-mentioned complete iris region (also referred to as corrosion).
  • the iris area the size of the iris area remains the same, but the holes, non-adjacent boundaries that may appear inside the iris area are eliminated.
  • Figure 17 illustrates a schematic diagram of a complete iris region resulting from the elimination of holes in the iris region, non-adjacent boundaries, in accordance with an embodiment of the present invention, whereby it can be seen that a complete iris region without any interference, holes, has been obtained.
  • FIG. 18 illustrates a block diagram of a structure of an iris extraction apparatus under non-uniform illumination conditions, including a filtering device 1802, a computing device 1804, and an iris extraction device 1806, in accordance with an embodiment of the present invention.
  • the filtering device 1802 is configured to filter out the uneven illumination intensity component from the image containing the iris
  • the computing device 1804 is configured to calculate the color difference in the filtered image of the filtering device 1802
  • the iris extraction device 1806 is configured to use the computing device.
  • the iris was extracted by the color difference calculated in 1804.
  • the filtering device 1802 includes a conversion unit for converting a color mode of the image such that the uneven illumination intensity component and the color component in the converted image are separated.
  • the conversion unit converts the image from the RGB color mode to the Lab color mode to obtain respective image components in the L channel, the a channel, and the b channel, wherein the image component in the L channel includes the unevenness
  • the illumination intensity component, the image component in the a channel and the image component in the b channel include the color component.
  • computing device 1804 includes a first computing unit for calculating a color difference using image components in the a channel and image components in the b channel.
  • the computing unit calculates the color difference using the image component in the a channel and the image component in the b channel according to the following formula: Where ⁇ E ab represents the chromatic aberration, ⁇ (a) represents the difference between two values in the a channel, and ⁇ (b) represents the difference between the two values in the b channel.
  • the iris extraction device 1806 includes: a first comparison unit configured to compare the color difference with a predetermined color difference threshold to obtain a first comparison result; and an iris extraction unit configured to perform the first comparison result according to the first comparison result Iris extraction.
  • the filtering device 1802 further includes a filtering unit for filtering out the uneven light intensity component from the image components in the L channel.
  • the filtering unit includes: a first computing module for connecting the L The image component in the track performs a logarithm operation; the second operation module is configured to perform fast Fourier transform on the log transformed image component; and the third operation module is configured to filter the fast Fourier transform by high-pass filtering a low frequency portion of the image component; a fourth operation module for performing an inverse transform of the high-pass filtered image component on the fast Fourier transform; and a fifth operation module for performing an exponential operation on the inverse transformed image component.
  • the iris extraction device 1806 further includes: a second calculation unit, configured to calculate a gray value corresponding to the image component in the L channel after filtering the uneven illumination intensity component; and a second comparison unit, configured to The gray value is compared with a predetermined gray threshold to obtain a second comparison result; a third calculating unit is configured to calculate a first intersection of the first comparison result and the second comparison result; and an iris extraction unit is configured to use the first Intersection for iris extraction.
  • a second calculation unit configured to calculate a gray value corresponding to the image component in the L channel after filtering the uneven illumination intensity component
  • a second comparison unit configured to The gray value is compared with a predetermined gray threshold to obtain a second comparison result
  • a third calculating unit is configured to calculate a first intersection of the first comparison result and the second comparison result
  • an iris extraction unit is configured to use the first Intersection for iris extraction.
  • the iris extraction device 1806 further includes: a fourth calculation unit configured to calculate a brightness of any one of the image components in the L channel after filtering the uneven illumination intensity component minus a predetermined number of the points around the point a difference between the average values of the brightness of the points; a third comparing unit for comparing the difference with a predetermined brightness difference threshold to obtain a third comparison result; and a fourth calculating unit for calculating the first intersection and the first a second intersection of the three comparison results; and an iris extraction unit for performing iris extraction based on the second intersection.
  • the fourth calculation unit calculates the difference according to the following formula:
  • (x, y) represents the coordinates of any point in the image component in the L channel after filtering out the uneven illumination intensity component
  • f(x, y) represents the luminance value of the point of the coordinate (x, y)
  • the iris extraction device 1806 further includes: a determining unit, configured to respectively determine whether each foreground pixel in the binary image corresponding to the second intersection has an adjacent foreground pixel, and if present, the foreground Pixel points are divided into the same target area in the same target area to obtain a plurality of target areas; and a fifth calculating unit is configured to separately calculate the number of foreground pixel points owned by the plurality of target areas; The element is used to determine the target area of the largest number as the iris area; and the iris extraction unit performs iris extraction in the iris area.
  • the iris extraction device 1806 further includes: an elimination unit for eliminating holes in the iris region, non-adjacent boundaries, to obtain a complete iris region; and an iris extraction unit for iris extraction in the intact iris region .
  • the elimination unit eliminates holes, non-adjacent boundaries in the iris region according to the following formula:
  • A represents the iris region
  • represents any point in A
  • B represents a collection of structural elements
  • b represents any point in B. Represents a two-dimensional integer grid.
  • the elimination unit also eliminates holes, non-adjacent boundaries in the iris region according to the following formula:
  • A represents the iris region
  • represents any point in A
  • B represents a collection of structural elements. Represents a two-dimensional integer grid.
  • the influence of the uneven illumination condition on the iris extraction can be avoided; the local mean calculation can be performed by removing the image component in the L channel. The position of the spectacle frame is obtained, thereby eliminating the effect of the spectacle frame on iris extraction.

Abstract

Disclosed are a method and device for extracting an iris image under the condition of non-uniform illumination. The method comprises: filtering a non-uniform illumination intensity component from an image containing an iris image (S102); calculating a colour difference in the filtered image (S104); and conducting iris image extraction according to the colour difference (S106). In the embodiments, by filtering a non-uniform illumination intensity component, the influence of the condition of non-uniform illumination on the iris image extraction can be avoided.

Description

不均匀光照条件下的虹膜提取方法及设备Iris extraction method and device under uneven illumination 技术领域Technical field
本发明的实施例一般涉及图像处理领域,更具体地,涉及一种不均匀光照条件下的虹膜提取方法及设备。Embodiments of the present invention generally relate to the field of image processing, and more particularly to an iris extraction method and apparatus under uneven illumination conditions.
背景技术Background technique
在当今信息化时代,如何准确鉴定一个人的身份,保护信息安全是一个必须解决的关键社会问题。为此,生物特征鉴别技术悄然兴起,并成为目前世界信息安全管理领域的前沿研究课题。In today's information age, how to accurately identify a person's identity and protect information security is a key social issue that must be resolved. To this end, biometric identification technology has quietly emerged and has become a frontier research topic in the world of information security management.
生物特征鉴别技术是指利用人体所固有的生理特征或行为特征来进行个人身份鉴定,其中虹膜身份识别技术是生物特征鉴别技术的一个重要分支。对于一个人来说,他眼珠的中心是黑色的瞳孔,而虹膜就是瞳孔外缘间的环形组织,其呈现出相互交错的类似于斑点、细丝、条纹、隐窝的纹理特征。虹膜在一个人的一生中几乎不会发生改变,不同人的虹膜是完全不一样的。虹膜身份识别技术是计算机图像处理技术和模式识别技术在个人身份识别领域的应用。由于其存在高稳定性和高准确性,同时还可以使人们摆脱记忆信用卡号、银行帐号、身份证号、网络登录号的繁琐,因此在银行、公安、机场、网络安全等行业领域得到广泛使用。Biometric identification technology refers to the use of the physiological characteristics or behavioral characteristics inherent in the human body for personal identification. The iris identification technology is an important branch of biometric identification technology. For one person, the center of his eye is a black pupil, and the iris is the annular tissue between the outer edges of the pupil, which presents interlaced texture features similar to spots, filaments, stripes, and crypts. The iris is hardly changed in one's life, and the iris of different people is completely different. Iris identification technology is the application of computer image processing technology and pattern recognition technology in the field of personal identification. Because of its high stability and high accuracy, it can also make people get rid of the cumbersome memory credit card number, bank account number, ID number, and network access number. Therefore, it is widely used in banking, public security, airport, network security and other industries. .
典型的虹膜身份识别主要由图像采集、图像预处理、特征编码、特征匹配四部分构成。在图像预处理过程中,虹膜提取尤其关键,其执行时间和精度将直接影响整个虹膜身份识别过程的识别速度和精度。传统的虹膜提取方法包括Hough变换圆检测的方法、Daugman的基于微分积分算子的方法以及Wildes等采用的边界检测方法。他们都将虹膜边界建模为包括内外两条边界的圆环形,通过提取圆环的内外圆来实现虹膜提取。The typical iris identification is mainly composed of image acquisition, image preprocessing, feature coding and feature matching. In the image preprocessing process, iris extraction is especially critical, and its execution time and accuracy will directly affect the recognition speed and accuracy of the entire iris identification process. Traditional iris extraction methods include Hough transform circle detection method, Daugman's method based on differential integral operator, and boundary detection method adopted by Wildes et al. They both model the iris boundary as a circular ring with two inner and outer boundaries, and extract the iris by extracting the inner and outer circles of the ring.
然而,在实际采集过程中,由于上眼睑总是遮挡虹膜的上半部分,因此虹膜通常不是一个完整的圆环形。以及,虹膜外边界常常 较为模糊,即虹膜与巩膜(巩膜即眼球外围的白色部分,是眼睛最外层的纤维膜)的过渡不明显,从而难以用上述传统的边界检测方法来检测外边界。尤其需要指出的是,传统方法费时费力,难以满足实时快速准确检测的要求。However, during the actual acquisition process, the iris is usually not a complete circular ring since the upper eyelid always blocks the upper half of the iris. And the outer boundary of the iris is often It is more ambiguous, that is, the transition between the iris and the sclera (the white part of the sclera, that is, the outer part of the eyeball, which is the outermost layer of the eye) is not obvious, so it is difficult to detect the outer boundary by the above-mentioned conventional boundary detection method. In particular, it should be pointed out that the traditional method is time-consuming and laborious, and it is difficult to meet the requirements of real-time fast and accurate detection.
为此,对于这种非理想采集的虹膜提取,许多研究者也开展了相关工作。例如,在发明专利申请CN103246871A中,发明人首先利用现有Daugman方法来提取图像的内边界,进而选定内边界左右两侧的矩形区域,使得矩形区域尽可能覆盖内边界之外的左右两侧的虹膜,然后对该矩形区域进行灰度变换,利用Canny算子对矩形区域进行边界检测。具体地,在对矩形区域进行边界检测时,根据预定的阈值,低于阈值的像素值被设定为0,否则被设定为1,从而得到边界点集合。由此可见,在该发明专利申请CN103246871A中采用了所谓阈值分割的常用图像处理基本操作手段,用以提取出感兴趣的目标。To this end, many researchers have also carried out related work on this non-ideal collection of iris extraction. For example, in the invention patent application CN103246871A, the inventor first uses the existing Daugman method to extract the inner boundary of the image, and then selects the rectangular regions on the left and right sides of the inner boundary such that the rectangular region covers the left and right sides outside the inner boundary as much as possible. The iris is then gray-scale transformed into the rectangular region, and the Canny operator is used to detect the boundary of the rectangular region. Specifically, when performing boundary detection on a rectangular area, a pixel value lower than the threshold is set to 0 according to a predetermined threshold, and otherwise set to 1, thereby obtaining a boundary point set. It can be seen that in the invention patent application CN103246871A, a common image processing basic operation means called threshold segmentation is used to extract the target of interest.
然而,仅仅采用单一阈值操作通常难以稳健地检测出虹膜边界。特别是对于双眼虹膜身份识别的情形,由于双眼通常处于不均匀光照条件下,这就使得一个单一的、固定的阈值无法同时适应双眼。However, it is often difficult to robustly detect iris boundaries using only a single threshold operation. Especially for the case of binocular iris recognition, since the eyes are usually in uneven illumination, this makes a single, fixed threshold not adaptable to both eyes at the same time.
发明内容Summary of the invention
鉴于上述不均匀光照条件影响虹膜提取的问题,本发明的实施例提出一种不均匀光照条件下的虹膜提取方法及设备。In view of the above problem that uneven illumination conditions affect iris extraction, embodiments of the present invention provide an iris extraction method and apparatus under uneven illumination conditions.
根据本发明的一个方面,提供了一种不均匀光照条件下的虹膜提取方法,包括:从包含虹膜的图像中滤除不均匀光照强度分量;在滤除后的所述图像中计算色差;以及根据所述色差进行虹膜提取。According to an aspect of the present invention, there is provided an iris extraction method under uneven illumination conditions, comprising: filtering out an uneven illumination intensity component from an image containing an iris; calculating a color difference in the filtered image; Iris extraction is performed according to the color difference.
在一个实施例中,所述从包含虹膜的图像中滤除不均匀光照强度分量包括:对所述图像的色彩模式进行转换,使得转换后的所述图像中的所述不均匀光照强度分量和色彩分量相分离。In one embodiment, the filtering out the uneven illumination intensity component from the image containing the iris comprises: converting a color mode of the image such that the uneven illumination intensity component in the converted image and The color components are separated.
在一个实施例中,将所述图像从RGB色彩模式转换到Lab色彩模式,得到L通道、a通道和b通道中的相应的图像分量,以使得 转换后的所述图像中的所述不均匀光照强度分量和色彩分量相分离,其中所述L通道中的图像分量包括所述不均匀光照强度分量,所述a通道中的图像分量和所述b通道中的图像分量包括所述色彩分量。In one embodiment, the image is converted from an RGB color mode to a Lab color mode to obtain corresponding image components in the L channel, the a channel, and the b channel such that The uneven illumination intensity component and the color component in the converted image are separated, wherein the image component in the L channel includes the uneven illumination intensity component, an image component in the a channel, and the The image component in the b channel includes the color component.
在一个实施例中,所述在滤除后的所述图像中计算色差包括:利用所述a通道中的图像分量和所述b通道中的图像分量计算所述色差。In one embodiment, the calculating the color difference in the filtered image comprises calculating the color difference using an image component in the a channel and an image component in the b channel.
在一个实施例中,根据如下公式,利用所述a通道中的图像分量和所述b通道中的图像分量计算所述色差:
Figure PCTCN2015070058-appb-000001
其中ΔEab表示所述色差,Δ(a)表示所述a通道内两个数值之差,Δ(b)表示所述b通道内两个数值之差。
In one embodiment, the color difference is calculated using image components in the a channel and image components in the b channel according to the following formula:
Figure PCTCN2015070058-appb-000001
Where ΔE ab represents the chromatic aberration, Δ(a) represents the difference between the two values in the a channel, and Δ(b) represents the difference between the two values in the b channel.
在一个实施例中,所述根据所述色差进行虹膜提取包括:将所述色差与预定的色差阈值进行比较,得到第一比较结果;以及根据所述第一比较结果进行虹膜提取。In one embodiment, the performing iris extraction according to the color difference comprises: comparing the color difference with a predetermined color difference threshold to obtain a first comparison result; and performing iris extraction according to the first comparison result.
在一个实施例中,所述从包含虹膜的图像中滤除不均匀光照强度分量还包括:从所述L通道中的图像分量中滤除所述不均匀光照强度分量。In one embodiment, the filtering out the uneven illumination intensity component from the image comprising the iris further comprises filtering the non-uniform illumination intensity component from the image components in the L channel.
在一个实施例中,从所述L通道中的图像分量中滤除所述不均匀光照强度分量包括:将所述L通道中的图像分量进行对数运算;将对数变换后的所述图像分量进行快速傅里叶变换;通过高通滤波滤除快速傅里叶变换后的所述图像分量中的低频部分;将高通滤波后的所述图像分量进行快速傅里叶变换的逆变换;以及将逆变换后的所述图像分量进行指数运算。In one embodiment, filtering the uneven illumination intensity component from the image components in the L channel comprises: performing a logarithmic operation on the image components in the L channel; the log transformed image Performing a fast Fourier transform on the component; filtering the low frequency portion of the image component after the fast Fourier transform by high-pass filtering; performing inverse transform of the high-pass filtered image component on the fast Fourier transform; The inversely transformed image component is subjected to an exponential operation.
在一个实施例中,所述根据所述色差进行虹膜提取还包括:将所述色差与预定的色差阈值进行比较,得到第一比较结果;计算滤除所述不均匀光照强度分量后的所述L通道中的图像分量对应的灰度值;将所述灰度值与预定的灰度阈值进行比较,得到第二比较结果;计算所述第一比较结果和所述第二比较结果的第一交集;以及根据所述第一交集进行虹膜提取。 In one embodiment, the performing iris extraction according to the color difference further comprises: comparing the color difference with a predetermined color difference threshold to obtain a first comparison result; and calculating the filtering after filtering the uneven light intensity component a gray value corresponding to the image component in the L channel; comparing the gray value with a predetermined gray threshold to obtain a second comparison result; calculating the first comparison result and the first of the second comparison result Intersection; and performing iris extraction based on the first intersection.
在一个实施例中,所述根据所述色差进行虹膜提取还包括:计算滤除所述不均匀光照强度分量后的所述L通道中的图像分量中任意一点的亮度减去所述点周围的预定数目的点的亮度的平均值的差值;将所述差值与预定的亮度差值阈值进行比较,得到第三比较结果;计算所述第一交集和所述第三比较结果的第二交集;以及根据所述第二交集进行虹膜提取。In one embodiment, the performing iris extraction according to the color difference further comprises: calculating a brightness of any point in the image component in the L channel after filtering the uneven light intensity component minus a point around the point a difference between the average values of the brightness of the predetermined number of points; comparing the difference with a predetermined brightness difference threshold to obtain a third comparison result; calculating the first intersection and the second of the third comparison result Intersection; and performing iris extraction based on the second intersection.
在一个实施例中,根据如下公式计算所述差值:
Figure PCTCN2015070058-appb-000002
其中,(x,y)表示滤除所述不均匀光照强度分量后的所述L通道中的图像分量中的任意一点的坐标,f(x,y)表示坐标(x,y)的点的亮度值,
Figure PCTCN2015070058-appb-000003
表示非负整数集合,n表示非负整数。
In one embodiment, the difference is calculated according to the following formula:
Figure PCTCN2015070058-appb-000002
Where (x, y) represents the coordinates of any point in the image component in the L channel after filtering out the uneven illumination intensity component, and f(x, y) represents the point of the coordinate (x, y) Brightness value,
Figure PCTCN2015070058-appb-000003
Represents a non-negative integer set and n represents a non-negative integer.
在一个实施例中,根据所述第二交集进行虹膜提取包括:分别判断所述第二交集对应的二值图像中的每一前景像素点是否存在毗邻的前景像素点,如果存在,则将所述前景像素点与所述毗邻的前景像素点划分在同一目标区域内从而得到多个目标区域;分别计算所述多个目标区域所拥有的前景像素点的数目;将所述数目最多的目标区域确定为虹膜区域;以及在所述虹膜区域中进行虹膜提取。In an embodiment, performing iris extraction according to the second intersection includes: determining, respectively, whether each foreground pixel in the binary image corresponding to the second intersection has an adjacent foreground pixel, and if so, Deriving a foreground pixel in the same target area as the adjacent foreground pixel to obtain a plurality of target areas; respectively calculating a number of foreground pixel points owned by the plurality of target areas; and the target area having the largest number Determined to be an iris region; and perform iris extraction in the iris region.
在一个实施例中,在所述虹膜区域中进行虹膜提取包括:消除所述虹膜区域中的孔洞、不毗邻的边界,得到完整的虹膜区域;以及在所述完整的虹膜区域中进行虹膜提取。In one embodiment, performing iris extraction in the iris region includes: eliminating holes in the iris region, non-adjacent boundaries, resulting in a complete iris region; and performing iris extraction in the intact iris region.
在一个实施例中,根据如下公式消除所述虹膜区域中的孔洞、不毗邻的边界:In one embodiment, the holes in the iris region, non-adjacent boundaries are eliminated according to the following formula:
Figure PCTCN2015070058-appb-000004
Figure PCTCN2015070058-appb-000004
Figure PCTCN2015070058-appb-000005
Figure PCTCN2015070058-appb-000005
其中A表示所述虹膜区域,α表示A中的任意一点,B表示结构元素集合,b表示B中的任意一点,
Figure PCTCN2015070058-appb-000006
表示二维整数网格。
Where A represents the iris region, α represents any point in A, B represents a collection of structural elements, and b represents any point in B.
Figure PCTCN2015070058-appb-000006
Represents a two-dimensional integer grid.
在一个实施例中,还根据如下公式消除虹膜区域中的孔洞、不 毗邻的边界根据:In one embodiment, the holes in the iris region are also eliminated according to the following formula, Adjacent borders are based on:
Figure PCTCN2015070058-appb-000007
Figure PCTCN2015070058-appb-000007
其中A表示所述虹膜区域,α表示A中的任意一点,B表示结构元素集合,
Figure PCTCN2015070058-appb-000008
表示二维整数网格。
Where A represents the iris region, α represents any point in A, and B represents a collection of structural elements.
Figure PCTCN2015070058-appb-000008
Represents a two-dimensional integer grid.
根据本发明的另一个方面,提供了一种不均匀光照条件下的虹膜提取设备,包括:滤除装置,用于从包含虹膜的图像中滤除不均匀光照强度分量;计算装置,用于在滤除后的所述图像中计算色差;以及虹膜提取装置,用于根据所述色差进行虹膜提取。According to another aspect of the present invention, there is provided an iris extraction apparatus under uneven illumination conditions, comprising: a filtering device for filtering out uneven illumination intensity components from an image containing an iris; and a computing device for Calculating the color difference in the filtered image; and an iris extracting means for performing iris extraction based on the color difference.
在一个实施例中,所述滤除装置包括:转换单元,用于对所述图像的色彩模式进行转换,使得转换后的所述图像中的所述不均匀光照强度分量和色彩分量相分离。In one embodiment, the filtering device includes a conversion unit for converting a color mode of the image such that the uneven illumination intensity component and the color component in the converted image are separated.
在一个实施例中,所述转换单元将所述图像从RGB色彩模式转换到Lab色彩模式,得到L通道、a通道和b通道中的相应的图像分量,以使得转换后的所述图像中的所述不均匀光照强度分量和色彩分量相分离,其中所述L通道中的图像分量包括所述不均匀光照强度分量,所述a通道中的图像分量和所述b通道中的图像分量包括所述色彩分量。In one embodiment, the conversion unit converts the image from an RGB color mode to a Lab color mode to obtain corresponding image components in the L channel, the a channel, and the b channel, such that in the converted image The uneven illumination intensity component and the color component are separated, wherein the image component in the L channel includes the uneven illumination intensity component, and the image component in the a channel and the image component in the b channel include The color component.
在一个实施例中,所述计算装置包括:第一计算单元,用于利用所述a通道中的图像分量和所述b通道中的图像分量计算所述色差。In one embodiment, the computing device includes a first computing unit for calculating the color difference using image components in the a channel and image components in the b channel.
在一个实施例中,所述第一计算单元根据如下公式,利用所述a通道中的图像分量和所述b通道中的图像分量计算所述色差:
Figure PCTCN2015070058-appb-000009
其中ΔEab表示所述色差,Δ(a)表示所述a通道内两个数值之差,Δ(b)表示所述b通道内两个数值之差。
In one embodiment, the first calculating unit calculates the color difference by using an image component in the a channel and an image component in the b channel according to the following formula:
Figure PCTCN2015070058-appb-000009
Where ΔE ab represents the chromatic aberration, Δ(a) represents the difference between the two values in the a channel, and Δ(b) represents the difference between the two values in the b channel.
在一个实施例中,所述虹膜提取装置包括:第一比较单元,用于将所述色差与预定的色差阈值进行比较,得到第一比较结果;以及虹膜提取单元,用于根据所述第一比较结果进行虹膜提取。In one embodiment, the iris extraction device includes: a first comparison unit configured to compare the color difference with a predetermined color difference threshold to obtain a first comparison result; and an iris extraction unit configured to The results were compared for iris extraction.
在一个实施例中,所述滤除装置还包括:滤除单元,用于从所述L通道中的图像分量中滤除所述不均匀光照强度分量。 In one embodiment, the filtering device further comprises: a filtering unit for filtering out the uneven illumination intensity component from image components in the L channel.
在一个实施例中,所述滤除单元包括:第一运算模块,用于将所述L通道中的图像分量进行对数运算;第二运算模块,用于将对数变换后的所述图像分量进行快速傅里叶变换;第三运算模块,用于通过高通滤波滤除快速傅里叶变换后的所述图像分量中的低频部分;第四运算模块,用于将高通滤波后的所述图像分量进行快速傅里叶变换的逆变换;以及第五运算模块,用于将逆变换后的所述图像分量进行指数运算。In one embodiment, the filtering unit includes: a first operation module, configured to perform a logarithm operation on image components in the L channel; and a second operation module, configured to convert the logarithmically transformed image The component performs a fast Fourier transform; a third operation module is configured to filter the low frequency portion of the image component after the fast Fourier transform by high-pass filtering; and a fourth operation module, configured to filter the high-pass The image component performs an inverse transform of the fast Fourier transform; and a fifth operation module is configured to perform an exponential operation on the inverse transformed image component.
在一个实施例中,所述虹膜提取装置包括:第一比较单元,用于将所述色差与预定的色差阈值进行比较,得到第一比较结果;第二计算单元,用于计算滤除所述不均匀光照强度分量后的所述L通道中的图像分量对应的灰度值;第二比较单元,用于将所述灰度值与预定的灰度阈值进行比较,得到第二比较结果;第三计算单元,用于计算所述第一比较结果和所述第二比较结果的第一交集;以及所述虹膜提取单元用于根据所述第一交集进行虹膜提取。In one embodiment, the iris extraction device includes: a first comparison unit, configured to compare the color difference with a predetermined color difference threshold to obtain a first comparison result; and a second calculation unit, configured to calculate the filtering a gray value corresponding to the image component in the L channel after the uneven light intensity component; a second comparing unit, configured to compare the gray value with a predetermined gray threshold to obtain a second comparison result; a third calculating unit, configured to calculate a first intersection of the first comparison result and the second comparison result; and the iris extraction unit is configured to perform iris extraction according to the first intersection.
在一个实施例中,所述虹膜提取装置还包括:第四计算单元,用于计算滤除所述不均匀光照强度分量后的所述L通道中的图像分量中任意一点的亮度减去所述点周围的预定数目的点的亮度的平均值的差值;第三比较单元,用于将所述差值与预定的亮度差值阈值进行比较,得到第三比较结果;第四计算单元,用于计算所述第一交集和所述第三比较结果的第二交集;以及所述虹膜提取单元用于根据所述第二交集进行虹膜提取。In one embodiment, the iris extraction device further includes: a fourth calculation unit, configured to calculate a brightness of any one of the image components in the L channel after filtering the uneven illumination intensity component, minus the a difference between the average values of the brightness of the predetermined number of points around the point; a third comparing unit for comparing the difference with a predetermined brightness difference threshold to obtain a third comparison result; the fourth calculating unit, Computing a second intersection of the first intersection and the third comparison result; and the iris extraction unit is configured to perform iris extraction according to the second intersection.
在一个实施例中,所述第四计算单元根据如下公式计算所述差值:In one embodiment, the fourth calculating unit calculates the difference according to the following formula:
Figure PCTCN2015070058-appb-000010
Figure PCTCN2015070058-appb-000010
其中,(x,y)表示滤除所述不均匀光照强度分量后的所述L通道中的图像分量中的任意一点的坐标,f(x,y)表示坐标(x,y)的点的亮度值,
Figure PCTCN2015070058-appb-000011
表示非负整数集合,n表示非负整数。
Where (x, y) represents the coordinates of any point in the image component in the L channel after filtering out the uneven illumination intensity component, and f(x, y) represents the point of the coordinate (x, y) Brightness value,
Figure PCTCN2015070058-appb-000011
Represents a non-negative integer set and n represents a non-negative integer.
在一个实施例中,所述虹膜提取装置还包括:判断单元,用于分 别判断所述第二交集对应的二值图像中的每一前景像素点是否存在毗邻的前景像素点,如果存在,则将所述前景像素点与所述毗邻的前景像素点划分在同一目标区域内从而得到多个目标区域;第五计算单元,用于分别计算所述多个目标区域所拥有的前景像素点的数目;确定单元,用于将所述数目最多的目标区域确定为虹膜区域;以及所述虹膜提取单元在所述虹膜区域中进行虹膜提取。In an embodiment, the iris extraction device further includes: a determining unit, configured to divide Determining whether each foreground pixel in the binary image corresponding to the second intersection has an adjacent foreground pixel, and if present, dividing the foreground pixel and the adjacent foreground pixel into the same target area a plurality of target regions are obtained thereby; a fifth calculating unit, configured to separately calculate a number of foreground pixel points owned by the plurality of target regions; and a determining unit configured to determine the target region with the largest number as an iris region; And the iris extraction unit performs iris extraction in the iris region.
在一个实施例中,所述虹膜提取装置还包括:消除单元,用于消除所述虹膜区域中的孔洞、不毗邻的边界,得到完整的虹膜区域;以及所述虹膜提取单元在所述完整的虹膜区域中进行虹膜提取。In one embodiment, the iris extraction device further includes: an elimination unit for eliminating holes in the iris region, non-adjacent boundaries, to obtain a complete iris region; and the iris extraction unit at the complete Iris extraction is performed in the iris region.
在一个实施例中,所述消除单元根据如下公式消除所述虹膜区域中的孔洞、不毗邻的边界:In one embodiment, the elimination unit eliminates holes, non-adjacent boundaries in the iris region according to the following formula:
Figure PCTCN2015070058-appb-000012
Figure PCTCN2015070058-appb-000012
Figure PCTCN2015070058-appb-000013
Figure PCTCN2015070058-appb-000013
其中A表示所述虹膜区域,α表示A中的任意一点,B表示结构元素集合,b表示B中的任意一点,
Figure PCTCN2015070058-appb-000014
表示二维整数网格。
Where A represents the iris region, α represents any point in A, B represents a collection of structural elements, and b represents any point in B.
Figure PCTCN2015070058-appb-000014
Represents a two-dimensional integer grid.
在一个实施例中,所述消除单元还根据如下公式消除虹膜区域中的孔洞、不毗邻的边界:In one embodiment, the elimination unit also eliminates holes, non-adjacent boundaries in the iris region according to the following formula:
Figure PCTCN2015070058-appb-000015
Figure PCTCN2015070058-appb-000015
其中A表示所述虹膜区域,α表示A中的任意一点,B表示结构元素集合,
Figure PCTCN2015070058-appb-000016
表示二维整数网格。
Where A represents the iris region, α represents any point in A, and B represents a collection of structural elements.
Figure PCTCN2015070058-appb-000016
Represents a two-dimensional integer grid.
通过下文详细描述将会理解,根据本发明的实施例,滤除不均匀光照强度分量后,逐步提取从而可以避免该不均匀光照条件对虹膜提取的影响。As will be understood from the detailed description below, according to an embodiment of the present invention, after filtering out the uneven illumination intensity component, the extraction is gradually performed to avoid the influence of the uneven illumination condition on the iris extraction.
附图说明DRAWINGS
通过参考附图阅读下文的详细描述,本发明的实施例的上述以及其它目的、特征和优点将变得易于理解。在附图中,以示例性而非限制性的方式示出了本发明的若干实施例,其中: The above and other objects, features and advantages of the embodiments of the present invention will become <RTIgt; In the drawings, several embodiments of the invention are illustrated in the
图1图示了根据本发明的实施例的不均匀光照条件下的虹膜提取方法的流程图;1 illustrates a flow chart of an iris extraction method under uneven illumination conditions, in accordance with an embodiment of the present invention;
图2图示了根据本发明的实施例的色彩模式转换前的原始RGB的示意图;2 illustrates a schematic diagram of raw RGB before color mode conversion, in accordance with an embodiment of the present invention;
图3图示了根据本发明的实施例的色彩模式转换后的L通道的示意图;3 illustrates a schematic diagram of a color mode converted L channel in accordance with an embodiment of the present invention;
图4图示了根据本发明的实施例的色彩模式转换后的a通道的示意图;4 illustrates a schematic diagram of a channel after color mode conversion, in accordance with an embodiment of the present invention;
图5图示了根据本发明的实施例的色彩模式转换后的b通道的示意图;FIG. 5 illustrates a schematic diagram of a b-channel after color mode conversion in accordance with an embodiment of the present invention; FIG.
图6图示了根据本发明的实施例a通道和b通道的色差的示意图;6 illustrates a schematic diagram of chromatic aberrations of a channel and b channel in accordance with an embodiment of the present invention;
图7图示了根据本发明的实施例的左眼区域色差效果的示意图;FIG. 7 illustrates a schematic diagram of a left eye region color difference effect according to an embodiment of the present invention; FIG.
图8图示了根据本发明的实施例的将图7中的色差与预定的色差阈值进行比较所得到的第一比较结果的示意图;8 illustrates a schematic diagram of a first comparison result obtained by comparing the color difference in FIG. 7 with a predetermined color difference threshold according to an embodiment of the present invention;
图9图示了根据本发明的实施例的滤除不均匀光照后的L通道中的图像分量的示意图;9 illustrates a schematic diagram of filtering image components in an L channel after uneven illumination, in accordance with an embodiment of the present invention;
图10图示了根据本发明的实施例的将图9对应的灰度值与预定的灰度阈值进行比较所得到的第二比较结果的示意图;FIG. 10 illustrates a schematic diagram of a second comparison result obtained by comparing the grayscale value corresponding to FIG. 9 with a predetermined grayscale threshold value according to an embodiment of the present invention; FIG.
图11图示了根据本发明的实施例的对图9对应的灰度值执行去除局部均值的操作后的示意图;11 illustrates a schematic diagram of performing an operation of removing a local mean for the gray value corresponding to FIG. 9 according to an embodiment of the present invention;
图12图示了根据本发明的实施例的将图11中的去除局部均值的操作后的亮度差值与预定的亮度差值阈值进行比较所得到的第三比较结果的示意图;12 is a diagram illustrating a third comparison result obtained by comparing a luminance difference value after the partial mean subtraction operation in FIG. 11 with a predetermined luminance difference threshold value according to an embodiment of the present invention;
图13图示了根据本发明的实施例的图8中的第一比较结果和图10中的第二比较结果的第一交集的示意图;Figure 13 illustrates a schematic diagram of a first intersection of the first comparison result in Figure 8 and the second comparison result in Figure 10, in accordance with an embodiment of the present invention;
图14图示了根据本发明的实施例的图13中的第一交集和图12中的第三比较结果的第二交集的示意图; 14 illustrates a schematic diagram of a second intersection of the first intersection of FIG. 13 and the third comparison result of FIG. 12, in accordance with an embodiment of the present invention;
图15图示了根据本发明的实施例的将各前景像素点分别划分在不同目标区域所得到的多个目标区域的示意图;15 illustrates a schematic diagram of a plurality of target regions obtained by dividing each foreground pixel into different target regions, respectively, according to an embodiment of the present invention;
图16图示了根据本发明的实施例的图15中的拥有的前景像素点的数目最多的虹膜区域的示意图;16 illustrates a schematic diagram of the iris region having the largest number of foreground pixel points in FIG. 15 in accordance with an embodiment of the present invention;
图17图示了根据本发明的实施例的消除虹膜区域中的孔洞、不毗邻的边界所得到完整的虹膜区域的示意图;17 illustrates a schematic view of a complete iris region obtained by eliminating holes, non-adjacent boundaries in an iris region, in accordance with an embodiment of the present invention;
图18图示了根据本发明的实施例的不均匀光照条件下的虹膜提取设备的结构框图。FIG. 18 illustrates a structural block diagram of an iris extraction apparatus under uneven illumination conditions according to an embodiment of the present invention.
在各个附图中,相同或对应的标号表示相同或对应的部分。In the various figures, the same or corresponding reference numerals indicate the same or corresponding parts.
具体实施方式detailed description
下面将参考附图中示出的若干示例性实施例来描述本发明的原理和精神。应当理解,给出这些实施例仅仅是为了使本领域技术人员能够更好地理解进而实现本发明,而并非以任何方式限制本发明的范围。The principles and spirit of the present invention are described below with reference to a few exemplary embodiments illustrated in the drawings. It is to be understood that the examples are given only to enable those skilled in the art to understand the invention, and not to limit the scope of the invention in any way.
图1图示了根据本发明的实施例的不均匀光照条件下的虹膜提取方法的流程图,其包括如下的步骤S102至步骤S106。1 illustrates a flow chart of an iris extraction method under uneven illumination conditions, including steps S102 through S106 as follows, in accordance with an embodiment of the present invention.
步骤S102,从包含虹膜的图像中滤除不均匀光照强度分量。Step S102, filtering out the uneven light intensity component from the image containing the iris.
步骤S104,在滤除后的图像中计算色差。In step S104, the color difference is calculated in the filtered image.
步骤S106,根据该色差进行虹膜提取。In step S106, iris extraction is performed according to the color difference.
鉴于在不均匀光照条件下,图像中的各点具有不均匀光照强度分量,在本实施例中,通过滤除该不均匀光照强度分量,从而可以避免该不均匀光照条件对虹膜提取的影响。In view of the uneven illumination condition, each point in the image has a non-uniform illumination intensity component. In the present embodiment, by filtering the uneven illumination intensity component, the influence of the uneven illumination condition on the iris extraction can be avoided.
由此可见,本发明的实施例能够改善图像内部各区域之间的非均匀光照。然而,相关技术中往往针对各图像帧之间的亮度差异,例如将各图像帧中的亮度各自提取出来并计算平均值,然后以该平均值作为各图像帧的平均亮度。这样,虽然可以改善各图像帧之间的亮度差异,但是并不能够改善图像内部各区域之间的非均匀光照。Thus, embodiments of the present invention are capable of improving non-uniform illumination between regions within an image. However, in the related art, the luminance difference between the image frames is often used, for example, the luminances in the respective image frames are each extracted and the average value is calculated, and then the average value is used as the average luminance of each image frame. Thus, although the difference in luminance between image frames can be improved, it is not possible to improve non-uniform illumination between regions within the image.
在一个实施例中,步骤S102中的滤除该不均匀光照强度分量可 以经由转换该图像的色彩模式而实现,例如,可以通过转换而使得在转换后的图像中的不均匀光照强度分量和色彩分量相分离,从而有利于滤除该不均匀光照强度分量。In an embodiment, filtering the uneven light intensity component in step S102 can be This is achieved by converting the color mode of the image, for example, by converting the uneven illumination intensity component and the color component in the converted image, thereby facilitating filtering out the uneven illumination intensity component.
在这一实施例中,可以将该图像从RGB色彩模式转换到Lab色彩模式(Lab经常用做CIE 1976(L*,a*,b*)色彩模式的非正式缩写,ISO 11664-4:2008(E)/CIE S 014-4/E:2007)。由于Lab色彩模式仅仅在其L通道的图像分量中存在该不均匀光照强度分量,因此使得在转换后的图像中的不均匀光照强度分量和色彩分量相分离得以实现。In this embodiment, the image can be converted from RGB color mode to Lab color mode (Lab is often used as an informal abbreviation for CIE 1976 (L*, a*, b*) color mode, ISO 11664-4: 2008 (E)/CIE S 014-4/E: 2007). Since the Lab color mode only has the uneven illumination intensity component in the image component of its L channel, the separation of the uneven illumination intensity component and the color component in the converted image is achieved.
图2至图5分别图示了上述色彩模式转换前后的图像。其中,在图2所图示的转换前的原始RGB图像中,可以明显看出左眼和右眼处于不均匀光照条件下。在图3所图示的转换后的L通道的图像中,该左眼和右眼处的不均匀光照情况仍然存在。但是,在图4和图5所分别图示的转换后的a通道和b通道的图像中,在该左眼和右眼处已不存在不均匀光照情况,而仅仅存在色差情况。2 to 5 respectively illustrate images before and after the above color mode conversion. Among them, in the original RGB image before the conversion illustrated in FIG. 2, it can be clearly seen that the left eye and the right eye are in uneven illumination conditions. In the image of the converted L channel illustrated in FIG. 3, uneven illumination conditions at the left and right eyes still exist. However, in the images of the converted a-channel and b-channel illustrated in FIGS. 4 and 5, respectively, there is no uneven illumination at the left and right eyes, and only the chromatic aberration exists.
本领域技术人员应当知晓,上述色彩模式的转换仅仅是示例性而非限制性的,其他任何能够实现不均匀光照强度分量的分离的实施例、甚至不分离而直接能够滤除该不均匀光照强度分量的实施例均应当纳入本发明的保护范围。It will be appreciated by those skilled in the art that the above-described color mode conversion is merely exemplary and not limiting, and any other embodiment capable of achieving a separation of uneven illumination intensity components, even without separation, can directly filter out the uneven illumination intensity. Embodiments of the components should all be included in the scope of protection of the present invention.
在一个实施例中,由于在Lab色彩模式中虹膜和皮肤之间存在明显的颜色差异,因此可以利用Lab色彩模式的a通道中的图像分量和b通道中的图像分量的色差来进行虹膜提取,即将该色差与预定的色差阈值进行比较,得到第一比较结果;以及根据该第一比较结果进行虹膜提取。其中,该色差可以根据如下公式进行计算并且其提取效果可以如图6所图示。In one embodiment, since there is a significant color difference between the iris and the skin in the Lab color mode, iris extraction can be performed using the image component in the a channel of the Lab color mode and the color difference of the image component in the b channel. That is, the color difference is compared with a predetermined color difference threshold to obtain a first comparison result; and iris extraction is performed according to the first comparison result. Wherein, the color difference can be calculated according to the following formula and the extraction effect thereof can be illustrated in FIG. 6.
Figure PCTCN2015070058-appb-000017
Figure PCTCN2015070058-appb-000017
其中ΔEab表示色差,Δ(a)表示a通道内两个数值之差,Δ(b)表示b通道内两个数值之差。Where ΔE ab represents the color difference, Δ(a) represents the difference between the two values in the a channel, and Δ(b) represents the difference between the two values in the b channel.
在一个实施例中,还可以从L通道中的图像分量中滤除该不均 匀光照强度分量,计算滤除该不均匀光照强度分量后的L通道中的图像分量对应的灰度值,以及根据该灰度值和上述色差来进行更精确地虹膜提取。在这一实施例中,根据该灰度值和上述色差来进行更精确地虹膜提取可以实现如下:通过将该灰度值与预定的灰度阈值进行比较,得到第二比较结果;计算该第二比较结果和上述第一比较结果的交集;以及根据该交集进行虹膜提取。出于使得本发明实施例的上下文描述更加清晰、连贯的目的,根据这一实施例的联合提取效果可以去参考图13。在图13所图示的图像中可以明显看出,除去眼镜框之外,其余大部分虹膜区域都已经呈现白色并且能够明确区分。In one embodiment, the unevenness can also be filtered out from the image components in the L channel. The uniform light intensity component is calculated, and the gray value corresponding to the image component in the L channel after filtering the uneven light intensity component is calculated, and the iris extraction is performed more accurately according to the gray value and the color difference. In this embodiment, performing more precise iris extraction based on the gray value and the color difference may be achieved by comparing the gray value with a predetermined gray threshold to obtain a second comparison result; And comparing the result of the comparison with the first comparison result; and extracting the iris according to the intersection. For the purpose of making the context description of the embodiment of the present invention clearer and more consistent, the joint extraction effect according to this embodiment can be referred to FIG. As is apparent from the image illustrated in Fig. 13, most of the iris regions have been whitened and can be clearly distinguished except for the eyeglass frame.
在这一实施例中,从L通道中的图像分量中滤除该不均匀光照强度分量可以通过对L通道中的图像分量进行同态滤波,即先后进行对数运算、快速傅里叶变换、高通滤波、快速傅里叶变换的逆变换以及指数运算。由于图像中的不均匀光照强度相对变化很小,因此其可以看作是图像的低频成份,因此通过高通滤波该不均匀光照强度分量得以被滤除。In this embodiment, filtering the uneven illumination intensity component from the image components in the L channel can perform homomorphic filtering on the image components in the L channel, that is, performing logarithmic operations, fast Fourier transforms, High-pass filtering, inverse transform of fast Fourier transform, and exponential operation. Since the relative variation of the uneven illumination intensity in the image is small, it can be regarded as a low frequency component of the image, so that the uneven illumination intensity component can be filtered out by high-pass filtering.
本领域技术人员应当知晓,上述涉及同态滤波的操作仅仅是示例性而非限制性的,其他任何能够实现不均匀光照强度分量的滤除的实施例均应当纳入本发明的保护范围。Those skilled in the art will appreciate that the above-described operations involving homomorphic filtering are merely exemplary and not limiting, and any other embodiment that enables filtering of non-uniform illumination intensity components should be included in the scope of the present invention.
此外,在虹膜提取过程中,由于眼镜框与虹膜距离相近并且通常颜色近似,因此影响虹膜提取的效果。为了解决该问题,在本发明的一个实施例中,可以执行去除局部均值的操作,即计算滤除不均匀光照后的L通道中的图像分量中的任意一点的亮度减去该点周围的预定数目的点的亮度的平均值的差值,具体来说,根据如下公式计算该差值:
Figure PCTCN2015070058-appb-000018
其中,(x,y)表示滤除所述不均匀光照强度分量后的所述L通道中的图像分量中的任意一点的坐标,f(x,y)表示坐标(x,y)的点的亮度值,
Figure PCTCN2015070058-appb-000019
表示非负整 数集合,n表示非负整数。将该差值与预定的亮度差值阈值进行比较,得到第三比较结果;计算第一交集和第三比较结果的第二交集;以及根据第二交集进行虹膜提取。这样,可以去除包括眼镜框在内的各种物体(例如眼睫毛)对虹膜提取的影响。
In addition, in the iris extraction process, since the eyeglass frame is close to the iris and generally has a similar color, the effect of iris extraction is affected. In order to solve this problem, in one embodiment of the present invention, the operation of removing the local mean may be performed, that is, calculating the brightness of any point in the image component in the L channel after filtering out the uneven illumination minus the reservation around the point. The difference between the average values of the brightness of the number of points, specifically, the difference is calculated according to the following formula:
Figure PCTCN2015070058-appb-000018
Where (x, y) represents the coordinates of any point in the image component in the L channel after filtering out the uneven illumination intensity component, and f(x, y) represents the point of the coordinate (x, y) Brightness value,
Figure PCTCN2015070058-appb-000019
Represents a non-negative integer set and n represents a non-negative integer. Comparing the difference with a predetermined brightness difference threshold to obtain a third comparison result; calculating a second intersection of the first intersection and the third comparison result; and performing iris extraction according to the second intersection. In this way, the influence of various objects including the eyeglass frame, such as eyelashes, on iris extraction can be removed.
下面结合图7至图14对其进行详细说明。为了简化说明起见,图7至图14仅仅以右眼为例,左眼的情况与右眼类似。This will be described in detail below with reference to FIGS. 7 to 14. For the sake of simplicity of explanation, FIGS. 7 to 14 only take the right eye as an example, and the case of the left eye is similar to the right eye.
本领域技术人员应当理解,在执行去除局部均值的操作之前,可以针对不同图像而通过人工选择是否已经戴眼镜,如果选择结果指示已经戴眼镜则使用图7至图14所描述的涉及第一比较结果、第二比较结果和第三比较结果的虹膜提取方法,否则仅仅使用涉及第一比较结果和第二比较结果的虹膜提取方法。同时,本领域技术人员还应当理解,由于上述涉及第一比较结果、第二比较结果和第三比较结果的虹膜提取方法的计算速度很快,因此更加优选的是无需上述人工选择步骤而针对所有图像都执行上述涉及第一比较结果、第二比较结果和第三比较结果的虹膜提取方法。上述这些方法均应当纳入本发明的保护范围。It will be understood by those skilled in the art that before performing the operation of removing the local mean, it is possible to manually select whether or not the glasses have been worn for different images, and if the result of the selection indicates that the glasses have been worn, the first comparison described with reference to FIGS. 7 to 14 is used. The result, the second comparison result, and the iris extraction method of the third comparison result, otherwise only the iris extraction method involving the first comparison result and the second comparison result is used. Meanwhile, those skilled in the art should also understand that since the above-mentioned iris extraction method involving the first comparison result, the second comparison result, and the third comparison result is fast, it is more preferable that all of the manual selection steps are not required. The image performs the above-described iris extraction method involving the first comparison result, the second comparison result, and the third comparison result. All of the above methods should be included in the scope of protection of the present invention.
图7图示了根据本发明的实施例的a通道中的图像分量和b通道中的图像分量之间的色差的示意图,上文已经详细描述了如何获取该色差,因此这里不再赘述。7 illustrates a schematic diagram of chromatic aberration between image components in a channel and image components in a b channel, which has been described in detail above, and thus will not be described again herein, in accordance with an embodiment of the present invention.
图8图示了根据本发明的实施例的将图7中的色差与预定的色差阈值进行比较所得到的第一比较结果的示意图,其中存在眼镜框以及扩大的虹膜边缘区域。8 illustrates a schematic diagram of a first comparison result obtained by comparing the color difference in FIG. 7 with a predetermined color difference threshold, in which there is a spectacle frame and an enlarged iris edge region, in accordance with an embodiment of the present invention.
图9图示了根据本发明的实施例的滤除不均匀光照后的L通道中的图像分量的示意图,上文已经详细描述了如何获取该L通道中的图像分量,因此这里不再赘述。FIG. 9 illustrates a schematic diagram of filtering image components in an L channel after uneven illumination according to an embodiment of the present invention, and how to acquire image components in the L channel has been described in detail above, and thus will not be described herein.
图10图示了根据本发明的实施例的将图9对应的灰度值与预定的灰度阈值进行比较所得到的第二比较结果的示意图。FIG. 10 illustrates a schematic diagram of a second comparison result obtained by comparing the grayscale value corresponding to FIG. 9 with a predetermined grayscale threshold value, in accordance with an embodiment of the present invention.
图11图示了根据本发明的实施例的对图9对应的灰度值执行去除局部均值的操作后的示意图,上文已经详细描述了如何去除局部 均值,因此这里不再赘述。FIG. 11 is a schematic diagram showing an operation of removing a local mean for the gray value corresponding to FIG. 9 according to an embodiment of the present invention, and how to remove the local portion has been described in detail above. Mean, so I won't go into details here.
图12图示了根据本发明的实施例的将图11中的去除局部均值的操作后的亮度差值与预定的亮度差值阈值进行比较所得到的第三比较结果的示意图,其中可以清晰地看到眼镜框的位置。12 is a diagram illustrating a third comparison result obtained by comparing the luminance difference value after the partial mean subtraction operation in FIG. 11 with a predetermined luminance difference threshold value, in which the clearing can be clearly performed, according to an embodiment of the present invention. See the position of the eyeglass frame.
图13图示了根据本发明的实施例的图8中的第一比较结果和图10中的第二比较结果的第一交集的示意图,其中如上所述,第二比较结果是将图9中的L通道中的图像分量滤除不均匀光照强度分量后所对应的灰度值与预定的灰度阈值进行比较而得到的。图13中除了眼镜框,几乎不存在其他无关区域。13 illustrates a schematic diagram of a first intersection of the first comparison result in FIG. 8 and the second comparison result in FIG. 10, wherein the second comparison result is as shown in FIG. 9 according to an embodiment of the present invention. The image component in the L channel is obtained by filtering the grayscale value corresponding to the uneven illumination intensity component and comparing it with a predetermined grayscale threshold. In Fig. 13, except for the eyeglass frame, there are almost no other unrelated areas.
图14图示了根据本发明的实施例的图13中的第一交集和图12中的第三比较结果的第二交集(即第一比较结果、第二比较结果和第三比较结果的交集)的示意图,图中主要剩下虹膜区域和很少的眼镜架残存片段。14 illustrates a second intersection of the first intersection in FIG. 13 and the third comparison result in FIG. 12 (ie, the intersection of the first comparison result, the second comparison result, and the third comparison result, according to an embodiment of the present invention). The schematic diagram of the image mainly contains the iris area and few remaining fragments of the spectacle frame.
在一个实施例中,根据第二交集进行虹膜提取可以实现如下:分别判断第二交集对应的二值图像中的每一前景像素点是否存在毗邻的前景像素点,如果存在,则将该前景像素点与该毗邻的前景像素点划分在同一目标区域内从而得到多个目标区域;分别计算该多个目标区域所拥有的前景像素点的数目;将数目最多的目标区域确定为虹膜区域;以及在该虹膜区域中进行虹膜提取。In one embodiment, performing iris extraction according to the second intersection may be implemented as follows: respectively determining whether each foreground pixel in the binary image corresponding to the second intersection has an adjacent foreground pixel, and if present, the foreground pixel a point is divided into the adjacent target pixel in the same target area to obtain a plurality of target areas; respectively calculating a number of foreground pixel points owned by the plurality of target areas; determining the most-numbered target area as an iris area; Iris extraction is performed in the iris region.
图15图示了根据本发明的实施例的将各前景像素点分别划分在不同目标区域所得到的多个目标区域的示意图。基于图15,上述判断步骤可以实现如下:逐个判断该二值图像中的每一前景像素点在其上方、下方、左方、右方、左上方、右上方、左下方和/或右下方是否存在毗邻的前景像素点,如果存在,则将该前景像素点与该毗邻的前景像素点划分在同一目标区域内,例如,针对每个目标区域采用统一的连通标记(例如图15中作为示例而采用数字1、2、3);如果不存在,则该前景像素点就是分离的前景像素点,并且通常是非虹膜区域中的残存的背景区域。FIG. 15 illustrates a schematic diagram of a plurality of target regions obtained by dividing respective foreground pixel points into different target regions, respectively, according to an embodiment of the present invention. Based on FIG. 15, the determining step may be implemented as follows: whether each foreground pixel in the binary image is determined one above the other, below, below, left, right, upper left, upper right, lower left, and/or lower right. There is an adjacent foreground pixel, if present, dividing the foreground pixel into the same target area in the same target area, for example, using a uniform connectivity flag for each target area (eg, as an example in FIG. 15 The numbers 1, 2, 3) are used; if not present, the foreground pixel is the separated foreground pixel and is typically the residual background area in the non-iris region.
图16图示了根据本发明的实施例的图15中的拥有的前景像素 点的数目最多的虹膜区域的示意图,其中位于中间部位的目标区域拥有的前景像素点的数目最多,因此将该目标区域确定虹膜区域。FIG. 16 illustrates the owned foreground pixels of FIG. 15 in accordance with an embodiment of the present invention. A schematic diagram of the iris region having the largest number of dots, wherein the target region located at the intermediate portion has the largest number of foreground pixel points, and thus the target region is determined to be the iris region.
在一个实施例中,还可以消除该虹膜区域中的孔洞、不毗邻的边界以便得到完整的虹膜区域;以及在该完整的虹膜区域中进行虹膜提取。具体来说,根据如下公式消除虹膜区域中的孔洞、不毗邻的边界:In one embodiment, holes in the iris region, non-adjacent boundaries can also be eliminated to obtain a complete iris region; and iris extraction is performed in the intact iris region. Specifically, the holes in the iris area and the non-adjacent boundaries are eliminated according to the following formula:
Figure PCTCN2015070058-appb-000020
Figure PCTCN2015070058-appb-000020
Figure PCTCN2015070058-appb-000021
Figure PCTCN2015070058-appb-000021
Figure PCTCN2015070058-appb-000022
Figure PCTCN2015070058-appb-000022
其中A表示虹膜区域,α表示A中的任意一点,B表示结构元素集合,b表示B中的任意一点,
Figure PCTCN2015070058-appb-000023
表示二维整数网格。
Where A represents the iris region, α represents any point in A, B represents a collection of structural elements, and b represents any point in B.
Figure PCTCN2015070058-appb-000023
Represents a two-dimensional integer grid.
下面针对这些公式如何消除虹膜区域中的孔洞、不毗邻的边界进行详细描述。The following is a detailed description of how these formulas eliminate holes and non-contiguous boundaries in the iris region.
假设A表示虹膜区域,B表示半径为R的结构元素集合(例如半径为3的圆),则上述公式可以形象地理解为将B的中心沿A的边缘点滑动,得到重叠区域;将该重叠区域与B进行并集运算,得到边缘被扩大的、孔洞、不毗邻的边界被消除的虹膜区域(也可被称为膨胀的虹膜区域)。然后,将B的中心沿该膨胀的虹膜区域的边缘点滑动,得到另一重叠区域;将该另一重叠区域与B进行交集取非运算,得到上述完整的虹膜区域(也可被称为腐蚀的虹膜区域),其虹膜区域的尺寸保持不变但虹膜区域内部可能出现的孔洞、不毗邻的边界被消除。Assuming that A represents the iris region and B represents a set of structural elements with a radius R (for example, a circle with a radius of 3), the above formula can be visually understood as sliding the center of B along the edge point of A to obtain an overlapping region; The region performs a union operation with B to obtain an iris region (also referred to as an expanded iris region) whose edges are enlarged, holes, and non-adjacent boundaries are eliminated. Then, the center of B is slid along the edge point of the expanded iris region to obtain another overlapping region; the other overlapping region is intersected with B to obtain the above-mentioned complete iris region (also referred to as corrosion). The iris area), the size of the iris area remains the same, but the holes, non-adjacent boundaries that may appear inside the iris area are eliminated.
图17图示了根据本发明的实施例的消除虹膜区域中的孔洞、不毗邻的边界所得到完整的虹膜区域的示意图,由此可见,已经能够获得没有任何干扰、孔洞的完整的虹膜区域。Figure 17 illustrates a schematic diagram of a complete iris region resulting from the elimination of holes in the iris region, non-adjacent boundaries, in accordance with an embodiment of the present invention, whereby it can be seen that a complete iris region without any interference, holes, has been obtained.
本领域技术人员应当知晓,上述涉及通过图像形态学后处理去除眼镜框残存片段的操作仅仅是示例性而非限制性的,其他任何能够实现去除眼镜框残存片段的实施例均应当纳入本发明的保护范 围。Those skilled in the art will appreciate that the above-described operations involving the removal of the remaining segments of the eyeglass frame by image morphology post-processing are merely exemplary and not limiting, and any other embodiments capable of achieving the removal of the remaining segments of the eyeglass frame should be incorporated into the present invention. Protection Wai.
本领域技术人员应当知晓,上述涉及判别眼镜框位置的操作仅仅是示例性而非限制性的,其他任何能够实现判别眼镜框位置的实施例均应当纳入本发明的保护范围。It should be appreciated by those skilled in the art that the above-described operations relating to determining the position of the eyeglass frame are merely exemplary and not limiting, and any other embodiment capable of realizing the position of the eyeglass frame should be included in the scope of protection of the present invention.
图18图示了根据本发明的实施例的不均匀光照条件下的虹膜提取设备的结构框图,其包括滤除装置1802、计算装置1804和虹膜提取装置1806。其中,滤除装置1802用于从包含虹膜的图像中滤除不均匀光照强度分量,计算装置1804用于在滤除装置1802滤除后的图像中计算色差,虹膜提取装置1806用于根据计算装置1804计算的色差进行虹膜提取。18 illustrates a block diagram of a structure of an iris extraction apparatus under non-uniform illumination conditions, including a filtering device 1802, a computing device 1804, and an iris extraction device 1806, in accordance with an embodiment of the present invention. Wherein the filtering device 1802 is configured to filter out the uneven illumination intensity component from the image containing the iris, the computing device 1804 is configured to calculate the color difference in the filtered image of the filtering device 1802, and the iris extraction device 1806 is configured to use the computing device. The iris was extracted by the color difference calculated in 1804.
在一个实施例中,滤除装置1802包括:转换单元,用于对所述图像的色彩模式进行转换,使得转换后的所述图像中的所述不均匀光照强度分量和色彩分量相分离。In one embodiment, the filtering device 1802 includes a conversion unit for converting a color mode of the image such that the uneven illumination intensity component and the color component in the converted image are separated.
在一个实施例中,转换单元将所述图像从RGB色彩模式转换到Lab色彩模式,以得到L通道、a通道和b通道中的相应的图像分量,其中L通道中的图像分量包括该不均匀光照强度分量,a通道中的图像分量和b通道中的图像分量包括该色彩分量。在一个实施例中,计算装置1804包括:第一计算单元,用于利用a通道中的图像分量和b通道中的图像分量计算色差。In one embodiment, the conversion unit converts the image from the RGB color mode to the Lab color mode to obtain respective image components in the L channel, the a channel, and the b channel, wherein the image component in the L channel includes the unevenness The illumination intensity component, the image component in the a channel and the image component in the b channel include the color component. In one embodiment, computing device 1804 includes a first computing unit for calculating a color difference using image components in the a channel and image components in the b channel.
在一个实施例中,计算单元根据如下公式,利用a通道中的图像分量和b通道中的图像分量计算色差:
Figure PCTCN2015070058-appb-000024
其中ΔEab表示所述色差,Δ(a)表示a通道内两个数值之差,Δ(b)表示b通道内两个数值之差。
In one embodiment, the computing unit calculates the color difference using the image component in the a channel and the image component in the b channel according to the following formula:
Figure PCTCN2015070058-appb-000024
Where ΔE ab represents the chromatic aberration, Δ(a) represents the difference between two values in the a channel, and Δ(b) represents the difference between the two values in the b channel.
在一个实施例中,所述虹膜提取装置1806包括:第一比较单元,用于将色差与预定的色差阈值进行比较,得到第一比较结果;以及虹膜提取单元,用于根据第一比较结果进行虹膜提取。In one embodiment, the iris extraction device 1806 includes: a first comparison unit configured to compare the color difference with a predetermined color difference threshold to obtain a first comparison result; and an iris extraction unit configured to perform the first comparison result according to the first comparison result Iris extraction.
在一个实施例中,滤除装置1802还包括:滤除单元,用于从L通道中的图像分量中滤除不均匀光照强度分量。In one embodiment, the filtering device 1802 further includes a filtering unit for filtering out the uneven light intensity component from the image components in the L channel.
在一个实施例中,滤除单元包括:第一运算模块,用于将L通 道中的图像分量进行对数运算;第二运算模块,用于将对数变换后的图像分量进行快速傅里叶变换;第三运算模块,用于通过高通滤波滤除快速傅里叶变换后的图像分量中的低频部分;第四运算模块,用于将高通滤波后的图像分量进行快速傅里叶变换的逆变换;以及第五运算模块,用于将逆变换后的图像分量进行指数运算。In one embodiment, the filtering unit includes: a first computing module for connecting the L The image component in the track performs a logarithm operation; the second operation module is configured to perform fast Fourier transform on the log transformed image component; and the third operation module is configured to filter the fast Fourier transform by high-pass filtering a low frequency portion of the image component; a fourth operation module for performing an inverse transform of the high-pass filtered image component on the fast Fourier transform; and a fifth operation module for performing an exponential operation on the inverse transformed image component.
在一个实施例中,虹膜提取装置1806还包括:第二计算单元,用于计算滤除不均匀光照强度分量后的L通道中的图像分量对应的灰度值;第二比较单元,用于将该灰度值与预定的灰度阈值进行比较,得到第二比较结果;第三计算单元,用于计算第一比较结果和第二比较结果的第一交集;以及虹膜提取单元用于根据第一交集进行虹膜提取。如上面已经详细说明的,在没有戴眼镜的情况下,虹膜提取装置1806的执行可以到此为止。In an embodiment, the iris extraction device 1806 further includes: a second calculation unit, configured to calculate a gray value corresponding to the image component in the L channel after filtering the uneven illumination intensity component; and a second comparison unit, configured to The gray value is compared with a predetermined gray threshold to obtain a second comparison result; a third calculating unit is configured to calculate a first intersection of the first comparison result and the second comparison result; and an iris extraction unit is configured to use the first Intersection for iris extraction. As has been explained in detail above, the execution of the iris extraction device 1806 can be performed here without wearing glasses.
在一个实施例中,虹膜提取装置1806还包括:第四计算单元,用于计算滤除不均匀光照强度分量后的L通道中的图像分量中任意一点的亮度减去该点周围的预定数目的点的亮度的平均值的差值;第三比较单元,用于将该差值与预定的亮度差值阈值进行比较,得到第三比较结果;第四计算单元,用于计算第一交集和第三比较结果的第二交集;以及虹膜提取单元用于根据第二交集进行虹膜提取。In one embodiment, the iris extraction device 1806 further includes: a fourth calculation unit configured to calculate a brightness of any one of the image components in the L channel after filtering the uneven illumination intensity component minus a predetermined number of the points around the point a difference between the average values of the brightness of the points; a third comparing unit for comparing the difference with a predetermined brightness difference threshold to obtain a third comparison result; and a fourth calculating unit for calculating the first intersection and the first a second intersection of the three comparison results; and an iris extraction unit for performing iris extraction based on the second intersection.
在一个实施例中,第四计算单元根据如下公式计算差值:In one embodiment, the fourth calculation unit calculates the difference according to the following formula:
Figure PCTCN2015070058-appb-000025
Figure PCTCN2015070058-appb-000025
其中,(x,y)表示滤除不均匀光照强度分量后的L通道中的图像分量中的任意一点的坐标,f(x,y)表示坐标(x,y)的点的亮度值,
Figure PCTCN2015070058-appb-000026
表示非负整数集合,n表示非负整数。
Where (x, y) represents the coordinates of any point in the image component in the L channel after filtering out the uneven illumination intensity component, and f(x, y) represents the luminance value of the point of the coordinate (x, y),
Figure PCTCN2015070058-appb-000026
Represents a non-negative integer set and n represents a non-negative integer.
在一个实施例中,虹膜提取装置1806还包括:判断单元,用于分别判断第二交集对应的二值图像中的每一前景像素点是否存在毗邻的前景像素点,如果存在,则将该前景像素点与该毗邻的前景像素点划分在同一目标区域内从而得到多个目标区域;第五计算单元,用于分别计算该多个目标区域所拥有的前景像素点的数目;确定单 元,用于将数目最多的目标区域确定为虹膜区域;以及虹膜提取单元在该虹膜区域中进行虹膜提取。In an embodiment, the iris extraction device 1806 further includes: a determining unit, configured to respectively determine whether each foreground pixel in the binary image corresponding to the second intersection has an adjacent foreground pixel, and if present, the foreground Pixel points are divided into the same target area in the same target area to obtain a plurality of target areas; and a fifth calculating unit is configured to separately calculate the number of foreground pixel points owned by the plurality of target areas; The element is used to determine the target area of the largest number as the iris area; and the iris extraction unit performs iris extraction in the iris area.
在一个实施例中,虹膜提取装置1806还包括:消除单元,用于消除虹膜区域中的孔洞、不毗邻的边界,得到完整的虹膜区域;以及虹膜提取单元在该完整的虹膜区域中进行虹膜提取。In one embodiment, the iris extraction device 1806 further includes: an elimination unit for eliminating holes in the iris region, non-adjacent boundaries, to obtain a complete iris region; and an iris extraction unit for iris extraction in the intact iris region .
在一个实施例中,消除单元根据如下公式消除虹膜区域中的孔洞、不毗邻的边界:In one embodiment, the elimination unit eliminates holes, non-adjacent boundaries in the iris region according to the following formula:
Figure PCTCN2015070058-appb-000027
Figure PCTCN2015070058-appb-000027
Figure PCTCN2015070058-appb-000028
Figure PCTCN2015070058-appb-000028
其中A表示虹膜区域,α表示A中的任意一点,B表示结构元素集合,b表示B中的任意一点,
Figure PCTCN2015070058-appb-000029
表示二维整数网格。
Where A represents the iris region, α represents any point in A, B represents a collection of structural elements, and b represents any point in B.
Figure PCTCN2015070058-appb-000029
Represents a two-dimensional integer grid.
在一个实施例中,消除单元还根据如下公式消除虹膜区域中的孔洞、不毗邻的边界:In one embodiment, the elimination unit also eliminates holes, non-adjacent boundaries in the iris region according to the following formula:
Figure PCTCN2015070058-appb-000030
Figure PCTCN2015070058-appb-000030
其中A表示虹膜区域,α表示A中的任意一点,B表示结构元素集合,
Figure PCTCN2015070058-appb-000031
表示二维整数网格。
Where A represents the iris region, α represents any point in A, and B represents a collection of structural elements.
Figure PCTCN2015070058-appb-000031
Represents a two-dimensional integer grid.
综上所述,根据本发明的上述实施例,通过滤除不同的不均匀光照强度分量,从而可以避免该不均匀光照条件对虹膜提取的影响;通过L通道中的图像分量去除局部均值计算可以得到眼镜架位置,从而消除眼镜架对虹膜提取的影响。In summary, according to the above embodiment of the present invention, by filtering out different uneven illumination intensity components, the influence of the uneven illumination condition on the iris extraction can be avoided; the local mean calculation can be performed by removing the image component in the L channel. The position of the spectacle frame is obtained, thereby eliminating the effect of the spectacle frame on iris extraction.
虽然已经参考若干具体实施例描述了本发明,但是应该理解,本发明并不限于所公开的具体实施例。本发明旨在涵盖所附权利要求的精神和范围内所包括的各种修改和等同布置。所附权利要求的范围符合最宽泛的解释,从而包含所有这样的修改及等同结构和功能。 Although the invention has been described with reference to a particular embodiment thereof, it is understood that the invention is not limited to the specific embodiments disclosed. The invention is intended to cover various modifications and equivalents The scope of the following claims is to be accorded

Claims (30)

  1. 一种不均匀光照条件下的虹膜提取方法,包括:An iris extraction method under uneven illumination conditions, comprising:
    从包含虹膜的图像中滤除不均匀光照强度分量;Filtering out uneven illumination intensity components from images containing irises;
    在滤除后的所述图像中计算色差;以及Calculating the color difference in the filtered image;
    根据所述色差进行虹膜提取。Iris extraction is performed according to the color difference.
  2. 根据权利要求1所述的方法,其中所述从包含虹膜的图像中滤除不均匀光照强度分量包括:The method of claim 1 wherein said filtering out the uneven illumination intensity component from the image comprising the iris comprises:
    对所述图像的色彩模式进行转换,使得转换后的所述图像中的所述不均匀光照强度分量和色彩分量相分离。The color mode of the image is converted such that the uneven illumination intensity component and the color component in the converted image are separated.
  3. 根据权利要求2所述的方法,其中将所述图像从RGB色彩模式转换到Lab色彩模式,得到L通道、a通道和b通道中的相应的图像分量,以使得转换后的所述图像中的所述不均匀光照强度分量和色彩分量相分离,其中所述L通道中的图像分量包括所述不均匀光照强度分量,所述a通道中的图像分量和所述b通道中的图像分量包括所述色彩分量。The method of claim 2, wherein converting the image from an RGB color mode to a Lab color mode results in respective image components in the L channel, the a channel, and the b channel such that the converted image is in the image The uneven illumination intensity component and the color component are separated, wherein the image component in the L channel includes the uneven illumination intensity component, and the image component in the a channel and the image component in the b channel include The color component.
  4. 根据权利要求3所述的方法,其中所述在滤除后的所述图像中计算色差包括:The method of claim 3 wherein said calculating the color difference in said filtered image comprises:
    利用所述a通道中的图像分量和所述b通道中的图像分量计算所述色差。The color difference is calculated using image components in the a channel and image components in the b channel.
  5. 根据权利要求4所述的方法,其中根据如下公式,利用所述a通道中的图像分量和所述b通道中的图像分量计算所述色差:The method of claim 4, wherein the color difference is calculated using image components in the a channel and image components in the b channel according to the following formula:
    Figure PCTCN2015070058-appb-100001
    Figure PCTCN2015070058-appb-100001
    其中ΔEab表示所述色差,Δ(a)表示所述a通道内两个数值之差,Δ(b)表示所述b通道内两个数值之差。Where ΔE ab represents the chromatic aberration, Δ(a) represents the difference between the two values in the a channel, and Δ(b) represents the difference between the two values in the b channel.
  6. 根据权利要求4或5所述的方法,其中所述根据所述色差进行虹膜提取包括:The method according to claim 4 or 5, wherein said performing iris extraction based on said color difference comprises:
    将所述色差与预定的色差阈值进行比较,得到第一比较结果;以及 Comparing the color difference with a predetermined color difference threshold to obtain a first comparison result;
    根据所述第一比较结果进行虹膜提取。Iris extraction is performed based on the first comparison result.
  7. 根据权利要求3所述的方法,其中所述从包含虹膜的图像中滤除不均匀光照强度分量还包括:The method of claim 3 wherein said filtering out the uneven illumination intensity component from the image comprising the iris further comprises:
    从所述L通道中的图像分量中滤除所述不均匀光照强度分量。The uneven illumination intensity component is filtered from image components in the L channel.
  8. 根据权利要求7所述的方法,其中从所述L通道中的图像分量中滤除所述不均匀光照强度分量包括:The method of claim 7 wherein filtering the non-uniform illumination intensity component from image components in the L channel comprises:
    将所述L通道中的图像分量进行对数运算;Performing a logarithmic operation on image components in the L channel;
    将对数变换后的所述图像分量进行快速傅里叶变换;Performing a fast Fourier transform on the log transformed image component;
    通过高通滤波滤除快速傅里叶变换后的所述图像分量中的低频部分;Filtering the low frequency portion of the image component after the fast Fourier transform by high pass filtering;
    将高通滤波后的所述图像分量进行快速傅里叶变换的逆变换;以及Performing an inverse transform of the high-pass filtered image component on the fast Fourier transform;
    将逆变换后的所述图像分量进行指数运算。The inversely transformed image component is subjected to an exponential operation.
  9. 根据权利要求7所述的方法,其中所述根据所述色差进行虹膜提取还包括:The method of claim 7, wherein said performing iris extraction based on said color difference further comprises:
    将所述色差与预定的色差阈值进行比较,得到第一比较结果;Comparing the color difference with a predetermined color difference threshold to obtain a first comparison result;
    计算滤除所述不均匀光照强度分量后的所述L通道中的图像分量对应的灰度值;Calculating a gray value corresponding to an image component in the L channel after filtering the uneven light intensity component;
    将所述灰度值与预定的灰度阈值进行比较,得到第二比较结果;Comparing the gray value with a predetermined gray threshold to obtain a second comparison result;
    计算所述第一比较结果和所述第二比较结果的第一交集;以及Calculating a first intersection of the first comparison result and the second comparison result;
    根据所述第一交集进行虹膜提取。Iris extraction is performed according to the first intersection.
  10. 根据权利要求9所述的方法,其中所述根据所述色差进行虹膜提取还包括:The method of claim 9, wherein said performing iris extraction based on said color difference further comprises:
    计算滤除所述不均匀光照强度分量后的所述L通道中的图像分量中任意一点的亮度减去所述点周围的预定数目的点的亮度的平均值的差值;Calculating a difference between an average value of luminances of the image components in the L channel after filtering the uneven illumination intensity component minus an average of luminances of a predetermined number of points around the point;
    将所述差值与预定的亮度差值阈值进行比较,得到第三比较结果;Comparing the difference with a predetermined brightness difference threshold to obtain a third comparison result;
    计算所述第一交集和所述第三比较结果的第二交集;以及 Calculating a second intersection of the first intersection and the third comparison result;
    根据所述第二交集进行虹膜提取。Iris extraction is performed according to the second intersection.
  11. 根据权利要求10所述的方法,其中根据如下公式计算所述差值:The method of claim 10 wherein said difference is calculated according to the following formula:
    Figure PCTCN2015070058-appb-100002
    Figure PCTCN2015070058-appb-100002
    其中,(x,y)表示滤除所述不均匀光照强度分量后的所述L通道中的图像分量中的任意一点的坐标,f(x,y)表示坐标(x,y)的点的亮度值,
    Figure PCTCN2015070058-appb-100003
    表示非负整数集合,n表示非负整数。
    Where (x, y) represents the coordinates of any point in the image component in the L channel after filtering out the uneven illumination intensity component, and f(x, y) represents the point of the coordinate (x, y) Brightness value,
    Figure PCTCN2015070058-appb-100003
    Represents a non-negative integer set and n represents a non-negative integer.
  12. 根据权利要求10或11所述的方法,根据所述第二交集进行虹膜提取包括:The method according to claim 10 or 11, wherein extracting the iris according to the second intersection comprises:
    分别判断所述第二交集对应的二值图像中的每一前景像素点是否存在毗邻的前景像素点,如果存在,则将所述前景像素点与所述毗邻的前景像素点划分在同一目标区域内从而得到多个目标区域;Determining, respectively, whether each foreground pixel in the binary image corresponding to the second intersection has an adjacent foreground pixel, and if present, dividing the foreground pixel and the adjacent foreground pixel into the same target area Thereby obtaining a plurality of target areas;
    分别计算所述多个目标区域所拥有的前景像素点的数目;Calculating, respectively, the number of foreground pixel points owned by the plurality of target regions;
    将所述数目最多的目标区域确定为虹膜区域;以及Determining the target area of the largest number as the iris area;
    在所述虹膜区域中进行虹膜提取。Iris extraction is performed in the iris region.
  13. 根据权利要求12所述的方法,在所述虹膜区域中进行虹膜提取包括:The method of claim 12, wherein performing iris extraction in the iris region comprises:
    消除所述虹膜区域中的孔洞、不毗邻的边界,得到完整的虹膜区域;以及Eliminating holes and non-adjacent boundaries in the iris region to obtain a complete iris region;
    在所述完整的虹膜区域中进行虹膜提取。Iris extraction is performed in the intact iris region.
  14. 根据权利要求13所述的方法,其中根据如下公式消除所述虹膜区域中的孔洞、不毗邻的边界:The method of claim 13 wherein the holes in the iris region, non-adjacent boundaries are eliminated according to the following formula:
    Figure PCTCN2015070058-appb-100004
    Figure PCTCN2015070058-appb-100004
    Figure PCTCN2015070058-appb-100005
    Figure PCTCN2015070058-appb-100005
    其中A表示所述虹膜区域,α表示A中的任意一点,B表示结构元素集合,b表示B中的任意一点,
    Figure PCTCN2015070058-appb-100006
    表示二维整数网格。
    Where A represents the iris region, α represents any point in A, B represents a collection of structural elements, and b represents any point in B.
    Figure PCTCN2015070058-appb-100006
    Represents a two-dimensional integer grid.
  15. 根据权利要求14所述的方法,其中还根据如下公式消除虹膜 区域中的孔洞、不毗邻的边界:The method of claim 14 wherein the iris is further eliminated according to the following formula Holes in the area, non-adjacent boundaries:
    Figure PCTCN2015070058-appb-100007
    Figure PCTCN2015070058-appb-100007
    其中A表示所述虹膜区域,α表示A中的任意一点,B表示结构元素集合,
    Figure PCTCN2015070058-appb-100008
    表示二维整数网格。
    Where A represents the iris region, α represents any point in A, and B represents a collection of structural elements.
    Figure PCTCN2015070058-appb-100008
    Represents a two-dimensional integer grid.
  16. 一种不均匀光照条件下的虹膜提取设备,包括:An iris extraction device under uneven illumination conditions, comprising:
    滤除装置,用于从包含虹膜的图像中滤除不均匀光照强度分量;a filtering device for filtering out uneven light intensity components from the image containing the iris;
    计算装置,用于在滤除后的所述图像中计算色差;以及a computing device for calculating a color difference in the filtered image;
    虹膜提取装置,用于根据所述色差进行虹膜提取。An iris extraction device for performing iris extraction according to the chromatic aberration.
  17. 根据权利要求16所述的设备,其中所述滤除装置包括:The apparatus of claim 16 wherein said filtering means comprises:
    转换单元,用于对所述图像的色彩模式进行转换,使得转换后的所述图像中的所述不均匀光照强度分量和色彩分量相分离。And a converting unit, configured to convert a color mode of the image, so that the uneven light intensity component and the color component in the converted image are separated.
  18. 根据权利要求17所述的设备,其中所述转换单元将所述图像从RGB色彩模式转换到Lab色彩模式,得到L通道、a通道和b通道中的相应的图像分量,以使得转换后的所述图像中的所述不均匀光照强度分量和色彩分量相分离,其中所述L通道中的图像分量包括所述不均匀光照强度分量,所述a通道中的图像分量和所述b通道中的图像分量包括所述色彩分量。The apparatus according to claim 17, wherein said converting unit converts said image from an RGB color mode to a Lab color mode, and obtains corresponding image components in the L channel, the a channel, and the b channel, so that the converted image The uneven illumination intensity component and the color component in the image are separated, wherein the image component in the L channel includes the uneven illumination intensity component, an image component in the a channel, and the b channel The image component includes the color component.
  19. 根据权利要求18所述的设备,其中所述计算装置包括:The device of claim 18 wherein said computing device comprises:
    第一计算单元,用于利用所述a通道中的图像分量和所述b通道中的图像分量计算所述色差。a first calculating unit, configured to calculate the color difference by using an image component in the a channel and an image component in the b channel.
  20. 根据权利要求19所述的设备,其中所述第一计算单元根据如下公式,利用所述a通道中的图像分量和所述b通道中的图像分量计算所述色差:The apparatus according to claim 19, wherein said first calculating unit calculates said color difference using an image component in said a channel and an image component in said b channel according to the following formula:
    Figure PCTCN2015070058-appb-100009
    Figure PCTCN2015070058-appb-100009
    其中ΔEab表示所述色差,Δ(a)表示所述a通道内两个数值之差,Δ(b)表示所述b通道内两个数值之差。Where ΔE ab represents the chromatic aberration, Δ(a) represents the difference between the two values in the a channel, and Δ(b) represents the difference between the two values in the b channel.
  21. 根据权利要求19或20所述的设备,其中所述虹膜提取装置包括:The apparatus according to claim 19 or 20, wherein said iris extraction means comprises:
    第一比较单元,用于将所述色差与预定的色差阈值进行比较,得 到第一比较结果;以及a first comparing unit, configured to compare the color difference with a predetermined color difference threshold To the first comparison result;
    虹膜提取单元,用于根据所述第一比较结果进行虹膜提取。An iris extraction unit configured to perform iris extraction according to the first comparison result.
  22. 根据权利要求18所述的设备,其中所述滤除装置还包括:The apparatus according to claim 18, wherein said filtering means further comprises:
    滤除单元,用于从所述L通道中的图像分量中滤除所述不均匀光照强度分量。And filtering a unit for filtering the uneven light intensity component from image components in the L channel.
  23. 根据权利要求22所述的设备,其中所述滤除单元包括:The apparatus of claim 22 wherein said filtering unit comprises:
    第一运算模块,用于将所述L通道中的图像分量进行对数运算;a first operation module, configured to perform logarithm operation on image components in the L channel;
    第二运算模块,用于将对数变换后的所述图像分量进行快速傅里叶变换;a second operation module, configured to perform fast Fourier transform on the log transformed image component;
    第三运算模块,用于通过高通滤波滤除快速傅里叶变换后的所述图像分量中的低频部分;a third operation module, configured to filter, by high-pass filtering, a low-frequency portion of the image component after fast Fourier transform;
    第四运算模块,用于将高通滤波后的所述图像分量进行快速傅里叶变换的逆变换;以及a fourth operation module, configured to perform inverse transform of the high-pass filtered image component by fast Fourier transform;
    第五运算模块,用于将逆变换后的所述图像分量进行指数运算。And a fifth operation module, configured to perform an exponential operation on the inverse transformed image component.
  24. 根据权利要求22所述的设备,其中所述虹膜提取装置包括:The apparatus according to claim 22, wherein said iris extraction means comprises:
    第一比较单元,用于将所述色差与预定的色差阈值进行比较,得到第一比较结果;a first comparing unit, configured to compare the color difference with a predetermined color difference threshold to obtain a first comparison result;
    第二计算单元,用于计算滤除所述不均匀光照强度分量后的所述L通道中的图像分量对应的灰度值;a second calculating unit, configured to calculate a gray value corresponding to the image component in the L channel after filtering the uneven light intensity component;
    第二比较单元,用于将所述灰度值与预定的灰度阈值进行比较,得到第二比较结果;a second comparing unit, configured to compare the gray value with a predetermined gray threshold to obtain a second comparison result;
    第三计算单元,用于计算所述第一比较结果和所述第二比较结果的第一交集;以及a third calculating unit, configured to calculate a first intersection of the first comparison result and the second comparison result;
    所述虹膜提取单元用于根据所述第一交集进行虹膜提取。The iris extraction unit is configured to perform iris extraction according to the first intersection.
  25. 根据权利要求24所述的设备,其中所述虹膜提取装置还包括:The apparatus according to claim 24, wherein said iris extracting means further comprises:
    第四计算单元,用于计算滤除所述不均匀光照强度分量后的所述L通道中的图像分量中任意一点的亮度减去所述点周围的预定数目的点的亮度的平均值的差值; a fourth calculating unit, configured to calculate a difference between an average value of brightness of the image component in the L channel after filtering the uneven light intensity component minus a brightness of a predetermined number of points around the point Value
    第三比较单元,用于将所述差值与预定的亮度差值阈值进行比较,得到第三比较结果;a third comparing unit, configured to compare the difference with a predetermined brightness difference threshold to obtain a third comparison result;
    第四计算单元,用于计算所述第一交集和所述第三比较结果的第二交集;以及a fourth calculating unit, configured to calculate a second intersection of the first intersection and the third comparison result;
    所述虹膜提取单元用于根据所述第二交集进行虹膜提取。The iris extraction unit is configured to perform iris extraction according to the second intersection.
  26. 根据权利要求25所述的设备,其中所述第四计算单元根据如下公式计算所述差值:The apparatus according to claim 25, wherein said fourth calculating unit calculates said difference value according to the following formula:
    Figure PCTCN2015070058-appb-100010
    Figure PCTCN2015070058-appb-100010
    其中,(x,y)表示滤除所述不均匀光照强度分量后的所述L通道中的图像分量中的任意一点的坐标,f(x,y)表示坐标(x,y)的点的亮度值,
    Figure PCTCN2015070058-appb-100011
    表示非负整数集合,n表示非负整数。
    Where (x, y) represents the coordinates of any point in the image component in the L channel after filtering out the uneven illumination intensity component, and f(x, y) represents the point of the coordinate (x, y) Brightness value,
    Figure PCTCN2015070058-appb-100011
    Represents a non-negative integer set and n represents a non-negative integer.
  27. 根据权利要求25或26所述的设备,所述虹膜提取装置还包括:The apparatus according to claim 25 or 26, wherein the iris extraction device further comprises:
    判断单元,用于分别判断所述第二交集对应的二值图像中的每一前景像素点是否存在毗邻的前景像素点,如果存在,则将所述前景像素点与所述毗邻的前景像素点划分在同一目标区域内从而得到多个目标区域;a determining unit, configured to determine, respectively, whether each foreground pixel in the binary image corresponding to the second intersection has an adjacent foreground pixel, and if present, the foreground pixel and the adjacent foreground pixel Divided into the same target area to obtain multiple target areas;
    第五计算单元,用于分别计算所述多个目标区域所拥有的前景像素点的数目;a fifth calculating unit, configured to separately calculate a number of foreground pixel points owned by the plurality of target areas;
    确定单元,用于将所述数目最多的目标区域确定为虹膜区域;以及a determining unit, configured to determine the target area of the largest number as an iris area;
    所述虹膜提取单元在所述虹膜区域中进行虹膜提取。The iris extraction unit performs iris extraction in the iris region.
  28. 根据权利要求27所述的设备,所述虹膜提取装置还包括:The apparatus according to claim 27, wherein the iris extraction device further comprises:
    消除单元,用于消除所述虹膜区域中的孔洞、不毗邻的边界,得到完整的虹膜区域;以及a eliminating unit for eliminating holes in the iris region, non-adjacent boundaries, and obtaining a complete iris region;
    所述虹膜提取单元在所述完整的虹膜区域中进行虹膜提取。The iris extraction unit performs iris extraction in the intact iris region.
  29. 根据权利要求28所述的设备,其中所述消除单元根据如下公式消除所述虹膜区域中的孔洞、不毗邻的边界: The apparatus according to claim 28, wherein said eliminating unit eliminates holes, non-adjacent boundaries in said iris region according to the following formula:
    Figure PCTCN2015070058-appb-100012
    Figure PCTCN2015070058-appb-100012
    Figure PCTCN2015070058-appb-100013
    Figure PCTCN2015070058-appb-100013
    其中A表示所述虹膜区域,α表示A中的任意一点,B表示结构元素集合,b表示B中的任意一点,
    Figure PCTCN2015070058-appb-100014
    表示二维整数网格。
    Where A represents the iris region, α represents any point in A, B represents a collection of structural elements, and b represents any point in B.
    Figure PCTCN2015070058-appb-100014
    Represents a two-dimensional integer grid.
  30. 根据权利要求29所述的设备,其中所述消除单元还根据如下公式消除虹膜区域中的孔洞、不毗邻的边界:The apparatus according to claim 29, wherein said eliminating unit further eliminates holes, non-adjacent boundaries in the iris region according to the following formula:
    Figure PCTCN2015070058-appb-100015
    Figure PCTCN2015070058-appb-100015
    其中A表示所述虹膜区域,α表示A中的任意一点,B表示结构元素集合,
    Figure PCTCN2015070058-appb-100016
    表示二维整数网格。
    Where A represents the iris region, α represents any point in A, and B represents a collection of structural elements.
    Figure PCTCN2015070058-appb-100016
    Represents a two-dimensional integer grid.
PCT/CN2015/070058 2014-01-09 2015-01-04 Method and device for extracting iris image image under condition of non-uniform illumination WO2015103953A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201410016888.8A CN104778729A (en) 2014-01-09 2014-01-09 Iris extraction method and equipment under uneven illumination condition
CN201410016888.8 2014-01-09

Publications (1)

Publication Number Publication Date
WO2015103953A1 true WO2015103953A1 (en) 2015-07-16

Family

ID=53523539

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2015/070058 WO2015103953A1 (en) 2014-01-09 2015-01-04 Method and device for extracting iris image image under condition of non-uniform illumination

Country Status (2)

Country Link
CN (1) CN104778729A (en)
WO (1) WO2015103953A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117094913A (en) * 2023-10-17 2023-11-21 深圳市城市交通规划设计研究中心股份有限公司 Asphalt pavement detection image quality improvement method under non-uniform illumination

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107463913B (en) * 2017-08-11 2020-07-14 苏炯龙 L ab-based image collagen fiber recognition and calculation system and method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1760887A (en) * 2004-10-11 2006-04-19 中国科学院自动化研究所 The robust features of iris image extracts and recognition methods
CN102542243A (en) * 2010-12-17 2012-07-04 北京理工大学 LBP (Local Binary Pattern) image and block encoding-based iris feature extracting method
CN102819850A (en) * 2012-08-16 2012-12-12 李勃 Method for detecting edge of color image on basis of local self-adaption color difference threshold

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1209732C (en) * 2002-10-09 2005-07-06 中国科学院自动化研究所 Iris identifying method based on ripple analysis and zero passage description
JP4537143B2 (en) * 2004-07-30 2010-09-01 キヤノン株式会社 Image processing apparatus and method, imaging apparatus, and program
JP5089405B2 (en) * 2008-01-17 2012-12-05 キヤノン株式会社 Image processing apparatus, image processing method, and imaging apparatus
CN101930543B (en) * 2010-08-27 2012-06-27 南京大学 Method for adjusting eye image in self-photographed video

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1760887A (en) * 2004-10-11 2006-04-19 中国科学院自动化研究所 The robust features of iris image extracts and recognition methods
CN102542243A (en) * 2010-12-17 2012-07-04 北京理工大学 LBP (Local Binary Pattern) image and block encoding-based iris feature extracting method
CN102819850A (en) * 2012-08-16 2012-12-12 李勃 Method for detecting edge of color image on basis of local self-adaption color difference threshold

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117094913A (en) * 2023-10-17 2023-11-21 深圳市城市交通规划设计研究中心股份有限公司 Asphalt pavement detection image quality improvement method under non-uniform illumination
CN117094913B (en) * 2023-10-17 2024-03-01 深圳市城市交通规划设计研究中心股份有限公司 Asphalt pavement detection image quality improvement method under non-uniform illumination

Also Published As

Publication number Publication date
CN104778729A (en) 2015-07-15

Similar Documents

Publication Publication Date Title
Vidya et al. Skin cancer detection using machine learning techniques
Guan et al. Accurate segmentation of partially overlapping cervical cells based on dynamic sparse contour searching and GVF snake model
Eadgahi et al. Localization of hard exudates in retinal fundus image by mathematical morphology operations
Fiorese et al. VirtualShave: automated hair removal from digital dermatoscopic images
Kamboj A color-based approach for melanoma skin cancer detection
Xiao et al. Retinal hemorrhage detection by rule-based and machine learning approach
CN109785321A (en) Meibomian gland method for extracting region based on deep learning and Gabor filter
Diaz et al. Glaucoma diagnosis by means of optic cup feature analysis in color fundus images
Pal et al. Mathematical morphology aided optic disk segmentation from retinal images
Madhankumar et al. Characterization of skin lesions
Brancati et al. Automatic segmentation of pigment deposits in retinal fundus images of Retinitis Pigmentosa
WO2015103953A1 (en) Method and device for extracting iris image image under condition of non-uniform illumination
Hashim et al. Optic disc boundary detection from digital fundus images
Zhou et al. A novel approach for red lesions detection using superpixel multi-feature classification in color fundus images
Roy et al. Detection of retinal microaneurysms using fractal analysis and feature extraction technique
Victor et al. Detection of skin cancer cells-A review
Thakur et al. Pre-processing of retinal images for removal of outliers
Reza et al. Automatic detection of optic disc in fundus images by curve operator
Elbalaoui et al. Segmentation of optic disc from fundus images
Akhade et al. Automatic optic disc detection in digital fundus images using image processing techniques
Helwan ITDS: Iris tumor detection system using image processing techniques
Tamilarasi et al. Template matching algorithm for exudates detection from retinal fundus images
Kamil et al. A novel approach for optic disc detection in RGB retinal fundus images
Patil et al. Automated macula proximity diagnosis for early finding of diabetic macular edema
Tulasigeri et al. An advanced thresholding algorithm for diagnosis of glaucoma in fundus images

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 15735094

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 15735094

Country of ref document: EP

Kind code of ref document: A1