JPH0410177A - 図形検出装置 - Google Patents
図形検出装置Info
- Publication number
- JPH0410177A JPH0410177A JP11410890A JP11410890A JPH0410177A JP H0410177 A JPH0410177 A JP H0410177A JP 11410890 A JP11410890 A JP 11410890A JP 11410890 A JP11410890 A JP 11410890A JP H0410177 A JPH0410177 A JP H0410177A
- Authority
- JP
- Japan
- Prior art keywords
- picture
- ternary
- pixel
- image
- straight line
- 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.)
- Granted
Links
Landscapes
- Image Analysis (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔産業上の利用分野]
本発明は、画像中から所望する直線成分を抽出する図形
検出装置に関する。
検出装置に関する。
〔従来の技術]
画像を用いた受動的な静止剛体の3次元形状計測手法と
して、エピポーラ画像を用いた3次元形状計測手法があ
げられる(Yamamoto、 “連続ステレオ画像
からの3次元情報の抽出”、信学論り。
して、エピポーラ画像を用いた3次元形状計測手法があ
げられる(Yamamoto、 “連続ステレオ画像
からの3次元情報の抽出”、信学論り。
Vol、J69−D、 No、11. pp、1631
4638.1986.)。エピポーラ画像とは、視軸に
垂直な方向に等速直進運動するカメラから連続的に静止
剛体を撮影して得られる画像群を時間軸方向に積み重ね
、これをエピポーラ線にそって切断した画像である。こ
のエピポーラ画像には通常はぼ平行に近い直線群が見ら
れる。この直線群は撮影された静止剛体の特徴点の見か
けの動きを表しており、この各直線の傾きを求めること
で、各特徴点までの距離を求めることができる。すなわ
ち撮影した静止剛体の3次元形状を計測することが可能
となる。
4638.1986.)。エピポーラ画像とは、視軸に
垂直な方向に等速直進運動するカメラから連続的に静止
剛体を撮影して得られる画像群を時間軸方向に積み重ね
、これをエピポーラ線にそって切断した画像である。こ
のエピポーラ画像には通常はぼ平行に近い直線群が見ら
れる。この直線群は撮影された静止剛体の特徴点の見か
けの動きを表しており、この各直線の傾きを求めること
で、各特徴点までの距離を求めることができる。すなわ
ち撮影した静止剛体の3次元形状を計測することが可能
となる。
このエピポーラ画像中から直線を検出する手法として、
従来、 KalmanfilterやRa+u+er
法などの直線追跡手法を用いていたが、直線追跡法はノ
イズや切れなどに弱いという欠点を有していた。
従来、 KalmanfilterやRa+u+er
法などの直線追跡手法を用いていたが、直線追跡法はノ
イズや切れなどに弱いという欠点を有していた。
一方、ノイズや切れに強い大局的な直線検出手法として
、 Hough変換があげられる。Hough変換に
よる直線検出手順は次の様なものである(B、Brow
n、“コンピュータビジョン”3日本コンピュータ協会
、 1987.)。
、 Hough変換があげられる。Hough変換に
よる直線検出手順は次の様なものである(B、Brow
n、“コンピュータビジョン”3日本コンピュータ協会
、 1987.)。
■°画像Eを構成する各画素の輝度値をF(x、y)(
x=O,、、、、M−1,y=o、 、、、、N−1)
とする。図形成分は画像E中の輝度変化が激しい画素や
等輝度の画素が連結して構成されるから9画像Eを一次
微分して輝度変化が大きい画素を「1」小さい画素を「
0」とする二値画像e°や、あるいは画像EにLapl
acian−Gaussianフィルタを通用して、零
交差点となる画素は「IJをそれ以外の画素は「0」と
する二値画像e”を求め、これらの二値画像にHoug
h変換を施していた。以後、これらの二値画像e”、e
”を総称して二値画像eと呼ぶことにする。二値画像e
の各画素をf (x、y)、(x=o、、、、、M−1
,y=o、、、、、N−1)とする。
x=O,、、、、M−1,y=o、 、、、、N−1)
とする。図形成分は画像E中の輝度変化が激しい画素や
等輝度の画素が連結して構成されるから9画像Eを一次
微分して輝度変化が大きい画素を「1」小さい画素を「
0」とする二値画像e°や、あるいは画像EにLapl
acian−Gaussianフィルタを通用して、零
交差点となる画素は「IJをそれ以外の画素は「0」と
する二値画像e”を求め、これらの二値画像にHoug
h変換を施していた。以後、これらの二値画像e”、e
”を総称して二値画像eと呼ぶことにする。二値画像e
の各画素をf (x、y)、(x=o、、、、、M−1
,y=o、、、、、N−1)とする。
■°直線を方程式y=mx+cで表し、パラメータm、
cで構成されるパラメータ空間を適当に量子化して配列
A (c、m)を構成し、その初期値を「0」とする。
cで構成されるパラメータ空間を適当に量子化して配列
A (c、m)を構成し、その初期値を「0」とする。
この配列Aをパラメータ配列と呼ぶことにする。
■゛二値画像eの画素値f (x、y)が「1」なら点
(x、y)を通る全ての直線のパラメータC1mについ
て、パラメータ配列A(c、m)の内容を「1」だけ加
算する。すなわちy=mx十〇を満たす全てのmとCと
について。
(x、y)を通る全ての直線のパラメータC1mについ
て、パラメータ配列A(c、m)の内容を「1」だけ加
算する。すなわちy=mx十〇を満たす全てのmとCと
について。
A (c、m):=A (c、m)+1とする。これを
すべてのf (x、y)、(x=O,、、、、M−1y
・帆、、、、N−1)について実行する。
すべてのf (x、y)、(x=O,、、、、M−1y
・帆、、、、N−1)について実行する。
■゛パラメータ配列A (c、m)上の極大点は二値画
像e上で一直線上にあるf (x、y) = 1である
画素の群れに対応するから、これを該画像における直線
として抽出する。
像e上で一直線上にあるf (x、y) = 1である
画素の群れに対応するから、これを該画像における直線
として抽出する。
ただしy=mx+cの形で直線をパラメータ表現すると
mが無限大となり得るから5通常はX CO5θ+y
sinθ8ρ の形で直線をパラメータ表現し、パラメータ配列A (
c、 m)の代わりにA(θ、ρ)を用いるのが一般的
である。この場合、Xとyとを固定するとρ−θ空間で
は正弦波曲線が得られるが、その他の点では上記手順と
変わらない。
mが無限大となり得るから5通常はX CO5θ+y
sinθ8ρ の形で直線をパラメータ表現し、パラメータ配列A (
c、 m)の代わりにA(θ、ρ)を用いるのが一般的
である。この場合、Xとyとを固定するとρ−θ空間で
は正弦波曲線が得られるが、その他の点では上記手順と
変わらない。
〔発明が解決しようとする課題]
上記手順からも明らかな様に、 Hough変換では
二値画像の上に直線をあてはめ、その直線に画素値が「
1jである点が何個台まれるかをパラメータ配列A(θ
、ρ)上で数えることで、直線検出を行っていた。した
がって、ノイズや切れの影響を受けにくり、大局的な直
線検出が可能である。
二値画像の上に直線をあてはめ、その直線に画素値が「
1jである点が何個台まれるかをパラメータ配列A(θ
、ρ)上で数えることで、直線検出を行っていた。した
がって、ノイズや切れの影響を受けにくり、大局的な直
線検出が可能である。
しかし第5図に示すように1画像中にほとんど平行に近
い密な直線群が存在する場合1点線で示す様な該直線群
を横切る直線が偽直線として抽出されてしまうという欠
点ををしていた。これは該直線群を偽直線 x cosθ +y sinθ18ρゝが横切る回数が
A(θ°、ρ“)に加算されるためρ−θ空間上に偽の
極大点が発生するためである。
い密な直線群が存在する場合1点線で示す様な該直線群
を横切る直線が偽直線として抽出されてしまうという欠
点ををしていた。これは該直線群を偽直線 x cosθ +y sinθ18ρゝが横切る回数が
A(θ°、ρ“)に加算されるためρ−θ空間上に偽の
極大点が発生するためである。
エピポーラ画像中に現れる直線成分は、被写体の形状や
表面の模様が複雑になり、特徴点が増えれば増えるほど
きわめて密な直線群となる。したがって従来の)tou
gh変換をそのまま該エピポーラ画像に適用すると、偽
の直線群を検出してしまい正確な3次元形状の計測は不
可能になるといった欠点を有していた。
表面の模様が複雑になり、特徴点が増えれば増えるほど
きわめて密な直線群となる。したがって従来の)tou
gh変換をそのまま該エピポーラ画像に適用すると、偽
の直線群を検出してしまい正確な3次元形状の計測は不
可能になるといった欠点を有していた。
本発明は上記の点を解決し、エピポーラ画像中からHo
ugh変換によって直線成分を検出し、3次元情報を抽
出することを目的としている。
ugh変換によって直線成分を検出し、3次元情報を抽
出することを目的としている。
(課題を解決するための手段)
本発明では9次のような手順で直線抽出を行う。
■ 二値画像eの各画素f (x、y)、 (X=0.
、、、、M1□y・O,、、、、N−1)に対して次の
ような三値化処理T : f (x、y) −’ g
(x+y)を行い、三値画像Gに変換する。
、、、、M1□y・O,、、、、N−1)に対して次の
ような三値化処理T : f (x、y) −’ g
(x+y)を行い、三値画像Gに変換する。
A(θ、ρ):=A(θ、ρ) + g (x、y)と
する。
する。
■ パラメータ配列A(θ、ρ)上の累積値の高い極大
点を直線として抽出する。
点を直線として抽出する。
上記手順からも明らかな様に1本発明では二値画像eを
三値画像Gに変換し、その画素+[g(x、y)をパラ
メータ配列A(θ、ρ)上に累加することを主な特徴と
している。
三値画像Gに変換し、その画素+[g(x、y)をパラ
メータ配列A(θ、ρ)上に累加することを主な特徴と
している。
・−−m−・−−−−m−−−−−・−・−・(1)た
だしg (x、y)+ (x=o、、、、、M−1,y
=o、、、、、N−1)は三値画像Gの画素である。
だしg (x、y)+ (x=o、、、、、M−1,y
=o、、、、、N−1)は三値画像Gの画素である。
■ パラメータ配列A(θ、ρ)を用意し、その初期値
を「0」とする。
を「0」とする。
■ 三値画像Gの画素値g(x、y)≠0ならば1点(
x、y)を通る全ての直線のパラメータθ。
x、y)を通る全ての直線のパラメータθ。
ρについて、パラメータ配列A(θ、ρ)に画素値g(
x、y)を加算する。すなわちx cosθ+y si
nθ=ρ −−−−−−−−・−・・−・・・−(2)
を満たす全てのρとθとについて 〔作用〕 二値画像eに即ち、第1図(A)図示の如き二値画像e
に代えて■の様な三値化処理を行うと。
x、y)を加算する。すなわちx cosθ+y si
nθ=ρ −−−−−−−−・−・・−・・・−(2)
を満たす全てのρとθとについて 〔作用〕 二値画像eに即ち、第1図(A)図示の如き二値画像e
に代えて■の様な三値化処理を行うと。
得られる三値画像Gは第1図(B)に示すような三値の
パターンをもつものとなる。したがって。
パターンをもつものとなる。したがって。
直線群を横切るような偽の直線に沿って画素値g(x、
y)を累加すると、必ずg(x+y)=1の画素を横切
るとg(x+y) = 1の画素を横切ることになる
から、偽の直線に対応する累積値は「0」か負になる。
y)を累加すると、必ずg(x+y)=1の画素を横切
るとg(x+y) = 1の画素を横切ることになる
から、偽の直線に対応する累積値は「0」か負になる。
一方、真の直線に沿って画素値g(x、y)を累加する
と、 g(x、y) =1の画素だけを横切るから、
累積値は大きな極大値となる。すなわち偽の直線に対応
する累積値は「0」近くに抑制され、一方、真の直線に
対応する累積値は大きな極大値となるから、真の直線だ
けが検出される。三値化処理Tによる三値化は、 f
(x、y) =1である画素の隣接画素に「−1」の値
を与えているが。
と、 g(x、y) =1の画素だけを横切るから、
累積値は大きな極大値となる。すなわち偽の直線に対応
する累積値は「0」近くに抑制され、一方、真の直線に
対応する累積値は大きな極大値となるから、真の直線だ
けが検出される。三値化処理Tによる三値化は、 f
(x、y) =1である画素の隣接画素に「−1」の値
を与えているが。
これによって二値画像e上でのf (x、y) = 1
である画素の連続性を、 Rough変換実施時のパ
ラメータ配列A(θ、ρ)上での累積値に反映すること
が可能になっている。
である画素の連続性を、 Rough変換実施時のパ
ラメータ配列A(θ、ρ)上での累積値に反映すること
が可能になっている。
従来の技術では、二値画像eに直接Rough変換を実
施するため、二値画像e上での画素の連続性をパラメー
タ配列A(θ、ρ)上での累積値に反映することができ
ず、したがって偽の直線を検出してしまうものであった
。
施するため、二値画像e上での画素の連続性をパラメー
タ配列A(θ、ρ)上での累積値に反映することができ
ず、したがって偽の直線を検出してしまうものであった
。
(実施例〕
第2図は第1図(B)図示の場合において真の直線に対
応する累積値が大であり、かつ偽の直線に対応する累積
値が小であることを示す。第3図は本発明の実施例を示
す。同図において、1は入力部、2は三値化処理部、2
−1は二値画像蓄積部、2−2は三値化部、3はHou
gh変換処理部。
応する累積値が大であり、かつ偽の直線に対応する累積
値が小であることを示す。第3図は本発明の実施例を示
す。同図において、1は入力部、2は三値化処理部、2
−1は二値画像蓄積部、2−2は三値化部、3はHou
gh変換処理部。
3−1は三値画像蓄積部、3−2はHough変換部。
3−3は加算部、3−4はパラメータ配列蓄積部。
3−5は極大値抽出部、4は出力部である。
パラメータ配列蓄積部3−4はBough変換実施時の
パラメータ配列A(θ、ρ)を蓄積する部分であり、初
期状態ではすべてのパラメータ配列の要素値は「OJで
ある。
パラメータ配列A(θ、ρ)を蓄積する部分であり、初
期状態ではすべてのパラメータ配列の要素値は「OJで
ある。
本装置を動作するには、まず直線成分を抽出したい二値
画像を入力部1に入力する。入力部1に入力された該二
値画像は二値画像蓄積部2−1に蓄積される。三値化部
2−2は該二値画像蓄積部2−1に蓄積された該二値画
像を(1)式で表される変換Tに従って画素毎に読みだ
して三値化し、結果を三値画像蓄積部3−1に蓄積する
。Hough変換部3−2は三値画像蓄積部3−1に蓄
積された三値画像を画素毎に読みだし2 もし該画素値
が「0」でなければ、(2)式を満足するすべてのパラ
メータの組(θ、ρ)と該画素の座標値(χy)とを順
次、加算部3−3に出力する。該加算部3−3では入力
された座標値(x、y)に対応する三値画像蓄積部3−
1中の画素値を読みだし。
画像を入力部1に入力する。入力部1に入力された該二
値画像は二値画像蓄積部2−1に蓄積される。三値化部
2−2は該二値画像蓄積部2−1に蓄積された該二値画
像を(1)式で表される変換Tに従って画素毎に読みだ
して三値化し、結果を三値画像蓄積部3−1に蓄積する
。Hough変換部3−2は三値画像蓄積部3−1に蓄
積された三値画像を画素毎に読みだし2 もし該画素値
が「0」でなければ、(2)式を満足するすべてのパラ
メータの組(θ、ρ)と該画素の座標値(χy)とを順
次、加算部3−3に出力する。該加算部3−3では入力
された座標値(x、y)に対応する三値画像蓄積部3−
1中の画素値を読みだし。
該画素値に、入力されたパラメータの&U(θρ)に対
応するパラメータ配列蓄積部3−4のパラメータ配列要
素を読みだして加算し、結果を該パラメータ配列蓄積部
3−4に蓄積する。Hough変換部3−2と加算部3
−3とが上記の動作を該三値画像蓄積部3−1内のすべ
ての画素について実行すると、極大値抽出部3−5は5
該パラメ一タ配列蓄積部3−4に蓄積されたパラメータ
配列要素の極大値を抽出し、該パラメータ配列要素に対
応するパラメータの組(θ、ρ)を出力部4に出力する
。
応するパラメータ配列蓄積部3−4のパラメータ配列要
素を読みだして加算し、結果を該パラメータ配列蓄積部
3−4に蓄積する。Hough変換部3−2と加算部3
−3とが上記の動作を該三値画像蓄積部3−1内のすべ
ての画素について実行すると、極大値抽出部3−5は5
該パラメ一タ配列蓄積部3−4に蓄積されたパラメータ
配列要素の極大値を抽出し、該パラメータ配列要素に対
応するパラメータの組(θ、ρ)を出力部4に出力する
。
この様な構成となっているから、二値画像中にほぼ平行
な密な直線群が存在しても、偽の直線を検出することな
く、真の直線群だけを検出することが可能である。
な密な直線群が存在しても、偽の直線を検出することな
く、真の直線群だけを検出することが可能である。
本実施例は直線を抽出する場合のものであるがこれは一
例であって、実際にはパラメータ表現されるすべての曲
線図形(たとえば第4図図示の様な曲線群)について同
様の処理構成で該曲線抽出が可能なことは明らかである
。
例であって、実際にはパラメータ表現されるすべての曲
線図形(たとえば第4図図示の様な曲線群)について同
様の処理構成で該曲線抽出が可能なことは明らかである
。
本発明によれば、密な直線群が存在するエピポーラ画像
に対してもHough変換の導入が可能となり、ノイズ
や切れなどの影響を受けにくい直線抽出、3次元形状計
測が可能となる。
に対してもHough変換の導入が可能となり、ノイズ
や切れなどの影響を受けにくい直線抽出、3次元形状計
測が可能となる。
第1図(A)(B )は本発明の原理説明図、第2図は
本発明の場合の利点を説明する説明図、第3図は本発明
の実施例、第4図は蜜な曲線群を表す図、第5図は従来
の問題点を説明する図を示す。 図中、1は入力部、2は三値化処理部、2−1は二値画
像蓄積部、2−2は三値化部、3はHough変換処理
部、3−1は三値画像蓄積部、3−2はHough変換
部、3−3は加算部、3−4はパラメータ配列蓄積部、
3−5は極大値抽出部4は出力部を表す。
本発明の場合の利点を説明する説明図、第3図は本発明
の実施例、第4図は蜜な曲線群を表す図、第5図は従来
の問題点を説明する図を示す。 図中、1は入力部、2は三値化処理部、2−1は二値画
像蓄積部、2−2は三値化部、3はHough変換処理
部、3−1は三値画像蓄積部、3−2はHough変換
部、3−3は加算部、3−4はパラメータ配列蓄積部、
3−5は極大値抽出部4は出力部を表す。
Claims (1)
- 【特許請求の範囲】 入力された「0」と「1」との二値から構成される二値
画像に対して、該二値画像中の「1」の画素を正値に、
「1」の画素に隣接する「0」の画素を負値に、それ以
外の画素を「0」とする三値画像を出力する三値化処理
部と、 該三値化処理部から出力された三値画像に Hough変換を行い線成分を抽出するHough変換
処理部とから構成され、 該二値画像中に存在する密な線群についての各画素の累
積値を抽出して、真の線成分だけを検出するようにした ことを特徴とする図形検出装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP11410890A JP2912676B2 (ja) | 1990-04-27 | 1990-04-27 | 図形検出装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP11410890A JP2912676B2 (ja) | 1990-04-27 | 1990-04-27 | 図形検出装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0410177A true JPH0410177A (ja) | 1992-01-14 |
| JP2912676B2 JP2912676B2 (ja) | 1999-06-28 |
Family
ID=14629326
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP11410890A Expired - Fee Related JP2912676B2 (ja) | 1990-04-27 | 1990-04-27 | 図形検出装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2912676B2 (ja) |
-
1990
- 1990-04-27 JP JP11410890A patent/JP2912676B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2912676B2 (ja) | 1999-06-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jeon et al. | Depth from a light field image with learning-based matching costs | |
| CN105427298B (zh) | 基于各向异性梯度尺度空间的遥感图像配准方法 | |
| CN104197861A (zh) | 基于结构光灰度向量的三维数字成像方法 | |
| JPH0935061A (ja) | 画像処理方法 | |
| Kimura et al. | A convolver-based real-time stereo machine (SAZAN) | |
| CN118537377A (zh) | 一种多镜头红外与可见光的图像配准方法 | |
| Yang et al. | Robust structured light 3D imaging with two fringe patterns using recurrent classification neural network | |
| Seetharaman et al. | A piecewise affine model for image registration in nonrigid motion analysis | |
| JPH04116481A (ja) | 磁場測定方法、およびその装置 | |
| JP4918615B2 (ja) | 対象個数検出装置および対象個数検出方法 | |
| JPH0410177A (ja) | 図形検出装置 | |
| US5008744A (en) | Method and apparatus for picture motion measurement using picture correlation as a function of displacement | |
| CN111696143A (zh) | 一种事件数据的配准方法与系统 | |
| JPH0875454A (ja) | 測距装置 | |
| CN106056599B (zh) | 一种基于物体深度数据的物体识别算法及装置 | |
| Papanikolaou et al. | Colour digital image correlation method for monitoring of cultural heritage objects with natural texture | |
| JP3074292B2 (ja) | テンプレートマッチング処理方法 | |
| Nomura et al. | Edge location to subpixel precision and analysis | |
| Gadelmawla et al. | On-line measurement of product dimensions using computer vision | |
| JP3060114B2 (ja) | 対応点探索方法 | |
| Tang et al. | Weak texture 3D reconstruction method based on laser-assisted scanning and moving platforms | |
| JPH02138677A (ja) | 画像特徴点検出方法 | |
| Rajendran et al. | A method for subpixel accuracy detection along the edges and ridges using colors for images | |
| JP2840748B2 (ja) | 平面検出装置 | |
| CN119323591A (zh) | 一种结构全局三维振型的非接触式识别方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090409 Year of fee payment: 10 |
|
| LAPS | Cancellation because of no payment of annual fees |