[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

JP7044016B2 - Line sensor camera calibration device and method - Google Patents

Line sensor camera calibration device and method Download PDF

Info

Publication number
JP7044016B2
JP7044016B2 JP2018166067A JP2018166067A JP7044016B2 JP 7044016 B2 JP7044016 B2 JP 7044016B2 JP 2018166067 A JP2018166067 A JP 2018166067A JP 2018166067 A JP2018166067 A JP 2018166067A JP 7044016 B2 JP7044016 B2 JP 7044016B2
Authority
JP
Japan
Prior art keywords
line sensor
sensor camera
dimensional
image
coordinates
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.)
Active
Application number
JP2018166067A
Other languages
Japanese (ja)
Other versions
JP2020038550A (en
Inventor
寛修 深井
勇介 渡部
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.)
Meidensha Corp
Original Assignee
Meidensha 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 Meidensha Corp filed Critical Meidensha Corp
Priority to JP2018166067A priority Critical patent/JP7044016B2/en
Publication of JP2020038550A publication Critical patent/JP2020038550A/en
Application granted granted Critical
Publication of JP7044016B2 publication Critical patent/JP7044016B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)
  • Studio Devices (AREA)
  • Image Analysis (AREA)

Description

本発明は、画像処理分野において、ラインセンサカメラのキャリブレーションを行う装置及び方法に関する。 The present invention relates to an apparatus and a method for calibrating a line sensor camera in the field of image processing.

下記非特許文献1では、カメラのキャリブレーション方法として、エリアカメラで平面マーカを様々な位置姿勢で撮影することにより、各撮影毎の外部パラメータ及びカメラの内部パラメータを求める提案を行っている。 In Non-Patent Document 1 below, as a camera calibration method, it is proposed to obtain an external parameter and an internal parameter of the camera for each shooting by shooting a plane marker in various positions and postures with an area camera.

また、下記特許文献1には、校正用パターンを計測対象物に投影し、計測対象物を撮影し、撮影画像を画像解析して特徴点を抽出し、取得した複数の特徴点に基づいて三次元計測に必要なパラメータのうちの少なくとも一つを推定するようにした校正装置、方法及びプログラムにおいて、全てのパラメータを未知数とすると不定性があり、値が一意に決まらないため、評価関数に正則化項をつけるなどして最適化を行うことが記載されている。 Further, in Patent Document 1 below, a calibration pattern is projected onto a measurement object, the measurement object is photographed, the photographed image is image-analyzed to extract feature points, and a tertiary feature point is extracted based on the acquired plurality of feature points. In a calibrator, method and program that estimates at least one of the parameters required for the original measurement, if all the parameters are unknown, there is indeterminacy and the values are not uniquely determined, so the evaluation function is regular. It is described that optimization is performed by adding a chemical term.

また、下記特許文献2には、対象画像から一部の画像領域を指定し、指定された画像領域からエッジを抽出し、エッジ位置の画像上の座標及び強度に基づいて歪み係数を算出するようにした撮影パラメータ測定方法及び装置並びに記録媒体において、各エッジにおける強度を、歪み係数を算出する処理における重み付けに用いることが記載されている。 Further, in Patent Document 2 below, a part of the image area is designated from the target image, an edge is extracted from the designated image area, and the distortion coefficient is calculated based on the coordinates and intensity of the edge position on the image. It is described that the intensity at each edge is used for weighting in the process of calculating the strain coefficient in the imaging parameter measuring method and apparatus and the recording medium.

また、下記特許文献3には、キャリブレーション用マーカを用いてラインセンサカメラの内部パラメータを校正するようにしたラインセンサカメラのキャリブレーション装置及び方法が記載されている。 Further, Patent Document 3 below describes a calibration device and method for a line sensor camera in which the internal parameters of the line sensor camera are calibrated using a calibration marker.

特開2015-36629号公報JP-A-2015-36629 特許第3661073号公報Japanese Patent No. 3661073 特開2016-218815号公報Japanese Unexamined Patent Publication No. 2016-218815

Z. Zhang. "A Flexible New Technique for Camera Calibration", IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol. 22, No. 11, 2000年, pp.1330-1334.Z. Zhang. "A Flexible New Technique for Camera Calibration", IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol. 22, No. 11, 2000, pp.1330-1334. Richard Szeliski,玉木 徹 他,“コンピュータビジョン‐アルゴリズムと応用‐”,共立出版,2013年, pp.51Richard Szeliski, Toru Tamaki et al., "Computer Vision-Algorithms and Applications-", Kyoritsu Shuppan, 2013, pp.51 C.M.ビショップ著,「パターン認識と機械学習 上」,シュプリンガー・ジャパン,2007年,pp.9C.M. Bishop, "Pattern Recognition and Machine Learning," Springer Japan, 2007, pp.9

非特許文献1に記載のキャリブレーション方法は、平面マーカを様々な位置姿勢で撮影し、内部パラメータと各撮影の外部パラメータを非線形最適化手法によって解く方法である。平面マーカを用いるためマーカの作成が容易なこと、複数回の撮影を行うため焦点距離などを高精度に推定できることといったメリットがある。これらの特徴から、非特許文献1の方法はカメラのキャリブレーション方法として現在最も用いられている。しかしながら、この方法はラインセンサカメラには対応していなかった。例えば、ラインセンサカメラにおいて、非特許文献1の方法のように平面マーカを用いると、焦点距離などの算出が不安定になるという問題があった。 The calibration method described in Non-Patent Document 1 is a method in which plane markers are photographed in various positions and orientations, and internal parameters and external parameters of each imaging are solved by a nonlinear optimization method. Since the planar marker is used, it is easy to create a marker, and since multiple shootings are performed, the focal length and the like can be estimated with high accuracy. Due to these characteristics, the method of Non-Patent Document 1 is currently most used as a camera calibration method. However, this method has not been compatible with line sensor cameras. For example, in a line sensor camera, if a plane marker is used as in the method of Non-Patent Document 1, there is a problem that the calculation of the focal length and the like becomes unstable.

また、特許文献1に記載のものは、評価関数に正則化項を付けるなどして最適化を行うことで安定的にパラメータを求める方法であるが、各パラメータをどの程度正則化項として重要視するかといった事前知識については特に考慮されていなかった。 Further, the method described in Patent Document 1 is a method of stably obtaining parameters by performing optimization such as adding a regularization term to the evaluation function, but how important each parameter is as a regularization term. No particular consideration was given to prior knowledge such as whether to do so.

また、特許文献2に記載のものは、評価関数に重み付けを取り入れる点で正則化に重み付けをする本発明の考え方と類似しているが、本発明とは、重み付けをする箇所および考え方が異なる。また、パラメータの値を直接的に安定的に推定する機能がないという問題があった。 Further, what is described in Patent Document 2 is similar to the idea of the present invention in which the regularization is weighted in that the weighting is incorporated into the evaluation function, but the place and the idea of weighting are different from the present invention. In addition, there is a problem that there is no function to directly and stably estimate the value of the parameter.

また、特許文献3に記載のものは、エリアカメラに比べ情報量の少ないラインセンサカメラでマーカ画像を複数枚撮影するだけでキャリブレーションが可能な装置だが、歪み係数が極端に大きくなったり焦点距離が実際の値と大きく異なったりということが撮影画像によっては起こり得るという問題があった。 Further, the device described in Patent Document 3 is a device capable of calibrating by simply taking a plurality of marker images with a line sensor camera having a smaller amount of information than an area camera, but the distortion coefficient becomes extremely large and the focal length becomes extremely large. However, there is a problem that it may occur depending on the captured image that the value is significantly different from the actual value.

このようなことから本発明は、ラインセンサカメラのキャリブレーションを安定して行うことが可能なラインセンサカメラのキャリブレーション装置及び方法を提供することを目的とする。 Therefore, it is an object of the present invention to provide a calibration device and method for a line sensor camera capable of stably calibrating the line sensor camera.

上記課題を解決する第1の発明に係るラインセンサカメラのキャリブレーション装置は、
ラインセンサカメラの撮像範囲となる2次元平面上での2次元平面座標が既知の複数の検出点を有する立体マーカと、
前記立体マーカを前記ラインセンサカメラで撮像した複数の画像に基づいて、前記ラインセンサカメラのキャリブレーション用の内部パラメータ及び外部パラメータを算出する処理部とを有し、
前記処理部は、
位置及び姿勢が各々異なる前記立体マーカの前記画像が各々入力される画像入力部と、
前記立体マーカの前記検出点の前記画像上の1次元画像座標を前記画像毎に算出する計測座標算出部と、
前記2次元平面座標から前記1次元画像座標へ変換する前記ラインセンサカメラの数式モデルに基づいて求めた変換行列について、前記変換行列における初期の内部パラメータとして既知の初期値を入力し、線形解法により前記変換行列から初期の外部パラメータを前記画像毎に算出する外部パラメータ算出部と、
前記数式モデルに基づいて求めたレンズ歪みを考慮した誤差関数に、前記内部パラメータを考慮した値を加算してなる評価関数について、前記評価関数に前記初期の内部パラメータ及び前記初期の外部パラメータを入力し、前記検出点の全数及び前記画像の全数についての前記誤差関数の総和に前記内部パラメータを考慮した値を加算した評価関数が最小となる前記内部パラメータ及び前記外部パラメータを、反復計算を用いた非線形解法により算出する非線形最小二乗部とを有する
ことを特徴とする。
The calibration device for the line sensor camera according to the first invention, which solves the above problems, is
A three-dimensional marker having a plurality of detection points whose two-dimensional plane coordinates on the two-dimensional plane which is the imaging range of the line sensor camera are known, and
It has a processing unit for calculating internal parameters and external parameters for calibration of the line sensor camera based on a plurality of images of the three-dimensional marker captured by the line sensor camera.
The processing unit
An image input unit in which the images of the three-dimensional markers having different positions and postures are input, and
A measurement coordinate calculation unit that calculates the one-dimensional image coordinates of the detection point of the three-dimensional marker on the image for each image, and the measurement coordinate calculation unit.
For the transformation matrix obtained based on the mathematical model of the line sensor camera that converts the two-dimensional plane coordinates to the one-dimensional image coordinates, input known initial values as the initial internal parameters in the transformation matrix, and use a linear solution method. An external parameter calculation unit that calculates initial external parameters from the transformation matrix for each image,
For an evaluation function obtained by adding a value considering the internal parameters to an error function considering the lens distortion obtained based on the mathematical model, the initial internal parameters and the initial external parameters are input to the evaluation function. Then, the internal parameter and the external parameter, which minimizes the evaluation function obtained by adding the value considering the internal parameter to the sum of the error functions for the total number of the detection points and the total number of the images, were repeatedly calculated. It is characterized by having a non-linear minimum squared part calculated by a non-linear solution method.

上記課題を解決する第2の発明に係るラインセンサカメラのキャリブレーション装置は、
上記第1の発明に係るラインセンサカメラのキャリブレーション装置において、
前記評価関数は、レンズの焦点距離の大きさ、レンズの主点座標の位置及びレンズの歪み係数の大きさにそれぞれ重み付けをして前記誤差関数に加算したものである
ことを特徴とする。
The calibration device for the line sensor camera according to the second invention, which solves the above problems, is
In the calibration device for the line sensor camera according to the first invention.
The evaluation function is characterized in that the magnitude of the focal length of the lens, the position of the principal point coordinates of the lens, and the magnitude of the distortion coefficient of the lens are each weighted and added to the error function.

上記課題を解決する第3の発明に係るラインセンサカメラのキャリブレーション装置は、
上記第1または第2の発明に係るラインセンサカメラのキャリブレーション装置において、
前記立体マーカは、L字型の形状であり、L字型に折れ曲がった内側の面に前記複数の検出点が配置されている
ことを特徴とする。
The calibration device for the line sensor camera according to the third invention, which solves the above problems, is
In the calibration device for the line sensor camera according to the first or second invention.
The three-dimensional marker has an L-shaped shape, and is characterized in that the plurality of detection points are arranged on an inner surface that is bent into an L-shape.

上記課題を解決する第4の発明に係るラインセンサカメラのキャリブレーション装置は、
上記第1から第3のいずれか一つの発明に係るラインセンサカメラのキャリブレーション装置において、
前記数式モデルは、ピンホールカメラモデルを用いたものである
ことを特徴とする。
The calibration device for the line sensor camera according to the fourth invention, which solves the above problems, is
In the calibration device for the line sensor camera according to any one of the first to third inventions.
The mathematical model is characterized by using a pinhole camera model.

上記課題を解決する第5の発明に係るラインセンサカメラのキャリブレーション方法は、
ラインセンサカメラの撮像範囲となる2次元平面上での2次元平面座標が既知の複数の検出点を有する立体マーカを前記ラインセンサカメラで撮像した複数の画像に基づいて、前記ラインセンサカメラのキャリブレーション用の内部パラメータ及び外部パラメータを算出するラインセンサカメラのキャリブレーション方法であって、
位置及び姿勢が各々異なる前記立体マーカの前記画像が各々入力される画像入力工程と、
前記立体マーカの前記検出点の前記画像上の1次元画像座標を前記画像毎に算出する計測座標算出工程と、
前記2次元平面座標から前記1次元画像座標へ変換する前記ラインセンサカメラの数式モデルに基づいて求めた変換行列について、前記変換行列における初期の内部パラメータとして既知の初期値を入力し、線形解法により前記変換行列から初期の外部パラメータを前記画像毎に算出する外部パラメータ算出工程と、
前記数式モデルに基づいて求めたレンズ歪みを考慮した誤差関数に、前記内部パラメータを考慮した値を加算してなる評価関数について、前記評価関数に前記初期の内部パラメータ及び前記初期の外部パラメータを入力し、前記検出点の全数及び前記画像の全数についての前記誤差関数の総和に前記内部パラメータを考慮した値を加算した評価関数が最小となる前記内部パラメータ及び前記外部パラメータを、反復計算を用いた非線形解法により算出する非線形最小二乗工程とを有する
ことを特徴とする。
The method for calibrating the line sensor camera according to the fifth invention, which solves the above problems, is
Calibration of the line sensor camera based on a plurality of images taken by the line sensor camera of a three-dimensional marker having a plurality of detection points whose two-dimensional plane coordinates are known on the two-dimensional plane which is the imaging range of the line sensor camera. A line sensor camera calibration method that calculates internal and external parameters for the camera.
An image input step in which the images of the three-dimensional markers having different positions and postures are input, and
A measurement coordinate calculation step of calculating the one-dimensional image coordinates of the detection point of the three-dimensional marker on the image for each image, and
For the transformation matrix obtained based on the mathematical model of the line sensor camera that converts the two-dimensional plane coordinates to the one-dimensional image coordinates, input known initial values as the initial internal parameters in the transformation matrix, and use a linear solution method. An external parameter calculation step of calculating the initial external parameters from the transformation matrix for each image, and
For an evaluation function obtained by adding a value considering the internal parameters to an error function considering the lens distortion obtained based on the mathematical model, the initial internal parameters and the initial external parameters are input to the evaluation function. Then, the internal parameter and the external parameter, which minimizes the evaluation function obtained by adding the value considering the internal parameter to the sum of the error functions for the total number of the detection points and the total number of the images, were repeatedly calculated. It is characterized by having a non-linear minimum square step calculated by a non-linear solution method.

上記課題を解決する第6の発明に係るラインセンサカメラのキャリブレーション方法は、
上記第5の発明に係るラインセンサカメラのキャリブレーション方法において、
前記評価関数として、レンズの焦点距離の大きさ、レンズの主点座標の位置及びレンズの歪み係数の大きさにそれぞれ重み付けをして前記誤差関数に加算したものを用いる
ことを特徴とする。
The method for calibrating the line sensor camera according to the sixth invention for solving the above problems is as follows.
In the calibration method of the line sensor camera according to the fifth aspect of the invention.
The evaluation function is characterized in that the size of the focal length of the lens, the position of the principal point coordinates of the lens, and the size of the distortion coefficient of the lens are each weighted and added to the error function.

