WO2010103850A1 - 画像識別子抽出装置 - Google Patents

画像識別子抽出装置 Download PDF

Info

Publication number
WO2010103850A1
WO2010103850A1 PCT/JP2010/001770 JP2010001770W WO2010103850A1 WO 2010103850 A1 WO2010103850 A1 WO 2010103850A1 JP 2010001770 W JP2010001770 W JP 2010001770W WO 2010103850 A1 WO2010103850 A1 WO 2010103850A1
Authority
WO
WIPO (PCT)
Prior art keywords
image identifier
region
extraction
dimension
partial
Prior art date
Application number
PCT/JP2010/001770
Other languages
English (en)
French (fr)
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 日本電気株式会社
Priority to JP2011503733A priority Critical patent/JP4935960B2/ja
Priority to EP10750606.5A priority patent/EP2407932B1/en
Priority to US13/256,088 priority patent/US8744193B2/en
Priority to KR1020117021100A priority patent/KR101335877B1/ko
Priority to CN201080011278.9A priority patent/CN102349093B/zh
Priority to KR1020137022201A priority patent/KR101411127B1/ko
Publication of WO2010103850A1 publication Critical patent/WO2010103850A1/ja

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/28Quantising the image, e.g. histogram thresholding for discrimination between background and foreground patterns
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/42Global feature extraction by analysis of the whole pattern, e.g. using frequency domain transformations or autocorrelation
    • G06V10/421Global feature extraction by analysis of the whole pattern, e.g. using frequency domain transformations or autocorrelation by analysing segments intersecting the pattern
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/50Extraction of image or video features by performing operations within image blocks; by using histograms, e.g. histogram of oriented gradients [HoG]; by summing image-intensity values; Projection analysis

