JPH0799327B2 - カメラの焦点距離較正方法及び装置 - Google Patents
カメラの焦点距離較正方法及び装置Info
- Publication number
- JPH0799327B2 JPH0799327B2 JP5008997A JP899793A JPH0799327B2 JP H0799327 B2 JPH0799327 B2 JP H0799327B2 JP 5008997 A JP5008997 A JP 5008997A JP 899793 A JP899793 A JP 899793A JP H0799327 B2 JPH0799327 B2 JP H0799327B2
- Authority
- JP
- Japan
- Prior art keywords
- focal length
- camera
- vanishing point
- calibration
- image
- 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 - Fee Related
Links
Landscapes
- Automatic Focus Adjustment (AREA)
- Closed-Circuit Television Systems (AREA)
- Testing, Inspecting, Measuring Of Stereoscopic Televisions And Televisions (AREA)
- Length Measuring Devices By Optical Means (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、画像に写っている物体
の３次元位置や方向を計算するにあたって必要な、カメ
ラの焦点距離較正方法に関し、特に最小の誤差で焦点距
離を得るための、カメラの焦点距離較正方法及び装置に
関する。
の３次元位置や方向を計算するにあたって必要な、カメ
ラの焦点距離較正方法に関し、特に最小の誤差で焦点距
離を得るための、カメラの焦点距離較正方法及び装置に
関する。
【０００２】
【従来の技術】従来、この種のカメラの焦点距離較正方
法は、例えば、「１９９１年１２月、日本ロボット学会
誌、第９巻 第７号 ８１３〜８２９頁」に示されるよ
うに、カメラから取得された画像によって、工業製品の
自動検査・組立や移動ロボットの制御などを目的とする
際に、基本となるカメラパラメータ（焦点距離，視点の
３次元位置，カメラ座標系の３次元方向など）のうち、
最も基本となるカメラの焦点距離の較正方法として用い
られている。このカメラの焦点距離較正方法について、
図４を参照して説明する。
法は、例えば、「１９９１年１２月、日本ロボット学会
誌、第９巻 第７号 ８１３〜８２９頁」に示されるよ
うに、カメラから取得された画像によって、工業製品の
自動検査・組立や移動ロボットの制御などを目的とする
際に、基本となるカメラパラメータ（焦点距離，視点の
３次元位置，カメラ座標系の３次元方向など）のうち、
最も基本となるカメラの焦点距離の較正方法として用い
られている。このカメラの焦点距離較正方法について、
図４を参照して説明する。
【０００３】図４にあるようなカメラの撮像モデル５を
考える。Ｚ軸がカメラの光軸と一致するようにＸＹＺ座
標系をとり、原点０をレンズの中心とみなし、平面Ｚ＝
ｆを画像面とする。ｆはレンズの中心と撮像面との距離
であり、「焦点距離」と呼ぶ。
考える。Ｚ軸がカメラの光軸と一致するようにＸＹＺ座
標系をとり、原点０をレンズの中心とみなし、平面Ｚ＝
ｆを画像面とする。ｆはレンズの中心と撮像面との距離
であり、「焦点距離」と呼ぶ。
【０００４】正方格子５１を描いた表示板を用意し、そ
の画像を入力する。正方格子５１は画像面上に投影さ
れ、投影像５２を生ずる。この投影像５２の直線部分の
エッジを画像処理によって検出し、直線を最小２乗法に
よって当てはめる。
の画像を入力する。正方格子５１は画像面上に投影さ
れ、投影像５２を生ずる。この投影像５２の直線部分の
エッジを画像処理によって検出し、直線を最小２乗法に
よって当てはめる。
【０００５】ここで、焦点距離ｆを求めるために、仮の
【０００６】
【外１】
【０００７】を用いて、消失点Ｐ1 ，Ｐ2 の
【０００８】
【外２】
【０００９】を最小２乗法によって推定する。ただし、
図４にあるように、画像上の点Ｐi の「Ｎベクトル」と
は、視点から点Ｐi へ向う単位ベクトルである。
図４にあるように、画像上の点Ｐi の「Ｎベクトル」と
は、視点から点Ｐi へ向う単位ベクトルである。
【００１０】焦点距離ｆは、
【００１１】
【数１】
【００１２】を解くことによって求まる。ただし、
【００１３】
【外３】
【００１４】である。
【００１５】つぎに焦点距離ｆの分散Ｖ［ｆ］を計算す
る。それは
る。それは
【００１６】
【数２】
【００１７】で求まる。ここで、
【００１８】
【外４】
【００１９】に対する共分散行列であり、
【００２０】
【数３】
【００２１】で定義される。ただし、
【００２２】
【外５】
【００２３】であり、Ｅ［・・・］は期待値を表す。
【００２４】以上を撮像角度を変えながら繰り返し、異
なる撮像角度に対する複数の焦点距離と各々の焦点距離
に対する分散を得る。測定をｎ回行ったとし、得られた
焦点距離と分散をｆi ，Ｖ［ｆi ］，ｉ＝１，・・・，
ｎとすると、最も確からしい
なる撮像角度に対する複数の焦点距離と各々の焦点距離
に対する分散を得る。測定をｎ回行ったとし、得られた
焦点距離と分散をｆi ，Ｖ［ｆi ］，ｉ＝１，・・・，
ｎとすると、最も確からしい
【００２５】
【外６】
【００２６】は、最小２乗法により次のように求まる。
【００２７】
【数４】
【００２８】である。
【００２９】
【発明が解決しようとする課題】以上に述べた従来技術
によると、カメラから入力された画像によりカメラの焦
点距離を求めるには、 １．式（１）により焦点距離ｆを求め、 ２．式（２）により分散Ｖ［ｆ］を求め、 ３．異なる撮像角度に対して手順１、２を行い、 ４．式（３）、（４）により最も確からしい
によると、カメラから入力された画像によりカメラの焦
点距離を求めるには、 １．式（１）により焦点距離ｆを求め、 ２．式（２）により分散Ｖ［ｆ］を求め、 ３．異なる撮像角度に対して手順１、２を行い、 ４．式（３）、（４）により最も確からしい
【００３０】
【外７】
【００３１】を得る、のように手順を踏む必要があり、
手間がかかる。また、求められた焦点距離がどの程度信
頼できるか明らかでない。さらに、較正用パターンとし
て正方格子を用いているが、これは信頼性の高い焦点距
離を得るためには最適なものではない。
手間がかかる。また、求められた焦点距離がどの程度信
頼できるか明らかでない。さらに、較正用パターンとし
て正方格子を用いているが、これは信頼性の高い焦点距
離を得るためには最適なものではない。
【００３２】本発明の目的は、較正用パターンとして斜
方格子を用い、１回の測定により、正方格子を用いる方
法よりも、信頼性の高い焦点距離を得ることのできる、
カメラの焦点距離較正方法を提供することにある。
方格子を用い、１回の測定により、正方格子を用いる方
法よりも、信頼性の高い焦点距離を得ることのできる、
カメラの焦点距離較正方法を提供することにある。
【００３３】また、本発明の他の目的は、カメラの焦点
距離をその信頼区間と共に得ることができる較正方法を
提供することにある。
距離をその信頼区間と共に得ることができる較正方法を
提供することにある。
【００３４】さらに、本発明の他の目的は、このような
方法を実施するカメラの焦点距離較正装置を提供するこ
とにある。
方法を実施するカメラの焦点距離較正装置を提供するこ
とにある。
【００３５】
【課題を解決するための手段】本発明は、平面上に描か
れた焦点距離較正用パターンをカメラで撮像して得られ
る画像からカメラの焦点距離を得る、カメラの焦点距離
較正方法において、平面上に描かれた２組の斜交した平
行線からなる較正用格子パターンをカメラで撮像し、入
力画像から、撮像された前記２組の平行線の消失点を計
算し、前記消失点から焦点距離を計算することを特徴と
する。
れた焦点距離較正用パターンをカメラで撮像して得られ
る画像からカメラの焦点距離を得る、カメラの焦点距離
較正方法において、平面上に描かれた２組の斜交した平
行線からなる較正用格子パターンをカメラで撮像し、入
力画像から、撮像された前記２組の平行線の消失点を計
算し、前記消失点から焦点距離を計算することを特徴と
する。
【００３６】また、本発明は、平面上に描かれた焦点距
離較正用パターンをカメラで撮像して得られる画像から
カメラの焦点距離を得る、カメラの焦点距離較正装置に
おいて、平面上に描かれた２組の斜交した平行線からな
る較正用格子パターンと、この較正用格子パターンを撮
像するカメラと、入力画像から、撮像された前記２組の
平行線の消失点を計算する消失点演算手段と、前記消失
点から焦点距離を計算する焦点距離演算手段とを備える
ことを特徴とする。
離較正用パターンをカメラで撮像して得られる画像から
カメラの焦点距離を得る、カメラの焦点距離較正装置に
おいて、平面上に描かれた２組の斜交した平行線からな
る較正用格子パターンと、この較正用格子パターンを撮
像するカメラと、入力画像から、撮像された前記２組の
平行線の消失点を計算する消失点演算手段と、前記消失
点から焦点距離を計算する焦点距離演算手段とを備える
ことを特徴とする。
【００３７】
【作用】本発明における、３次元空間中の２次元平面上
に描かれた格子パターンから焦点距離を求める方法と、
信頼性の高い焦点距離を得るための格子パターンの構成
と撮像方法について、図５，図６，図７，図８，図９を
参照して述べる。
に描かれた格子パターンから焦点距離を求める方法と、
信頼性の高い焦点距離を得るための格子パターンの構成
と撮像方法について、図５，図６，図７，図８，図９を
参照して述べる。
【００３８】３次元空間に置かれた平面上に、Ｎ本の平
行線Ｌi ，ｉ＝１，・・・，Ｎと、Ｎ′本の平行線
Ｌi ′，ｉ＝１，・・・，Ｎ′がある。このとき、画像
上のＬi，Ｌi ′の射影をｌi ，ｌi ′とする。
行線Ｌi ，ｉ＝１，・・・，Ｎと、Ｎ′本の平行線
Ｌi ′，ｉ＝１，・・・，Ｎ′がある。このとき、画像
上のＬi，Ｌi ′の射影をｌi ，ｌi ′とする。
【００３９】ｌi を延長してできる交点、すなわち消失
点の
点の
【００４０】
【外８】
【００４１】とする。図５に示すように、直線ｌi の
【００４２】
【外９】
【００４３】とすると、
【００４４】
【外１０】
【００４５】は次式により求めることができる。
【００４６】
【数５】
【００４７】ただしｗi は直線ｌi を当てはめたエッジ
の長さ（単位は画素）、
の長さ（単位は画素）、
【００４８】
【外１１】
【００４９】はその方向を示す単位ベクトルであり、
【００５０】
【外１２】
【００５１】はその中心点のＮベクトルである。また、
ｆは焦点距離の推定値であり、適当な値を用いればよ
い。
ｆは焦点距離の推定値であり、適当な値を用いればよ
い。
【００５２】ここで、直線ｌi の消失点の
【００５３】
【外１３】
【００５４】を次の「くりこみ法」によって計算する。
【００５５】「くりこみ法」： １．ｃ＝１、Ｗi ＝１，ｉ＝１，・・・，Ｎとおく。
【００５６】２．
【００５７】
【外１４】
【００５８】を次のように計算する。
【００５９】
【数６】
【００６０】そして、その最小固有値λm 、対応する
【００６１】
【外１５】
【００６２】を計算する。
【００６３】３．ｃ、Ｗi を次のように更新する。
【００６４】
【数７】
【００６５】４．更新が収束していなければ２へ戻る。
収束していれば
収束していれば
【００６６】
【外１６】
【００６７】である。
【００６８】そして、
【００６９】
【外１７】
【００７０】の最大固有値λu 、対応する
【００７１】
【外１８】
【００７２】中間固有値λv 、対応する
【００７３】
【外１９】
【００７４】を計算する。
【００７５】５．
【００７６】
【外２０】
【００７７】を次式により計算する。
【００７８】
【数８】
【００７９】以上から、直線ｌi の
【００８０】
【外２１】
【００８１】が求まる。
【００８２】上では、「くりこみ法」を用いる方法を述
べたが、これとは別に、「モーメント行列」を用いる方
法もある。「モーメント行列」Ｍは、
べたが、これとは別に、「モーメント行列」を用いる方
法もある。「モーメント行列」Ｍは、
【００８３】
【数９】
【００８４】で定義される。ここで、重みＷi は
【００８５】
【数１０】
【００８６】で与えられ、消失点の
【００８７】
【外２２】
【００８８】は、この「モーメント行列」を対角化した
ときの最小固有値に対応する固有ベクトルになる。この
方法で注意すべきことは、重みの式に求めるべき
ときの最小固有値に対応する固有ベクトルになる。この
方法で注意すべきことは、重みの式に求めるべき
【００８９】
【外２３】
【００９０】が含まれていることである。ただし、実際
の計算においては、重みに含まれる
の計算においては、重みに含まれる
【００９１】
【外２４】
【００９２】として、全ての重みを１として計算した推
定値を用いれば十分である。このときにも、「くりこみ
法」の場合と同様に、最大固有値λu と対応する
定値を用いれば十分である。このときにも、「くりこみ
法」の場合と同様に、最大固有値λu と対応する
【００９３】
【外２５】
【００９４】中間固有値λv と対応する
【００９５】
【外２６】
【００９６】を計算し、式（９）によって共分散行列を
計算する。
計算する。
【００９７】図６を参照し、「くりこみ法」または「モ
ーメント行列」を用いる方法で得られる３つの固有ベク
トル、
ーメント行列」を用いる方法で得られる３つの固有ベク
トル、
【００９８】
【外２７】
【００９９】について述べる。まず、最大固有値に対応
する
する
【０１００】
【外２８】
【０１０１】は、直線ｌi の仮想的な中心線ｌG のＮベ
クトルであり、中間固有値に対応する
クトルであり、中間固有値に対応する
【０１０２】
【外２９】
【０１０３】に等しい。
【０１０４】さて、以上により直線ｌi の
【０１０５】
【外３０】
【０１０６】と、その
【０１０７】
【外３１】
【０１０８】が求められた。これと同じ操作を直線
ｌi ′に対しても行い、
ｌi ′に対しても行い、
【０１０９】
【外３２】
【０１１０】を求める。
【０１１１】次に焦点距離の求め方について述べる。図
７に示すように、３次元空間中で、斜方格子の交わる角
度をαとする。求められた２つの消失点の
７に示すように、３次元空間中で、斜方格子の交わる角
度をαとする。求められた２つの消失点の
【０１１２】
【外３３】
【０１１３】を
【０１１４】
【数１１】
【０１１５】とする。はじめに３つの量Ａ、Ｂ、Ｃを定
義する。
義する。
【０１１６】
【数１２】
【０１１７】次に、２次方程式
【０１１８】
【数１３】
【０１１９】の二つの解を計算し、そのうち
【０１２０】
【数１４】
【０１２１】であるものを選ぶ。すると、真の
【０１２２】
【外３４】
【０１２３】は次式により求まる：
【０１２４】
【数１５】
【０１２５】焦点距離ｆの分散Ｖ［ｆ］は、
【０１２６】
【外３５】
【０１２７】を用いて、
【０１２８】
【数１６】
【０１２９】で与えられ、焦点距離の水準（１００−
ａ）％の信頼区間は
ａ）％の信頼区間は
【０１３０】
【数１７】
【０１３１】で与えられる。λa は標準正規分布のａ％
点である。（例えばλ5 ＝１．９６）。また、κは「解
像度定数」と呼ばれる定数であり、
点である。（例えばλ5 ＝１．９６）。また、κは「解
像度定数」と呼ばれる定数であり、
【０１３２】
【数１８】
【０１３３】で表される量である。εは「画像精度」と
呼ばれる定数であり、各エッジ点の平均的ズレ（標準偏
差）を画素の単位で表したものである。γは「エッジ密
度」すなわち単位長さ当たりのエッジ点の数（単位は画
素の逆数）である。κは通常κ＝１として差し支えな
い。
呼ばれる定数であり、各エッジ点の平均的ズレ（標準偏
差）を画素の単位で表したものである。γは「エッジ密
度」すなわち単位長さ当たりのエッジ点の数（単位は画
素の逆数）である。κは通常κ＝１として差し支えな
い。
【０１３４】次に、焦点距離の分散Ｖ［ｆ］を最も小さ
くするような、平行線の配置について述べる。φi を
くするような、平行線の配置について述べる。φi を
【０１３５】
【外３６】
【０１３６】の間の角度とし、θi をｉ番目のエッジの
中心点の
中心点の
【０１３７】
【外３７】
【０１３８】と消失点の
【０１３９】
【外３８】
【０１４０】のなす角度とする。このとき、消失点の
【０１４１】
【外３９】
【０１４２】は次のように書くことができる。
【０１４３】
【数１９】
【０１４４】上式をＶ［ｆ］の表式に代入すると、
【０１４５】
【数２０】
【０１４６】を得る。各エッジの中点が画像原点の近傍
にあれば、θi ，θi ′はそれぞれ図８にあるθ，θ′
で近似できる。したがって、
にあれば、θi ，θi ′はそれぞれ図８にあるθ，θ′
で近似できる。したがって、
【０１４７】
【数２１】
【０１４８】となる。
【０１４９】図８を参照して、空間中におかれた格子パ
ターンの配置について考える。ｒを視点と格子パターン
の中心の距離、ｌを格子パターンの大きさ、ｄを各々の
格子の間隔とする。
ターンの配置について考える。ｒを視点と格子パターン
の中心の距離、ｌを格子パターンの大きさ、ｄを各々の
格子の間隔とする。
【０１５０】
【外４０】
【０１５１】は２つの消失点に対応するＮベクトルであ
り、これは空間中のそれぞれの平行線の方向ベクトルに
なっているので、θとθ′で、カメラによる撮像方向
と、較正用パターンが描かれている平面の関係を規定で
きる。
り、これは空間中のそれぞれの平行線の方向ベクトルに
なっているので、θとθ′で、カメラによる撮像方向
と、較正用パターンが描かれている平面の関係を規定で
きる。
【０１５２】このとき、焦点距離の分散Ｖ［ｆ］は
【０１５３】
【数２２】
【０１５４】と表される。ただし、
【０１５５】
【数２３】
【０１５６】である。
【０１５７】この分散Ｖ［ｆ］は、
【０１５８】
【数２４】
【０１５９】のとき最小となる。図９に示すように、カ
メラにより取得された画像上での、格子の交わる角度を
γとすると、
メラにより取得された画像上での、格子の交わる角度を
γとすると、
【０１６０】
【数２５】
【０１６１】が成り立つ。
【０１６２】以上から分散Ｖ［ｆ］を最小にする角度α
は、
は、
【０１６３】
【数２６】
【０１６４】と求まり、これに対応する格子パターンの
最適配置は、
最適配置は、
【０１６５】
【数２７】
【０１６６】となる。また、このときγ＝９０度であ
る。
る。
【０１６７】したがって、焦点距離の分散を最小にす
る、すなわち最も信頼性の高い焦点距離を得るための較
正用パターンとして、２組の平行線の交わる角度が式
（３１）のαであるものを用いる。さらに、較正用パタ
ーンに対するカメラの撮像方向としては、図８に示す
θ，θ′が、式（３２）の値を持つようにすればよい。
しかし、このように配置するよりも、撮像した画像中の
格子線の交わる角度γが、上に示したように、９０度に
なるようにパターンを配置するほうが簡単である。
る、すなわち最も信頼性の高い焦点距離を得るための較
正用パターンとして、２組の平行線の交わる角度が式
（３１）のαであるものを用いる。さらに、較正用パタ
ーンに対するカメラの撮像方向としては、図８に示す
θ，θ′が、式（３２）の値を持つようにすればよい。
しかし、このように配置するよりも、撮像した画像中の
格子線の交わる角度γが、上に示したように、９０度に
なるようにパターンを配置するほうが簡単である。
【０１６８】
【実施例】次に、本発明の実施例について図面を参照し
て説明する。図２に、本発明における較正用格子パター
ンの例を示す。図２を参照すると、較正用格子パターン
の例は、向かい合う辺の中点を結ぶ線分を有する菱形で
ある。この較正用格子パターンは、焦点距離の分散を最
小にするために、式（３１）に示されるように２ｓｉｎ
-1√（３／１０）度の角度を有するものである。
て説明する。図２に、本発明における較正用格子パター
ンの例を示す。図２を参照すると、較正用格子パターン
の例は、向かい合う辺の中点を結ぶ線分を有する菱形で
ある。この較正用格子パターンは、焦点距離の分散を最
小にするために、式（３１）に示されるように２ｓｉｎ
-1√（３／１０）度の角度を有するものである。
【０１６９】図３に、入力画像４を示す。図３を参照す
ると、入力画像４は、９０度の角度を有するように撮像
される。
ると、入力画像４は、９０度の角度を有するように撮像
される。
【０１７０】図１は、本発明の実施例を示すブロック図
である。図１を参照すると、本発明の実施例は、データ
入力装置１と、プログラム制御により動作するデータ処
理装置２と、ディスプレイ装置等の出力装置３とから構
成される。
である。図１を参照すると、本発明の実施例は、データ
入力装置１と、プログラム制御により動作するデータ処
理装置２と、ディスプレイ装置等の出力装置３とから構
成される。
【０１７１】データ入力装置１は、較正用格子パターン
１１と、その較正用格子パターン１１を撮像する手段で
あり、かつ、焦点距離を求めるべきカメラ１２とを備え
る。
１１と、その較正用格子パターン１１を撮像する手段で
あり、かつ、焦点距離を求めるべきカメラ１２とを備え
る。
【０１７２】データ処理装置２は、カメラ１２により入
力された画像から、画像処理によって格子線のエッジを
求めるエッジ検出手段２１と、消失点のＮベクトル及び
そのＮベクトルの分散を計算する消失点算出手段２２
と、消失点のＮベクトルから焦点距離を計算する焦点距
離算出手段２３と、消失点のＮベクトルと、Ｎベクトル
の分散と、焦点距離の値から焦点距離の分散を計算する
分散算出手段２４と、焦点距離の分散から焦点距離の信
頼区間を計算する信頼区間算出手段２５とを備える。
力された画像から、画像処理によって格子線のエッジを
求めるエッジ検出手段２１と、消失点のＮベクトル及び
そのＮベクトルの分散を計算する消失点算出手段２２
と、消失点のＮベクトルから焦点距離を計算する焦点距
離算出手段２３と、消失点のＮベクトルと、Ｎベクトル
の分散と、焦点距離の値から焦点距離の分散を計算する
分散算出手段２４と、焦点距離の分散から焦点距離の信
頼区間を計算する信頼区間算出手段２５とを備える。
【０１７３】次に、図１、図２、及び図３を参照して、
本実施例の動作について説明する。
本実施例の動作について説明する。
【０１７４】カメラ１２によって取得された較正用斜方
格子パターン１１の入力画像４は、出力装置３とエッジ
検出手段２１に供給される。
格子パターン１１の入力画像４は、出力装置３とエッジ
検出手段２１に供給される。
【０１７５】エッジ検出手段２１では、入力画像４から
画像処理により、直線のエッジを検出し、各エッジの方
向ベクトルと、中心点のＮベクトルと、画素単位の長さ
が計算される。
画像処理により、直線のエッジを検出し、各エッジの方
向ベクトルと、中心点のＮベクトルと、画素単位の長さ
が計算される。
【０１７６】消失点算出手段２２では、式（６）にした
がって、仮の焦点距離ｆを用いて直線の
がって、仮の焦点距離ｆを用いて直線の
【０１７７】
【外４１】
【０１７８】を計算し、「くりこみ法」を用いて、最小
固有値λm と対応する
固有値λm と対応する
【０１７９】
【外４２】
【０１８０】最大固有値λu と対応する
【０１８１】
【外４３】
【０１８２】中間固有値λv と対応する
【０１８３】
【外４４】
【０１８４】を計算し、式（９）によって
【０１８５】
【外４５】
【０１８６】を計算する。
【０１８７】焦点距離算出手段２３では、式（１２）〜
（１９）により、
（１９）により、
【０１８８】
【外４６】
【０１８９】が計算される。
【０１９０】分散算出手段２４では式（２０）により、
焦点距離算出手段２３で算出された
焦点距離算出手段２３で算出された
【０１９１】
【外４７】
【０１９２】を計算する。
【０１９３】信頼区間算出手段２５では、式（２１）に
より焦点距離算出手段２３で算出された
より焦点距離算出手段２３で算出された
【０１９４】
【外４８】
【０１９５】の信頼区間が計算される。
【０１９６】以上の実施例では、画像処理によって格子
線のエッジを抽出したが、これは入力画像を見て人間が
行っても良い。また、消失点の算出には「くりこみ法」
を用いたが、これは式（１０）、（１１）に示した「モ
ーメント行列」の方法を用いても実施することができ
る。
線のエッジを抽出したが、これは入力画像を見て人間が
行っても良い。また、消失点の算出には「くりこみ法」
を用いたが、これは式（１０）、（１１）に示した「モ
ーメント行列」の方法を用いても実施することができ
る。
【０１９７】
【発明の効果】以上説明したように、本発明によるカメ
ラの焦点距離較正方法及び装置では、平面上に描かれた
斜方格子を用いることによって、従来のように複数の測
定をすることなく、１回の測定によって、正方格子を用
いるよりも信頼性の高い焦点距離と、さらに焦点距離の
信頼区間を得ることができる。
ラの焦点距離較正方法及び装置では、平面上に描かれた
斜方格子を用いることによって、従来のように複数の測
定をすることなく、１回の測定によって、正方格子を用
いるよりも信頼性の高い焦点距離と、さらに焦点距離の
信頼区間を得ることができる。
【図１】本発明の実施例を示すブロック図である。
【図２】焦点距離較正用格子パターンの例を示す説明図
である。
である。
【図３】格子パターンの入力画像を示す説明図である。
【図４】撮像モデルを示す説明図である。
【図５】エッジのＮベクトルと方向ベクトルを示す説明
図である。
図である。
【図６】３つの固有ベクトルを示す説明図である。
【図７】角度αを示す説明図である。
【図８】格子パターンの配置を示す説明図である。
【図９】角度γを示す説明図である。
１ データ入力装置 ２ データ処理装置 ３ 出力装置 １１ 較正用斜方格子パタン １２ カメラ ２１ エッジ検出手段 ２２ 消失点算出手段 ２３ 焦点距離算出手段 ２４ 分散算出手段 ２５ 信頼区間算出手段
フロントページの続き (51)Int.Cl.6 識別記号 庁内整理番号 ＦＩ 技術表示箇所 Ｈ０４Ｎ 17/00 Ｋ
Claims (10)
- 【請求項１】平面上に描かれた焦点距離較正用パターン
をカメラで撮像して得られる画像からカメラの焦点距離
を得る、カメラの焦点距離較正方法において、 平面上に描かれた２組の斜交した平行線からなる較正用
格子パターンをカメラで撮像し、 入力画像から、撮像された前記２組の平行線の消失点を
計算し、 前記消失点から焦点距離を計算することを特徴とするカ
メラの焦点距離較正方法。 - 【請求項２】前記較正用斜方格子パターンにおける２組
の平行線が２ｓｉｎ-1√（３／１０）の角度を有し、カ
メラによる入力画像において２組の平行線が直角に交わ
ることを特徴とする請求項１記載のカメラの焦点距離較
正方法。 - 【請求項３】「くりこみ法」により消失点を計算するこ
とを特徴とする、請求項１又は２記載のカメラの焦点距
離較正方法。 - 【請求項４】「モーメント行列の方法」により消失点を
計算することを特徴とする、請求項１又は２記載のカメ
ラの焦点距離較正方法。 - 【請求項５】平面上に描かれた焦点距離較正用パターン
をカメラで撮像して得られる画像からカメラの焦点距離
を得る、カメラの焦点距離較正装置において、 平面上に描かれた２組の斜交した平行線からなる較正用
格子パターンと、 この較正用格子パターンを撮像するカメラと、 入力画像から、撮像された前記２組の平行線の消失点を
計算する消失点演算手段と、 前記消失点から焦点距離を計算する焦点距離演算手段と
を備えることを特徴とするカメラの焦点距離較正装置。 - 【請求項６】前記焦点距離演算手段が、 前記消失点演算手段により求められた消失点から、焦点
距離を計算する焦点距離算出手段と、 焦点距離の分散を計算する分散算出手段と、 焦点距離の分散から焦点距離の信頼区間を計算する信頼
区間算出手段とを有することを特徴とする請求項５記載
のカメラの焦点距離較正装置。 - 【請求項７】前記較正用斜方格子パターンにおける２組
の平行線が２ｓｉｎ-1√（３／１０）の角度を有し、カ
メラによる入力画像において２組の平行線が直角に交わ
ること特徴とする請求項５又は６記載のカメラの焦点距
離較正装置。 - 【請求項８】前記消失点演算手段が、 画像処理により入力画像から格子線を抽出するエッジ検
出手段と、 前記エッジ検出手段により得られる格子線から消失点を
求める消失点演算手段とを有することを特徴とする請求
項５又は６又は７記載のカメラの焦点距離較正装置。 - 【請求項９】前記消失点演算手段が、「くりこみ法」に
より消失点を計算することを特徴とする、請求項５又は
６又は７記載のカメラの焦点距離較正装置。 - 【請求項１０】前記消失点演算手段が、「モーメント行
列の方法」により消失点を計算することを特徴とする、
請求項５又は６又は７記載のカメラの焦点距離較正装
置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5008997A JPH0799327B2 (ja) | 1993-01-22 | 1993-01-22 | カメラの焦点距離較正方法及び装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5008997A JPH0799327B2 (ja) | 1993-01-22 | 1993-01-22 | カメラの焦点距離較正方法及び装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH06221832A JPH06221832A (ja) | 1994-08-12 |
| JPH0799327B2 true JPH0799327B2 (ja) | 1995-10-25 |
Family
ID=11708333
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP5008997A Expired - Fee Related JPH0799327B2 (ja) | 1993-01-22 | 1993-01-22 | カメラの焦点距離較正方法及び装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0799327B2 (ja) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100711206B1 (ko) * | 2005-05-04 | 2007-04-24 | 한국과학기술연구원 | 체크보드 형태의 교차점 자동 추출 방법 |
| JP2010048553A (ja) * | 2006-12-19 | 2010-03-04 | Panasonic Corp | 複眼測距装置の検査方法およびそれに用いるチャート |
| US8102516B2 (en) | 2008-02-26 | 2012-01-24 | Panasonic Corporation | Test method for compound-eye distance measuring apparatus, test apparatus, and chart used for the same |
| WO2021113810A1 (en) * | 2019-12-06 | 2021-06-10 | Mashgin, Inc. | System and method for identifying items |
-
1993
- 1993-01-22 JP JP5008997A patent/JPH0799327B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH06221832A (ja) | 1994-08-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20200096317A1 (en) | Three-dimensional measurement apparatus, processing method, and non-transitory computer-readable storage medium | |
| KR100386090B1 (ko) | 동심원 패턴을 이용한 카메라 내부변수 보정시스템 및카메라 보정방법 | |
| CN102472609A (zh) | 位置和姿势校准方法及设备 | |
| JP2011179907A (ja) | 位置姿勢計測装置、位置姿勢計測方法およびプログラム | |
| KR20130138247A (ko) | 신속 3ｄ 모델링 | |
| WO2020031950A1 (ja) | 計測校正装置、計測校正方法、及びプログラム | |
| CN102289803A (zh) | 图像处理设备、图像处理方法及程序 | |
| CN107680139A (zh) | 一种远心双目立体视觉测量系统的通用性标定方法 | |
| WO2023273108A1 (zh) | 单目测距方法、装置及智能装置 | |
| CN104794718A (zh) | 一种单图像ct机房监控摄像机标定的方法 | |
| Svoboda et al. | Matching in catadioptric images with appropriate windows, and outliers removal | |
| CN111105467A (zh) | 一种图像标定方法、装置及电子设备 | |
| JP3696336B2 (ja) | カメラのキャリブレーション方法 | |
| CN113658279A (zh) | 相机内参和外参估算方法、装置、计算机设备及存储介质 | |
| JPH0799327B2 (ja) | カメラの焦点距離較正方法及び装置 | |
| CN116563370A (zh) | 基于单目计算机视觉的测距方法和测速方法 | |
| JP2005275789A (ja) | 三次元構造抽出方法 | |
| JP2011022084A (ja) | 三次元姿勢測定装置および三次元姿勢測定方法 | |
| JP3207023B2 (ja) | 画像計測装置の校正方法 | |
| JP2003331263A (ja) | カメラ画像−実空間座標対応付けマップ生成方法、カメラ画像−実空間座標対応付けマップ生成装置、プログラムおよび記録媒体 | |
| JPH04264207A (ja) | 多視点ステレオ画像計測方法 | |
| JPH076769B2 (ja) | 物体計測装置 | |
| JPH0835828A (ja) | ３次元計測装置のキャリブレーション方法 | |
| JP4332629B2 (ja) | 量子化誤差を軽減する距離計測方法及び距離計測システム | |
| CN117611680A (zh) | 光平面标定方法、装置和电子设备 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Effective date: 20031211 Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20031212 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081226 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091226 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 6 Free format text: PAYMENT UNTIL: 20091226 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101226 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 7 Free format text: PAYMENT UNTIL: 20101226 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 8 Free format text: PAYMENT UNTIL: 20111226 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 8 Free format text: PAYMENT UNTIL: 20111226 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 9 Free format text: PAYMENT UNTIL: 20121226 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131226 Year of fee payment: 10 |
|
| LAPS | Cancellation because of no payment of annual fees |