上記課題を解決する第7の発明に係るラインセンサカメラのキャリブレーション方法は、
上記第5または第6の発明に係るラインセンサカメラのキャリブレーション方法において、
前記立体マーカとして、L字型の形状であって、L字型に折れ曲がった内側の面に前記複数の検出点が配置されているものを用いる
ことを特徴とする。
The method for calibrating the line sensor camera according to the seventh invention, which solves the above problems, is
In the calibration method of the line sensor camera according to the fifth or sixth invention.
The three-dimensional marker is characterized by using an L-shaped marker in which the plurality of detection points are arranged on an inner surface bent into an L-shape.

上記課題を解決する第8の発明に係るラインセンサカメラのキャリブレーション方法は、
上記第5から第7のいずれか一つの発明に係るラインセンサカメラのキャリブレーション方法において、
前記数式モデルとして、ピンホールカメラモデルを用いる
ことを特徴とする。
The method for calibrating the line sensor camera according to the eighth invention for solving the above problems is as follows.
In the calibration method of the line sensor camera according to any one of the fifth to seventh inventions.
A pinhole camera model is used as the mathematical model.

本発明に係るラインセンサカメラのキャリブレーション装置及び方法によれば、ラインセンサカメラのキャリブレーションを安定して行うことができる。 According to the calibration device and method for the line sensor camera according to the present invention, the line sensor camera can be calibrated in a stable manner.

本発明に係るラインセンサカメラのキャリブレーション装置の実施形態の一例を示す概略構成図である。It is a schematic block diagram which shows an example of embodiment of the calibration apparatus of the line sensor camera which concerns on this invention. 図1に示した装置で実施するラインセンサカメラのキャリブレーション方法を説明するフローチャートである。It is a flowchart explaining the calibration method of the line sensor camera carried out by the apparatus shown in FIG. L字マーカを撮影して取得した画像の例である。This is an example of an image obtained by photographing an L-shaped marker. L字マーカを撮影して取得した画像の他の例である。This is another example of an image obtained by photographing an L-shaped marker.

カメラを用いた計測を行う際、カメラの焦点距離、主点座標、歪み係数などといったカメラパラメータを事前に求めることは必須であり、これを求めることをカメラキャリブレーションという。なお、以降、カメラキャリブレーションを、単に、キャリブレーションと呼ぶ。 When making measurements using a camera, it is essential to obtain camera parameters such as the focal length, principal point coordinates, and distortion coefficient of the camera in advance, and obtaining this is called camera calibration. Hereinafter, camera calibration is simply referred to as calibration.

本発明では、非特許文献1の方法のように、ラインセンサカメラのキャリブレーションに、マーカの複数回撮影と非線形最適化手法によるカメラパラメータ算出を導入する。しかしながら、非特許文献1の方法はエリアカメラを前提としているため、この方法をそのままラインセンサカメラのキャリブレーションに導入することはできない。例えば、ラインセンサカメラでは、次元数が少ないため、回転行列の直交性を活用できず、又、ラインセンサカメラがキャリブレーション用の平面マーカをチェスボード状に撮影できないため、焦点距離と外部パラメータ算出が不安定になることが問題となる。そこで、本発明では、キャリブレーション用として、立体マーカを用いている。 In the present invention, as in the method of Non-Patent Document 1, the camera parameter calculation by the multiple shooting of the marker and the nonlinear optimization method is introduced into the calibration of the line sensor camera. However, since the method of Non-Patent Document 1 is premised on an area camera, this method cannot be directly introduced into the calibration of a line sensor camera. For example, in a line sensor camera, the number of dimensions is small, so the orthogonality of the rotation matrix cannot be utilized, and since the line sensor camera cannot take a chessboard-like image of a plane marker for calibration, the focal length and external parameters are calculated. The problem is that it becomes unstable. Therefore, in the present invention, a three-dimensional marker is used for calibration.

又、本発明では、カメラモデルとして、ピンホールカメラモデルを用いる。但し、本発明では、ラインセンサカメラを使うため、通常のエリアカメラでは撮像範囲の3次元空間が2次元画像になるのに対して、ラインセンサカメラでは撮像範囲の2次元平面が1次元画像となる。そして、撮像範囲の2次元平面上に立体マーカを設置することで、求める外部パラメータを、回転が1軸、平行移動が2次元とすることができる。 Further, in the present invention, a pinhole camera model is used as the camera model. However, in the present invention, since the line sensor camera is used, the three-dimensional space of the imaging range becomes a two-dimensional image in the normal area camera, whereas the two-dimensional plane of the imaging range becomes the one-dimensional image in the line sensor camera. Become. Then, by installing the three-dimensional marker on the two-dimensional plane of the imaging range, the external parameters to be obtained can be set to one axis for rotation and two dimensions for translation.

ここで、本発明の装置構成を図1に示すと共に、図1に示す構成における数式モデルを下記式(1)に示す。なお、図1に示す本発明の装置構成の詳細については後述する。 Here, the apparatus configuration of the present invention is shown in FIG. 1, and the mathematical model in the configuration shown in FIG. 1 is shown in the following equation (1). The details of the apparatus configuration of the present invention shown in FIG. 1 will be described later.

又、下記の式(1)~(17)の式中において、ボールド体(太字)ではないアルファベットはスカラーであり、ボールド体(太字)の小文字のアルファベットはベクトルであり、ボールド体(太字)の大文字のアルファベットは行列である。なお、文書中においては、ボールド体(太字)を使用できないので、通常の書式とし、該当するアルファベットの直後に“(ボールド体)”と付記している。 Further, in the following equations (1) to (17), the alphabet that is not bold (bold) is a scalar, and the lowercase alphabet of bold (bold) is a vector, and the bold alphabet (bold). The uppercase alphabet is a matrix. In addition, since bold type (bold type) cannot be used in the document, the normal format is used and "(bold type)" is added immediately after the corresponding alphabet.

図1において、ラインセンサカメラ10は1ラインしか撮像を行わないため、撮像範囲は3次元空間上では2次元平面11となる。この2次元平面11上のL字マーカ20のオブジェクト座標系(2次元平面座標系)をX軸、Y軸で表現し、ラインセンサカメラ10のカメラ座標系(1次元画像座標系)を、ラインセンサカメラ10の光軸方向の軸wに直交し、撮像範囲の2次元平面11と並行なu軸で表現している。又、オブジェクト座標系とカメラ座標系間の回転行列をR(ボールド体)(回転軸は1軸でその角度はθ)、並進ベクトルをt(ボールド体)(2次元ベクトルで、その要素はt1、t2)とする。 In FIG. 1, since the line sensor camera 10 captures only one line, the imaging range is the two-dimensional plane 11 in the three-dimensional space. The object coordinate system (two-dimensional plane coordinate system) of the L-shaped marker 20 on the two-dimensional plane 11 is represented by the X-axis and the Y-axis, and the camera coordinate system (one-dimensional image coordinate system) of the line sensor camera 10 is represented by a line. It is represented by a u-axis orthogonal to the axis w in the optical axis direction of the sensor camera 10 and parallel to the two-dimensional plane 11 of the imaging range. The rotation matrix between the object coordinate system and the camera coordinate system is R (bold body) (the axis of rotation is one axis and its angle is θ), and the translation vector is t (bold body) (two-dimensional vector whose elements are t). 1 , t 2 ).

Figure 0007044016000001
Figure 0007044016000001

ここで、sはスケーリング係数、fは焦点距離、cは主点座標、uは実際に観測されたカメラ座標系における位置(ピクセル座標における計測座標)、X、Yはオブジェクト座標系における位置(L字マーカ20の実空間(2次元平面11)上における真値座標)である。 Here, s is a scaling coefficient, f is a focal distance, c is a principal point coordinate, u is a position in an actually observed camera coordinate system (measured coordinates in pixel coordinates), and X and Y are positions in an object coordinate system (L). It is the true value coordinates in the real space (two-dimensional plane 11) of the character marker 20.

