WO1998020454A1 - Appareil d'extraction de structure - Google Patents

Appareil d'extraction de structure Download PDF

Info

Publication number
WO1998020454A1
WO1998020454A1 PCT/JP1996/003212 JP9603212W WO9820454A1 WO 1998020454 A1 WO1998020454 A1 WO 1998020454A1 JP 9603212 W JP9603212 W JP 9603212W WO 9820454 A1 WO9820454 A1 WO 9820454A1
Authority
WO
WIPO (PCT)
Prior art keywords
pattern
dimensional
pattern data
fourier
free
Prior art date
Legal status (The legal status 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 status listed.)
Ceased
Application number
PCT/JP1996/003212
Other languages
English (en)
French (fr)
Inventor
Koji Kobayashi
Hiroshi Nakajima
Takafumi Aoki
Masayuki Kawamata
Tatsuo Higuchi
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Azbil Corp
Original Assignee
Azbil Corp
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 Azbil Corp filed Critical Azbil Corp
Priority to PCT/JP1996/003212 priority Critical patent/WO1998020454A1/ja
Priority to US09/297,537 priority patent/US6195460B1/en
Priority to DE69631845T priority patent/DE69631845T2/de
Priority to JP08536385A priority patent/JP3035654B2/ja
Priority to EP96935527A priority patent/EP1011073B1/en
Publication of WO1998020454A1 publication Critical patent/WO1998020454A1/ja
Priority to KR1019997003884A priority patent/KR100295586B1/ko
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; 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/431Frequency domain transformation; Autocorrelation