Definitions

  • the present invention relates to an apparatus for extracting an image identifier, which is a feature amount for identifying an image (determining identity).
  • the image identifier is an image feature amount for identifying an image (determining identity). An image identifier extracted from one image is compared with an image identifier extracted from another image, and based on the comparison result, an identity scale (generally, similarity or distance indicating the degree to which the two images are identical) Can be calculated. Further, it is possible to determine whether or not two images are the same by comparing the calculated identity scale with a certain threshold value.
  • the two images are the same is not limited to the case where the two images are the same at the level of the image signal (the pixel values of the pixels constituting the image), but the conversion of the compression format (format) of the image, Image size / aspect ratio conversion, image tone adjustment, various image filtering (sharpening, smoothing, etc.), local processing of images (telop overlay, clipping, etc.), image recapturing
  • image recapturing The case where one image is a duplicated image of the other image by various modification processes such as. If an image identifier is used, for example, a copy of a moving image that is an image or a collection of images can be detected. Therefore, the image identifier can be applied to an illegal copy detection system for images or moving images.
  • FIG. 18 is a diagram illustrating an image identifier extraction method described in Patent Document 1.
  • This image identifier is a feature vector of a plurality of dimensions (16 dimensions in FIG. 18).
  • An average luminance value is calculated from 32 rectangular areas 244 (of which 16 rectangular areas are drawn in FIG. 18) at predetermined positions in the image 240, and between the paired rectangular areas (see FIG. 18, a pair of rectangular regions are connected by a dotted line 248), and a difference in average luminance value is calculated to obtain a 16-dimensional difference vector 250.
  • a composite vector is generated by vector conversion for the difference vector 250, and a 16-dimensional quantization index vector obtained by quantizing each dimension of the composite vector is used as an image identifier.
  • An image identifier composed of feature vectors of a plurality of dimensions has a higher level of identification capability that is the degree to which different images can be identified because the smaller the correlation between dimensions, the greater the amount of information that the feature vector has (smaller redundancy). Become.
  • the correlation between dimensions is the degree of similarity of occurrence of dimension features, and mathematically, for example, the phase between probability variables when occurrence of feature quantities of each dimension is a random variable. It is a value that can be calculated as the number of relations or mutual information. For this reason, it is desirable that the image identifier composed of feature vectors of a plurality of dimensions is designed so that the correlation between dimensions becomes small.
  • the image signal (pixel value of pixels constituting the image) has a correlation between the local regions of the image.
  • a specific image pattern repeatedly appears (especially when it appears in a regular cycle) (for example, an image of a building window arranged in a lattice pattern, see FIG. 19A) or In the case of an image composed of a specific texture (see FIG. 19B), the correlation between local regions of the image increases.
  • an image identifier composed of a feature vector composed of feature amounts extracted from a plurality of local regions of an image is an image with a large correlation between local regions of the image. Since the shape of the local region from which the feature value is extracted in the dimension is the same (in the example of Patent Document 1, a rectangular region having the same shape), the correlation between the dimensions of the extracted feature value becomes large. Therefore, there is a first problem that the identification ability of the image identifier (feature vector) is lowered.
  • the same shape means that the size and the angle (tilt or posture) of the region are the same.
  • Patent Document 1 describes an image in which a specific image pattern repeatedly appears (see FIG. 19A) or an image composed of a specific texture (see FIG. 19B). Image identifiers such as those described have poor identification capabilities.
  • An object of the present invention is that the image identifier extracted from an image having a large correlation between local regions of the image or an image in which signals are concentrated at a specific frequency is reduced in the ability to identify different images.
  • An object of the present invention is to provide an image identifier extraction device that solves the problem.
  • An image identifier extraction device includes a combination of shapes of two partial areas that form a pair and a relative positional relationship between the two partial areas that form a pair.
  • Extracting means for extracting a region feature amount from each partial region of the image according to a plurality of partial region pairs in the image including one or more partial region pairs different from the region pair, and region features for each of the extracted partial regions
  • generating means for generating an image identifier used for identifying the image by using a method different from that of the other one or more partial region pairs for at least one partial region pair based on the quantity.
  • the present invention is configured as described above, it is possible to increase the identification capability that is the degree to which different images can be identified. In particular, this effect is remarkable for an image having a large correlation between local regions of the image.
  • FIG. 1 It is a figure which shows the comparison / quantization method information classified by dimension used in the 5th Embodiment and 6th Embodiment of this invention. It is a figure which shows the comparison / quantization method information classified by dimension used in the 5th Embodiment and 6th Embodiment of this invention. It is a figure which shows the comparison / quantization method information classified by dimension used in the 5th Embodiment and 6th Embodiment of this invention. It is a figure which shows the comparison / quantization method information classified by dimension used in the 5th Embodiment and 6th Embodiment of this invention. It is a figure which shows the extraction method of the image identifier described in patent document 1. FIG.
  • FIG. 1 It is a figure which shows the example of the image from which the correlation between local areas becomes large. It is a block diagram of the collation means which collates a quantization index vector. It is a flowchart which shows the process example of the collation means which collates a quantization index vector. It is a flowchart which shows another example of a process of the collation means which collates a quantization index vector. It is a flowchart which shows another example of a process of the collation means which collates a quantization index vector. It is a figure which shows an example of the index provided with respect to 1024 blocks formed by dividing
  • the image identifier extraction apparatus applies a feature vector (more specifically, a quantization index vector) composed of a plurality of dimensions to an input image.
  • a feature vector (more specifically, a quantization index vector) composed of a plurality of dimensions to an input image.
  • This is a system for outputting as an image identifier, and comprises a dimension determining means 1, an extraction area acquiring means 2, an area feature quantity calculating means 3, and a comparing means 4.
  • the dimension determining unit 1 determines the dimension of the feature vector to be extracted next and supplies it to the extraction region acquiring unit 2.
  • the dimension determining unit 1 sequentially supplies the dimension of the feature vector to be extracted, and the constituent elements after the extraction region acquiring unit 2 extract the feature amount corresponding to the supplied dimension.
  • the dimension determining unit 1 may sequentially supply the first dimension to the Nth dimension to the extraction region acquiring unit 2. If all dimensions of the feature vector are finally supplied, the order of the supplied dimensions may be arbitrary. Multiple dimensions may be supplied in parallel.
  • the extraction area acquisition means 2 is supplied with the dimension-specific extraction area information as an input separately from the dimension from the dimension determination means 1.
  • the dimension-specific extraction area information is information that indicates a pair of a first extraction area and a second extraction area that are associated with each dimension of the feature vector and that extract the feature quantity of that dimension. is there.
  • the first and second extraction regions have the following characteristics as essential conditions.
  • An essential condition for the first and second extraction regions is that the relative positions of the extraction region pairs are different among the dimensions, and that the combinations of the shapes of the extraction region pairs are different between the dimensions.
  • FIG. 2 shows an example of a pair of extraction areas for each dimension indicated by the dimension-specific extraction information that satisfies the above essential conditions.
  • the combination of the shapes of the pairs of extraction regions between dimensions is different.
  • the different shapes are congruent shapes with different angles (for example, the first extraction region of the first dimension in FIG. 2 and the first extraction region of the seventh dimension), or similar shapes with different sizes ( For example, the first dimension second extraction area and the ninth dimension second extraction area in FIG. 2 are also included.
  • the minimum condition is that at least one dimensional pair with a different combination of shapes of the extracted region pairs exists in all dimensions of the feature vector.
  • the more the dimensions of the extraction region pairs (combinations) are different from each other the more desirable. This is because the more the dimensions of the pair of extraction regions (combinations) are different from each other, the smaller the correlation between more dimensions of the feature vector and the higher the discrimination ability.
  • the shape (combination) of a pair of extraction regions may be different between all dimensions of the feature vector.
  • the first extraction region and the second extraction region in a certain dimension do not need to have the same shape as in the ninth dimension in FIG. 2, but are different in shape from other dimensions in FIG. May be. If the shapes of the first extraction region and the second extraction region in each dimension are different, the correlation between the feature amounts extracted from the first extraction region and the second extraction region is small, and the discrimination ability is high. This is desirable. In addition, since the possibility that the first extraction region and the second extraction region become a frequency blind spot with respect to the same frequency at the same time is reduced, the discrimination capability is increased.
  • each extraction area is arbitrary.
  • an arbitrary complicated shape such as the second extraction region of the sixth dimension in FIG. 2 may be used.
  • it may be a line segment or a curve, for example, as in the seventh dimension or the tenth dimension in FIG.
  • the extraction area is composed of a plurality of small areas that are not continuous, such as the first extraction area in the eighth dimension, the first and second extraction areas in the eleventh dimension, and the first extraction area in the twelfth dimension. It may be configured.
  • an extraction region having an arbitrarily complicated shape the correlation between dimensions of feature amounts extracted therefrom can be reduced, and the discrimination ability can be increased.
  • the first extraction region and the second extraction region may partially overlap. Further, either one of the extraction region pairs may be included in the other. In this way, by allowing duplication of pairs of extraction regions, more patterns of extraction region pairs (relative positions and distances) can be obtained, and therefore, the number of patterns that can reduce the correlation between dimensions can be increased. This increases the possibility of higher identification ability.
  • the extraction regions may partially overlap between the dimensions as in each dimension shown in FIG. 2.
  • the patterns of extraction area pairs that can be taken are limited.
  • FIG. 2 by allowing duplication in the extraction areas between dimensions, more patterns of extraction area pairs can be obtained, so that the number of patterns that can reduce the correlation between dimensions can be increased. , The possibility of higher discrimination ability increases. However, if there are too many extracted regions between dimensions, the correlation between dimensions increases, and the discrimination ability decreases, which is not desirable.
  • the combinations of the shapes of the extraction regions between the dimensions (dimension pairs) are different from each other.
  • the combinations of the shapes of the extraction regions between the dimensions (dimension pairs) may be different from each other.
  • the combinations of the shapes of the extraction regions between the dimensions (dimension pairs) may be different from each other.
  • the combinations of the shape of the extraction area pairs between the dimensions (dimension pairs) between dimensions (dimension pairs) of a certain ratio or more. May be different from each other.
  • at least one dimension (dimension pair) the combination of the shapes of the extraction regions between the dimensions (dimension pairs) is mutually determined. It may be different.
  • the dimension pair in which the extraction regions exist in the vicinity of each other has a distance between the two-dimensional extraction regions (the first extraction region and the second extraction region) calculated by a certain predetermined method.
  • the distance between the two-dimensional extraction regions (the first extraction region and the second extraction region) is, for example, the sum of all the pixels included in each dimension extraction region (the first extraction region and the second extraction region). You may calculate as an average value of distance between. Further, for example, it may be calculated as the minimum value of the distance between all the pixels included in each dimension extraction region (first extraction region and second extraction region).
  • centroid coordinates of the extraction areas of each dimension are calculated (the centroid coordinates of the first extraction area and the centroid coordinates of the second extraction area), and the distances between them are calculated (a total of four distances are calculated).
  • the distance between the two-dimensional extraction regions may be defined as an average value or a minimum value thereof. However, the distance between the two-dimensional extraction regions is not limited to these calculation methods.
  • the following effects can be obtained by making the combinations of the shapes of the pairs of extraction regions between the dimensions in which the extraction regions are close to each other (dimension pairs) different from each other.
  • the correlation between the extracted feature amounts tends to increase between the dimensions in which the extraction regions are close to each other. If the combination of the shapes of the pairs of extraction regions between these dimensions is the same, the correlation between the feature amounts is further increased, the redundancy of the entire feature vector is increased, and the discrimination ability is reduced. Therefore, by canceling the increase in the correlation between the extracted feature amounts due to the presence of the extraction regions in the vicinity, the combinations of the shapes of the pairs of extraction regions between those dimensions are made different, It is possible to reduce the correlation between feature amounts, reduce the redundancy of the entire feature vector, and prevent the recognition ability from being lowered.
  • the distance between the extraction regions between the dimensions (dimension pairs) between dimensions (dimension pairs) having the same combination of the shapes of the pairs of extraction regions is larger.
  • the distance between the dimensions (dimension pairs) in which the combinations of the shapes of the extraction region pairs are the same may be larger than a predetermined threshold.
  • the distance between the dimensions (dimension pairs) is set to be larger than a predetermined threshold value. May be.
  • the distance between the dimensions (dimension pairs) between dimensions (dimension pairs) of a certain ratio or more is You may make it become larger than a predetermined threshold value.
  • the combinations of the shapes of the extraction areas between the dimensions (dimension pairs) are different from each other.
  • between the dimensions (dimension pairs) adjacent to each other in the extraction area means that a part of the extraction area is in contact with each other.
  • the combinations of the shapes of the extraction regions between the dimensions (dimension pairs) may be different from each other.
  • the combination of the shapes of pairs of extraction regions between the dimensions (dimension pairs) may be different between all dimensions (dimension pairs) adjacent to the extraction region.
  • the combinations of the shapes of the extraction areas between the dimensions (dimension pairs) are different from each other. You may do it. Further, among at least one dimension (dimension pair) among the dimensions (dimension pairs) adjacent to each other in the extraction area, the combinations of the shapes of the extraction area pairs between the dimensions (dimension pairs) are different from each other. May be.
  • the following effects can be obtained by making the combinations of the shapes of the pairs of extraction regions between the dimensions where the extraction regions are adjacent to each other (dimension pairs) different from each other.
  • the correlation between extracted feature amounts tends to increase between dimensions in which extraction regions are adjacent to each other. This is because the distance between the two-dimensional extraction regions is close between the dimensions where the extraction regions are adjacent to each other (because the correlation between the close regions in the image is high). If the combination of the shapes of the pairs of extraction regions between these dimensions is the same, the correlation between the feature amounts is further increased, the redundancy of the entire feature vector is increased, and the discrimination ability is reduced.
  • the extraction regions between the dimensions (dimension pairs) are not adjacent to each other between the dimensions (dimension pairs) having the same combination of shapes of the pairs of extraction regions.
  • the extraction areas between those dimensions (dimension pairs) may not be adjacent to each other.
  • the extraction areas between those dimensions (dimension pairs) may not be adjacent to each other.
  • the extracted areas between those dimensions (dimension pairs) are adjacent to each other between dimensions (dimension pairs) that exceed a certain ratio. You may make it not.
  • the extracted areas between those dimensions (dimension pairs) are not adjacent to each other. It may be.
  • the combinations of the shapes of the extraction regions between the dimensions (dimension pairs) are different from each other.
  • the combinations of the shapes of the extraction areas between the dimensions (dimension pairs) may be different from each other.
  • the dimensions of the pairs of extraction areas between those dimensions (dimension pairs) are more than a certain ratio.
  • the combinations may be different from each other.
  • the dimensions (dimension pairs) where the extraction areas partially overlap at least one dimension (dimension pair), the combination of the shapes of the extraction areas between those dimensions (dimension pairs) They may be different from each other.
  • the following effects can be obtained by making the combinations of the shapes of the pairs of extraction regions between the dimensions in which the extraction regions partially overlap (dimension pairs) different from each other.
  • the correlation between the extracted feature amounts tends to increase between dimensions where the extraction regions partially overlap. This is because a region used for calculation of a feature amount is partially shared between dimensions where the extraction regions partially overlap. If the combination of the shapes of the pairs of extraction regions between these dimensions is the same, the correlation between the feature amounts is further increased, the redundancy of the entire feature vector is increased, and the discrimination ability is reduced. Therefore, by canceling the correlation between the extracted feature quantities due to the partial overlap of the extracted areas, the combination of the shape of the extracted areas between these dimensions is made different. Correlation can be reduced, redundancy as a whole feature vector can be reduced, and degradation of identification ability can be prevented.
  • the feature quantity (area feature quantity) extracted from it is affected. There is a high possibility that it will be different from the feature amount extracted from the area before falsification. If the combinations of the shapes of pairs of extracted regions between these dimensions are the same, the feature values (region feature values) extracted from the extracted regions of each dimension are affected in the same way, and the feature values are changed at the same time.
  • the extraction regions between the dimensions (dimension pairs) do not overlap between the dimensions (dimension pairs) in which the combinations of the shapes of the pairs of extraction regions are the same.
  • the extraction areas between those dimensions (dimension pairs) may not overlap.
  • the extraction regions between the dimensions (dimension pairs) may not overlap between all the dimensions (dimension pairs) having the same combination of shapes of the extraction regions.
  • the extracted areas between those dimensions (dimension pairs) overlap between dimensions (dimension pairs) that exceed a certain percentage. You may make it not.
  • the extraction areas between the dimensions (dimension pairs) do not overlap in at least one dimension (dimension pair). It may be.
  • the extraction area should be such that when the extraction areas of all dimensions are integrated, the area where the feature amount is not extracted in the image becomes small (that is, covers almost the entire screen of the image). .
  • the extraction areas of all dimensions are integrated, the area from which no feature value is extracted in the image becomes small (that is, covers almost the entire screen of the image). Since more information contained in can be reflected in the feature amount, the discrimination ability can be increased.
  • the extraction regions of all dimensions are integrated, it is desirable that the extraction regions have no bias and are obtained uniformly from the entire image. However, if there is a high probability that local processing such as telop superimposition is performed on a specific area, it is desirable that the extraction area is acquired avoiding that area. In addition, since there are generally no image characteristic portions in the peripheral area such as the edge of the image, it is desirable that the extraction area is acquired avoiding the peripheral area.
  • the size and relative position (distance, direction) of the extraction region follow a constant distribution (for example, a uniform distribution).
  • a constant distribution for example, a uniform distribution.
  • the relative position (distance, direction) follows a uniform distribution, so there is no bias in the distance and direction, and there is no concentration in a specific distance or direction, so more diversity can be taken. Because. Also, the closer the relative position, the greater the correlation between the regions, so in order to cancel it, it is desirable that the closer the relative position, the larger the difference in shape.
  • the dimension-specific extraction area information may be information in any format as long as the first extraction area and the second extraction area for each dimension can be uniquely identified.
  • the dimension-specific extraction area information indicates the same extraction area for images of any size and aspect ratio.
  • the information must be in a format that can be obtained.
  • the dimension-specific extraction area information describes the position and shape of the extraction area of an image having a certain size and aspect ratio (for example, an image having a horizontal width of 320 pixels and a vertical width of 240 pixels). Also good.
  • the image is first resized to the specified size and aspect ratio, and then the position of the extraction region described in the dimension-specific extraction region information -
  • the extraction area may be specified according to the shape.
  • the extraction area may be specified by converting the position / shape of the extraction area described in the dimension-specific extraction area information in accordance with an input image having an arbitrary size and aspect ratio. .
  • Information indicating each extraction area included in the dimension-specific extraction area information constitutes an extraction area for an image having a certain size and aspect ratio (for example, an image having a horizontal width of 320 pixels and a vertical width of 240 pixels), for example. It may be information describing a set of coordinate values of all pixels.
  • the information indicating each extraction area included in the dimension-specific extraction area information may be, for example, information in which the position / shape of the extraction area is described as a parameter for an image having a certain size and aspect ratio. For example, when the shape of the extraction region is a rectangle, information describing the coordinate values of the four corners of the rectangle may be used. For example, when the shape of the extraction region is a circle, the coordinate value and radius value of the center of the circle may be used.
  • the pseudo random number seed (seed) is used as the extraction area information for each dimension, starting from the seed within the extraction area acquiring means 2, generating a pseudo random number, and generating an extraction area having a different shape according to the random number.
  • the four corners of a quadrangle are determined according to a random number.
  • the generated extraction area Since the extraction area is determined based on a random number, the generated extraction area has a different shape for each dimension. If the seeds of the pseudo random numbers are the same, the same random number sequence is generated every time (for any image), so the same extraction region is reproduced for different images.
  • the extraction region acquisition unit 2 acquires information indicating the first extraction region and the second extraction region corresponding to the dimension supplied from the dimension determination unit 1 from the extraction region information supplied by dimension supplied as an input, and performs extraction. This is supplied to the area representative value calculation means 3.
  • the region feature quantity calculation unit 3 is supplied with an image from which an image identifier is to be extracted as an input. Is done.
  • the area feature quantity calculation means 3 includes a first area feature quantity calculation means 31 and a second area feature quantity calculation means 32.
  • the area feature quantity calculation means 3 uses the first area feature quantity calculation means 31 to indicate the first extraction area supplied from the extraction area acquisition means 2 for each dimension from the image supplied as input. Based on the above, the feature amount of the first extraction region is calculated as the first region feature amount and supplied to the comparison means 4.
  • the area feature quantity calculation unit 3 uses the second area feature quantity calculation unit 32 to calculate the second extraction area supplied from the extraction area acquisition unit 2 for each dimension from the image supplied as input. Based on the indicated information, the feature amount of the second extraction region is calculated as the second region feature amount and supplied to the comparison unit 4.
  • the region feature amount calculation unit 3 may select the dimension by size as necessary. Resizes the image to the specified size and aspect ratio of the extraction area information.
  • the region feature amount calculation means 3 calculates the region feature amount of each extraction region using the pixel value of the pixel group included in each extraction region.
  • the pixel value is a value of a signal possessed by each pixel of the image, and is a scalar amount or a vector amount.
  • the pixel value is a luminance value (scalar amount).
  • the pixel value is a vector quantity representing a color component.
  • the color image is an RGB image
  • the pixel value is a three-dimensional vector amount of an R component, a G component, and a B component.
  • the color image is a YCbCr image
  • the pixel value is a three-dimensional vector amount of a Y component, a Cb component, and a Cr component.
  • the calculation method in the extraction region of the dimension is constant (the same calculation method for any input image). As long as it is).
  • the region feature amount to be calculated may be a scalar amount or a vector amount.
  • the pixel value is a scalar quantity such as a luminance value
  • the area feature quantity is calculated as the average value, median value, mode value, maximum value, minimum value, etc. of the pixel values included in the extraction area. (Both are scalar quantities).
  • the pixel value at the position of P% from the lower order of the permutation sorted in ascending order is, for example, Y (floor (N ⁇ P / 100)), and this value is calculated as the region feature amount of the extraction region.
  • Floor () is a function for truncating after the decimal point.
  • the region feature amount calculated by applying this formula (Y (floor (N ⁇ P / 100))) to the luminance value of the pixel included in the extraction region is referred to as a “percentile luminance value feature amount”. I will call it.
  • the area feature quantity may be calculated by the method described above after first converting them into a scalar quantity by an arbitrary method.
  • the region feature amount may be calculated by the method described above after first converting them into a luminance value that is a scalar amount.
  • an average vector of pixel values included in the extraction region may be used as the region feature amount.
  • an arbitrary calculation such as edge detection or template matching may be performed on the extraction region, and the calculation result may be used as the region feature amount.
  • the region feature amount may be a two-dimensional vector quantity representing the edge direction (gradient direction).
  • the edge direction may be a scalar amount that represents the degree of similarity with a certain template.
  • a color distribution included in the extraction region, an edge direction distribution, and a histogram representing the edge intensity distribution may be calculated as the region feature amount (all are vector amounts).
  • various feature quantities defined in the international standard ISO / IEC 15938-3 that is, Dominant Color, Color Layout, Scalable Color, Color Structure, EdgeHistogram, HomogeneousTexture, Strain. Shape 3D, Parametric Motion, Motion Activity and the like may be used.
  • the comparison unit 4 compares, for each dimension, the first region feature amount supplied from the region feature amount calculation unit 3 and the second region feature amount, and the quantization result obtained by quantizing the comparison result Output the index.
  • the comparison unit 4 outputs a quantization index for each dimension, so that finally a quantization index vector composed of a plurality of dimensions of the quantization index is output.
  • the comparison unit 4 can arbitrarily quantize the first region feature amount and the second region feature amount.
  • the number of quantization indexes per dimension is also arbitrary.
  • the comparison unit 4 compares the magnitudes, and if the first area feature quantity is larger, the quantization index is +1, otherwise In this case, the quantization index may be set to ⁇ 1, and quantization may be performed into binary quantization indexes of +1 and ⁇ 1.
  • the quantization index Qn of dimension n can be calculated by the following equation.
  • FIG. 3 shows a more detailed configuration diagram of the comparison unit 4 when the comparison unit 4 performs the comparison / quantization based on the above-described equation 1.
  • the comparison means 4 is composed of a magnitude comparison means 41 and a quantization means 42.
  • the size comparison means 41 compares the values of the first region feature value and the second region feature value.
  • the comparison result is supplied to the quantization means 42. That is, the magnitude comparison means 41 compares the magnitudes of Vn1 and Vn2, and quantifies information indicating whether the comparison result is Vn1> Vn2 or Vn1 ⁇ Vn2 as a magnitude comparison result.
  • the generating means 42 To the generating means 42.
  • the quantization means 42 performs quantization according to the expression 1 based on the magnitude comparison result supplied from the magnitude comparison means 41, and outputs a quantization index. That is, when information indicating that the comparison result is Vn1> Vn2 is supplied, the quantization unit 42 is supplied with information indicating that the quantization index is +1 and the comparison result is Vn1 ⁇ Vn2. , The quantization index is set to -1, and the quantization index is output.
  • comparison / quantization method A the comparison / quantization method based on Equation 1 is referred to as comparison / quantization method A.
  • the comparison unit 4 determines the first region feature amount and the first feature amount. 2 is regarded as a quantization index 0 indicating that there is no difference.
  • the size is compared, and if the first region feature is larger, the quantum The quantization index may be +1, otherwise the quantization index may be ⁇ 1, and the quantization index may be quantized to +1, 0, ⁇ 1.
  • the quantization index Qn of dimension n can be calculated by the following equation. .
  • FIG. 4 shows a more detailed configuration diagram of the comparison unit 4 when the comparison unit 4 performs the comparison / quantization based on the above-described Expression 2.
  • the comparison unit 4 includes a difference value calculation unit 43 and a quantization unit 44.
  • the quantization means 44 is supplied with a threshold value, which is information (quantization boundary information) indicating a quantization boundary, which is defined in advance, as an input.
  • the difference value calculation unit 43 calculates a difference value between the value of the first region feature value and the value of the second region feature value.
  • the calculated difference value is supplied to the quantization means 44. That is, the difference value calculation unit 43 calculates Vn1 ⁇ Vn2 and supplies the value to the quantization unit 44.
  • the quantization means 44 is based on the difference value supplied from the difference value calculation means 43 and a threshold value that is information (quantization boundary information) indicating a predetermined quantization boundary supplied as an input. Quantization is performed according to 2 and a quantization index is output. That is, the quantizing means 42 is based on the value of Vn1 ⁇ Vn2 supplied from the difference value calculating means 41 and the threshold value th supplied as input, in the case of
  • the quantization index is +1,
  • comparison / quantization method B the comparison / quantization method based on Equation 2 is referred to as comparison / quantization method B.
  • the quantization is performed into three values based on the difference value, but it may be quantized into a larger number (of levels) of quantization indexes according to the magnitude of the difference value.
  • the comparison means 4 has the configuration shown in FIG. 4, and the quantization means 44 has a plurality of threshold values as information (quantization boundary information) indicating the boundaries of quantization of each level specified in advance. Are supplied as inputs.
  • a comparison / quantization method that quantizes a quantization index of a plurality of levels of four or more levels based on the difference value and a plurality of threshold values supplied as input is referred to as a comparison / quantization method C. To.
  • the difference between the first region feature value and the second region feature value is small (below the prescribed threshold value)
  • a quantization index indicating that there is no difference is introduced.
  • the feature quantity (quantization index) of the pair of extracted areas with a small difference in area feature quantity is more stable, that is, more robust to various modification processes and noise. can do. Therefore, the difference in features between local regions is small overall, stable even for flat images with little overall change (for example, images of blue sky), that is, robust to various modification processes and noise, An image identifier (quantization index vector) can be output.
  • the comparison unit 4 may first convert the vector quantity into a scalar quantity by an arbitrary method and then perform quantization by the method described above ( This comparison / quantization method will be referred to as comparison / quantization method D). Further, for example, a difference vector that is a difference between the vector of the first extraction region and the vector of the second extraction region may be calculated, and the quantization index may be calculated by vector quantization of the difference vector. In this case, for example, a representative vector (centroid vector or the like) for each predetermined quantization index is supplied, and the quantization index having the largest similarity (the smallest distance) between the representative vector and the difference vector is used.
  • comparison / quantization method E this comparison / quantization method will be referred to as comparison / quantization method E.
  • this comparison / quantization method when the norm of the difference vector is equal to or smaller than a predetermined threshold, the difference between the first area feature quantity and the second area feature quantity is not found.
  • a quantization index 0 indicating that there is no difference
  • a quantization index indicating that there is no difference may be introduced.
  • the quantization index vectors output in the present invention are collated (the quantization index vector extracted from one image is compared with the quantization index vector extracted from another image, the images are identical).
  • the identity measure can be compared with a threshold to determine the identity of the image.
  • the comparison unit 4 calculates the quantization index based on Equation 2
  • the identity measure can be calculated as follows.
  • the quantization index vectors of two images are compared between corresponding dimensions, and the number of dimensions whose “quantization index is not 0” is calculated (this value is assumed to be A).
  • the number of dimensions with the same quantization index is calculated in dimensions where “both quantization indexes are not 0” (this value is B).
  • the similarity is calculated as B / A.
  • A 0 (that is, when the quantization index is 0 for all dimensions)
  • the similarity is set to a specified numerical value (for example, 0.5).
  • the extraction area acquisition means 2 acquires information indicating the first extraction area and the second extraction area of the dimension n from the dimension-specific extraction area information supplied as input, and the area feature amount calculation means 3 (Step A2).
  • the area feature quantity calculating means 3 calculates a first area feature quantity and a second area feature quantity of dimension n from the image supplied as input, and supplies them to the comparison means 4 (step A3). ).
  • the comparison unit 4 compares the first region feature value of dimension n with the second region feature value, quantizes the comparison result, and outputs a quantization index (step A4).
  • step A5 it is determined whether or not the output of the quantization index has been completed for all dimensions (that is, whether n ⁇ N is true or false) (step A5).
  • the process ends.
  • the output of the quantization index is not completed for all dimensions (that is, when n ⁇ N is true)
  • the process proceeds to step A6.
  • extraction processing is performed in order from dimension 1 to dimension N, but the order is not limited to this and may be arbitrary.
  • extraction processing for a plurality of dimensions may be performed in parallel.
  • the first effect is that it is possible to increase the identification ability, which is the degree to which different images can be identified, of image identifiers composed of feature vectors composed of a plurality of dimensions. In particular, this effect is remarkable for an image having a large correlation between local regions of the image.
  • the second effect is that the discrimination ability does not deteriorate even for an image in which signals are concentrated at a specific frequency.
  • the shape of the region from which features are extracted differs between dimensions (the shape of the region is diverse), so even for images where signals are concentrated at a specific frequency, This is because there is no difference in feature quantity between (many) extraction area pairs (dimensions), and it is difficult for occurrence of a reduction in discrimination ability.
  • the second embodiment of the present invention is different in that the comparison means 4 in the first embodiment shown in FIG. 1 is replaced with a comparison means 4A shown in detail in FIG. Other than the comparison unit 4A is the same as that of the first embodiment, and thus the description thereof is omitted here.
  • the comparison unit 4 ⁇ / b> A includes a difference value calculation unit 43, a quantization boundary determination unit 45, and a quantization unit 44.
  • the difference value calculation unit 43 calculates a difference value between the first region feature amount and the second region feature amount supplied from the region feature amount calculation unit 3 for each dimension, and the quantization boundary determination unit 45 To the quantizing means 44.
  • the difference value is, for example, a scalar obtained by subtracting the second region feature amount (or vice versa) from the first region feature amount when the region feature amount is a scalar amount (for example, an average value of luminance values). Amount. Further, when the region feature amount is a vector amount, for example, each vector may be converted into a scalar amount by an arbitrary method, and then the difference value of the scalar amount may be obtained. When the region feature amount is a vector amount, a difference vector between the first region feature amount and the second region feature amount may be used as a difference value (vector amount).
  • the quantization boundary determination unit 45 determines the quantization boundary based on the distribution of the difference values of all the dimensions. Then, the information on the determined quantization boundary is supplied to the quantization means 44.
  • the distribution of difference values in all dimensions is the frequency (probability) of occurrence with respect to the difference value (or difference vector).
  • the determination of the quantization boundary means that, when the difference value is quantized, a parameter to be assigned exclusively to the quantization index is determined without omission.
  • the difference value is a scalar quantity, for example, a value range (that is, a threshold value) for each quantization index (quantization level) is determined, and the value range (threshold value) is supplied to the quantization unit 43 as quantization boundary information.
  • the difference value is a vector quantity
  • a parameter for performing vector quantization for example, a representative vector (such as a centroid vector) of each quantization index is determined, and quantization means is used as information on the quantization boundary. 44.
  • the quantization range (threshold value) may be determined so that the ratios of the respective quantization indexes to all dimensions are uniform.
  • a constant ⁇ is used.
  • the quantization index is +1
  • the center point of the difference value distribution (the point at which the integrated values of the left and right distributions are equal) is determined as the quantization threshold ⁇ so that the ratios of the quantization indexes +1 and ⁇ 1 are equal. Also good.
  • the difference value is a vector quantity
  • the ratios of the respective quantization indexes to all dimensions are made uniform based on the distribution of the difference vectors of all dimensions.
  • a vector space region allocated to each quantization index may be determined, or a representative vector (such as a centroid vector) of each quantization index when performing vector quantization may be determined.
  • a representative vector such as a centroid vector
  • the entropy can be increased by equalizing the ratio of the quantization index to all dimensions (that is, eliminating the bias of the quantization index). can do.
  • the quantization boundary determination unit 45 determines the quantization boundary so that the ratio of the quantization index to all dimensions is uniform, and the quantization unit 44 performs quantization based on the comparison / quantization method. This is referred to as comparison / quantization method F.
  • the quantization boundary determination unit 45 performs the difference when the difference value is a scalar quantity and the ternary quantization is performed according to the above-described equation 2 (quantization index is +1, 0, ⁇ 1).
  • Threshold value th when quantizing to quantization index 0 indicating that there is no is determined based on the distribution of difference values in all dimensions.
  • the threshold th may be supplied to the quantization means 44 (in the comparison means 4 of FIG. 4 of the first embodiment, this threshold th is defined in advance).
  • the absolute values of the difference values of all dimensions are calculated, the absolute values of the calculated difference values are sorted, and a specified ratio (note that the specified ratio is supplied as an input, for example, from the upper or lower order)
  • the comparison / quantization method will be referred to as comparison / quantization method G).
  • the threshold th may be determined so that the ratios of the quantization indexes of +1, 0, and -1 approach each other instead of the prescribed ratio (this comparison / quantization method is compared with the comparison / quantization method H). I will call it).
  • the comparison / quantization method H corresponds to a specific example of the comparison / quantization method F in the case of following Formula 2.
  • Floor () is a function for truncating after the decimal point.
  • the method in the present embodiment can be compared with the case in which the comparison unit 4 in the first embodiment has the configuration of FIG.
  • a predetermined threshold th is supplied as an input, whereas the above-described method in the second embodiment is performed by the quantization boundary determining unit 45.
  • the threshold th is adaptively calculated for the image. As described above, the threshold th is fixed in the first embodiment, and the threshold th is adaptively calculated for the image in the second embodiment.
  • the threshold th By calculating the threshold th adaptively for the image, the dimension value of the feature vector is biased toward a specific quantization index (when the specific quantization index is Since the appearance probability is high (especially for an image with few undulations), the discrimination ability can be increased. For example, when the fixed threshold th in the first embodiment is used, an image with few undulations has the majority (or all dimensions) of the feature vector having a quantization index of 0. On the other hand, when the adaptive threshold value th in the second embodiment is used, the threshold value th is automatically adjusted to a small value for an image with few undulations. There is no situation where the index becomes zero.
  • the quantization unit 44 performs quantization for each dimension based on the difference value for each dimension supplied from the difference value calculation unit 43 and the quantization boundary information supplied from the quantization boundary determination unit 45. , Output quantization index.
  • the quantization means 44 is meaningless if the quantization boundary information outputted from the quantization boundary determination means 45 is ignored, and the quantization boundary determination means 45 determines the quantization boundary. It is necessary to follow the quantization method assumed in the above.
  • the extraction region acquisition unit 2 acquires information indicating the first extraction region and the second extraction region of the dimension n from the extraction region information by dimension supplied as input, and the extraction region representative value calculation unit 3 (step B2).
  • the extracted area representative value calculation unit 3 calculates the first area feature quantity and the second area feature quantity of the dimension n from the image supplied as input, and supplies them to the difference value calculation unit 43. (Step B3).
  • the difference value calculation unit 43 calculates a difference value between the first region feature amount and the second region feature amount of the dimension n, and supplies the difference value to the quantization boundary determination unit 45 and the quantization unit 44. (Step B4).
  • step B5 it is determined whether or not the processing up to the calculation of difference values for all dimensions has been completed (that is, whether n ⁇ N is true or false) (step B5).
  • the process proceeds to step B7. If the processing for all dimensions has not been completed (that is, if n ⁇ N is true), the process proceeds to step B6.
  • extraction processing is performed in order from dimension 1 to dimension N, but the order is not limited to this and may be arbitrary.
  • the quantization boundary determination means 45 performs quantization based on the distribution of the difference values of all dimensions.
  • the boundary is determined, and information on the determined quantization boundary is supplied to the quantization means 44 (step B7).
  • the quantization means 44 performs quantization based on the difference value of the dimension n and the quantization boundary supplied from the quantization boundary determination means 45, and outputs a quantization index (step B9).
  • step B10 it is determined whether or not the output of quantization indexes for all dimensions has been completed (that is, whether n ⁇ N is true or false) (step B10).
  • the process ends.
  • the output of the quantization index for all dimensions is not completed (that is, when n ⁇ N is true)
  • the process proceeds to step B11.
  • quantization processing is performed in order from dimension 1 to dimension N, but the order is not limited to this and may be arbitrary.
  • the quantization boundary is calculated adaptively (dynamically) with respect to the image. Different.
  • the dimension value of the feature vector is set to a specific quantum for a specific image (for example, a flat image with few undulations).
  • the index is biased toward a specific index (the probability of occurrence of a specific quantization index is high) (entropy is low), and a problem arises in that the discrimination ability of these images is reduced.
  • the quantization boundary is adaptively (dynamically) calculated for the image, so that the dimension value of the feature vector for any image becomes Since it can suppress that it is biased to a specific quantization index (the appearance probability of a specific quantization index is high), the identification capability can be increased.
  • a region feature quantity calculation method acquisition unit 5 is added to the configuration of the first embodiment shown in FIG. Is different in that it is replaced with a region feature amount calculation unit 3A having first and second region feature amount calculation units 31A and 32A. Since the other configuration is the same as that of the first embodiment, the description thereof is omitted here. In addition, although it demonstrated as a combination with 1st Embodiment here, the combination with 2nd Embodiment may be sufficient.
  • the area feature quantity calculation method acquisition means 5 is supplied with the dimensions from the dimension determining means 1 and the dimension-specific area feature quantity calculation method information.
  • the area feature value calculation method information for each dimension is information indicating a calculation method of the area feature value in the dimension, which is associated with each dimension of the feature vector defined in advance. Is a necessary condition.
  • the difference in the region feature amount calculation method includes a case where different parameters (threshold value, etc.) are applied to the same procedure.
  • the region feature amount calculation method is, for example, various methods described in the description of the region feature amount calculation means 3 of the first embodiment, and parameters associated therewith.
  • the region feature amount calculation method for each dimension indicated by the region feature amount calculation method information for each dimension is such that at least one pair of different dimension of the region feature amount calculation method exists in all dimensions of the feature vector. It is a condition. It is desirable that there are more dimensions in which the region feature quantity calculation methods are different from each other. This is because as the area feature quantity calculation method has more different dimensions, the correlation between more dimensions of the feature vector becomes smaller and the discrimination ability becomes higher. For example, the region feature amount calculation method may be different between all dimensions of the feature vector.
  • the format of the information indicating the region feature value calculation method for each dimension may be any format as long as the method for calculating the region feature value is uniquely specified.
  • FIG. 9 shows an example of a region feature amount calculation method for each dimension.
  • the region feature quantity calculation method differs between dimensions.
  • the feature quantity of the scalar quantity and the vector quantity may be mixed (the first, third, fifth, sixth, eighth, ninth, tenth and twelfth dimensions are the scalar quantity, the second quantity). 4, 7 and 11 dimensions are vector quantities).
  • the region feature quantity calculation method between dimensions (dimension pairs) is different between dimensions (dimension pairs) in which extraction regions exist in the vicinity of each other.
  • the region feature quantity calculation method between dimensions (dimension pairs) may be different between dimensions (dimension pairs) in which the extraction areas exist in the vicinity of each other.
  • the region feature quantity calculation method between the dimensions (dimension pairs) may be different between all the dimensions (dimension pairs) in which the extraction areas exist in the vicinity.
  • the area feature value calculation method between the dimensions (dimension pairs) is different between dimensions (dimension pairs) of a certain ratio or more. It may be.
  • at least one dimension (dimension pair) is set so that the region feature quantity calculation method between those dimensions (dimension pairs) is different. Also good.
  • the dimension pair in which the extraction regions exist in the vicinity of each other has a distance between the two-dimensional extraction regions (the first extraction region and the second extraction region) calculated by a certain predetermined method.
  • the distance between the two-dimensional extraction regions (the first extraction region and the second extraction region) is, for example, the sum of all the pixels included in each dimension extraction region (the first extraction region and the second extraction region). You may calculate as an average value of distance between. Further, for example, it may be calculated as the minimum value of the distance between all the pixels included in each dimension extraction region (first extraction region and second extraction region).
  • centroid coordinates of the extraction areas of each dimension are calculated (the centroid coordinates of the first extraction area and the centroid coordinates of the second extraction area), and the distances between them are calculated (a total of four distances are calculated).
  • the distance between the two-dimensional extraction regions may be defined as an average value or a minimum value thereof. However, the distance between the two-dimensional extraction regions is not limited to these calculation methods.
  • the following effects can be obtained by making the region feature amount calculation method between dimensions different between the dimensions (dimension pairs) in which the extraction areas are close to each other.
  • the correlation between the extracted feature amounts tends to increase between the dimensions in which the extraction regions are close to each other. If the region feature quantity calculation method between these dimensions is the same, the correlation between the feature quantities is further increased, the redundancy of the entire feature vector is increased, and the discrimination ability is lowered. Therefore, by canceling the increase in the correlation between the extracted feature quantities due to the presence of the extracted areas in the vicinity, the region feature quantity calculation method between these dimensions is different, so that Correlation can be reduced, redundancy as a whole feature vector can be reduced, and degradation of identification ability can be prevented.
  • the distance between the extracted regions between dimensions (dimension pairs) between dimensions (dimension pairs) having the same region feature quantity calculation method is larger.
  • the distance between these dimensions (dimension pairs) may be larger than a predetermined threshold.
  • the distance between these dimensions (dimension pairs) may be larger than a predetermined threshold.
  • at least one dimension (dimension pair) has a distance between those dimensions (dimension pairs) from a predetermined threshold. May also be increased.
  • the region feature quantity calculation method between dimensions (dimension pairs) be different between dimensions (dimension pairs) adjacent to each other.
  • between the dimensions (dimension pairs) adjacent to each other in the extraction area means that a part of the extraction area is in contact with each other.
  • the region feature quantity calculation method between the dimensions (dimension pairs) between the adjacent dimensions (dimension pairs) may be different.
  • the region feature quantity calculation method between the dimensions (dimension pairs) may be different between all the dimensions (dimension pairs) adjacent to the extraction area.
  • the area feature value calculation method between the dimensions (dimension pairs) between the dimensions (dimension pairs) exceeding a certain predetermined ratio may be different. Good.
  • the region feature quantity calculation method between the dimensions (dimension pairs) may be different between at least one dimension (dimension pair).
  • the following effects can be obtained by making the region feature quantity calculation method between dimensions different between adjacent dimensions (dimension pairs).
  • the correlation between extracted feature amounts tends to increase between dimensions in which extraction regions are adjacent to each other. This is because the distance between the two-dimensional extraction regions is close between the dimensions where the extraction regions are adjacent to each other (because the correlation between the close regions in the image is high). If the region feature quantity calculation method between these dimensions is the same, the correlation between the feature quantities is further increased, the redundancy of the entire feature vector is increased, and the discrimination ability is lowered. Therefore, by canceling the correlation between the extracted feature quantities due to the existence of the extracted areas adjacent to each other, the feature quantity calculation method between the dimensions is made different so that the feature quantities are different. Can be reduced, the redundancy of the entire feature vector can be reduced, and the discrimination ability can be prevented from lowering.
  • the extracted regions between dimensions (dimension pairs) are not adjacent to each other between dimensions (dimension pairs) having the same region feature value calculation method.
  • the extraction areas between those dimensions (dimension pairs) may not be adjacent to each other.
  • the extraction regions between these dimensions (dimension pairs) may not be adjacent to each other between all dimensions (dimension pairs) having the same region feature quantity calculation method.
  • the extracted areas between these dimensions (dimension pairs) should not be adjacent to each other between dimensions (dimension pairs) of a certain ratio or more. May be.
  • the extraction areas between these dimensions (dimension pairs) should not be adjacent to each other between at least one dimension (dimension pair). Good.
  • the region feature quantity calculation method between dimensions (dimension pairs) differs between dimensions (dimension pairs) where the extraction areas partially overlap.
  • the region feature amount calculation method between dimensions (dimension pairs) may be different between dimensions (dimension pairs) in which extraction regions partially overlap.
  • the region feature amount calculation method between the dimensions (dimension pairs) may be different between all the dimensions (dimension pairs) where the extraction areas partially overlap.
  • the area feature value calculation method between the dimensions (dimension pairs) is different between dimensions (dimension pairs) of a predetermined ratio or more. You may do it.
  • the region feature quantity calculation method between the dimensions (dimension pairs) is different between at least one dimension (dimension pair). May be.
  • the following effects can be obtained by making the region feature quantity calculation method between the dimensions in which the extracted regions partially overlap (dimension pairs) different.
  • the correlation between the extracted feature amounts tends to increase between dimensions where the extraction regions partially overlap. This is because a region used for calculation of a feature amount is partially shared between dimensions where the extraction regions partially overlap. If the region feature quantity calculation method between these dimensions is the same, the correlation between the feature quantities is further increased, the redundancy of the entire feature vector is increased, and the discrimination ability is lowered. Therefore, the correlation of feature quantities is made different by using different region feature quantity calculation methods between these dimensions so as to cancel the increase in correlation between extracted feature quantities due to partial overlap of extracted areas. Can be reduced, the redundancy of the entire feature vector can be reduced, and the degradation of the discrimination ability can be prevented.
  • the region feature quantity calculation method between dimensions (dimension pairs) in which extracted regions partially overlap each other different, the following effects can be obtained. It is possible to enhance the resistance against image alteration (especially malicious alteration) by making the region feature value calculation method different between the dimensions (dimension pairs) where the extracted areas partially overlap. There is an effect that. Malicious tampering with an image is performed, for example, by changing pixel values in a partial area of the image, for example, for the purpose of degrading the accuracy of identity determination based on an image identifier (feature vector) extracted from the image. .
  • the feature quantity (area feature quantity) extracted from it is affected. There is a high possibility that it will be different from the feature amount extracted from the area before falsification. If the area feature value calculation method between these dimensions is the same, the feature value (area feature value) extracted from the extraction area of each dimension is affected in the same way, and the feature value may be changed at the same time. , The robustness is low. If the area feature value calculation method between these dimensions is different, even if tampering is performed on the overlapping area, the feature value is similarly affected, and the possibility that the feature value will be changed at the same time is reduced.
  • the extracted regions between dimensions (dimension pairs) do not overlap between dimensions (dimension pairs) having the same area feature quantity calculation method.
  • the extracted regions between dimensions (dimension pairs) may not overlap between dimensions (dimension pairs) having the same region feature quantity calculation method.
  • the extraction regions between the dimensions (dimension pairs) may not overlap between all dimensions (dimension pairs) having the same region feature quantity calculation method.
  • the extracted areas between those dimensions (dimension pairs) should not overlap. May be.
  • at least one dimension (dimension pair) between the dimensions (dimension pairs) may not be overlapped. Good.
  • the combinations of the shapes of the extraction region pairs are different between dimensions (dimension pairs) in which the region feature quantity extraction method is the same.
  • the region feature quantity extraction method is different between dimensions (dimension pairs) in which the combinations of the shapes of the extraction area pairs are the same.
  • the combination of the shape of the pair of extraction regions may be different between dimensions (dimension pairs) having the same region feature extraction method.
  • the region feature quantity extraction method may be different between dimensions (dimension pairs) in which the combinations of the shapes of the pairs of extraction regions are the same.
  • the combination of the shape of the pair of extraction regions may be different between all dimensions (dimension pairs) having the same region feature quantity extraction method.
  • the region feature extraction method may be different between all dimensions (dimension pairs) in which the combination of the shapes of the extraction region pairs is the same. Further, for example, among the dimensions (dimension pairs) having the same region feature quantity extraction method, the combinations of the shapes of the extracted area pairs may be different between dimensions (dimension pairs) of a predetermined ratio or more. . Further, the region feature quantity extraction method may be different between dimensions (dimension pairs) that are equal to or higher than a predetermined ratio among dimensions (dimension pairs) in which the combinations of the shapes of the pairs of extraction areas are the same. Further, for example, among the dimensions (dimension pairs) having the same region feature quantity extraction method, the combination of the shapes of the extraction region pairs may be different between at least one dimension (dimension pair). In addition, the region feature quantity extraction method may be different between at least one dimension (dimension pair) among dimensions (dimension pairs) having the same combination of shapes of pairs of extraction areas.
  • the area feature quantity calculation method acquisition means 5 acquires information indicating the area feature quantity calculation method corresponding to the dimension supplied from the dimension determination means 1 from the dimension-specific area feature quantity calculation method information supplied as input. It is supplied to the feature amount calculation means 3A.
  • the area feature quantity calculation means 3A is based on information indicating the first extraction area and the second extraction area supplied from the extraction area acquisition means 2 for each dimension from the image supplied as input. According to the information indicating the region feature amount calculation method supplied from the calculation method acquisition unit 5, the feature amount of the first extraction region and the feature amount of the second extraction region are respectively set to the first region feature amount and the second region feature amount. And is supplied to the comparison means 4.
  • the dimension of information indicating the extraction area to be supplied and the dimension of information indicating the area feature quantity calculation method are synchronized.
  • the dimension (number) of the feature vector is represented by n, and there are a total of N dimensions from 1 to N.
  • the extraction region acquisition unit 2 acquires information indicating the first extraction region and the second extraction region of the dimension n from the dimension-specific extraction region information supplied as input, and the region feature amount calculation unit 3A. (Step C2).
  • the region feature amount calculation method acquisition unit 5 acquires information indicating the region feature amount calculation method corresponding to the dimension n from the dimension-specific region feature amount calculation method information supplied as an input, and the region feature amount calculation unit Supply to 3A (step C3).
  • the area feature quantity calculating means 3A calculates a first area feature quantity and a second area feature quantity of dimension n from the image supplied as input, and supplies them to the comparison means 4 (step C4). ).
  • the comparison unit 4 compares the first area feature quantity of the dimension n with the second area feature quantity, quantizes the comparison result, and outputs a quantization index (step C5).
  • it is determined whether or not the output of the quantization index has been completed for all dimensions step C6). If the output of the quantization index is completed for all dimensions, the process ends. If the output of the quantization index has not been completed for all dimensions, the process proceeds to step C7.
  • extraction processing is performed in order from dimension 1 to dimension N, but the order is not limited to this and may be arbitrary.
  • extraction processing for a plurality of dimensions may be performed in parallel.
  • step C2 and step C3 may be reversed.
  • the comparison method acquisition means 6 is added to the configuration of the first embodiment shown in FIG. 1, and the comparison means 4 is replaced with the comparison means 4B. It is different in point. Since the other configuration is the same as that of the first embodiment, the description thereof is omitted here. In addition, although it demonstrated as a combination with 1st Embodiment here, the combination with 2nd Embodiment and 3rd Embodiment may be sufficient.
  • the comparison method acquisition means 6 is supplied with the dimensions from the dimension determination means 1 and the comparison method information for each dimension.
  • the comparison / quantization method information for each dimension is information indicating a method for performing quantization by comparing region feature quantities in the dimensions associated with each dimension of the feature vector defined in advance. It is an essential condition that the comparison and quantization methods are different.
  • different comparison / quantization methods include the case where different parameters (threshold, number of quantization indexes, etc.) are applied to the same procedure.
  • the comparison / quantization method refers to, for example, various comparison / quantization methods described in the description of the comparison unit 4 of the first embodiment, parameters associated therewith (threshold value, number of quantization indexes, etc.), These are various comparison / quantization methods described in the description of the comparison means 4A of the second embodiment, and accompanying parameters (threshold, number of quantization indexes, etc.).
  • the comparison / quantization method for each dimension indicated by the comparison / quantization method information for each dimension is such that at least one pair of different dimensions of the comparison / quantization method exists in all dimensions of the feature vector. It is a condition. The more dimensions the comparison and quantization methods are different from each other, the better. This is because as the comparison / quantization method has more different dimensions, the correlation between more dimensions of the feature vector becomes smaller and the discrimination ability becomes higher. For example, the comparison / quantization method may be different between all dimensions of the feature vector.
  • the format of the information indicating the comparison / quantization method for each dimension may be any format as long as the method of comparing and quantizing the region feature quantities is uniquely specified.
  • FIG. 12 shows an example of a comparison / quantization method for each dimension.
  • the comparison / quantization method differs between dimensions.
  • different parameters may be set by the same comparison / quantization method as in the third, fifth, and twelfth dimensions.
  • the example of the comparison / quantization method for each dimension shown in FIG. 12 corresponds to the example of the region feature value calculation method for each dimension shown in FIG.
  • the comparison / quantization method of scalar quantities is shown
  • the comparison / quantization method of vector quantities is shown as an example for region feature quantities of vector quantities.
  • the comparison / quantization method between dimensions (dimension pairs) differs between dimensions (dimension pairs) in which extraction regions exist in the vicinity of each other.
  • the comparison / quantization method between dimensions (dimension pairs) may be different between dimensions (dimension pairs) in which extraction regions exist in the vicinity of each other.
  • the comparison / quantization method between these dimensions (dimension pairs) may be different between all dimensions (dimension pairs) in which the extraction regions exist in the vicinity of each other.
  • the comparison / quantization method between these dimensions (dimension pairs) is different between dimensions (dimension pairs) of a certain ratio or more. It may be.
  • at least one dimension (dimension pair) is different in the comparison / quantization method between those dimensions (dimension pairs). Also good.
  • the dimension pair in which the extraction regions exist in the vicinity of each other has a distance between the two-dimensional extraction regions (the first extraction region and the second extraction region) calculated by a certain predetermined method.
  • the distance between the two-dimensional extraction regions (the first extraction region and the second extraction region) is, for example, the sum of all the pixels included in each dimension extraction region (the first extraction region and the second extraction region). You may calculate as an average value of distance between. Further, for example, it may be calculated as the minimum value of the distance between all the pixels included in each dimension extraction region (first extraction region and second extraction region).
  • centroid coordinates of the extraction areas of each dimension are calculated (the centroid coordinates of the first extraction area and the centroid coordinates of the second extraction area), and the distances between them are calculated (a total of four distances are calculated).
  • the distance between the two-dimensional extraction regions may be defined as an average value or a minimum value thereof. However, the distance between the two-dimensional extraction regions is not limited to these calculation methods.
  • the comparison / quantization method between dimensions By making the comparison / quantization method between dimensions different between dimensions (dimension pairs) in which the extraction regions exist in the vicinity, the following effects can be obtained.
  • the correlation between the extracted feature amounts tends to increase between the dimensions in which the extraction regions are close to each other. If the comparison / quantization method between these dimensions is the same, the correlation between the feature amounts is further increased, the redundancy of the entire feature vector is increased, and the discrimination ability is lowered. Therefore, by canceling the correlation between the extracted feature quantities due to the presence of the extracted regions in the vicinity, the comparison / quantization method between these dimensions is different, so that Correlation can be reduced, redundancy as a whole feature vector can be reduced, and degradation of identification ability can be prevented.
  • the distance between the extracted regions between dimensions (dimension pairs) between dimensions (dimension pairs) having the same comparison / quantization method is larger.
  • the distance between these dimensions (dimension pairs) may be larger than a predetermined threshold.
  • the distance between these dimensions (dimension pairs) may be larger than a predetermined threshold.
  • at least one dimension (dimension pair) has a distance between those dimensions (dimension pairs) from a predetermined threshold. May also be increased.
  • the comparison / quantization method between dimensions (dimension pairs) differs between dimensions (dimension pairs) adjacent to each other in the extraction region.
  • between the dimensions (dimension pairs) adjacent to each other in the extraction area means that a part of the extraction area is in contact with each other.
  • the comparison / quantization method between dimensions (dimension pairs) may be different between dimensions (dimension pairs) adjacent to the extraction region.
  • the comparison / quantization method between these dimensions (dimension pairs) may be different between all the dimensions (dimension pairs) adjacent to the extraction region.
  • the comparison / quantization method between these dimensions (dimension pairs) is different between dimensions (dimension pairs) of a predetermined ratio or more. Good.
  • the comparison / quantization method between the dimensions (dimension pairs) may be different between at least one dimension (dimension pair) among the dimensions (dimension pairs) adjacent to the extraction region.
  • the following effects can be obtained by making the comparison / quantization methods between dimensions different between dimensions (dimension pairs) where the extraction regions are adjacent to each other.
  • the correlation between extracted feature amounts tends to increase between dimensions in which extraction regions are adjacent to each other. This is because the distance between the two-dimensional extraction regions is close between the dimensions where the extraction regions are adjacent to each other (because the correlation between the close regions in the image is high). If the comparison / quantization method between these dimensions is the same, the correlation between the feature amounts is further increased, the redundancy of the entire feature vector is increased, and the discrimination ability is lowered. Therefore, by canceling the correlation between the extracted feature quantities due to the existence of the extracted regions adjacent to each other, the comparison / quantization method between these dimensions is made different so that the feature quantities are different. Can be reduced, the redundancy of the entire feature vector can be reduced, and the discrimination ability can be prevented from lowering.
  • the extraction regions between dimensions (dimension pairs) are not adjacent to each other between dimensions (dimension pairs) having the same comparison / quantization method.
  • the extraction regions between those dimensions (dimension pairs) may not be adjacent to each other.
  • the extraction regions between those dimensions (dimension pairs) may not be adjacent to each other.
  • between dimensions (dimension pairs) that have the same comparison / quantization method (dimension pairs) do not allow the extraction areas between those dimensions (dimension pairs) to be adjacent to each other. May be.
  • at least one dimension (dimension pair) may not be adjacent to the extraction area between those dimensions (dimension pair). Good.
  • the comparison / quantization method between dimensions (dimension pairs) differs between dimensions (dimension pairs) in which the extraction regions partially overlap.
  • a comparison / quantization method between dimensions (dimension pairs) may be different between dimensions (dimension pairs) in which extraction regions partially overlap.
  • the comparison / quantization method between these dimensions (dimension pairs) may be different between all dimensions (dimension pairs) in which the extraction regions partially overlap.
  • the comparison / quantization method between the dimensions (dimension pairs) is different between dimensions (dimension pairs) of a predetermined ratio or more. You may do it.
  • the comparison / quantization method between the dimensions (dimension pairs) is different between at least one dimension (dimension pair). May be.
  • the correlation between the extracted feature amounts tends to increase between dimensions where the extraction regions partially overlap. This is because a region used for calculation of a feature amount is partially shared between dimensions where the extraction regions partially overlap. If the comparison / quantization method between these dimensions is the same, the correlation between the feature amounts is further increased, the redundancy of the entire feature vector is increased, and the discrimination ability is lowered. Therefore, the correlation of feature quantities is different by making the comparison / quantization method between these dimensions different so that the correlation of extracted feature quantities increases due to partial overlap of extracted regions. Can be reduced, the redundancy of the entire feature vector can be reduced, and the degradation of the discrimination ability can be prevented.
  • the feature quantity (area feature quantity) extracted from it is affected.
  • the feature values (region feature values) extracted from the extraction regions of each dimension are affected in the same way, and the feature values may be changed at the same time. , The robustness is low. If the comparison and quantization methods between these dimensions are different, even if tampering is performed on the overlapping area, the feature amount will be affected in the same way and the possibility that the feature amount will be changed at the same time is reduced.
  • the extraction regions between dimensions (dimension pairs) do not overlap between dimensions (dimension pairs) having the same comparison / quantization method.
  • the extraction regions between those dimensions (dimension pairs) may not overlap.
  • the extraction regions between these dimensions (dimension pairs) may not overlap.
  • the extraction areas between those dimensions (dimension pairs) should not overlap between dimensions (dimension pairs) that exceed a certain percentage. May be.
  • at least one dimension (dimension pair) between the dimensions (dimension pairs) may not overlap. Good.
  • the combinations of the shapes of the extraction region pairs be different between dimensions (dimension pairs) having the same comparison / quantization method.
  • the comparison / quantization method is different between dimensions (dimension pairs) in which the combinations of the shapes of pairs of extraction regions are the same.
  • the combination of the shapes of the extraction region pairs may be different between dimensions (dimension pairs) having the same comparison / quantization method.
  • the comparison / quantization method may be different between dimensions (dimension pairs) in which the combinations of the shapes of the extraction region pairs are the same. In this case, for example, the combination of the shape of the pair of extraction regions may be different between all dimensions (dimension pairs) having the same comparison / quantization method.
  • the comparison / quantization method may be different between all dimensions (dimension pairs) having the same combination of shapes of the pairs of extraction regions. Further, for example, among the dimensions (dimension pairs) having the same comparison / quantization method, the combinations of the shapes of the pairs of the extraction regions may be different between dimensions (dimension pairs) of a certain ratio or more. . Further, the comparison / quantization method may be different between dimensions (dimension pairs) that are equal to or greater than a predetermined ratio among dimensions (dimension pairs) having the same combination of shapes of pairs of extraction regions. Further, for example, among the dimensions (dimension pairs) having the same comparison / quantization method, the combination of the shapes of the extraction region pairs may be different between at least one dimension (dimension pair). In addition, the comparison / quantization method may be different between at least one dimension (dimension pair) among dimensions (dimension pairs) having the same combination of shapes of pairs of extraction regions.
  • the region feature quantity calculation method is different between dimensions (dimension pairs) having the same comparison / quantization method.
  • the comparison / quantization method is different between dimensions (dimension pairs) having the same region feature quantity calculation method.
  • the region feature quantity calculation method may be different between dimensions (dimension pairs) having the same comparison / quantization method.
  • the comparison / quantization method may be different between dimensions (dimension pairs) having the same region feature amount calculation method. In this case, for example, the region feature quantity calculation method may be different between all dimensions (dimension pairs) having the same comparison / quantization method. Further, the comparison / quantization method may be different between all dimensions (dimension pairs) having the same region feature quantity calculation method.
  • the region feature quantity calculation method may be different between dimensions (dimension pairs) of a predetermined ratio or more.
  • the comparison / quantization method may be different between dimensions (dimension pairs) of a certain ratio or more among dimensions (dimension pairs) having the same region feature quantity calculation method.
  • the region feature quantity calculation method may be different between at least one dimension (dimension pair) among dimensions (dimension pairs) having the same comparison / quantization method.
  • the comparison / quantization method may be different between at least one dimension (dimension pair) among dimensions (dimension pair) having the same region feature quantity calculation method.
  • the comparison method acquisition unit 6 acquires information indicating the comparison / quantization method corresponding to the dimension supplied from the dimension determination unit 1 from the dimension-specific comparison / quantization method information supplied as an input, and sends the information to the comparison unit 4B. Supply.
  • the comparison unit 4B compares, for each dimension, the first region feature amount and the second region feature amount supplied from the region feature amount calculation unit 3 and a comparison / quantization method supplied from the comparison method acquisition unit 6. Are compared and quantized according to the information indicating the output, and a quantization index is output.
  • the comparison unit 4B may be configured to include both the comparison unit 4 of the first embodiment and the comparison unit 4B of the second embodiment as necessary. .
  • the comparison means 4B needs to synchronize the dimension of the supplied region feature quantity with the dimension of the information indicating the comparison / quantization method.
  • the extraction area acquisition means 2 acquires information indicating the first extraction area and the second extraction area of the dimension n from the dimension-specific extraction area information supplied as input, and the area feature amount calculation means 3 (Step D2).
  • the comparison method acquisition unit 6 acquires information indicating the comparison / quantization method corresponding to the dimension n from the dimension-specific comparison / quantization method information supplied as an input, and supplies the information to the comparison unit 4B (step). D3).
  • the area feature quantity calculating means 3 calculates the first area feature quantity of the dimension n and the second area feature quantity from the image supplied as input, and supplies them to the comparison means 4B (step D4). ).
  • the comparison unit 4B compares the first area feature quantity of the dimension n and the second area feature quantity, quantizes the comparison result, and outputs a quantization index (step D5).
  • it is determined whether or not the output of the quantization index has been completed for all dimensions step D6). If the output of the quantization index is completed for all dimensions, the process ends. If the output of the quantization index has not been completed for all dimensions, the process proceeds to step D7.
  • extraction processing is performed in order from dimension 1 to dimension N, but the order is not limited to this and may be arbitrary. In addition to this processing procedure, extraction processing for a plurality of dimensions may be performed in parallel. Furthermore, the order of step D2 and step D3 may be reversed, and step D3 may be executed immediately before step D5.
  • the feature vector to be extracted has 300 dimensions (from the first dimension to the 300th dimension).
  • each dimension extraction area is composed of quadrangular shapes.
  • FIG. 14 shows dimension-specific extraction region information supplied as an input to the extraction region acquisition means 2 in the fifth embodiment.
  • FIG. 14 shows the XY coordinate values of the four corners of the extraction area (first extraction area and second extraction area) for each dimension with respect to an image size of 320 pixels wide ⁇ 240 pixels long, which is a prescribed image size. Indicates.
  • the extraction area of the first dimension is configured by a rectangle having four corners with coordinate values (262.000,163.000), coordinate values (178.068,230.967), coordinate values (184.594,67.411), and coordinate values (100.662,135.378).
  • a first extraction area and a first rectangle composed of four corners with coordinate values (161.000,133.000), coordinate values (156.027,132.477), coordinate values (164.240,102.170), and coordinate values (159.268,101.647) It consists of an extraction area.
  • An extraction area for each dimension is an area surrounded by the coordinate values of the four corners of an image normalized to an image size of 320 pixels wide by 240 pixels wide. It is a set of pixels with integer coordinate values included in. However, negative coordinate values included in the region surrounded by the coordinate values at the four corners are not included in the extraction region.
  • FIG. 15 shows the dimension-specific region feature value calculation method information supplied as an input to the region feature value calculation method acquisition unit 5 in the fifth embodiment.
  • the average value of the luminance values of the pixel groups included in each extraction region is This is a region feature amount.
  • FIG. 17 shows dimension-specific comparison / quantization method information supplied as an input to the comparison method acquisition unit 6 in the fifth embodiment.
  • the comparison / quantization method B or the comparison / quantization method G is used for each dimension, and the parameter values are different for each dimension.
  • the feature vector to be extracted has 300 dimensions (from the first dimension to the 300th dimension).
  • the information shown in FIG. 14 is used as dimension-specific extraction region information supplied as an input to the extraction region acquisition means 2 as in the fifth embodiment.
  • the information shown in FIG. 17 is used as the dimension-specific comparison / quantization method information supplied as an input to the comparison method acquisition unit 6 as in the fifth embodiment.
  • FIG. 16 shows the dimension-specific region feature value calculation method information supplied as an input to the region feature value calculation method acquisition unit 5 in the sixth embodiment.
  • the average value of the luminance values of the pixel groups included in the extraction region (the first extraction region and the second extraction region) or the percentile luminance value feature amount is used for each dimension. Even when the same percentile luminance value feature value is used, the feature value is different for each dimension.
  • the first dimension is an average value of luminance values of pixels included in the extraction region.
  • the fourth dimension is a percentile luminance value feature amount, Y (floor (N ⁇ 20.0 / 100).
  • the eighth dimension is a percentile luminance value feature amount, Y (floor (N ⁇ 80.0 / 100). 100).
  • the feature vector to be extracted has 325 dimensions (from the first dimension to the 325th dimension).
  • each area is configured by a combination of 1024 blocks that are obtained by dividing an image in the vertical direction 32 and the horizontal direction 32.
  • an index starting from 0 is assigned to each block in order from the upper left, and an area is described using this index.
  • the rectangular area is expressed as ab using the index a of the upper left block and the index b of the lower right block.
  • a rectangle composed of four blocks with indexes 0, 1, 32, and 33 is described as 0-33.
  • 2-67 is an area formed by connecting a rectangle defined by 0-33 and a rectangle defined by 2-67, that is, block numbers 0, 1, 2, 3, 32, An area composed of 33, 34, 35, 66, and 67 is shown.
  • FIG. 25 shows a region corresponding to each dimension of the seventh embodiment by this notation.
  • the above-described 325 dimensions are described for each type of region, divided into FIGS. 25-a, 25-b, 25-c, 25-d, 25-e, 25-f, and 25-g. is doing.
  • the region type is a grouping (typing) of regions having similar region patterns determined by a combination of relative positions and shapes between the first and second extraction regions.
  • the shapes of the first and second extraction regions are both rectangles consisting of 4 vertical blocks and 2 horizontal blocks, or rectangles consisting of 2 vertical blocks and 4 horizontal blocks.
  • the relative positional relationship between the first and second extraction regions is viewed, they exist at adjacent positions so that the long sides of the rectangle overlap each other.
  • the shape of the first and second extraction regions is a shape in which two squares each having two vertical and horizontal blocks are arranged on a diagonal of 45 degrees or 135 degrees so as to share one vertex.
  • the two squares constituting the second area are present at positions immediately adjacent to the left and below the upper left square of the first area. .
  • the shape of the first and second extraction regions is a square composed of 10 blocks vertically and horizontally. Further, when looking at the relative positional relationship between the first and second extraction regions, they exist at positions separated by an integral multiple of 10 blocks both vertically and horizontally.
  • the shape of the first and second extraction regions is a square made up of 6 blocks vertically and horizontally. Further, when looking at the relative positional relationship between the first and second extraction regions, they exist at positions separated by an integral multiple of 6 blocks both vertically and horizontally.
  • first and second extraction areas two areas obtained by dividing a square area into a central square and two outside thereof are defined as first and second extraction areas.
  • the shape of the region is such that the second extraction region is a central square, and the first square is a shape obtained by hollowing out the second extraction region from the entire square. Further, when the relative positional relationship of the regions is viewed, the second extraction region exists at the position of the center hole of the first extraction region.
  • the shape of the region is as follows: the first extraction region is 6 blocks long, the rectangle is 10 blocks wide, the second extraction region is 10 blocks long, It is a 6-block rectangle. Further, when the relative positional relationship between the first and second extraction regions is viewed, the center positions are arranged to coincide with each other.
  • a rectangle consisting of 4 vertical blocks and 12 horizontal blocks or a rectangle consisting of 12 vertical blocks and 4 horizontal blocks is divided into 3 equal parts.
  • the shape of the region is such that the first extraction region is a shape in which two squares composed of 4 blocks in length and width are arranged 4 blocks apart vertically or horizontally, and the second extraction region is a square composed of 4 blocks in length and width. is there.
  • a second extraction area exists between the first extraction areas.
  • region types in FIGS. 25-a, 25-b, 25-c, 25-d, 25-e, 25-f, and 25-g are respectively referred to as region type a, region type b, These will be referred to as region type c, region type d, region type e, region type f, and region type g.
  • an average of luminance values is calculated as a region feature amount, and a feature amount in each dimension is calculated.
  • a value extracted by the above-described various extraction methods such as median or maximum value instead of the average of luminance values may be obtained as the region feature amount.
  • the index corresponding to the threshold is 37.
  • indexes corresponding to threshold values can be obtained for other region types. This is shown in FIG. In this way, when the threshold value is determined for each region type, the probability of occurrence of 0, 1, and ⁇ 1 in each dimension can be made uniform as compared with the case where the threshold value is determined as a whole, and the discrimination ability is improved.
  • the quantization may be performed by the other various quantization methods described above.
  • a representative value for example, an average value of luminance values of pixels in a block
  • an area feature amount is extracted therefrom. You may do it.
  • region type has symmetry as a whole. For this reason, even when the right and left of the image are inverted or the image is inverted up and down, by appropriately changing the correspondence and sign of the dimensions, the feature amount extracted from the image that is horizontally or vertically inverted is also used. The image feature amount can be restored. For this reason, it becomes possible to collate with an image that is reversed left and right or up and down.
  • collating means for collating the quantization index vector output in the present invention will be described with reference to a block diagram.
  • FIG. 20 there is shown a block diagram of a collating unit 100 that collates the quantized index vectors output in the present invention.
  • the dimension determining means 101 is connected to the quantized value acquiring means 102 and 103 and outputs the determined dimension information.
  • the quantized value acquiring unit 102 acquires the quantized index value of the dimension input from the dimension determining unit 101 from the first quantized index vector, and outputs it to the scale calculating unit 104 as the first quantized index value.
  • the quantization value acquisition unit 103 acquires the dimension quantization index value input from the dimension determination unit 101 from the second quantization index vector, and outputs it to the scale calculation unit 104 as the second quantization index value.
  • the scale calculation means 104 calculates and outputs a scale representing identity from the first and second quantization index values output from the quantized value acquisition means 102 and 103, respectively.
  • the collation unit 100 receives a first quantization index vector that is a quantization index vector extracted from the first image and a second quantization that is a quantization index vector extracted from the second image. An index vector is input. The input first and second quantization index vectors are input to the quantized value acquisition means 102 and 103, respectively.
  • Dimension information output from the dimension determining means 101 is also input to the quantized value acquisition means 102 and 103.
  • the dimension determining unit 101 sequentially outputs information specifying each dimension of the quantization index vector which is an N-dimensional vector.
  • the order of output does not necessarily increase one by one from 1 to N, and any order may be used as long as the dimensions from 1 to N are specified without excess or deficiency.
  • the quantized value acquisition means 102 and 103 acquire the quantization index value of the dimension specified by the dimension information output from the dimension determination means 101 from the input quantization index vector. Then, the acquired quantization index value is output to the scale calculation means 104.
  • the scale calculation unit 104 compares the first quantization index value output from the quantization value acquisition unit 102 with the second quantization index value. This comparison is performed for each dimension, and a similarity measure (or distance measure) between the first and second quantization index vectors is calculated as an identity measure.
  • the obtained identity scale value is compared with a predetermined threshold value to determine identity.
  • identity scale is a scale representing similarity, it is determined that they are the same when the scale value is equal to or greater than a threshold value.
  • identity scale is a scale representing a distance, it is determined that they are the same when the scale value is equal to or smaller than a threshold value.
  • FIG. 21 is a flowchart showing the operation of the verification unit 100.
  • the dimension (number) of the quantization index vector is represented by n, and there are a total of N dimensions from 1 to N.
  • a variable for calculating the similarity is represented by B.
  • the quantized value acquisition means 102 and 103 the first quantized index value and the second quantized index value of the dimension n are obtained from the first quantized index vector and the second quantized index vector. It is acquired and supplied to the scale calculation means 104 (step S102).
  • the scale calculation means 104 calculates the similarity ⁇ B between the feature amounts corresponding to the respective quantization indexes from the first quantization index value and the second quantization index value (step S104). .
  • the representative value of the feature value before quantization may be calculated from the quantization index, and a value that increases as the difference between the representative values decreases may be used as ⁇ B.
  • a table that can subtract the value of ⁇ B by the combination of the quantization index values is held, and this table is obtained from the combination of the quantization index values.
  • the value of ⁇ B may be directly obtained by using it.
  • step S106 the value of ⁇ B is added to the variable B (step S106).
  • the value of ⁇ B is 0, instead of adding 0 to the variable B, it may be controlled not to add.
  • step S108 it is checked whether or not the dimension number n has reached the number of dimensions N (step S108). If not, the process proceeds to step S112. If the dimension number n has been reached, the value of the variable B at that time is determined as an identity measure. (Scale representing similarity) is output (step S110), and the process is terminated.
  • extraction processing is performed in order from dimensions 1 to N, but the order is not limited to this and may be arbitrary.
  • FIG. 22 is another flowchart showing the operation of the verification unit 100. Also in the flowchart of FIG. 22, it is assumed that the dimension (number) of the quantization index vector is represented by n, and there are a total of N dimensions from 1 to N. A variable for calculating the distance scale is represented by C.
  • steps S100, S104, S106, and S110 are replaced with steps S200, S204, S206, and S210, respectively.
  • variable C is set to zero.
  • the representative value of the feature value before quantization may be calculated from the quantization index, and a value that becomes smaller as the difference between the representative values is smaller may be used as ⁇ C.
  • a table that can subtract the value of ⁇ C by the combination of the quantization index values is held, and this table is obtained from the combination of the quantization index values. The value of ⁇ C may be directly obtained by using it.
  • step S206 the value of ⁇ C is added to the variable C.
  • the value of ⁇ C is 0, instead of adding 0 to the variable C, it may be controlled not to add.
  • step S210 the value of the variable C at that time is output as an identity scale (a scale representing distance), and the process ends.
  • step S108 the process proceeds to step S210.
  • extraction processing is performed in order from dimensions 1 to N, but the order is not limited to this and may be arbitrary.
  • FIG. 23 is another flowchart showing the operation of the verification unit 100. Also in the flowchart of FIG. 23, it is assumed that the dimension (number) of the quantization index vector is represented by n, and there are a total of N dimensions from 1 to N. In addition, a variable for calculating the similarity is represented by B, and a variable for counting dimensions where “both quantization indexes are not 0” is represented by A.
  • a and B are set to 0 (step S300), and the process proceeds to step S102.
  • Step S102 is the same as that in FIG. 21, and after the end, the process proceeds to step S314.
  • step S314 the scale calculation means 104 checks whether both the first quantization index value and the second quantization index value are zero. When both are 0, the process proceeds to step S108, and when either one is not 0, the value of the variable A is increased by 1 (step S316), and the process proceeds to step S104.
  • steps S104, S106, S108, and S112 is the same as that in FIG. If the dimension number n reaches the dimension number N in step S108, the process proceeds to step S310.
  • extraction processing is performed in order from dimensions 1 to N, but the order is not limited to this and may be arbitrary.
  • the image identifier extraction device of the present invention can be realized by a computer and a program, as well as by hardware.
  • the program is provided by being recorded on a computer-readable recording medium such as a magnetic disk or a semiconductor memory, and is read by the computer at the time of starting up the computer, etc.
  • dimensional determination means extraction area acquisition means, area feature quantity calculation means, comparison means, area feature quantity calculation method acquisition means, and comparison method acquisition means.
  • the present invention relates to a patent application of Japanese Patent Application No. 2009-061023 filed on March 13, 2009 in Japan, and Japanese Patent Application No. 2009- filed on April 14, 2009 in Japan.
  • the benefit of the priority claim based on the patent application No. 097862 is enjoyed, and all the contents described in the patent application are included in this specification.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Analysis (AREA)

Abstract

 この画像識別子抽出装置は、抽出手段と生成手段とを備える。抽出手段は、対をなす2つの部分領域の形状の組み合わせと、対をなす2つの部分領域の相対的な位置関係との双方が、他の少なくとも1つの部分領域対と相違する1以上の部分領域対を含む、画像中の複数の部分領域対に従って、画像の各部分領域から領域特徴量を抽出する。生成手段は、上記抽出した各部分領域毎の領域特徴量に基づいて、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる方法を使用して、上記画像の識別に用いる画像識別子を生成する。

Description

画像識別子抽出装置
 本発明は、画像を識別する(同一性を判定する)ための特徴量である画像識別子を抽出する装置に関する。
 画像識別子は、画像を識別する(同一性を判定する)ための画像特徴量である。ある画像から抽出した画像識別子と、別の画像から抽出した画像識別子とを比較し、その比較結果から、2つの画像が同一である度合いを示す同一性尺度(一般的には、類似度または距離という)を算出することができる。また、算出した同一性尺度をある閾値と比較することにより、2つの画像が同一であるか否かを判定することができる。ここで「2つの画像が同一」とは、画像信号(画像を構成する画素の画素値)のレベルで2つの画像が同一である場合だけに限らず、画像の圧縮形式(フォーマット)の変換、画像のサイズ・アスペクト比の変換、画像の色調の調整、画像への各種フィルタ処理(鮮鋭化、平滑化など)、画像への局所的な加工(テロップ重畳、切抜きなど)、画像の再キャプチャリング、などの各種改変処理によって、一方の画像が他方の画像の複製された画像である場合も含む。画像識別子を用いれば、例えば、画像、または画像の集合体である動画像の複製を検知できるため、画像または動画像の違法コピー検知システムなどに応用することができる。
 画像識別子の一例が、特許文献1に記載されている。図18は、特許文献1に記載されている画像識別子の抽出方法を示す図である。この画像識別子は、複数の次元(図18では16次元)の特徴ベクトルである。画像240内のあらかじめ定められた位置の32個の長方形領域244(図18ではそのうち16個の長方形領域が描かれている)からそれぞれ平均輝度値を算出し、対となる長方形領域の間(図18では対となる長方形領域を点線248で結んでいる)で平均輝度値の差を算出し、16次元の差ベクトル250を求める。差ベクトル250に対してベクトル変換により合成ベクトルを生成し、合成ベクトルの各次元を量子化して得られた16次元の量子化インデックスベクトルを画像識別子とする。
特表平8-500471号公報
 複数の次元の特徴ベクトルで構成される画像識別子は、次元間の相関が小さいほど、特徴ベクトルが持つ情報量が大きい(冗長性が小さい)ので、異なる画像を識別できる度合いである識別能力が高くなる。反対に、特徴ベクトルの次元間の相関が大きいと、特徴ベクトルが持つ情報量が小さい(冗長性が大きい)ので、識別能力が低くなる。ここで次元間の相関とは、次元の特徴量の生起の類似性の度合いであり、数学的には、例えば、各次元の特徴量の生起を確率変数とした場合の、確率変数間の相関係数や、相互情報量として算出できる値である。このため、複数の次元の特徴ベクトルで構成される画像識別子は、次元間の相関が小さくなるように設計されていることが望ましい。
 画像信号(画像を構成する画素の画素値)は、画像の局所領域間において相関がある。一般的に、局所領域間の距離が近いほど、相関は大きくなる。特に、例えば、ある特定の画像パターンが繰り返し出現する(特にそれが規則正しい周期で出現する場合に)画像(例えば格子状に配置されたビルの窓の画像など、図19(A)を参照)や、ある特定のテクスチャで構成されている画像(図19(B)を参照)などは、画像の局所領域間の相関が大きくなる。
[第1の問題点]
 特許文献1に記載されているような、画像の複数の局所領域から抽出した特徴量から成る特徴ベクトルで構成されている画像識別子は、画像の局所領域間の相関が大きい画像に対して、各次元において特徴量を抽出する局所領域の形状が同一であるため(特許文献1の例では同一の形状の長方形領域)、抽出される特徴量の次元間の相関が大きくなる。そのため、画像識別子(特徴ベクトル)の識別能力が低くなる、という第1の問題点がある。ここで形状が同一とは、領域の大きさや角度(傾き或いは姿勢)も含めて同一であるということである。
 例えば、ある特定の画像パターンが繰り返し出現する画像(図19(A)参照)や、ある特定のテクスチャで構成されている画像(図19(B)参照)などに対しては、特許文献1で記載されているような画像識別子は、識別能力が低くなる。
[第2の問題点]
 特許文献1に記載されている画像識別子の第2の問題点は、特徴量(特徴ベクトル)を算出するための各次元の領域の形状(大きさ、角度も含めて)が同一の長方形であるため、長方形の辺の長さと同じ、あるいは、その整数分の1の周期を持つ周波数成分を検知できないという、周波数上の盲点が存在するということである。その理由は、この特定の周波数の信号成分について領域内で平均をとると、信号成分の大小によらず0となってしまい、その周波数成分の信号を全く検知できなくなるためである。より具体的には、長方形の辺の長さと同じ周期を持つ周波数をf0とすると,周波数nf0(n=1,2,3,…)の成分が検知できなくなる。このため、直流成分とこの周波数成分に信号が集中している画像に対しては、画素値の平均値は直流成分と同じになってしまい、領域間で値の差がなくなる結果、領域間の平均画素値の差として抽出される特徴量は全て0になってしまい、識別できなくなる(識別能力が著しく低下する)。実際には、周波数nf0(n=1,2,3,…)の成分のみではなく、その近傍の一定の周波数領域に対しては同様に検知困難となるため、上記特定周波数に信号成分が集中していなくても、その周波数帯の信号成分が使えないことにより、識別能力が低下する。この問題を軽減するには、周波数f0の値を大きくし、上記検知困難な周波数帯に陥る信号電力を下げることが考えられる。しかしながら、周波数f0の値を大きくすることは、領域の大きさを小さくすることを意味し、特徴量の頑健性(各種改変処理やノイズに対して特徴量が変化しない度合い)の低下につながる。例えば、領域が小さくなることで、多少の位置ずれに対しても、特徴量の値が大きく変化することになり、特徴量の頑健性が下がる。このように、同一の長方形領域を用いる場合には、識別能力をあげた上で頑健性を確保することが極めて難しい。
[発明の目的]
 本発明の目的は、画像の局所領域間の相関が大きい画像や特定の周波数に信号が集中している画像から抽出される画像識別子は異なる画像を識別できる度合いである識別能力が低下する、という課題を解決する画像識別子抽出装置を提供することである。
 本発明の一形態にかかる画像識別子抽出装置は、対をなす2つの部分領域の形状の組み合わせと、対をなす2つの部分領域の相対的な位置関係との双方が、他の少なくとも1つの部分領域対と相違する1以上の部分領域対を含む、画像中の複数の部分領域対に従って、画像の各部分領域から領域特徴量を抽出する抽出手段と、上記抽出した各部分領域毎の領域特徴量に基づいて、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる方法を使用して、上記画像の識別に用いる画像識別子を生成する生成手段とを備える。
 本発明は上述したように構成されているため、画像識別子の、異なる画像を識別できる度合いである識別能力を高くすることができる。特に、画像の局所領域間の相関が大きい画像に対して、この効果は顕著である。
 また本発明によれば、特定の周波数に信号が集中している画像に対しても、識別能力が低下しないという効果がある。
本発明の第1の実施の形態のブロック図である。 次元別抽出情報が示す次元ごとの抽出領域の対の例を示す図である。 本発明の第1の実施の形態における比較手段の一例を示すブロック図である。 本発明の第1の実施の形態における比較手段の別の例を示すブロック図である。 本発明の第1の実施の形態の処理の流れを示すフローチャートである。 本発明の第2の実施の形態の要部ブロック図である。 本発明の第2の実施の形態の処理の流れを示すフローチャートである。 本発明の第3の実施の形態のブロック図である。 次元ごとの領域特徴量算出方法の例を示す図である。 本発明の第3の実施の形態の処理の流れを示すフローチャートである。 本発明の第4の実施の形態のブロック図である。 次元ごとの比較・量子化方法の例を示す図である。 本発明の第4の実施の形態の処理の流れを示すフローチャートである。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別抽出領域情報を示す図である。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別抽出領域情報を示す図である。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別抽出領域情報を示す図である。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別抽出領域情報を示す図である。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別抽出領域情報を示す図である。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別抽出領域情報を示す図である。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別抽出領域情報を示す図である。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別抽出領域情報を示す図である。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別抽出領域情報を示す図である。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別抽出領域情報を示す図である。 本発明の第5の実施の形態で使用する次元別領域特徴量算出方法情報を示す図である。 本発明の第5の実施の形態で使用する次元別領域特徴量算出方法情報を示す図である。 本発明の第5の実施の形態で使用する次元別領域特徴量算出方法情報を示す図である。 本発明の第5の実施の形態で使用する次元別領域特徴量算出方法情報を示す図である。 本発明の第5の実施の形態で使用する次元別領域特徴量算出方法情報を示す図である。 本発明の第6の実施の形態で使用する次元別領域特徴量算出方法情報を示す図である。 本発明の第6の実施の形態で使用する次元別領域特徴量算出方法情報を示す図である。 本発明の第6の実施の形態で使用する次元別領域特徴量算出方法情報を示す図である。 本発明の第6の実施の形態で使用する次元別領域特徴量算出方法情報を示す図である。 本発明の第6の実施の形態で使用する次元別領域特徴量算出方法情報を示す図である。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別比較・量子化方法情報を示す図である。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別比較・量子化方法情報を示す図である。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別比較・量子化方法情報を示す図である。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別比較・量子化方法情報を示す図である。 本発明の第5の実施の形態および第6の実施の形態で使用する次元別比較・量子化方法情報を示す図である。 特許文献1に記載されている画像識別子の抽出方法を示す図である。 局所領域間の相関が大きくなる画像の例を示す図である。 量子化インデックスベクトルを照合する照合手段のブロック図である。 量子化インデックスベクトルを照合する照合手段の処理例を示すフローチャートである。 量子化インデックスベクトルを照合する照合手段の別の処理例を示すフローチャートである。 量子化インデックスベクトルを照合する照合手段の更に別の処理例を示すフローチャートである。 画像を縦方向32、横方向32に分割してできる1024個のブロックに対して付与するインデックスの一例を示す図である。 本発明の第7の実施の形態の各次元に対応する領域のうち或る1つのタイプに属する領域を示す図である。 本発明の第7の実施の形態の各次元に対応する領域のうち或る1つのタイプに属する領域を示す図である。 本発明の第7の実施の形態の各次元に対応する領域のうち或る1つのタイプに属する領域を示す図である。 本発明の第7の実施の形態の各次元に対応する領域のうち或る1つのタイプに属する領域を示す図である。 本発明の第7の実施の形態の各次元に対応する領域のうち或る1つのタイプに属する領域を示す図である。 本発明の第7の実施の形態の各次元に対応する領域のうち或る1つのタイプに属する領域を示す図である。 本発明の第7の実施の形態の各次元に対応する領域のうち或る1つのタイプに属する領域を示す図である。 各次元の領域タイプと次元数、閾値に対応するインデックスとの関係を示す図である。 領域タイプaの次元の第1、第2の抽出領域の一例を示す図である。 領域タイプbの次元の第1、第2の抽出領域の一例を示す図である。 領域タイプcの次元の第1、第2の抽出領域の一例を示す図である。 領域タイプdの次元の第1、第2の抽出領域の一例を示す図である。 領域タイプeの次元の第1、第2の抽出領域の一例を示す図である。 領域タイプfの次元の第1、第2の抽出領域の一例を示す図である。 領域タイプgの次元の第1、第2の抽出領域の一例を示す図である。
[第1の実施の形態]
[第1の実施の形態の構成]
 次に、本発明の第1の実施の形態について図面を参照して詳細に説明する。
 図1を参照すると、本発明の第1の実施の形態に係る画像識別子抽出装置は、入力された画像に対して、複数の次元から成る特徴ベクトル(より具体的には量子化インデックスベクトル)を画像識別子として出力するシステムであり、次元決定手段1と、抽出領域取得手段2と、領域特徴量算出手段3と、比較手段4と、から構成されている。
 次元決定手段1は、次に抽出する特徴ベクトルの次元を決定し、抽出領域取得手段2へ供給する。次元決定手段1は、順次、抽出する特徴ベクトルの次元を供給し、抽出領域取得手段2以降の構成要素は、供給された次元に対応する特徴量を抽出する。例えば、特徴ベクトルがN次元から構成される場合、次元決定手段1は第1次元から第N次元までを順に抽出領域取得手段2へ供給してもよい。最終的に特徴ベクトルの全ての次元が供給されれば、供給する次元の順番は任意でよい。複数の次元が並列に供給されてもよい。
 抽出領域取得手段2には、次元決定手段1からの次元とは別に、入力として次元別抽出領域情報が供給される。
 次元別抽出領域情報は、あらかじめ規定された、特徴ベクトルの次元ごとに対応付けられた、その次元の特徴量を抽出するための第1の抽出領域と第2の抽出領域の対を示す情報である。第1および第2の抽出領域は必須条件として、以下の特徴を有する。
[第1および第2の抽出領域の必須条件]
 第1および第2の抽出領域の必須条件は、次元間で抽出領域対の相対的な位置が異なることに加えて、次元間で抽出領域対の形状の組み合わせが異なることである。
 上記必須条件を満たす、次元別抽出情報が示す次元ごとの抽出領域の対の例を図2に示す。図18に示した画像識別子の抽出領域とは異なり、次元間の抽出領域の対の形状の組み合わせが異なる。ここで異なる形状とは、角度の異なる合同な形状や(例えば、図2の第1次元の第2の抽出領域と、第7次元の第1の抽出領域)、大きさの異なる相似な形状(例えば、図2の第1次元の第2の抽出領域と、第9次元の第2の抽出領域)も含む。なお、特徴ベクトルの全次元の中に、抽出領域の対の形状の組み合わせの異なる次元のペアが、少なくとも1つ存在することが最低条件である。抽出領域の対の形状(の組み合わせ)が相互に異なる次元が多いほど、望ましい。これは、抽出領域の対の形状(の組み合わせ)が相互に異なる次元が多いほど、特徴ベクトルのより多くの次元間で相関が小さくなり、識別能力が高くなるからである。例えば、特徴ベクトルの全ての次元間で、抽出領域の対の形状(の組み合わせ)が相互に異なっていてもよい。
 ある次元における第1の抽出領域と、第2の抽出領域とは、図2の第9次元のように、同じ形状である必要はなく、図2の他の次元のように、形状が異なっていてもよい。各次元での第1の抽出領域と第2の抽出領域の形状が異なっていると、第1の抽出領域と第2の抽出領域から抽出される特徴量の相関が小さくなり、識別能力が高くなるため、望ましい。また、第1の抽出領域と第2の抽出領域が同時に同じ周波数に関して周波数的な盲点となる可能性が低くなるため、識別能力が高くなる。
 各々の抽出領域の形状は任意である。例えば、図2の第6次元の第2の抽出領域のような、任意の複雑な形状であっても構わない。画像の複数の画素で構成されるものであれば、例えば、図2の第7次元や第10次元のように、線分や曲線であっても構わない。また例えば、第8次元の第1の抽出領域、第11次元の第1と第2の抽出領域、第12次元の第1の抽出領域のように、抽出領域が、連続しない複数の小領域から構成されるものであってもよい。このように、任意の複雑な形状の抽出領域を含むことによって、そこから抽出される特徴量の次元間の相関を小さくすることができ、識別能力を高くすることができる。
 また、例えば、図2の第5次元のように、第1の抽出領域と第2の抽出領域の一部が重複していてもよい。また抽出領域対のいずれか一方が、もう一方の中に内包されていてもよい。このように、抽出領域の対に重複を許容することにより、より多くの抽出領域対のパターン(相対的位置・距離)を取れるため、次元間の相関を小さくすることができるパターンを増やすことができ、識別能力をより高くする可能性が増える。
 また、図18に示した画像識別子の抽出領域とは異なり、図2に示した各次元のように、次元間で抽出領域が一部重複していてもよい。図18に示した画像識別子の抽出領域のように、次元間で抽出領域を排他的に取ると、取れる抽出領域対のパターンが限られてしまう。図2に示したように、次元間での抽出領域に重複を許容することにより、より多くの抽出領域対のパターンを取れるため、次元間の相関を小さくすることができるパターンを増やすことができ、識別能力をより高くする可能性が増える。ただし、次元間での抽出領域の重複が多すぎると、次元間の相関が大きくなってしまい、識別能力が低くなるため、望ましくない。
また、抽出領域が互いに近傍に存在する次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせは互いに異なることが望ましい。例えば、抽出領域が互いに近傍に存在する次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせを互いに異なるようにしてもよい。この場合、抽出領域が互いに近傍に存在する全ての次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせを互いに異なるようにしてもよい。また、抽出領域が互いに近傍に存在する次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせを互いに異なるようにしてもよい。また、抽出領域が互いに近傍に存在する次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせを互いに異なるようにしてもよい。
ここで、抽出領域が互いに近傍に存在する次元ペアは、ある所定の方法によって算出される、2つの次元の抽出領域(第1の抽出領域と第2の抽出領域)の間の距離が、ある所定の閾値よりも小さくなるような次元のペアのことである。2つの次元の抽出領域(第1の抽出領域と第2の抽出領域)の距離は、例えば、それぞれの次元の抽出領域(第1の抽出領域と第2の抽出領域)に含まれる全画素の間の距離の平均値として算出してもよい。また例えば、それぞれの次元の抽出領域(第1の抽出領域と第2の抽出領域)に含まれる全画素の間の距離の最小値として算出してもよい。また例えば、それぞれの次元の抽出領域の重心座標を算出し(第1の抽出領域の重心座標と第2の抽出領域の重心座標)、それらの間の距離を算出し(計4つの距離が計算できる)、それらの平均値や最小値として、2つの次元の抽出領域の距離を定義してもよい。ただし、2つの次元の抽出領域の距離は、これらの算出方法に限られない。
抽出領域が互いに近傍に存在する次元間(次元ペア)で、それらの次元間の抽出領域の対の形状の組み合わせを互いに異なるようにすることで、以下の効果がある。抽出領域が互いに近傍に存在する次元間では、抽出される特徴量の相関が大きくなる傾向になる。それらの次元間の抽出領域の対の形状の組み合わせが同じであると、特徴量の相関はさらに大きくなり、特徴ベクトル全体としての冗長性が高くなり、識別能力が低下する。従って、抽出領域が互いに近傍に存在することにより、抽出される特徴量の相関が大きくなるのを打ち消すように、それらの次元間の抽出領域の対の形状の組み合わせを異なるようにすることにより、特徴量の相関を小さくでき、特徴ベクトル全体としての冗長性を小さくし、識別能力の低下を防ぐことができる。
また逆の観点から、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域の間の距離が大きいほうが望ましい。例えば、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)では、それらの次元間(次元ペア)の間の距離がある所定の閾値よりも大きくなるようにしてもよい。この場合、例えば、抽出領域の対の形状の組み合わせが同一である全ての次元間(次元ペア)において、それらの次元間(次元ペア)の間の距離がある所定の閾値よりも大きくなるようにしてもよい。また、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の間の距離がある所定の閾値よりも大きくなるようにしてもよい。また、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の間の距離がある所定の閾値よりも大きくなるようにしてもよい。
また、抽出領域が隣接する次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせは互いに異なることが望ましい。ここで、抽出領域が隣接する次元間(次元ペア)とは、抽出領域の一部の周囲が互いに接している、ということである。例えば、抽出領域が隣接する次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせを互いに異なるようにしてもよい。この場合、抽出領域が隣接する全ての次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせを互いに異なるようにしてもよい。また、抽出領域が隣接する次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせを互いに異なるようにしてもよい。また、抽出領域が隣接する次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせを互いに異なるようにしてもよい。
抽出領域が互いに隣接する次元間(次元ペア)で、それらの次元間の抽出領域の対の形状の組み合わせを互いに異なるようにすることで、以下の効果がある。抽出領域が互いに隣接する次元間では、抽出される特徴量の相関が大きくなる傾向になる。これは抽出領域が互いに隣接する次元間では、2つの次元の抽出領域の間の距離が近いためである(画像内の近い領域の間は相関が高いため)。それらの次元間の抽出領域の対の形状の組み合わせが同じであると、特徴量の相関はさらに大きくなり、特徴ベクトル全体としての冗長性が高くなり、識別能力が低下する。従って、抽出領域が互いに隣接して存在することにより、抽出される特徴量の相関が大きくなるのを打ち消すように、それらの次元間の抽出領域の対の形状の組み合わせを異なるようにすることにより、特徴量の相関を小さくでき、特徴ベクトル全体としての冗長性を小さくし、識別能力の低下を防ぐことができる。
また逆の観点から、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域が隣接しないようにすることが望ましい。例えば、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域が隣接しないようにしてもよい。この場合、例えば、抽出領域の対の形状の組み合わせが同一である全ての次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が隣接しないようにしてもよい。また、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が隣接しないようにしてもよい。また、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が隣接しないようにしてもよい。
また、抽出領域が一部重複している次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせは互いに異なることが望ましい。例えば、抽出領域が一部重複している次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせを互いに異なるようにしてもよい。この場合、抽出領域が一部重複している全ての次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせを互いに異なるようにしてもよい。また、抽出領域が一部重複している次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせを互いに異なるようにしてもよい。また、抽出領域が一部重複している次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域の対の形状の組み合わせを互いに異なるようにしてもよい。
抽出領域が一部重複している次元間(次元ペア)で、それらの次元間の抽出領域の対の形状の組み合わせを互いに異なるようにすることで、以下の効果がある。抽出領域が一部重複している次元間では、抽出される特徴量の相関が大きくなる傾向になる。これは抽出領域が一部重複する次元間では、特徴量の算出に用いられる領域を一部共有しているためである。それらの次元間の抽出領域の対の形状の組み合わせが同じであると、特徴量の相関はさらに大きくなり、特徴ベクトル全体としての冗長性が高くなり、識別能力が低下する。従って、抽出領域が一部重複することにより、抽出される特徴量の相関が大きくなるのを打ち消すように、それらの次元間の抽出領域の形状の組み合わせを異なるようにすることにより、特徴量の相関を小さくでき、特徴ベクトル全体としての冗長性を小さくし、識別能力の低下を防ぐことができる。
また、抽出領域が一部重複している次元間(次元ペア)で、それらの次元間の抽出領域の対の形状の組み合わせを互いに異なるようにすることで、さらに別の以下の効果がある。抽出領域が一部重複している次元間(次元ペア)で、それらの次元間の抽出領域の対の形状の組み合わせを互いに異なるようにすることで、画像の改竄(特に悪意のある改竄)に対する耐性を強化できる、という効果がある。画像に対する悪意のある改竄は、その画像から抽出される画像識別子(特徴ベクトル)による同一性判定の精度を劣化させる目的で、例えば画像の一部の領域の画素値を変更するなどして行われる。抽出領域が一部重複している次元間では、その画像の重複領域を改竄、すなわちその重複領域の画素値を変更した場合に、そこから抽出される特徴量(領域特徴量)が影響され、改竄前の領域から抽出した特徴量と異なってしまう可能性が高くなる。それらの次元間の抽出領域の対の形状の組み合わせが同じであると、各々の次元の抽出領域から抽出される特徴量(領域特徴量)は同様な影響を受け、同時に特徴量が変更されてしまう可能性が高くなるため、頑健性は低い。それらの次元間の抽出領域の対の形状の組み合わせが異なっていると、重複領域に改竄が行われた場合にも、特徴量が同様な影響を受け、同時に特徴量が変更されてしまう可能性を低くすることができるため、頑健性を確保することができる。すなわち、重複領域を改竄することにより、その重複領域を共有している複数の次元の特徴量を同時に変更してしまうことがより困難となる。従って、抽出領域が一部重複している次元間(次元ペア)で、それらの次元間の抽出領域の対の形状の組み合わせを互いに異なるようにすることで、画像の改竄に対する耐性を強化できる。
また逆の観点から、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域が重複しないようにすることが望ましい。例えば、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域が重複しないようにしてもよい。この場合、例えば、抽出領域の対の形状の組み合わせが同一である全ての次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が重複しないようにしてもよい。また、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が重複しないようにしてもよい。また、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が重複しないようにしてもよい。
 また、全ての次元の抽出領域を統合した場合に、画像内で特徴量が抽出されない領域が小さくなるような(すなわち、画像のほぼ全画面をカバーする)抽出領域の取り方であることが望ましい。これは、図18のように、画像内で特徴量が抽出されない領域が多く含まれていると、画像信号(画像を構成する画素の画素値)に含まれる多くの情報を使用しないことになり、識別能力が高くならないためである。全ての次元の抽出領域を統合した場合に、画像内で特徴量が抽出されない領域が小さくなるような(すなわち、画像のほぼ全画面をカバーする)抽出領域の取り方であることにより、画像信号に含まれるより多くの情報を特徴量に反映できるため、識別能力を高くすることができる。また、全ての次元の抽出領域を統合した場合に、抽出領域に偏りがなく、画像全体から満遍なく取得されていることが望ましい。ただし、ある特定の領域にテロップ重畳などの局所的な加工が施される確率が高い場合は、その領域を避けて抽出領域が取得されていることが望ましい。また、画像の縁などの周辺領域には画像の特徴部分が一般的に存在しないことが多いため、周辺領域を避けて抽出領域が取得されていることが望ましい。
 その他、抽出領域の大きさ、相対的位置(距離、方向)が一定の分布(例えば一様分布)に従うことが望ましい。その理由は、相対的位置(距離、方向)が一様分布に従うことによって、距離や方向に対して偏りがなく、特定の距離や方向に集中することがないため、より多くの多様性がとれるためである。また、相対的位置が近いほど、その領域間の相関が大きくなるため、それを打ち消すために、相対的位置が近いものほどより形状の差が大きいほうが望ましい。
 次元別抽出領域情報は、次元ごとの第1の抽出領域と第2の抽出領域とが一意に特定できる情報であれば、どのような形式の情報であっても構わない。また抽出領域は、如何なるサイズやアスペクト比の画像に対しても、常に同じ領域である必要があるため、次元別抽出領域情報は、如何なるサイズやアスペクト比の画像に対しても、同じ抽出領域を取得できる形式の情報である必要がある。例えば、次元別抽出領域情報は、ある規定のサイズとアスペクト比の画像(例えば、横幅320画素×縦幅240画素の画像)に対して、その抽出領域の位置・形状を記述したものであってもよい。この場合、ある任意のサイズとアスペクト比で入力された画像に対しては、まず画像をその規定のサイズとアスペクト比にリサイズしてから、次元別抽出領域情報に記述されている抽出領域の位置・形状に従って、抽出領域を特定すればよい。あるいは逆に、入力された画像の任意のサイズとアスペクト比の画像に合わせて、次元別抽出領域情報に記述されている抽出領域の位置・形状を変換して、抽出領域を特定してもよい。
 次元別抽出領域情報に含まれる各々の抽出領域を示す情報は、例えば、ある規定のサイズとアスペクト比の画像(例えば、横幅320画素×縦幅240画素の画像)に対して、抽出領域を構成する全ての画素の座標値の集合を記述した情報であってもよい。また次元別抽出領域情報に含まれる各々の抽出領域を示す情報は、例えば、ある規定のサイズとアスペクト比の画像に対して、抽出領域の位置・形状をパラメータ記述した情報であってもよい。例えば抽出領域の形が四角形である場合は、四角形の四隅の座標値を記述した情報であってもよい。また例えば抽出領域の形が円である場合は、円の中心の座標値と半径の値としてもよい。
 また、擬似乱数の種(シード)を次元別抽出領域情報として、抽出領域取得手段2の内部でその種からスタートして擬似乱数を発生させて、乱数に従って異なる形状の抽出領域を生成していく(例えば乱数に従って四角形の四隅を決定していくなど)、という方法も採用することができる。具体的には、例えば以下の手順で、次元別抽出領域を取得することができる。
 (1)擬似乱数の種(シード)が次元別抽出領域情報として供給される。
 (2)次元n=1とする。
 (3)擬似乱数を発生させ、次元nの第1の抽出領域の四角形の四隅を決定する。
 (4)擬似乱数を発生させ、次元nの第2の抽出領域の四角形の四隅を決定する。
 (5)次元n=n+1として、(3)へ戻る。
 乱数に基づいて抽出領域を決定しているので、生成される抽出領域は次元毎に異なる形状になる。また、擬似乱数のシードが同じであれば、毎回(どの画像に対しても)同じ乱数列が発生されるため、異なる画像に対しても同じ抽出領域が再現される。
 抽出領域取得手段2は、入力として供給される次元別抽出領域情報から、次元決定手段1から供給される次元に対応する第1の抽出領域と第2の抽出領域を示す情報を取得し、抽出領域代表値算出手段3へ供給する。
 領域特徴量算出手段3には、抽出領域取得手段2からの入力(第1の抽出領域と第2の抽出領域を示す情報)とは別に、入力として、画像識別子の抽出対象となる画像が供給される。領域特徴量算出手段3は、第1の領域特徴量算出手段31と第2の領域特徴量算出手段32とを有する。領域特徴量算出手段3は、第1の領域特徴量算出手段31を用いて、入力として供給される画像から、次元ごとに、抽出領域取得手段2から供給される第1の抽出領域を示す情報に基づき、第1の抽出領域の特徴量を第1の領域特徴量として算出し、比較手段4へ供給する。また、領域特徴量算出手段3は、第2の領域特徴量算出手段32を用いて、入力として供給される画像から、次元ごとに、抽出領域取得手段2から供給される第2の抽出領域を示す情報に基づき、第2の抽出領域の特徴量を第2の領域特徴量として算出し、比較手段4へ供給する。
 なお、第1の抽出領域と第2の抽出領域を示す情報に基づいて、入力される画像に対するそれぞれの抽出領域を特定するためには、必要に応じて領域特徴量算出手段3は、次元別抽出領域情報の規定のサイズとアスペクト比に画像をリサイズする。
 領域特徴量算出手段3は、それぞれの抽出領域に含まれる画素群の画素値を用いて、それぞれの抽出領域の領域特徴量を算出する。ここで画素値とは、画像の各画素が持つ信号の値であり、スカラー量またはベクトル量である。例えば、画像が輝度画像の場合は、画素値は輝度値(スカラー量)である。また例えば、画像がカラー画像の場合は、画素値は色成分を表すベクトル量である。例えばカラー画像がRGB画像である場合は、画素値はR成分、G成分、B成分の3次元のベクトル量である。また例えばカラー画像がYCbCr画像である場合は、画素値はY成分、Cb成分、Cr成分の3次元のベクトル量である。
 抽出領域の領域特徴量を算出する方法は、その次元の抽出領域(第1の抽出領域と第2の抽出領域)における算出方法が一定である(どの入力画像に対しても同じ算出方法である)限りは、任意の方法でよい。
 また、算出する領域特徴量は、スカラー量でもよいし、ベクトル量であってもよい。例えば、画素値が輝度値などのスカラー量である場合、領域特徴量を、その抽出領域に含まれる画素値の、平均値、メディアン値、最頻値、最大値、最小値、などと算出してもよい(いずれもスカラー量である)。また例えば、抽出領域に含まれる画素値をソートし、分布(ソートされた順列)の上位または下位から規定の割合の位置にある画素値を、領域特徴量として算出してもよい(これもスカラー量である)。より具体的に、規定の割合として、百分率でP%とした場合(例えばP=25%)を例に挙げて説明する。抽出領域に含まれる計N個の画素の画素値(輝度値)を昇順にソートし、昇順にソートされた画素値(輝度値)の集合をY(i)={Y(0)、Y(1)、Y(2)、…、Y(N-1)}と表す。ここで、昇順にソートされた順列の下位からP%の位置にある画素値は、例えば、Y(floor(N×P/100))となり、この値を抽出領域の領域特徴量として算出する。なお、floor()は、小数点以下の切り捨てを行う関数である。ここで、抽出領域に含まれる画素の輝度値に対して、この式(Y(floor(N×P/100)))を適用して算出された領域特徴量を、「パーセンタイル輝度値特徴量」と呼ぶことにする。
 また例えば、画素値が色成分などのベクトル量の場合は、まずそれらを任意の方法でスカラー量に変換してから、上述した方法によって領域特徴量を算出してもよい。例えば、画素値がRGB成分の3次元のベクトル量である場合は、まずそれらをスカラー量である輝度値に変換してから、上述した方法によって領域特徴量を算出してもよい。また画素値がベクトル量の場合は、例えば、その抽出領域に含まれる画素値の平均ベクトルを領域特徴量としてもよい。
 また例えば、抽出領域に対してエッジ検出や、テンプレートマッチングなどの任意の演算(微分演算、フィルタ演算)を行い、その演算結果を領域特徴量としてもよい。例えば、エッジの方向(勾配の方向)を表す2次元のベクトル量であってもよい。また例えば、あるテンプレートとの類似度などを表すスカラー量であってもよい。
 また例えば、抽出領域に含まれる色分布や、エッジの方向分布、エッジの強度分布を表すヒストグラムを、領域特徴量として算出してもよい(いずれもベクトル量である)。
 また例えば、国際標準規格ISO/IEC 15938-3に規定されている各種特徴量、すなわち、Dominant Color、Color Layout、Scalable Color、Color Structure、Edge Histogram、Homogeneous Texture、Texture Browsing、Region Shape、Contour Shape、Shape 3D、Parametric Motion、Motion Activityなどであってもよい。
 比較手段4は、次元ごとに、領域特徴量算出手段3から供給される第1の領域特徴量と、第2の領域特徴量とを比較し、比較した結果を量子化して得られた量子化インデックスを出力する。比較手段4が、次元ごとに、量子化インデックスを出力することで、最終的に、複数の次元の量子化インデックスから成る量子化インデックスベクトルが出力されることになる。
 比較手段4が、第1の領域特徴量と、第2の領域特徴量とを比較して、量子化する方法は、任意である。また、1つの次元当たりの量子化インデックスの数も任意である。
 比較手段4は、例えば、領域特徴量がスカラー量である場合(例えば輝度値の平均値)、その大小を比較して第1の領域特徴量のほうが大きい場合は量子化インデックスを+1、それ以外の場合は量子化インデックスを-1とする、のようにして+1と-1の2値の量子化インデックスに量子化してもよい。ここで、次元nの第1の領域特徴量をVn1、第2の領域特徴量をVn2とすると、次元nの量子化インデックスQnは、次式で算出することができる。
[式1]
Qn=+1  (Vn1>Vn2 の場合)
   -1  (Vn1≦Vn2 の場合)
 ここで、比較手段4が、上述の式1に基づいた比較・量子化を行う場合における、比較手段4のより詳細な構成図を図3に示す。
 図3を参照すると、比較手段4は、大小比較手段41と、量子化手段42と、から構成されている。
 大小比較手段41は、第1の領域特徴量と、第2の領域特徴量とが供給されると、第1の領域特徴量の値と第2の領域特徴量の値との大小を比較し、その比較結果を量子化手段42へ供給する。すなわち、大小比較手段41は、Vn1とVn2の大小を比較し、比較結果が、Vn1>Vn2であるか、Vn1≦Vn2であるか、のいずれであるかを示す情報を、大小比較結果として量子化手段42へ供給する。
 量子化手段42は、大小比較手段41から供給される大小比較結果に基づいて、式1に従って量子化を行い、量子化インデックスを出力する。すなわち量子化手段42は、比較結果がVn1>Vn2であることを示す情報が供給される場合は、量子化インデックスを+1、比較結果がVn1≦Vn2であることを示す情報が供給される場合は、量子化インデックスを-1、として量子化インデックスを出力する。
 なお、この式1に基づいた比較・量子化方法を比較・量子化方法Aと呼ぶことにする。
 また、比較手段4は、例えば、領域特徴量がスカラー量である場合(例えば輝度値の平均値)、差分値の絶対値がある規定の閾値以下の場合は、第1の領域特徴量と第2の領域特徴量との差がないものをみなし、差がないことを示す量子化インデックス0とし、それ以外の場合は、その大小を比較して第1の領域特徴量のほうが大きい場合は量子化インデックスを+1、それ以外の場合は量子化インデックスを-1とする、のようにして+1、0、-1の3値の量子化インデックスに量子化してもよい。ここで、次元nの第1の領域特徴量をVn1、第2の領域特徴量をVn2とし、規定の閾値をthとすると、次元nの量子化インデックスQnは、次式で算出することができる。
[式2]
Qn=+1  (|Vn1-Vn2|>th かつ Vn1>Vn2 の場合)
   0   (|Vn1-Vn2|≦th の場合)
   -1  (|Vn1-Vn2|>th かつ Vn1≦Vn2 の場合)
 ここで、比較手段4が、上述の式2に基づいた比較・量子化を行う場合における、比較手段4のより詳細な構成図を図4に示す。
 図4を参照すると、比較手段4は、差分値算出手段43と、量子化手段44と、から構成されている。量子化手段44には、あらかじめ規定された、量子化の境界を表す情報(量子化境界情報)である閾値が、入力として供給される。
 差分値算出手段43は、第1の領域特徴量と、第2の領域特徴量とが供給されると、第1の領域特徴量の値と第2の領域特徴量の値との差分値を算出し、算出した差分値を量子化手段44へ供給する。すなわち、差分値算出手段43は、Vn1-Vn2を算出し、その値を量子化手段44へ供給する。
 量子化手段44は、差分値算出手段43から供給される差分値と、入力として供給されるあらかじめ規定された量子化の境界を表す情報(量子化境界情報)である閾値とに基づいて、式2に従って量子化を行い、量子化インデックスを出力する。すなわち量子化手段42は、差分値算出手段41から供給されるVn1-Vn2の値と、入力として供給される閾値thとに基づいて、|Vn1-Vn2|>th かつ Vn1-Vn2>0 の場合は量子化インデックスを+1、|Vn1-Vn2|>th かつ Vn1-Vn2≦0 の場合は量子化インデックスを-1、|Vn1-Vn2|≦th の場合は量子化インデックスを0、として量子化インデックスを出力する。
 なお、この式2に基づいた比較・量子化方法を比較・量子化方法Bと呼ぶことにする。
 また、ここでは差分値に基づいて3値に量子化しているが、差分値の大きさに応じて、より多数(のレベルの)の量子化インデックスに量子化してもよい。この場合も、比較手段4は、図4に示した構成をとり、量子化手段44には、あらかじめ規定された、各レベルの量子化の境界を表す情報(量子化境界情報)として複数の閾値が、入力として供給される。なお、この差分値と、入力として供給される複数の閾値とに基づいて、4レベル以上の複数のレベルの量子化インデックスに量子化する比較・量子化方法を比較・量子化方法Cと呼ぶことにする。
 このように、第1の領域特徴量と第2の領域特徴量との差が小さい(規定の閾値以下の)ときに、差がないものとして、差がないことを表す量子化インデックスを導入することで、式1の方法に比べて、領域特徴量の差が小さい抽出領域の対の次元の特徴量(量子化インデックス)をより安定に、すなわち各種改変処理やノイズに対してより頑健に、することができる。そのため、局所領域間の特徴の差が全体的に少ない、全体的に変化の少ない平坦な画像(例えば青空の画像)に対しても安定した、すなわち各種改変処理やノイズに対しても頑健な、画像識別子(量子化インデックスベクトル)を出力することができる。
 また、比較手段4は、例えば、領域特徴量がベクトル量である場合は、ベクトル量をまずそれらを任意の方法でスカラー量に変換してから、上述した方法によって量子化を行ってもよい(この比較・量子化方法を比較・量子化方法Dと呼ぶことにする)。また例えば、第1の抽出領域のベクトルから第2の抽出領域のベクトルとの差分である差分ベクトルを算出し、差分ベクトルをベクトル量子化して量子化インデックスを算出してもよい。この場合は、例えば、あらかじめ規定された量子化インデックスごとの代表ベクトル(重心ベクトルなど)が供給され、それら代表ベクトルと差分ベクトルとの類似度が最も大きく(距離が最も小さく)なる量子化インデックスに分類してもよい(この比較・量子化方法を比較・量子化方法Eと呼ぶことにする)。また、上述の式2によるスカラー量の量子化と同様に、差分ベクトルのノルムがある規定の閾値以下の場合は、第1の領域特徴量と第2の領域特徴量との差がないものをみなし、差がないことを示す量子化インデックス0として、差がないことを表す量子化インデックスを導入してもよい。
 なお、本発明で出力される量子化インデックスベクトルを照合する際(ある画像から抽出した量子化インデックスベクトルと、別の画像から抽出した量子化インデックスベクトルとを比較して、それらの画像が同一であるか否かを判定する際)は、量子化インデックスが一致する次元数(類似度)、あるいは量子化インデックスが非一致である次元数(ハミング距離)を同一性尺度として算出し、算出された同一性尺度をある閾値と比較して、画像の同一性の判定を行うことができる。また、比較手段4において、量子化インデックスが式2に基づいて算出された場合は、以下のように同一性尺度(類似度)を算出することができる。まず、2つの画像の量子化インデックスベクトルを対応する次元どうしで比較して、「共に量子化インデックスが0」ではない次元の数を算出する(この値をAとする)。次に、「共に量子化インデックスが0」ではない次元において、量子化インデックスが一致する次元の数を算出する(この値をBとする)。そして、類似度をB/Aとして算出する。ここでA=0の場合(すなわち、全ての次元が共に量子化インデックスが0となる場合)は、類似度を規定の数値(例えば0.5)とする。
[第1の実施の形態の動作]
 次に、図5のフローチャートを参照して、第1の実施の形態における画像識別子抽出装置の動作を説明する。図5のフローチャートでは、特徴ベクトルの次元(の番号)をnで表し、次元は1からNまでの合計N次元あるものとする。
 まず、次元決定手段1は、抽出する特徴ベクトルの最初の次元として、次元1を決定し(n=1)、抽出領域取得手段2へ供給する(ステップA1)。
 次に、抽出領域取得手段2は、入力として供給される次元別抽出領域情報から、次元nの第1の抽出領域と第2の抽出領域とを示す情報を取得し、領域特徴量算出手段3へ供給する(ステップA2)。
 次に、領域特徴量算出手段3は、入力として供給される画像から、次元nの第1の領域特徴量と、第2の領域特徴量とを算出し、比較手段4へ供給する(ステップA3)。
 次に、比較手段4は、次元nの第1の領域特徴量と第2の領域特徴量とを比較し、比較した結果を量子化して、量子化インデックスを出力する(ステップA4)。
 次に、全ての次元に対して量子化インデックスの出力が終了したか否かを判定(すなわちn<Nが真であるか偽であるかを判定)する(ステップA5)。全ての次元に対して量子化インデックスの出力が終了した場合(すなわちn<Nが偽である場合)は処理を終了する。全ての次元に対して量子化インデックスの出力が終了していない場合(すなわちn<Nが真である場合)は、ステップA6へ移行する。ステップA6では、次元決定手段1が、抽出する特徴ベクトルの次元として、次の次元を決定し(n=n+1)、抽出領域取得手段2へ供給する。そして、再度ステップA2へ移行する。
 なお、ここでは、次元1から次元Nまで順番に抽出処理を行っているが、順番はこれに限らず任意でよい。またこの処理手順に限らず、複数の次元に対する抽出処理を並列に行うようにしてもよい。
[第1の実施の形態の効果]
 次に、本発明の第1の実施の形態の効果について説明する。
 第1の効果は、複数の次元から成る特徴ベクトルで構成される画像識別子の、異なる画像を識別できる度合いである識別能力を高くすることができることである。特に、画像の局所領域間の相関が大きい画像に対して、この効果は顕著である。
 その理由は、次元間で特徴量を抽出する領域の形状が異なる(領域の形状に多様性がある)ことにより、次元間の相関を小さくできるからである。
 第2の効果は、特定の周波数に信号が集中している画像に対しても、識別能力が低下することがないことである。
 その理由は、次元間で特徴量を抽出する領域の形状が異なる(領域の形状に多様性がある)ことにより、ある特定の周波数に信号が集中している画像に対しても、同時に全ての(多くの)抽出領域の対(次元)の間で特徴量の差が無くなり識別能力が低下するようなことが発生しにくくなるからである。
[第2の実施の形態]
[第2の実施の形態の構成]
 次に、本発明の第2の実施の形態について図面を参照して詳細に説明する。
 本発明の第2の実施の形態は、図1に示した第1の実施の形態における比較手段4が、図6に詳細を示す比較手段4Aに置き換わる点において、異なる。比較手段4A以外に関しては、第1の実施の形態と同様であるため、ここでは説明を省略する。
 図6を参照すると、比較手段4Aは、差分値算出手段43と、量子化境界決定手段45と、量子化手段44と、から構成されている。
 差分値算出手段43は、次元ごとに、領域特徴量算出手段3から供給される第1の領域特徴量と、第2の領域特徴量との差分値を算出し、量子化境界決定手段45と、量子化手段44とへ供給する。
 差分値は、領域特徴量がスカラー量の場合(例えば輝度値の平均値)は、例えば、第1の領域特徴量から第2の領域特徴量を(あるいはその逆)減算して得られたスカラー量である。また、領域特徴量がベクトル量の場合は、例えば、それぞれのベクトルを任意の方法でスカラー量に変換してから、スカラー量の差分値を求めてもよい。また、領域特徴量がベクトル量の場合は、第1の領域特徴量と第2の領域特徴量との差分ベクトルを、差分値(ベクトル量)としてもよい。
 量子化境界決定手段45は、差分値算出手段43から供給される特徴ベクトルの全ての次元の差分値が供給されると、全ての次元の差分値の分布に基づいて、量子化の境界を決定し、決定した量子化境界の情報を量子化手段44へ供給する。ここで全ての次元の差分値の分布とは、差分値(あるいは差分ベクトル)に対する生起の頻度(確率)である。
 また量子化の境界を決定するとは、差分値を量子化する際に、漏れなく、かつ排他的に量子化インデックスに割り当てるためのパラメータを決定する、ということである。差分値がスカラー量である場合は、例えば、各量子化インデックス(量子化レベル)に対する値域(すなわち閾値)を決定し、その値域(閾値)を量子化境界の情報として量子化手段43へ供給する。また差分値がベクトル量である場合は、例えばベクトル量子化を行うためのパラメータ、例えば、各量子化インデックスの代表ベクトル(重心ベクトルなど)を決定し、それを量子化境界の情報として量子化手段44へ供給する。
 量子化境界決定手段45は、差分値がスカラー量の場合であって、M値の量子化を行う場合(M=2、3、…など)に、全ての次元の差分値の分布に基づいて、それぞれの量子化インデックスの全次元に対する割合が均等になるように、量子化の値域(閾値)を決定してもよい。
 例えば、前記式1の変形として、定数αを用いて、Vn1+α>Vn2の場合は量子化インデックス+1、Vn1+α≦Vnの場合は量子化インデックス-1とする2値の量子化(M=2)の場合に、量子化インデックスの+1と-1の割合が均等になるように、差分値の分布の中央の点(左右の分布の積分値が等しくなる点)を量子化の閾値αとして決定してもよい。また差分値がベクトル量である場合も同様に、M値の量子化を行う場合に、全ての次元の差分ベクトルの分布に基づいて、それぞれの量子化インデックスの全次元に対する割合が均等になるように、各量子化インデックスに割り当てられるベクトル空間の領域を決定したり、ベクトル量子化を行う際の各量子化インデックスの代表ベクトル(重心ベクトルなど)を決定してもよい。このように、どの画像に対しても、全次元に対する量子化インデックスの割合を均等にすることで(すなわち、量子化インデックスの偏りを無くす)、エントロピーを高くすることができるため、識別能力を高くすることができる。
 なお、量子化境界決定手段45が、量子化インデックスの全次元に対する割合が均等になるように量子化の境界を決定し、それに基づいて量子化手段44が量子化を行う比較・量子化方法を、比較・量子化方法Fと呼ぶことにする。
 また例えば、量子化境界決定手段45は、差分値がスカラー量の場合であって、上述の式2による3値の量子化を行う場合に(量子化インデックスが+1、0、-1)、差分がないことを示す量子化インデックス0に量子化する際の閾値th(この閾値以下の場合に量子化インデックスを0とする)を、全ての次元の差分値の分布に基づいて決定し、決定した閾値thを量子化手段44へ供給してもよい(第1の実施の形態の図4の比較手段4では、この閾値thはあらかじめ規定されているものである)。例えば、全ての次元の差分値の絶対値を算出し、算出した差分値の絶対値をソートして、その上位または下位から、ある規定の割合(なおこの規定の割合は、例えば、入力として供給されるとする)の点を閾値thとしてもよい(この比較・量子化方法を比較・量子化方法Gと呼ぶことにする)。またここで規定の割合ではなく、+1、0、-1の量子化インデックスの割合が均等に近づくように、閾値thを決定してもよい(この比較・量子化方法を比較・量子化方法Hと呼ぶことにする)。比較・量子化方法Hは、式2に従った場合の、比較・量子化方法Fの具体例に相当する。
 比較・量子化方法Gのより具体的な方法を、規定の割合として、百分率でP%とした場合(例えばP=25%)を例に挙げて説明する。全ての次元(次元数=Nとする)の差分値の絶対値を、昇順にソートし、昇順にソートされた差分値の絶対値の集合をD(i)={D(0)、D(1)、D(2)、…、D(N-1)}と表す。ここで、昇順にソートされた順列の下位からP%の位置にある値は、例えば、D(floor(N×P/100))となり、閾値th=D(floor(N×P/100))となる。なお、floor()は、小数点以下の切り捨てを行う関数である。
 本実施の形態における方法は、第1の実施の形態における、比較手段4が図4の構成をとる場合と対比することができる。第1の実施の形態における図4の構成では、あらかじめ規定された閾値thが入力として供給されるのに対して、第2の実施の形態における上述の方法は、量子化境界決定手段45において、全ての次元の差分値の分布に基づいて、画像に対して適応的に閾値thが算出される。このように第1の実施の形態では閾値thが固定化されており、第2の実施の形態では閾値thが画像に適応的に算出される。画像に適応的に閾値thが算出されることで、閾値thが固定化されている場合と比較して、特徴ベクトルの次元の値が、特定の量子化インデックスに偏る(特定の量子化インデックスの出現確率が高い)ことを抑えることができるため(特に起伏の少ない画像に対してなど)、識別能力を高くすることができる。例えば、第1の実施の形態における固定化された閾値thを用いた場合、起伏の少ない画像は、特徴ベクトルの大多数の次元(または全ての次元)が量子化インデックス0になってしまうのに対して、第2の実施の形態における適応的な閾値thを用いると、起伏の少ない画像に対しては閾値thが小さい値に自動的に調整されるため、特徴ベクトルの大多数の次元が量子化インデックス0になるような事態が発生しない。
 量子化手段44は、次元ごとに、差分値算出手段43から供給される次元ごとの差分値と、量子化境界決定手段45から供給される量子化境界の情報とに基づいて、量子化を行い、量子化インデックスを出力する。
 なお、量子化手段44は、量子化境界決定手段45が出力した量子化境界の情報を無視した量子化を行っては意味がなくなるため、量子化境界決定手段45で量子化境界を決定した際に想定していた量子化方法に従う必要がある。
[第2の実施の形態の動作]
 次に、図7のフローチャートを参照して、第2の実施の形態における画像識別子抽出装置の動作を説明する。図7のフローチャートでは、特徴ベクトルの次元(の番号)をnで表し、次元は1からNまでの合計N次元あるものとする。
 まず、次元決定手段1は、抽出する特徴ベクトルの最初の次元として、次元1を決定し(n=1)、抽出領域取得手段2へ供給する(ステップB1)。
 次に、抽出領域取得手段2は、入力として供給される次元別抽出領域情報から、次元nの第1の抽出領域と第2の抽出領域とを示す情報を取得し、抽出領域代表値算出手段3へ供給する(ステップB2)。
 次に、抽出領域代表値算出手段3は、入力として供給される画像から、次元nの第1の領域特徴量と、第2の領域特徴量とを算出し、差分値算出手段43へ供給する(ステップB3)。
 次に、差分値算出手段43は、次元nの第1の領域特徴量と第2の領域特徴量との差分値を算出し、量子化境界決定手段45と、量子化手段44とへ供給する(ステップB4)。
 次に、全ての次元に対する差分値の算出までの処理が終了したか否かを判定(すなわちn<Nが真であるか偽であるかを判定)する(ステップB5)。全ての次元に対する差分値算出までの処理を終了した場合(すなわちn<Nが偽である場合)はステップB7へ移行する。全ての次元に対する処理が終了していない場合(すなわちn<Nが真である場合)は、ステップB6へ移行する。ステップB6では、次元決定手段1が、抽出する特徴ベクトルの次元として、次の次元を決定し(n=n+1)、抽出領域取得手段2へ供給する。そして、再度ステップB2へ移行する。
 なお、ここでは、次元1から次元Nまで順番に抽出処理を行っているが、順番はこれに限らず任意でよい。
 次に、量子化境界決定手段45は、差分値算出手段43から供給される特徴ベクトルの全ての次元の差分値が供給されると、全ての次元の差分値の分布に基づいて、量子化の境界を決定し、決定した量子化境界の情報を量子化手段44へ供給する(ステップB7)。
 次にステップB8では、量子化を行う(量子化インデックスを算出する)特徴ベクトルの最初の次元として、次元1をセット(n=1)する。
 次に、量子化手段44は、次元nの差分値と、量子化境界決定手段45から供給される量子化境界とに基づいて、量子化を行い、量子化インデックスを出力する(ステップB9)。
 次に、全ての次元に対する量子化インデックスの出力が終了したか否かを判定(すなわちn<Nが真であるか偽であるかを判定)する(ステップB10)。全ての次元に対する量子化インデックスの出力を終了した場合(すなわちn<Nが偽である場合)は処理を終了する。全ての次元に対する量子化インデックスの出力が終了していない場合(すなわちn<Nが真である場合)は、ステップB11へ移行する。ステップB11では、量子化を行う(量子化インデックスを算出する)特徴ベクトルの次元として、次の次元をセットする(n=n+1)。そして、再度ステップB9へ移行する。
 なお、ここでは、次元1から次元Nまで順番に量子化処理を行っているが、順番はこれに限らず任意でよい。
[第2の実施の形態の効果]
 第2の実施の形態では、量子化の境界が固定されている第1の実施の形態と比較して、量子化の境界が画像に対して適応的に(動的に)算出される点が異なる。第1の実施の形態のように、量子化の境界が固定化されていると、特定の画像(例えば起伏の少ない平坦な画像など)に対して、特徴ベクトルの次元の値が、特定の量子化インデックスに偏る(特定の量子化インデックスの出現確率が高い)という事態が発生し(エントロピーが低くなる)、これらの画像に対して識別能力が低下するという問題が発生する。一方で第2の実施の形態のように、量子化の境界が画像に対して適応的に(動的に)算出されることにより、どの画像に対しても、特徴ベクトルの次元の値が、特定の量子化インデックスに偏る(特定の量子化インデックスの出現確率が高い)ことを抑えることができるため、識別能力を高くすることができる。
[第3の実施の形態]
[第3の実施の形態の構成]
 次に、本発明の第3の実施の形態について図面を参照して詳細に説明する。
 図8を参照すると、本発明の第3の実施の形態は、図1に示した第1の実施の形態の構成に、領域特徴量算出方法取得手段5が追加され、領域特徴量算出手段3が、第1および第2の領域特徴量算出手段31Aおよび32Aを有する領域特徴量算出手段3Aに置き換わる点で異なる。なお、それ以外の構成に関しては、第1の実施の形態の構成と同様であるため、ここでは説明を省略する。なお、ここでは、第1の実施の形態との組み合わせとして説明しているが、第2の実施の形態との組み合わせであってもよい。
 領域特徴量算出方法取得手段5には、次元決定手段1からの次元と、次元別領域特徴量算出方法情報とが供給される。
 次元別領域特徴量算出方法情報は、あらかじめ規定された、特徴ベクトルの次元ごとに対応付けられた、その次元での領域特徴量の算出方法を示す情報であり、次元間で領域特徴量算出方法が異なることが必須条件である。なおここで、領域特徴量算出方法が異なるとは、同一の手順に対して異なるパラメータ(閾値など)を適用する場合も含む。
 ここで領域特徴量算出方法とは、例えば、第1の実施の形態の領域特徴量算出手段3の説明で記述した各種方法、またそれに伴うパラメータなどである。
 なお次元別領域特徴量算出方法情報が示す次元ごとの領域特徴量算出方法は、特徴ベクトルの全次元の中に、領域特徴量算出方法の異なる次元のペアが、少なくとも1つ存在することが最低条件である。領域特徴量算出方法が相互に異なる次元が多いほど、望ましい。これは、領域特徴量算出方法が相互に異なる次元が多いほど、特徴ベクトルのより多くの次元間で相関が小さくなり、識別能力が高くなるからである。例えば、特徴ベクトルの全ての次元間で、領域特徴量算出方法が相互に異なっていてもよい。
 なお、次元ごとの領域特徴量算出方法を示す情報の形式は、領域特徴量を算出する方法が一意に特定される限りは、任意の形式であってよい。
 図9に、次元ごとの領域特徴量算出方法の例を示す。図9に示すように、次元間で領域特徴量算出方法が異なる。また図9に示した例のように、スカラー量とベクトル量の特徴量が混在していてもよい(第1、3、5、6、8、9、10、12次元はスカラー量、第2、4、7、11次元はベクトル量)。
また、抽出領域が互いに近傍に存在する次元間(次元ペア)では、それらの次元間(次元ペア)の領域特徴量算出方法が異なることが望ましい。例えば、抽出領域が互いに近傍に存在する次元間(次元ペア)では、それらの次元間(次元ペア)の領域特徴量算出方法が異なるようにしてもよい。この場合、抽出領域が互いに近傍に存在する全ての次元間(次元ペア)において、それらの次元間(次元ペア)の領域特徴量算出方法が異なるようにしてもよい。また、抽出領域が互いに近傍に存在する次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の領域特徴量算出方法が異なるようにしてもよい。また、抽出領域が互いに近傍に存在する次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の領域特徴量算出方法が異なるようにしてもよい。
ここで、抽出領域が互いに近傍に存在する次元ペアは、ある所定の方法によって算出される、2つの次元の抽出領域(第1の抽出領域と第2の抽出領域)の間の距離が、ある所定の閾値よりも小さくなるような次元のペアのことである。2つの次元の抽出領域(第1の抽出領域と第2の抽出領域)の距離は、例えば、それぞれの次元の抽出領域(第1の抽出領域と第2の抽出領域)に含まれる全画素の間の距離の平均値として算出してもよい。また例えば、それぞれの次元の抽出領域(第1の抽出領域と第2の抽出領域)に含まれる全画素の間の距離の最小値として算出してもよい。また例えば、それぞれの次元の抽出領域の重心座標を算出し(第1の抽出領域の重心座標と第2の抽出領域の重心座標)、それらの間の距離を算出し(計4つの距離が計算できる)、それらの平均値や最小値として、2つの次元の抽出領域の距離を定義してもよい。ただし、2つの次元の抽出領域の距離は、これらの算出方法に限られない。
抽出領域が互いに近傍に存在する次元間(次元ペア)で、それらの次元間の領域特徴量算出方法が異なるようにすることで、以下の効果がある。抽出領域が互いに近傍に存在する次元間では、抽出される特徴量の相関が大きくなる傾向になる。それらの次元間の領域特徴量算出方法が同じであると、特徴量の相関はさらに大きくなり、特徴ベクトル全体としての冗長性が高くなり、識別能力が低下する。従って、抽出領域が互いに近傍に存在することにより、抽出される特徴量の相関が大きくなるのを打ち消すように、それらの次元間の領域特徴量算出方法が異なるようにすることにより、特徴量の相関を小さくでき、特徴ベクトル全体としての冗長性を小さくし、識別能力の低下を防ぐことができる。
また逆の観点から、領域特徴量算出方法が同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域の間の距離が大きいほうが望ましい。例えば、領域特徴量算出方法が同一である次元間(次元ペア)では、それらの次元間(次元ペア)の間の距離がある所定の閾値よりも大きくなるようにしてもよい。この場合、例えば、領域特徴量算出方法が同一である全ての次元間(次元ペア)において、それらの次元間(次元ペア)の間の距離がある所定の閾値よりも大きくなるようにしてもよい。また、領域特徴量算出方法が同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の間の距離がある所定の閾値よりも大きくなるようにしてもよい。また、領域特徴量算出方法が同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の間の距離がある所定の閾値よりも大きくなるようにしてもよい。
また、抽出領域が隣接する次元間(次元ペア)では、それらの次元間(次元ペア)の領域特徴量算出方法が異なることが望ましい。ここで、抽出領域が隣接する次元間(次元ペア)とは、抽出領域の一部の周囲が互いに接している、ということである。例えば、抽出領域が隣接する次元間(次元ペア)では、それらの次元間(次元ペア)の領域特徴量算出方法が異なるようにしてもよい。この場合、抽出領域が隣接する全ての次元間(次元ペア)において、それらの次元間(次元ペア)の領域特徴量算出方法が異なるようにしてもよい。また、抽出領域が隣接する次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の領域特徴量算出方法が異なるようにしてもよい。また、抽出領域が隣接する次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の領域特徴量算出方法が異なるようにしてもよい。
抽出領域が互いに隣接する次元間(次元ペア)で、それらの次元間の領域特徴量算出方法が異なるようにすることで、以下の効果がある。抽出領域が互いに隣接する次元間では、抽出される特徴量の相関が大きくなる傾向になる。これは抽出領域が互いに隣接する次元間では、2つの次元の抽出領域の間の距離が近いためである(画像内の近い領域の間は相関が高いため)。それらの次元間の領域特徴量算出方法が同じであると、特徴量の相関はさらに大きくなり、特徴ベクトル全体としての冗長性が高くなり、識別能力が低下する。従って、抽出領域が互いに隣接して存在することにより、抽出される特徴量の相関が大きくなるのを打ち消すように、それらの次元間の領域特徴量算出方法が異なるようにすることにより、特徴量の相関を小さくでき、特徴ベクトル全体としての冗長性を小さくし、識別能力の低下を防ぐことができる。
また逆の観点から、領域特徴量算出方法が同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域が隣接しないようにすることが望ましい。例えば、領域特徴量算出方法が同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域が隣接しないようにしてもよい。この場合、例えば、領域特徴量算出方法が同一である全ての次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が隣接しないようにしてもよい。また、領域特徴量算出方法が同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が隣接しないようにしてもよい。また、領域特徴量算出方法が同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が隣接しないようにしてもよい。
また、抽出領域が一部重複している次元間(次元ペア)では、それらの次元間(次元ペア)の領域特徴量算出方法が異なることが望ましい。例えば、抽出領域が一部重複している次元間(次元ペア)では、それらの次元間(次元ペア)の領域特徴量算出方法が異なるようにしてもよい。この場合、抽出領域が一部重複している全ての次元間(次元ペア)において、それらの次元間(次元ペア)の領域特徴量算出方法が異なるようにしてもよい。また、抽出領域が一部重複している次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の領域特徴量算出方法が異なるようにしてもよい。また、抽出領域が一部重複している次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の領域特徴量算出方法が異なるようにしてもよい。
抽出領域が一部重複している次元間(次元ペア)で、それらの次元間の領域特徴量算出方法が異なるようにすることで、以下の効果がある。抽出領域が一部重複している次元間では、抽出される特徴量の相関が大きくなる傾向になる。これは抽出領域が一部重複する次元間では、特徴量の算出に用いられる領域を一部共有しているためである。それらの次元間の領域特徴量算出方法が同じであると、特徴量の相関はさらに大きくなり、特徴ベクトル全体としての冗長性が高くなり、識別能力が低下する。従って、抽出領域が一部重複することにより、抽出される特徴量の相関が大きくなるのを打ち消すように、それらの次元間の領域特徴量算出方法が異なるようにすることにより、特徴量の相関を小さくでき、特徴ベクトル全体としての冗長性を小さくし、識別能力の低下を防ぐことができる。
また、抽出領域が一部重複している次元間(次元ペア)で、それらの次元間の領域特徴量算出方法が異なるようにすることで、さらに別の以下の効果がある。抽出領域が一部重複している次元間(次元ペア)で、それらの次元間の領域特徴量算出方法が異なるようにすることで、画像の改竄(特に悪意のある改竄)に対する耐性を強化できる、という効果がある。画像に対する悪意のある改竄は、その画像から抽出される画像識別子(特徴ベクトル)による同一性判定の精度を劣化させる目的で、例えば画像の一部の領域の画素値を変更するなどして行われる。抽出領域が一部重複している次元間では、その画像の重複領域を改竄、すなわちその重複領域の画素値を変更した場合に、そこから抽出される特徴量(領域特徴量)が影響され、改竄前の領域から抽出した特徴量と異なってしまう可能性が高くなる。それらの次元間の領域特徴量算出方法が同じであると、各々の次元の抽出領域から抽出される特徴量(領域特徴量)は同様な影響を受け、同時に特徴量が変更されてしまう可能性が高くなるため、頑健性は低い。それらの次元間の領域特徴量算出方法が異なっていると、重複領域に改竄が行われた場合にも、特徴量が同様な影響を受け、同時に特徴量が変更されてしまう可能性を低くすることができるため、頑健性を確保することができる。すなわち、重複領域を改竄することにより、その重複領域を共有している複数の次元の特徴量を同時に変更してしまうことがより困難となる。従って、抽出領域が一部重複している次元間(次元ペア)で、それらの次元間の領域特徴量算出方法が異なるようにすることで、画像の改竄に対する耐性を強化できる。
また逆の観点から、領域特徴量算出方法が同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域が重複しないようにすることが望ましい。例えば、領域特徴量算出方法が同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域が重複しないようにしてもよい。この場合、例えば、領域特徴量算出方法が同一である全ての次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が重複しないようにしてもよい。また、領域特徴量算出方法が同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が重複しないようにしてもよい。また、領域特徴量算出方法が同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が重複しないようにしてもよい。
また同様な観点から、領域特徴量抽出方法が同一である次元間(次元ペア)では、抽出領域の対の形状の組み合わせが異なるようにすることが望ましい。また、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)では、領域特徴量抽出方法が異なるようにすることが望ましい。例えば、領域特徴量抽出方法が同一である次元間(次元ペア)では、抽出領域の対の形状の組み合わせが異なるようにしてもよい。また、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)では、領域特徴量抽出方法が異なるようにしてもよい。この場合、例えば、領域特徴量抽出方法が同一である全ての次元間(次元ペア)において、抽出領域の対の形状の組み合わせが異なるようにしてもよい。また、抽出領域の対の形状の組み合わせが同一である全ての次元間(次元ペア)において、領域特徴量抽出方法が異なるようにしてもよい。また例えば、領域特徴量抽出方法が同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、抽出領域の対の形状の組み合わせが異なるようにしてもよい。また、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、領域特徴量抽出方法が異なるようにしてもよい。また例えば、領域特徴量抽出方法が同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、抽出領域の対の形状の組み合わせが異なるようにしてもよい。また、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、領域特徴量抽出方法が異なるようにしてもよい。
 領域特徴量算出方法取得手段5は、入力として供給される次元別領域特徴量算出方法情報から、次元決定手段1から供給される次元に対応する領域特徴量算出方法を示す情報を取得し、領域特徴量算出手段3Aへ供給する。
 領域特徴量算出手段3Aは、入力として供給される画像から、次元ごとに、抽出領域取得手段2から供給される第1の抽出領域と第2の抽出領域とを示す情報に基づき、領域特徴量算出方法取得手段5から供給される領域特徴量算出方法を示す情報に従って、第1の抽出領域の特徴量と、第2の抽出領域の特徴量とを、それぞれ第1の領域特徴量と第2の領域特徴量として算出し、比較手段4へ供給する。
 領域特徴量算出手段3Aでは、供給される抽出領域を示す情報の次元と、領域特徴量算出方法を示す情報の次元との同期が取れている必要がある。
[第3の実施の形態の動作]
 次に、図10のフローチャートを参照して、第3の実施の形態における画像識別子抽出装置の動作を説明する。図10のフローチャートでは、特徴ベクトルの次元(の番号)をnで表し、次元は1からNまでの合計N次元あるものとする。
 まず、次元決定手段1は、抽出する特徴ベクトルの最初の次元として、次元1を決定し(n=1)、抽出領域取得手段2と領域特徴量算出方法取得手段5とへ供給する(ステップC1)。次に、抽出領域取得手段2は、入力として供給される次元別抽出領域情報から、次元nの第1の抽出領域と第2の抽出領域とを示す情報を取得し、領域特徴量算出手段3Aへ供給する(ステップC2)。
 次に、領域特徴量算出方法取得手段5は、入力として供給される次元別領域特徴量算出方法情報から、次元nに対応する領域特徴量算出方法を示す情報を取得し、領域特徴量算出手段3Aへ供給する(ステップC3)。
 次に、領域特徴量算出手段3Aは、入力として供給される画像から、次元nの第1の領域特徴量と、第2の領域特徴量とを算出し、比較手段4へ供給する(ステップC4)。次に、比較手段4は、次元nの第1の領域特徴量と第2の領域特徴量とを比較し、比較した結果を量子化して、量子化インデックスを出力する(ステップC5)。次に、全ての次元に対して量子化インデックスの出力が終了したか否かを判定する(ステップC6)。全ての次元に対して量子化インデックスの出力が終了した場合は処理を終了する。全ての次元に対して量子化インデックスの出力が終了していない場合は、ステップC7へ移行する。ステップC7では、次元決定手段1が、抽出する特徴ベクトルの次元として、次の次元を決定し(n=n+1)、抽出領域取得手段2と領域特徴量算出方法取得手段5とへ供給する。そして、再度ステップC2へ移行する。
 なお、ここでは、次元1から次元Nまで順番に抽出処理を行っているが、順番はこれに限らず任意でよい。またこの処理手順に限らず、複数の次元に対する抽出処理を並列に行うようにしてもよい。さらに、ステップC2とステップC3の順序を逆にしてもよい。
[第3の実施の形態の効果]
 第1の実施の形態の効果に加えて、異なる画像を識別できる度合いである識別能力を更に高くすることができる。
 その理由は、次元間で領域特徴量算出方法が異なる(領域特徴量算出方法に多様性がある)ことにより、次元間の相関をより小さくできるからである。
[第4の実施の形態]
[第4の実施の形態の構成]
 次に、本発明の第4の実施の形態について図面を参照して詳細に説明する。
 図11を参照すると、本発明の第4の実施の形態は、図1に示した第1の実施の形態の構成に、比較方法取得手段6が追加され、比較手段4が比較手段4Bに置き換わる点で異なる。なお、それ以外の構成に関しては、第1の実施の形態の構成と同様であるため、ここでは説明を省略する。なお、ここでは、第1の実施の形態との組み合わせとして説明しているが、第2の実施の形態および第3の実施の形態との組み合わせであってもよい。
 比較方法取得手段6には、次元決定手段1からの次元と、次元別比較方法情報とが供給される。
 次元別比較・量子化方法情報は、あらかじめ規定された、特徴ベクトルの次元ごとに対応付けられた、その次元での領域特徴量を比較して量子化を行う方法を示す情報であり、次元間で比較・量子化方法が異なることが必須条件である。なおここで、比較・量子化方法が異なるとは、同一の手順に対して異なるパラメータ(閾値、量子化インデックス数など)を適用する場合も含む。
 ここで比較・量子化方法とは、例えば第1の実施の形態の比較手段4の説明で記述した各種比較・量子化の方法、またそれに伴うパラメータ(閾値、量子化インデックス数など)や、第2の実施の形態の比較手段4Aの説明で記述した各種比較・量子化の方法、またそれに伴うパラメータ(閾値、量子化インデックス数など)などである。
 なお次元別比較・量子化方法情報が示す次元ごとの比較・量子化方法は、特徴ベクトルの全次元の中に、比較・量子化方法の異なる次元のペアが、少なくとも1つ存在することが最低条件である。比較・量子化方法が相互に異なる次元が多いほど、望ましい。これは、比較・量子化方法が相互に異なる次元が多いほど、特徴ベクトルのより多くの次元間で相関が小さくなり、識別能力が高くなるからである。例えば、特徴ベクトルの全ての次元間で、比較・量子化方法が相互に異なっていてもよい。
 なお、次元ごとの比較・量子化方法を示す情報の形式は、領域特徴量を比較して量子化する方法が一意に特定される限りは、任意の形式であってよい。
 図12に、次元ごとの比較・量子化方法の例を示す。図12に示すように、次元間で比較・量子化方法が異なる。また、第3、5、12次元のように、同じ比較・量子化方法で、異なるパラメータ(閾値th)を設定してもよい。なお、図12に示した、次元ごとの比較・量子化方法の例は、図9に示した、次元ごとの領域特徴量算出方法の例と対応させており、スカラー量の領域特徴量に対してはスカラー量の比較・量子化方法を、ベクトル量の領域特徴量に対してはベクトル量の比較・量子化方法を例として示した。