上記式(1)は、L字マーカ20のオブジェクト座標(2次元平面座標)からラインセンサカメラ10のカメラ座標(1次元画像座標)へ変換する数式モデルとなっている。但し、上記式(1)では、ラインセンサカメラ10のレンズの歪みが考慮されていない。そこで、下記式(2)のように(x,y)を定義すると、下記式(3)~(5)のように歪みを表すことができる。なお、k1~k3は半径歪み係数である。また、レンズの歪みについては、非特許文献2(特に、51頁)を参考にしている。 The above equation (1) is a mathematical model for converting the object coordinates (two-dimensional plane coordinates) of the L-shaped marker 20 into the camera coordinates (one-dimensional image coordinates) of the line sensor camera 10. However, in the above equation (1), the distortion of the lens of the line sensor camera 10 is not taken into consideration. Therefore, if (x, y) is defined as in the following equation (2), the strain can be expressed as in the following equations (3) to (5). Note that k 1 to k 3 are radial distortion coefficients. Further, regarding the distortion of the lens, Non-Patent Document 2 (particularly, page 51) is referred to.

Figure 0007044016000002
Figure 0007044016000002

上記式(1)~(5)について展開したいが、式(5)にx”をそのまま代入できないため、次式(5)’を式(4)に代入する。式(1)~(3)についても式(4)に代入すると下記式(6)となる。 I would like to expand the above equations (1) to (5), but since x "cannot be substituted into equation (5) as it is, the following equation (5)'is substituted into equation (4). Equations (1) to (3) Is also substituted into the equation (4) to obtain the following equation (6).

Figure 0007044016000003
Figure 0007044016000003

このカメラモデルでキャリブレーションを行う場合には、単位焦点距離面での誤差最小化を行うことになり、下記式(7)が最小となるように、レーベンバーグマーカート法により、内部及び外部パラメータを求めることとなる。 When calibrating with this camera model, the error in the unit focal length plane will be minimized, and the internal and external parameters will be minimized by the Levenberg-Marquardt method so that the following equation (7) is minimized. Will be sought.

Figure 0007044016000004
Figure 0007044016000004

上記式(7)において、Mは画像の撮影回数、Nは各撮影において取得される特徴点数(白黒帯23の検出点数)である。上記式(7)は、上記式(6)に基づく誤差関数を、撮像した画像の全数M及び検出点の全数Nについて総和したものとなる。 In the above equation (7), M is the number of times the image is taken, and N is the number of feature points (the number of detected points of the black and white band 23) acquired in each picture. The above equation (7) is the sum of the error functions based on the above equation (6) for the total number M of the captured images and the total number N of the detection points.

ここで、立体マーカを用いて実際にキャリブレーションを行ってみると立体マーカを用いるだけでは解が不安定になることを完全には防げないことが分かった。ここで不安定というのは「焦点距離はレンズに定められた初期値から大きくは変化しないはずだ」、「レンズの歪み係数k1,k2,k3は極端に大きな値にはならず、またレンズの歪み係数k1,k2,k3は一般的にk1>k2>k3となることが多い」というカメラキャリブレーションを行う熟練者であれば知っている一般的な知識から大きく外れることを指す。 Here, when the calibration was actually performed using the steric marker, it was found that the solution could not be completely prevented from becoming unstable just by using the steric marker. Unstable here means that "the focal length should not change significantly from the initial value set for the lens", "the distortion coefficients k 1 , k 2 , k 3 of the lens do not become extremely large values, and In addition, the distortion coefficients k 1 , k 2 , and k 3 of the lens are generally k 1 > k 2 > k 3 "from the general knowledge known to those who are experts in camera calibration. It means that it deviates greatly.

すなわち、上記式(7)が最小になるパラメータを最適化処理で求めると、パラメータが非現実的な値になる可能性があった。そこで本実施例では、上記式(6)に基づく誤差関数に、「焦点距離fはレンズの規定値f0から大きくは変わらない」、「主点座標cは(カメラによるが)画像中心c0の近くである」、「レンズの歪み係数k1,k2,k3は極端に大きな数字にはならない」、「レンズの歪み係数の大きさはk1>k2>k3となることが多い」等の専門家の知識を導入した評価関数(下記式(8))が最小となるように、レーベンバーグマーカート法により内部及び外部パラメータを求めることとした。 That is, if the parameter that minimizes the above equation (7) is obtained by the optimization process, the parameter may have an unrealistic value. Therefore, in this embodiment, the error function based on the above equation (6) is that "the focal length f does not change significantly from the specified value f 0 of the lens" and "the principal point coordinate c is the image center c 0 (depending on the camera)". "The distortion coefficient of the lens k 1 , k 2 , k 3 is not an extremely large number", "The magnitude of the distortion coefficient of the lens is k 1 > k 2 > k 3 " It was decided to obtain the internal and external parameters by the Lebenberg Marcart method so that the evaluation function (the following equation (8)) that introduced the knowledge of experts such as "many" is minimized.

Figure 0007044016000005
Figure 0007044016000005

上記式(8)において、wfは焦点距離の正則化項の重み係数、wcは主点座標の正則化項の重み係数、c0は主点座標の既定値、wk1は歪み係数k1の正則化項の重み係数、wk2は歪み係数k2の正則化項の重み係数、wk3は歪み係数k3の正則化項の重み係数である。 In the above equation (8), w f is the weighting coefficient of the regularization term of the focal distance, w c is the weighting coefficient of the regularization term of the main point coordinates, c 0 is the default value of the main point coordinates, and w k 1 is the distortion coefficient k. 1 is the weighting coefficient of the regularized term, w k2 is the weighting coefficient of the regularized term of the strain coefficient k 2 , and w k 3 is the weighting coefficient of the regularized term of the strain coefficient k 3 .

上記式(8)を用いることで、誤差が小さいだけでなく、焦点距離の大きさ、主点座標の位置、歪み係数の大きさ等の内部パラメータを考慮したキャリブレーションを行うことができる。また各内部パラメータ(ここでは、焦点距離の大きさ、主点座標の位置、歪み係数の大きさ)に重み付けをすることで、焦点距離の大きさ、主点座標の位置、歪み係数の大きさのうち、どのパラメータを重要視するかを必要に応じて設定することができ、これにより、設計者や専門家の知識を導入することができる。なお、上記式(8)ではL2ノルムと呼ばれる二乗の指標を用いたが、L1ノルムと呼ばれる絶対値指標を使用することも有効である。 By using the above equation (8), not only the error is small, but also calibration can be performed in consideration of internal parameters such as the size of the focal length, the position of the principal point coordinates, and the size of the strain coefficient. In addition, by weighting each internal parameter (here, the size of the focal length, the position of the principal point coordinates, the size of the strain coefficient), the size of the focal length, the position of the principal point coordinates, and the size of the strain coefficient are weighted. Of these, which parameter is important can be set as needed, which can introduce the knowledge of designers and experts. Although the square index called the L2 norm is used in the above equation (8), it is also effective to use the absolute value index called the L1 norm.

次に、上述した内部及び外部パラメータが全て校正済み(キャリブレーション済み)の場合に、計測座標uが得られた時のX、Yについて求める。これは、上記式(6)をX、Yの形になるように解けばよいが、上記式(6)にはX、Y両方が含まれているため、このままでは解くことができない。そこで、L字マーカ20の検出点では、X、Yの2次元平面11上の真値座標がどちらかは0であることを利用し、真値座標が0となる方向で得られたパラメータによって算出される値も0と仮定して解く。その結果、上記式(6)は下記式(9)、(10)となり、解析的にX、Yの値を求めることができた。ここで、x’やx”は上記式(3)~(5)で定義された値である。 Next, when all the above-mentioned internal and external parameters have been calibrated (calibrated), X and Y when the measurement coordinates u are obtained are obtained. This may be solved by solving the above equation (6) in the form of X and Y, but since the above equation (6) includes both X and Y, it cannot be solved as it is. Therefore, at the detection point of the L-shaped marker 20, the fact that either of the true value coordinates on the two-dimensional plane 11 of X and Y is 0 is utilized, and the parameters obtained in the direction in which the true value coordinates become 0 are used. The calculated value is also assumed to be 0 and solved. As a result, the above equation (6) became the following equations (9) and (10), and the values of X and Y could be obtained analytically. Here, x'and x'are the values defined by the above equations (3) to (5).

