JPH09166421A - 表面三次元形状の測定方法 - Google Patents

表面三次元形状の測定方法

Info

Publication number
JPH09166421A
JPH09166421A JP7327595A JP32759595A JPH09166421A JP H09166421 A JPH09166421 A JP H09166421A JP 7327595 A JP7327595 A JP 7327595A JP 32759595 A JP32759595 A JP 32759595A JP H09166421 A JPH09166421 A JP H09166421A
Authority
JP
Japan
Prior art keywords
shape
measured
light
estimated
estimated shape
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP7327595A
Other languages
English (en)
Inventor
Yoshiyuki Sonda
嘉之 尊田
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.)
AGC Inc
Original Assignee
Asahi Glass 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 Asahi Glass Co Ltd filed Critical Asahi Glass Co Ltd
Priority to JP7327595A priority Critical patent/JPH09166421A/ja
Publication of JPH09166421A publication Critical patent/JPH09166421A/ja
Pending legal-status Critical Current

Links

Landscapes

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

Abstract

(57)【要約】 【課題】透明板状体の光学歪を、人間が視覚する歪に見
合ったものとして評価する方法を得る。 【解決手段】視対象物上の実点から発せられた光線群の
瞳に入射する方向から、透明板状体を介して視対象物を
観察した場合の虚像点群を決め、この虚像点群によって
視対象物の虚像の三次元形状を形成し(14)、得られ
た虚像の三次元形状の奥行き幅から、透明板状体の透視
歪の大きさを評価する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、光を反射する面を
有する被測定物の外観検査方法に関し、特に、建築用、
車両用等のガラス板やブラウン管等の表面三次元形状の
測定方法に関する。
【0002】
【従来の技術】自動車用、建築用等のガラス板やブラウ
ン管等の成形後の表面三次元形状の検査には、その被測
定物の所定形状ごとの検査ゲージに被測定物を載置し
て、接触式で検査する方法が一般的であった。しかし、
この方法では検査ゲージが型式ごとに必要であり、載置
を含めた検査に要する時間が多く、さらにはゲージへの
接触によって被測定物の変形の恐れがあった。
【0003】そこで、被測定物の三次元形状の検査や測
定には、上記の接触式よりも非接触式のものが好ましい
との見方ができた。この非接触式の三次元形状の測定方
法としては、例えば、複数の点状の散乱光源から発せら
れた光を被測定物の表面へ照射し、その反射光の受光装
置への入射方向を手がかりに表面の傾きを算出すること
で、被測定物の表面三次元形状を得る方法(特開平3−
44504)などが提案されている。
【0004】
【発明が解決しようとする課題】上記の測定方法は、具
体的には光源から発せられた光の被測定物からの反射光
によって複数の局所平面を形成しこれらを局所曲面とみ
なし、これらの局所曲面を接続して局所曲面の集合体と
して三次元形状を求めるものである。
【0005】しかし、この方法は次に示す欠点があっ
た。 (1)局所平面を局所曲面とみなすため、曲率が大きな
被測定物については、測定誤差が大きい。 (2)被測定物を局所的に表現するにあたって、各々の
局所曲面間の境界における接続性の考慮が困難である。 (3)初点から順次接続していくため、初点から離れる
にしたがって誤差が蓄積される。
【0006】本発明の目的は、従来技術の前述の欠点を
解消する新規な表面三次元形状の測定方法を提供するこ
とにある。
【0007】
【課題を解決するための手段】本発明は、光を反射する
面を有する被測定物の表面に実質的に点状の散乱光源か
ら発せられた光を照射し、その反射光線を受光装置で受
光して反射光の受光装置への入射方向に延びる直線を特
定する工程Aと、被測定物の推定形状を関数の線形結合
によって設定する工程Bと、前記直線と推定形状の面と
の交点を定める工程Cと、複数の光源位置から発せられ
た光に関し工程AおよびCを経て複数の交点を定める工
程Dと、これら複数の交点における各々の光源から受光
装置までの光路長の変分を用いた評価指標が規定値以下
になるまで推定形状を補正し、最終的な推定形状を被測
定物の測定形状とする工程Eと、を有する表面三次元形
状の測定方法を提供する。
【0008】
【発明の実施の形態】以下、図面に基づいて本発明をさ
らに詳細に説明する。図1は、本発明の表面三次元形状
を測定する工程の一例を説明する流れ図である。また図
2は、本発明における被測定物の三次元形状の測定状態
を示す概略斜視図である。
【0009】光源2から発せられた光は、被測定物1で
反射されて受光装置3で受光される。この受光データに
基づいて、演算装置4にて演算され、被測定物1の三次
元形状が測定される。なお、受光装置3としては、カメ
ラ、ビデオカメラ、スティルカメラ、フォトセンサをマ
トリックス状に配置したもの等、反射光の入射方向を特
定しうるものである。
【0010】具体的には、点状の散乱光源2から発せら
れた光が被測定物1面上で反射して、受光装置3でその
入射方向を測定する(11)。別途、被測定物の形状を
関数の線形結合によって推定、表現しておき、これを初
期推定形状として演算装置4に設定する(12)。11
で得られた入射方向に延びる直線と上記の推定形状との
交点を求め(13)、さらに複数の光源位置から発せら
れる光に関するこの交点を求める。
【0011】これらの各交点における光源から受光装置
までの光路長の変分(以下簡単に交点の変分、交点にお
ける変分、光路長の変分という)の二乗和を評価指標と
し、この評価指標が規定値以下になるかの判定を行う
(14)。規定値以下であれば、その推定形状が被測定
物の求める表面三次元形状の測定結果となる。
【0012】評価指標が規定値を超えている場合、被測
定物の表面三次元形状が、推定した形状から掛け離れて
いることを示す。したがって、被測定物の表面三次元形
状を得るためには、推定形状を修正する必要がある。こ
の修正を行って(15)、先に測定した光の入射方向に
延びる直線と修正後の推定形状との交点を求め、評価指
標が規定値になるまで上記の工程を繰り返すことによっ
て、被測定物の表面三次元形状が得られる。
【0013】図3は光学における変分原理を説明する概
略断面図である。光源2から発せられて被測定物1で反
射した光は、反射像5aとして受光装置の撮像面上の点
6に撮像され、受光装置へは7の方向から入射したと認
識される。8は受光装置の視点位置である。このとき、
光線経路2→5a→8の光路長は近傍のいかなる光線経
路2→5b→8の光路長よりも短くなる。すなわち、5
aにおいて光路長の変分が0となる。
【0014】図4は本発明の原理を説明する概略断面図
である。受光装置の視点位置8を通り、2aの位置にあ
る散乱光源から発せられた光の被測定物による反射光の
入射方向へ延びる直線10aと、被測定物の推定形状9
の面との交点を11aとする。この交点11aは、光源
2aの座標、推定形状9(推定形状上の座標)、視点位
置8の座標がわかっているため、容易に求めうる。
【0015】推定形状が被測定物の表面三次元形状と一
致していれば、11aにおける光路長の変分は0とな
る。(しかし、光路長の変分が0になることは、推定形
状が被測定物の表面三次元形状と一致する必要条件であ
って、十分条件ではない。)ここでは、被測定物の推定
形状は、被測定物の実際の表面三次元形状にできるだけ
近い形状を予想して設定するのが通常であるため、基本
的には、推定形状が被測定物の表面三次元形状と異なっ
ていれば光路長の変分は0とはならない。
【0016】ほかの位置2bに散乱光源がある場合につ
いても同様であり、2bの位置にある光源から発せられ
た光の被測定物による反射光の入射方向へ延びる直線1
0bと、被測定物の推定形状9の面との交点11bの変
分を求める。このように多数の光源の位置2a、2b・
・・・について、交点11a、11b・・・・の変分を
求め、これらの変分が0であれば、被測定物の形状は推
定形状に一致しているといえる。
【0017】仮に0でなくても、きわめて小さい値であ
れば、推定形状が実際の被測定物の形状に比べて大きく
違わないことになる。すなわち、実際の被測定物に比べ
ての推定形状の違いが許容できる範囲であれば、変分が
0にならなくても、そのときの推定形状を被測定物の測
定形状とみなしうる。この場合、その判断基準(推定形
状の被測定物形状への一致の程度)として、上記の各交
点における変分の二乗和を評価指標とし、あらかじめ評
価指標の上限となる規定値を定めておき、評価指標をこ
の規定値以下にする推定形状を、被測定物の測定形状と
する。
【0018】この規定値としては、推定形状が実際の形
状に一致する点に鑑みれば、0が好ましい。この値0に
ついても、受光装置の解像度等を考慮して、小数点以下
数桁で切り捨てることによって、完全に0の値でなくて
も0とみなすこともできる。
【0019】上記の規定値以下にならない場合、推定形
状は実際の被測定物の表面三次元形状と異なっているこ
とになる。そのため、推定した形状を実際の形状に近づ
けるために推定形状に補正を加える必要がある。
【0020】その補正の方法の一例として、最適化手法
がある。これは、非線形関数である複数の交点の変分の
二乗和の零点を求めるものであり、ニュートン法をはじ
め、準ニュートン法、最急降下法、共役勾配法、方向集
合法等が例示できる。このような最適化手法によって、
複数の交点の変分の二乗和、すなわち評価指標が規定値
以下、好ましくは0になるまで推定形状を補正し、評価
指標を規定値以下、好ましくは0にする推定形状を、被
測定物の測定形状とする。
【0021】本発明において用いる評価指標としては、
上記の交点の変分の二乗和のほか、変分の絶対値の和
等、適宜の値を採用できるが、通常は変分の二乗和を用
いる。
【0022】本発明において、推定形状は関数の線形結
合によって表現される。この関数としては、初等n次関
数、B−スプライン関数、三角関数等、表面を表現でき
る関数であれば特に制限はない。そして、これらの関数
に(線形結合)係数をかけて線形結合させることによっ
て、面形状を表現できる。さらに、この関数としてB−
スプライン関数、三角関数などの基底関数の使用が好ま
しく、最適化計算の収束が速い点から、B−スプライン
関数の使用が特に好ましい。
【0023】本発明における散乱光源は実質的に点状の
ものであり、例えばLED、粗面状のスクリーン上にレ
ーザ光を照射したもの、蛍光灯等があげられる。本発明
において、被測定物に照射される光が発せられる位置は
複数設定される。そのためには、複数の光源、例えばL
EDを複数個配列し、順次点滅させることによって、被
測定物に照射される光が発せられる位置を複数にするこ
とができる。ほかに、レーザ光等を走査させる等、ひと
つの光源の位置を移動させることによっても、被測定物
に照射される光が発せられる位置を複数にすることがで
きる。
【0024】上記の三次元形状の測定は、被測定物が建
築用や自動車用等の曲げガラス板である場合に、ガラス
板の曲げ成形後の形状が所定(例えば設計上のCADデ
ータ等)の形状であるかを検査する場合に、好ましく用
いられる。建築用や自動車用の曲げガラス板は、窓の開
口部に取り付けられるものであり、その曲げ精度が悪い
場合には取り付け寸法精度が悪くなり、水密性や開口部
への納まり具合等の劣化を招く。そのため、成形後のガ
ラス板の検査が重要になり、上記方法の採用によってそ
の検査を非接触で行いうる。
【0025】さらに、この検査によって所定の形状に比
べて曲率が大きい、端縁や中央部等局所的に曲率が小さ
い、等の形状不良が見いだされた場合に、曲げ成形装置
にその結果をフィードバックすることにより、所定の形
状のガラス板が得られる。例えば、形状不良のあり方に
よって、加熱炉の温度を局所的にまたは全体的に上昇さ
せたり下降させる等、成形条件を調整することによっ
て、曲げ成形の精度を上げうる。
【0026】これをさらに発展させれば、ガラス板の曲
げ加工後の曲げ成形装置における下流側に、本発明の三
次元形状の測定方法を用いた装置を配することによっ
て、成形されるガラス板の曲げ成形精度を成形後ただち
に把握できる。この場合、本発明の測定方法は、ガラス
板に非接触で行いうる。このため、成形直後であって
も、容易に三次元形状の測定を行いうる。こうして、成
形装置に対して成形条件設定の補正をリアルタイムで行
うことができ、効率のよいガラス板の曲げ成形が実現で
きる。
【0027】上記の曲げガラス板のほか、平板状のガラ
ス板の表面三次元形状、すなわち表面の平滑性、凹凸の
有無の検査や、ブラウン管等の表面形状の検査もでき
る。ガラス板のほかにも表面が光を反射しうるものであ
れば、樹脂ガラス板と呼ばれる有機透明樹脂板、金属
体、樹脂成形品等の表面三次元形状の測定もできる。
【0028】
【実施例】以下、本発明における三次元形状の測定の手
法を用いた例を示す。本例では、図5のように座標系を
設定した。原点は受光装置の視点位置8であり、z軸は
視点位置8から被測定物に下ろした垂線である。
【0029】本例において、図1に示した実際の被測定
物1に散乱光源2から照射されて、その反射光を受光装
置3で受光するかわりに、被測定物の形状を数式で与え
て数値計算によって、本発明における手法が被測定物の
三次元形状を正確に測定できるかを確認した。
【0030】具体的には、まず、以下の式で被測定物の
面形状を表した。なお、G=36.0(m)、h=0.
0284756(m)、s=3.0(m)とした。この
形状は図6に示すような形状であった。
【0031】
【数1】
【0032】次に、直交格子上に配置された複数の散乱
光源について、上記の式によって表される形状の面での
反射光の入射方向に延びる直線を得た。一方、被測定物
の推定形状として、図7に示すような被測定物の初期推
定形状を設定した。このときの推定形状は、B−スプラ
イン関数を基底関数に選びそれらの線形結合により設定
しており、B−スプライン関数の線形結合における線形
結合係数をすべて同じ値とした。
【0033】その後、複数の散乱光源各々について、反
射光の受光装置への入射方向に延びる直線と推定形状面
との交点を定め、これら複数の交点の変分の二乗和を評
価指標とした。この評価指標が0になるまで、推定形状
の面を表す式に共役勾配法によって、B−スプライン関
数の線形結合係数の値を変化させる補正を加え、評価指
標が0になるときの推定形状を図8に示す形状として得
た。
【0034】こうして得られた推定形状を測定形状と
し、図8と図6とを比較すると、測定形状は被測定物の
形状をよく再現していることが確認できた。
【0035】
【発明の効果】以上のように本発明によれば、ガラス板
など粗面でない被測定物の表面形状を非接触で測定する
にあたり、推定形状を関数の線形結合により表現してい
るため測定形状の接続性に優れ、曲率の大きな被測定物
に対しても精度を落とすことなく測定できる。
【0036】特に、初期拘束点の近傍から表面形状を順
次推定するのではなく、最適化手法により、関数の線形
結合における線形結合係数を適宜変化させ、表面形状全
体を同時に修正すれば、測定精度を表面内で均一に、か
つ高く保つことができる。
【図面の簡単な説明】
【図1】本発明における光学歪の評価に用いる視対象物
の虚像を得る方法の一例を説明する流れ図
【図2】本発明における透明板状体と視対象物の配置構
成を示す概略斜視図
【図3】光学における変分原理を説明する概略断面図
【図4】本発明の原理を説明する概略断面図
【図5】本発明における実施例に用いた座標系を示す概
略斜視図
【図6】本発明における実施例に用いた被測定物の形状
を示す概略斜視図
【図7】本発明における実施例に用いた推定形状を示す
概略斜視図
【図8】本発明における実施例に用いた測定形状を示す
概略斜視図
【符号の説明】 1:被測定物 2:点状の散乱光源 3:受光装置 4:演算装置

Claims (4)

    【特許請求の範囲】
  1. 【請求項1】光を反射する面を有する被測定物の表面に
    実質的に点状の散乱光源から発せられた光を照射し、そ
    の反射光線を受光装置で受光して反射光の受光装置への
    入射方向に延びる直線を特定する工程Aと、被測定物の
    推定形状を関数の線形結合によって設定する工程Bと、
    前記直線と推定形状の面との交点を定める工程Cと、複
    数の光源位置から発せられた光に関し工程AおよびCを
    経て複数の交点を定める工程Dと、これら複数の交点に
    おける各々の光源から受光装置までの光路長の変分を用
    いた評価指標が規定値以下になるまで推定形状を補正
    し、最終的な推定形状を被測定物の測定形状とする工程
    Eと、を有する表面三次元形状の測定方法。
  2. 【請求項2】工程Eにおいて、前記評価指標を変分の二
    乗和とし、この変分の二乗和が0になるまで推定形状を
    補正することを特徴とする請求項1の表面三次元形状の
    測定方法。
  3. 【請求項3】工程Eにおいて、最適化手法によって推定
    形状を補正することを特徴とする請求項1または2の表
    面三次元形状の測定方法。
  4. 【請求項4】工程Eにおいて、前記推定形状を表す関数
    の線形結合の線形結合係数を、最適化手法によって補正
    することを特徴とする請求項3の表面三次元形状の測定
    方法。
JP7327595A 1995-12-15 1995-12-15 表面三次元形状の測定方法 Pending JPH09166421A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP7327595A JPH09166421A (ja) 1995-12-15 1995-12-15 表面三次元形状の測定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP7327595A JPH09166421A (ja) 1995-12-15 1995-12-15 表面三次元形状の測定方法

Publications (1)

Publication Number Publication Date
JPH09166421A true JPH09166421A (ja) 1997-06-24

Family

ID=18200820

Family Applications (1)

Application Number Title Priority Date Filing Date
JP7327595A Pending JPH09166421A (ja) 1995-12-15 1995-12-15 表面三次元形状の測定方法

Country Status (1)

Country Link
JP (1) JPH09166421A (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001052196A (ja) * 1999-08-06 2001-02-23 Mitsubishi Electric Inf Technol Center America Inc 文字枠により囲まれた距離フィールドのディテール指向型階層表現を生成する方法及びデータを記憶するメモリ
AU766662B2 (en) * 2000-05-16 2003-10-23 Canon Kabushiki Kaisha Rapid detection of whether paths of arbitrary complexity intersect
JP2006071372A (ja) * 2004-08-31 2006-03-16 Sunway Burekkusu:Kk 路面断面プロファイルの測定方法
JP2007315968A (ja) * 2006-05-26 2007-12-06 Univ Of Tokyo モーションキャプチャを用いた運動学パラメータの推定法及び装置
JP2021089265A (ja) * 2019-11-26 2021-06-10 みどり精密工業株式会社 炉壁の形状・損耗測定装置及び炉壁の形状・損耗測定方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001052196A (ja) * 1999-08-06 2001-02-23 Mitsubishi Electric Inf Technol Center America Inc 文字枠により囲まれた距離フィールドのディテール指向型階層表現を生成する方法及びデータを記憶するメモリ
AU766662B2 (en) * 2000-05-16 2003-10-23 Canon Kabushiki Kaisha Rapid detection of whether paths of arbitrary complexity intersect
JP2006071372A (ja) * 2004-08-31 2006-03-16 Sunway Burekkusu:Kk 路面断面プロファイルの測定方法
JP2007315968A (ja) * 2006-05-26 2007-12-06 Univ Of Tokyo モーションキャプチャを用いた運動学パラメータの推定法及び装置
JP2021089265A (ja) * 2019-11-26 2021-06-10 みどり精密工業株式会社 炉壁の形状・損耗測定装置及び炉壁の形状・損耗測定方法

Similar Documents

Publication Publication Date Title
JP4924426B2 (ja) 形状検査方法および装置
US8111907B2 (en) Method for repeatable optical determination of object geometry dimensions and deviations
EP2102587B1 (en) Method of automated quantitative analysis of distortion shaped vehicle glass by reflected optical imaging
EP2252856B1 (en) Method of determination of glass surface shape and optical distortion by reflected optical imaging
JP6519265B2 (ja) 画像処理方法
Liu et al. An improved online dimensional measurement method of large hot cylindrical forging
JP4998711B2 (ja) 面歪の測定装置及び方法
CN1158684A (zh) 自动图像监视器校准系统中坐标转换的方法与装置
WO2016145582A1 (zh) 相位偏移校准方法、3d形状检测的方法、系统及投影系统
KR19980063676A (ko) 대향하는 부분들 사이의 간극 및 부정합을 측정하기 위한시스템
US9797833B2 (en) Method for determining the refractive power of a transparent object, and corresponding device
JP4613626B2 (ja) 鏡面形状測定方法および装置並びに検査方法および装置
JPH09166421A (ja) 表面三次元形状の測定方法
US20130162816A1 (en) Device for measuring the shape of a mirror or of a specular surface
JP6951469B2 (ja) 光学測定装置をキャリブレーションするための方法
JP2572286B2 (ja) 3次元形状寸法計測装置
JP2009192483A (ja) 三次元形状計測方法および三次元形状計測装置
JPH0623992B2 (ja) 板ガラスの検査方法
JPH0660807B2 (ja) 光の中心位置の高精度計測方法
JP4806767B2 (ja) 歪み測定方法及び装置
JP2961140B2 (ja) 画像処理方法
JP2022030132A (ja) 目地寸法の計測方法
JP2870908B2 (ja) 透視歪の測定方法及びその装置
JP2006145231A (ja) 表面形状測定方法及び表面形状測定装置
EP1400802A1 (en) Method and arrangement for detecting and evaluating surface irregularities

Legal Events

Date Code Title Description
FPAY Renewal fee payment (prs date is renewal date of database)

Year of fee payment: 7

Free format text: PAYMENT UNTIL: 20080727

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080727

Year of fee payment: 7

FPAY Renewal fee payment (prs date is renewal date of database)

Year of fee payment: 8

Free format text: PAYMENT UNTIL: 20090727

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090727

Year of fee payment: 8

FPAY Renewal fee payment (prs date is renewal date of database)

Year of fee payment: 9

Free format text: PAYMENT UNTIL: 20100727

FPAY Renewal fee payment (prs date is renewal date of database)

Year of fee payment: 9

Free format text: PAYMENT UNTIL: 20100727

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110727

Year of fee payment: 10

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120727

Year of fee payment: 11

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130727

Year of fee payment: 12