また、抽出領域が互いに近傍に存在する次元間(次元ペア)では、それらの次元間(次元ペア)の比較・量子化方法が異なることが望ましい。例えば、抽出領域が互いに近傍に存在する次元間(次元ペア)では、それらの次元間(次元ペア)の比較・量子化方法が異なるようにしてもよい。この場合、抽出領域が互いに近傍に存在する全ての次元間(次元ペア)において、それらの次元間(次元ペア)の比較・量子化方法が異なるようにしてもよい。また、抽出領域が互いに近傍に存在する次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の比較・量子化方法が異なるようにしてもよい。また、抽出領域が互いに近傍に存在する次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の比較・量子化方法が異なるようにしてもよい。
ここで、抽出領域が互いに近傍に存在する次元ペアは、ある所定の方法によって算出される、2つの次元の抽出領域(第1の抽出領域と第2の抽出領域)の間の距離が、ある所定の閾値よりも小さくなるような次元のペアのことである。2つの次元の抽出領域(第1の抽出領域と第2の抽出領域)の距離は、例えば、それぞれの次元の抽出領域(第1の抽出領域と第2の抽出領域)に含まれる全画素の間の距離の平均値として算出してもよい。また例えば、それぞれの次元の抽出領域(第1の抽出領域と第2の抽出領域)に含まれる全画素の間の距離の最小値として算出してもよい。また例えば、それぞれの次元の抽出領域の重心座標を算出し(第1の抽出領域の重心座標と第2の抽出領域の重心座標)、それらの間の距離を算出し(計4つの距離が計算できる)、それらの平均値や最小値として、2つの次元の抽出領域の距離を定義してもよい。ただし、2つの次元の抽出領域の距離は、これらの算出方法に限られない。
抽出領域が互いに近傍に存在する次元間(次元ペア)で、それらの次元間の比較・量子化方法が異なるようにすることで、以下の効果がある。抽出領域が互いに近傍に存在する次元間では、抽出される特徴量の相関が大きくなる傾向になる。それらの次元間の比較・量子化方法が同じであると、特徴量の相関はさらに大きくなり、特徴ベクトル全体としての冗長性が高くなり、識別能力が低下する。従って、抽出領域が互いに近傍に存在することにより、抽出される特徴量の相関が大きくなるのを打ち消すように、それらの次元間の比較・量子化方法が異なるようにすることにより、特徴量の相関を小さくでき、特徴ベクトル全体としての冗長性を小さくし、識別能力の低下を防ぐことができる。
また逆の観点から、比較・量子化方法が同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域の間の距離が大きいほうが望ましい。例えば、比較・量子化方法が同一である次元間(次元ペア)では、それらの次元間(次元ペア)の間の距離がある所定の閾値よりも大きくなるようにしてもよい。この場合、例えば、比較・量子化方法が同一である全ての次元間(次元ペア)において、それらの次元間(次元ペア)の間の距離がある所定の閾値よりも大きくなるようにしてもよい。また、比較・量子化方法が同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の間の距離がある所定の閾値よりも大きくなるようにしてもよい。また、比較・量子化方法が同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の間の距離がある所定の閾値よりも大きくなるようにしてもよい。
また、抽出領域が隣接する次元間(次元ペア)では、それらの次元間(次元ペア)の比較・量子化方法が異なることが望ましい。ここで、抽出領域が隣接する次元間(次元ペア)とは、抽出領域の一部の周囲が互いに接している、ということである。例えば、抽出領域が隣接する次元間(次元ペア)では、それらの次元間(次元ペア)の比較・量子化方法が異なるようにしてもよい。この場合、抽出領域が隣接する全ての次元間(次元ペア)において、それらの次元間(次元ペア)の比較・量子化方法が異なるようにしてもよい。また、抽出領域が隣接する次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の比較・量子化方法が異なるようにしてもよい。また、抽出領域が隣接する次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の比較・量子化方法が異なるようにしてもよい。
抽出領域が互いに隣接する次元間(次元ペア)で、それらの次元間の比較・量子化方法が異なるようにすることで、以下の効果がある。抽出領域が互いに隣接する次元間では、抽出される特徴量の相関が大きくなる傾向になる。これは抽出領域が互いに隣接する次元間では、2つの次元の抽出領域の間の距離が近いためである(画像内の近い領域の間は相関が高いため)。それらの次元間の比較・量子化方法が同じであると、特徴量の相関はさらに大きくなり、特徴ベクトル全体としての冗長性が高くなり、識別能力が低下する。従って、抽出領域が互いに隣接して存在することにより、抽出される特徴量の相関が大きくなるのを打ち消すように、それらの次元間の比較・量子化方法が異なるようにすることにより、特徴量の相関を小さくでき、特徴ベクトル全体としての冗長性を小さくし、識別能力の低下を防ぐことができる。
また逆の観点から、比較・量子化方法が同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域が隣接しないようにすることが望ましい。例えば、比較・量子化方法が同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域が隣接しないようにしてもよい。この場合、例えば、比較・量子化方法が同一である全ての次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が隣接しないようにしてもよい。また、比較・量子化方法が同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が隣接しないようにしてもよい。また、比較・量子化方法が同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が隣接しないようにしてもよい。
また、抽出領域が一部重複している次元間(次元ペア)では、それらの次元間(次元ペア)の比較・量子化方法が異なることが望ましい。例えば、抽出領域が一部重複している次元間(次元ペア)では、それらの次元間(次元ペア)の比較・量子化方法が異なるようにしてもよい。この場合、抽出領域が一部重複している全ての次元間(次元ペア)において、それらの次元間(次元ペア)の比較・量子化方法が異なるようにしてもよい。また、抽出領域が一部重複している次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の比較・量子化方法が異なるようにしてもよい。また、抽出領域が一部重複している次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の比較・量子化方法が異なるようにしてもよい。
抽出領域が一部重複している次元間(次元ペア)で、それらの次元間の比較・量子化方法が異なるようにすることで、以下の効果がある。抽出領域が一部重複している次元間では、抽出される特徴量の相関が大きくなる傾向になる。これは抽出領域が一部重複する次元間では、特徴量の算出に用いられる領域を一部共有しているためである。それらの次元間の比較・量子化方法が同じであると、特徴量の相関はさらに大きくなり、特徴ベクトル全体としての冗長性が高くなり、識別能力が低下する。従って、抽出領域が一部重複することにより、抽出される特徴量の相関が大きくなるのを打ち消すように、それらの次元間の比較・量子化方法が異なるようにすることにより、特徴量の相関を小さくでき、特徴ベクトル全体としての冗長性を小さくし、識別能力の低下を防ぐことができる。
また、抽出領域が一部重複している次元間(次元ペア)で、それらの次元間の比較・量子化方法が異なるようにすることで、さらに別の以下の効果がある。抽出領域が一部重複している次元間(次元ペア)で、それらの次元間の比較・量子化方法が異なるようにすることで、画像の改竄(特に悪意のある改竄)に対する耐性を強化できる、という効果がある。画像に対する悪意のある改竄は、その画像から抽出される画像識別子(特徴ベクトル)による同一性判定の精度を劣化させる目的で、例えば画像の一部の領域の画素値を変更するなどして行われる。抽出領域が一部重複している次元間では、その画像の重複領域を改竄、すなわちその重複領域の画素値を変更した場合に、そこから抽出される特徴量(領域特徴量)が影響され、改竄前の領域から抽出した特徴量と異なってしまう可能性が高くなる。それらの次元間の比較・量子化方法が同じであると、各々の次元の抽出領域から抽出される特徴量(領域特徴量)は同様な影響を受け、同時に特徴量が変更されてしまう可能性が高くなるため、頑健性は低い。それらの次元間の比較・量子化方法が異なっていると、重複領域に改竄が行われた場合にも、特徴量が同様な影響を受け、同時に特徴量が変更されてしまう可能性を低くすることができるため、頑健性を確保することができる。すなわち、重複領域を改竄することにより、その重複領域を共有している複数の次元の特徴量を同時に変更してしまうことがより困難となる。従って、抽出領域が一部重複している次元間(次元ペア)で、それらの次元間の比較・量子化方法が異なるようにすることで、画像の改竄に対する耐性を強化できる。
また逆の観点から、比較・量子化方法が同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域が重複しないようにすることが望ましい。例えば、比較・量子化方法が同一である次元間(次元ペア)では、それらの次元間(次元ペア)の抽出領域が重複しないようにしてもよい。この場合、例えば、比較・量子化方法が同一である全ての次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が重複しないようにしてもよい。また、比較・量子化方法が同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が重複しないようにしてもよい。また、比較・量子化方法が同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、それらの次元間(次元ペア)の抽出領域が重複しないようにしてもよい。
また同様な観点から、比較・量子化方法が同一である次元間(次元ペア)では、抽出領域の対の形状の組み合わせが異なるようにすることが望ましい。また、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)では、比較・量子化方法が異なるようにすることが望ましい。例えば、比較・量子化方法が同一である次元間(次元ペア)では、抽出領域の対の形状の組み合わせが異なるようにしてもよい。また、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)では、比較・量子化方法が異なるようにしてもよい。この場合、例えば、比較・量子化方法が同一である全ての次元間(次元ペア)において、抽出領域の対の形状の組み合わせが異なるようにしてもよい。また、抽出領域の対の形状の組み合わせが同一である全ての次元間(次元ペア)において、比較・量子化方法が異なるようにしてもよい。また例えば、比較・量子化方法が同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、抽出領域の対の形状の組み合わせが異なるようにしてもよい。また、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、比較・量子化方法が異なるようにしてもよい。また例えば、比較・量子化方法が同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、抽出領域の対の形状の組み合わせが異なるようにしてもよい。また、抽出領域の対の形状の組み合わせが同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、比較・量子化方法が異なるようにしてもよい。
また同様な観点から、比較・量子化方法が同一である次元間(次元ペア)では、領域特徴量算出方法が異なるようにすることが望ましい。また、領域特徴量算出方法が同一である次元間(次元ペア)では、比較・量子化方法が異なるようにすることが望ましい。例えば、比較・量子化方法が同一である次元間(次元ペア)では、領域特徴量算出方法が異なるようにしてもよい。また、領域特徴量算出方法が同一である次元間(次元ペア)では、比較・量子化方法が異なるようにしてもよい。この場合、例えば、比較・量子化方法が同一である全ての次元間(次元ペア)において、領域特徴量算出方法が異なるようにしてもよい。また、領域特徴量算出方法が同一である全ての次元間(次元ペア)において、比較・量子化方法が異なるようにしてもよい。また例えば、比較・量子化方法が同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、領域特徴量算出方法が異なるようにしてもよい。また、領域特徴量算出方法が同一である次元間(次元ペア)のうち、ある所定の割合以上の次元間(次元ペア)において、比較・量子化方法が異なるようにしてもよい。また例えば、比較・量子化方法が同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、領域特徴量算出方法が異なるようにしてもよい。また、領域特徴量算出方法が同一である次元間(次元ペア)のうち、少なくとも1以上の次元間(次元ペア)において、比較・量子化方法が異なるようにしてもよい。
 比較方法取得手段6は、入力として供給される次元別比較・量子化方法情報から、次元決定手段1から供給される次元に対応する比較・量子化方法を示す情報を取得し、比較手段4Bへ供給する。
 比較手段4Bは、次元ごとに、領域特徴量算出手段3から供給される第1の領域特徴量と、第2の領域特徴量とを、比較方法取得手段6から供給される比較・量子化方法を示す情報に従って、比較・量子化して、量子化インデックスを出力する。比較手段4Bは、比較・量子化方法によって、必要に応じて、第1の実施の形態の比較手段4と、第2の実施の形態の比較手段4Bの両方を内包した構成となる場合もある。
 比較手段4Bでは、供給される領域特徴量の次元と、比較・量子化方法を示す情報の次元の同期が取れている必要がある。
[第4の実施の形態の動作]
 次に、図13のフローチャートを参照して、第4の実施の形態における画像識別子抽出装置の動作を説明する。図13のフローチャートでは、特徴ベクトルの次元(の番号)をnで表し、次元は1からNまでの合計N次元あるものとする。
 まず、次元決定手段1は、抽出する特徴ベクトルの最初の次元として、次元1を決定し(n=1)、抽出領域取得手段2と比較方法取得手段6とへ供給する(ステップD1)。次に、抽出領域取得手段2は、入力として供給される次元別抽出領域情報から、次元nの第1の抽出領域と第2の抽出領域とを示す情報を取得し、領域特徴量算出手段3へ供給する(ステップD2)。
 次に、比較方法取得手段6は、入力として供給される次元別比較・量子化方法情報から、次元nに対応する比較・量子化方法を示す情報を取得し、比較手段4Bへ供給する(ステップD3)。
 次に、領域特徴量算出手段3は、入力として供給される画像から、次元nの第1の領域特徴量と、第2の領域特徴量とを算出し、比較手段4Bへ供給する(ステップD4)。次に、比較手段4Bは、次元nの第1の領域特徴量と第2の領域特徴量とを比較し、比較した結果を量子化して、量子化インデックスを出力する(ステップD5)。次に、全ての次元に対して量子化インデックスの出力が終了したか否かを判定する(ステップD6)。全ての次元に対して量子化インデックスの出力が終了した場合は処理を終了する。全ての次元に対して量子化インデックスの出力が終了していない場合は、ステップD7へ移行する。ステップD7では、次元決定手段1が、抽出する特徴ベクトルの次元として、次の次元を決定し(n=n+1)、抽出領域取得手段2と比較方法取得手段6とへ供給する。そして、再度ステップD2へ移行する。
 なお、ここでは、次元1から次元Nまで順番に抽出処理を行っているが、順番はこれに限らず任意でよい。またこの処理手順に限らず、複数の次元に対する抽出処理を並列に行うようにしてもよい。さらに、ステップD2とステップD3との順序を逆にしてもよく、ステップD3をステップD5の直前に実行するようにしてもよい。
[第4の実施の形態の効果]
 第1の実施の形態の効果に加えて、異なる画像を識別できる度合いである識別能力を更に高くすることができる。
 その理由は、次元間で比較・量子化方法が異なる(比較・量子化方法に多様性がある)ことにより、次元間の相関をより小さくできるからである。
 次に、本発明における第5~第7の実施の形態について説明する。
[第5の実施の形態]
 第5の実施の形態では、抽出する特徴ベクトルの次元数は300次元(第1次元から第300次元)である。
 第5の実施の形態では、次元ごとの抽出領域(第1の抽出領域と第2の抽出領域)は、様々な形状の四角形から構成される。第5の実施の形態において、抽出領域取得手段2に入力として供給される次元別抽出領域情報を図14に示す。図14は、規定の画像サイズである、横幅320画素×縦幅240画素の画像サイズに対する、次元ごとの抽出領域(第1の抽出領域と第2の抽出領域)の四角形の四隅のXY座標値を示す。例えば、第1次元の抽出領域は、座標値(262.000,163.000)、座標値(178.068,230.967)、座標値(184.594,67.411)、座標値(100.662,135.378)を四隅とする四角形で構成される第1の抽出領域と、座標値(161.000,133.000)、座標値(156.027,132.477)、座標値(164.240,102.170)、座標値(159.268,101.647)を四隅とする四角形で構成される第1の抽出領域とで構成される。
 次元ごとの抽出領域(第1の抽出領域と第2の抽出領域)は、横幅320画素×縦幅240画素の画像サイズに正規化された画像に対して、この四隅の座標値で囲まれる領域の中に含まれる整数値の座標値の画素の集合となる。ただし、四隅の座標値で囲まれる領域の中に含まれる負の座標値は、抽出領域に含まない。
 第5の実施の形態において、領域特徴量算出方法取得手段5に入力として供給される次元別領域特徴量算出方法情報を図15に示す。第5の実施の形態では、全ての次元に対して、それぞれの抽出領域(第1の抽出領域と第2の抽出領域)に含まれる画素群の輝度値の平均値が、それぞれの抽出領域の領域特徴量となる。
 第5の実施の形態において、比較方法取得手段6に入力として供給される次元別比較・量子化方法情報を図17に示す。第5の実施の形態では、次元ごとに、比較・量子化方法Bまたは比較・量子化方法Gが用いられ、次元ごとにそのパラメータの値も異なる。例えば、第1次元は、比較・量子化方法Gで、閾値th=D(floor(300×5.0/100))である。また、例えば第2次元は、比較・量子化方法Gで、閾値th=D(floor(300×10.0/100))である。また、例えば第9次元は、比較・量子化方法Bで、閾値th=3.0である。