Figure 0007044016000006
Figure 0007044016000006

次に、外部パラメータを撮影毎に求める方法について述べる。並進と回転を表す変換行列を下記式(11)とすると、上記式(1)に基づいて、変換行列は下記式(12)で表すことができる。 Next, a method of obtaining external parameters for each shooting will be described. Assuming that the transformation matrix representing translation and rotation is the following equation (11), the transformation matrix can be expressed by the following equation (12) based on the above equation (1).

Figure 0007044016000007
Figure 0007044016000007

これを展開、整理すると、下記式(13)となる。 When this is expanded and organized, the following equation (13) is obtained.

Figure 0007044016000008
Figure 0007044016000008

ここで、スケールが次数を1つ減らせる性質を利用し、3本以上の同時方程式を立てることで変換行列を解くことができる。実際には、L字マーカ20の座標は3点以上あるが、この場合、方程式と未知数の数が一致しない。そこで、最小二乗法を用いる。二乗誤差をE、各計測点で生じる誤差ベクトルをe(ボールド体)、マーカの座標数をnとし、下記式(14)で表す行列を用いると、二乗誤差Eを下記式(15)で表すことができる。なお、“T”は、転置行列を表す記号である。 Here, the transformation matrix can be solved by constructing three or more simultaneous equations by utilizing the property that the scale can reduce the order by one. Actually, the coordinates of the L-shaped marker 20 are three or more points, but in this case, the equation and the number of unknowns do not match. Therefore, the least squares method is used. The squared error is E, the error vector generated at each measurement point is e (bold body), the number of coordinates of the marker is n, and the squared error E is expressed by the following equation (15) using the matrix expressed by the following equation (14). be able to. Note that " T " is a symbol representing a transposed matrix.

Figure 0007044016000009
Figure 0007044016000009

二乗誤差Eを最小化するためのベクトルは、BTB(ボールド体)の最小の固有値に対応する固有ベクトルとなる。そこで、解となる固有ベクトルを求めるため、特異値分解を用いる。これは下記式(16)となる。ここで、U(ボールド体)は左特異行列、V(ボールド体)は右特異行列、Σ(ボールド体)は特異値行列であり、I(ボールド体)は単位行列である。これにより、BTB(ボールド体)は下記式(17)で表すことができる。 The vector for minimizing the square error E is the eigenvector corresponding to the minimum eigenvalue of BTB (bold body). Therefore, singular value decomposition is used to obtain the eigenvector to be the solution. This is given by the following equation (16). Here, U (bold form) is a left singular matrix, V (bold form) is a right singular matrix, Σ (bold form) is a singular value matrix, and I (bold form) is an identity matrix. As a result, the BTB (bold body) can be expressed by the following equation (17).

Figure 0007044016000010
Figure 0007044016000010

この結果より、固有ベクトルが右特異行列V(ボールド体)の右特異ベクトルと同じになり、固有値が特異値の二乗と同じになることが分かる。このことから、特異値分解を使うことで、変換行列を求めることができる。 From this result, it can be seen that the eigenvector is the same as the right singular vector of the right singular matrix V (bold), and the eigenvalue is the same as the square of the singular value. From this, the transformation matrix can be obtained by using the singular value decomposition.

実際の外部パラメータとするには、スケール不定及びt2>0のため、t2で正規化を行い、回転行列の性質より、(cosθ)2+(sinθ)2=1を使用して、回転と並進のパラメータを求める。 To make it an actual external parameter, the scale is indefinite and t 2 > 0, so normalization is performed at t 2 , and due to the nature of the rotation matrix, (cos θ) 2 + (sin θ) 2 = 1 is used for rotation. And find the translation parameters.

このように、ラインセンサカメラであっても、計測座標u、真値座標X、Y、焦点距離f、主点座標cが分かっていれば、線形解法で解を求めることができる。実際には、焦点距離f、主点座標cは最初には分からないため、適当な値を入力することで、初期の外部パラメータを求める。そして、最終的に、この初期の外部パラメータを初期値として、式(8)をレーベンバーグマーカート法によって解くことで、最終的な内部及び外部パラメータを得ることができる。 As described above, even with the line sensor camera, if the measurement coordinates u, the true value coordinates X and Y, the focal length f, and the principal point coordinates c are known, the solution can be obtained by the linear solution method. Actually, since the focal length f and the principal point coordinates c are not known at first, the initial external parameters are obtained by inputting appropriate values. Finally, the final internal and external parameters can be obtained by solving the equation (8) by the Levenberg-Marquart method with this initial external parameter as the initial value.

次に、上述した計算を行う本発明に係るラインセンサカメラのキャリブレーション装置及び方法について、図1及び図2を参照して以下に説明する。 Next, the calibration device and method of the line sensor camera according to the present invention for performing the above calculation will be described below with reference to FIGS. 1 and 2.

[実施例1]
本実施例の装置は、ラインセンサカメラ10のキャリブレーション装置であり、図1に示すように、L字マーカ20と、処理システム30とを有し、処理システム30は、記憶部31、画像入力部32、計測座標算出部33、外部パラメータ算出部34及び非線形最小二乗部35を有している。
[Example 1]
The device of this embodiment is a calibration device of the line sensor camera 10, and as shown in FIG. 1, has an L-shaped marker 20 and a processing system 30, and the processing system 30 has a storage unit 31 and an image input. It has a unit 32, a measurement coordinate calculation unit 33, an external parameter calculation unit 34, and a non-linear minimum square unit 35.

ラインセンサカメラ10は、CCD(Charge-Coupled Device)などの撮像素子が一次元の線状に配置され、1次元の線状に撮像対象物の画像を撮像するものである。 In the line sensor camera 10, an image pickup device such as a CCD (Charge-Coupled Device) is arranged in a one-dimensional linear shape, and an image of an object to be imaged is captured in the one-dimensional linear shape.

L字マーカ20は、L字型の形状の立体マーカであり、その中央21で90°に折れ曲がり、折れ曲がった内側の面22に白黒帯23を有する剛体からなる。この白黒帯23は、中央21を中心に対称的に配置されると共に、内側の面22の長手方向に沿って交互に白黒となるように配置されており、白黒帯23の白黒の幅が既知、つまり、白黒帯23による複数の検出点が2次元平面11上における2次元平面座標で既知となっている。上述した非特許文献1の方法のように、平面マーカを用いる場合には焦点距離が不安定になるが、本実施例では、L字マーカ20をキャリブレーションに用いており、焦点距離を正確に算出可能となる。 The L-shaped marker 20 is a three-dimensional marker having an L-shaped shape, and is formed of a rigid body that is bent at 90 ° at the center 21 and has a black-and-white band 23 on the bent inner surface 22. The black-and-white band 23 is arranged symmetrically about the center 21 and alternately arranged in black and white along the longitudinal direction of the inner surface 22, and the black-and-white width of the black-and-white band 23 is known. That is, a plurality of detection points by the black-and-white band 23 are known in the two-dimensional plane coordinates on the two-dimensional plane 11. When a plane marker is used as in the method of Non-Patent Document 1 described above, the focal length becomes unstable, but in this embodiment, the L-shaped marker 20 is used for calibration, and the focal length is accurately set. It can be calculated.

なお、本実施例では、L字型の立体マーカを用いているが、複数の検出点に奥行きを持たせることができ、かつ、全ての検出点の2次元平面11上のX、Yの真値座標のどちらかを必ず0とすることができれば、L字型の立体マーカに限らない。例えば、直方体や直角三角柱において、直角を挟む2辺の外面に図1に示したような白黒帯23を設けても良い。 Although an L-shaped three-dimensional marker is used in this embodiment, it is possible to give depth to a plurality of detection points, and the true X and Y on the two-dimensional plane 11 of all the detection points are true. If either of the value coordinates can be set to 0, the marker is not limited to the L-shaped solid marker. For example, in a rectangular parallelepiped or a right-angled triangular prism, a black-and-white band 23 as shown in FIG. 1 may be provided on the outer surfaces of two sides sandwiching the right angle.

処理システム30において、記憶部31では、後述するように、ラインセンサカメラ10からの画像データ、各撮影でのL字マーカ20の白黒帯23の画像上の計測座標データ(計測座標u)、レンズ初期値、L字マーカ20の白黒帯23の真値座標X、Y、ラインセンサカメラ10の内部パラメータ(焦点距離f、主点座標c、歪み係数k1~k3)及び各撮影でのラインセンサカメラ10の外部パラメータ(回転θ、並進t1、t2)を保管する。記憶部31としては、RAM(Random Access memory)、HDD(Hard Disk Drive)、SSD(Solid State Drive)などの記憶装置が使用可能である。 In the processing system 30, in the storage unit 31, as will be described later, image data from the line sensor camera 10, measurement coordinate data (measurement coordinates u) on the image of the black and white band 23 of the L-shaped marker 20 in each shooting, and a lens. Initial values, true value coordinates X and Y of the black and white band 23 of the L-shaped marker 20, internal parameters of the line sensor camera 10 (focal length f, principal point coordinates c, distortion coefficient k 1 to k 3 ), and lines in each shooting. The external parameters (rotation θ, translation t 1 , t 2 ) of the sensor camera 10 are stored. As the storage unit 31, a storage device such as a RAM (Random Access memory), an HDD (Hard Disk Drive), or an SSD (Solid State Drive) can be used.

又、画像入力部32では、ラインセンサカメラ10で取得したL字マーカ20の画像データが、ラインセンサカメラ10から入力されて、入力された画像データを記憶部31へ保管している。 Further, in the image input unit 32, the image data of the L-shaped marker 20 acquired by the line sensor camera 10 is input from the line sensor camera 10, and the input image data is stored in the storage unit 31.

又、計測座標算出部33では、ラインセンサカメラ10で撮像した画像データからL字マーカ20の白黒帯23(検出点)を検出し、各撮影でのL字マーカ20の白黒帯23(検出点)の画像上の計測座標u(1次元画像座標)を求め、計測座標データとして記憶部31へ保管する。なお、L字マーカ20の白黒帯23(検出点)の画像上の座標については、2値化処理やエッジ検出といった画像処理を行うことにより検出している。 Further, the measurement coordinate calculation unit 33 detects the black-and-white band 23 (detection point) of the L-shaped marker 20 from the image data captured by the line sensor camera 10, and the black-and-white band 23 (detection point) of the L-shaped marker 20 in each shooting. ), The measurement coordinates u (one-dimensional image coordinates) on the image are obtained and stored in the storage unit 31 as the measurement coordinate data. The coordinates on the image of the black and white band 23 (detection point) of the L-shaped marker 20 are detected by performing image processing such as binarization processing and edge detection.

又、外部パラメータ算出部34では、計測座標算出部33で算出した各撮影での計測座標u、レンズ初期値、L字マーカ20の白黒帯23(検出点)の真値座標X、Yから、上述した式(11)~(17)を用いて、各撮影でのラインセンサカメラ10の外部パラメータを求め、記憶部31へ保管している。 Further, in the external parameter calculation unit 34, the measurement coordinates u in each shooting calculated by the measurement coordinate calculation unit 33, the initial lens value, and the true value coordinates X and Y of the black and white band 23 (detection point) of the L-shaped marker 20 are used. Using the above-mentioned equations (11) to (17), the external parameters of the line sensor camera 10 in each shooting are obtained and stored in the storage unit 31.

又、非線形最小二乗部35では、計測座標算出部33で算出した各撮影での計測座標u、レンズ初期値、L字マーカ20の白黒帯23(検出点)の真値座標X、Y、外部パラメータ算出部34で求めた各撮影でのラインセンサカメラ10の外部パラメータを用い、最終的なラインセンサカメラ10の内部パラメータ及び各撮影でのラインセンサカメラ10の外部パラメータについて、上述した式(8)をレーベンバーグマーカート法により算出し、記憶部31へ保管している。 Further, in the non-linear minimum squared unit 35, the measured coordinates u in each shooting calculated by the measured coordinate calculation unit 33, the initial value of the lens, the true value coordinates X, Y of the black and white band 23 (detection point) of the L-shaped marker 20, and the outside. Using the external parameters of the line sensor camera 10 for each shooting obtained by the parameter calculation unit 34, the above-mentioned equation (8) is used for the final internal parameters of the line sensor camera 10 and the external parameters of the line sensor camera 10 for each shooting. ) Is calculated by the Levenberg-Markart method and stored in the storage unit 31.

次に、図1と共に図2を参照して、本発明に係るラインセンサカメラのキャリブレーション方法を説明する。 Next, the calibration method of the line sensor camera according to the present invention will be described with reference to FIG. 2 together with FIG.

(ステップS1)
画像入力工程において、位置を固定したL字マーカ20に対し、様々な位置、姿勢(w軸周りの回転角θおよびu軸及びv軸に平行な向きの位置t1,t2の少なくともひとつが異なり、かつ、立体マーカ全体が撮影できる位置、姿勢)になるようにラインセンサカメラ10を設置し、それぞれの位置、姿勢での画像を撮像し、記憶部31へ保管する(画像入力部32)。撮像した画像の一例を図3A,図3Bに示す。L字マーカ20は、その位置を固定しているので、設置誤差は含まれないようになっている。
(Step S1)
In the image input process, at least one of various positions and postures (rotation angle θ around the w-axis and positions t1 and t2 parallel to the u-axis and v-axis is different from the L-shaped marker 20 whose position is fixed. In addition, the line sensor camera 10 is installed so that the entire three-dimensional marker can be photographed (position and posture), images at each position and posture are captured, and stored in the storage unit 31 (image input unit 32). An example of the captured image is shown in FIGS. 3A and 3B. Since the position of the L-shaped marker 20 is fixed, the installation error is not included.

(ステップS2)
計測座標算出工程において、撮影した画像毎に、L字マーカ20の白黒帯23(検出点)の検出処理を行い、L字マーカ20の白黒帯23(検出点)の画像上の計測座標u(1次元画像座標)を算出し、計測座標データとして記憶部31へ保管する(計測座標算出部33)。
(Step S2)
In the measurement coordinate calculation step, the black-and-white band 23 (detection point) of the L-shaped marker 20 is detected for each captured image, and the measurement coordinates u (detection point) of the black-and-white band 23 (detection point) of the L-shaped marker 20 are performed. (One-dimensional image coordinates) is calculated and stored as measurement coordinate data in the storage unit 31 (measurement coordinate calculation unit 33).

(ステップS3)
外部パラメータ算出工程において、ラインセンサカメラ10の初期の内部パラメータとなるレンズ初期値として、焦点距離fをメーカの公表値、主点座標cを総画素数の半分、歪み係数k1~k3を全て0とし、撮影した画像毎に、上述した式(11)~(17)を用い、線形解法により初期の外部パラメータを算出し、記憶部31へ保管する(外部パラメータ算出部34)。
(Step S3)
In the external parameter calculation process, the focal length f is the manufacturer's published value, the principal point coordinate c is half the total number of pixels, and the distortion coefficients k 1 to k 3 are set as the lens initial values that are the initial internal parameters of the line sensor camera 10. All are set to 0, and the initial external parameters are calculated by the linear solution method using the above-mentioned equations (11) to (17) for each captured image and stored in the storage unit 31 (external parameter calculation unit 34).

(ステップS4)
非線形最小二乗工程において、上述した初期の内部パラメータ及びステップS3で算出した初期の外部パラメータを初期値として、撮影した画像毎に、上述した式(8)を用い、レーベンバーグマーカート法で最終的な内部及び外部パラメータを算出し、記憶部31へ保管する(非線形最小二乗部35)。
(Step S4)
In the nonlinear least squares step, the above-mentioned initial internal parameters and the initial external parameters calculated in step S3 are used as initial values, and the above-mentioned equation (8) is used for each photographed image, and the final result is obtained by the Levenberg-Marquart method. Internal and external parameters are calculated and stored in the storage unit 31 (nonlinear least squares unit 35).

本実施例は、1台のラインセンサカメラ10と1つのL字マーカ20を用い、以上の手順により、ラインセンサカメラ10の焦点距離f、主点座標c、歪み係数k1~k3といった内部パラメータを求めると共に、撮影毎のL字マーカ20の回転θ、平行移動t1、t2といった外部パラメータを求めることが可能となる。 In this embodiment, one line sensor camera 10 and one L-shaped marker 20 are used, and the focal length f, the principal point coordinates c, and the distortion coefficients k 1 to k 3 of the line sensor camera 10 are internally determined by the above procedure. In addition to obtaining the parameters, it is possible to obtain external parameters such as the rotation θ of the L-shaped marker 20 and the parallel movements t 1 and t 2 for each shooting.

以上説明したように、本実施例では、誤差関数だけでなく内部パラメータの値も考慮した評価関数を最小化する非線形解法を用いることで、レンズの歪みを含む内部パラメータと外部パラメータとを安定的かつ同時に算出することが可能となる。また、内部パラメータの値を考慮する際に重み係数を導入することにより、どの内部パラメータについて重要視するかといった専門家の知識を導入することができ、より好適にキャリブレーションを行うことが可能となる。 As described above, in this embodiment, by using a nonlinear solution method that minimizes the evaluation function that considers not only the error function but also the value of the internal parameter, the internal parameter including the distortion of the lens and the external parameter are stable. Moreover, it is possible to calculate at the same time. In addition, by introducing a weighting coefficient when considering the value of the internal parameter, it is possible to introduce expert knowledge such as which internal parameter is important, and it is possible to perform calibration more appropriately. Become.

10 ラインセンサカメラ
20 L字マーカ
23 白黒帯
30 処理システム
31 記憶部
32 画像入力部
33 計測座標算出部
34 外部パラメータ算出部
35 非線形最小二乗部
10 Line sensor camera 20 L-shaped marker 23 Black and white band 30 Processing system 31 Storage unit 32 Image input unit 33 Measurement coordinate calculation unit 34 External parameter calculation unit 35 Non-linear least squares unit

Claims (8)

ラインセンサカメラの撮像範囲となる2次元平面上での2次元平面座標が既知の複数の検出点を有する立体マーカと、
前記立体マーカを前記ラインセンサカメラで撮像した複数の画像に基づいて、前記ラインセンサカメラのキャリブレーション用の内部パラメータ及び外部パラメータを算出する処理部とを有し、
前記処理部は、
位置及び姿勢が各々異なる前記立体マーカの前記画像が各々入力される画像入力部と、
前記立体マーカの前記検出点の前記画像上の1次元画像座標を前記画像毎に算出する計測座標算出部と、
前記2次元平面座標から前記1次元画像座標へ変換する前記ラインセンサカメラの数式モデルに基づいて求めた変換行列について、前記変換行列における初期の内部パラメータとして既知の初期値を入力し、線形解法により前記変換行列から初期の外部パラメータを前記画像毎に算出する外部パラメータ算出部と、
前記数式モデルに基づいて求めたレンズ歪みを考慮した誤差関数に、前記内部パラメータを考慮した値を加算してなる評価関数について、前記評価関数に前記初期の内部パラメータ及び前記初期の外部パラメータを入力し、前記検出点の全数及び前記画像の全数についての前記誤差関数の総和に前記内部パラメータを考慮した値を加算した評価関数が最小となる前記内部パラメータ及び前記外部パラメータを、反復計算を用いた非線形解法により算出する非線形最小二乗部とを有する
ことを特徴とするラインセンサカメラのキャリブレーション装置。
A three-dimensional marker having a plurality of detection points whose two-dimensional plane coordinates on the two-dimensional plane which is the imaging range of the line sensor camera are known, and
It has a processing unit for calculating internal parameters and external parameters for calibration of the line sensor camera based on a plurality of images of the three-dimensional marker captured by the line sensor camera.
The processing unit
An image input unit in which the images of the three-dimensional markers having different positions and postures are input, and
A measurement coordinate calculation unit that calculates the one-dimensional image coordinates of the detection point of the three-dimensional marker on the image for each image, and the measurement coordinate calculation unit.
For the transformation matrix obtained based on the mathematical model of the line sensor camera that converts the two-dimensional plane coordinates to the one-dimensional image coordinates, input known initial values as the initial internal parameters in the transformation matrix, and use a linear solution method. An external parameter calculation unit that calculates initial external parameters from the transformation matrix for each image,
For the evaluation function obtained by adding the value considering the internal parameters to the error function considering the lens distortion obtained based on the mathematical model, the initial internal parameters and the initial external parameters are input to the evaluation function. Then, the internal parameter and the external parameter, which minimizes the evaluation function obtained by adding the value considering the internal parameter to the sum of the error functions for the total number of the detection points and the total number of the images, were repeatedly calculated. A line sensor camera calibration device characterized by having a non-linear minimum squared portion calculated by a non-linear solution method.
請求項1に記載のラインセンサカメラのキャリブレーション装置において、
前記評価関数は、レンズの焦点距離の大きさ、レンズの主点座標の位置及びレンズの歪み係数の大きさにそれぞれ重み付けをして前記誤差関数に加算したものである
ことを特徴とするラインセンサカメラのキャリブレーション装置。
In the calibration device for the line sensor camera according to claim 1,
The evaluation function is a line sensor characterized in that the magnitude of the focal length of the lens, the position of the principal point coordinates of the lens, and the magnitude of the distortion coefficient of the lens are each weighted and added to the error function. Camera calibration device.
請求項1または請求項2に記載のラインセンサカメラのキャリブレーション装置において、
前記立体マーカは、L字型の形状であり、L字型に折れ曲がった内側の面に前記複数の検出点が配置されている
ことを特徴とするラインセンサカメラのキャリブレーション装置。
In the calibration device for the line sensor camera according to claim 1 or 2.
The three-dimensional marker has an L-shaped shape, and is a calibration device for a line sensor camera, characterized in that the plurality of detection points are arranged on an inner surface bent into an L-shape.
請求項1から請求項3のいずれか1項に記載のラインセンサカメラのキャリブレーション装置において、
前記数式モデルは、ピンホールカメラモデルを用いたものである
ことを特徴とするラインセンサカメラのキャリブレーション装置。
In the calibration device for the line sensor camera according to any one of claims 1 to 3.
The mathematical model is a calibration device for a line sensor camera, which is characterized by using a pinhole camera model.
ラインセンサカメラの撮像範囲となる2次元平面上での2次元平面座標が既知の複数の検出点を有する立体マーカを前記ラインセンサカメラで撮像した複数の画像に基づいて、前記ラインセンサカメラのキャリブレーション用の内部パラメータ及び外部パラメータを算出するラインセンサカメラのキャリブレーション方法であって、
位置及び姿勢が各々異なる前記立体マーカの前記画像が各々入力される画像入力工程と、
前記立体マーカの前記検出点の前記画像上の1次元画像座標を前記画像毎に算出する計測座標算出工程と、
前記2次元平面座標から前記1次元画像座標へ変換する前記ラインセンサカメラの数式モデルに基づいて求めた変換行列について、前記変換行列における初期の内部パラメータとして既知の初期値を入力し、線形解法により前記変換行列から初期の外部パラメータを前記画像毎に算出する外部パラメータ算出工程と、
前記数式モデルに基づいて求めたレンズ歪みを考慮した誤差関数に、前記内部パラメータを考慮した値を加算してなる評価関数について、前記評価関数に前記初期の内部パラメータ及び前記初期の外部パラメータを入力し、前記検出点の全数及び前記画像の全数についての前記誤差関数の総和に前記内部パラメータを考慮した値を加算した評価関数が最小となる前記内部パラメータ及び前記外部パラメータを、反復計算を用いた非線形解法により算出する非線形最小二乗工程とを有する
ことを特徴とするラインセンサカメラのキャリブレーション方法。
Calibration of the line sensor camera based on a plurality of images taken by the line sensor camera of a three-dimensional marker having a plurality of detection points whose two-dimensional plane coordinates are known on the two-dimensional plane which is the imaging range of the line sensor camera. A line sensor camera calibration method that calculates internal and external parameters for the camera.
An image input step in which the images of the three-dimensional markers having different positions and postures are input, and
A measurement coordinate calculation step of calculating the one-dimensional image coordinates of the detection point of the three-dimensional marker on the image for each image, and
For the transformation matrix obtained based on the mathematical model of the line sensor camera that converts the two-dimensional plane coordinates to the one-dimensional image coordinates, input known initial values as the initial internal parameters in the transformation matrix, and use a linear solution method. An external parameter calculation step of calculating the initial external parameters from the transformation matrix for each image, and
For the evaluation function obtained by adding the value considering the internal parameters to the error function considering the lens distortion obtained based on the mathematical model, the initial internal parameters and the initial external parameters are input to the evaluation function. Then, the internal parameter and the external parameter, which minimizes the evaluation function obtained by adding the value considering the internal parameter to the sum of the error functions for the total number of the detection points and the total number of the images, were repeatedly calculated. A line sensor camera calibration method comprising a non-linear minimum square step calculated by a non-linear solution method.
請求項5に記載のラインセンサカメラのキャリブレーション方法において、
前記評価関数として、レンズの焦点距離の大きさ、レンズの主点座標の位置及びレンズの歪み係数の大きさにそれぞれ重み付けをして前記誤差関数に加算したものを用いる
ことを特徴とするラインセンサカメラのキャリブレーション方法。
In the calibration method of the line sensor camera according to claim 5,
As the evaluation function, a line sensor characterized in that the size of the focal length of the lens, the position of the principal point coordinates of the lens, and the size of the distortion coefficient of the lens are each weighted and added to the error function. How to calibrate the camera.
請求項5または請求項6に記載のラインセンサカメラのキャリブレーション方法において、
前記立体マーカとして、L字型の形状であって、L字型に折れ曲がった内側の面に前記複数の検出点が配置されているものを用いる
ことを特徴とするラインセンサカメラのキャリブレーション方法。
In the calibration method of the line sensor camera according to claim 5 or 6.
A method for calibrating a line sensor camera, which comprises using as the three-dimensional marker an L-shaped marker in which the plurality of detection points are arranged on an inner surface bent into an L-shape.
請求項5から請求項7のいずれか1項に記載のラインセンサカメラのキャリブレーション方法において、
前記数式モデルとして、ピンホールカメラモデルを用いる
ことを特徴とするラインセンサカメラのキャリブレーション方法。
The method for calibrating a line sensor camera according to any one of claims 5 to 7.
A method for calibrating a line sensor camera, which comprises using a pinhole camera model as the mathematical model.
JP2018166067A 2018-09-05 2018-09-05 Line sensor camera calibration device and method Active JP7044016B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018166067A JP7044016B2 (en) 2018-09-05 2018-09-05 Line sensor camera calibration device and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018166067A JP7044016B2 (en) 2018-09-05 2018-09-05 Line sensor camera calibration device and method

Publications (2)

Publication Number Publication Date
JP2020038550A JP2020038550A (en) 2020-03-12
JP7044016B2 true JP7044016B2 (en) 2022-03-30

Family

ID=69738564

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018166067A Active JP7044016B2 (en) 2018-09-05 2018-09-05 Line sensor camera calibration device and method

Country Status (1)

Country Link
JP (1) JP7044016B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111667539B (en) * 2020-06-08 2023-08-29 武汉唯理科技有限公司 Camera calibration and plane measurement method
CN112082512B (en) * 2020-09-08 2023-04-14 深圳广成创新技术有限公司 Calibration optimization method and device for phase measurement deflection technique and computer equipment
CN112634375B (en) * 2020-12-21 2022-08-05 杭州东信北邮信息技术有限公司 Plane calibration and three-dimensional reconstruction method in AI intelligent detection
JP7517170B2 (en) 2021-01-22 2024-07-17 株式会社明電舎 Calibration device and calibration method for line sensor camera

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009288165A (en) 2008-05-30 2009-12-10 Jfe Steel Corp Calibration method of line sensor camera, calibration device and program
JP2016218815A (en) 2015-05-22 2016-12-22 株式会社明電舎 Calibration device and method for line sensor camera
JP2017161245A (en) 2016-03-07 2017-09-14 株式会社明電舎 Stereo calibration device, and stereo calibration method, for line sensor cameras

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6645140B2 (en) * 2015-11-27 2020-02-12 株式会社明電舎 Image calibration apparatus and calibration method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009288165A (en) 2008-05-30 2009-12-10 Jfe Steel Corp Calibration method of line sensor camera, calibration device and program
JP2016218815A (en) 2015-05-22 2016-12-22 株式会社明電舎 Calibration device and method for line sensor camera
JP2017161245A (en) 2016-03-07 2017-09-14 株式会社明電舎 Stereo calibration device, and stereo calibration method, for line sensor cameras

Also Published As

Publication number Publication date
JP2020038550A (en) 2020-03-12

Similar Documents

Publication Publication Date Title
CN110689579B (en) Rapid monocular vision pose measurement method and measurement system based on cooperative target
JP7044016B2 (en) Line sensor camera calibration device and method
JP6641729B2 (en) Line sensor camera calibration apparatus and method
EP2795579B1 (en) Image registration method and system robust to noise
CN108510551B (en) Method and system for calibrating camera parameters under long-distance large-field-of-view condition
EP2926543B1 (en) A method of calibrating a camera and a system therefor
JP6079333B2 (en) Calibration apparatus, method and program
JP6217227B2 (en) Calibration apparatus, method and program
WO2016042779A1 (en) Triangulation device, triangulation method, and recording medium recording program therefor
US10771776B2 (en) Apparatus and method for generating a camera model for an imaging system
CN113920205B (en) Calibration method of non-coaxial camera
CN112465877B (en) Kalman filtering visual tracking stabilization method based on motion state estimation
CN114820739B (en) Multispectral camera-oriented image rapid registration method and device
CN116433737A (en) Method and device for registering laser radar point cloud and image and intelligent terminal
CN113706635A (en) Long-focus camera calibration method based on point feature and line feature fusion
JP7121936B2 (en) Camera calibration information acquisition device, image processing device, camera calibration information acquisition method and program
JP2023538706A (en) Fusion-based digital image correlation framework for performing distortion measurements
CN108416811B (en) Camera self-calibration method and device
JP7517170B2 (en) Calibration device and calibration method for line sensor camera
Claus et al. A Plumbline Constraint for the Rational Function Lens Distortion Model.
CN116051634A (en) Visual positioning method, terminal and storage medium
CN112927299B (en) Calibration method and device and electronic equipment
WO2023166617A1 (en) Camera parameter estimation device, camera parameter estimation method, and computer-readable recording medium
CN114758015B (en) Collimator outgoing beam modeling and camera internal parameter calibration method based on general imaging model
CN116801100B (en) Calibration focal length verification method and device for binocular camera module in automatic driving system

Legal Events

Date Code Title Description
RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20190529

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20190529

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210203

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20220215

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220228

R150 Certificate of patent or registration of utility model

Ref document number: 7044016

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150