Definitions

  • an N-dimensional pattern for example, sound (one-dimensional), planar image (two-dimensional), three-dimensional (three-dimensional)] is collated based on spatial frequency characteristics, and a registered pattern and an illuminated pattern are compared.
  • the present invention relates to a pattern extracting device for extracting a difference from a moving pattern and a moving pattern.
  • the difference between two similar patterns has been extracted by a human visual check. That is, one pattern is set as a reference pattern, and the difference is extracted by comparing the reference pattern with the other pattern with human eyes. Extraction of common patterns (moving patterns) at different positions in two similar patterns (whole patterns) is also performed by human visual check. That is, one pattern is set as a reference pattern, and the reference pattern and the other pattern are compared with human eyes to extract a moving pattern.
  • human visual checks can only deal with differences between two similar patterns or when the movement pattern is clear. In other words, if the entire pattern is composed of complex patterns, or if the differences or movement patterns are very small, it takes time to extract and accurate checks Can not do. In addition, it is difficult to detect all movement patterns on a platform where there are multiple movement patterns.
  • An object of the present invention is to provide a pattern extracting device capable of extracting a difference between similar patterns and a moving pattern in a short time and accurately.
  • the present invention provides a registered Fourier N-dimensional pattern data by performing an N-dimensional discrete Fourier transform on the N-dimensional pattern data of a registered pattern, and the N-dimensional pattern data of a collation pattern.
  • N-dimensional discrete Fourier transform is applied to the data to create a matching Fourier N-dimensional pattern data, and the registered Fourier N-dimensional pattern data and the matching Fourier N-dimensional pattern
  • the first pattern processing means performs N-dimensional discrete Fourier transform and N-dimensional discrete inverse Fourier on the combined free-dimensional N-dimensional pattern data obtained by combining them. Either one of the Fourier transforms is performed, the correlation peak in the correlation component area appearing in the synthesized Fourier N-dimensional pattern data subjected to the Fourier transform is obtained, and the obtained correlation peak is obtained.
  • an N-dimensional discrete free-run is performed on the synthesized free-range data whose part is masked and a part of which is masked.
  • the synthesized Fourier N-dimensional pattern data subjected to the linear transformation and the registered Fourier N-dimensional pattern data are re-synthesized, and the re-combined Fourier N-dimensional pattern data obtained as a result is recombined.
  • N-dimensional discrete inverse Fourier transform is performed.
  • the N-dimensional pattern data of the registered pattern is subjected to the N-dimensional discrete Fourier transform to generate the registered free N-dimensional pattern data, and the N-dimensional pattern data of the collation pattern is generated.
  • the N-dimensional discrete Fourier transform is applied to the input data to create a matching Fourier N-dimensional pattern data.
  • the registered Fourier N-dimensional pattern data and the matching Fourier N-dimensional pattern data are synthesized, and the first pattern is obtained from the synthesized Fourier processing N-dimensional pattern data obtained by this.
  • the N-dimensional discrete Fourier transform or the N-dimensional discrete inverse Fourier transform is performed by the Fourier transform means, and the correlation component appearing in the synthesized Fourier N-dimensional pattern data subjected to the Fourier transform is obtained.
  • the correlation peak in the end is obtained. Then, the periphery including and including the correlation peak is masked, and the synthesized Fourier N-dimensional pattern data in which a part of the correlation peak is masked is subjected to N-dimensional pattern processing by the first pattern processing means. If a discrete Fourier transform has been applied, an N-dimensional inverse discrete Fourier transform is applied, and the N-dimensional inverse discrete Fourier transform is applied in the first pattern processing means. If so, an N-dimensional discrete Fourier transform is performed.
  • the synthesized Fourier N-dimensional pattern data subjected to the Fourier transform and the registered Fourier N-dimensional pattern data are re-synthesized, and the re-synthesized Fourier ⁇ N-dimensional N-dimensional discrete inverse Fourier transform is applied to the pattern data.
  • the recombined hood subjected to the inverse Fourier transform obtained by this is In the N-dimensional pattern data, the outline of the difference is extracted, and the outline of the movement pattern is extracted. It is possible to know where the difference and the movement pattern are and where they exist.
  • the present invention further combines the registered Fourier N-dimensional pattern data and the illuminated Fourier N-dimensional pattern data, and obtains an amplitude with respect to the Taisei Fourier N-dimensional pattern data obtained thereby.
  • the second pattern processing means After performing suppression processing (1 og processing, T processing, etc.), one of ⁇ -dimensional discrete Fourier transform and ⁇ -dimensional discrete inverse Fourier transform is performed, and the second pattern processing means performs Performs amplitude restoration processing (inverse function processing such as 1 og processing and processing) on the synthesized Fourier ⁇ -dimensional pattern data that has been subjected to the Rier transform, and then synthesizes the data after the amplitude restoration processing has been performed.
  • the free N-dimensional pattern data and the registered Fourier N-dimensional pattern data are re-synthesized, and the re-composed Fourier N-dimensional pattern data obtained by this is N-dimensional discrete inverse Fourier Perform the conversion.
  • the present invention further synthesizes registered Fourier N-dimensional pattern data and collated Fourier N-dimensional pattern data, and converts the resulting data into Taisei Fourier N-dimensional pattern data.
  • amplitude suppression processing (10 g processing and processing, etc.) on either of them
  • one of N-dimensional discrete free-time transform and N-dimensional discrete inverse-free transform is performed, and the second pattern is applied.
  • the reconstructed Fourier obtained by recombining the Fourier transformed N-dimensional pattern data and the registered Fourier N-dimensional pattern data subjected to the Fourier transform by the processing means. Performs N-dimensional discrete inverse Fourier transform on N-dimensional pattern data.
  • the present invention further provides an N-dimensional discrete Fourier transform to the N-dimensional pattern data of the registered pattern, and then performs amplitude suppression processing (1 og processing, ⁇ processing, and the like), thereby obtaining a registered free space.
  • N-dimensional discrete Fourier transform to the N-dimensional pattern data of the registered pattern, and then performs amplitude suppression processing (1 og processing, ⁇ processing, and the like), thereby obtaining a registered free space.
  • the amplitude reconstruction processing (10 g processing and the inverse function processing such as "processing" is performed on the regenerated free-flowing N-dimensional pattern data.
  • N-dimensional discrete inverse Fourier transform is performed.
  • the present invention further provides an N-dimensional discrete file transformation of the N-dimensional pattern data of the registration pattern, and then performs an amplitude suppression process (1 og process, T process, etc.), thereby making the registration file available.
  • the matching Fourier N-dimensional pattern data is created, and the combined Fourier N-dimensional pattern data subjected to Fourier transform by the second pattern processing means and the registered Fourier N-dimensional pattern data are re-synthesized. Then, an N-dimensional discrete inverse Fourier transform is applied to the recombined Fourier N-dimensional pattern data obtained by this.
  • FIG. 1A to 1G are diagrams for explaining a process of extracting a difference between a registered pattern and a collation pattern in the pattern extraction device shown in FIG.
  • FIG. 2 is a block diagram of a pattern extracting apparatus according to an embodiment of the present invention.
  • FIG. 3 is a flowchart for explaining the registration (registration) pattern registration operation in the pattern extraction device shown in FIG.
  • FIG. 4 is a flowchart for explaining an operation of extracting a difference between a registered pattern and a matching pattern in the pattern extracting apparatus shown in FIG.
  • FIG. 5 is a flowchart, following FIG. 4, for explaining the operation of extracting the difference between the registered pattern and the matching pattern.
  • 6A to 6G are diagrams for explaining a process of extracting a moving pattern in the pattern extraction device shown in FIG.
  • FIG. 7 is a flowchart for explaining the movement pattern extraction operation in the pattern extraction device shown in FIG.
  • FIG. 8 is a flowchart following FIG. 7 for explaining the movement pattern extraction operation. It is.
  • FIG. 9 is a functional block diagram of the pattern extraction algorithm corresponding to the flowcharts shown in FIGS.
  • FIG. 10 is a functional block diagram of the pattern extraction algorithm corresponding to the flowcharts shown in FIGS. 5. BEST MODE FOR CARRYING OUT THE INVENTION
  • FIG. 2 is a block diagram of a pattern extracting apparatus according to an embodiment of the present invention, and illustrates a case where two-dimensional pattern data composed of image data is collated.
  • 10 is an operation unit
  • 20 is a control unit
  • the operation unit 10 is a numeric keypad 10-1 and a display (LCD: Liquid Crystal Display) 10-2.
  • CCD Charge Coupled Device
  • the control unit 20 includes a control unit 20 — 1 having a CPU (Central Processing Unit), a ROM (Ead Only Memory) 20 — 2, and a RAM (Random Access Memory) 20 — 3, a disk (HD) 20-4, a frame memory (FM) 20-5, an external connection (I / F) 20-6, and a Fourier transform section (FFT) 20-7 and a pattern extraction program is stored in R 0 M 20 -2.
  • a control unit 20 — 1 having a CPU (Central Processing Unit), a ROM (Ead Only Memory) 20 — 2, and a RAM (Random Access Memory) 20 — 3, a disk (HD) 20-4, a frame memory (FM) 20-5, an external connection (I / F) 20-6, and a Fourier transform section (FFT) 20-7 and a pattern extraction program is stored in R 0 M 20 -2.
  • a control unit 20 — 1 having a CPU (Central Processing Unit), a ROM (Ead Only Memory) 20
  • the reference pattern (registered pattern) is registered as shown in Fig. 3. That is, before starting the pattern matching, the user inputs the ID number assigned to the reference pattern using the numeric keypad 10-1 (step S301), and inputs the ID number to the CCD camera. Place the registration pattern at a predetermined position in the field of view of the camera 10-3. As a result, the original image of the registered pattern is converted to A
  • control unit 20 It is given to the control unit 20 as a grayscale image (image data: two-dimensional pattern data) with 0 pixels and 256 gradations.
  • the control unit 20-1 captures the image data of the registration pattern given by the operation unit 10 via the frame memory 20-5 (step S302), and captures the captured image data.
  • the image data of the registered pattern (see Fig. 1A) is sent to the Fourier transform unit 20-7 to be subjected to a two-dimensional discrete Fourier transform (two-dimensional DFT: two-dimensional discrete Fourier transform). Step S304).
  • the image data of the registered pattern shown in FIG. 1A becomes Fourier image data (registered Fourier image data) F A as shown in FIG. 1B.
  • the control unit 20-1 receives the Fourier image data F A as the original image data of the registered pattern and inputs the data into the hard disk 20-4.
  • the file is filed in correspondence with the 1D number (step S305).
  • the two-dimensional discrete Fourier transform is described in, for example, “Introduction to Computer Image Processing, Japan Industrial Technology Center, pp. 44-45” (Reference 1).
  • the difference between the registered pattern and the matching pattern is extracted as shown in FIG. That is, the user inputs the ID number assigned to the reference pattern using the numeric keypad 10-1 (step S401), and enters the field of view of the CCD camera 10-1. Place the matching pattern at a predetermined position in the range. As a result, the original image data of the matching pattern is
  • grayscale image image data: two-dimensional pattern data
  • control unit 20-1 receives the registration pattern filed in the disk 20-14. Then, the Fourier image data FA of the registered pattern corresponding to the ID number is read (step S).
  • the control unit 20-1 captures the image data of the collation pattern provided from the operation unit 10 via the frame memory 20-5 (step S403), and captures the image data.
  • the image data of the matching pattern (see Fig. 1C) is sent to the Fourier transform unit 20-7 to perform a two-dimensional discrete Fourier transform (two-dimensional DFT) (step S405).
  • the image data of the collation pattern shown in FIG. 1C becomes Fourier image data (collation Fourier image data) FB as shown in FIG. 1D.
  • the control unit 20-1 receives the Fourier image data FB of the collation pattern obtained in step S405 and the registered pattern read out in step S402.
  • the free image data FA is synthesized (step S406) to obtain the synthesized free image data.
  • the matching Fourier image data is assumed to be A e0
  • the registered free image data is assumed to be B ⁇ . — Represented by 0).
  • A, B, ⁇ , and ⁇ are also functions of frequency (four ⁇ ) space (u, V).
  • a 'B' e J ( e — 0) A 'B' cos (0- ⁇ ) + j-A-B-sin ( ⁇
  • Step S After obtaining the image data, perform amplitude suppression processing using the phase-only correlation method (Step S).
  • 10 g processing is performed as the amplitude suppression processing. That is, an arithmetic expression of a synthetic full one Li E image data described above A 'B' e j - take the log of (0 3 ⁇ 4?), Log ( A 'B) - ej' ⁇ - in ⁇ and the child Thus, the amplitude ⁇ ⁇ ⁇ ⁇ ⁇ is suppressed to log (A ⁇ B) (A'B> log (A'B)).
  • the above-mentioned phase-only correlation method is a cross-correlation modified so as to focus on a spatial phase change of an image, and is a synthetic Fourier in which amplitude information is suppressed and limited to only phase information. Image data is required.
  • the influence of the illuminance difference between when the registered pattern is collected and when the collation pattern is collected is reduced. That is, by performing the amplitude suppression processing, the spectrum intensity of each pixel is suppressed, and there is no extraordinary value, so that more information is effective.
  • 1 og processing is performed as amplitude suppression processing, but ⁇ processing may be performed. Further, the processing is not limited to the 10 g processing and the ⁇ processing, and any processing may be used as long as the amplitude can be suppressed.
  • a constant value for all amplitudes with amplitude suppression For example, when the value is set to 1, that is, when only the phase is used, there is an advantage that the amount of calculation can be reduced and an amount of data can be reduced as compared with log processing or “processing”.
  • control unit 20-1 After performing the amplitude suppression processing in step, the control unit 20-1 sends the synthesized Fourier image data subjected to the amplitude suppression processing to the Fourier transformation unit 20-7, and performs the second 2 A dimensional DFT is performed (step S408), whereby the synthesized free image data on which the amplitude suppression processing has been performed is a synthesized free image as shown in Fig. 1E. It becomes image data.
  • the control unit 20-1 captures the synthesized Fourier image data obtained in step S408, and obtains a predetermined image including the central part from the synthesized Fourier image data.
  • the pixel with the highest intensity (correlation peak) in the rear is extracted (step S409). In this case, a correlation peak appears near the center of the correlation component identifier.
  • the control unit 20-1 masks the periphery including the correlation peak extracted in step S409 (step S410). That is, as shown in FIG. 1F, the area SO surrounded by the white dotted line is masked with respect to the combined free-flow image data shown in FIG. 1E. Then, the synthesized Fourier image data in which this region SO is masked is subjected to a two-dimensional discrete inverse Fourier transform (two-dimensional IDFT) (step S). 411), the amplitude restoration processing is performed on the synthesized Fourier image data on which the two-dimensional IDFT has been performed (step S412).
  • the amplitude restoration process is a process in which the inverse function of the function performed in the amplitude suppression process in step S407 is performed on the amplitude. If 2 and 1 og e A, then e A.
  • control unit 20-1 transmits the synthesized Fourier image data subjected to the amplitude restoration processing in step S412 and the registered file read out in step S402.
  • the re-combined Fourier image data FA is recombined with the re-combined Fourier image data (step S4113), and the two-dimensional IDFT is performed on the recombined Fourier image data (step S413).
  • S 4 14) to obtain recombined free-view image data as shown in FIG. 1G.
  • the re-synthesis in step S413 means registration pattern B and collation pattern A.
  • the contour of the pattern that exists only in the matching pattern appears at the corresponding position. .
  • the pattern of the car is superimposed on a part of the matching pattern.
  • the pattern of this car appears in Fig. 1G as a difference between the verification pattern and the registered pattern.
  • the control unit 20-1 extracts the pattern of the vehicle appearing in FIG. 1G as a pattern existing only in the matching pattern (step S4 15) o
  • the registered pattern does not have the car pattern superimposed thereon, and the car pattern is superimposed only on the matching pattern.
  • the case has been described.
  • the case where the car pattern is superimposed on both the registered pattern and the collation pattern, and the position of the car pattern is moving will be described below.
  • control unit 20-1 performs the same processing as the flowchart shown in FIGS. 4 and 5, and extracts the movement pattern through the processing steps shown in FIG.
  • FIG. 6E which is a process corresponding to FIG. 1E
  • a correlation value P 1 indicating a background and a correlation value P 2 indicating a car appear.
  • a mask is applied to the periphery including the correlation peak P 1 (see FIG. 6F).
  • the control unit 20-1 performs a two-dimensional IDFT on the masked composite free image data (step S411), and performs the composite Fourier subjected to the two-dimensional IDFT.
  • An amplitude restoration process is performed on the image data (step S412).
  • the synthesized Fourier image data subjected to the amplitude restoration processing is read out in step S402.
  • the registered registered free image data FA is recombined (step S4113), and a two-dimensional IDFT is performed on the recombined free image data obtained thereby (step S4).
  • Step S4 14) obtains recombined Fourier image data as shown in FIG. 6G.
  • the contours of the moving pattern existing in both the registered pattern and the matching pattern are represented by: Appears at the corresponding position in the matching pattern.
  • the position of the vehicle in the registered pattern shown in FIG. 6A is moving.
  • This car pattern appears in Fig. 6G as a moving pattern that exists at different positions in the matching pattern and the registered pattern.
  • the control unit 20-1 extracts the pattern of the vehicle shown in FIG. 6G as a moving pattern existing at different positions of the registered pattern and the matching pattern (step). Step S 4 16).
  • the two-dimensional IDFT is performed in the Fourier transform units 20-7, but may be performed in the CPU 20-1.
  • the two-dimensional DFT is performed in step S408 shown in FIG. 4, but the two-dimensional DFT may be performed instead of the two-dimensional DFT. That is, instead of performing two-dimensional DFT on the synthesized Fourier image data on which the amplitude suppression processing has been performed, two-dimensional IDFT may be performed.
  • two-dimensional IDFT may be performed in step D408.
  • two-dimensional DFT is performed in step S411.
  • the matching accuracy of the two-dimensional DFT and the two-dimensional IDFT does not change quantitatively.
  • the two-dimensional IDFT is described in Reference 1 above.
  • the 2-dimensional DFT in scan STEP S 4 0 8 may be subjected to amplitude suppression processing and then synthesized. That is, as shown in FIG. 7, the step S407 of FIG. 4 is eliminated, and the step S702 of reading out the registered Fourier image data FA and the image data of the collation pattern are eliminated. A step S703 for performing the amplitude suppression processing is provided between the step S704 for performing the input.
  • a step S706 of applying two-dimensional DFT to the image data of the collation pattern and a collation Fourier after the amplitude suppression processing are combined with the image data FB and the registered Fourier image data FA.
  • a step S707 for performing amplitude suppression processing is provided between the step S708 and the step S708.
  • the combined Fourier image data subjected to the two-dimensional IDFT and the Fourier image data of the registered pattern read out in step 720 A step S714 for performing amplitude restoration processing is provided after the step S711 for recombining the signals.
  • steps S701, S709 to S712, and S717 to S717 are the steps S401 and S410 shown in FIGS. 9 to S 4 12 and S 17 to S 4 17 are the same as those described above, and a description thereof will be omitted.
  • the suppression rate of the amplitude of the synthesized free image data at this time is smaller than the case where the amplitude is suppressed after the synthesized free image data shown in FIGS. Therefore, the method of performing amplitude suppression processing on the synthesized Fourier image data shown in Fig. 4 is more effective than the method of performing amplitude suppression processing on the synthesized Fourier image data shown in Fig. 7 and generating the synthesized Fourier image data.
  • the position of the difference and the matching accuracy of the movement pattern are improved. It should be noted that even when the amplitude suppression processing shown in FIGS. 7 and 8 is performed before the combined Fourier image data is obtained, two-dimensional! ) Instead of FT, two-dimensional IDFT may be performed.
  • the two-dimensional pattern extraction processing has been described.
  • the three-dimensional pattern extraction processing can be performed in the same manner, regardless of the two-dimensional or three-dimensional pattern.
  • the multidimensional pattern extraction processing can be performed in the same manner.
  • the amplitude suppression processing is performed, the amplitude suppression processing need not always be performed. Further, when all amplitudes are set to 1 in the amplitude suppression processing, that is, when the phase is limited, the amplitude restoration processing may not be performed.
  • Fig. 9 shows a functional block diagram of the pattern extraction algorithm corresponding to the flow charts shown in Figs. 4 and 5, and Fig. 10 shows the flow chart shown in Figs. 7 and 8.
  • a functional block diagram of the pattern extraction algorithm corresponding to Fig. 4 is shown.
  • each functional block has the same step number as each step of the flowchart, and corresponds to the appended step number. The function of the step Each has.
  • the registered pattern and the matching pattern are matched based on the spatial frequency characteristics, and the difference between similar patterns and the movement pattern are extracted as the matching result. This makes it possible to perform quality inspection, abnormality detection (analysis), and moving object detection in a short time and accurately.

Landscapes

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

Description

明 細 書 パターン抽出装置
1 . 技術分野
こ の発明は、 空間周波数特性に基づいて N次元のパターン 〔例えば、 音声 ( 1 次元) 、 平面画像 ( 2次元) 、 立体 ( 3次元) 〕 の照合を行い、 登録パタ ー ンと 照台パター ンと の相違点や移動パター ンを抽出するパター ン抽出装置に関する。
2 . 背景技術
従来、 2 つの類似したパター ンの相違点の抽出は、 人間の目視によるチェ ッ ク によ って行われている。 すなわち、 一方のパターンを基準パターンと し、 この基 準パターン と他方のパターンとを人間の目で比較して相違点を抽出している。 また、 2つの類似したパターン (全体パター ン) 中のそれぞれ異なる位置に存 在する共通パター ン (移動パターン) の抽出も、 人間の目視によるチェ ッ ク によ つて行われている。 すなわち、 一方のパター ンを基準パター ンと し、 この基準パ ターンと他方のパターンとを人間の目で比較して、 移動パタ一ンを抽出している。 しかしながら、 人間の目視によるチ ェ ッ クでは、 2つの類似したパター ンの相 違点や移動パターンが明らかなものである場合に しか対応できない。 すなわち、 全体パタ一ンが複雑なパタ一ンで構成されていたり、 あるいは相違点や移動パ夕 一ンが微小な場合には、 抽出するのに時間がかかると と もに正確なチェ ッ クが行 えない。 また、 移動パター ンが複数存在する場台には、 全ての移動パターンを検 出する こ とが難しい。
本発明の目的は、 短時間でしかも正確に、 類似パターン間の相違点や移動パタ ー ンを抽出する こ とのできるパターン抽出装置を提供する こ とにある。
3 . 発明の開示
本発明は、 登録パターンの N次元パター ンデ一タに N次元離散的フ一 リ ヱ変換 を施して登録フー リ エ N次元パター ンデータを作成し、 照合バタ一 ンの N次元パ 夕一ンデータに N次元離散的フ一リ ェ変換を施して照合フー リ エ N次元パター ン デ一タを作成し、 登録フー リ ェ N次元パターンデータ と照合フー リ ェ N次元パタ ー ンデータ とを合成し、 これによつて得られる合成フ ー リ ェ N次元パター ンデー 夕に対して第 1 のパターン処理手段において N次元離散的フ一 リ ェ変換および N 次元離散的逆フー リ エ変換の何れか一方を施し、 こ のフー リ エ変換の施された合 成フー リ エ N次元パター ンデータに出現する相関成分エ リ ア内の相関ピークを求 め、 この求めた相関ピークを含みその周辺をマスク し、 これによ つてその一部が マス ク された合成フ一 リ ヱ N次元パター ンデ一タに対し、 第 1 のパターン処理手 段において N次元離散的フ一 リ ェ変換が施されている場合には N次元離散的逆フ — リ ヱ変換を施し、 第 1 のパターン処理手段において N次元離散的逆フー リ エ変 換が施されている場合には N次元離散的フー リ ヱ変換を施し、 こ のフー リ ヱ変換 の施された合成フ一 リ ェ N次元パター ンデータ と登録フー リ ェ N次元パター ンデ —タ とを再合成し、 これによつて得られる再合成フー リ エ N次元パター ンデータ に対して N次元離散的逆フ一 リ ェ変換を施す。
本発明によれば、 登録パタ一ンの N次元パターンデータに N次元離散的フ一 リ ェ変換が施されて登録フ 一 リ ェ N次元パター ンデータが作成され、 照合パター ン の N次元パタ 一 ンデータ に N次元離散的フ一 リ ェ変換が施されて照合フー リ エ N 次元パターンデータが作成される。 そ して、 登録フー リ エ N次元パターンデータ と照合フ一 リ エ N次元パター ンデータ とが合成され、 これによつて得られる合成 フー リ 工 N次元パタ ー ンデータ に対して第 1 のパタ ー ン処理手段において N次元 離散的フー リ ヱ変換あるいは N次元離散的逆フ一 リ ヱ変換が施され、 このフー リ ェ変換の施された合成フー リ ェ N次元パター ンデータに出現する相関成分ヱ リ 了 内の相関ピークが求められる。 そ して、 この相関ピークを含みその周辺がマ ス ク され、 これによつてその一部がマス ク された合成フ一 リ エ N次元パターンデータ に対し、 第 1 のパターン処理手段において N次元離散的フ一 リ ヱ変換が施されて いる場合には N次元離散的逆フー リ エ変換が施され、 第 1 のパターン処理手段に おいて N次元離散的逆フー リ エ変換が施されている場合には N次元離散的フー リ ェ変換が施される。 そして、 こ のフー リ エ変換の施された合成フー リ エ N次元パ ターンデータ と登録フー リ ェ N次元パター ンデータ とが再合成され、 これによつ て得られる再合成フー リ ヱ N次元パターンデータに対して N次元離散的逆フ一 リ ェ変換が施される。 これによつて得られる逆フー リ ェ変換の施された再合成フー リ エ N次元パター ンデータ には、 相違点の輪郭が抽出され、 また移動パター ンの 輪郭が抽出され、 相違点や移動パター ンが何で合っても、 どこに存在するかが分 かる。
また、 本発明は、 さらに登録フー リ エ N次元パター ンデータ と照台フー リ エ N 次元パター ンデータ とを合成し、 これによつて得られる台成フー リ エ N次元パタ — ンデータ に対して振幅抑制処理 ( 1 o g処理や T処理等) を行ったうえ Ν次元 離散的フー リ エ変換および Ν次元離散的逆フー リ エ変換の何れか一方を施し、 第 2のパタ一ン処理手段によってフー リ エ変換の施された合成フ一 リ エ Ν次元パタ ー ンデータ に対して振幅復元処理 ( 1 o g処理や 処理等の逆関数処理) を行つ たうえ、 この振幅復元処理の行われた合成フ一 リ ェ N次元パターンデータ と登録 フー リ ェ N次元パター ンデータ とを再合成し、 これによつて得られる再台成フ ー リェ N次元パタ一ンデータに対して N次元離散的逆フー リ ェ変換を施す。
また、 本発明は、 さ らに登録フー リ エ N次元パタ ー ンデータ と照合フー リ エ N 次元パター ンデータ とを合成し、 これによつて得られる台成フー リ ェ N次元パ夕 — ンデータに対して振幅抑制処理 ( 1 0 g処理や 処理等) を行ったうえ N次元 離散的フ 一 リ ュ変換および N次元離散的逆フ ー リ ュ変換の何れか一方を施し、 第 2のパター ン処理手段によってフ ー リ ヱ変換の施された台成フ ー リ ェ N次元パタ — ンデータ と登録フー リ エ N次元パター ンデータ とを再合成し、 こ れによ つて得 られる再合成フー リ エ N次元パター ンデータに対して N次元離散的逆フ一 リ エ変 換を施す。
また、 本発明は、 さらに登録パター ンの N次元パター ンデータに N次元離散的 フー リ ェ変換を施してから振幅抑制処理 ( 1 o g処理や^処理等) を行う こ と に より登録フ ー リ ェ N次元パター ンデータを作成し、 照合パター ンの N次元パター ンデータに N次元離散的フ一 リ エ変換を施してから振幅抑制処理 ( l o g処理や 処理等) を行う こ とによ り照合フー リ ヱ N次元パター ンデータを作成し、 第 2 のパター ン処理手段によってフ ー リ ヱ変換の施された台成フー リ ェ N次元パター ンデ一夕 と登録フー リ エ N次元パタ ー ンデ一夕 とを再合成し、 これによ つ て得ら れる再台成フ ー リ ヱ N次元パター ンデータ に対して振幅復元処理 ( 1 0 g処理や ; "処理等の逆関数処理) をを行ってから N次元離散的逆フー リ エ変換を施す。 また、 本発明は、 さらに登録パター ンの N次元パター ンデータに N次元離散的 フ一 リ ヱ変換を施してから振幅抑制処理 ( 1 o g処理や T処理等) を行う ことに より登録フ一 リ ェ Ν次元パター ンデータを作成し、 照合パター ンの Ν次元パター ンデータに Ν次元離散的フー リ ェ変換を施してから振幅抑制処理 ( 1 o g処理や ^処理等) を行う ことにより照合フー リ エ N次元パターンデータを作成し、 第 2 のパター ン処理手段によってフー リ ェ変換の施された合成フ一リ エ N次元パ夕一 ンデータ と登録フー リ エ N次元パタ一 ンデータ とを再合成し、 これによ つて得ら れる再合成フー リ ェ N次元パター ンデータに対して N次元離散的逆フー リ ェ変換 を施す。
4 . 図面の簡単な説明
図 1 A〜 1 Gは、 図 2 に示すパタ一ン抽出装置における登録パタ一ンと照合パ ター ン との相違点の抽出過程を説明する図である。
図 2 は、 本発明の一実施例を示すパターン抽出装置のプロ ッ ク構成図である。 図 3 は、 図 2 に示すパター ン抽出装置における基準 (登録) パタ一ンの登録動 作を説明するためのフローチヤ一 トである。
図 4 は、 図 2 に示すパター ン抽出装置における登録パター ンと照合パター ン と の相違点の抽出動作を説明するためのフ ローチ ャー トである。
図 5 は、 図 4 に続いて、 登録パター ン と照合パター ンと の相違点の抽出動作を 説明するためのフローチヤ一 トである。
図 6 A ~ 6 Gは、 図 2 に示すパター ン抽出装置における移動パター ンの抽出過 程を説明する図である。
図 7 は、 図 2 に示すパターン抽出装置における移動パターンの抽出動作を説明 するためのフ ローチヤ一トである。
図 8 は、 図 7 に続いて、 移動パターンの抽出動作を説明するためのフローチヤ — トである。 である。
図 9 は、 図 4 および図 5 に示したフ ローチ ヤ 一 トに対応するパター ン抽出アル ゴリ ズムの機能ブロ ッ ク図である。
図 1 0 は、 図 7 および図 8 に示したフローチ ャー ト に対応するパター ン抽出ァ ルゴリ ズムの機能プロ ッ ク図である。 5. 発明を実施するための最良の形態
以下、 本発明を図面に基づき詳細に説明する。
図 2 は本発明の一実施例を示すパター ン抽出装置のブロ ッ ク構成図であり、 画 像データからなる 2次元パター ンデータを照合する場合を説明する。 図 2 におい て、 1 0 は操作部、 2 0 はコ ン ト ロール部であり、 操作部 1 0 はテ ンキー 1 0 — 1, デ ィ ス プ レイ ( L C D : Liquid Crystal Display) 1 0 — 2 およ び C C D (Charge Coupled Device) カ メ ラ 1 0 — 3 を備えている。 コ ン ト ロール部 2 0 は、 C P U (Central Processing Unit) を有してなる制御部 2 0 — 1 と、 R O M (E ead Only Memory) 2 0 — 2 と、 R A M (Random Access Memory) 2 0 — 3 と、 ー ドディ ス ク ( H D ) 2 0 — 4 と、 フ レーム メ モ リ ( F M) 2 0 — 5 と、 外部接 続部 ( I / F ) 2 0 — 6 と、 フー リ エ変換部 ( F F T ) 2 0 — 7 とを備えてなり、 R 0 M 2 0 - 2 にはパター ン抽出プロ グラ ムが格納されている。
〔基準パタ一ンの登録〕
このパター ン抽出装置において基準パター ン (登録パター ン) は図 3 に示すよ うにして登録される。 すなわち、 パター ン照合の開始前に、 利用者は、 テ ンキー 1 0 - 1 を用いて基準パター ンに割り当てられた I Dナ ンバを入力のうえ (ス テ ップ S 3 0 1 ) 、 C C Dカ メ ラ 1 0 — 3 の視野範囲の所定の位置に登録パター ン を置く。 これによ り、 登録パタ一ンの原画像が、 A / D変換により、 3 2 0 X 4
0 0画素, 2 5 6階調の濃淡画像 (画像データ : 2次元パター ンデータ) と して、 コ ン ト ロール部 2 0へ与えられる。
制御部 2 0 — 1 は、 この操作部 1 0 よ り与えられる登録パター ンの画像データ をフ レーム メ モ リ 2 0 — 5 を介して取り込み (ス テ ッ プ S 3 0 2 ) 、 取り込んだ 登録パターンの画像データ (図 1 A参照) をフー リ エ変換部 2 0 — 7 へ送って 2 次元離散的フー リ ヱ変換 ( 2次元 D F T : two-dimensional Discrete Fourier T ransform) を施す (ステ ッ プ S 3 0 4 ) 。 これによ り、 図 1 Aに示された登録パ タ一 ンの画像データは、 図 1 Bに示されるよ うなフー リ エ画像データ (登録フー リ ェ画像データ) F Aとなる。 制御部 2 0 — 1 は、 このフー リ ェ画像デ一タ F A を登録パター ンの原画像データ と して、 ハー ドディ スク 2 0 — 4内に入力された
1 Dナ ンバと対応させてフ ア イ ル化する (ス テ ッ プ S 3 0 5 ) 。 なお、 2次元離散的フー リ エ変換については、 例えば 「コ ン ピュータ画像処理 入門、 日本工業技術セ ンタ ー編、 P P. 4 4 — 4 5」 (文献 1 ) 等に説明されて いる o
〔相違点の抽出〕
このパターン抽出装置において、 登録パターン と照合パター ン との相違点の抽 出は、 図 4 に示すよ う に して行われる。 すなわち、 利用者は、 テ ンキー 1 0 — 1 を用いて基準パター ンに割り当てられた I Dナ ンバを入力のうえ ( ス テ ッ プ S 4 0 1 ) 、 C C Dカ メ ラ 1 0— 1 の視野範囲の所定の位置に照合バター ンを置く。 これによ り、 登録パター ンの場合と同様にして、 照合パター ンの原画像データが、
3 2 0 X 4 0 0 画素, 2 5 6階調の濃淡画像 (画像データ : 2次元パタ ー ンデ一 タ) と して、 コ ン ト ロール部 2 0 へ与えられる。
制御部 2 0 — 1 は、 テ ンキー 1 0 — 1 を介して I Dナ ンバが与えられる と、 ノ、 一ドディ ス ク 2 0 一 4内にフ ァ イ ル化されている登録パタ ー ンから、 その I Dナ ンバに対応する登録パター ンのフー リ ェ画像データ F Aを読み出す (ス テ ッ プ S
4 0 2〉 。 次に、 制御部 2 0 — 1 は、 操作部 1 0 よ り与えられる照合パタ一 ンの 画像データをフ レームメ モ リ 2 0 — 5を介して取り込み (ステ ッ プ S 4 0 3 ) 、 取り込んだ照合パター ンの画像データ (図 1 C参照) をフー リヱ変換部 2 0 - 7 へ送って 2次元離散的フー リ エ変換 ( 2次元 D F T) を施す (ステ ッ プ S 4 0 5 ) 。 これによ り、 図 1 Cに示された照合パター ンの画像データは、 図 1 Dに示される よ う なフー リ エ画像データ (照合フー リ エ画像データ) F B となる。
次に、 制御部 2 0 — 1 は、 ス テ ッ プ S 4 0 5で得られた照合パター ンのフー リ ェ画像データ F B と ステ ッ プ S 4 0 2で読み出した登録パター ンのフ一 リ ェ画像 データ F Aとを合成し (ス テ ッ プ S 4 0 6 ) 、 合成フ一リ ェ画像データを得る。 こ こで、 合成フ一リ ェ画像データ は、 照合フー リ ェ画像データ を A · e 』 0 と し、 登録フ一 リ エ画像データを B · と した場合、 A ' B · e ' (0— 0 )で表される。 但し、 A, B, Θ, φヒも周波数 (フー リヱ) 空間 ( u, V ) の関数とする。 そ して、 A · B · e ' W¾6 >は、
A ' B ' e J ( e — 0 ) = A ' B ' c o s ( 0 - ^ ) + j - A - B - s i n ( θ
- Φ ) · · · ( 1 ) と して表され、 Α · e 』 0 = a , + j /3 , 、 Β · e ' φ = a 2 + ) β 2 とすると、
A = ( α ,2+ ^ ,2) 1/2
Β = ( α 2 2 + β 2 2) 1/2
θ = t a τι ' { β ι / a , )
φ = \. Ά η ~ ι { β ι / 2 )
となる。 この ( 1 ) 式を計算する こ とにより合成フー リ エ画像データを得る。 なお、
A , B . e "0—¾6 ) = A , B . e . e—
= Α · e 1 θ · Β · e ' j φ
= ( t + i β , ) · ( a 2 — j >S 2 )
= ( a i · 2 + β , · β 2 ) + i ( a 2 · S > - a >
• β z ) · · · ( 2 ) と して、 式 ( 2 〉 を用いて合成フー リ エ画像データを求めるようにしてもよい。 制御部 2 0 — 1 は、 このよ うに して合成フ一 リ ヱ画像データを得た後、 位相限 定相関法 (Phase Only Correlation) を用いて振幅抑制処理を行う (ス テ ッ プ S
4 0 7 ) 。 本実施例では、 振幅抑制処理と して、 1 0 g処理を行う。 すなわち、 前述した合成フ一リ ェ画像データの演算式である A ' B ' e j (0 -?¾ )の l o gを とり、 l o g ( A ' B ) - e j ' θ - φ とするこ とによ り、 振幅である Α · Βを l o g ( A · B ) に抑制する ( A ' B > l o g ( A ' B ) ) 。 上述した位相限定 相関法は、 画像の空間的な位相変化に着目するよ う に修正した相互相関 (cross correlation) であ り、 振幅情報が抑制されて位相情報のみに限定された合成フー リ ェ画像デ一タが求められる。
振幅抑制処理を施した合成フ一 リ ェ画像データでは登録パターンの採取時と照 合パター ンの採取時の照度差による影響が小さ く なる。 すなわち、 振幅抑制処理 を行う こ とによ り、 各画素のスぺク ト ラ ム強度が抑圧され、 飛び抜けた値がなく なり、 よ り多く の情報が有効となる。
なお、 この実施例では、 振幅抑制処理と して 1 o g処理を行う ものと したが、 ^処理を行う よ うにしてもよい。 また、 1 0 g処理や^処理に限らず、 振幅を抑 制するこ とができればどのよ うな処理でもよい。 振幅抑制で全ての振幅を一定値、 例えば 1 にする と、 すなわち位相のみにする と、 l o g処理や "処理等に比べ、 計算量を減らすこ とができるという利点とデータが少なく なるという利点がある。 ス テ ッ プ S 4 0 7 で振幅抑制処理を行つた後、 制御部 2 0 — 1 は、 その振幅抑 制処理を行つた合成フー リ ェ画像データをフ ー リ ェ変換部 2 0 — 7 へ送り、 第 2 回目の 2次元 D F Tを施す (ステ ッ プ S 4 0 8 ) 。 これによ り、 振幅抑制処理の 行われた合成フ一リ ェ画像デ一タは、 図 1 E に示されるよ うな合成フ一リ ェ画像 データ とな る。
次に、 制御部 2 0 — 1 は、 ス テ ッ プ S 4 0 8で得られた合成フー リ エ画像デ一 夕を取り込み、 この合成フ一 リ ヱ画像データ よ り中心部分を含む所定の相関成分 ェ リ ァの各画素の相関成分の強度 (振幅) をスキ ャ ン し、 各画素の相関成分の強 度の ヒ ス ト グラ ムを求め、 この ヒ ス ト グラ ムよ り相関成分ヱ リ ア内の最も強度の 高い画素 (相関ピーク) を抽出する (ス テ ッ プ S 4 0 9 ) 。 この場合、 相関成分 ェ リ ァの中心付近に相関ピークが現れる。
次に、 制御部 2 0 — 1 は、 ス テ ッ プ S 4 0 9 で抽出した相関ピークを含みその 周辺をマスクする (ステッ プ S 4 1 0 ) 。 すなわち、 図 1 Fに示すよう に、 図 1 Eに示される合成フ ー リ ヱ画像データに対し、 白い点線で囲んだ領域 S O をマ ス クする。 そ して、 この領域 S O がマ ス ク された合成フー リ エ画像データに対して 2次元離散的逆フー リ ヱ変換 ( 2次元 I D F T : two-dimensional Inverse DFT) を施し (ス テ ッ プ S 4 1 1 ) 、 こ の 2次元 I D F Tの施された合成フー リ エ画像 データ に対して振幅復元処理を行う (ス テ ッ プ S 4 1 2 ) 。 こ こで、 振幅復元処 理は、 ス テ ッ プ S 4 0 7 での振幅抑制処理で行われた関数の逆関数を振幅に対し て行う処理のこ とを言い、 と した場合には A 2 と し、 1 o g e Aと した場合 には e A とする。
次に、 制御部 2 0 - 1 は、 ステ ッ プ S 4 1 2 にて振幅復元処理の行われた合成 フー リ エ画像データ とス テ ッ プ S 4 0 2で読み出されている登録フー リ エ画像デ ータ F Aとを再合成し (ステ ッ プ S 4 1 3 ) 、 これによつて得られる再合成フー リ ェ画像デ一タに対して 2次元 I D F Tを施し (ス テ ッ プ S 4 1 4 ) 、 図 1 Gに 示すよ うな再合成フ ーリ ェ画像データを得る。
こ こで、 ス テ ッ プ S 4 1 3 での再合成とは、 登録パター ン Bと照合パタ ー ン A 9 とを合成したものから照合パター ン Aを抽出する意味である。 具体的には次の 2 通り考えられる。 この 2つはどち ら も等価である。
①ス テ ッ プ S 4 0 6での合成を A · Β · e J (0 -i¾)と した場合、 ス テ ッ プ S 4 1 3での再合成は、 (Α · B ) Z B · e ' 1 ίθ ' φ' + φι Ά - e j 6> とする。
②ス テ ッ プ S 4 0 6での合成を A · B · e j (0— 0 )と した場合、 ス テ ッ プ S 4 1 3での再合成は、 (A · B ) ZB · e 1 'φ ' (θ -φ , > =Α · e J0 とする。
図 1 Gを見ても分かるよ う に、 こ の 2次元 I D F Tの施された再合成フー リ エ 画像データには、 照合パタ一ンのみに存在するパター ンの輪郭がそれと対応する 位置に現れる。 すなわち、 図 1 Cに示す照合パター ンでは、 その一部に自動車の パター ンが重畳されている。 この自動車のパター ンが照合パター ン と登録パター ソとの相違点と して図 1 Gに現れる。 制御部 2 0 — 1 は、 図 1 Gに現れている自 動車のパター ンを、 照合パター ンのみに存在するパター ンと して抽出する (ス テ ップ S 4 1 5 ) o
[移動パターンの抽出〕
上述した登録パター ンと照合パターンとの相違点の抽出の場合、 登録パター ン には自動車のパタ一ンが重畳されておらず、 照合パタ一ンにのみ自動車のパタ一 ンが重畳されている場合について説明した。 これに対し、 登録パター ンにも照合 パタ ー ンに も自動車のパタ ー ンが重畳されており、 こ の自動車のパタ ー ンの位置 が移動している場合について以下に説明する。
この場合にも、 制御部 2 0 — 1 は、 図 4および図 5 に示したフ ローチ ャー ト と 同様の処理を行い、 図 6 に示すような処理過程を経て移動パターンを抽出する。 但し、 この場合、 図 1 Eに対応する処理過程である図 6 Eにおいては、 背景を示 す相関値 P 1 と自動車を示す相関値 P 2が出現する。 この場合、 相関値 P 1 の方 が大き く、 相関ピーク と して抽出されるので、 こ の相関ピーク P 1 を含みその周 辺にマ ス クが施される (図 6 F参照) 。
制御部 2 0 — 1 は、 このマス ク された合成フ一 リ ェ画像データに対して 2次元 I D F Tを施し (ス テ ッ プ S 4 1 1 ) 、 この 2次元 I D F Tの施された合成フー リ ェ画像データに対して振幅復元処理を行う (ス テ ッ プ S 4 1 2 ) 。 そ して、 こ の振幅復元処理の行われた合成フー リ ェ画像デ一タ とステ ッ プ S 4 0 2で読み出 された登録フ一 リ ェ画像データ F Aとを再合成し (ステツ プ S 4 1 3 ) 、 これに よつて得られる再合成フー リ ェ画像デ一タに対して 2次元 I D F Tを施し (ス テ ップ S 4 1 4 ) 、 図 6 Gに示すよ うな再合成フー リ ェ画像データを得る。
図 6 Gを見ても分かるように、 この 2次元 I D F Tの施された再合成フー リ エ 画像データ には、 登録パタ ー ン と照合パター ンと の両方に存在する移動パター ン の輪郭が、 照合パター ンのそれと対応する位置に現れる。 すなわち、 図 6 Cに示 す照合パター ンでは、 図 6 Aに示す登録パター ン中の自動車の位置が移動してい る。 こ の自動車のパター ンが照合パター ンと登録パター ンのそれぞれ異なる位置 に存在する移動パター ン と して図 6 Gに現れる。 制御部 2 0 ― 1 は、 図 6 Gに現 れている自動車のパター ンを、 登録パター ン と照合パター ンのそれぞれ異なる位 置に存在する移動パタ一 ンと して抽出する ( ス テ ッ プ S 4 1 6 ) 。
なお、 本実施例では、 2次元 I D F Tをフー リ ヱ変換部 2 0 — 7 において行う ものと したが、 C P U 2 0 — 1 内で行う ものと してもよい。 また、 本実施例では、 図 4 に示したス テ ッ プ S 4 0 8 にて 2次元 D F Tを行う ようにしたが、 2次元 D F Tではな く 2次元 I D F Tを行うよ う にしてもよい。 すなわち、 振幅抑制処理 の施された合成フー リ ェ画像データに対して 2次元 D F Tを行うのに代えて、 2 次元 I D F Tを行うようにしてもよい。 但し、 ス テ ッ プ D 4 0 8で 2次元 I D F Tを行った場合、 ス テ ッ プ S 4 1 1 では 2次元 D F Tを行う。 2次元 D F Tと 2 次元 I D F T とは、 定量的にみて照合精度は変わらない。 2次元 I D F Tについ ては、 先の文献 1 に説明されている。
また、 本実施例では、 合成後のフ一リ ヱ画像デー夕に対してス テ ッ プ S 4 0 7 で振幅抑制処理を施した後に、 ス テ ッ プ S 4 0 8 で 2次元 D F Tを行うよ うにし たが、 合成前の登録フー リ エ画像データ F Aおよび照合フー リ エ画像データ F B にそれぞれ振幅抑制処理を行った後に合成するよ う にしてもよい。 すなわち、 図 7 に示すよ う に、 図 4のステ ッ プ S 4 0 7をな く し、 登録フー リ エ画像データ F Aを読み出すス テ ッ プ S 7 0 2 と照合パター ンの画像データ の入力を行う ス テ ツ プ S 7 0 4 との間に振幅抑制処理を行う ステ ッ プ S 7 0 3 を設ける。 さ らに、 照 合パタ―ンの画像データに 2次元 D F Tを施すス テ ッ プ S 7 0 6 と振幅抑制処理 後の照合フー リ ヱ画像データ F B と登録フー リ ェ画像データ F Aとの合成を行う ステ ッ プ S 7 0 8 との間に振幅抑制処理を行う ス テ ッ プ S 7 0 7 を設ける。 但し、 この場合、 図 8 に示すよ うに、 2次元 I D F Tの施された合成フー リ エ画像デー タ とス テ ッ プ 7 0 2で読み出された登録パタ一 ンのフー リ ェ画像データ とを再合 成するス テ ッ プ S 7 1 3 の後に、 振幅復元処理を行う ス テ ッ プ S 7 1 4 を設ける。 このよ う に した場合、 ス テ ッ プ S 7 0 3 および S 7 0 7での振幅抑制処理によ つて、 振幅抑制処理の施された登録フー リ ヱ画像データおよび照合フ一 リ ヱ画像 データが得られ、 これらのフ一 リ ヱ画像データが合成されて合成フー リ ヱ画像デ —タが得られる。 なお、 ス テ ッ プ S 7 0 1, S 7 0 9〜 S 7 1 2, S 7 1 7 ~ S 7 1 7 は図 4 および図 5 に示すス テ ッ プ S 4 0 1, S 4 0 9 ~ S 4 1 2 , S 1 7 ~ S 4 1 7 と同一であり、 その説明を省略する。
この時の合成フ一 リ ヱ画像データ の振幅の抑制率は、 図 4 および図 5 に示す合 成フー リ ェ画像データ と してから振幅抑制処理を行う場合に対して小さい。 した がって、 図 4 に示す合成フー リ ェ画像データ と してから振幅抑制処理を行う方が、 図 7 に示す振幅抑制処理を行ってから合成フー リ エ画像データ とする方法に比べ て、 その相違点の位置や移動パター ンの照合精度がア ッ プする。 なお、 図 7 およ び図 8 に示す振幅抑制処理を行つてから合成フー リ エ画像データ とする場合にも、 合成フ一 リ ェ画像データに対して 2次元!) F Tではな く、 2次元 I D F Tを行う ように してもよい。
また、 上述した実施例においては、 2次元パター ンの抽出処理について説明し たが、 3次元パタ一ンの抽出処理についても同様にして行う ことが可能であり、 2次元, 3次元に拘わらず多次元のパター ンの抽出処理を同様にして行う こ とが できる。 また、 振幅抑制処理を行う ものと したが、 振幅抑制処理は必ずし も行わ なく てもよい。 また、 振幅抑制処理で全ての振幅を 1 と した場合、 すなわち位相 限定を行った場合には、 振幅復元処理を行わな く てもよい。
図 9 は図 4 および図 5 に示したフ ローチ ャー ト に対応するパター ン抽出アルゴ リ ズムの機能ブロ ッ ク図を示し、 図 1 0 は図 7 および図 8 に示したフ ローチ ヤ一 ト に対応するパタ ー ン抽出アルゴ リ ズムの機能ブロ ッ ク図を示す。 図 9 および図 1 0 において、 各機能ブロ ッ ク にはフ ローチ ヤ一 ト の各ス テ ッ プと同一のス テ ツ プ番号が付されており、 付されたステ ッ プ番号に対応するステッ プの機能をそれ ぞれ有する。
以上説明したよ うに本発明は、 登録パターンと照合パターンとが空間周波数特 性に基づいて照合され、 この照合結果と して類似パターン間の相違点や移動パタ ― ンの抽出を行う こ とが可能となり、 短時間でしかも正確に品質検査や異常検知 (解析) , 移動物体の検出を行う こ とができるよ う になる。

Claims

請 求 の 範 囲
1 . 登録パター ンの N次元パタ一 ンデータに N次元離散的フ一リ ヱ変換を施して 登録フー リ エ N次元パター ンデータを作成する登録フ一 リ エパター ンデータ作成 手段と、
照台パタ一ンの N次元パターンデータに N次元離散的フー リ ェ変換を施して照 合フ一 リ ェ N次元パターンデータを作成する照合フ一 リ ェパターンデータ作成手 段と、
前記登録フ一 リ ヱパターンデータ作成手段によつて作成された登録フ一 リ ェ N 次元パター ンデータ と前記照合フ一 リ エパター ンデータ作成手段によ つて作成さ れた照合フ一リ ェ N次元パター ンデータ とを合成して得られる台成フ一リ エ N次 元パターンデータに対して N次元離散的フー リ ェ変換および N次元離散的逆フ一 リ エ変換の何れか一方を施す第 1 のパター ン処理手段と、
前記第 1 のパターン処理手段によってフー リ ェ変換が施された合成フー リ エ N 次元パターンデータに出現する相関成分ェ リ ァ内の相関ピ一クを求め、 求めた相 関ピー クを含む所定の領域をマ ス クするマス ク処理手段と、
前記マ ス ク処理手段によ って所定領域がマ ス ク された台成フー リ エ N次元バタ ー ンデータに対し、 前記第 1 のパターン処理手段において N次元離散的フ一リ ェ 変換が施されている場合には N次元離散的逆フー リ ヱ変換を施し、 前記第 1 のパ ターン処理手段において N次元離散的逆フー リ ェ変換が施されている場合には N 次元離散的フ一 リ ェ変換を施す第 2 のパターン処理手段と、
前記第 2 のパター ン処理手段によ つてフー リ ェ変換の施された合成フー リ N 次元パター ンデータ と前記登録フー リ ヱパター ンデータ作成手段によつて作成さ れた登録フー リ エ N次元パター ンデータ とを再合成して得られる再合成フー リ ェ N次元パター ンデータ に対して N次元離散的逆フー リ ヱ変換を施す第 3 のパター ン処理手段と
を備えたこ とを特徴とするパターン抽出装置。
2 . 前記第 1 のパター ン処理手段は、
前記登録フー リ ヱパターンデータ作成手段によつて作成された登録フ一 リ ェ N 次元パタ一 ンデータ と前記照合フー リ ェパタ一ンデータ作成手段によつて作成さ れた照合フー リ ェ N次元パター ンデータ とを合成するパター ンデータ合成手段と、 前記パター ンデータ合成手段によって得られる合成フー リ ェ N次元パター ンデ -夕に対して振幅抑制処理を行う振幅抑制処理手段と、
前記振幅抑制処理手段によつて振幅抑制された合成フ一 リ N次元パターンデ 一夕に対して N次元離散的フー リ ェ変換および N次元離散的逆フー リ ュ変換の何 れか一方を施す第 1 のフ一 リ エ変換手段と
を備え、
前記第 3 のパターン処理手段は、
前記第 2 のパターン処理手段によってフー リ ェ変換された合成フー リ ヱ N次元 パターンデ—タに対して振幅復元処理を行う振幅復元処理手段と、
前記振幅抑制処理手段によ つて振幅復元された合成フ一 リ ェ N次元パターンデ 一夕 と前記登録フ一 リ ヱパター ンデータ作成手段によ つて作成された登録フ一 リ ェ N次元パター ンデータ とを再合成するパターンデータ再合成手段と、
前記パタ一ンデータ再合成手段によつて得られた再合成フ一 リ ヱ N次元パタ一 ンデータに対して N次元離散的逆フー リ ヱ変換を施す第 2 のフー リ ヱ変換手段と を備えたこ とを特徴とする請求の範囲第 1 項記載のパターン抽出装置。
3 . 前記第 1 のパター ン処理手段は、
前記登録フー リ エパター ンデータ作成手段によ つて作成された登録フ一 リ N 次元パター ンデータ と前記照合フ一 リ ヱパター ンデータ作成手段によつて作成さ れた照合フー リ ェ N次元パター ンデータ とを合成するパターンデ一タ合成手段と、 前記パター ンデータ合成手段によって得られた合成フー リ エ N次元パター ンデ —タに対して振幅抑制処理を行う振幅抑制処理手段と、
前記振幅抑制処理手段によ つて振幅抑制された合成フー リ N次元パターンデ —タに対して N次元離散的フー リ ェ変換および N次元離散的逆フー リ ェ変換の何 れか一方を施す第 1 のフー リ ヱ変換手段と
を備え、
前記第 3 のパターン処理手段は、
前記第 2 のパターン処理手段によ つてフー リ ェ変換された合成フー リ ェ N次元 パター ンデ一タ と前記登録フ ー リ エパター ンデータ作成手段によつて作成された 登録フ ー リ ェ N次元パター ンデータ とを再合成するパター ンデータ再合成手段と、 前記パターンデータ再合成手段によつて得られた再合成フ 一 リ ヱ N次元パター ンデータに対して N次元離散的逆フ ー リ ヱ変換を施す第 2 のフー リ ヱ変換手段を を備えたこ とを特徵とする請求項第 1 項記載のパターン抽出装置。
4 . 前記登録フ ー リ ヱパタ一 ンデータ作成手段は、
登録パター ンの N次元パター ンデータ に N次元離散的フー リ ヱ変換を施す第 1 のフ ー リ ュ変換手段と、
前記第 1 の フ 一 リ ヱ変換手段によつてフー リ ヱ変換された登録パターンの N次 元パタ一 ンデータに対し振幅抑制処理を行つて登録フー リ ェ N次元パタ一ンデー タを作成する第 1 の振幅抑制処理手段と
を備え、
前記照合フ ー リ ヱパタ ー ンデータ作成手段は、
照合パター ンの N次元パター ンデータに N次元離散的フ一リ ェ変換を施す第 2 のフ ー リ α 変換手段と、
前記第 2 のフ 一 リ ヱ変換手段によ つてフー リ ヱ変換された照合パタ一 ンの N次 元パター ンデータに対し振幅抑制処理を行つて照合フ一リ ヱ N次元パター ンデー タを作成する第 2 の振幅抑制処理手段と
を備え、
前記第 3 のパターン処理手段は、
前記第 2 のパターン処理手段によってフー リ ェ変換の施された合成フ 一 リ ェ N 次元パターンデータ と前記登録フ 一 リ ェパターンデータ作成手段によつて作成さ れた登録フー リ エ N次元パター ンデータ とを再合成するパター ンデータ再合成手 段と、
前記パターンデータ再合成手段によつて得られた再合成フ 一 リ ェ N次元パター ンデータに対して振幅復元処理を行う振幅復元処理手段と、
前記振幅復元処理手段によって振幅復元された再合成フ ー リ エ N次元パター ン データに対して N次元離散的逆フ ー リ ェ変換を施す第 3 の フ ー リ ェ変換手段と を備えたこ とを特徴とする請求項第 1 項記載のパターン抽出装置。
5 . 前記登録フ 一 リ エパター ンデータ作成手段は、
登録パターンの N次元パター ンデータに N次元離散的フー リ ェ変換を施す第 1 のフー リ ヱ変換手段と、
前記第 1 のフ一 リ ヱ変換手段によ ってフー リ ヱ変換された登録パタ ー ンの N次 元パタ一ンデータに対し振幅抑制処理を行って登録フー リ ヱ N次元パタ一ンデー タを作成する第 1 の振幅抑制処理手段と
を備え、 前記照合フ一 リ ェパター ンデータ作成手段は、
照合パター ンの N次元パター ンデータに N次元離散的フ一 リ ヱ変換を施す第 2 のフ一 リ 変換手段と、
前記第 2 のフ ー リ ェ変換手段によつてフ 一 リ ェ変換された照合パターンの N次 元パター ンデータに対し振幅抑制処理を行つて照台フー リ エ N次元パター ンデー タを作成する第 2 の振幅抑制処理手段と
を備え、
前記第 3 のパターン処理手段は、
前記第 2 のパターン処理手段によ ってフ一 リ ェ変換された合成フ一 リ ヱ N次元 パター ンデータ と前記登録フー リ エパター ンデータ作成手段によ つて作成された 登録フ ー リ ェ N次元パタ一ンデータ とを再合成するパター ンデータ再合成手段と、 前記パター ンデータ再合成手段によ つて得られた再合成フ一 リ ェ N次元パタ一 ンデータ に対して N次元離散的逆フ一 リ ヱ変換を施す第 3 のフー リ エ変換手段と を備えたこ とを特徴とする請求項第 1 項記載のパターン抽出装置。
6 . 前記マス ク処理手段は、
前記第 1 のパタ一ン処理手段によつてフ一 リ ヱ変換が施された合成フ一 リ エ N 次元パター ンデータの相関成分の強度のヒ ス ト グラ ムに基づいて相関成分ヱ リ ァ 内の相関ピークを求める相関ビーク算出手段と、
前記相関ピーク算出手段によつて求められた相関ピークを含む所定の領域をマ ス クするマ ス ク手段と
を備えたこ とを特徴とする請求項第 1 項記載のパターン抽出装置。
7 . 前記第 3 のパター ン処理手段によつて逆フー リ ェ変換された再合成フ一 リ ェ N次元パター ンデータから照台パター ンにのみに存在する相違パター ンを抽出す るパター ン抽出手段をさ らに備えたこ とを特徴とする請求項第 1項記載のパター ン抽出装置。
8 . 前記第 3 のパター ン処理手段によ つて逆フー リ ヱ変換された再合成フー リ エ N次元パター ンデータから照合パター ンと登録パター ンのそれぞれ異なる位置に 存在する移動パター ンを抽出するパター ン抽出手段をさ らに備えたこ とを特徴と する請求項第 1 項記載のパターン抽出装置。
PCT/JP1996/003212 1996-11-01 1996-11-01 Appareil d'extraction de structure Ceased WO1998020454A1 (fr)

Priority Applications (6)

Application Number Priority Date Filing Date Title
PCT/JP1996/003212 WO1998020454A1 (fr) 1996-11-01 1996-11-01 Appareil d'extraction de structure
US09/297,537 US6195460B1 (en) 1996-11-01 1996-11-01 Pattern extraction apparatus
DE69631845T DE69631845T2 (de) 1996-11-01 1996-11-01 Musterextrahierungsgerät
JP08536385A JP3035654B2 (ja) 1996-11-01 1996-11-01 パターン抽出装置
EP96935527A EP1011073B1 (en) 1996-11-01 1996-11-01 Pattern extraction apparatus
KR1019997003884A KR100295586B1 (ko) 1996-11-01 1999-04-30 패턴 추출장치

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP1996/003212 WO1998020454A1 (fr) 1996-11-01 1996-11-01 Appareil d'extraction de structure

Publications (1)

Publication Number Publication Date
WO1998020454A1 true WO1998020454A1 (fr) 1998-05-14

Family

ID=14154065

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP1996/003212 Ceased WO1998020454A1 (fr) 1996-11-01 1996-11-01 Appareil d'extraction de structure

Country Status (6)

Country Link
US (1) US6195460B1 (ja)
EP (1) EP1011073B1 (ja)
JP (1) JP3035654B2 (ja)
KR (1) KR100295586B1 (ja)
DE (1) DE69631845T2 (ja)
WO (1) WO1998020454A1 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3468182B2 (ja) 1999-12-21 2003-11-17 日本電気株式会社 画像マッチング装置および画像マッチング方法
WO2008023726A1 (fr) * 2006-08-22 2008-02-28 Yamatake Corporation Appareil de radar et procédé de mesure de la distance
JP2008089402A (ja) * 2006-10-02 2008-04-17 Konica Minolta Holdings Inc 情報処理システム、プログラムおよび情報処理方法
JP2010230578A (ja) * 2009-03-27 2010-10-14 Fujifilm Corp 偏芯量測定方法
KR101255952B1 (ko) 2011-07-21 2013-04-23 주식회사 앤비젼 패턴층이 형성된 기판의 간섭 현상을 이용한 패턴검사방법 및 패턴검사장치

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3846851B2 (ja) * 2001-02-01 2006-11-15 松下電器産業株式会社 画像のマッチング処理方法及びその装置
WO2005093654A2 (en) 2004-03-25 2005-10-06 Fatih Ozluturk Method and apparatus to correct digital image blur due to motion of subject or imaging device
US10721405B2 (en) 2004-03-25 2020-07-21 Clear Imaging Research, Llc Method and apparatus for implementing a digital graduated filter for an imaging apparatus
US9826159B2 (en) 2004-03-25 2017-11-21 Clear Imaging Research, Llc Method and apparatus for implementing a digital graduated filter for an imaging apparatus
WO2005103610A1 (ja) * 2004-04-22 2005-11-03 The University Of Electro-Communications 微小変位計測法及び装置
US20060034531A1 (en) * 2004-05-10 2006-02-16 Seiko Epson Corporation Block noise level evaluation method for compressed images and control method of imaging device utilizing the evaluation method
JP4525286B2 (ja) * 2004-10-14 2010-08-18 沖電気工業株式会社 生体情報認証装置および認証方法
JP4385139B2 (ja) 2006-02-01 2009-12-16 国立大学法人電気通信大学 変位検出方法、及び、変位検出装置、変位検出プログラム、並びに、位相特異点マッチング処理方法、位相特異点マッチング処理プログラム
JP5630957B2 (ja) * 2006-05-19 2014-11-26 日産化学工業株式会社 ハイパーブランチポリマー及びその製造方法
JP5178561B2 (ja) 2009-02-06 2013-04-10 Hoya株式会社 パターン検査方法、パターン検査装置、フォトマスク製造方法、およびパターン転写方法
CN102521590B (zh) * 2011-04-29 2013-12-04 北京大学 一种基于方向的左右掌纹识别方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02206294A (ja) * 1989-02-06 1990-08-16 Nippon Telegr & Teleph Corp <Ntt> 画像マッチング装置
JPH0816785A (ja) * 1994-06-27 1996-01-19 Nec Corp 画像認識装置および画像認識方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1570070A (en) * 1976-01-07 1980-06-25 Emi Ltd Arrangement for correlating pictures
US5619596A (en) * 1993-10-06 1997-04-08 Seiko Instruments Inc. Method and apparatus for optical pattern recognition
JPH07108526A (ja) 1993-10-12 1995-04-25 Nippon Steel Corp スライシングマシン
US5684890A (en) * 1994-02-28 1997-11-04 Nec Corporation Three-dimensional reference image segmenting method and apparatus
JPH07254062A (ja) 1994-03-14 1995-10-03 Yamatake Honeywell Co Ltd 指紋照合装置
US5915034A (en) * 1995-05-02 1999-06-22 Yamatake-Honeywell, Co., Ltd. Pattern collation apparatus based on spatial frequency characteristics
US5886744A (en) * 1995-09-08 1999-03-23 Intel Corporation Method and apparatus for filtering jitter from motion estimation video data
JPH09147115A (ja) * 1995-11-20 1997-06-06 Hamamatsu Photonics Kk 人物照合装置
US6067367A (en) * 1996-10-31 2000-05-23 Yamatake-Honeywell Co., Ltd. Moving direction measuring device and tracking apparatus

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02206294A (ja) * 1989-02-06 1990-08-16 Nippon Telegr & Teleph Corp <Ntt> 画像マッチング装置
JPH0816785A (ja) * 1994-06-27 1996-01-19 Nec Corp 画像認識装置および画像認識方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP1011073A4 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3468182B2 (ja) 1999-12-21 2003-11-17 日本電気株式会社 画像マッチング装置および画像マッチング方法
WO2008023726A1 (fr) * 2006-08-22 2008-02-28 Yamatake Corporation Appareil de radar et procédé de mesure de la distance
JP2008089402A (ja) * 2006-10-02 2008-04-17 Konica Minolta Holdings Inc 情報処理システム、プログラムおよび情報処理方法
JP2010230578A (ja) * 2009-03-27 2010-10-14 Fujifilm Corp 偏芯量測定方法
KR101255952B1 (ko) 2011-07-21 2013-04-23 주식회사 앤비젼 패턴층이 형성된 기판의 간섭 현상을 이용한 패턴검사방법 및 패턴검사장치

Also Published As

Publication number Publication date
DE69631845T2 (de) 2005-01-05
EP1011073A1 (en) 2000-06-21
JP3035654B2 (ja) 2000-04-24
DE69631845D1 (de) 2004-04-15
KR20000053002A (ko) 2000-08-25
EP1011073B1 (en) 2004-03-10
EP1011073A4 (en) 2000-06-21
KR100295586B1 (ko) 2001-08-07
US6195460B1 (en) 2001-02-27

Similar Documents

Publication Publication Date Title
WO1998020454A1 (fr) Appareil d&#39;extraction de structure
CN101061489B (zh) 用来处理视频数据的装置和方法
JPWO1998020454A1 (ja) パターン抽出装置
Uddin et al. Deepfake face detection via multi-level discrete wavelet transform and vision transformer
Cox Template matching and measures of match in image processing
JP2021120814A (ja) 学習プログラム、学習方法、及び情報処理装置
Wong et al. Local statistics for generative image detection
WO2019123917A1 (ja) 画像照合装置
JPH08248045A (ja) 速度測定方法および速度測定装置
Chen et al. Stereoscopic ranging by matching image modulations
CN119700062A (zh) 一种非接触心率检测方法、系统及相关设备
CN112807000A (zh) 鲁棒性脑电信号的生成方法及装置
US12425792B2 (en) Video processing device and method
Elbouz et al. Correlation based efficient face recognition and color change detection
JP2008225970A (ja) 画像合成装置およびそれを用いた画像照合装置ならびに画像合成方法
Popescu et al. VIBROTOOL—Software tool for change detection and diagnosis in vibration signals
JPH10134196A (ja) 距離測定装置
Bhowmik Improved single-sensor-based modal identification using singular spectrum analysis
Sitanggang et al. Image Edge Detection for Batak Ulos Motif Recognition using Canny Operators
JPH11296679A (ja) 画像特徴抽出方法および装置並びに画像特徴抽出プロ グラムを記録したコンピュータ読み取り可能な記録媒 体
Deravi Biometrics standards
JP3197193B2 (ja) 2次元信号フィルタリング装置
Kumar et al. Image disparity estimation based on fractional dual-tree complex wavelet transform: A multi-scale approach
CN121837127A (zh) 一种基于多尺度边缘特征融合的骨折智能检测方法
JP3058951B2 (ja) 画像パターン認識装置

Legal Events

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

Ref document number: 96180491.2

Country of ref document: CN

AK Designated states

Kind code of ref document: A1

Designated state(s): CN JP KR US

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE CH DE DK ES FI FR GB GR IE IT LU MC NL PT SE

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

Ref document number: 09297537

Country of ref document: US

Ref document number: 1019997003884

Country of ref document: KR

WWE Wipo information: entry into national phase

Ref document number: 1996935527

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 1996935527

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 1019997003884

Country of ref document: KR

WWG Wipo information: grant in national office

Ref document number: 1019997003884

Country of ref document: KR

WWG Wipo information: grant in national office

Ref document number: 1996935527

Country of ref document: EP