JPS61211783A - Image signal processing system - Google Patents
Image signal processing systemInfo
- Publication number
- JPS61211783A JPS61211783A JP60052921A JP5292185A JPS61211783A JP S61211783 A JPS61211783 A JP S61211783A JP 60052921 A JP60052921 A JP 60052921A JP 5292185 A JP5292185 A JP 5292185A JP S61211783 A JPS61211783 A JP S61211783A
- Authority
- JP
- Japan
- Prior art keywords
- pattern
- image
- processing
- spatial frequency
- power spectrum
- 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.)
- Pending
Links
Landscapes
- Character Discrimination (AREA)
Abstract
Description
【発明の詳細な説明】
(産業上の利用分野)
本発明は、OCR等のパタン認識装置に適用する画像信
号処理方式に関する。DETAILED DESCRIPTION OF THE INVENTION (Field of Industrial Application) The present invention relates to an image signal processing method applied to a pattern recognition device such as OCR.
(従来技術とその問題点)
OCR等のパタン認識装置の認識方式は、大別してアナ
ログ的マツチング法(テンプレートマツチング法)及び
構造解析法とに分ける事ができる。(Prior art and its problems) The recognition methods of pattern recognition devices such as OCR can be broadly divided into analog matching methods (template matching methods) and structural analysis methods.
前者は、標準パタンと未知パタンとの重ね合わせによっ
てパタンを認識する方式であり、汎用性に優れ、特に活
字の読み取シには、優利であるとされている。The former method recognizes a pattern by superimposing a standard pattern and an unknown pattern, and is said to be highly versatile and particularly useful for reading printed characters.
しかしこの方式には、重ね合わせの際に位置ずれがあっ
たり、回転、伸縮などのパタンの変形があると、標準パ
タンとの重なりあいが少くなり、正しく認識されないと
いう欠点があった。そこで従来は、位置ずれに対しては
、標準のパタンに人為的にぼかしを加える、或いは未知
のパタンを上下左右に移動させて、最も重なりの大きい
位置を探す等の処理を行なって認識率を向上させていた
。However, this method has the drawback that if there is a positional shift during superimposition or if the pattern is deformed such as rotation or expansion/contraction, the overlap with the standard pattern will be reduced and the pattern will not be recognized correctly. Conventionally, to deal with positional deviations, the recognition rate was improved by artificially blurring the standard pattern, or by moving the unknown pattern up, down, left and right, and searching for the position with the greatest overlap. It was improving.
しかしながら、このような方式では、未知パタンの位置
の微小な変位については、認識率の改善がなされるもの
の、大きく位置がずれているパタンの認識には、効果が
なく、そのためさらに位置合わせ等の前処理が必要とな
る欠点があった。また、人為的にぼかしを加える方法で
は、空間分解能が落ちる欠点がある。However, although this method improves the recognition rate for small displacements in the position of unknown patterns, it is ineffective for recognizing patterns with large displacements, and therefore requires further alignment, etc. There was a drawback that pretreatment was required. Furthermore, the method of artificially adding blur has the disadvantage that spatial resolution decreases.
“プロシーデイングズオブジアイイーイーイー。“The Proceedings of the Ieeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee.
1970.58巻&l(“ProceedinJs o
f theよりEE。1970.58 Volume & l (“ProceedinJs o
EE from f the.
197(L vol s 8″)に掲載のジー、ジー、
レンダリス゛(’G、 G、 Lenaariθ)ら″
による文献“デイフラクションーパタンサンプリングフ
オーオートマテイツクパタンリコグニション″(“Di
ffraction −Patl;ern Sampl
inJ for Automatic Pattern
RecoJlnition″)には、ホログラフィ−を
利用してパタンのフランホーファー回折像を作り、との
回折像によってパタン認識を行なう方法が述べられてい
る。この方法は、パタンに2次元フーリエ変換を施して
、パタンをその空間周波数のパワースペクトルへと変換
する事に相当し、もとのパタンを並進移動させても、そ
のパワースペクトルは不変であるために、パタンの位置
ずれに対して強い方式である。しかしながら、この方式
では、呈示されるパタンか光を透過するフィルム等であ
る必要があるために、実用化には難点があった。G, G, published in 197 (L vol s 8″)
Renderis゛('G, G, Lenaariθ) et al.''
The document “Diffraction-Pattern Sampling Flow Automation Pattern Recognition” (“Di
ffraction-Patl;ern Sample
inJ for Automatic Pattern
RecoJlnition'') describes a method in which a Franhofer diffraction image of a pattern is created using holography, and the pattern is recognized from the diffraction image.This method involves performing two-dimensional Fourier transformation on the pattern. , which corresponds to converting a pattern into a power spectrum of its spatial frequency, and since the power spectrum remains unchanged even if the original pattern is translated, it is a method that is robust against pattern displacement. However, this method has been difficult to put into practical use because the pattern to be presented needs to be a film or the like that transmits light.
(発明の目的)
本発明の目的は、アナログ的マツチング法に関するこの
ような従来の欠点を取り除き、パタンの位置ずれに対し
て強いパタン認識のだめの高速の画像信号処理方式を提
供することにある。(Object of the Invention) An object of the present invention is to eliminate such conventional drawbacks of the analog matching method and to provide a high-speed image signal processing method for pattern recognition that is resistant to pattern displacement.
(′発明の構成)
パタンの位置ずれに対して強くするという要求は、入力
画像を互いにオーバーラツプする又はオーバーラツプし
ない小領域に分割し、各小領域ごとに入力画像をその空
間周波数パワースペクトル又は自己相関関数のパタンへ
変換する処理を行なう事により満たされる。とうして得
られたパタンのセットと同様の処理を予め施しておいだ
標準パタンとてアナログ的マツチング法により重なり合
いを調べれば、パタ/の位置ずれに対して強い高速なパ
タン認識が実現できる。('Structure of the Invention) The requirement to be strong against pattern displacement is to divide the input image into small regions that overlap or do not overlap each other, and to analyze the spatial frequency power spectrum or autocorrelation of the input image for each small region. This is satisfied by converting it into a function pattern. By checking the overlap using an analog matching method using a standard pattern that has been previously subjected to the same processing as the set of patterns obtained, high-speed pattern recognition that is resistant to positional deviations of patterns can be realized.
(構成の詳細な説明)
第1図は本発明方式を利用した、パタン認識方式の流れ
図である。テレビカメラもしくはイメージスキャナーか
ら入力されたパタンの画像情報は、画素ごとに明暗度が
数値化されて、画像メモリに収納される。本方式による
処理では、この画像情報に対してプロセッサー等を用い
て以下のような変換を行なう。(Detailed Description of Configuration) FIG. 1 is a flowchart of a pattern recognition method using the method of the present invention. Pattern image information input from a television camera or image scanner is converted into numerical values for the brightness of each pixel and stored in an image memory. In the processing according to this method, the following conversion is performed on this image information using a processor or the like.
今、入力画像の画素数をNXNとする。まず、この入力
画像を、画素数n、 xn、からなる、n2×n2コの
小領域に分割する。ここで互いにオーバーラツプがない
ように画像を分割する場合には、N−n。Now, assume that the number of pixels of the input image is NXN. First, this input image is divided into n2×n2 small regions each having n and xn pixels. If the images are to be divided so that there is no overlap, N-n.
xn、であるが、互いにオーバーラツプがあるように分
割してもかまわない。次に、これらの小領域に分割され
た各画像に対して、適当な処理を施して呈示画像を並進
移動させても、変換されたパタンか不変になるようにす
る。例えば、画像の空間周波数パワースペクトルは画像
を並進移動させても不変であるから、空間周波数パワー
スペクトルを求めるような処理を施せば良い。あるいは
画像の自己相関関数もまた画像の並進移動に関して不変
であるから自己相関関数を求めるような処理でもかまわ
ない。xn, but they may be divided so that they overlap each other. Next, each image divided into these small regions is subjected to appropriate processing so that even if the presented image is translated, the transformed pattern remains unchanged. For example, since the spatial frequency power spectrum of an image remains unchanged even when the image is translated, processing to obtain the spatial frequency power spectrum may be performed. Alternatively, since the autocorrelation function of the image also remains unchanged with respect to the translational movement of the image, a process of determining the autocorrelation function may be used.
ココテ各小領域を添字(X、 Y ) (X、 Y =
l 〜n、 )で表わし、各小領域中での画素の位置
を添字(X。Subscript each small area (X, Y) (X, Y =
The position of a pixel in each small area is represented by a subscript (X).
7)(X、7=1〜n、)で表わす。また(X、Y)で
表わされる小領域中の(x、y)の位置にある画素の値
を工(X、3”、XIY)で表わす1.すると、上に述
べた空間周波数パワースペクトルの・母タンは、次式に
より計算される。7) Represented by (X, 7=1 to n,). Also, the value of the pixel at the position (x, y) in the small area represented by (X, Y) is expressed as (X, 3", XIY)1. Then, the spatial frequency power spectrum described above・Mother tongue is calculated by the following formula.
・・・・・・・・・ (1)
P (R,、R2; X、 Y )−F(R11R2:
X、’ Y) ” 1.、(2)(1)式ハ、I (x
、7 ; xt Y) の添字(x+y)に関する2
次元の離散値ツーIJ 工変換を表わす。(1)式の計
算に関しては、高速フーリエ変換のマルゴリズムを利用
する事が可能であり、その場合には、処理時間がより短
縮され得る。(2)式のPが、空間周波数・9ワースベ
クトルである〇
また、先述の処理に於いて、自己相関関数を求める場合
には、これは次式により計算される07+d、 :X、
Y)
もしも、今季領域の境界での効果が無視できるとすると
、PもしくはCは、入力画像・母タンの並進移動に対し
て不変なノぞタンになるから、このような処理によって
入力画像の位置ずれに関して非常に安定なパタンか得ら
れる事になる(この条件は入力画像中の各パタンの大き
さが、小領域の大きさく nt Xns )に比べて充
分小さい場合に満足される(第2図))。従って、この
ような処理を施した・母タンを次のパタン判定ユニット
へ送って、あらかじめ同様の処理を施しておいた標準パ
タンとのテンプレートマツチングを行なえば、パタンの
(n、程度以下の大きさの)位置ずれに対して強い・母
タン認識方式が実現される事になる。・・・・・・・・・ (1) P (R,, R2; X, Y)-F(R11R2:
X, ' Y) ” 1., (2) (1) Formula c, I (x
, 7 ; xt Y) 2 regarding the subscript (x+y)
Represents the discrete value-to-IJ transformation of the dimension. Regarding the calculation of equation (1), it is possible to use fast Fourier transform algorithm, and in that case, the processing time can be further reduced. P in equation (2) is the spatial frequency/9worth vector. Also, in the above process, when calculating the autocorrelation function, it is calculated by the following equation: 07+d, :X,
Y) If the effect at the boundary of the current season region can be ignored, P or C will be a nozotan that remains unchanged with respect to the translational movement of the input image/mother tongue, so by such processing, the input image A very stable pattern with respect to positional shift can be obtained (this condition is satisfied when the size of each pattern in the input image is sufficiently small compared to the size of the small area nt figure)). Therefore, by sending the mother tongue that has undergone such processing to the next pattern judgment unit and performing template matching with a standard pattern that has undergone similar processing in advance, it is possible to A mother tongue recognition method that is resistant to positional deviations (in terms of size) will be realized.
入力画は中のパタンか小領域の大きさに比べである程度
大きい場合には、各小領域での処理に於いて境界での効
果が無視できなくなるために、P(R11R29XI
Y)あるイハC(d+ s d2 :X、Y) は厳
密には、入力・9タンの並進移動に関して不変ではなく
なる。しかしながら、この場合にもやはり、Pあるいは
Cは、入力画像そのものと比較すると、並進移動に関し
て安定なパタンになっている。If the input image is relatively large compared to the size of the inner pattern or small area, the effect at the boundary cannot be ignored in the processing of each small area, so P(R11R29XI
Y) Strictly speaking, a certain IhaC(d+s d2 :X, Y) is no longer invariant with respect to the translational movement of the input/9 tongue. However, in this case as well, P or C has a stable pattern with respect to translational movement when compared with the input image itself.
第3図を用いて、その理由を説明する(本発明が主に対
象とするパタンは2次元のパタンであるが、ここでは、
簡単のため1次元の信号・9タンを用いて説明する)。The reason for this will be explained using FIG. 3 (the patterns that the present invention mainly targets are two-dimensional patterns;
For simplicity, we will explain using a one-dimensional signal, 9 tangs).
良く知られているように、1悠激な変化のある信号に対
して、ある幅、nl の範囲で移動平均をとると、信
号の高周波成分が取り除かれなだらかな変化をするよう
になる。このように適当な重み関数を用いて・母タンの
移動平均をとってからテングレートマツチングを行なう
方式は・ザタンを人為的にぼかす方法に対応しており、
位置ずれに対して強い方式になる。As is well known, when a moving average is taken over a certain width (nl) for a signal that changes drastically, the high frequency components of the signal are removed and the signal changes gradually. This method of using an appropriate weighting function to take the moving average of the mother tongue and then performing ten rate matching corresponds to the method of artificially blurring the mother tongue.
This method is resistant to misalignment.
このような移動平均を用いる方式と本発明による処理方
式との最大の違いは、前者が、・9タンの位置変動に強
くするため、パタンの高周波成分を捨てているのに対し
て、後者に於いては・9タンの高周波成分の情報を捨て
る事なく、・9タンの位置変動に関する安定性を実現し
ている点にある。これは、前者がn、の幅の入力信号か
ら単に、その平均というスカラー量を計算するだけであ
るのに対して、後者では、同じ入力信号から、その空間
周波数パワースペクトル(もしくは、自己相関関数)の
パタンというベクトル量を計算しているからである。従
って、前者では許容できる位置ずれの大きさは、分解能
とトレードオフの関係にあるから、許容できる位置ずれ
の大きさは分解能により制限を受ける事になるが、後者
に於いては、これらの二つの量は、独立であるから、許
容できる位置ずれの大きさは自由に設定できる事になる
。The biggest difference between the method using such a moving average and the processing method according to the present invention is that the former discards high frequency components of the pattern in order to make it resistant to positional fluctuations, whereas the latter discards the high frequency components of the pattern. In this case, ・Stability regarding the positional fluctuation of the 9-tongue is achieved without discarding the information on the high-frequency components of the 9-tongue. The former simply calculates a scalar quantity called the average from an input signal with a width of n, whereas the latter calculates its spatial frequency power spectrum (or autocorrelation function) from the same input signal. This is because we are calculating a vector quantity called the pattern of ). Therefore, in the former case, the allowable size of positional deviation is in a trade-off relationship with resolution, so the allowable size of positional deviation is limited by the resolution, but in the latter case, these two Since the two quantities are independent, the allowable magnitude of positional deviation can be set freely.
レンダリス(Lendari日)らがホログラフィ−を
用いて実現したように、入力画像全体の空間周波数パワ
ースペクトルを求める処理をプロセッサー等を用いて実
現する方法によっても、本発明による方式と同様に、画
像の高周波成分な落とする事なく、・臂タンの位置変動
に対して強い・ヤタン認識が実現できると考えられるが
、その場合には、計算量が膨大となるので実用的とはe
えない〇今 N2この画素からなる画1象の空間周波数
パワースペクトル又は自己相関関数を求めるのに、N
この計算(かけ算)が必要だとする。本発明のように、
小頭域に分割して処理を行なうと、111Xl’l、、
ΣNxn1 この計算量で済み、計算量の減少は、N
X n、 / N αl / n;となる。高速フーリ
エ変換を用いると、この評側よりも、計算量の減少は小
さくなるが、やはり大幅な減少が実現できる。従って本
発明による方式を用いる事により、パタンの位置変動に
強い、高速のパタン認識が実現できる。As was achieved by Lendari et al. using holography, a method of using a processor or the like to obtain the spatial frequency power spectrum of the entire input image can also be used to obtain image data in the same way as the method according to the present invention. It is conceivable that it is possible to achieve yataan recognition that is strong against changes in the position of the armpit tongue without dropping high-frequency components, but in that case, the amount of calculation would be enormous, so it is not practical.
No 〇 Now N2 To find the spatial frequency power spectrum or autocorrelation function of this image consisting of pixels, N
Suppose we need this calculation (multiplication). As in the present invention,
When processing is performed by dividing into small head regions, 111Xl'l,...
ΣNxn1 This amount of calculation is enough, and the reduction in calculation amount is N
X n, / N αl / n; If fast Fourier transform is used, the reduction in calculation amount will be smaller than on this evaluation side, but a significant reduction can still be achieved. Therefore, by using the method according to the present invention, high-speed pattern recognition that is resistant to pattern position fluctuations can be realized.
(実施例)
第4図は、本発明による処理をマルチ・グロセツサー構
成による並列処理方式により実現した実施例を示す流れ
図である。テレビカメラもしくはイメージスキャナーか
ら入力されたパタンの画像情報は画素ごとに明暗団が数
値化されて、画1象メモリに収納される。次に入力画1
M (lI!Ii素数なN×Nとする)を、画素数n1
Xn、からなる”2 xntこの小頭域に分割する。こ
れらの小頭域に分割された各画1家は、それぞれ専用の
プロセッサーへと送られ、各小頭域ごとに、その空間周
波数・9ワースペクトルの・9タンへと変換される。こ
の時、その処理に、高速フーリエ変換を用いても良い。(Embodiment) FIG. 4 is a flow chart showing an embodiment in which the processing according to the present invention is realized by a parallel processing method using a multi-grocer configuration. Image information of a pattern input from a television camera or an image scanner is converted into numerical values of chiaroscuro for each pixel and stored in the image memory. Next, input image 1
M (Let lI!Ii be a prime number N×N), the number of pixels n1
Xn, is divided into "2xnt" sub-regions. Each stroke divided into these sub-regions is sent to a dedicated processor, and for each sub-region, its spatial frequency, It is converted into 9-tan of the 9-war spectrum.At this time, fast Fourier transform may be used for the processing.
このようにして得られた、小領域ごとの空間周波数7f
ワースベクトルのノ4タンのセットは、次のパタン判定
ユニットへと送られ、そこに於いて、予め同様の処理を
施しておいて標準パタンを用いてテンプレートマツチン
グによシ、パタン判定が行なわれる0本発明による処理
は、分割された小領域ごとに、独立に、実行できるので
、第4図のような、並列処理化が可能になるのである。Spatial frequency 7f for each small area obtained in this way
The set of four patterns of worth vectors is sent to the next pattern judgment unit, where the same processing is performed in advance and pattern judgment is performed by template matching using a standard pattern. Since the processing according to the present invention can be executed independently for each divided small area, parallel processing as shown in FIG. 4 is possible.
勿論、処理を並列化することで、処理スピードは増加す
る。Of course, processing speed increases by parallelizing processing.
(発明の効果)
本方式による処理を行なってノfタンを変換しておいて
から、ノ量タンマツチングを行なうと、呈示・母タンの
位置ずれに強く、かつ高速に実行できる・ぐタン認識方
式が実現できる。従来のような・fタンの移動平均を取
って、・9タンに人為的に、ぼかしを加える方式と比べ
ると本方式は呈示・臂タンの高周波成分の情報を捨てる
事なく、位置ずれに対して強くできるので分解能を落と
す事なく、太きな位置ずれにも対応できるという特長を
有する。(Effects of the Invention) By performing the processing using this method to convert the nof tongue and then performing the no amount tongue matching, the tongue recognition method is resistant to misalignment of the presentation and mother tongue and can be executed at high speed. can be realized. Compared to the conventional method of taking the moving average of the f tongue and artificially adding blur to the nine tongue, this method does not discard information on the high frequency components of the presented arm tongue, and is effective against positional shifts. It has the advantage of being able to handle large displacements without compromising resolution.
また、入力画像全体を、その空間周波数パワースペクト
ルのパタンへと変換する方式と比べて、本発明の方式で
は、計算量の大幅な削減が実現でき、高速な処理が可能
である。まだ、処理が分割された、各小領域ごとに独立
に実行できるので、実施例に示したように、マルチプロ
セッサー化して高速化を計るのに適した方式もできる。Furthermore, compared to a method that converts the entire input image into a pattern of its spatial frequency power spectrum, the method of the present invention can achieve a significant reduction in the amount of calculation and can perform high-speed processing. However, since processing can be executed independently for each divided small area, it is possible to use a method suitable for increasing speed by using multiple processors as shown in the embodiment.
第1図は本発明による画像信号処理を用いた・臂タン認
識方式の流れ図、第2図は本発明による処理で並進移動
によっても変換された・9タンが全く変化しないような
zfJンの例を示す図、第3図は移動平均操作と本発明
による処理との比較を示す図、第4図はマルチプロセッ
サー構成による本発明の一実施例を示す流れ図である。
代理人 弁理士 本 庄 伸 介
第3図
X+rh
十多動’If−少勺: I(x) j(x)= (、
xI(x)dx(昭二二I(X) P(R,X)
第4図Figure 1 is a flowchart of the arm tongue recognition method using image signal processing according to the present invention, and Figure 2 is an example of a zfJn that is converted by translational movement using the processing according to the present invention and does not change at all. FIG. 3 is a diagram showing a comparison between the moving average operation and the processing according to the present invention, and FIG. 4 is a flowchart showing an embodiment of the present invention using a multiprocessor configuration. Agent Patent Attorney Shinsuke Honjo Figure 3
xI(x)dx(Sho22 I(X) P(R,X) Fig. 4
Claims (1)
プしない小領域に分割し、これら各小領域ごとに、前記
入力画像を空間周波数パワースペクトル又は自己相関関
数のパタンへ変換する事を特徴とする画像信号処理方式
。An image signal processing method characterized by dividing an input image into small regions that overlap or do not overlap each other, and converting the input image into a pattern of a spatial frequency power spectrum or an autocorrelation function for each of these small regions. .
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP60052921A JPS61211783A (en) | 1985-03-16 | 1985-03-16 | Image signal processing system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP60052921A JPS61211783A (en) | 1985-03-16 | 1985-03-16 | Image signal processing system |
Publications (1)
Publication Number | Publication Date |
---|---|
JPS61211783A true JPS61211783A (en) | 1986-09-19 |
Family
ID=12928293
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP60052921A Pending JPS61211783A (en) | 1985-03-16 | 1985-03-16 | Image signal processing system |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPS61211783A (en) |
-
1985
- 1985-03-16 JP JP60052921A patent/JPS61211783A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Giachetti | Matching techniques to compute image motion | |
US5311305A (en) | Technique for edge/corner detection/tracking in image frames | |
US20060159311A1 (en) | Dominant motion analysis | |
JP4527885B2 (en) | Rotation correction and duplicate image identification by Fourier transform correlation | |
US20140126839A1 (en) | Defect detection using joint alignment and defect extraction | |
US11244475B2 (en) | Determining a pose of an object in the surroundings of the object by means of multi-task learning | |
JP2010086267A (en) | Image processing apparatus for vehicle | |
CN116433728A (en) | DeepSORT target tracking method for shake blur scene | |
US4985765A (en) | Method and apparatus for picture motion measurement whereby two pictures are correlated as a function of selective displacement | |
JPS61211783A (en) | Image signal processing system | |
RU2718211C1 (en) | Device for combining infrared images | |
JPS63211474A (en) | Hierarchization structural template matching method | |
Techmer | Real-time motion based vehicle segmentation in traffic lanes | |
CN113921412A (en) | Method, device and equipment for calculating chip period in wafer | |
EP0366165A2 (en) | Picture motion measurement | |
JP3074292B2 (en) | Template matching processing method | |
JPS61241885A (en) | Recognizing method for picture signal | |
CN104318580A (en) | Pattern search method and device | |
JPS62205488A (en) | Recognizing method for image signal | |
JP3331275B2 (en) | Projection parameter extraction processing method | |
JPH05225335A (en) | Object area segmenting device | |
JPH08279048A (en) | Image search device | |
Shivakumara et al. | Fourier transform based column-block and row-block matching procedure for document image mosaicing | |
JP2918667B2 (en) | Mobile vehicle environment recognition device | |
JP4037503B2 (en) | Image feature extraction method and apparatus |