[第6の実施の形態]
 第6の実施の形態は、第5の実施の形態と同じく、抽出する特徴ベクトルの次元数は300次元(第1次元から第300次元)である。また第6の実施の形態では、抽出領域取得手段2に入力として供給される次元別抽出領域情報として、第5の実施の形態と同じく図14に示す情報を使用する。さらに第6の実施の形態では、比較方法取得手段6に入力として供給される次元別比較・量子化方法情報として、第5の実施の形態と同じく図17に示す情報を使用する。
 第6の実施の形態において、領域特徴量算出方法取得手段5に入力として供給される次元別領域特徴量算出方法情報を図16に示す。第6の実施の形態では、次元ごとに、抽出領域(第1の抽出領域と第2の抽出領域)に含まれる画素群の輝度値の平均値、または、パーセンタイル輝度値特徴量が用いられ、同じパーセンタイル輝度値特徴量を用いる場合でも、次元ごとにその特徴量は異なる。例えば、第1次元は、抽出領域に含まれる画素の輝度値の平均値である。また、例えば第4次元は、パーセンタイル輝度値特徴量で、Y(floor(N×20.0/100)である。また、第8次元は、パーセンタイル輝度値特徴量で、Y(floor(N×80.0/100)である。
[第7の実施の形態]
 第7の実施の形態は、抽出する特徴ベクトルの次元数は325次元(第1次元から第325次元)である。第7の実施の形態の場合は、各領域は、画像を縦方向32、横方向32に分割してできる1024個のブロックの組み合わせによって構成されている。ここで、各ブロックに対して、図24に示すように、左上から順に0から始まるインデックスを付与し、このインデックスを用いて領域を記述する。具体的には、長方形領域を、その左上のブロックのインデックスaと右下のブロックのインデックスbを用いてa-bのように表現する。例えば、インデックス0、1、32、33の4つのブロックからなる長方形は、0-33のように記述する。また、このようにしてできる長方形を記号“|”によって繋げた場合は、その記号の前後の長方形を連結してできる領域を表現するものとする。例えば、0-33|2-67は、0-33で定義される長方形と、2-67で定義される長方形を連結してできる領域、すなわち、ブロック番号0、1、2、3、32、33、34、35、66、67によって構成される領域を表している。
 この表記によって第7の実施の形態の各次元に対応する領域を示したものが図25である。図では、領域のタイプ別に図25-a、図25-b、図25-c、図25-d、図25-e、図25-f、図25-gに分けて上述の325次元を記述している。ここで、領域のタイプとは、第1、第2の抽出領域間の相対位置や形状の組み合わせによって定まる領域パターンが似たもの同士でグループ化(類型化)したものである。
 具体的には、図25-aの場合は、図27-aに一例を示すように、縦横4ブロックからなる正方形を縦方向か横方向に2等分してできる2つの領域を第1、第2の抽出領域とした場合に相当する。このため、第1、第2の抽出領域の形状は、ともに縦4ブロック、横2ブロックからなる長方形、あるいは縦2ブロック、横4ブロックからなる長方形である。また、第1、第2の抽出領域の相対的な位置関係を見ると、長方形の長い辺同士が重なるように隣接する位置に存在する。
 図25-bの場合は、図27-bに一例を示すように、縦横8ブロックからなる正方形を縦横2等分してできる4つの正方形のうち、左上と右下、右上と左下をそれぞれ組み合わせてできる2つの領域を第1、第2の抽出領域とした場合に相当する。このため、第1、第2の抽出領域の形状は、ともに縦横2ブロックからなる正方形を1つの頂点を共有するように45度あるいは135度の対角線上に2つ配置した形状となっている。また、領域の相対的な位置関係を見ると、第2の領域を構成する2つの正方形が、第1の領域の左上の正方形のすぐ左と下に隣接する位置に第2の領域が存在する。
 図25-cの場合は、図27-cに一例を示すように、第1、第2の抽出領域の形状は、ともに縦横10ブロックからなる正方形である。また、第1、第2の抽出領域の相対的な位置関係を見ると、縦横ともに10ブロックの整数倍だけ離れた位置に存在する。
 図25-dの場合は、図27-dに一例を示すように、第1、第2の抽出領域の形状は、ともに縦横6ブロックからなる正方形である。また第1、第2の抽出領域の相対的な位置関係を見ると、縦横ともに6ブロックの整数倍だけ離れた位置に存在する。
 図25-eの場合は、図27-eに一例を示すように、正方形領域を中心部分の正方形とその外側の2つに分けてできる2つの領域を第1、第2の抽出領域とした場合に相当する。このため、領域の形状は、第2の抽出領域が中心部分の正方形、第1の正方形は全体の正方形から第2の抽出領域をくりぬいた形状である。また、領域の相対的な位置関係を見ると、第1の抽出領域の中央の穴の位置に第2の抽出領域が存在する。
 図25-fの場合は、図27-fに一例を示すように、領域の形状は、第1の抽出領域は縦6ブロック、横10ブロックの長方形、第2の抽出領域は縦10ブロック、横6ブロックの長方形である。また、第1、第2の抽出領域の相対的な位置関係を見ると、中心位置が一致するように配置されている。
 図25-gの場合には、図27-gに一例を示すように、縦4ブロック、横12ブロックからなる長方形、あるいは縦12ブロック、横4ブロックからなる長方形を、長い辺を3等分してできる中央の正方形とそれ以外の2領域を第1、第2の抽出領域とした場合に相当する。このため、領域の形状は、第1の抽出領域は縦横4ブロックからなる正方形を2つ、縦か横に4ブロック離れて配置した形状で、第2の抽出領域は縦横4ブロックからなる正方形である。また、領域の相対的な位置関係を見ると、第1の抽出領域の間に第2の抽出領域が存在する。
 以後、図25-a、図25-b、図25-c、図25-d、図25-e、図25-f、図25-gの領域タイプを、それぞれ領域タイプa、領域タイプb、領域タイプc、領域タイプd、領域タイプe、領域タイプf、領域タイプgと呼ぶことにする。
 第7の実施の形態では、図25で示した各領域において、領域特徴量として輝度値の平均を算出し、各次元の特徴量を算出する。もちろん、輝度値の平均のかわりにメディアンや最大値など、前述の様々な抽出方法によって抽出した値を領域特徴量として求めるようにしてもよい。
 各次元の特徴量の量子化では、上述の領域のタイプ別に閾値を定め、量子化を行うようにする。例えば、式2に従って特徴量を3値に量子化する場合には、領域のタイプ別に、0、1、-1の生起の割合が均等になるように量子化の閾値thを決定し、量子化を行うようにする。具体的には、段落0096で記述した方法をP=33.333%、Nを領域タイプ別の次元数として領域タイプ別に適用し、閾値thを求める。例えば、領域タイプaの場合にはN=113となるため、th=D(floor(113×33.333/100))=D(37)により閾値を算出する。ここで、D(i)(i=0、1、…、N-1)は、領域タイプaに該当する第1次元から第113次元の差分値の絶対値を昇順にソートした集合になる。この場合は閾値に対応するインデックスが37となる。同様に、他の領域タイプに対しても、閾値に対応するインデックスを求めることができる。これを示したのが図26である。このように領域タイプ別に閾値を求める方が、全体で閾値を決める場合に比べて各次元での0、1、-1の発生確率を均一化できるようになり、識別能力が向上する。もちろん、前述の他の様々な量子化方法によって量子化するようにしてもよい。
 なお、第7の実施の形態の場合には、図24で示したブロックごとに代表値(例えば、ブロック内の画素の輝度値の平均値)を先に算出し、それから領域特徴量を抽出するようにしてもよい。これにより、領域内の全画素から直接領域特徴量を抽出する場合よりも高速に抽出できるようになる。また、各領域タイプの抽出領域は、全体として対称性を有する。このため、画像の右と左を反転させたり、上下を反転させたりした場合でも、次元の対応関係と符号を適切に変更することによって、左右または上下反転した画像から抽出された特徴量からもとの画像の特徴量を復元できる。このため、左右あるいは上下を反転させた画像とも照合することができるようになる。
[照合手段の実施の形態]
 次に、本発明で出力される量子化インデックスベクトルを照合する照合手段についてブロック図を用いて説明する。
 図20を参照すると、本発明で出力される量子化インデックスベクトルを照合する照合手段100のブロック図が示されており、次元決定手段101、量子化値取得手段102、103、尺度算出手段104とからなる。
 次元決定手段101は量子化値取得手段102、103へ接続され、決定された次元情報を出力する。量子化値取得手段102は、第1の量子化インデックスベクトルから、次元決定手段101から入力される次元の量子化インデックス値を取得し、第1の量子化インデックス値として尺度算出手段104へ出力する。量子化値取得手段103は、第2の量子化インデックスベクトルから、次元決定手段101から入力される次元の量子化インデックス値を取得し、第2の量子化インデックス値として尺度算出手段104へ出力する。尺度算出手段104は、量子化値取得手段102、103からそれぞれ出力される第1、第2の量子化インデックス値から同一性を表す尺度を算出し、出力する。
 次に、図20の照合手段100の動作について説明する。
 まず、照合手段100へは、第1の画像から抽出される量子化インデックスベクトルである第1の量子化インデックスベクトルと、第2の画像から抽出される量子化インデックスベクトルである第2の量子化インデックスベクトルとが入力される。入力された第1、第2の量子化インデックスベクトルは、それぞれ量子化値取得手段102、103へ入力される。
 量子化値取得手段102、103へは、次元決定手段101から出力される次元情報も入力される。次元決定手段101では、N次元ベクトルである量子化インデックスベクトルの各次元を指定する情報を順次出力する。出力する順序は必ずしも1からNまで1つずつ増えていく必要はなく、1からNまでの次元が過不足なく指定される順序であれば、どのような順序であってもよい。
 量子化値取得手段102、103では、入力された量子化インデックスベクトルから、次元決定手段101から出力される次元情報で指定される次元の量子化インデックス値を取得する。そして、取得した量子化インデックス値を尺度算出手段104へ出力する。
 尺度算出手段104では、量子化値取得手段102から出力される第1の量子化インデックス値と第2の量子化インデックス値とを比較する。この比較を各次元に対して行い、第1、第2の量子化インデックスベクトル間の類似尺度(あるいは距離尺度)を同一性尺度として算出する。
 得られた同一性尺度値は予め定めた閾値と比較し、同一性の判定を行う。同一性尺度が類似度をあらわす尺度である場合には、尺度値が閾値以上の場合に同一と判定する。一方、同一性尺度が距離をあらわす尺度である場合には、尺度値が閾値以下の場合に同一と判定する。
 次に、フローチャートを用いて図20の照合手段100の動作を説明する。まず、同一性尺度として類似度を用いる場合の動作について説明する。
 図21は、照合手段100の動作を示すフローチャートである。図21のフローチャートでは、量子化インデックスベクトルの次元(の番号)をnで表し、次元は1からNまでの合計N次元あるものとする。また、類似度を算出する変数をBで表すこととする。
 まず、次元決定手段101は、照合する量子化インデックスベクトルの最初の次元として、次元1を決定し(n=1)、量子化値取得手段102、103へ供給するとともに、尺度算出手段104において変数Bを0にセットする。(ステップS100)。
 次に、量子化値取得手段102、103において、第1の量子化インデックスベクトル、第2の量子化インデックスベクトルから、次元nの第1の量子化インデックス値と第2の量子化インデックス値とを取得し、尺度算出手段104へ供給する(ステップS102)。
 次に、尺度算出手段104において、第1の量子化インデックス値と第2の量子化インデックス値とから、それぞれの量子化インデックスに対応する特徴量の間の類似度ΔBを算出する(ステップS104)。例えば、量子化インデックスが一致する場合にはΔB=1とし、それ以外の場合はΔB=0とする。あるいは、量子化インデックスから量子化前の特徴量の代表値を算出し、代表値間の差分が小さいほど大きくなる値をΔBとして用いてもよい。この際、特徴量の代表値を算出して差分を求めるかわりに、量子化インデックス値の組み合わせによってΔBの値を引くことができるテーブルを保持しておき、量子化インデックス値の組み合わせからこのテーブルを用いてΔBの値を直接求めるようになっていてもよい。
 次に、ΔBの値は変数Bに加算される(ステップS106)。この際、ΔBの値が0の場合には、変数Bに0を加算するかわりに、加算しないように制御してもよい。
 次に、次元の番号nが次元数Nに到達したかどうかを調べ(ステップS108)、到達しない場合はステップS112へ移行し、到達した場合には、そのときの変数Bの値を同一性尺度(類似度を表す尺度)として出力し(ステップS110)、処理を終了する。
 ステップ112では、次元決定手段101が、取得する量子化インデックスの次元として、n=n+1によって次の次元を決定し、量子化値取得手段102、103へ供給する。そして、再度ステップS102へ移行する。
 なお、ここでは、次元1からNまで順番に抽出処理を行っているが、順番はこれに限らず任意でよい。
 次に、同一性尺度として距離を用いる場合の動作について説明する。
 図22は、照合手段100の動作を示す別のフローチャートである。図22のフローチャートでも、量子化インデックスベクトルの次元(の番号)をnで表し、次元は1からNまでの合計N次元あるものとする。また、距離尺度を算出する変数をCで表すこととする。
 基本的なフローは、図21の場合と同じであるが、ステップS100、S104、S106、S110がそれぞれステップS200、S204、S206、S210に置き換わっている点が異なる。
 まず、ステップS200では、次元決定手段101において、照合する量子化インデックスベクトルの最初の次元として、次元1を決定し(n=1)、量子化値取得手段102、103へ供給するとともに、尺度算出手段104において変数Cを0にセットする。
 ステップS204では、尺度算出手段104において、第1の量子化インデックス値と第2の量子化インデックス値とから、それぞれの量子化インデックスに対応する特徴量の距離ΔCを算出する。例えば、量子化インデックスが一致する場合にはΔC=0とし、それ以外の場合はΔC=1とする。あるいは、量子化インデックスから量子化前の特徴量の代表値を算出し、代表値間の差分が小さいほど小さくなる値をΔCとして用いてもよい。この際、特徴量の代表値を算出して差分を求めるかわりに、量子化インデックス値の組み合わせによってΔCの値を引くことができるテーブルを保持しておき、量子化インデックス値の組み合わせからこのテーブルを用いてΔCの値を直接求めるようになっていてもよい。
 ステップS206では、ΔCの値は変数Cに加算される。この際、ΔCの値が0の場合には、変数Cに0を加算するかわりに、加算しないように制御してもよい。
 ステップS210では、そのときの変数Cの値を同一性尺度(距離を表す尺度)として出力し、処理を終了する。
 それ以外のステップについては、図21の場合と同様である。ただし、ステップS108で次元の番号nが次元数Nに到達した場合にはステップS210へ移行する。
 なお、ここでは、次元1からNまで順番に抽出処理を行っているが、順番はこれに限らず任意でよい。
 次に、第1の量子化インデックス値と第2の量子化インデックス値とで、「共に量子化インデックスが0」である次元を除外し、同一性尺度として類似度を用いる場合の動作について説明する。
 図23は、照合手段100の動作を示す別のフローチャートである。図23のフローチャートでも、量子化インデックスベクトルの次元(の番号)をnで表し、次元は1からNまでの合計N次元あるものとする。また、類似度を算出する変数をBで表すこととし、「共に量子化インデックスが0」ではない次元をカウントするための変数をAで表すこととする。
 まず、次元決定手段101は、照合する量子化インデックスベクトルの最初の次元として、次元1を決定し(n=1)、量子化値取得手段102、103へ供給するとともに、尺度算出手段104において変数A、Bを0にセットし(ステップS300)、ステップS102へ移行する。
 ステップS102は図21の場合と同様であり、終了後、ステップS314へ移行する。
 ステップS314では、尺度算出手段104において、第1の量子化インデックス値と第2の量子化インデックス値とがともに0であるかどうかを調べる。ともに0である場合には、ステップS108へ移行し、どちらか一方が0でない場合には、変数Aの値をひとつ増やし(ステップS316)、ステップS104へ移行する。
 ステップS104、S106、S108、S112の処理は図21の場合と同様である。ステップS108で次元の番号nが次元数Nに到達した場合には、ステップS310へ移行する。
 ステップS310では、尺度算出手段104において、B/Aの値を算出し、同一性尺度として出力し、処理を終了する。ただし、A=0の場合には、規定の値(例えば0.5)を出力する。
 なお、ここでは、次元1からNまで順番に抽出処理を行っているが、順番はこれに限らず任意でよい。
 以上本発明の実施の形態について説明したが、本発明は、上述した実施形態に限定されるものではない。本発明の構成や詳細には、本発明の範囲内で当業者が理解しうる様々な変更をすることができる。また、本発明の画像識別子抽出装置は、その有する機能をハードウェア的に実現することは勿論、コンピュータとプログラムとで実現することができる。プログラムは、磁気ディスクや半導体メモリ等のコンピュータ可読記録媒体に記録されて提供され、コンピュータの立ち上げ時などにコンピュータに読み取られ、そのコンピュータの動作を制御することにより、そのコンピュータを前述した各実施の形態における次元決定手段、抽出領域取得手段、領域特徴量算出手段、比較手段、領域特徴量算出方法取得手段、比較方法取得手段として機能させる。
 なお、本発明は、日本国にて2009年3月13日に特許出願された特願2009-061023の特許出願、および日本国にて2009年4月14日に特許出願された特願2009-097862の特許出願に基づく優先権主張の利益を享受するものであり、当該特許出願に記載された内容は、全て本明細書に含まれるものとする。
1…次元決定手段
2…抽出領域取得手段
3、3A…領域特徴量算出手段
 31、31A…第1の領域特徴量算出手段
 32、32A…第2の領域特徴量算出手段
4、4B…比較手段
 41…大小比較手段
 42、44…量子化手段
 43…差分値算出手段
 45…量子化境界決定手段
5…領域特徴量算出方法取得手段
6…比較方法取得手段

Claims (45)

  1.  対をなす2つの部分領域の形状の組み合わせと、対をなす2つの部分領域の相対的な位置関係との双方が、他の少なくとも1つの部分領域対と相違する1以上の部分領域対を含む、画像中の複数の部分領域対に従って、画像の各部分領域から領域特徴量を抽出する抽出手段と、
     前記抽出した各部分領域毎の領域特徴量に基づいて、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる方法を使用して、前記画像の識別に用いる画像識別子を生成する生成手段と
    を備えることを特徴とする画像識別子抽出装置。
  2.  対をなす2つの部分領域の形状の組み合わせと、対をなす2つの部分領域の相対的な位置関係との双方が同じ部分領域対では、絶対的な位置が相違する
    ことを特徴とする請求項1に記載の画像識別子抽出装置。
  3.  前記生成手段は、部分領域対に対応する画像識別子要素を算出し、該画像識別子要素の集合を前記画像識別子とする
    ことを特徴とする請求項1または2に記載の画像識別子抽出装置。
  4.  前記生成手段は、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる方法を使用して、部分領域対に対応する画像識別子要素を算出し、該画像識別子要素の集合を前記画像識別子とする
    ことを特徴とする請求項1乃至3の何れか1項に記載の画像識別子抽出装置。
  5.  前記抽出手段は、対をなす2つの部分領域の双方から領域特徴量を抽出し、
     前記生成手段は、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる方法を使用して、前記対をなす2つの部分領域の双方から抽出された領域特徴量を用いて画像識別子要素を算出し、該画像識別子要素の集合を前記画像識別子とする
    ことを特徴とする請求項1乃至4の何れか1項に記載の画像識別子抽出装置。
  6.  前記抽出手段は、対をなす2つの部分領域の双方から領域特徴量を抽出し、
     前記生成手段は、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる方法を使用して、前記対をなす2つの部分領域の双方から抽出された領域特徴量の差分値から画像識別子要素を算出し、該画像識別子要素の集合を前記画像識別子とする
    ことを特徴とする請求項1乃至5の何れか1項に記載の画像識別子抽出装置。
  7.  前記抽出手段は、対をなす2つの部分領域の双方から領域特徴量を抽出し、
     前記生成手段は、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる方法を使用して、前記対をなす2つの部分領域の双方から抽出された領域特徴量の差分値を量子化して画像識別子要素を算出し、該画像識別子要素の集合を前記画像識別子とする
    ことを特徴とする請求項1乃至6の何れか1項に記載の画像識別子抽出装置。
  8.  前記抽出手段は、対をなす2つの部分領域の双方から領域特徴量を抽出し、
     前記生成手段は、少なくとも1つの部分領域対については他の1以上の部分領域対と異なるパラメータを使用して、前記対をなす2つの部分領域の双方から抽出された領域特徴量の差分値を量子化して画像識別子要素を算出し、該画像識別子要素の集合を前記画像識別子とする
    ことを特徴とする請求項1乃至7の何れか1項に記載の画像識別子抽出装置。
  9.  前記抽出手段は、対をなす2つの部分領域の双方から領域特徴量を抽出し、
     前記生成手段は、少なくとも1つの部分領域対については他の1以上の部分領域対と異なるパラメータである量子化閾値を使用して、前記対をなす2つの部分領域の双方から抽出された領域特徴量の差分値を量子化して画像識別子要素を算出し、該画像識別子要素の集合を前記画像識別子とする
    ことを特徴とする請求項1乃至8の何れか1項に記載の画像識別子抽出装置。
  10.  前記生成手段は、部分領域対を複数のタイプに分類し、タイプごとにそのタイプに属する部分領域対の量子化に用いるパラメータを決定する
    ことを特徴とする請求項7乃至9の何れか1項に記載の画像識別子抽出装置。
  11.  前記生成手段は、同じタイプに属する部分領域対ごとの双方の領域特徴量の差分値の分布に基づいて、同じタイプに属する部分領域対の量子化に用いるパラメータを決定する
    ことを特徴とする請求項10に記載の画像識別子抽出装置。
  12.  前記生成手段は、異なる量子化値に分類される部分領域対の割合がタイプ内でほぼ均等になるように、量子化に用いるパラメータを決定する
    ことを特徴とする請求項7または11に記載の画像識別子抽出装置。
  13.  前記生成手段は、3値に量子化する
    ことを特徴とする請求項10乃至12の何れか1項に記載の画像識別子抽出装置。
  14.  前記3値は、第1の量子化値、該第1の量子化値より小さな値の第2の量子化値、該第2の量子化値より小さな値の第3の量子化値から構成され、第1と第2の量子化値の差と、第2と第3の量子化値の差とが等しい
    ことを特徴とする請求項13に記載の画像識別子抽出装置。
  15.  前記生成手段は、部分領域対を構成する2つの部分領域のうちの第1の部分領域から抽出した領域特徴量をV1、第2の部分領域から抽出した領域特徴量をV2とし、量子化に用いるパラメータをthとするとき、V1からV2を減算した値がthより大きければ第1の量子化値を生成し、V1からV2を減算した値が-thより小さければ第3の量子化値を生成し、V1からV2を減算した値の絶対値がth以下であれば第2の量子化値を生成する
    ことを特徴とする請求項13または14に記載の画像識別子抽出装置。
  16.  前記生成手段は、同じタイプに属する部分領域対ごとの双方の領域特徴量の差分値の絶対値をソートして、その上位または下位から予め定められた割合の位置にある値を、当該タイプに属する部分領域対の量子化に用いるパラメータthとして用いる
    ことを特徴とする請求項15に記載の画像識別子抽出装置。
  17.  前記抽出手段は、画像の各部分領域の画素値の代表値を、当該部分領域の領域特徴量として抽出する
    ことを特徴とする請求項1乃至16の何れか1項に記載の画像識別子抽出装置。
  18.  前記抽出手段は、画像の各部分領域の平均輝度値を、当該部分領域の領域特徴量として抽出する
    ことを特徴とする請求項1乃至17の何れか1項に記載の画像識別子抽出装置。
  19.  部分領域が隣接または一部重複する部分領域対どうしでは、対をなす2つの部分領域の形状の組み合わせと、対をなす2つの部分領域の相対的な位置関係との双方が相違する
    ことを特徴とする請求項1乃至18の何れか1項に記載の画像識別子抽出装置。
  20.  前記抽出手段は、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる算出方法により前記領域特徴量を抽出する
    ことを特徴とする請求項1乃至19の何れか1項に記載の画像識別子抽出装置。
  21.  請求項1乃至20の何れかに記載の画像識別子抽出装置から生成された画像識別子を用いて照合を行うことを特徴とする照合装置。
  22.  請求項1乃至20の何れかに記載の画像識別子抽出装置から生成された画像識別子を用いて識別を行うことを特徴とする識別装置。
  23.  対をなす2つの部分領域の形状の組み合わせと、対をなす2つの部分領域の相対的な位置関係との双方が、他の少なくとも1つの部分領域対と相違する1以上の部分領域対を含む、画像中の複数の部分領域対に従って、画像の各部分領域から領域特徴量を抽出し、
     前記抽出した各部分領域毎の領域特徴量に基づいて、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる方法を使用して、前記画像の識別に用いる画像識別子を生成する
    ことを特徴とする画像識別子抽出方法。
  24.  対をなす2つの部分領域の形状の組み合わせと、対をなす2つの部分領域の相対的な位置関係との双方が同じ部分領域対では、絶対的な位置が相違する
    ことを特徴とする請求項23に記載の画像識別子抽出方法。
  25.  前記生成手段は、部分領域対に対応する画像識別子要素を算出し、該画像識別子要素の集合を前記画像識別子とする
    ことを特徴とする請求項23または24に記載の画像識別子抽出方法。
  26.  前記画像識別子の生成では、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる方法を使用して、部分領域対に対応する画像識別子要素を算出し、該画像識別子要素の集合を前記画像識別子とする
    ことを特徴とする請求項23乃至25の何れか1項に記載の画像識別子抽出方法。
  27.  前記領域特徴量の抽出では、対をなす2つの部分領域の双方から領域特徴量を抽出し、
     前記生成手段は、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる方法を使用して、前記対をなす2つの部分領域の双方から抽出された領域特徴量を用いて画像識別子要素を算出し、該画像識別子要素の集合を前記画像識別子とする
    ことを特徴とする請求項23乃至26の何れか1項に記載の画像識別子抽出方法。
  28.  前記領域特徴量の抽出では、対をなす2つの部分領域の双方から領域特徴量を抽出し、
     前記画像識別子の生成では、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる方法を使用して、前記対をなす2つの部分領域の双方から抽出された領域特徴量の差分値から画像識別子要素を算出し、該画像識別子要素の集合を前記画像識別子とする
    ことを特徴とする請求項23乃至27の何れか1項に記載の画像識別子抽出方法。
  29.  前記領域特徴量の抽出では、対をなす2つの部分領域の双方から領域特徴量を抽出し、
     前記画像識別子の生成では、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる方法を使用して、前記対をなす2つの部分領域の双方から抽出された領域特徴量の差分値を量子化して画像識別子要素を算出し、該画像識別子要素の集合を前記画像識別子とする
    ことを特徴とする請求項23乃至28の何れか1項に記載の画像識別子抽出方法。
  30.  前記領域特徴量の抽出では、対をなす2つの部分領域の双方から領域特徴量を抽出し、
     前記画像識別子の生成では、少なくとも1つの部分領域対については他の1以上の部分領域対と異なるパラメータを使用して、前記対をなす2つの部分領域の双方から抽出された領域特徴量の差分値を量子化して画像識別子要素を算出し、該画像識別子要素の集合を前記画像識別子とする
    ことを特徴とする請求項23乃至29の何れか1項に記載の画像識別子抽出方法。
  31.  前記領域特徴量の抽出では、対をなす2つの部分領域の双方から領域特徴量を抽出し、
     前記画像識別子の生成では、少なくとも1つの部分領域対については他の1以上の部分領域対と異なるパラメータである量子化閾値を使用して、前記対をなす2つの部分領域の双方から抽出された領域特徴量の差分値を量子化して画像識別子要素を算出し、該画像識別子要素の集合を前記画像識別子とする
    ことを特徴とする請求項23乃至30の何れか1項に記載の画像識別子抽出方法。
  32.  前記画像識別子の生成では、部分領域対を複数のタイプに分類し、タイプごとにそのタイプに属する部分領域対の量子化に用いるパラメータを決定する
    ことを特徴とする請求項29乃至31の何れか1項に記載の画像識別子抽出方法。
  33.  前記画像識別子の生成では、同じタイプに属する部分領域対ごとの双方の領域特徴量の差分値の分布に基づいて、同じタイプに属する部分領域対の量子化に用いるパラメータを決定する
    ことを特徴とする請求項32に記載の画像識別子抽出方法。
  34.  前記画像識別子の生成では、異なる量子化値に分類される部分領域対の割合がタイプ内でほぼ均等になるように、量子化に用いるパラメータを決定する
    ことを特徴とする請求項29または33に記載の画像識別子抽出方法。
  35.  前記画像識別子の生成では、3値に量子化する
    ことを特徴とする請求項32乃至34の何れか1項に記載の画像識別子抽出方法。
  36.  前記3値は、第1の量子化値、該第1の量子化値より小さな値の第2の量子化値、該第2の量子化値より小さな値の第3の量子化値から構成され、第1と第2の量子化値の差と、第2と第3の量子化値の差とが等しい
    ことを特徴とする請求項35に記載の画像識別子抽出方法。
  37.  前記画像識別子の生成では、部分領域対を構成する2つの部分領域のうちの第1の部分領域から抽出した領域特徴量をV1、第2の部分領域から抽出した領域特徴量をV2とし、量子化に用いるパラメータをthとするとき、V1からV2を減算した値がthより大きければ第1の量子化値を生成し、V1からV2を減算した値が-thより小さければ第3の量子化値を生成し、V1からV2を減算した値の絶対値がth以下であれば第2の量子化値を生成する
    ことを特徴とする請求項35または36に記載の画像識別子抽出方法。
  38.  前記画像識別子の生成では、同じタイプに属する部分領域対ごとの双方の領域特徴量の差分値の絶対値をソートして、その上位または下位から予め定められた割合の位置にある値を、当該タイプに属する部分領域対の量子化に用いるパラメータthとして用いる
    ことを特徴とする請求項37に記載の画像識別子抽出方法。
  39.  前記領域特徴量の抽出では、画像の各部分領域の画素値の代表値を、当該部分領域の領域特徴量として抽出する
    ことを特徴とする請求項23乃至38の何れか1項に記載の画像識別子抽出方法。
  40.  前記領域特徴量の抽出では、画像の各部分領域の平均輝度値を、当該部分領域の領域特徴量として抽出する
    ことを特徴とする請求項23乃至39の何れか1項に記載の画像識別子抽出方法。
  41.  部分領域が隣接または一部重複する部分領域対どうしでは、対をなす2つの部分領域の形状の組み合わせと、対をなす2つの部分領域の相対的な位置関係との双方が相違する
    ことを特徴とする請求項23乃至40の何れか1項に記載の画像識別子抽出方法。
  42.  前記領域特徴量の抽出では、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる算出方法により前記領域特徴量を抽出する
    ことを特徴とする請求項23乃至41の何れか1項に記載の画像識別子抽出方法。
  43.  請求項23乃至42の何れかに記載の画像識別子抽出方法により生成された画像識別子を用いて照合を行うことを特徴とする照合方法。
  44.  請求項23乃至42の何れかに記載の画像識別子抽出方法により生成された画像識別子を用いて識別を行うことを特徴とする識別方法。
  45.  コンピュータを、
     対をなす2つの部分領域の形状の組み合わせと、対をなす2つの部分領域の相対的な位置関係との双方が、他の少なくとも1つの部分領域対と相違する1以上の部分領域対を含む、画像中の複数の部分領域対に従って、画像の各部分領域から領域特徴量を抽出する抽出手段と、
     前記抽出した各部分領域毎の領域特徴量に基づいて、少なくとも1つの部分領域対については他の1以上の部分領域対と異なる方法を使用して、前記画像の識別に用いる画像識別子を生成する生成手段と
    して機能させるためのプログラム。
PCT/JP2010/001770 2009-03-13 2010-03-12 画像識別子抽出装置 WO2010103850A1 (ja)

Priority Applications (6)

Application Number Priority Date Filing Date Title
JP2011503733A JP4935960B2 (ja) 2009-03-13 2010-03-12 画像識別子抽出装置
EP10750606.5A EP2407932B1 (en) 2009-03-13 2010-03-12 Image identifier extraction device
US13/256,088 US8744193B2 (en) 2009-03-13 2010-03-12 Image signature extraction device
KR1020117021100A KR101335877B1 (ko) 2009-03-13 2010-03-12 화상 식별자 추출 장치
CN201080011278.9A CN102349093B (zh) 2009-03-13 2010-03-12 图像签名提取设备
KR1020137022201A KR101411127B1 (ko) 2009-03-13 2010-03-12 화상 식별자 추출 장치

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2009-061023 2009-03-13
JP2009061023 2009-03-13
JP2009097862 2009-04-14
JP2009-097862 2009-04-14

Publications (1)

Publication Number Publication Date
WO2010103850A1 true WO2010103850A1 (ja) 2010-09-16

Family

ID=42728147

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2010/001770 WO2010103850A1 (ja) 2009-03-13 2010-03-12 画像識別子抽出装置

Country Status (6)

Country Link
US (1) US8744193B2 (ja)
EP (1) EP2407932B1 (ja)
JP (4) JP4935960B2 (ja)
KR (2) KR101335877B1 (ja)
CN (1) CN102349093B (ja)
WO (1) WO2010103850A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2420973B1 (en) * 2009-04-14 2019-11-13 Nec Corporation Image identifier extraction device
US9870517B2 (en) * 2011-12-08 2018-01-16 Excalibur Ip, Llc Image object retrieval
JP6456031B2 (ja) * 2014-03-25 2019-01-23 キヤノン株式会社 画像認識装置、画像認識方法およびプログラム
WO2016004330A1 (en) * 2014-07-03 2016-01-07 Oim Squared Inc. Interactive content generation

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004192555A (ja) * 2002-12-13 2004-07-08 Fuji Xerox Co Ltd 情報管理方法、情報管理装置及び情報管理プログラム
JP2005031901A (ja) * 2003-07-10 2005-02-03 Nippon Conlux Co Ltd 硬貨識別装置および方法
JP2006268825A (ja) * 2005-02-28 2006-10-05 Toshiba Corp オブジェクト検出装置、学習装置、オブジェクト検出システム、方法、およびプログラム
JP2007293438A (ja) * 2006-04-21 2007-11-08 Omron Corp 特徴量取得装置

Family Cites Families (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2585606B2 (ja) 1987-07-06 1997-02-26 株式会社日立製作所 画像パターン検索方法
JPH02268388A (ja) * 1989-04-10 1990-11-02 Hitachi Ltd 文字認識方法
JPH0368272A (ja) 1989-08-05 1991-03-25 Ricoh Co Ltd 画像データ処理装置
US5436653A (en) * 1992-04-30 1995-07-25 The Arbitron Company Method and system for recognition of broadcast segments
DE69435133D1 (de) 1993-12-10 2008-10-09 Ricoh Kk Bildextrahierungsverfahren und Gerät, und Bilderkennungsverfahren und Gerät, zum Extrahieren und Erkennen eines bestimmten Bildes aus einem Bildeingabesignal
JP4458429B2 (ja) 1993-12-10 2010-04-28 株式会社リコー カラー画像認識方法および装置
JP3847832B2 (ja) * 1995-03-31 2006-11-22 キヤノン株式会社 画像情報処理方法及びその装置、並びにその制御方法
JP2002109596A (ja) * 2000-09-28 2002-04-12 Nippon Conlux Co Ltd 貨幣識別方法及び装置
JP3819236B2 (ja) 2000-12-13 2006-09-06 株式会社東芝 パターン認識方法及びパターン認識を行わせるプログラムを記憶したコンピュータ読み取り可能な記憶媒体
JP4190739B2 (ja) * 2001-02-19 2008-12-03 日本電気株式会社 画像特徴量生成装置および方法、ならびに画像特徴量生成プログラムを記録した記憶媒体
US20040012600A1 (en) * 2002-03-22 2004-01-22 Deering Michael F. Scalable high performance 3d graphics
JP2008529414A (ja) * 2005-01-28 2008-07-31 ユークリッド・ディスカバリーズ・エルエルシー ビデオデータを処理する装置および方法
US20060193520A1 (en) * 2005-02-28 2006-08-31 Takeshi Mita Object detection apparatus, learning apparatus, object detection system, object detection method and object detection program
JP2006351001A (ja) 2005-05-19 2006-12-28 Nippon Telegr & Teleph Corp <Ntt> コンテンツ特徴量抽出方法及び装置及びコンテンツ同一性判定方法及び装置
WO2007129374A1 (ja) * 2006-04-26 2007-11-15 Mitsubishi Denki Kabushiki Kaisha 物体検出装置及びエレベータの物体検出装置
EP2036354A1 (en) * 2006-06-20 2009-03-18 Koninklijke Philips Electronics N.V. Generating fingerprints of video signals
JP2008009938A (ja) 2006-06-30 2008-01-17 Omron Corp 動画像データ処理装置、動画像データ処理方法、動画像データ処理プログラム、およびこれを記録した記録媒体
EP2315161A1 (en) 2006-10-11 2011-04-27 Mitsubishi Electric Information Technology Centre Europe B.V. Image descriptor for image recognition
EP2390838B1 (en) * 2009-01-23 2019-10-02 Nec Corporation Image identifier extracting apparatus
US8170341B2 (en) * 2009-01-23 2012-05-01 Nec Corporation Image signature extraction device
EP2333718B1 (en) * 2009-01-29 2013-08-28 Nec Corporation Feature amount selecting device
WO2010103848A1 (ja) * 2009-03-13 2010-09-16 日本電気株式会社 画像識別子照合装置
WO2010146786A1 (ja) * 2009-06-16 2010-12-23 日本電気株式会社 画像識別子照合装置
JP5573131B2 (ja) * 2009-12-01 2014-08-20 日本電気株式会社 映像識別子抽出装置および方法、映像識別子照合装置および方法、ならびにプログラム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004192555A (ja) * 2002-12-13 2004-07-08 Fuji Xerox Co Ltd 情報管理方法、情報管理装置及び情報管理プログラム
JP2005031901A (ja) * 2003-07-10 2005-02-03 Nippon Conlux Co Ltd 硬貨識別装置および方法
JP2006268825A (ja) * 2005-02-28 2006-10-05 Toshiba Corp オブジェクト検出装置、学習装置、オブジェクト検出システム、方法、およびプログラム
JP2007293438A (ja) * 2006-04-21 2007-11-08 Omron Corp 特徴量取得装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DAIDI ZHONG ET AL.: "ADVANCES IN INFORMATION SYSTEMS; [LECTURE NOTES IN COMPUTER SCIENCE", 28 June 2007, SPRINGER BERLIN HEIDELBERG, article "Face Image Retrieval System Using Combination of Subimages", pages: 49 - 60
KOTA IWAMOTO ET AL.: "Response to the Call for Proposals on Video Signature Tools", 87. MPEG MEETING; 2-2-2009 - 6-2-2009; LAUSANNE; (MOTION PICTURE EXPERT GROUP OR ISO/IEC JTC1/SC29/WG11, 29 January 2009 (2009-01-29)
See also references of EP2407932A4

Also Published As

Publication number Publication date
EP2407932A4 (en) 2012-08-29
US8744193B2 (en) 2014-06-03
EP2407932A1 (en) 2012-01-18
US20120014609A1 (en) 2012-01-19
JP5598487B2 (ja) 2014-10-01
CN102349093B (zh) 2014-10-22
KR20110123771A (ko) 2011-11-15
KR20130101156A (ko) 2013-09-12
JP4935960B2 (ja) 2012-05-23
KR101335877B1 (ko) 2013-12-02
JP2012099156A (ja) 2012-05-24
JP5598486B2 (ja) 2014-10-01
EP2407932B1 (en) 2019-11-27
CN102349093A (zh) 2012-02-08
JP5454600B2 (ja) 2014-03-26
JPWO2010103850A1 (ja) 2012-09-13
KR101411127B1 (ko) 2014-06-25
JP2012123831A (ja) 2012-06-28
JP2012099157A (ja) 2012-05-24

Similar Documents

Publication Publication Date Title
JP5045817B2 (ja) 画像識別子抽出装置
JP4735783B2 (ja) 画像識別子照合装置
JP5573863B2 (ja) 画像識別子抽出装置
JP4735786B2 (ja) 画像識別子照合装置
JP5454600B2 (ja) 画像識別子抽出装置
JP5686120B2 (ja) 画像識別子抽出装置

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 201080011278.9

Country of ref document: CN

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

Ref document number: 10750606

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2011503733

Country of ref document: JP

ENP Entry into the national phase

Ref document number: 20117021100

Country of ref document: KR

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 13256088

Country of ref document: US

Ref document number: 2010750606

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE