JPH0527152B2 - - Google Patents

Info

Publication number
JPH0527152B2
JPH0527152B2 JP59013371A JP1337184A JPH0527152B2 JP H0527152 B2 JPH0527152 B2 JP H0527152B2 JP 59013371 A JP59013371 A JP 59013371A JP 1337184 A JP1337184 A JP 1337184A JP H0527152 B2 JPH0527152 B2 JP H0527152B2
Authority
JP
Japan
Prior art keywords
local
area
local area
absolute value
position correction
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.)
Expired - Lifetime
Application number
JP59013371A
Other languages
Japanese (ja)
Other versions
JPS60159972A (en
Inventor
Mutsumi Watanabe
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.)
Toshiba Corp
Original Assignee
Tokyo Shibaura Electric Co Ltd
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 Tokyo Shibaura Electric Co Ltd filed Critical Tokyo Shibaura Electric Co Ltd
Priority to JP59013371A priority Critical patent/JPS60159972A/en
Publication of JPS60159972A publication Critical patent/JPS60159972A/en
Publication of JPH0527152B2 publication Critical patent/JPH0527152B2/ja
Granted legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/80Geometric correction

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Description

【発明の詳細な説明】 [発明の技術分野] 本発明は、相関係数と2枚の画像の画素毎の濃
度差という2つの指標を用いて、2枚の画像間の
位置ズレを高密度に補正できる、実用性の高い画
像の精密位置合わせ装置に関する。
Detailed Description of the Invention [Technical Field of the Invention] The present invention uses two indicators, a correlation coefficient and a density difference for each pixel of two images, to calculate the positional deviation between two images with high density. The present invention relates to a highly practical precision image positioning device that can perform image correction.

[発明の技術的背景] 衛星画像処理、或いは医用画像処理などの分野
において、時間差をもつて撮像入力された複数枚
の画像より新たに1枚の画像を合成しようとする
場合、各画像間の位置合わせが必要となる。この
位置ズレの特性が既知であり、画像全体について
一様とみなせ、かつ平行移動・回転・拡大・縮小
等の1次の座標変換式で表わされる場合には、比
較的簡単な信号処理で補正可能である。
[Technical Background of the Invention] In fields such as satellite image processing or medical image processing, when attempting to synthesize a new image from multiple images captured and input with a time difference, it is necessary to Positioning is required. If the characteristics of this positional shift are known, can be considered uniform for the entire image, and can be expressed by a linear coordinate transformation formula such as parallel translation, rotation, enlargement, or reduction, it can be corrected using relatively simple signal processing. It is possible.

[背景技術の問題点] だが、前記衛星画像、或いは医用画像の場合に
は位置ズレの原因は幾つかの要素が複合したもの
である場合が多く、位置ズレの量が局所的に変化
するため、単一の1次座標変換で画像全体の位置
合わせを行なうのは困難である。
[Problems with the background art] However, in the case of the above-mentioned satellite images or medical images, the cause of positional deviation is often a combination of several factors, and the amount of positional deviation changes locally. , it is difficult to align the entire image with a single linear coordinate transformation.

このため、従来は位置を合わせるべき画像を幾
つかの小領域に分割し、この各修復対象領域部分
を移動させつつ、合わせられるべき画像との相関
係数を計算し、この値が最大になつた部分を位置
の合つたところとしていた。
For this reason, conventionally, the image to be aligned is divided into several small areas, and each of these areas to be repaired is moved while calculating the correlation coefficient with the image to be aligned, and this value is maximized. The correct position was the part where the

しかしこの方法では、位置ズレが微小、或いは
画像のコントラスト差が小さい場合には相関係数
値の変化量が非常に小さくなり感度が悪く、画像
にノイズが重畳した場合などこの変化量が隠れて
しまい誤つた位置に収束する場合がある、などの
問題点があつたため高精度の位置合わせは困難で
あり、実用性に乏しかつた。
However, with this method, when the positional shift is minute or the contrast difference between images is small, the amount of change in the correlation coefficient becomes very small and the sensitivity is poor, and when noise is superimposed on the image, this amount of change may be hidden. Because of problems such as convergence sometimes converging on the wrong position, highly accurate alignment was difficult and impractical.

[発明の目的] 本発明はこのような問題点を克服するよう考慮
してなされたもので、その目的とするところは、
相関係数と2枚の画像の画素毎の濃度差という2
つの指標を用いて2枚の画像間の位置ズレを高精
度に補正する画像の精密位置合わせ装置を提供す
ることにある。
[Object of the invention] The present invention has been made in consideration of overcoming such problems, and its purpose is to:
2. Correlation coefficient and pixel-by-pixel density difference between two images
An object of the present invention is to provide a precision image positioning device that uses two indices to correct positional deviation between two images with high precision.

[発明の概要] 第1図に本発明の原理図を示す。[Summary of the invention] FIG. 1 shows a diagram of the principle of the present invention.

本発明は、2枚の画像をそれぞれ位置ズレ量が
均一と考えられる複数個の局所領域にP11,P12
P21,…,Q11,Q12,Q21,…分割した後、これら
の各局所領域毎にまず近傍領域内における相関係
数を各々求め、この値が最大となる位置に画像面
内で移動を繰り返す粗位置補正を行なつた後、前
記近傍領域内における相関係数の変化の最大値が
既め与えられた閾値以下になつた時点で今度は前
記近傍領域における2枚の画像の画素毎の濃度差
の評価関数が最大(最小)となる方向に移動を反
復して、精密位置補正を実行するものである。上
記移動方向を最大値方向とするか最小値方向とす
るかは上記評価関数の作り方に依り決定される。
In the present invention, two images are divided into a plurality of local regions where the amount of positional deviation is considered to be uniform, P 11 , P 12 ,
After dividing P 21 , ..., Q 11 , Q 12 , Q 21 , ..., first calculate the correlation coefficient in the neighboring region for each local region, and then set the correlation coefficient in the image plane at the position where this value is maximum. After performing coarse position correction by repeating movement, when the maximum value of the change in the correlation coefficient in the neighboring area becomes equal to or less than a predetermined threshold, the pixels of the two images in the neighboring area are Precise position correction is executed by repeatedly moving in the direction in which the evaluation function of each density difference is maximized (minimum). Whether the moving direction is the maximum value direction or the minimum value direction is determined depending on how the evaluation function is created.

最初は位置ズレが大きくてコントラスト差が無
視できないため、正規化された指標である相関係
数値を用いて上記粗位置補正を行ない、位置がほ
とんど合つた後はコントラストがほとんど同一と
考えられるため、2枚の画像の画素毎の濃度差を
も考慮した上記評価関数を指標とした上記精密位
置補正を実行することにより、相関係数値のみを
用いた場合に比べて前記位置ズレが微小、コント
ラスト差が小、或いは多少異なる画像の位置合わ
せの場合にも誤つた位置への収束防止が期待で
き、その実用的効果は多大である。
Initially, the positional deviation is large and the contrast difference cannot be ignored, so the above coarse positional correction is performed using the correlation coefficient, which is a normalized index, and after the positions are almost aligned, the contrast is considered to be almost the same. By executing the above-mentioned precise position correction using the above-mentioned evaluation function as an index, which also takes into account the density difference for each pixel of the two images, the positional deviation is minute and the contrast difference is smaller than when only the correlation coefficient is used. Even in the case of alignment of images with small or slightly different values, it can be expected to prevent convergence to an erroneous position, and this has a great practical effect.

[発明の実施例] 以下、本発明の一実施例につき図面を参照して
説明する。
[Embodiment of the Invention] An embodiment of the present invention will be described below with reference to the drawings.

第2図に全体の構成図を示す。 Figure 2 shows the overall configuration.

撮像系1より画像メモリ2,3に入力された2
枚の画像の位置合わせを行なう場合、まず局所領
域分割回路4により各々複数個の局所領域に分割
した後位置補正対象領域指定回路5によつて位置
合わせを行なう局所領域対を選択し、局所領域メ
モリ6,7に格納する。
2 input from the imaging system 1 to the image memories 2 and 3
When aligning two images, the local area dividing circuit 4 first divides each into a plurality of local areas, and then the position correction target area specifying circuit 5 selects a pair of local areas to be aligned. Store in memories 6 and 7.

まず粗位置補正部100では、上記局所領域メ
モリ6及び近傍領域切り出し回路8により切り出
された上記局所領域メモリ7内の両画像に対して
相関係数計算回路9により相関係数値を求め相関
係数値メモリ10に格納する。ここで位置合わせ
すべき一方の局所領域Aの画素毎の濃度をaij(i
=1〜m,j=1〜n)、その平均値をμAとし、
他方の局所領域Bの画素毎の濃度をbij(i=1〜
m,j=1〜n)、その平均値をμBとすると、こ
れら局所領域対の相関係数CABは以下の式で表さ
れる。
First, in the coarse position correction section 100, a correlation coefficient is calculated by a correlation coefficient calculating circuit 9 for both images in the local area memory 7 that have been cut out by the local area memory 6 and the neighboring area cutting circuit 8, and a correlation coefficient is calculated. It is stored in the memory 10. Here, the density of each pixel of one local area A to be aligned is a ij (i
= 1 to m, j = 1 to n), the average value is μ A ,
The density of each pixel in the other local area B is expressed as b ij (i=1~
m, j=1 to n), and its average value is μ B , then the correlation coefficient C AB of these local area pairs is expressed by the following formula.

この操作を前記近傍領域切り出し回路8により指
定される全近傍領域について行なつた後、これら
の相関係数値群のうち最大のものを最大値検出回
路11により求め、移動位置決定回路12により
この相関係数の最大となる方向に画像内で位置移
動を行なう。同時に最大変化量計算回路13によ
り、前記近傍領域内の相関係数値の最大変化量と
予め与えられた変化量閾値14とを比較器15に
より比較し、上記最大変化量が上記閾値以下にな
る(つまり相関係数の変化がほとんど無くなる)
迄前記粗位置補正をくり返す。
After performing this operation for all the neighboring regions specified by the neighboring region extraction circuit 8, the maximum value detection circuit 11 obtains the largest one among these correlation coefficient groups, and the movement position determining circuit 12 calculates this correlation value. The position within the image is moved in the direction that maximizes the relationship coefficient. At the same time, the maximum change amount calculation circuit 13 compares the maximum change amount of the correlation coefficient in the neighborhood region with a predetermined change amount threshold value 14 using a comparator 15, and the maximum change amount becomes equal to or less than the threshold value ( In other words, there is almost no change in the correlation coefficient)
The above coarse position correction is repeated until

閾値以下になつた場合は、粗/精密切り換え回
路16により精密位置補正部200にうつる。
If the value falls below the threshold, the coarse/fine switching circuit 16 transfers the information to the fine position correction section 200.

本実施例では評価関数として絶対値誤差を用い
ることとする。この場合移動方向は最小値の方向
となる。
In this embodiment, an absolute value error is used as the evaluation function. In this case, the direction of movement is the direction of the minimum value.

上記精密位置補正部200では、前記局所領域
メモリ6,7内の画像間の絶対値誤差を絶対値誤
差計算回路17により求め絶対値誤差メモリ18
に格納する。
In the precise position correction section 200, the absolute value error between the images in the local area memories 6 and 7 is calculated by the absolute value error calculation circuit 17, and the absolute value error memory 18
Store in.

この操作を前記近傍領域切り出し回路8により
指定される全近傍領域について行なつた後、これ
らの絶対値誤差のうち最小のものを最小値検出回
路19により求め移動位置決定回路20によりこ
の絶対値誤差の最小となる方向に画像内での位置
移動を行なう。
After performing this operation for all the neighboring regions specified by the neighboring region extraction circuit 8, the minimum value detection circuit 19 determines the smallest absolute value error among these absolute value errors, and the movement position determining circuit 20 calculates this absolute value error. The position within the image is moved in the direction that minimizes .

同時に上記近傍領域内における絶対値誤差の最
小値と予め与えられた絶対値誤差閾値21を比較
器22により比較し精密位置補正終了判定回路2
3により上記最小値が上記閾値以下になつた場合
は位置が合つたと判断して精密位置補正を終了
し、前記位置補正対象領域指定回路5により新し
い位置補正の対象となる局所領域対を指定し、再
び上記の操作を繰り返す。
At the same time, a comparator 22 compares the minimum value of the absolute value error in the above-mentioned neighborhood area with a predetermined absolute value error threshold 21, and the precision position correction completion determination circuit 2
If the minimum value becomes equal to or less than the threshold value in step 3, it is determined that the position has been matched, and the precise position correction is terminated, and the position correction target area designation circuit 5 designates a pair of local areas to be a new position correction target. and repeat the above operation again.

[発明の効果] 以上述べたように、まず粗位置補正部により相
関係数値を指標とする位置補正を行ない、この相
関係数値の変化がほとんど無くなつた時点で2枚
の画像の画素毎の濃度差の評価関数を指標とする
精密位置補正を行なうことによつて2枚の画像間
の簡易かつ高精度な位置合わせが実現でき、その
実用的効果は多大である。
[Effects of the Invention] As described above, first, the coarse position correction section performs position correction using the correlation coefficient as an index, and when the change in the correlation coefficient has almost disappeared, the position correction for each pixel of the two images is performed. By performing precise position correction using the density difference evaluation function as an index, it is possible to achieve simple and highly accurate position alignment between two images, which has great practical effects.

尚、本発明は上述した実施例に限定されるもの
ではない。例えば、前記実施例では精密位置補正
評価関数として絶対値誤差を用いたが、より一般
的に2枚の画像の画素毎の濃度差と相関係数とを
組み合わせた評価関数を採ることにより、より高
精度な位置補正が実現できる。
Note that the present invention is not limited to the embodiments described above. For example, in the above embodiment, the absolute value error was used as the precision position correction evaluation function, but more generally, by adopting an evaluation function that combines the density difference for each pixel of two images and the correlation coefficient, it is possible to Highly accurate position correction can be achieved.

又、粗位置補正部において、相関係数を指標と
する代りに、画像中の特徴点をそれぞれ抽出し、
これらの位置を合わせるfeature watchingの手
法を用いることにより、更に高速な位置合わせが
期待できる。
In addition, in the coarse position correction section, instead of using the correlation coefficient as an index, each feature point in the image is extracted,
By using a feature watching method to align these positions, even faster alignment can be expected.

要するに、本発明はその主旨を逸脱しない範囲
で種々変形して実施することができる。
In short, the present invention can be implemented with various modifications without departing from the spirit thereof.

【図面の簡単な説明】[Brief explanation of the drawing]

第1図は本発明の原理図、第2図は本発明の一
実施例の構成図である。 1……撮像系、2,3……画像メモリ、6,7
……局所領域メモリ、100……粗位置補正部、
200……精密位置補正部。
FIG. 1 is a principle diagram of the present invention, and FIG. 2 is a block diagram of an embodiment of the present invention. 1... Imaging system, 2, 3... Image memory, 6, 7
... Local area memory, 100 ... Coarse position correction section,
200...Precise position correction section.

Claims (1)

【特許請求の範囲】 1 2枚の画像を格納した画像メモリと、 前記2枚の画像を夫々複数個の局所領域に分割
する局所領域分割手段と、 前記2枚の画像夫々の複数個の局所領域から位
置合わせすべき1局所領域を夫々指定する位置補
正対象領域指定手段と、 この位置補正対象領域指定手段より指定された
一方の局所領域を格納すると共に他方の局所領域
をその近傍領域を含めて格納する局所領域メモリ
と、 この局所領域メモリの前記他方の局所領域を前
記近傍領域内で順次移動させて切り出す近傍領域
切り出し手段と、 前記一方の局所領域と前記近傍領域切り出し手
段より順次切り出される他方の局所領域との相関
係数値を夫々計算する相関係数計算手段、この相
関係数計算手段より計算された最大の相関係数値
に対応する他方の局所領域と前記一方の局所領域
とを位置合わせする手段、前記相関係数計算手段
より夫々計算された相関係数値の最大変化量が予
め与えられた閾値以下の場合に該位置合わせを終
了させる手段、前記最大変化量が前記閾値より大
の場合に前記最大の相関係数値に対応する他方の
局所領域の前記近傍領域内での移動切り出し、
夫々の相関係数値の計算及び位置合わせを繰り返
させる手段を備えた粗位置補正手段と、 この粗位置補正手段により位置合わせされ前記
近傍領域内切り出し手段より順次切り出される他
方の局所領域と前記一方の局所領域との画素毎の
濃度の絶対値誤差を夫々計算する絶対値誤差計算
手段、この絶対値誤差計算手段より計算された最
小の絶対値誤差に対応する他方の局所領域と前記
一方の局所領域とを位置合わせする手段、前記最
小の絶対値誤差が予め与えられた閾値以下の場合
に該位置合わせを終了させる手段、前記最小の絶
対値誤差が前記閾値より大の場合に前記最小の絶
対値誤差に対応する他方の局所領域の前記近傍領
域内での移動切り出し、夫々の絶対値誤差の計算
及び位置合わせを繰り返させる手段を備えた精密
位置補正手段と、 を具備したことを特徴とする画像の精密位置合わ
せ装置。
[Scope of Claims] 1. An image memory storing two images; local region dividing means for dividing each of the two images into a plurality of local regions; and a plurality of local regions of each of the two images. a position correction target area specifying means for specifying one local area to be aligned from the area; and a position correction target area specifying means for storing one local area specified by the position correction target area specifying means, and storing the other local area including its neighboring area. a local area memory for storing the other local area; a neighboring area cutting means for sequentially moving and cutting the other local area of the local area memory within the neighboring area; and a local area memory for sequentially cutting out the one local area and the neighboring area cutting means. Correlation coefficient calculation means for calculating a correlation coefficient with the other local area, and positioning the other local area and said one local area corresponding to the maximum correlation coefficient calculated by the correlation coefficient calculation means. means for aligning, means for terminating the alignment when the maximum amount of change in the correlation coefficients calculated by the correlation coefficient calculation means is less than or equal to a predetermined threshold; and means for terminating the alignment when the maximum amount of change is greater than the threshold moving the other local region corresponding to the maximum correlation coefficient within the neighboring region in the case;
a coarse position correction means having a means for repeating calculation and alignment of respective correlation coefficients, and another local area aligned by the coarse position correction means and successively cut out by the neighboring area cut-out means and one of the local areas. Absolute value error calculation means for calculating the absolute value error of density for each pixel with respect to the local region, the other local region corresponding to the minimum absolute value error calculated by the absolute value error calculation means, and said one local region. means for terminating the alignment when the minimum absolute value error is less than or equal to a predetermined threshold; and means for terminating the alignment when the minimum absolute value error is greater than the threshold; An image characterized by comprising: precision position correction means having means for repeating movement cutting out of the other local region corresponding to the error within the neighboring region, calculation of each absolute value error, and alignment; Precision alignment device.
JP59013371A 1984-01-30 1984-01-30 Precise picture positioning device Granted JPS60159972A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP59013371A JPS60159972A (en) 1984-01-30 1984-01-30 Precise picture positioning device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP59013371A JPS60159972A (en) 1984-01-30 1984-01-30 Precise picture positioning device

Publications (2)

Publication Number Publication Date
JPS60159972A JPS60159972A (en) 1985-08-21
JPH0527152B2 true JPH0527152B2 (en) 1993-04-20

Family

ID=11831225

Family Applications (1)

Application Number Title Priority Date Filing Date
JP59013371A Granted JPS60159972A (en) 1984-01-30 1984-01-30 Precise picture positioning device

Country Status (1)

Country Link
JP (1) JPS60159972A (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0695340B2 (en) * 1985-11-05 1994-11-24 株式会社東芝 Image matching method
JP6344039B2 (en) 2014-04-28 2018-06-20 富士通株式会社 Image display device, image display method, and program

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5134369A (en) * 1974-09-18 1976-03-24 Honda Motor Co Ltd Bureekideisuku no seiho
JPS54149433A (en) * 1978-05-16 1979-11-22 Mitsubishi Electric Corp Pattern collating unit
JPS57137978A (en) * 1981-02-20 1982-08-25 Toshiba Corp Pattern detecting device
JPS5886659A (en) * 1981-11-18 1983-05-24 Toshiba Corp Positioning device for picture
JPS58161082A (en) * 1982-03-19 1983-09-24 Fujitsu Ltd Collating system of seal impression

Also Published As

Publication number Publication date
JPS60159972A (en) 1985-08-21

Similar Documents

Publication Publication Date Title
Gruen et al. Geometrically constrained multiphoto matching
US7471809B2 (en) Method, apparatus, and program for processing stereo image
US6438507B1 (en) Data processing method and processing device
US6847392B1 (en) Three-dimensional structure estimation apparatus
US6961481B2 (en) Method and apparatus for image processing using sub-pixel differencing
EP0686942A2 (en) Stereo matching method and disparity measuring method
CN117237406B (en) Robot vision tracking method
CN115824190A (en) Vision and GPS-based target ship fusion positioning method
KR20230148999A (en) High-resolution image stitching method for semiconductor defect inspection
US7912321B1 (en) Image registration with uncertainty analysis
CN114690226A (en) Monocular vision distance measurement method and system based on carrier phase difference technology assistance
JPH08292014A (en) Measuring method of pattern position and device thereof
JP3066173B2 (en) Pattern matching method
JPH0527152B2 (en)
WO2002012830A1 (en) Height measurement apparatus
JPH04293178A (en) Configuration match system and its method
JP3066137B2 (en) Pattern matching method
US20260112054A1 (en) Method and apparatus for estimating rotation of video, and electronic device and storage medium
CN116916170A (en) Image acquisition device and image acquisition method
CN119146876A (en) Mobile and standing combined three-dimensional scanning method, device and equipment
US20040013299A1 (en) System and method for contrast enhanced registration with complex polynomial interpolation
JPS59218573A (en) Image alignment method
CN114739382B (en) Robot and mapping method, device and storage medium thereof
Ely et al. Model supported positioning
JPS6280768A (en) Stereo image processing method