JP6503568B2 - Dual image inspection system - Google Patents
Dual image inspection system Download PDFInfo
- Publication number
- JP6503568B2 JP6503568B2 JP2014250399A JP2014250399A JP6503568B2 JP 6503568 B2 JP6503568 B2 JP 6503568B2 JP 2014250399 A JP2014250399 A JP 2014250399A JP 2014250399 A JP2014250399 A JP 2014250399A JP 6503568 B2 JP6503568 B2 JP 6503568B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- double
- intensity distribution
- light
- light source
- 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
Links
- 238000007689 inspection Methods 0.000 title claims description 39
- 230000009977 dual effect Effects 0.000 title claims description 28
- 239000011521 glass Substances 0.000 claims description 39
- 238000012545 processing Methods 0.000 claims description 28
- 230000004044 response Effects 0.000 claims description 24
- 230000000007 visual effect Effects 0.000 claims description 18
- 238000005311 autocorrelation function Methods 0.000 claims description 14
- 238000002474 experimental method Methods 0.000 claims description 10
- 230000003287 optical effect Effects 0.000 claims description 10
- 238000000034 method Methods 0.000 claims description 9
- 238000003384 imaging method Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000005340 laminated glass Substances 0.000 description 9
- 230000007246 mechanism Effects 0.000 description 8
- 238000005259 measurement Methods 0.000 description 7
- 239000005357 flat glass Substances 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 210000001747 pupil Anatomy 0.000 description 2
- 230000003252 repetitive effect Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- -1 for example Substances 0.000 description 1
- 239000005341 toughened glass Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Landscapes
- Testing Of Optical Devices Or Fibers (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Description
本発明は、例えば合わせガラスなどの被検体ガラスで発生する二重像を検査する二重像検査システムに関する。 The present invention relates to a dual image inspection system for inspecting dual images generated on an object glass, such as, for example, laminated glass.
近年、乗用車のフロントガラスには、従来の強化ガラスに替わり、合わせガラスが使用されている。このフロントガラスをスクリーンとする虚像視ヘッドアップディスプレイ(Head-Up Display:HUD)が普及し始めている。虚像視HUDでは、プロジェクタからの投射画像が、合わせガラスの2つの界面で反射するので、一般的には、二重像が発生する。そこで、貼り合わされた2枚のガラスの間にある中間層に対して適当な楔角を与えることにより二重像を補正することが、例えば特許文献1で提案されている。
In recent years, laminated glass has been used as a windshield for passenger cars in place of conventional tempered glass. A virtual image head-up display (HUD) using this windshield as a screen is beginning to spread. In virtual vision HUD, a projected image from a projector is reflected at two interfaces of laminated glass, so a double image generally occurs. Therefore, for example,
特許文献1の記載によれば、曲面ガラスによって発生する二重像の視角と、楔角を有する平板ガラスによって発生する二重像の視角とが、互いに等しいと設定することにより、二重像を解消することの可能な条件を導出することができる、ようである。しかしながら、二重像の視角が0°になるように楔角を設定したり、合わせガラスを平板ガラスに置き換えたりした場合であっても、二重像が発生することがあることを本出願人は発見した。これは、特許文献1に記載の理論では説明ができない。この問題を検証するためには、発生した二重像を精密に検査することが必要であるが、従来の検査システムでは、二重像を精密に検査することが困難であった。
According to the description of
本発明はかかる問題点に鑑みてなされたもので、その目的は、二重像を精密に検査することの可能な二重像検査システムを提供することにある。 The present invention has been made in view of such problems, and an object thereof is to provide a dual image inspection system capable of precisely examining dual images.
本発明の一実施の形態としての二重像検査システムは、光源ユニット、カメラユニットおよび処理部を備える。光源ユニットは、光源と、光源に照明されるスリットと、スリットを介して光源から発せられた線状光を投影光として被検体ガラスに投影する投影光学系とを含んで構成される。カメラユニットは、光源ユニットから発せられ、被検体ガラスの界面で繰り返し反射された光を集光する対物レンズと、対物レンズにより集光された光を受光するイメージセンサとを含んで構成される。処理部は、カメラユニットにより得られた画像データを処理するようになっている。具体的には、処理部は、画像データに対して所定のデータ変換処理を行うことにより、イメージセンサの受光面に生ずる干渉縞パターンに対応した二重像の線像強度分布を求めるようになっている。処理部は、さらに、得られた線像強度分布に基づき、被検体ガラスにおける投影光の反射の際に生じた二重像の視角を求めるようになっている。 A dual image inspection system according to an embodiment of the present invention includes a light source unit, a camera unit, and a processing unit. The light source unit includes a light source, a slit illuminated by the light source, and a projection optical system that projects linear light emitted from the light source through the slit as projection light onto the subject glass. The camera unit is configured to include an objective lens for collecting light emitted from the light source unit and repeatedly reflected at the interface of the object glass , and an image sensor for receiving the light collected by the objective lens. The processing unit is configured to process image data obtained by the camera unit. Specifically, the processing unit performs predetermined data conversion processing on the image data to obtain a line image intensity distribution of a double image corresponding to the interference fringe pattern generated on the light receiving surface of the image sensor. ing. The processing unit is further configured to obtain the viewing angle of the double image generated upon reflection of the projection light on the subject glass based on the obtained linear image intensity distribution.
本発明の一実施の形態としての二重像検査システムによれば、二重像を精密に検査することができる。 According to the dual image inspection system as one embodiment of the present invention, double images can be precisely inspected.
以下、本発明の実施の形態について、図面を参照して詳細に説明する。以下の説明は本発明の一具体例であって、本発明は以下の態様に限定されるものではない。また、本発明は、各図に示す各構成要素の配置や寸法、寸法比などについても、それらに限定されるものではない。 Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings. The following description is one specific example of the present invention, and the present invention is not limited to the following embodiments. Further, the present invention is not limited to the arrangement, dimensions, dimensional ratio, and the like of each component shown in each drawing.
[構成]
図1は、本発明の一実施の形態に係る二重像検査システム1の全体構成例を模式的に表したものである。二重像検査システム1は、被検体ガラス100の二重像の視角θを検査するものであり、光源ユニット10、カメラユニット20および画像処理部30を備える。光源ユニット10は、本発明の「光源ユニット」の一具体例に対応する。カメラユニット20は、本発明の「カメラユニット」の一具体例に対応する。画像処理部30は、本発明の「処理部」の一具体例に対応する。被検体ガラス100は、例えば、合わせガラスであり、光源ユニット10から発せられる投影光Lを波面分割により分割可能な向きで配置されている。被検体ガラス100は、例えば、投影光Lが60°で入射するように配置された合わせガラスである。被検体ガラス100は、本発明の「被検体ガラス」の一具体例に対応する。
[Constitution]
FIG. 1 schematically shows an example of the entire configuration of a double-
光源ユニット10は、線状の投影光Lを出射するものである。光源ユニット10は、例えば、LED光源11、スリット12、レンズ13および投影レンズ14を、LED光源11から発せられた光の光路上に、この順に有している。LED光源11は、本発明の「光源」の一具体例に対応する。スリット12は、本発明の「スリット」の一具体例に対応する。投影レンズ14は、本発明の「投影光学系」の一具体例に対応する。
The
LED光源11は、インコヒーレント光を発する発光ダイオードであり、可視領域の波長の光を出射するようになっている。LED光源11は、例えば、白色発光ダイオードである。LED光源11は、スリット12を均一に照明するようになっている。スリット12は、LED光源11から発せられた点状光を、線状光に成形するようになっている。また、スリット12は、LED光源11から発せられた光を、位相のそろった光に変換するようになっている。つまり、スリット12は、波面のそろった線状光を出射するようになっている。スリット12の開口幅は、例えば、0.01mmとなっている。スリット12は、スリット12の長手方向と平行な線分がLED光源11から発せられた光の光軸と直交するとともに被検体ガラス100の表面に対して平行となるように、配置されている。つまり、スリット12は、図1の紙面に垂直な方向に延在している。レンズ13は、LED光源11から発せられた光を集光するものである。レンズ13は、例えば、凸レンズ13aおよび凹レンズ13bを、LED光源11から発せられた光の光路上に、この順に有する。凸レンズ13aは、例えば、コンデンサレンズである。投影レンズ14は、スリット12を介してLED光源11から発せられた線状光を投影光Lとして被検体ガラス100に投影するようになっている。具体的には、投影レンズ14は、スリット12で生成され、レンズ13で集光された線状光を、投影光Lとして被検体ガラス100に投影するようになっている。
The LED light source 11 is a light emitting diode that emits incoherent light, and emits light of a wavelength in the visible region. The LED light source 11 is, for example, a white light emitting diode. The LED light source 11 illuminates the
光源ユニット10は、LED光源11から発せられる光を凸レンズ13aの焦平面に沿って移動させる機構(第1の機構)をさらに有していてもよい。光源ユニット10は、第1の機構を用いて、LED光源11から発せられる光を凸レンズ13aの焦平面に沿って走査することにより、被検体ガラス100の広範囲で投影光Lを走査することができる。第1の機構は、例えば、LED光源11と第1の機構とを互いに連結する光ファイバと、光ファイバの光出射部を移動させるステージとを有する。レンズ13は、凹レンズ13bを2つ有していてもよい。この場合に、光源ユニット10は、凸レンズ13aと1つの凹レンズ13bとの間隔を変化させる機構(第2の機構)を有していてもよい。光源ユニット10は、第2の機構を用いて、凸レンズ13aと1つの凹レンズ13bとの間隔を変化させ、LED光源11から発せられた光が被検体ガラス100の界面で繰り返し反射されることにより生成される虚像Ivの位置を変化させることができる。
The
カメラユニット20は、光源ユニット10から発せられ、被検体ガラス100で反射された投影光Lを受光し、画像データDを生成するようになっている。カメラユニット20は、例えば、対物レンズ21、イメージセンサ22および画像処理部23を有する。対物レンズ21は、光源ユニット10から発せられ、被検体ガラス100で反射された投影光Lを集光するようになっている。対物レンズ21は、対物レンズ21の焦点距離をfとし、対物レンズ21の撮像倍率をmとするときに、(1+1/m)の位置に結像するようになっている。対物レンズ21のレイリーの2点分解能は、対物レンズ21の開口数の逆数に比例し、例えば、0.1’となっている。イメージセンサ22は、対物レンズ21で集光された光を受光面22Aで受光し、画像データDを生成するようになっている。イメージセンサ22は、例えば、CCD(Charge Coupled Device)イメージセンサである。対物レンズ21およびイメージセンサ22は、人の目の角度分解能(例えば、1’)よりも高い角度分解能を有することが好ましい。
The
画像処理部30は、カメラユニット20により得られた画像データDを処理し、検体ガラス100における投影光Lの反射の際に生じた二重像の視角θを求めるようになっている。画像処理部23は、画像データDに対して所定のデータ変換処理を行うことにより、イメージセンサ22の受光面22Aに生ずる、干渉縞パターンに対応した二重像の線像強度分布を求めるようになっている。画像処理部23は、さらに、得られた線像強度分布に基づき、被検体ガラス100における投影光Lの反射の際に生じた二重像の視角θを求めるようになっている。
The
画像処理部23は、二重像強度分布のモデル関数として、Youngの干渉実験装置200における2つのピンホール230aに等価な、フーリエ結像論上の空間周波数応答を用いることにより上記線像強度分布を導出し、導出した上記線像強度分布に基づいて二重像の視角θを求めるようになっている。画像処理部23は、上記モデル関数が示す線像強度分布のプロファイルを、画像データDに基づいて導出された線像強度分布のプロファイルにフィッティングさせることにより、二重像の視角θを求めるようになっている。
The image processing unit 23 uses the spatial frequency response on Fourier imaging theory equivalent to the two
画像処理部23は、上記モデル関数が示す二重像の線像強度分布の自己相関関数の積分値に基づいて、二重像の視角θを求めるようになっている。具体的には、画像処理部23は、上記モデル関数が示す二重像の線像強度分布の自己相関関数の積分値と、二重像の視角θとの相関を示す既知の関数を用いて、二重像の視角θを求めるようになっている。上記の既知の関数は、所定の視角θの範囲において一価の近似関数であり、例えば、後述するように、3次の多項式となっている。 The image processing unit 23 obtains the visual angle θ of the double image based on the integral value of the autocorrelation function of the linear image intensity distribution of the double image indicated by the model function. Specifically, the image processing unit 23 uses a known function indicating the correlation between the integral value of the autocorrelation function of the line image intensity distribution of the double image indicated by the model function and the visual angle θ of the double image. The viewing angle θ of the double image is to be obtained. The above-mentioned known function is a monovalent approximation function in the range of a predetermined viewing angle θ, and is, for example, a third-order polynomial as described later.
以下で、二重像強度分布のモデル関数、および近似関数について、詳細に説明する。 The model function of the dual image intensity distribution and the approximation function will be described in detail below.
[二重像強度分布のモデル関数]
まず、二重像強度分布のモデル関数について説明する。二重像分解能力は、一般には、レイリーの2点分解能で計測される。レイリーの2点分解能は、検出系のレンズの開口直径が光の波長と同程度の大きさになると、光の波動性によって、開口による回折効果によって分解能が制限されることを記述している。例えば、焦点距離f、開口半径wの望遠鏡では、レイリーの2点分解能εは、以下の式(1)で表される。従って、視角θ=1’を解像するときの条件は、ε=fθとすることにより導出される。
First, a model function of dual image intensity distribution will be described. Dual image resolution capabilities are generally measured at Rayleigh's two point resolution. Rayleigh's two-point resolution describes that when the aperture diameter of the lens of the detection system becomes as large as the wavelength of light, the wave nature of light limits the resolution by the diffractive effect of the aperture. For example, in a telescope with a focal length f and an aperture radius w, Rayleigh's two-point resolution ε is expressed by the following equation (1). Therefore, the condition for resolving the viewing angle θ = 1 ′ is derived by setting ε = fθ.
これにより、焦点距離fに無関係に、レイリーの2点分解能εに対応する開口半径wが求められる。つまり、開口半径wは、以下の式(2)で表される。
しかし、実際に合わせガラスで発生する二重像を計測する対物レンズ21の開口半径wを1.153mmよりも大きくしても、有効開口径は、1.15mm程度に制限されている。そのため、レイリーの2点分解能εは、1’以下にはならない。本出願人は、このような複雑な空間周波数応答を説明するために、二重像強度分布のモデル式を、図2のYoungの干渉実験装置200における干渉像の式を元に考案した。
However, even if the aperture radius w of the
Youngの干渉実験装置200は、図2に示したように、光源210から発せられる光の光路上に、レンズ220、ピンホールスクリーン230、レンズ240およびスクリーン250を、光源210側からこの順に備えている。ピンホールスクリーン230には、2つのピンホール230aが設けられている。干渉実験装置200では、2つのピンホール230aを通過した2つの光線がスクリーン250上で干渉し合い、スクリーン250上に干渉縞を発生させる。しかし、干渉実験装置200において、ピンホールスクリーン230が除去されると、スクリーン250上に干渉縞が発生しない。従って、2つのピンホール230aの機能の1つとして、波面分割があることが理解される。
As shown in FIG. 2, Young's
一方、本出願人が考案したモデル式を用いた二重像検査装置1でも、被検体ガラス100に光が照射されることにより、受光面22Aに二重像が発生する。ところが、二重像検査装置1には、Youngの干渉実験装置200における2つのピンホール230aに相当する構造物は存在しない。二重像検査装置1では、スリット12、被検体ガラス100内の繰り返し反射による波面分割、および対物レンズ21の開口半径が組み合わさることにより、Youngの干渉実験装置200における2つのピンホール230aに等価な、フーリエ結像論上の空間周波数応答が二重像強度分布のモデル関数として実現されている。
On the other hand, even in the double
二重像検査装置1における空間周波数応答は、以下の式(3)で表される。また、二重像検査装置1における空間周波数応答は、図3に示した二重像の干渉縞モデル300で説明される。干渉像モデル300は、LED光源11から発せられる光の光路上に、スリット12、レンズ13、投影レンズ14、被検体ガラス100、仮想ピンホールスクリーン310、対物レンズ21および受光面22Aを、LED光源11側から順に有する。仮想ピンホールスクリーン310は、2つのピンホール230aに対応する2つの仮想ピンホールH1を有する。
The spatial frequency response in the dual
ここで、J1は、1次の第1種ベッセル関数である。λは、LED光源11から発せられた光の平均波長である。zは、像面距離である。δは、仮想ピンホールH1の直径である。{x,y}は、受光面22A上の座標である。x0は、二重像の像高である。h'は、互いに隣接する2つの仮想ピンホールH1の間隔hの複数倍である。一方の仮想ピンホールH1には、LED光源11から発せられ、被検体ガラス100を透過した光L1が入射する。他方の仮想ピンホールH1には、被検体ガラス100の裏面で反射された後、受光面22A側に射出された光L2が入射する。μ12は、干渉縞の可視度であり、0から1の間の値を取る。α12は、光L2の位相であり、−π〜πの範囲の値を取る。
Here, J 1 is a first-order Bessel function of the first kind. λ is an average wavelength of the light emitted from the LED light source 11. z is an image plane distance. δ is the diameter of the virtual pinhole H1. {X, y} are coordinates on the
ところで、広く知られているYoungの干渉縞の式では、右辺の点像は1個となっている。これは、2つのピンホール230aで回折される光の干渉縞の観測位置が無限遠となっているからである。一方で、二重像検査装置1における空間周波数応答では、LED光源11の虚像Ivの位置が有限となっている。その結果、受光面22Aにおいて、二重像に干渉縞が重畳されている様子が観測される。また、干渉像型モデル300では、Youngの干渉実験装置200における2つのピンホール230aに対応するものとして、2つの仮想ピンホールH1が設けられている。2つの仮想ピンホールH1の直径δは、以下の式(4)で規定されている。被検体ガラス100の厚さtを4mm、被検体ガラス100に対する光L1の入射角ψを60°としたとき、ピンホールH1の直径δは、2.77mmとなる。このとき、下記の条件を適用することにより、二重像検査装置1における空間周波数応答から得られる二重像の線像強度分布が、Youngの干渉実験装置200により得られた二重像の線像強度分布に良好に近似することを本出願人は見出した。
By the way, in the widely known Young's interference fringe equation, the point image on the right side is one. This is because the observation position of the interference fringes of the light diffracted by the two
−条件―
スリット12の幅a:0.01mm
対物レンズ21の開口半径w:12mm
2つの仮想ピンホールH1の間隔h:2.77mm
-Condition-
Width a of slit 12: 0.01 mm
Aperture radius w of objective lens 21: 12 mm
Distance h between two virtual pinholes H1: 2.77 mm
図4は、二重像検査装置1における空間周波数応答を用いて算出した二重像の線像強度分布を等高線で表したものである。図4では、二重像の線像強度分布が被検体ガラス100の傾き方向に干渉縞が重畳している。図5は、二重像の線像強度分布の垂直断面の計測値および計算値をグラフで表したものである。図5において、計測値は点で表されており、計算値は実線で表されている。計測値は、イメージセンサ22から得られた画像データに基づいて導出された強度分布の垂直断面の値である。計算値は、二重像検査装置1における空間周波数応答を用いて算出した二重像の線像強度分布の値である。なお、二重像の視角θは約0.75’に仮定されている。図5の横軸は、像面(受光面22A)の空間座標であり、その単位は、イメージセンサ22の画素ピッチであり、具体的には、4.65μmである。図5に示したように、二重像検査装置1における空間周波数応答から得られる二重像の線像強度分布のプロファイルが、イメージセンサ22から得られた画像データに基づいて導出された線像強度分布のプロファイルに良好に近似することがわかる。なお、図6、図7、図8、図9、図10に、二重像の視角θが0.866444’、0.649833’、0.433222’、0.216611’、0’に仮定されているときに得られた二重像の線像強度分布を示す。図10において、二重像が分離されていることがわかる。
FIG. 4 is a diagram showing the line image intensity distribution of the double image calculated using the spatial frequency response in the double
二重像検査装置1における空間周波数応答は、以下のようにして導出される。小瀬輝次著の「フーリエ結像論」(共立出版 2013年復刊)に記載の複ピンホール光学系のOTF(Optical Transfer Function)の概念を踏襲して、ガラス表面とガラス裏面とによる多重繰り返し反射光のうち、表面反射光が入射する仮想ピンホールH1と、一次の裏面反射光が二重像の視角θ分だけ傾いて入射する仮想ピンホールH1を考える。このとき、仮想ピンホールH1の直径δを、表面反射光の光軸と、一次の裏面反射光の光軸との距離と等しくする。その上で、対物レンズ21の開口半径wが、仮想ピンホールH1の間隔hの複数個分となっていると規定する。最後に、複数の仮想ピンホールH1の自己相関関数から、二重像検査装置1における空間周波数応答を導出する。このようにして求めた空間周波数応答のフーリエ変換は、干渉縞の重畳した二重像の線像強度分布になるはずである。
The spatial frequency response in the dual
仮想ピンホールH1の直径δの大きさが、点像サイズを決定する。レイリーの2点分解の意味では、点像サイズの1/2に対応する角度分解能が限界値を与えるはずである。仮想ピンホールH1の直径δが2.77mmであるとき、対物レンズ21の焦点距離fを135mmとすると、点像のエアリー直径dは、1.22λz/δ=65.4μmとなる。ここで、像面距離zを、(1+1/m)fと近似した。図5では、画素ピッチは4.65μmであるから、エアリー直径dは、14画素ピッチに対応し、95番目の画素から109番目の画素までの範囲に相当している。一方、予想される二重像の視角θは0.75’であるから、像高さfθは29.4μmである。ここで、繰り返し反射光により干渉縞を重畳させないときの二重像強度分布を図11に示す。焦点距離を50mmにした場合で、繰り返し反射光により干渉縞を重畳させないときの二重像強度分布を図12に示す。これにより、長焦点距離の方が、角度分解能を上げることができる点で有利あることがわかる。ただし、図5では、多光束干渉効果、コンボリューションおよび収差が考慮されており、一方で、図11、図12では、これらが考慮されていない。
The size of the diameter δ of the virtual pinhole H1 determines the point image size. In the sense of Rayleigh's two-point resolution, an angular resolution corresponding to one-half of the point spread should provide a limit. Assuming that the focal distance f of the
上述したように、式(3)で示した空間周波数応答を、計測値にフィッティングさせることは可能である。しかし、式(3)で示した空間周波数応答では、パラメータの自由度が多いので、計測値にフィッティングさせるには、膨大な計算が必要になる。そこで、本出願人は、上式(3)で示した空間周波数応答を計測値にフィッティングさせる代わりに、以下に示す近似関数を用いて、二重像の視角θを導出する手法を見出した。 As described above, it is possible to fit the spatial frequency response shown in equation (3) to the measurement value. However, in the spatial frequency response shown in the equation (3), there are many degrees of freedom of parameters, and therefore, a great deal of calculation is required to fit the measurement value. Therefore, the present applicant has found a method for deriving the viewing angle θ of a double image by using an approximation function shown below, instead of fitting the spatial frequency response shown in the above equation (3) to the measurement value.
[近似関数]
次に、近似関数について説明する。図13は、上式(3)で示した空間周波数応答を用いて得られた二重像の線像強度分布の自己相関関数の積分値と、二重像の視角θとの関係を表したものである。図13には、二重像の視角θを横軸に、上記自己相関関数の積分値を縦軸にしたときのグラフが示されている。図13から、上式(3)で示した空間周波数応答が示す二重像の線像強度分布の自己相関関数の積分値に基づいて、二重像の視角θを求めることができる。図13に示したグラフは、二重像の視角θが大きくなるにつれて、上記自己相関関数の積分値も単調に増加する一価の近似関数で表される。この近似関数は、例えば、図13に示したように、Y=aX3+bX+cという3次の多項式で表される。二重像の視角θが、少なくとも、0°から0.9’までの範囲においては、この近似関数は、一価の近似関数となっている。従って、この既知の関数である近似関数を用いて、1’以下の二重像の視角θを求めることができる。
[Approximate function]
Next, the approximation function will be described. FIG. 13 shows the relationship between the integral value of the autocorrelation function of the line image intensity distribution of the double image obtained using the spatial frequency response shown by the above equation (3) and the visual angle θ of the double image. It is a thing. FIG. 13 shows a graph in which the viewing angle θ of the double image is on the horizontal axis and the integral value of the autocorrelation function is on the vertical axis. The visual angle θ of the double image can be determined from FIG. 13 based on the integral value of the autocorrelation function of the linear image intensity distribution of the double image indicated by the spatial frequency response indicated by the above equation (3). The graph shown in FIG. 13 is represented by a single-valued approximation function in which the integral value of the autocorrelation function monotonously increases as the viewing angle θ of the double image increases. For example, as shown in FIG. 13, this approximate function is expressed by a third-order polynomial of Y = aX 3 + bX + c. In the range of the viewing angle θ of the double image at least from 0 ° to 0.9 ′, this approximate function is a single-valued approximate function. Therefore, the viewing angle θ of the double image of 1 ′ or less can be determined using the approximation function that is this known function.
[処理手順]
次に、二重像検査装置1における処理手順の一例について説明する。図14は、二重像検査装置1における処理手順の一例を表したものである。投影光Lが光源ユニット10から発せられると、被検体ガラス100で反射され、受光面22Aに入射する。受光面22Aに入射した光は受光面22Aで受光され、画像データDが生成される。画像データDは、画像処理部30に取り込まれる(ステップS101)。画像データDは、画像処理部30によって、以下のように処理される。具体的には、まず、窓関数を用いて、画像データDのノイズが除去される(ステップS102)。次に、画像データDが規格化されたのち、複素フーリエ変換される(ステップS103)。これにより、イメージセンサ22の受光面22Aに生ずる、干渉縞パターンに対応した二重像の線像強度分布に対応する複素瞳関数が求められる。次に、得られた複素瞳関数の絶対値の2乗に対して逆フーリエ変換が行われる(ステップS104,105)。これにより、線像強度分布の自己相関関数が得られる。Bracewellの等価幅計算が行われる(ステップS106)。つまり、得られた自己相関関数の積分値が求められる。最後に、自己相関関数の積分値と、二重像の視角θとの相関を示す既知の関数を用いて、二重像の視角θが求められる(ステップS107)。なお、ステップS104〜S107を実施する代わりに、上記モデル関数が示す線像強度分布のプロファイルを、画像データDに基づいて導出された線像強度分布のプロファイルにフィッティングさせることにより、二重像の視角θが求められてもよい。
[Procedure]
Next, an example of the processing procedure in the double-
[効果]
次に、特許文献1の記載の発明と対比しつつ、二重像検査システム1の効果について説明する。
[effect]
Next, the effects of the double-
特許文献1の記載によれば、入射角ψを固定し、曲面ガラスによって発生する二重像の視角η1(式(5))と、楔角αを有する平板ガラスによって発生する二重像の視角η2(式(6))とが、互いに等しいと設定することにより、二重像を解消することの可能なガラスの厚さt、ガラスの曲率半径Rcおよび楔角αを、式(7)から導出することができる、ようである。
しかしながら、二重像の視角が0°になるように楔角αを設定したり、合わせガラスを平板ガラスに置き換えたりした場合であっても、二重像が発生することがあることを本出願人は発見した。これは、特許文献1では、式(6)の導出の際に、虚像距離を無限大に仮定しているためで、虚像距離が有限である場合には、平板ガラスの楔角αが無くても二重像は有限になるからである。虚像距離をd、表面反射光の光軸と裏面反射光の光軸との間隔をhとするとき、以下の式(8)が得られる。従って、有限虚像距離における二重像の解消に必要な楔角αは、以下の式(9)で表される。式(9)において、d→∞とすれば、式(7)に帰着する。
However, even if the depression angle α is set so that the visual angle of the double image is 0 ° or the laminated glass is replaced with flat glass, the double image may be generated. People have found. This is because, in
一方、二重像検査システム1では、二重像の線像強度分布のモデル関数として、Youngの干渉実験装置200における2つのピンホール230aに等価な、フーリエ結像論上の空間周波数応答が用いられる。これにより、Youngの干渉実験装置200における2つのピンホール230aに相当する構造物が存在しないにもかかわらず、イメージセンサ22の受光面22Aに生ずる、干渉縞パターンに対応した二重像の線像強度分布を導出することができる。さらに、得られた線像強度分布に基づき、被検体ガラス100における投影光Lの反射の際に生じた二重像の視角θが求められる。このように、二重像検査システム1では、幾何光学的なアプローチから二重像の視角θを導出するのではなく、波動光学的なアプローチから二重像の視角θを導出している。これにより、二重像を精密に検査することができる。
On the other hand, in the double
また、二重像検査システム1では、上式(3)で示した空間周波数応答を計測値にフィッティングさせる代わりに、上記モデル関数が示す二重像の線像強度分布の自己相関関数の積分値に基づいて、二重像の視角θが求められる。これにより、式(3)で示した空間周波数応答を、計測値にフィッティングさせる際に要する膨大な計算をすることなく、軽い計算量で、二重像の視角θを導出することができる。
Further, in the double
1…二重像検査システム、10…光源ユニット、11…LED光源、12…スリット、13…レンズ、13a…凸レンズ、13b…凹レンズ、14…投影レンズ、20…カメラユニット、21…対物レンズ、22…イメージセンサ、30…画像処理装置、100…被検体ガラス、200…干渉実験装置、210…光源、220…レンズ、230…ピンホールスクリーン、230a…ピンホール、24…レンズ、250…スクリーン、300…干渉縞モデル、310…仮想ピンホールスクリーン、a…スリット幅、D…画像データ、f…焦点距離、H1…仮想ピンホール、h…互いに隣接する2つの仮想ピンホールの間隔、Iv…虚像、J1…1次の第1種ベッセル関数、L…投影光、L1,L2…光、w…対物レンズの開口半径、λ…LED光源から発せられた光の平均波長、{x,y}…受光面上の座標、x0…二重像の像高、z…像面距離、α…楔角、α12…光L2の位相、δ…仮想ピンホールの直径、θ…二重像の視角、μ12…干渉縞の可視度。
DESCRIPTION OF
Claims (7)
前記光源ユニットから発せられ、前記被検体ガラスの界面で繰り返し反射された光を集光する対物レンズと、前記対物レンズにより集光された光を受光するイメージセンサとを含むカメラユニットと、
前記カメラユニットにより得られた画像データを処理する処理部と
を備え、
前記処理部は、
前記画像データに対して所定のデータ変換処理を行うことにより、前記イメージセンサの受光面に生ずる干渉縞パターンに対応した二重像の線像強度分布を求め、
得られた前記線像強度分布に基づき、前記被検体ガラスにおける前記投影光の反射の際に生じた二重像の視角を求める
二重像検査システム。 A light source unit including a light source, a slit illuminated to the light source, and a projection optical system for projecting linear light emitted from the light source via the slit as projection light onto a subject glass;
A camera unit including an objective lens for collecting light emitted from the light source unit and reflected repeatedly at the interface of the object glass , and an image sensor for receiving the light collected by the objective lens;
A processing unit that processes image data obtained by the camera unit;
The processing unit is
By performing predetermined data conversion processing on the image data, a line image intensity distribution of a double image corresponding to the interference fringe pattern generated on the light receiving surface of the image sensor is obtained.
The double image inspection system which calculates | requires the visual angle of the double image which arose at the time of reflection of the said projection light in the said object glass based on the obtained said linear image intensity distribution.
請求項1に記載の二重像検査システム。 The processing unit obtains the line image intensity distribution by using a spatial frequency response on Fourier imaging theory equivalent to a double pin hole in Young's interference experiment as a model function of a line image intensity distribution of a double image. Item 2. A dual image inspection system according to item 1.
請求項2に記載の二重像検査システム。
λ:前記光源から発せられた光の平均波長
z:像面距離
δ:前記ピンホールに対応する仮想ピンホールの直径
{x,y}:前記受光面上の座標
x0:前記二重像の像高
h':互いに隣接する2つの前記仮想ピンホールの間隔hの複数倍
μ12:干渉縞の可視度
α12は、前記被検体ガラスの表面および裏面で繰り返し反射された後、前記受光面側に射出された光の位相 The dual image inspection system according to claim 2, wherein the model function is expressed by the following equation (1).
請求項3に記載の二重像検査システム。 The processing unit obtains a viewing angle of the double image by fitting a profile of a linear image intensity distribution represented by the model function to a profile of a linear image intensity distribution derived based on the image data. Dual image inspection system as described in.
請求項3に記載の二重像検査システム。 The double image inspection system according to claim 3, wherein the processing unit obtains a viewing angle of the double image based on an integral value of an autocorrelation function of a line image intensity distribution of the double image indicated by the model function.
請求項5に記載の二重像検査システム。 The processing unit uses the known function that indicates the correlation between the integrated value of the autocorrelation function of the line image intensity distribution of the double image indicated by the model function and the viewing angle of the double image. The dual image inspection system according to claim 5, wherein a visual angle of
請求項6に記載の二重像検査システム。 The dual image inspection system according to claim 6, wherein the known function is a cubic polynomial.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014250399A JP6503568B2 (en) | 2014-12-10 | 2014-12-10 | Dual image inspection system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014250399A JP6503568B2 (en) | 2014-12-10 | 2014-12-10 | Dual image inspection system |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2016114362A JP2016114362A (en) | 2016-06-23 |
JP6503568B2 true JP6503568B2 (en) | 2019-04-24 |
Family
ID=56139832
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2014250399A Active JP6503568B2 (en) | 2014-12-10 | 2014-12-10 | Dual image inspection system |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6503568B2 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111323209A (en) * | 2020-03-13 | 2020-06-23 | 江苏泽景汽车电子股份有限公司 | HUD stray light automatic test system and test method |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102016114485A1 (en) | 2016-08-04 | 2018-02-08 | Isra Surface Vision Gmbh | Device and method for determining a double image angle and / or a viewing angle |
KR101968689B1 (en) * | 2017-04-18 | 2019-04-16 | 자동차부품연구원 | Measurement system and method for ghost images on head-up display |
CN107576483B (en) * | 2017-09-13 | 2019-03-19 | 上海伟世通汽车电子系统有限公司 | Camera calibrated device and method thereof for automobile HUD Meter Test |
CN107702734B (en) * | 2017-09-15 | 2020-06-16 | 大连理工大学 | Spatial phase modulation type laser interference measuring instrument and method |
CN109855845B (en) * | 2019-03-27 | 2022-05-24 | 广东技术师范大学 | Binocular eye lens measurement vehicle-mounted HUD virtual image distance and correction method |
CN110006634B (en) * | 2019-04-15 | 2020-12-18 | 北京京东方光电科技有限公司 | Viewing field angle measuring method, viewing field angle measuring device, display method and display equipment |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3707331A1 (en) * | 1987-03-07 | 1988-09-15 | Zeiss Carl Fa | INTERFEROMETER FOR MEASURING OPTICAL PHASE DIFFERENCES |
JPH04346048A (en) * | 1991-05-23 | 1992-12-01 | Ricoh Co Ltd | Method and device for receiving light in scanning line position detection device |
JPH07218233A (en) * | 1994-02-08 | 1995-08-18 | Ricoh Co Ltd | Method and apparatus for measuring surface to be detected |
JPH09218132A (en) * | 1996-02-15 | 1997-08-19 | Dainippon Screen Mfg Co Ltd | Inspection and measuring device for optical element |
-
2014
- 2014-12-10 JP JP2014250399A patent/JP6503568B2/en active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111323209A (en) * | 2020-03-13 | 2020-06-23 | 江苏泽景汽车电子股份有限公司 | HUD stray light automatic test system and test method |
Also Published As
Publication number | Publication date |
---|---|
JP2016114362A (en) | 2016-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6503568B2 (en) | Dual image inspection system | |
JP5337050B2 (en) | A method for automatically and quantitatively analyzing the distortion of molded glass for vehicles using reflected optical images. | |
EP2805670B1 (en) | Improvements in and relating to ophthalmoscopes | |
JP6494205B2 (en) | Wavefront measuring method, shape measuring method, optical element manufacturing method, optical device manufacturing method, program, wavefront measuring apparatus | |
JP6529830B2 (en) | Calculation method of height map of transparent solid with inclined surface or curved surface | |
US8045181B2 (en) | Inspection system and method with multi-image phase shift analysis | |
JP6112909B2 (en) | Shape measuring device and shape measuring method using Shack-Hartmann sensor | |
US9239237B2 (en) | Optical alignment apparatus and methodology for a video based metrology tool | |
JP2015505039A (en) | Non-contact surface shape evaluation using modulated light | |
JP2018515747A5 (en) | ||
JP5748414B2 (en) | Cylindrical shape measurement method | |
JP2018515747A (en) | Method and apparatus for determining surface data and / or measurement data relating to the surface of an at least partly transparent object | |
TW201346215A (en) | Image forming optical system, imaging apparatus, profile measuring apparatus, structure manufacturing system and structure manufacturing method | |
AU2010246485B2 (en) | Triangulation camera device and triangulation imaging method | |
JP6574674B2 (en) | Coating inspection device and inspection method | |
TWI553291B (en) | System for measuring transparent object by fringe projection | |
JP7076586B2 (en) | Imaging device | |
JP6429503B2 (en) | Measuring device, measuring method, optical element processing apparatus, and optical element | |
JP2013019786A (en) | Inside tube wall inspection apparatus and inside tube wall inspection method | |
US20220334068A1 (en) | Method for evaluating the optical quality of a delineated region of a glazing | |
JP6196841B2 (en) | Transmitted wavefront measuring apparatus and transmitted wavefront measuring method | |
WO2016013518A1 (en) | Image-forming optical system, illumination device, observation device, and wavefront recovery element | |
JP5325481B2 (en) | Measuring method of optical element and manufacturing method of optical element | |
US11536939B2 (en) | Method for digitally correcting an optical image of a sample by means of a microscope, and microscope | |
KR102180648B1 (en) | Apparatus and method for 3-dimensional tomographic inspection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20170926 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20180629 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20180717 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20180914 |
|
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: 20190219 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20190228 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6503568 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |