JP2022007906A - Three-dimensional measuring device - Google Patents
Three-dimensional measuring device Download PDFInfo
- Publication number
- JP2022007906A JP2022007906A JP2020198371A JP2020198371A JP2022007906A JP 2022007906 A JP2022007906 A JP 2022007906A JP 2020198371 A JP2020198371 A JP 2020198371A JP 2020198371 A JP2020198371 A JP 2020198371A JP 2022007906 A JP2022007906 A JP 2022007906A
- Authority
- JP
- Japan
- Prior art keywords
- pattern
- dimensional
- light
- light source
- image pickup
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000005259 measurement Methods 0.000 claims abstract description 109
- 230000010355 oscillation Effects 0.000 claims abstract description 21
- 238000000034 method Methods 0.000 claims description 62
- 230000010363 phase shift Effects 0.000 claims description 40
- 230000000737 periodic effect Effects 0.000 claims description 25
- 238000005070 sampling Methods 0.000 claims description 5
- 238000003384 imaging method Methods 0.000 abstract description 25
- 230000009467 reduction Effects 0.000 abstract description 3
- 239000013598 vector Substances 0.000 description 103
- 230000003287 optical effect Effects 0.000 description 30
- 230000005484 gravity Effects 0.000 description 24
- 238000010586 diagram Methods 0.000 description 20
- 239000000470 constituent Substances 0.000 description 18
- 239000002184 metal Substances 0.000 description 14
- 229910052751 metal Inorganic materials 0.000 description 14
- 239000004065 semiconductor Substances 0.000 description 13
- 239000000758 substrate Substances 0.000 description 10
- 239000011159 matrix material Substances 0.000 description 7
- 238000013461 design Methods 0.000 description 6
- 238000007689 inspection Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 230000002093 peripheral effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000001839 endoscopy Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000000691 measurement method Methods 0.000 description 3
- 239000000523 sample Substances 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 230000007480 spreading Effects 0.000 description 3
- 101710195649 Gamma-crystallin M1 Proteins 0.000 description 2
- 101710195647 Gamma-crystallin M2 Proteins 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005530 etching Methods 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 230000001678 irradiating effect Effects 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 229910001218 Gallium arsenide Inorganic materials 0.000 description 1
- 229910021417 amorphous silicon Inorganic materials 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 150000004767 nitrides Chemical class 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Semiconductor Lasers (AREA)
Abstract
Description
本開示は、三次元計測装置に関する。 The present disclosure relates to a three-dimensional measuring device.
従来の三次元計測法として、例えば特許文献1に記載の手法がある。この特許文献1の手法では、ランダムなドットパターンを被計測物に照射し、2台のカメラで同位置のドットパターンをそれぞれ撮像する。そして、2つのドットパターンの視差に基づいて、三角測量の原理で被計測物の三次元計測を実施する。
As a conventional three-dimensional measurement method, for example, there is a method described in
また、例えば特許文献2に記載の手法は、位相シフト法を用いた計測手法である。この特許文献2の手法では、格子パターンが投影される基準面を有する基準平板を用意し、当該基準平板をステージによって法線方向に平行移動させる。基準面に投影された格子パターンの画像と、被計測物に投影された格子パターンの画像とを撮像し、格子パターンの位相と空間座標とを対応付けるテーブルを用いて被計測物の空間座標を算出する。
Further, for example, the method described in
上述した特許文献1の手法では、光源としてプロジェクタを用いており、特許文献2の手法では、光源としてLEDアレイを用いている。このため、三次元計測装置が比較的大型化してしまうという問題があった。撮像装置としては、例えば1mm角以下の超小型カメラも開発されているため、三次元計測装置を全体として小型化するためには、光源の小型化が重要となっている。三次元計測装置の全体を小型化できれば、例えば口腔検査、内視鏡検査、管の内部や壁の隙間といった狭小箇所の検査、家具や装置等の床下からの検査といった用途への適用や、ハンディタイプの三次元計測装置の構築が可能となると考えられる。また、三次元計測装置への光源の適用にあたっては、計測精度の向上の観点から、出力される光のノイズや歪みが抑えられた光源であることが好ましい。
In the method of
本開示は、上記課題の解決のためになされたものであり、装置の小型化により適用範囲の拡大が図られ、且つ計測精度の向上が図られる三次元計測装置を提供することを目的とする。 The present disclosure has been made to solve the above problems, and an object of the present disclosure is to provide a three-dimensional measuring device capable of expanding the applicable range by downsizing the device and improving the measurement accuracy. ..
本開示の一側面に係る三次元計測装置は、被計測物に所定パターンを有する計測光を照射する一又は複数の光源部と、計測光が照射された被計測物を撮像する一又は複数の撮像部と、撮像部による撮像結果に基づいて被計測物の三次元形状を計測する計測部と、を備え、光源部は、M点発振のS-iPMSELによって構成されている。 The three-dimensional measuring device according to one aspect of the present disclosure includes one or a plurality of light source units that irradiate an object to be measured with measurement light having a predetermined pattern, and one or a plurality of light sources that image the object to be measured irradiated with the measurement light. It includes an image pickup unit and a measurement unit that measures the three-dimensional shape of the object to be measured based on the image pickup result by the image pickup unit, and the light source unit is configured by S-iPMSEL of M point oscillation.
この三次元計測装置では、光源部がM点発振のS-iPMSELによって構成されている。S-iPMSELは、基本層と、基本層とは屈折率が異なる複数の異屈折率領域とを有する位相変調層を有し、各異屈折率領域の重心位置が出力光像に応じて仮想的な正方格子の格子点位置からずれている。S-iPMSELは、例えば針先ほどのサイズで構成され、位相変調層が設けられた基板の主面に垂直な方向或いは傾斜した方向に2次元的なパターンの光像を出力できる。したがって、S-iPMSELを光源とすることで、三次元計測装置の全体の小型化を実現でき、装置の適用範囲の拡大が図られる。また、M点発振のS-iPMSELを用いることで、所望の2次元的なパターンの光像とは異なる0次光(位相変調されない回折波成分)の出力をなくすことができる。これにより、0次光によるノイズや歪みのないパターンの計測光を被計測物に照射することが可能となり、計測精度の向上が図られる。 In this three-dimensional measuring device, the light source unit is configured by S-iPMSEL of M point oscillation. The S-iPMSEL has a phase modulation layer having a basic layer and a plurality of different refractive index regions having different refractive indexes from the basic layer, and the position of the center of gravity of each different refractive index region is virtual according to the output light image. It deviates from the grid point position of a square grid. The S-iPMSEL can output a two-dimensional pattern of light images in a direction perpendicular to or inclined to the main surface of a substrate provided with a phase modulation layer, for example, having a size similar to that of a needle tip. Therefore, by using the S-iPMSEL as a light source, the overall size of the three-dimensional measuring device can be reduced, and the applicable range of the device can be expanded. Further, by using S-iPMSEL of M-point oscillation, it is possible to eliminate the output of 0th-order light (diffraction wave component that is not phase-modulated) different from the optical image of a desired two-dimensional pattern. As a result, it becomes possible to irradiate the object to be measured with the measurement light having a pattern without noise or distortion due to the 0th-order light, and the measurement accuracy can be improved.
三次元計測装置は、単体の光源部と、複数の撮像部と、を含み、計測光の前記所定パターンは、ドットパターン、ストライプパターン、及び格子パターンのいずれかからなる周期パターンであり、計測部は、周期パターンを用いたアクティブステレオ法に基づいて被計測物の三次元形状を計測してもよい。この場合、テクスチャの少ない画像を用いた三次元計測、暗部での三次元計測が可能となる。 The three-dimensional measuring device includes a single light source unit and a plurality of imaging units, and the predetermined pattern of the measurement light is a periodic pattern including any of a dot pattern, a stripe pattern, and a grid pattern, and is a measurement unit. May measure the three-dimensional shape of the object to be measured based on the active stereo method using a periodic pattern. In this case, three-dimensional measurement using an image with few textures and three-dimensional measurement in a dark part become possible.
三次元計測装置は、単体の前記光源部と、複数の撮像部と、を含み、計測光の所定パターンは、ランダムドットパターンであり、計測部は、ランダムドットパターンを用いたアクティブステレオ法に基づいて被計測物の三次元形状を計測する。この場合、周期的なドットパターンではなく、ランダムドットパターンを用いることで、ドットパターンの同一点を異なる撮像部で撮像した際の誤認識を抑制できる。 The three-dimensional measuring device includes the single light source unit and a plurality of imaging units, a predetermined pattern of measurement light is a random dot pattern, and the measuring unit is based on an active stereo method using a random dot pattern. Measure the three-dimensional shape of the object to be measured. In this case, by using a random dot pattern instead of a periodic dot pattern, it is possible to suppress erroneous recognition when the same point of the dot pattern is imaged by different imaging units.
三次元計測装置は、単体の光源部と、複数の撮像部と、を含み、計測光の所定パターンは、均一な密度を有するパターンであり、計測部は、均一な密度を有するパターンを用いたアクティブステレオ法に基づいて被計測物の三次元形状を計測してもよい。S-iPMSELからの出射光は、レーザ光であるため、散乱光中にスペックルが生じ得る。したがって、均一な密度を有するパターンを用いた場合であっても、計測光のパターンにランダムドットパターンが形成される。このランダムドットパターンを用いることで、ドットパターンの同一点を異なる撮像部で撮像した際の誤認識を抑制できる。 The three-dimensional measuring device includes a single light source unit and a plurality of imaging units, and a predetermined pattern of measurement light is a pattern having a uniform density, and the measuring unit uses a pattern having a uniform density. The three-dimensional shape of the object to be measured may be measured based on the active stereo method. Since the light emitted from the S-iPMSEL is a laser light, speckle may occur in the scattered light. Therefore, even when a pattern having a uniform density is used, a random dot pattern is formed in the pattern of the measurement light. By using this random dot pattern, it is possible to suppress erroneous recognition when the same point of the dot pattern is imaged by different imaging units.
三次元計測装置は、複数の光源部と、単体の撮像部と、を含み、計測光の所定パターンは、グレイコードパターンであり、計測部は、グレイコードパターンを用いた三角測量法に基づいて被計測物の三次元形状を計測してもよい。グレイコードパターンのパターン数は、撮像部の画素数に対して少数でよいため、グレイコードパターンを有する計測光の照射は、少数の光源部によって実現できる。グレイコードを用いる場合、隣接する画素のハミング距離が1となり、ビット列を復元する際にビットエラーが生じたとしても、その誤差が1に収まる。すなわち、グレイコードでは、ノイズに強い符号が得られる。 The three-dimensional measuring device includes a plurality of light source units and a single image pickup unit, and a predetermined pattern of measurement light is a Gray code pattern, and the measurement unit is based on a triangulation method using a Gray code pattern. The three-dimensional shape of the object to be measured may be measured. Since the number of patterns of the Gray code pattern may be small with respect to the number of pixels of the image pickup unit, irradiation of the measurement light having the Gray code pattern can be realized by a small number of light source units. When the Gray code is used, the Hamming distance of adjacent pixels is 1, and even if a bit error occurs when restoring the bit string, the error is within 1. That is, in the Gray code, a code that is resistant to noise can be obtained.
三次元計測装置は、複数の光源部と、単体の撮像部と、を含み、計測光の所定パターンは、正弦波状のストライプパターンであり、計測部は、正弦波状のストライプパターンを用いた位相シフト法に基づいて被計測物の三次元形状を計測してもよい。この場合、計測した位相を高さ換算することで、正弦波状のストライプパターンのピッチよりも小さい間隔で被計測物の高さを計測することができる。 The three-dimensional measuring device includes a plurality of light source units and a single image pickup unit, and a predetermined pattern of measurement light is a sinusoidal stripe pattern, and the measurement unit is a phase shift using a sinusoidal stripe pattern. The three-dimensional shape of the object to be measured may be measured based on the method. In this case, by converting the measured phase into height, the height of the object to be measured can be measured at intervals smaller than the pitch of the sinusoidal stripe pattern.
複数の光源部は、互いに周期の異なる正弦波状のストライプパターンをそれぞれ出力してもよい。位相シフト法では、位相2πにおける不連続性が課題となっている。これに対し、互いに周期の異なる正弦波状のストライプパターンを用いることで、位相2πにおける不連続性を改善することが可能となり、少数のパターンで高精度の三次元計測を実現できる。 The plurality of light source units may output sinusoidal stripe patterns having different periods from each other. In the phase shift method, the discontinuity at the phase 2π is a problem. On the other hand, by using sinusoidal stripe patterns having different periods from each other, it is possible to improve the discontinuity in the phase 2π, and it is possible to realize highly accurate three-dimensional measurement with a small number of patterns.
三次元計測装置は、複数の光源部と、単体の撮像部と、を含み、計測光の所定パターンは、正弦波状のストライプパターンであり、計測部は、正弦波状のストライプパターンを用いたサンプリングモアレ法に基づいて被計測物の三次元形状を計測してもよい。この場合、更に少数のパターンで高精度の三次元計測を実現できる。 The three-dimensional measuring device includes a plurality of light source units and a single image pickup unit, and a predetermined pattern of measurement light is a sinusoidal stripe pattern, and the measurement unit is a sampling moire using a sinusoidal stripe pattern. The three-dimensional shape of the object to be measured may be measured based on the method. In this case, highly accurate three-dimensional measurement can be realized with a smaller number of patterns.
三次元計測装置は、複数の光源部と、単体の撮像部と、を含み、計測光の所定パターンは、正弦波状のストライプパターンと、ランダムドットパターンとを重畳させた重畳パターンであり、計測部は、重畳パターンを用いた位相シフト法に基づいて被計測物の三次元形状を計測してもよい。位相シフト法では、位相2πにおける不連続性が課題となっている。これに対し、ランダムドットパターンを用いることで、位相2πにおける不連続性を改善することが可能となり、少数のパターンで高精度の三次元計測を実現できる。 The three-dimensional measuring device includes a plurality of light source units and a single image pickup unit, and a predetermined pattern of measurement light is a superimposed pattern in which a sinusoidal stripe pattern and a random dot pattern are superimposed. May measure the three-dimensional shape of the object to be measured based on the phase shift method using the superimposition pattern. In the phase shift method, the discontinuity at the phase 2π is a problem. On the other hand, by using the random dot pattern, it is possible to improve the discontinuity in the phase 2π, and it is possible to realize highly accurate three-dimensional measurement with a small number of patterns.
三次元計測装置は、複数の光源部と、単体の撮像部と、を含み、計測光の所定パターンは、正弦波状のストライプパターンと、グレイコードパターンとを含み、計測部は、正弦波状のストライプパターンを用いた位相シフト法、及びグレイコードパターンを用いた三角測量法に基づいて被計測物の三次元形状を計測してもよい。位相シフト法では、位相2πにおける不連続性が課題となっている。これに対し、グレイコードを用いることで、位相2πにおける不連続性を改善することが可能となり、少数のパターンで高精度の三次元計測を実現できる。 The three-dimensional measuring device includes a plurality of light source units and a single image pickup unit, a predetermined pattern of measurement light includes a sinusoidal stripe pattern and a Gray code pattern, and the measurement unit includes sinusoidal stripes. The three-dimensional shape of the object to be measured may be measured based on a phase shift method using a pattern and a triangular survey method using a Gray code pattern. In the phase shift method, the discontinuity at the phase 2π is a problem. On the other hand, by using the Gray code, it is possible to improve the discontinuity in the phase 2π, and it is possible to realize highly accurate three-dimensional measurement with a small number of patterns.
三次元計測装置は、複数の光源部と、複数の撮像部と、を含み、計測光の所定パターンは、正弦波状のストライプパターンであり、計測部は、正弦波状のストライプパターンを用いた位相シフト法及びアクティブステレオ法に基づいて被計測物の三次元形状を計測してもよい。この場合、計測した位相を高さ換算することで、正弦波状のストライプパターンのピッチよりも小さい間隔で被計測物の高さを計測することができる。また、位相シフト法では、位相2πにおける不連続性が課題となっている。これに対し、複数の撮像部を用いたアクティブステレオ法を組み合わせることで、位相2πにおける不連続性を改善することが可能となり、少数のパターンで高精度の三次元計測を実現できる。 The three-dimensional measuring device includes a plurality of light source units and a plurality of imaging units, and a predetermined pattern of measurement light is a sinusoidal stripe pattern, and the measurement unit is a phase shift using a sinusoidal stripe pattern. The three-dimensional shape of the object to be measured may be measured based on the method and the active stereo method. In this case, by converting the measured phase into height, the height of the object to be measured can be measured at intervals smaller than the pitch of the sinusoidal stripe pattern. Further, in the phase shift method, the discontinuity at the phase 2π is a problem. On the other hand, by combining the active stereo method using a plurality of imaging units, it is possible to improve the discontinuity in the phase 2π, and it is possible to realize highly accurate three-dimensional measurement with a small number of patterns.
光源部及び撮像部は、立体物の表面に配置されていてもよい。この場合、光源部及び撮像部が配置された立体物を三次元計測装置のプローブとして構成できる。立体物を用いることで、光源部及び撮像部のそれぞれの組を互いに異なる方向に向けることができるため、被計測物の三次元形状計測を広い立体角で実施することが可能となる。また、例えば口腔検査、内視鏡検査、管の内部や壁の隙間といった狭小箇所の検査、家具や装置等の床下からの検査といった用途への適用や、ハンディタイプの三次元計測装置の構築が容易となる。 The light source unit and the image pickup unit may be arranged on the surface of a three-dimensional object. In this case, a three-dimensional object in which the light source unit and the imaging unit are arranged can be configured as a probe of the three-dimensional measuring device. By using a three-dimensional object, each pair of the light source unit and the image pickup unit can be oriented in different directions, so that the three-dimensional shape measurement of the object to be measured can be performed at a wide solid angle. In addition, for example, it can be applied to applications such as oral examination, endoscopy, inspection of narrow places such as the inside of pipes and gaps in walls, inspection from under the floor of furniture and equipment, and construction of handy type 3D measuring equipment. It will be easy.
本開示によれば、装置の小型化により適用範囲の拡大が図られ、且つ計測精度の向上が図られる。 According to the present disclosure, the application range can be expanded and the measurement accuracy can be improved by downsizing the apparatus.
以下、図面を参照しながら、本開示の一側面に係る三次元計測装置の好適な実施形態について詳細に説明する。 Hereinafter, preferred embodiments of the three-dimensional measuring device according to one aspect of the present disclosure will be described in detail with reference to the drawings.
本実施形態に係る三次元計測装置101は、被計測物SAに所定パターンを有する計測光105を照射する一又は複数の光源部102と、計測光105が照射された被計測物SAを撮像する一又は複数の撮像部103と、撮像部103による撮像結果に基づいて被計測物SAの三次元形状を計測する計測部104と、を備えて構成されている(図15等参照)。また、光源部102は、M点発振のS-iPMSEL(Static-integrable Phase Modulating Surface Emitting Lasers)1によって構成されている。
The three-dimensional measuring device 101 according to the present embodiment captures one or a plurality of
三次元計測装置101では、針先ほどのサイズで構成されるS-iPMSEL1を用いて光源部102を構成することにより、装置全体を小型化でき、装置の適用範囲の拡大が図られる。また、三次元計測装置101では、M点発振のS-iPMSEL1を用いることで、所望の2次元的なパターンの光像とは異なる0次光(位相変調されない回折波成分)の出力をなくすことができる。これにより、0次光によるノイズや歪みのないパターンの計測光105を被計測物SAに照射することが可能となり、計測精度の向上が図られる。
[M点発振のS-iPMSEL]
In the three-dimensional measuring device 101, by configuring the
[S-iPMSEL of M point oscillation]
まず、M点発振のS-iPMSEL1について説明する。図1は、S-iPMSELの構成を示す一部断面斜視図である。図2は、S-iPMSELの積層構造を示す断面図である。図1では、S-iPMSEL1の中心においてS-iPMSEL1の厚さ方向に延びる軸をZ軸とするXYZ直交座標系を定義している。 First, S-iPMSEL1 of M point oscillation will be described. FIG. 1 is a partial cross-sectional perspective view showing the configuration of S-iPMSEL. FIG. 2 is a cross-sectional view showing a laminated structure of S-iPMSEL. FIG. 1 defines an XYZ Cartesian coordinate system whose Z axis is an axis extending in the thickness direction of the S-iPMSEL1 at the center of the S-iPMSEL1.
S-iPMSEL1は、XY面内方向において定在波を形成し、位相制御された平面波をZ軸方向に出力するレーザ光源である。S-iPMSEL1は、半導体基板10の主面10aに垂直な方向(すなわちZ軸方向)又はこれに対して傾斜した方向、或いはその両方を含む二次元的な任意形状の光像を出力する。
The S-iPMSEL1 is a laser light source that forms a standing wave in the XY in-plane direction and outputs a phase-controlled plane wave in the Z-axis direction. The S-iPMSEL1 outputs a two-dimensional arbitrary-shaped optical image including a direction perpendicular to the
図1及び図2に示されるように、S-iPMSEL1は、半導体基板10上に設けられた発光部としての活性層12と、活性層12を挟む一対のクラッド層11,13と、クラッド層13上に設けられたコンタクト層14と、を備えている。これらの半導体基板10、クラッド層11,13、及びコンタクト層14は、例えばGaAs系半導体、InP系半導体、もしくは窒化物系半導体といった化合物半導体によって構成されている。クラッド層11のエネルギーバンドギャップ、及びクラッド層13のエネルギーバンドギャップは、活性層12のエネルギーバンドギャップよりも大きくなっている。半導体基板10及び各層11~14の厚さ方向は、Z軸方向と一致している。
As shown in FIGS. 1 and 2, the S-
S-iPMSEL1は、活性層12と光学的に結合された位相変調層15を更に備えている。本実施形態では、位相変調層15は、活性層12とクラッド層13との間に設けられている。位相変調層15の厚さ方向は、Z軸方向と一致している。位相変調層15は、クラッド層11と活性層12との間に設けられていてもよい。活性層12とクラッド層13との間、及び活性層12とクラッド層11との間のうち少なくとも一方には、必要に応じて光ガイド層が設けられてもよい。光ガイド層は、キャリアを活性層12に効率的に閉じ込めるためのキャリア障壁層を含んでいてもよい。
The S-
位相変調層15は、第1屈折率媒質からなる基本層15aと、第1屈折率媒質とは屈折率の異なる第2屈折率媒質からなり、基本層15a内に存在する複数の異屈折率領域15bとを含んで構成されている。複数の異屈折率領域15bは、略周期構造を含んでいる。モードの等価屈折率をnとした場合、位相変調層15によって選択される波長λ0(=(√2)a×n、aは格子間隔)は、活性層12の発光波長範囲内に含まれている。位相変調層15は、活性層12の発光波長のうちの波長λ0近傍のバンド端波長を選択して、外部に出力することができる。位相変調層15内に入射したレーザ光は、位相変調層15内において異屈折率領域15bの配置に応じた所定のモードを形成し、所望のパターンを有するレーザビームとして、S-iPMSEL1の表面から外部に出射される。
The
S-iPMSEL1は、コンタクト層14上に設けられた電極16と、半導体基板10の裏面10b上に設けられた電極17とを更に備えている。電極16は、コンタクト層14とオーミック接触しており、電極17は、半導体基板10とオーミック接触している。電極17は、開口17aを有している。電極16は、コンタクト層14の中央領域に設けられている。コンタクト層14上における電極16以外の部分は、保護膜18(図2を参照)によって覆われている。電極16と接触していないコンタクト層14は、電流範囲の限定のために除去されていてもよい。半導体基板10の裏面10bのうち、電極17以外の部分は、開口17a内を含めて反射防止膜19によって覆われている。開口17a以外の領域にある反射防止膜19は、除去されていてもよい。
The S-
S-iPMSEL1では、電極16と電極17との間に駆動電流が供給されると、活性層12内において電子と正孔の再結合が生じ、活性層12が発光する。この発光に寄与する電子、正孔、及び活性層12で発生した光は、クラッド層11及びクラッド層13の間に効率的に閉じ込められる。
In S-iPMSEL1, when a driving current is supplied between the
活性層12から出射された光は、位相変調層15の内部に入射し、位相変調層15の内部の格子構造に応じた所定のモードを形成する。位相変調層15から出射したレーザ光は、裏面10bから開口17aを通ってS-iPMSEL1の外部へ直接的に出力される。或いは、位相変調層15から出射したレーザ光は、電極16において反射したのち、裏面10bから開口17aを通ってS-iPMSEL1の外部へ出力される。このとき、レーザ光に含まれる信号光(計測光105)は、主面10aに垂直な方向又はこれに対して傾斜した方向を含む二次元的な任意方向へ出射する。所望の光像を形成するのは、この信号光である。信号光は、主としてレーザ光の1次光及び-1次光である。本実施形態の位相変調層15からは、レーザ光の0次光が出力されないようになっている。
The light emitted from the
図3は、位相変調層15の平面図である。同図に示すように、位相変調層15は、第1屈折率媒質からなる基本層15aと、第1屈折率媒質とは屈折率の異なる第2屈折率媒質からなる複数の異屈折率領域15bとを含んでいる。図3では、位相変調層15に対し、XY面内における仮想的な正方格子を設定している。正方格子の一辺は、X軸と平行であり、他辺はY軸と平行である。正方格子の格子点Oを中心とする正方形状の単位構成領域Rは、X軸に沿った複数列及びY軸に沿った複数行にわたって二次元状に設定されている。各単位構成領域RのXY座標をそれぞれの単位構成領域Rの重心位置で規定すると、これらの重心位置は、仮想的な正方格子の格子点Oに一致する。複数の異屈折率領域15bは、各単位構成領域R内に例えば1つずつ設けられる。異屈折率領域15bの平面形状は、例えば円形状である。格子点Oは、異屈折率領域15bの外部に位置してもよく、異屈折率領域15bの内部に含まれていてもよい。
FIG. 3 is a plan view of the
1つの単位構成領域R内に占める異屈折率領域15bの面積Sの比率は、フィリングファクタ(FF)と称される。正方格子の格子間隔をaとすると、異屈折率領域15bのフィリングファクタFFは、S/a2として与えられる。Sは、XY平面における異屈折率領域15bの面積である。例えば異屈折率領域15bの形状が真円形状の場合には、フィリングファクタFFは、真円の直径dを用いてS=π(d/2)2として与えられる。異屈折率領域15bの形状が正方形の場合には、フィリングファクタFFは、正方形の一辺の長さLAを用いてS=LA2として与えられる。
The ratio of the area S of the different
図4は、単位構成領域Rを拡大して示す図である。同図に示すように、異屈折率領域15bのそれぞれは、重心Gを有する。ここでは、格子点Oから重心Gに向かうベクトルとX軸とのなす角度をφ(x,y)とする。xは、X軸におけるx番目の格子点の位置、yは、Y軸におけるy番目の格子点の位置を示す。回転角度φが0°である場合、格子点Oと重心Gとを結ぶベクトルの向きは、X軸の正方向と一致する。また、格子点Oと重心Gとを結ぶベクトルの長さをr(x,y)とする。一例では、r(x,y)は、x、yによらず、位相変調層15の全体にわたって一定である。
FIG. 4 is an enlarged view showing the unit constituent area R. As shown in the figure, each of the different
図3に示したように、格子点Oと重心Gとを結ぶベクトルの向き、すなわち、異屈折率領域15bの重心Gの格子点O周りの回転角度φは、所望の光像に応じた位相パターンに従って、各格子点O毎に個別に設定される。位相パターン、すなわち、回転角度分布φ(x,y)は、x,yの値で決まる位置毎に特定の値を有するが、必ずしも特定の関数で表わされるとは限らない。回転角度分布φ(x,y)は、所望の光像をフーリエ変換して得られる複素振幅分布のうち、位相分布を抽出したものから決定される。所望の光像から複素振幅分布を求める際には、ホログラム生成の計算時に一般的に用いられるGerchberg-Saxton(GS)法のような繰り返しアルゴリズムを適用することによって、ビームパターンの再現性を向上させることが可能である。
As shown in FIG. 3, the direction of the vector connecting the grid point O and the center of gravity G, that is, the rotation angle φ around the grid point O of the center of gravity G in the different
図5は、位相変調層の特定領域内に屈折率略周期構造を適用した例を示す平面図である。図5に示す例では、正方形の内側領域RINの内部に、目的となるビームパターンを出射するための略周期構造(例えば図3のに示した構造)が形成されている。一方、内側領域RINを囲む外側領域ROUTには、正方格子の格子点位置に、重心位置が一致する真円形の異屈折率領域が配置されている。外側領域ROUTにおけるフィリングファクタFFは、例えば12%に設定されている。内側領域RINの内部及び外側領域ROUT内では、仮想的に設定される正方格子の格子間隔は、同一(=a)である。この構造の場合、外側領域ROUT内にも光が分布するため、内側領域RINの周辺部での光強度の急激な変化によって生じる高周波ノイズ(いわゆる窓関数ノイズ)の発生を抑制できる。また、面内方向への光漏れを抑制することができ、閾値電流の低減が期待できる。 FIG. 5 is a plan view showing an example in which a refractive index substantially periodic structure is applied within a specific region of the phase modulation layer. In the example shown in FIG. 5, a substantially periodic structure (for example, the structure shown in FIG. 3) for emitting a target beam pattern is formed inside the inner region RIN of the square. On the other hand, in the outer region ROUT surrounding the inner region RIN, a perfect circular different refractive index region having the same center of gravity position is arranged at the grid point position of the square lattice. The filling factor FF in the outer region ROUT is set to, for example, 12%. Within the inner region RIN and the outer region ROUT, the grid spacing of the square grid virtually set is the same (= a). In the case of this structure, since the light is also distributed in the outer region ROUT, it is possible to suppress the generation of high frequency noise (so-called window function noise) caused by a sudden change in the light intensity in the peripheral portion of the inner region RIN. In addition, light leakage in the in-plane direction can be suppressed, and a reduction in the threshold current can be expected.
図6は、S-iPMSEL1の出力ビームパターンが結像して得られる光像と、位相変調層15における回転角度分布φ(x,y)との関係を説明する図である。出力ビームパターンの中心Qは、半導体基板10の主面10aに対して垂直な軸線上に位置するとは限らないが、垂直な軸線上に配置させることもできる。図6では、説明の便宜のため、中心Qが主面10aに対して垂直な軸線上にあるものとする。図6には、中心Qを原点とする4つの象限が示されている。図6の例では、第3象限に文字「A」が現れ、第1象限に文字「A」を180度回転したものが現れている。出力ビームパターンが回転対称な光像(例えば、十字、丸、二重丸など)である場合には、重なって一つの光像として観察される。なお、図4に示したように、S-iPMSEL1における異屈折率領域15bの重心Gを格子点O回りの円周方向にずらす場合には、図6に示すように、第1象限の出力ビームパターンと第3象限の出力ビームパターンとの間に強度差はないが、後述の図14のように、S-iPMSEL1における異屈折率領域15bの重心Gを格子点Oを通る直線上にずらす場合には、第1象限の出力ビームパターンと第3象限の出力ビームパターンとの間に強度差を持たせることが可能である。
FIG. 6 is a diagram illustrating the relationship between the optical image obtained by forming an image of the output beam pattern of S-iPMSEL1 and the rotation angle distribution φ (x, y) in the
S-iPMSEL1の出力ビームパターンの光像は、スポット、ドット、直線、十字架、線画、格子パターン、写真、縞状パターン、CG(コンピュータグラフィクス)、及び文字のうち少なくとも1つを含んでいる。所望の光像を得るためには、以下の手順によって、位相変調層15における異屈折率領域15bの回転角度分布φ(x、y)を決定する。
The optical image of the output beam pattern of S-iPMSEL1 includes at least one of spots, dots, straight lines, crosses, line art, grid patterns, photographs, striped patterns, CG (computer graphics), and characters. In order to obtain a desired optical image, the rotation angle distribution φ (x, y) of the different
第1の前提条件として、法線方向に一致するZ軸と、複数の異屈折率領域15bを含む位相変調層15の一方の面に一致したX-Y平面と、によって規定されるXYZ直交座標系において、正方形状を有するM1(1以上の整数)×N1(1以上の整数)個の単位構成領域Rにより構成される仮想的な正方格子をX-Y平面上に設定する。
As a first precondition, the XYZ Cartesian coordinates defined by the Z axis matching the normal direction and the XY plane matching one surface of the
第2の前提条件として、XYZ直交座標系における座標(ξ,η,ζ)は、図7に示すように、動径の長さrと、Z軸からの傾き角θtiltと、X-Y平面上で特定されるX軸からの回転角θrotと、で規定される球面座標(r,θrot,θtilt)に対して、以下の式(1)~式(3)で示された関係を満たしているものとする。図7は、球面座標(r,θrot,θtilt)からXYZ直交座標系における座標(ξ,η,ζ)への座標変換を説明するための図であり、座標(ξ,η,ζ)により、実空間であるXYZ直交座標系において設定される所定平面上の設計上の光像が表現される。 As a second precondition, the coordinates (ξ, η, ζ) in the XYZ Cartesian coordinate system are, as shown in FIG. 7, the length r of the driving diameter, the inclination angle θtilt from the Z axis, and the XY plane. The relationship shown by the following equations (1) to (3) is satisfied with respect to the rotation angle θrot from the X axis specified above and the spherical coordinates (r, θrot, θtilt) defined by. It is assumed that there is. FIG. 7 is a diagram for explaining coordinate conversion from spherical coordinates (r, θrot, θtilt) to coordinates (ξ, η, ζ) in the XYZ Cartesian coordinate system, and is a diagram for explaining the coordinate conversion (ξ, η, ζ). A design optical image on a predetermined plane set in the XYZ Cartesian coordinate system, which is a real space, is represented.
S-iPMSEL1から出力される光像に相当するビームパターンを角度θtilt及びθrotで規定される方向に向かう輝点の集合とするとき、角度θtiltおよびθrotは、以下の式(4)で規定される規格化波数であってX軸に対応したKx軸上の座標値kxと、以下の式(5)で規定される規格化波数であってY軸に対応すると共にKx軸に直交するKy軸上の座標値kyに換算されるものとする。規格化波数は、仮想的な正方格子の格子間隔に相当する波数2π/aを1.0として規格化された波数を意味する。このとき、Kx軸およびKy軸により規定される波数空間において、光像に相当するビームパターンを含む特定の波数範囲は、それぞれが正方形状のM2(1以上の整数)×N2(1以上の整数)個の画像領域FRで構成される。なお、整数M2は、整数M1と一致する必要はない。同様に、整数N2は、整数N1と一致する必要もない。また、式(4)および式(5)は、例えばY. Kurosaka et al.," Effects of non-lasing band intwo-dimensionalphotonic-crystal lasers clarified using omnidirectional bandstructure,"Opt. Express 20, 21773-21783 (2012)に開示されている。
a:仮想的な正方格子の格子定数
λ:S-iPMSEL1の発振波長
When the beam pattern corresponding to the optical image output from S-iPMSEL1 is a set of bright spots directed in the directions defined by the angles θtilt and θrot, the angles θtilt and θrot are defined by the following equation (4). The coordinate value kx on the Kx axis, which is the normalized wave number and corresponds to the X axis, and the Ky axis, which is the normalized wave number defined by the following equation (5) and corresponds to the Y axis and is orthogonal to the Kx axis. It shall be converted into the coordinate value ky of. The normalized wave number means a wave number standardized with a wave number of 2π / a corresponding to the grid spacing of a virtual square grid as 1.0. At this time, in the wavenumber space defined by the Kx axis and the Ky axis, the specific wavenumber range including the beam pattern corresponding to the optical image is each square-shaped M2 (integer of 1 or more) × N2 (integer of 1 or more). ) Consists of image area FR. The integer M2 does not have to match the integer M1. Similarly, the integer N2 does not have to match the integer N1. Equations (4) and (5) also include, for example, Y. Kurosaka et al., "Effects of non-lasing band intwo-dimensionalphotonic-crystal lasers clarified using omnidirectional bandstructure," Opt. Express 20, 21773-21783 (2012). ).
a: Lattice constant of a virtual square lattice
λ: Oscillation wavelength of S-iPMSEL1
第3の前提条件として、波数空間において、Kx軸方向の座標成分kx(0以上M2-1以下の整数)とKy軸方向の座標成分ky(0以上N2-1以下の整数)とで特定される画像領域FR(kx,ky)それぞれを、X軸方向の座標成分x(0以上M1-1以下の整数)とY軸方向の座標成分y(0以上N1-1以下の整数)とで特定されるX-Y平面上の単位構成領域R(x,y)に二次元逆離散フーリエ変換することで得られる複素振幅F(x,y)は、jを虚数単位として、以下の式(6)で与えられる。複素振幅F(x,y)は、振幅項をA(x,y)とすると共に位相項をP(x,y)とするとき、以下の式(7)により規定される。第4の前提条件として、単位構成領域R(x,y)は、X軸およびY軸にそれぞれ平行であって単位構成領域R(x,y)の中心となる格子点O(x,y)において直交するs軸およびt軸で規定される。
上記第1~第4の前提条件の下、位相変調層15は、以下の第5条件及び第6条件を満たすように構成される。すなわち、第5条件は、単位構成領域R(x,y)内において、重心Gが格子点O(x,y)から離れた状態で配置されていることで満たされる。第6条件は、格子点O(x,y)から対応する重心Gまでの線分長r2(x,y)がM1個×N1個の単位構成領域Rそれぞれにおいて共通の値に設定された状態で、格子点O(x,y)と対応する重心Gとを結ぶ線分と、s軸と、の成す角度φ(x,y)が、
φ(x,y)=C×P(x,y)+B
C:比例定数であって例えば180°/π
B:任意の定数であって例えば0
となる関係を満たすように、対応する異屈折率領域15bが単位構成領域R(x,y)内に配置されることで満たされる。
Under the first to fourth preconditions, the
φ (x, y) = C × P (x, y) + B
C: Proportional constant, for example 180 ° / π
B: An arbitrary constant, for example 0
The corresponding different
次に、S-iPMSEL1のM点発振について説明する。S-iPMSEL1のM点発振のためには、仮想的な正方格子の格子間隔a、活性層12の発光波長λ、及びモードの等価屈折率nが、λ=(√2)n×aといった条件を満たすとよい。図8は、M点発振のS-iPMSELの位相変調層に関する逆格子空間を示す平面図である。図中の点Pは、逆格子点を表している。図中の矢印B1は、基本逆格子ベクトルを表しており、矢印K1,K2,K3,及びK4は、4つの面内波数ベクトルを表している。面内波数ベクトルK1~K4は、回転角度分布φ(x,y)による波数拡がりSPをそれぞれ有している。
Next, the M point oscillation of S-iPMSEL1 will be described. For the M-point oscillation of S-iPMSEL1, the conditions such that the lattice spacing a of the virtual square lattice, the emission wavelength λ of the
なお、波数拡がりSPの形状及び大きさは、上述したΓ点発振の場合と同様である。M点発振のS-iPMSEL1では、面内波数ベクトルK1~K4の大きさ(すなわち面内方向の定在波の大きさ)は、基本逆格子ベクトルB1の大きさよりも小さくなっている。したがって、面内波数ベクトルK1~K4と基本逆格子ベクトルB1とのベクトル和が0にはならず、回折によって面内方向の波数が0となり得ないため、面垂直方向(Z軸方向)への回折は生じない。このままでは、M点発振のS-iPMSEL1では、面垂直方向(Z軸方向)への0次光、Z軸方向に対して傾斜した方向への1次光及び-1次光が出力しない。 The shape and size of the wave number spreading SP are the same as in the case of the above-mentioned Γ point oscillation. In the M-point oscillation S-iPMSEL1, the magnitude of the in-plane wave vector K1 to K4 (that is, the magnitude of the standing wave in the in-plane direction) is smaller than the magnitude of the basic reciprocal lattice vector B1. Therefore, the vector sum of the in-plane wave vector K1 to K4 and the basic reciprocal lattice vector B1 cannot be 0, and the wave number in the in-plane direction cannot be 0 due to diffraction. No diffraction occurs. As it is, the S-iPMSEL1 of M point oscillation does not output the 0th-order light in the plane vertical direction (Z-axis direction), the 1st-order light in the direction inclined with respect to the Z-axis direction, and the -1st-order light.
本実施形態では、M点発振のS-iPMSEL1において次のような工夫を位相変調層15に施すことにより、0次光を出力させずに、1次光及び-1次光の一部を出力させることができる。具体的には、図9に示すように、面内波数ベクトルK1~K4に対し、ある一定の大きさ及び向きを有する回折ベクトルVを加えることにより、面内波数ベクトルK1~K4のうち少なくとも1つ(図では面内波数ベクトルK3)の大きさを2π/λよりも小さくする。言い換えると、回折ベクトルVが加えられた後の面内波数ベクトルK1~K4のうち少なくとも1つ(面内波数ベクトルK3)を、半径2π/λの円状領域(ライトライン)LL内に収める。
In the present embodiment, by applying the following device to the
図9において破線で示される面内波数ベクトルK1~K4は、回折ベクトルVの加算前を表しており、実線で示される面内波数ベクトルK1~K4は、回折ベクトルVの加算後を表している。ライトラインLLは、全反射条件に対応しており、ライトラインLL内に収まる大きさの波数ベクトルは、面垂直方向(Z軸方向)の成分を有することとなる。一例では、回折ベクトルVの方向は、Γ-M1軸又はΓ-M2軸に沿っている。回折ベクトルVの大きさは、2π/(√2)a-2π/λから2π/(√2)a+2π/λの範囲内となっており、一例として、2π/(√2)aとなっている。 In-plane wave vector K1 to K4 shown by a broken line in FIG. 9 represent before addition of the diffraction vector V, and in-plane wave vector K1 to K4 shown by a solid line represent after addition of the diffraction vector V. .. The light line LL corresponds to the total reflection condition, and the wave vector having a size within the light line LL has a component in the plane vertical direction (Z-axis direction). In one example, the direction of the diffraction vector V is along the Γ-M1 axis or the Γ-M2 axis. The magnitude of the diffraction vector V is in the range of 2π / (√2) a-2π / λ to 2π / (√2) a + 2π / λ, and is 2π / (√2) a as an example. There is.
続いて、面内波数ベクトルK1~K4のうち、少なくとも1つをライトラインLL内に収めるための回折ベクトルVの大きさ及び向きについて検討する。下記の数式(8)~(11)は、回折ベクトルVが加えられる前の面内波数ベクトルK1~K4を示す。
波数ベクトルの広がりΔkx及びΔkyは、下記の数式(12)及び(13)をそれぞれ満たす。面内波数ベクトルのx軸方向の広がりの最大値Δkxmax及びy軸方向の広がりの最大値Δkymaxは、設計の光像の角度広がりにより規定される。
The wave vector spreads Δkx and Δky satisfy the following equations (12) and (13), respectively. The maximum value Δkxmax of the spread in the x-axis direction and the maximum value Δkymax of the spread in the y-axis direction of the in-plane wave vector are defined by the angular spread of the optical image of the design.
回折ベクトルVを下記の数式(14)のように表したとき、回折ベクトルVが加えられた後の面内波数ベクトルK1~K4は下記の数式(15)~(18)となる。
数式(15)~(18)において波数ベクトルK1~K4のいずれかがライトラインLL内に収まることを考慮すると、下記の数式(19)の関係が成り立つ。
すなわち、数式(19)を満たす回折ベクトルVを加えることにより、波数ベクトルK1~K4のいずれかがライトラインLL内に収まり、1次光及び-1次光の一部が出力される。
Considering that any of the wave vector K1 to K4 fits within the light line LL in the formulas (15) to (18), the following relation of the formula (19) is established.
That is, by adding the diffraction vector V satisfying the equation (19), any of the wave vector K1 to K4 is contained in the light line LL, and a part of the primary light and the -1st order light is output.
ライトラインLLの大きさ(半径)を2π/λとしたのは、以下の理由による。図10は、ライトラインLLの周辺構造を模式的に説明するための図である。同図では、Z軸方向に垂直な方向から見たデバイスと空気との境界を示している。真空中の光の波数ベクトルの大きさは2π/λとなるが、図10のようにデバイス媒質中を光が伝搬するときには、屈折率nの媒質内の波数ベクトルKaの大きさは2πn/λとなる。このとき、デバイスと空気の境界を光が伝搬するためには、境界に平行な波数成分が連続している必要がある(波数保存則)。 The size (radius) of the light line LL is set to 2π / λ for the following reasons. FIG. 10 is a diagram for schematically explaining the peripheral structure of the light line LL. The figure shows the boundary between the device and air as seen from the direction perpendicular to the Z-axis direction. The magnitude of the wave vector of light in vacuum is 2π / λ, but when light propagates in the device medium as shown in FIG. 10, the magnitude of the wave vector Ka in the medium having a refractive index n is 2πn / λ. Will be. At this time, in order for light to propagate along the boundary between the device and air, the wavenumber component parallel to the boundary must be continuous (wavenumber conservation law).
図10において、波数ベクトルKaとZ軸とが角度θをなす場合、面内に投影した波数ベクトル(すなわち面内波数ベクトル)Kbの長さは、(2πn/λ)sinθとなる。一方で、一般には媒質の屈折率n>1の関係から、媒質内の面内波数ベクトルKbが2π/λより大きくなる角度では、波数保存則が成立しなくなる。このとき、光は全反射し、空気側に取り出すことができなくなる。この全反射条件に対応する波数ベクトルの大きさがライトラインLLの大きさ、すなわち、2π/λとなる。 In FIG. 10, when the wave vector Ka and the Z axis form an angle θ, the length of the wave vector (that is, the in-plane wave vector) Kb projected in the plane is (2πn / λ) sin θ. On the other hand, in general, due to the relationship of the refractive index n> 1 of the medium, the wavenumber conservation law does not hold at an angle where the in-plane wave vector Kb in the medium becomes larger than 2π / λ. At this time, the light is totally reflected and cannot be taken out to the air side. The magnitude of the wave vector corresponding to this total reflection condition is the magnitude of the light line LL, that is, 2π / λ.
面内波数ベクトルK1~K4に回折ベクトルVを加える具体的な方式の一例として、光像に応じた位相分布である回転角度分布φ1(x,y)(第1の位相分布)に対し、光像とは無関係の回転角度分布φ2(x,y)(第2の位相分布)を重畳する方式が考えられる。この場合、位相変調層15の回転角度分布φ(x,y)は、φ(x,y)=φ1(x,y)+φ2(x,y)として表される。φ1(x,y)は、前に述べたように光像をフーリエ変換したときの複素振幅の位相に相当する。また、φ2(x,y)は、上記の数式(19)を満たす回折ベクトルVを加えるための回転角度分布である。
As an example of a specific method of adding the diffraction vector V to the in-plane wave vector K1 to K4, light is applied to the rotation angle distribution φ1 (x, y) (first phase distribution) which is the phase distribution according to the optical image. A method of superimposing a rotation angle distribution φ2 (x, y) (second phase distribution) unrelated to the image can be considered. In this case, the rotation angle distribution φ (x, y) of the
図11は、回転角度分布φ2(x,y)の一例を概念的に示す図である。同図の例では、第1の位相値φAと、第1の位相値φAとは異なる値の第2の位相値φBとが市松模様に配列されている。一例では、位相値φAは、0(rad)であり、位相値φBは、π(rad)である。この場合、第1の位相値φAと、第2の位相値φBとがπずつ変化する。このような位相値の配列によって、Γ-M1軸又はΓ-M2軸に沿う回折ベクトルVを好適に実現することができる。市松模様の配列の場合、V=(±π/a,±π/a)となり、回折ベクトルVと図8の波数ベクトルK1~K4とが丁度相殺される。なお、回折ベクトルVの角度分布θ2(x,y)は、回折ベクトルV(Vx,Vy)と位置ベクトルr(x,y)との内積で表される。すなわち、回折ベクトルVの角度分布θ2(x,y)は、θ2(x,y)=V・r=Vxx+Vyyで表される。 FIG. 11 is a diagram conceptually showing an example of the rotation angle distribution φ2 (x, y). In the example of the figure, the first phase value φA and the second phase value φB having a value different from the first phase value φA are arranged in a checkered pattern. In one example, the phase value φA is 0 (rad) and the phase value φB is π (rad). In this case, the first phase value φA and the second phase value φB change by π. With such an array of phase values, a diffraction vector V along the Γ-M1 axis or the Γ-M2 axis can be suitably realized. In the case of the checkered pattern arrangement, V = (± π / a, ± π / a), and the diffraction vector V and the wave vector K1 to K4 in FIG. 8 are just canceled out. The angular distribution θ2 (x, y) of the diffraction vector V is represented by the inner product of the diffraction vector V (Vx, Vy) and the position vector r (x, y). That is, the angular distribution θ2 (x, y) of the diffraction vector V is represented by θ2 (x, y) = V · r = Vxx + Vyy.
上記実施形態において、光像の角度広がりに基づく波数広がりが、波数空間上の或る点を中心とする半径Δkの円に含まれる場合、次のように簡略に考えることもできる。4方向の面内波数ベクトルK1~K4に回折ベクトルVを加えることにより、4方向の面内波数ベクトルK1~K4のうち少なくとも1つの大きさを2π/λ(ライトラインLL)よりも小さくする。このことは、4方向の面内波数ベクトルK1~K4から波数拡がりΔkを除いたものに対して回折ベクトルVを加えることにより、4方向の面内波数ベクトルK1~K4のうち少なくとも1つの大きさを、2π/λから波数拡がりΔkを差し引いた値{(2π/λ)-Δk}より小さくする、と考えてよい。 In the above embodiment, when the wavenumber spread based on the angular spread of the light image is included in a circle having a radius Δk centered on a certain point in the wavenumber space, it can be simply considered as follows. By adding the diffraction vector V to the in-plane wave vector K1 to K4 in the four directions, the magnitude of at least one of the in-plane wave vectors K1 to K4 in the four directions is made smaller than 2π / λ (light line LL). This means that by adding the diffraction vector V to the wave vector K1 to K4 in the four directions minus the wave number spread Δk, the magnitude of at least one of the wave vectors K1 to K4 in the four directions is large. Can be considered to be smaller than the value {(2π / λ) −Δk} obtained by subtracting the wavenumber spread Δk from 2π / λ.
図12は、上記の状態を概念的に示す図である。同図に示すように、波数拡がりΔkを除いた面内波数ベクトルK1~K4に対して回折ベクトルVを加えると、面内波数ベクトルK1~K4のうち少なくとも1つの大きさが{(2π/λ)-Δk}よりも小さくなる。図12において、領域LL2は、半径が{(2π/λ)-Δk}の円状の領域である。図12において、破線で示される面内波数ベクトルK1~K4は、回折ベクトルVの加算前を表しており、実線で示される面内波数ベクトルK1~K4は、回折ベクトルVの加算後を表している。領域LL2は、波数拡がりΔkを考慮した全反射条件に対応しており、領域LL2内に収まる大きさの波数ベクトルは、面垂直方向(Z軸方向)にも伝搬することとなる。 FIG. 12 is a diagram conceptually showing the above state. As shown in the figure, when the diffraction vector V is added to the in-plane wave vector K1 to K4 excluding the wave number spread Δk, the magnitude of at least one of the in-plane wave vectors K1 to K4 becomes {(2π / λ). ) -Δk}. In FIG. 12, the region LL2 is a circular region having a radius of {(2π / λ) −Δk}. In FIG. 12, the in-plane wave vector K1 to K4 shown by the broken line represents before the addition of the diffraction vector V, and the in-plane wave vector K1 to K4 shown by the solid line represent after the addition of the diffraction vector V. There is. The region LL2 corresponds to the total reflection condition in consideration of the wavenumber spread Δk, and the wave vector having a size within the region LL2 also propagates in the plane vertical direction (Z-axis direction).
本形態において、面内波数ベクトルK1~K4のうち少なくとも1つを領域LL2内に収めるための回折ベクトルVの大きさ及び向きを説明する。下記の数式(20)~(23)は、回折ベクトルVが加えられる前の面内波数ベクトルK1~K4を示す。
ここで、回折ベクトルVを前述した数式(14)のように表したとき、回折ベクトルVが加えられた後の面内波数ベクトルK1~K4は、下記の数式(24)~(27)となる。
数式(24)~(27)において、面内波数ベクトルK1~K4のいずれかが領域LL2内に収まることを考慮すると、下記の数式(28)の関係が成り立つ。すなわち、数式(28)を満たす回折ベクトルVを加えることにより、波数拡がりΔkを除いた面内波数ベクトルK1~K4のいずれかが領域LL2内に収まる。このような場合であっても、0次光を出力させずに、1次光及び-1次光の一部を出力させることができる。
図13は、変形例に係る位相変調層の平面図である。また、図14は、変形例に係る位相変調層における異屈折率領域の位置関係を示す図である。図13及び図14に示すように、変形例に係る位相変調層15の各異屈折率領域15bの重心Gは、直線D上に配置されている。直線Dは、各単位構成領域Rに対応する格子点Oを通り、正方格子の各辺に対して傾斜する直線である。つまり、直線Dは、X軸及びY軸の双方に対して傾斜する直線である。正方格子の一辺(X軸)に対する直線Dの傾斜角は、θである。
FIG. 13 is a plan view of the phase modulation layer according to the modified example. Further, FIG. 14 is a diagram showing the positional relationship of the different refractive index regions in the phase modulation layer according to the modified example. As shown in FIGS. 13 and 14, the center of gravity G of each different
傾斜角θは、位相変調層15B内において一定である。傾斜角θは、0°<θ<90°を満たし、一例ではθ=45°である。或いは、傾斜角θは、180°<θ<270°を満たし、一例ではθ=225°である。傾斜角θが0°<θ<90°または180°<θ<270°を満たす場合、直線Dは、X軸及びY軸によって規定される座標平面の第1象限から第3象限にわたって延びる。傾斜角θは、90°<θ<180°を満たし、一例ではθ=135°である。或いは、傾斜角θは、270°<θ<360°を満たし、一例ではθ=315°である。傾斜角θが90°<θ<180°または270°<θ<360°を満たす場合、直線Dは、X軸及びY軸によって規定される座標平面の第2象限から第4象限にわたって延びる。このように、傾斜角θは、0°、90°、180°及び270°を除く角度となっている。
The tilt angle θ is constant in the
ここで、格子点Oと重心Gとの距離をr(x,y)とする。xは、X軸におけるx番目の格子点の位置であり、yは、Y軸におけるy番目の格子点の位置である。距離r(x,y)が正の値である場合、重心Gは、第1象限(または第2象限)に位置する。距離r(x,y)が負の値である場合、重心Gは、第3象限(または第4象限)に位置する。距離r(x,y)が0である場合、格子点Oと重心Gとが互いに一致する。傾斜角度は、45°、135°、225°、275°が好適である。これらの傾斜角度では、M点の定在波を形成する4つの波数ベクトル(例えば、面内波数ベクトル(±π/a、±π/a))の中の2つのみが位相変調され、その他の2つが位相変調されないため、安定した定在波を形成することができる。 Here, the distance between the grid point O and the center of gravity G is r (x, y). x is the position of the x-th grid point on the X-axis, and y is the position of the y-th grid point on the Y-axis. When the distance r (x, y) is a positive value, the center of gravity G is located in the first quadrant (or the second quadrant). When the distance r (x, y) is a negative value, the center of gravity G is located in the third quadrant (or the fourth quadrant). When the distance r (x, y) is 0, the grid point O and the center of gravity G coincide with each other. The inclination angle is preferably 45 °, 135 °, 225 °, 275 °. At these tilt angles, only two of the four wave vectors (eg, in-plane wave vectors (± π / a, ± π / a)) that form the standing wave at point M are phase-modulated, and the others. Since the two are not phase-modulated, a stable standing wave can be formed.
各異屈折率領域の重心Gと各単位構成領域Rに対応する格子点Oとの距離r(x,y)は、所望の光像に応じた位相パターンに従って各異屈折率領域15b毎に個別に設定される。位相パターン、すなわち距離r(x,y)の分布は、x,yの値で決まる位置毎に特定の値を有するが、必ずしも特定の関数で表わされるとは限らない。距離r(x,y)の分布は、所望の光像を逆フーリエ変換して得られる複素振幅分布のうち位相分布を抽出したものから決定される。
The distance r (x, y) between the center of gravity G of each different refractive index region and the lattice point O corresponding to each unit constituent region R is individually set for each different
すなわち、図14に示すように、或る座標(x,y)における位相P(x,y)がP0である場合には、距離r(x,y)を0と設定し、位相P(x,y)がπ+P0である場合には、距離r(x,y)を最大値R0に設定し、位相P(x,y)が-π+P0である場合には、距離r(x,y)を最小値-R0に設定する。そして、その中間の位相P(x,y)に対しては、r(x,y)={P(x,y)-P0}×R0/πとなるように距離r(x,y)をとる。初期位相P0は、任意に設定することができる。 That is, as shown in FIG. 14, when the phase P (x, y) at a certain coordinate (x, y) is P 0 , the distance r (x, y) is set to 0, and the phase P ( When x, y) is π + P 0 , the distance r (x, y) is set to the maximum value R 0 , and when the phase P (x, y) is −π + P 0 , the distance r (x, y) is set. , Y) is set to the minimum value -R 0 . Then, for the intermediate phase P (x, y), the distance r (x, y) is such that r (x, y) = {P (x, y) −P 0 } × R 0 / π. ). The initial phase P 0 can be set arbitrarily.
仮想的な正方格子の格子間隔をaとすると、r(x,y)の最大値R0は、例えば下記式(29)の範囲内となる。所望の光像から複素振幅分布を求める際には、ホログラム生成の計算時に一般的に用いられるGerchberg-Saxton(GS)法のような繰り返しアルゴリズムを適用することによって、ビームパターンの再現性を向上させることが可能である。
本形態においては、位相変調層15の異屈折率領域15bの距離r(x,y)の分布を決定することにより、所望の光像を得ることができる。前述の実施形態と同様の第1~第4の前提条件の下、位相変調層15は、以下の条件を満たすよう構成される。すなわち、格子点O(x,y)から対応する異屈折率領域15bの重心Gまでの距離r(x,y)が、
r(x,y)=C×(P(x,y)-P0)
C:比例定数で例えばR0/π
P0:任意の定数であって例えば0
となる関係を満たすように、対応する異屈折率領域15bが単位構成領域R(x,y)内に配置される。
In this embodiment, a desired optical image can be obtained by determining the distribution of the distance r (x, y) in the different
r (x, y) = C × (P (x, y) -P 0 )
C: Proportional constant, for example R 0 / π
P 0 : An arbitrary constant, for example, 0
The corresponding different
すなわち、距離r(x,y)は、或る座標(x,y)における位相P(x,y)がP0である場合には0に設定され、位相P(x,y)がπ+P0である場合には最大値R0に設定され、位相P(x,y)が-π+P0である場合には最小値-R0に設定される。所望の光像を得たい場合、当該光像を逆フーリエ変換して、その複素振幅の位相P(x,y)に応じた距離r(x,y)の分布を複数の異屈折率領域15bに与えるとよい。位相P(x,y)と距離r(x,y)とは、互いに比例してもよい。
That is, the distance r (x, y) is set to 0 when the phase P (x, y) at a certain coordinate (x, y) is P 0 , and the phase P (x, y) is π + P 0 . If is, the maximum value is set to R 0 , and if the phase P (x, y) is −π + P 0 , the minimum value is set to −R 0 . When a desired optical image is desired, the optical image is subjected to inverse Fourier transform, and the distribution of the distance r (x, y) corresponding to the phase P (x, y) of the complex amplitude is distributed in a plurality of different
本形態においても、前述した実施形態と同様に、仮想的な正方格子の格子間隔aと活性層12の発光波長λとがM点発振の条件を満たす。さらに、位相変調層15において逆格子空間を考えるとき、距離r(x,y)の分布による波数拡がりをそれぞれ含む4方向の面内波数ベクトルのうち少なくとも1つの大きさは、2π/λ(ライトライン)よりも小さくすることができる。
Also in this embodiment, similarly to the above-described embodiment, the lattice spacing a of the virtual square lattice and the emission wavelength λ of the
本形態においては、M点で発振するS-iPMSEL1において次のような工夫を位相変調層15に施すことにより、0次光をライトライン内に出力させずに、1次光及び-1次光の一部を出力する。具体的には、図9に示したように、面内波数ベクトルK1~K4に対してある一定の大きさ及び向きを有する回折ベクトルVを加えることにより、面内波数ベクトルK1~K4のうち少なくとも1つの大きさを、2π/λよりも小さくする。すなわち、回折ベクトルVが加えられた後の面内波数ベクトルK1~K4のうち少なくとも1つを半径2π/λの円状領域(ライトライン)LL内に収める。前述した数式(19)を満たす回折ベクトルVを加えることにより、面内波数ベクトルK1~K4のいずれかがライトラインLL内に収まり、1次光及び-1次光の一部が出力される。
In this embodiment, the
或いは、図12に示したように、4方向の面内波数ベクトルK1~K4から波数拡がりΔkを除いたもの(すなわちM点発振の正方格子PCSELにおける4方向の面内波数ベクトル)に対して回折ベクトルVを加えることにより、4方向の面内波数ベクトルK1~K4のうち、少なくとも1つの大きさを2π/λから波数拡がりΔkを差し引いた値{(2π/λ)-Δk}より小さくしてもよい。すなわち、前述した数式(28)を満たす回折ベクトルVを加えることにより、面内波数ベクトルK1~K4のいずれかが領域LL2内に収まり、1次光及び-1次光の一部が出力される。 Alternatively, as shown in FIG. 12, the wave vector K1 to K4 in four directions minus the wave number spread Δk (that is, the in-plane wave vector in four directions in the square grid PCSEL of M point oscillation) is diffracted. By adding the vector V, at least one of the wave vector K1 to K4 in the four directions is made smaller than the value {(2π / λ) −Δk} obtained by subtracting the wavenumber expansion Δk from 2π / λ. May be good. That is, by adding the diffraction vector V satisfying the above-mentioned mathematical formula (28), any one of the in-plane wave vector K1 to K4 is contained in the region LL2, and a part of the primary light and the -1st order light is output. ..
面内波数ベクトルK1~K4に回折ベクトルVを加える具体的な方式の一例として、光像に応じた位相分布である距離分布r1(x,y)(第1の位相分布)に対し、光像とは無関係の距離分布r2(x,y)(第2の位相分布)を重畳する方式が考えられる。この場合、位相変調層15の距離分布r(x,y)は、
r(x,y)=r1(x,y)+r2(x,y)
として表される。r1(x,y)は、前述したように、光像をフーリエ変換したときの複素振幅の位相に相当する。r2(x,y)は、上記の数式(19)或いは数式(28)を満たす回折ベクトルVを加えるための距離分布である。なお、距離分布r2(x,y)の具体例は、図11と同様である。
[三次元計測装置の第1実施形態]
As an example of a specific method of adding the diffraction vector V to the in-plane wave vector K1 to K4, the optical image is obtained with respect to the distance distribution r1 (x, y) (first phase distribution) which is the phase distribution according to the optical image. A method of superimposing a distance distribution r2 (x, y) (second phase distribution) unrelated to the above can be considered. In this case, the distance distribution r (x, y) of the
r (x, y) = r1 (x, y) + r2 (x, y)
It is expressed as. As described above, r1 (x, y) corresponds to the phase of the complex amplitude when the optical image is Fourier transformed. r2 (x, y) is a distance distribution for adding a diffraction vector V that satisfies the above formula (19) or formula (28). A specific example of the distance distribution r2 (x, y) is the same as in FIG.
[First Embodiment of 3D Measuring Device]
図15は、第1実施形態に係る三次元計測装置101Aの構成を示す概略的な図である。同図に示すように、三次元計測装置101Aは、単体の光源部102と、複数(一対)の撮像部103と、計測部104とを含んで構成されている。光源部102は、上述したM点発振のS-iPMSEL1によって構成されている。光源部102から出射される計測光105は、ステージ106上に載置された被計測物SAの表面の一定の領域に照射される。ステージ106は、2次元方向又は3次元方向に走査可能な走査ステージであってもよい。なお、計測光105の照射範囲が被計測物SAの測定範囲に対して十分に広い場合、ステージ106の配置を省略してもよい。
FIG. 15 is a schematic diagram showing the configuration of the three-
本実施形態では、計測光105の所定パターンは、ドットパターン、ストライプパターン、及び格子パターンのいずれかからなる周期パターンW1となっている。図16の例では、計測光105の周期パターンW1は、100×100ピクセルの画像領域で示す周期的なドットパターンとなっている。このドットパターンでは、マトリクス状にドットが配列されており、ドット周期は、縦横共に5ピクセル周期となっている。また、図17は、周期的パターンの遠視野像の一例を示す図である。図17(a)は、40×40ドット、図17(b)は、60×60ドット、図17(c)は、80×80ドット、図17(d)は、120×120ドットの遠視野像である。光源部102の駆動条件は、電流0.5A、パルス幅50ns、パルス間隔5μs、温度25℃である。図の中心は、計測光105の面垂直方向の中心であり、図中のスケールバーは、15°に対応している。これらの図で示される遠視野像は、平面スクリーンにおいて各ドットがマトリクス状に配列するように設計されており、中心から離れた部分の配列の歪みは、測定系の光学系に起因するものである。
In the present embodiment, the predetermined pattern of the
撮像部103は、光源部102から出射される計測光105に対して感度を有する装置によって構成されている。撮像部103としては、例えばCCD(Charge Coupled Device)カメラ、CMOS(ComplementaryMOS)カメラ、その他の二次元イメージセンサなどを用いることができる。撮像部103は、計測光105が照射された状態の被計測物SAを撮像し、撮像結果を示す出力信号を計測部104に出力する。
The
計測部104は、例えばプロセッサ、メモリ等を含んで構成されるコンピュータシステムによって構成されている。計測部104は、各種の制御機能をプロセッサによって実行する。コンピュータシステムとしては、例えばパーソナルコンピュータ、マイクロコンピュータ、クラウドサーバ、スマートデバイス(スマートフォン、タブレット端末など)などが挙げられる。計測部104は、PLC(programmable logic controller)によって構成されていてもよく、FPGA(Field-programmable gate array)等の集積回路によって構成されていてもよい。
The measuring
計測部104は、撮像部103と通信可能に接続されており、撮像部103から入力される出力信号に基づいて、被計測物SAの三次元形状計測を実施する。本実施形態では、計測部104は、周期パターンW1を用いたアクティブステレオ法に基づいて被計測物SAの三次元形状を計測する。ここでは、一例として、平行等位ステレオの原理に基づく3次元形状計測法を示す。一対の撮像部103,103による視差をDとし、一対の撮像部103,103間の距離をb、一対の撮像部103,103の焦点距離をf、一対の撮像部103,103から被計測物SAまでの距離をZとすると、視差Dは、D=(f/Z)bで与えられる。撮像部103,103間の距離b及び撮像部103,103の焦点距離は、いずれも固有の値であるため、視差Dを求めることにより、被計測物SAまでの距離Zを求めることができる。
The
本実施形態では、周期パターンW1を有する計測光105が被計測物SAに照射される。このとき、撮像部103,103のそれぞれで撮像した周期パターンW1の同一点を計測部104で判別することが可能となる。また、パッシブステレオ法において課題となっていた、テクスチャの少ない画像を用いた三次元計測、暗部での三次元計測が可能となる。周期的なドットに代表される周期パターンW1を用いることで、計測光105のパターン密度の偏りが抑制され、計測光105の照明位置による計測精度のムラを抑えることが可能となる。
In the present embodiment, the
本実施形態において、周期パターンW1に代えて、例えば図18に示すようなランダムドットパターンW2を用いてもよい。ランダムドットパターンW2は、図16に示したドットパターンの各ドットを、格子点の位置から基本周期領域(隣接する格子点間の中点に垂直な線分で囲まれた矩形領域)の範囲で2次元的にランダムにシフトさせたパターンとなっている。一例として、格子点に位置するドットにそれぞれ乱数φ(ix,iy)を割り当て、当該乱数φに基づいて、各ドットを格子点の位置からシフトさせてもよい。 In the present embodiment, for example, a random dot pattern W2 as shown in FIG. 18 may be used instead of the periodic pattern W1. The random dot pattern W2 sets each dot of the dot pattern shown in FIG. 16 in the range from the position of the lattice point to the basic period region (a rectangular area surrounded by a line segment perpendicular to the midpoint between adjacent lattice points). The pattern is two-dimensionally randomly shifted. As an example, a random number φ (ix, ii) may be assigned to each dot located at the grid point, and each dot may be shifted from the position of the grid point based on the random number φ.
この場合、ランダムドットパターンW2が疑似的な周期性を有することとなるため、計測光105のパターン密度の偏りが抑制され、計測光105の照明位置による計測精度のムラを抑えることが可能となる。また、周期的なドットパターンではなく、ランダムドットパターンW2を用いることで、ドットパターンの同一点を異なる撮像部103で撮像した際の誤認識を抑制できる。したがって、視差Dの計測精度を向上でき、三次元形状計測の精度を高めることが可能となる。
In this case, since the random dot pattern W2 has a pseudo periodicity, the deviation of the pattern density of the
また、本実施形態において、周期パターンW1に代えて、図19に示すような均一な密度を有するパターンW3を用いてもよい。S-iPMSEL1からの出射光は、レーザ光であるため、散乱光中にスペックルが生じ得る。また、位相計算において、意図しないスペックル状ノイズの混入が生じる場合もある。したがって、均一な密度を有するパターンW3を用いた場合であっても、計測光105のパターンにランダムドットパターンが形成される。このランダムドットパターンを用いることで、ドットパターンの同一点を異なる撮像部103で撮像した際の誤認識を抑制できる。
Further, in the present embodiment, instead of the periodic pattern W1, a pattern W3 having a uniform density as shown in FIG. 19 may be used. Since the light emitted from the S-iPMSEL1 is a laser light, speckles may occur in the scattered light. In addition, unintended speckle-like noise may be mixed in the phase calculation. Therefore, even when the pattern W3 having a uniform density is used, a random dot pattern is formed in the pattern of the
図20は、均一な密度を有するパターンW3のFFP(Far Field Pattern)の一例を示す図である。同図の例では、計測光105のパルス幅を50ns、繰り返し間隔5μsとし、常温下でFFPを観測している。また、明るさ+40%、コントラスト-40%の色調補正を加えている。同図において、均一な密度を有するパターンW3を用いた場合であっても、計測光105のパターンにランダムドットパターンが形成されることが確認できる。
FIG. 20 is a diagram showing an example of an FFP (Far Field Pattern) of a pattern W3 having a uniform density. In the example of the figure, the pulse width of the
なお、図15の例では、三次元計測装置101Aは、単体の光源部102を備えているが、三次元計測装置101Aは、複数の光源部102を備えていてもよい。この場合、各光源部102からの計測光105を被計測物SAの異なる領域に照射することで、ステージ106の走査を行うことなく計測領域を拡大することができる。この構成を採用する場合には、ステージ106の配置を省略してもよい。
[三次元計測装置の第2実施形態]
In the example of FIG. 15, the three-
[Second Embodiment of 3D Measuring Device]
図21は、第2実施形態に係る三次元計測装置の構成を示す概略的な図である。同図に示すように、第2実施形態に係る三次元計測装置101Bは、複数の光源部102と、単体の撮像部103と、計測部104とを含んで構成されている。撮像部103、光源部102、及び計測部104の構成は、第1実施形態と同様である。本実施形態では、計測光105の所定パターンは、グレイコードパターンであり、計測部104は、グレイコードパターンを用いた三角測量法に基づいて被計測物SAの三次元形状を計測する。
FIG. 21 is a schematic diagram showing the configuration of the three-dimensional measuring device according to the second embodiment. As shown in the figure, the three-
図22は、グレイコードパターンの一例を示す図である。同図の例では撮像部103の画素をNx×Nyとし、X方向について示している。X方向の画素位置n(nは0~Nx-1の整数)をMx桁の2進数とすると、グレイコードパターンW4は、対象の数の2進数表現と、当該対象の数の2進数表現を1ビット右にシフトし、その先頭に0を付した数との排他的論理和で表される。すなわち、対象の数をnとすると、グレイコードパターンW4は、n^(n>>1)との論理式で与えられる。図22の例では、4ビット(4パターン)の場合のグレイコードパターンW4a~W4dを示している。グレイコードパターンW4の生成には、例えばOpenCVなどを用いることができる。
FIG. 22 is a diagram showing an example of a Gray code pattern. In the example of the figure, the pixels of the
グレイコードでは、隣接する画素のハミング距離が1となる。ハミング距離とは、桁数が同じ2つの値を比べたときに、対応する位置にある異なった値の桁の個数を指す。したがって、ハミング距離が1であるグレイコードでは、ビット列を復元する際にビットエラーが生じた場合でも誤差は1に収まる。単純なバイナリコードでは、上位ビットにエラーが生じた場合の位置の誤差が大きくなるが、グレイコードでは、ノイズに強い符号が得られる。 In Gray code, the Hamming distance of adjacent pixels is 1. The Hamming distance refers to the number of digits of different values at the corresponding positions when comparing two values with the same number of digits. Therefore, in the Gray code where the Hamming distance is 1, the error is within 1 even if a bit error occurs when restoring the bit string. In simple binary code, the position error becomes large when an error occurs in the high-order bit, but in Gray code, a code that is resistant to noise can be obtained.
グレイコードを用いる場合、光源部102の配置数は、2進数の各桁に対応するパターンの数であればよい。すなわち、グレイコードパターンW4a~W4dは、最上位ビットから最下位ビットまでの各桁の各画素の0,1が互いに異なるように設定された複数の縞状のパターンによって構成される。光源部102において最上位ビットのグレイコードパターンW4aから最下位ビットのグレイコードパターンW4dまでの各パターンを順に切り替えながら撮像部103で撮像を行う場合、Mx回の撮像で値Xが得られる。この値Xに基づいて、X番目の画素の位置を計測していることが分かる。Y方向についても同様に、グレイコードパターンW4a~W4dを順に切り替えながら撮像部103で撮像を行うことで、My回の撮像で値Yが得られる。この値Yに基づいて、Y番目の画素の位置を計測していることが分かる。
When the Gray code is used, the number of arrangements of the
被計測物SAの表面の色による誤認識を避けるため、図22に示したグレイコードパターンW4a~W4dとは、白黒が反転したグレイコードパターンを合わせて用いてもよい。この場合には、光源部102の配置数を2Mx+2Myとすればよい。
In order to avoid erroneous recognition due to the color of the surface of the object to be measured SA, the gray code patterns W4a to W4d shown in FIG. 22 may be used in combination with the gray code pattern in which black and white are reversed. In this case, the number of arrangements of the
本実施形態において、グレイコードパターンW4に代えて、例えば図23に示すように、正弦波状のストライプパターンW5を用いてもよい。図23に示す正弦波状のストライプパターンW5は、100×100ピクセルの画像領域で示す周期的なストライプパターンとなっている。正弦波状のストライプパターンW5の周期は、20ピクセル周期となっている。計測部104は、正弦波状のストライプパターンW5を用いた位相シフト法に基づいて被計測物SAの三次元形状を計測する。この形態では、例えば格子ピッチに対して1周期分を等分した位相シフト(位置ずれ)がそれぞれ与えられた複数の正弦波状のストライプパターンW5が用いられる。位相シフトのパターンは、位相が2π/N(Nは整数)ずつずれたものを用意すればよい。
In the present embodiment, instead of the Gray code pattern W4, for example, as shown in FIG. 23, a sinusoidal stripe pattern W5 may be used. The sinusoidal stripe pattern W5 shown in FIG. 23 is a periodic stripe pattern shown in an image area of 100 × 100 pixels. The period of the sinusoidal stripe pattern W5 is a period of 20 pixels. The measuring
ここでは、異なる位相シフトを有する4つの正弦波状のストライプパターンW5を用いる場合を例示する。4の正弦波状のストライプパターンW5を有する計測光105の光強度をそれぞれI0~I3とし、撮像部103の画素を(x,y)とすると、被計測物SAの表面での光強度I0~I3は、下記式(30)~(33)で表される。Ia(x,y)は、格子模様の振幅、Ib(x,y)は、背景強度、θ(x,y)は、初期位相である。
初期位相θは、tanθ=-(I3-I1)/(I2-I0)によって求めることができる。正弦波状のストライプパターンW5の位相シフト数がNである場合、初期位相θは、下記式(34)により求めることができる。
このような位相シフト法を用いる場合、計測した位相を高さ換算することで、正弦波状のストライプパターンW5のピッチよりも小さい間隔で被計測物SAの高さを計測することができる。三次元計測装置101Bの構成にあたっては、正弦波状のストライプパターンW5におけるストライプと平行な方向に光源部102を配列してもよい。この場合、光源部102の位置ずれに起因する位相シフトを無くすことが可能となり、複数の正弦波状のストライプパターンW5のそれぞれにおける初期位相のずれを解消できる。
When such a phase shift method is used, the height of the object to be measured SA can be measured at intervals smaller than the pitch of the sinusoidal stripe pattern W5 by converting the measured phase into height. In the configuration of the three-
本実施形態において、光源部102を互いに直交する二軸方向に配列してもよい。この場合、軸毎に計測光105のオン・オフを切り替えることで、被計測物SAの高さプロファイルを2軸で取得することができる。また、正弦波状のストライプパターンW5に代えて、例えば図24に示すように、互いに直交する二軸方向について正弦波状に変化するマトリクスパターンW6を用いてもよい。このようなマトリクスパターンW6を用いる場合、被計測物SAの高さプロファイルを2軸方向に同時に計測することができる。
In the present embodiment, the
本実施形態において、複数の光源部102は、互いに周期の異なる正弦波状のストライプパターンW5をそれぞれ出力してもよい。上述した位相シフト法では、位相2πにおける不連続性が課題となっている。これに対し、互いに周期の異なる正弦波状のストライプパターンW5を用いる場合には、例えば図25に示すように、全ての周波数で一致する座標を選択することにより、位相2πにおける不連続性を改善することが可能となり、少数のパターンで高精度の三次元計測を実現できる。位相2πにおける不連続性を改善することで、三次元形状計測の計測レンジの拡張や、凹凸の顕著な被計測物SAの高精度な計測を実現できる。
In the present embodiment, the plurality of
本実施形態において、計測部104は、正弦波状のストライプパターンW5を用いたサンプリングモアレ法に基づいて被計測物SAの三次元形状を計測してもよい。サンプリングモアレ法では、被計測物SAの表面に投影された正弦波状のストライプパターンW5の格子が被計測物SAの高さに応じて変形することを利用する。ここでは、撮像部103で撮像された画像において、基準面の高さでの位相シフト数Nの1つの正弦波パターンの縞間隔がカメラのN画素と対応するように予め調整する。ここでは、位相シフト数Nを4とする。1つの正弦波パターンを照射し、撮像部103の画素をN=4画素毎にサンプリングすることで、図26(a)に示すように、4画素毎に撮像された(撮像画素間の3画素が間引きされた)4つのパターンP1~P4が得られる。これらのパターンP1~P4間では、撮像画素が1画素ずつシフトしており、撮像画素の輝度値を線形補完することで、図26(b)に示すように、互いに位相がシフトしたモアレ縞パターンM1~M4が得られる。これらのモアレ縞パターンM1~M4を用いて上述した位相シフト法を適用することにより、正弦波状のストライプパターンW5のピッチよりも小さい間隔で被計測物SAの高さを計測することができる。このような手法によれば、前述の位相シフト法に比べて、照射する正弦波パターンの数を少なくすることができ、光源部102をコンパクトにすることができる。
In the present embodiment, the measuring
本実施形態において、正弦波状のストライプパターンW5に代えて、例えば図27に示すように、正弦波状のストライプパターンW5と、ランダムドットパターンW2とを重畳させた重畳パターンW7を用いてもよい。このような重畳パターンW7を用いることで、正弦波状のストライプパターンW5のピッチよりも小さい間隔で被計測物SAの高さを計測することができる。また、ランダムドットパターンを組み合わせることで、位相2πにおける不連続性を改善することが可能となり、少数のパターンで高精度の三次元計測を実現できる。重畳パターンは、図28に示すように、正弦波状に変化するマトリクスパターンW6と、ランダムドットパターンW2とを重畳させた重畳パターンW8であってもよい。この場合、前述の効果に加え、被計測物SAの高さプロファイルを2軸方向に同時に計測することができる。 In the present embodiment, instead of the sinusoidal stripe pattern W5, for example, as shown in FIG. 27, a superimposition pattern W7 in which the sinusoidal stripe pattern W5 and the random dot pattern W2 are superimposed may be used. By using such a superposed pattern W7, the height of the object to be measured SA can be measured at intervals smaller than the pitch of the sinusoidal stripe pattern W5. Further, by combining random dot patterns, it is possible to improve the discontinuity in the phase 2π, and it is possible to realize highly accurate three-dimensional measurement with a small number of patterns. As shown in FIG. 28, the superimposition pattern may be a superimposition pattern W8 in which a matrix pattern W6 that changes in a sinusoidal shape and a random dot pattern W2 are superposed. In this case, in addition to the above-mentioned effect, the height profile of the object to be measured SA can be measured simultaneously in the biaxial directions.
本実施形態において、正弦波状のストライプパターンW5及びグレイコードパターンW4の双方を用いてもよい。この場合、計測部104は、正弦波状のストライプパターンW5を用いた位相シフト法、及びグレイコードパターンW4を用いた三角測量法に基づいて被計測物SAの三次元形状を計測する。この場合、グレイコードパターンW4を用いた三角測量法によってピクセルレベルの計測を行い、正弦波状のストライプパターンW5を用いた位相シフト法によってサブピクセルレベルの計測を行うことができる。また、グレイコードを用いることで、位相2πにおける不連続性を改善することが可能となり、少数のパターンで高精度の三次元計測を実現できる。
[三次元計測装置の第3実施形態]
In this embodiment, both the sinusoidal stripe pattern W5 and the Gray code pattern W4 may be used. In this case, the measuring
[Third Embodiment of a three-dimensional measuring device]
図29は、第3実施形態に係る三次元計測装置の構成を示す概略的な図である。同図に示すように、第3実施形態に係る三次元計測装置101Cは、複数の光源部102と、複数(一対)の撮像部103と、計測部104とを含んで構成されている。撮像部103、光源部102、及び計測部104の構成は、第1実施形態と同様である。本実施形態では、計測光105の所定パターンは、正弦波状のストライプパターンW5であり、計測部104は、正弦波状のストライプパターンW5を用いた位相シフト法及びアクティブステレオ法に基づいて被計測物SAの三次元形状を計測する。
FIG. 29 is a schematic diagram showing the configuration of the three-dimensional measuring device according to the third embodiment. As shown in the figure, the three-
本実施形態では、計測した位相を高さ換算することで、正弦波状のストライプパターンW5のピッチよりも小さい間隔で被計測物SAの高さを計測することができる。また、複数の撮像部103を用いたアクティブステレオ法を組み合わせることで、位相2πにおける不連続性を改善することが可能となり、少数のパターンで高精度の三次元計測を実現できる。なお、アクティブステレオ法を用いる場合、前述したドットパターンを切り替えながら用いてもよい。
[S-iPMSELによる正弦波状のストライプパターンの位相シフト]
In the present embodiment, the height of the object to be measured SA can be measured at intervals smaller than the pitch of the sinusoidal stripe pattern W5 by converting the measured phase into height. Further, by combining the active stereo method using a plurality of
[Phase shift of sinusoidal stripe pattern by S-iPMSEL]
S-iPMSEL1では、設計した1次光のほか、出射面の法線について対称な-1次光が出力される(図6参照)。このため、正弦波状のストライプパターンW5の位相シフトを行う場合、±1次光が重なった正弦波を考えると、1次光と-1次光との間でストライプがシフトする向きが反転し、パターンが設計からずれてしまうおそれがある。説明の簡単化のため、X軸方向へのストライプのシフトを考えると、1次光の複素振幅は、下記式(35)で表される。-1次光の複素振幅は、面法線に対して1次光と対称となる位置に出射するものであり、下記式(36)で表される。式中、k(=kx、ky,kz)は波数ベクトル(大きさ2π/λ)、λは波長、ωは光の各周波数、Δθは位相シフト、a1は1次光振幅(理想的な位相分布に対し、実際の孔配置の位相分布に起因する成分)、Kは正弦波の縞の波数(=2π/Λ(Λは正弦波の周期))、θは正弦波の位相シフト量、(x,y,z)は投影ビームの座標である。
このとき、1次光及び-1次光の振幅に基づいて、合成振幅Aを下記式(37)により求めることができる。
実際の光強度は、合成振幅Aの二乗に比例するため、下記式(38)により求めることができる。
基本光波の周期は、正弦波の周期に比べて十分に小さい(λ<<Λ)。したがって、基本光波の波数kは、正弦波の縞の波数Kに比べて十分に大きい(k>>K)。このため、上記式(38)において、kの変化に対応する項を平均化してもよいと考えられる。この場合、±1次光を重畳した光の強度Iは、下記式(39)により近似できる。
Since the actual light intensity is proportional to the square of the combined amplitude A, it can be obtained by the following equation (38).
The period of the fundamental light wave is sufficiently smaller than the period of the sine wave (λ << Λ). Therefore, the wave number k of the fundamental light wave is sufficiently larger than the wave number K of the sine wave stripe (k >> K). Therefore, in the above equation (38), it is considered that the terms corresponding to the changes in k may be averaged. In this case, the intensity I of the light on which ± primary light is superimposed can be approximated by the following equation (39).
これらの式から、正弦波状のストライプパターンW5の位相シフトを行う場合、ストライプは、±1次光が重なったとしても、正弦波形状を維持したままシフトすることが分かる。また、±1次光が重畳した正弦波パターンでは、設計パターン(1次光振幅)に対して実際に得られる光強度(±1次光の振幅の和の二乗)の縞の間隔が半減し、位相シフト量が1周期に対して2倍となることが分かる。このため、例えば最終的に得られる光強度においてπ/2の位相シフトを実現させる場合、1次光振幅の位相シフト量の設計値をπ/4とすればよい。2軸方向に周期を持つ正弦波状のマトリクスパターンW6についても同様である。 From these equations, it can be seen that when the phase shift of the sinusoidal stripe pattern W5 is performed, the stripes shift while maintaining the sinusoidal shape even if the ± primary light overlaps. In addition, in the sine wave pattern in which ± primary light is superimposed, the interval between the stripes of the light intensity (the square of the sum of the amplitudes of ± primary light) actually obtained is halved with respect to the design pattern (primary light amplitude). It can be seen that the phase shift amount is doubled for one cycle. Therefore, for example, when a phase shift of π / 2 is realized in the finally obtained light intensity, the design value of the phase shift amount of the primary light amplitude may be set to π / 4. The same applies to the sinusoidal matrix pattern W6 having a period in the biaxial direction.
図4に示したように、S-iPMSEL1における異屈折率領域15bの重心Gを格子点O回りの円周方向にずらすことによって形成した場合、1次光及び-1次光の振幅aは、互いに等しい値となる。一方、図14に示したように、S-iPMSEL1における異屈折率領域15bの重心Gを、格子点Oを通り且つ正方格子の各辺に対して傾斜する直線D上にずらすことによって形成した場合、1次光及び-1次光の振幅aは、互いに異なる値となる。いずれの場合においても、±1次光が重畳した正弦波パターンを用いることができる。
As shown in FIG. 4, when the center of gravity G of the different
1次光及び-1次光が非対称なパターンである場合には、1次光と-1次光とが重なってしまうと、設計したパターンが得られないという問題がある。このような問題の一例としては、1次光の輝点1点当たりの構造が非対称な拡がりを持ち、設計パターンが不明瞭になってしまうことが挙げられる。この場合、1次光の出射領域を立体角がπとなる領域に限定すればよい。例えば1次光の出射領域を第1象限及び第2象限に限定する場合、-1次光の出射領域は、第4象限及び第3象限となるため、1次光と-1次光とが重なってしまうことを回避できる。これにより、1次光と-1次光との重なりによる輝点の拡がりを抑制できる。位相シフト法によって格子パターンをシフトさせる場合、-1次光のシフト方向は、1次光のシフト方向に対して反転する。したがって、1次光及び-1次光のそれぞれの出射領域に合わせて位相シフト演算で得られる位相を反転させることが好適である。一方、上記の問題が生じない場合でも、1次光及び-1次光を重畳せずにノイズの少ない像が得られる場合には、±1次光の投影領域を重畳せずにそれぞれ用いてもよい。
[光源部及び撮像部の配置例]
When the primary light and the -1st order light have an asymmetric pattern, there is a problem that the designed pattern cannot be obtained if the primary light and the -1st order light overlap. As an example of such a problem, the structure per bright spot of the primary light has an asymmetric spread, and the design pattern becomes unclear. In this case, the emission region of the primary light may be limited to the region where the solid angle is π. For example, when the emission region of the primary light is limited to the first quadrant and the second quadrant, the emission region of the -1st quadrant is the 4th quadrant and the 3rd quadrant. It is possible to avoid overlapping. As a result, it is possible to suppress the spread of bright spots due to the overlap of the primary light and the -1st order light. When the grid pattern is shifted by the phase shift method, the shift direction of the -1st order light is inverted with respect to the shift direction of the primary light. Therefore, it is preferable to invert the phase obtained by the phase shift operation according to the respective emission regions of the primary light and the -1st order light. On the other hand, even if the above problem does not occur, if an image with less noise can be obtained without superimposing the primary light and the -1st order light, use the ± primary light projection area without superimposing them. May be good.
[Example of arrangement of light source unit and image pickup unit]
図30は、光源部及び撮像部の配置例を示す概略的な斜視図である。図30に示すように、三次元計測装置101を構成するにあたって、光源部102及び撮像部103は、立体物111の表面に配置されていてもよい。立体物111は、三次元計測装置101A~101Cのプローブに相当する部分を構成する。立体物111は、例えば金属或いは樹脂によって円筒形状に形成されている。立体物111は、剛性を有していてもよく、可撓性を有していてもよい。立体物111は、内部空間を有していてもよい。
FIG. 30 is a schematic perspective view showing an example of arrangement of the light source unit and the image pickup unit. As shown in FIG. 30, in constructing the three-dimensional measuring device 101, the
光源部102及び撮像部103は、円筒形状の立体物111の周面111aにおいて、周方向にそれぞれ一定の間隔(ここでは45°の位相角)をもって配置されている。図30の例では、立体物111の先端側から基端側にわたって、一方の撮像部103の群、光源部102の群、他方の撮像部103の群が一定の間隔をもって配置されている。立体物111を長手方向から見た場合、一方の撮像部103、光源部102、他方の撮像部103が一列に並んでおり、これらの組が被計測物SAに対する一つの計測領域を構成している。立体物111が内部空間を有する場合、光源部102及び撮像部103に対する配線等が当該内部空間に収容されていてもよい。なお、光源部102と撮像部103との配置間隔は、必ずしも等間隔でなくてもよい。また、被計測物SAに対する測定範囲が網羅される場合には、単体の光源部102と単体の撮像部103とが立体物111に配置されていてもよい。
The
図31は、光源部及び撮像部の別の配置例を示す概略的な斜視図である。図31の例では、立体物121は、球形状をなしており、例えば円筒形状の支持部122の先端部分に設けられている。光源部102及び撮像部103は、球形状の立体物121の球面121aにおいて、経度方向にそれぞれ一定の間隔(ここでは45°の位相角)をもって配置されている。一方の撮像部103の群、光源部102の群、他方の撮像部103の群は、立体物121の緯度方向にわたって一定の間隔をもって配置されている。そして、立体物121の経度方向に並ぶ一方の撮像部103、光源部102、他方の撮像部103の組が被計測物SAに対する一つの計測領域を構成している。立体物121が内部空間を有する場合、光源部102及び撮像部103に対する配線等が当該内部空間に収容されていてもよい。図30の場合と同様、光源部102と撮像部103との配置間隔は、必ずしも等間隔でなくてもよい。また、被計測物SAに対する測定範囲が網羅される場合には、単体の光源部102と単体の撮像部103とが立体物121に配置されていてもよい。
FIG. 31 is a schematic perspective view showing another arrangement example of the light source unit and the image pickup unit. In the example of FIG. 31, the three-
以上のような構成によれば、光源部102及び撮像部103が配置された立体物111,121を三次元計測装置101のプローブとして構成できる。このような立体物111,121を用いることで、光源部102及び撮像部103のそれぞれの組を互いに異なる方向に向けることができるため、被計測物SAの三次元形状計測を広い立体角で実施することが可能となる。また、例えば口腔検査、内視鏡検査、管の内部や壁の隙間といった狭小箇所の検査、家具や装置等の床下からの検査といった用途への適用や、ハンディタイプの三次元計測装置の構築が容易となる。
According to the above configuration, the three-
図23に、正弦波状のストライプパターンW5を示したが、かかるストライプパターンを形成するにあたっては、隣接するパターン間のノイズ(輝度揺らぎ)を低減することが重要となる。隣接するパターン間のノイズは、例えば位相シフト法を適用する際の位置の揺らぎの要因となり、測定精度に影響を与えることが考えられる。そこで、隣接するパターン間のノイズ低減を考慮したストライプパターンの形成を実現するにあたっては、例えば図32に示すように、1次元の多点パターンを出射するS-iPMSEL1と、1次元レンズ51とを組み合わせた構成を採用し得る。
FIG. 23 shows a sine and cosine-shaped stripe pattern W5, but in forming such a stripe pattern, it is important to reduce noise (luminance fluctuation) between adjacent patterns. Noise between adjacent patterns may cause, for example, position fluctuation when applying the phase shift method, and may affect the measurement accuracy. Therefore, in order to realize the formation of a stripe pattern in consideration of noise reduction between adjacent patterns, for example, as shown in FIG. 32, the S-iPMSEL1 that emits a one-dimensional multipoint pattern and the one-
図32の例では、1次元レンズ51は、1次元凹レンズ52である。1次元凹レンズ52の媒質は、例えばガラスである。1次元凹レンズ52の一方面52aは平坦面であり、他方面52bは凹面となっている。1次元凹レンズ52は、一方面52aをS-iPMSEL1を向けた状態でS-iPMSEL1の表面(レーザ光の出射面)に対して配置されている。1次元凹レンズ52は、S-iPMSEL1の表面に対して結合され、S-iPMSEL1と一体化していてもよい。1次元凹レンズ52のレンズ位相は、下記式(40)によって求められる。下記式(40)において、φはレンズ位相、λはレンズ媒質中のレーザ光の波長、fは焦点距離である。
S-iPMSEL1からの多点パターンのレーザ光Laは、図32及び図33(a)の例では、X方向に所定の間隔をもって配列されている。図32の例では、1次元凹レンズ52は、凹面がX軸方向に延在するように配置されている。1次元凹レンズ52を通過した多点パターンのレーザ光Laは、X軸方向には変化せず、Y軸方向にのみ拡張する。したがって、多点パターンのレーザ光Laを1次元凹レンズ52に通すことにより、図33(b)に示すように、Y方向に拡張したライン状のレーザ光LbがX軸方向に並ぶストライプパターンW11が得られる。
The multi-point pattern laser beam La from S-iPMSEL1 is arranged at predetermined intervals in the X direction in the examples of FIGS. 32 and 33 (a). In the example of FIG. 32, the one-dimensional
ストライプパターンをより正弦波状に近づける場合、例えば図34(a)に示すように、多点パターンのレーザ光Laを形成し、各レーザ光の輝度がX軸方向について正弦波状となるように制御する。このような多点パターンのレーザ光Laを1次元凹レンズ52に通すことにより、図34(b)に示すストライプパターンW12では、Y方向に拡張したライン状のレーザ光LbがX軸方向に並ぶと共に、各レーザ光Lbの輝度がX軸方向について正弦波状に変化した状態となる。
When the stripe pattern is made closer to a sinusoidal shape, for example, as shown in FIG. 34 (a), a multipoint pattern laser beam La is formed, and the brightness of each laser beam is controlled to be sinusoidal in the X-axis direction. .. By passing the laser beam La of such a multi-point pattern through the one-dimensional
図34(a)及び図35(a)では、多点パターンのレーザ光LaがX軸方向に一直線状に並んでいるが、各レーザ光Laは、必ずしも一直線状に並んでいなくてもよく、Y軸方向に周期的或いはランダムにずれていてもよい。1次元レンズ51は、多点パターンのレーザ光Laを1次元方向に拡張できるものであればよく、1次元凹レンズ52に限られず、ラインジェネレータとして機能するパウエルレンズやラインマンレンズであってもよい。1次元レンズ51は、例えばフレネルレンズ、マイクロレンズ、メタレンズといったフラットレンズであってもよい。
In FIGS. 34 (a) and 35 (a), the laser beams La of the multipoint pattern are aligned in a straight line in the X-axis direction, but the laser light Las do not necessarily have to be aligned in a straight line. , May be periodically or randomly deviated in the Y-axis direction. The one-
メタレンズを用いる場合、例えば図35(a)に示すような共鳴型のメタレンズ構造53Aを採用してもよく、例えば図35(b)に示すような屈折率変調型のメタレンズ構造53Bを採用してもよい。図35(a)に示すように、共鳴型のメタレンズ構造53Aを採用する場合、メタレンズ構造53Aの構成材料は、下地となる層よりも高い屈折率を有する材料である。例えば下地となる層(例えば反射防止膜19)がSiNである場合、メタレンズ構造53Aの構成材料としては、アモルファスシリコンを用いることができる。メタレンズ構造53Aを構成する単位格子の高さ及び直径は、上述した式(40)で求まるレンズ位相に基づいて設定される。
When a metal lens is used, for example, a resonance type
図35(b)に示すように、屈折率変調型のメタレンズ構造53Bを採用する場合、S-iPMSEL1の表面のエッチングによってメタレンズ構造53Bを形成できる。例えばS-iPMSEL1の表面において、最表面となる層(例えば反射防止膜19)からその下層(例えば半導体基板10)の途中に至る孔部54をエッチングによって形成することで、屈折率変調型のメタレンズ構造53Bを形成できる。メタレンズ構造53Bを構成する各孔部54の深さ及び直径は、上述した式(40)で求まるレンズ位相に基づいて設定される。
As shown in FIG. 35B, when the refractive index modulation type
1…S-iPMSEL、101(101A~101C)…三次元計測装置、102…光源部、103…撮像部、104…計測部、105…計測光、SA…被計測物、W1…周期パターン、W2…ランダムドットパターン、W3…均一な密度を有するパターン、W4…グレイコードパターン、W5,W6…正弦波状のストライプパターン、W7,W8…重畳パターン、W11,W12…正弦波状のストライプパターン、111,121…立体物。 1 ... S-iPMSEL, 101 (101A to 101C) ... 3D measuring device, 102 ... light source unit, 103 ... imaging unit, 104 ... measuring unit, 105 ... measured light, SA ... object to be measured, W1 ... periodic pattern, W2 ... Random dot pattern, W3 ... Pattern with uniform density, W4 ... Gray code pattern, W5, W6 ... Sine wave-shaped stripe pattern, W7, W8 ... Superimposition pattern, W11, W12 ... Sine wave-shaped stripe pattern, 111, 121 … Three-dimensional object.
Claims (12)
前記計測光が照射された前記被計測物を撮像する一又は複数の撮像部と、
前記撮像部による撮像結果に基づいて前記被計測物の三次元形状を計測する計測部と、を備え、
前記光源部は、M点発振のS-iPMSELによって構成されている三次元計測装置。 One or more light source units that irradiate the object to be measured with measurement light having a predetermined pattern,
One or more image pickup units that image the object to be measured irradiated with the measurement light, and
A measurement unit that measures the three-dimensional shape of the object to be measured based on the image pickup result by the image pickup unit is provided.
The light source unit is a three-dimensional measuring device configured by S-iPMSEL of M point oscillation.
前記計測光の前記所定パターンは、ドットパターン、ストライプパターン、及び格子パターンのいずれかからなる周期パターンであり、
前記計測部は、前記周期パターンを用いたアクティブステレオ法に基づいて前記被計測物の三次元形状を計測する請求項1記載の三次元計測装置。 A single light source unit and a plurality of image pickup units are included.
The predetermined pattern of the measurement light is a periodic pattern including any one of a dot pattern, a stripe pattern, and a grid pattern.
The three-dimensional measuring device according to claim 1, wherein the measuring unit measures the three-dimensional shape of the object to be measured based on the active stereo method using the periodic pattern.
前記計測光の前記所定パターンは、ランダムドットパターンであり、
前記計測部は、前記ランダムドットパターンを用いたアクティブステレオ法に基づいて前記被計測物の三次元形状を計測する請求項1記載の三次元計測装置。 A single light source unit and a plurality of image pickup units are included.
The predetermined pattern of the measurement light is a random dot pattern.
The three-dimensional measuring device according to claim 1, wherein the measuring unit measures a three-dimensional shape of the object to be measured based on an active stereo method using the random dot pattern.
前記計測光の前記所定パターンは、均一な密度を有するパターンであり、
前記計測部は、前記均一な密度を有するパターンを用いたアクティブステレオ法に基づいて前記被計測物の三次元形状を計測する請求項1記載の三次元計測装置。 A single light source unit and a plurality of image pickup units are included.
The predetermined pattern of the measurement light is a pattern having a uniform density.
The three-dimensional measuring device according to claim 1, wherein the measuring unit measures a three-dimensional shape of the object to be measured based on an active stereo method using a pattern having a uniform density.
前記計測光の前記所定パターンは、グレイコードパターンであり、
前記計測部は、前記グレイコードパターンを用いた三角測量法に基づいて前記被計測物の三次元形状を計測する請求項1記載の三次元計測装置。 A plurality of the light source units and a single image pickup unit are included.
The predetermined pattern of the measurement light is a Gray code pattern.
The three-dimensional measuring device according to claim 1, wherein the measuring unit measures a three-dimensional shape of the object to be measured based on a triangulation method using the Gray code pattern.
前記計測光の前記所定パターンは、正弦波状のストライプパターンであり、
前記計測部は、前記正弦波状のストライプパターンを用いた位相シフト法に基づいて前記被計測物の三次元形状を計測する請求項1記載の三次元計測装置。 A plurality of the light source units and a single image pickup unit are included.
The predetermined pattern of the measurement light is a sinusoidal stripe pattern.
The three-dimensional measuring device according to claim 1, wherein the measuring unit measures the three-dimensional shape of the object to be measured based on a phase shift method using the sinusoidal stripe pattern.
前記計測光の前記所定パターンは、正弦波状のストライプパターンであり、
前記計測部は、前記正弦波状のストライプパターンを用いたサンプリングモアレ法に基づいて前記被計測物の三次元形状を計測する請求項1記載の三次元計測装置。 A plurality of the light source units and a single image pickup unit are included.
The predetermined pattern of the measurement light is a sinusoidal stripe pattern.
The three-dimensional measuring device according to claim 1, wherein the measuring unit measures a three-dimensional shape of the object to be measured based on a sampling moire method using the sinusoidal stripe pattern.
前記計測光の前記所定パターンは、正弦波状のストライプパターンと、ランダムドットパターンとを重畳させた重畳パターンであり、
前記計測部は、前記重畳パターンを用いた位相シフト法に基づいて前記被計測物の三次元形状を計測する請求項1記載の三次元計測装置。 A plurality of the light source units and a single image pickup unit are included.
The predetermined pattern of the measurement light is a superposition pattern in which a sinusoidal stripe pattern and a random dot pattern are superposed.
The three-dimensional measuring device according to claim 1, wherein the measuring unit measures the three-dimensional shape of the object to be measured based on a phase shift method using the superimposed pattern.
前記計測光の前記所定パターンは、正弦波状のストライプパターンと、グレイコードパターンとを含み、
前記計測部は、前記正弦波状のストライプパターンを用いた位相シフト法、及び前記グレイコードパターンを用いた三角測量法に基づいて前記被計測物の三次元形状を計測する請求項1記載の三次元計測装置。 A plurality of the light source units and a single image pickup unit are included.
The predetermined pattern of the measurement light includes a sinusoidal stripe pattern and a Gray code pattern.
The three-dimensional shape according to claim 1, wherein the measuring unit measures a three-dimensional shape of the object to be measured based on a phase shift method using the sinusoidal stripe pattern and a triangulation method using the Gray code pattern. Measuring device.
前記計測光の前記所定パターンは、正弦波状のストライプパターンであり、
前記計測部は、正弦波状のストライプパターンを用いた位相シフト法及びアクティブステレオ法に基づいて前記被計測物の三次元形状を計測する請求項1記載の三次元計測装置。 A plurality of the light source units and the plurality of image pickup units are included.
The predetermined pattern of the measurement light is a sinusoidal stripe pattern.
The three-dimensional measuring device according to claim 1, wherein the measuring unit measures the three-dimensional shape of the object to be measured based on a phase shift method using a sinusoidal stripe pattern and an active stereo method.
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US17/175,059 US20210262787A1 (en) | 2020-02-21 | 2021-02-12 | Three-dimensional measurement device |
CN202211213682.5A CN115388809A (en) | 2020-02-21 | 2021-02-18 | Three-dimensional measuring device |
CN202110187715.2A CN113295110A (en) | 2020-02-21 | 2021-02-18 | Three-dimensional measuring device |
JP2022129305A JP2022161985A (en) | 2020-02-21 | 2022-08-15 | Three-dimensional measurement device and light source device |
US18/195,508 US20230273015A1 (en) | 2020-02-21 | 2023-05-10 | Three-dimensional measurement device, light source device, and light receiving and emitting module |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020028446 | 2020-02-21 | ||
JP2020028446 | 2020-02-21 | ||
JP2020142715 | 2020-08-26 | ||
JP2020142715 | 2020-08-26 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2022129305A Division JP2022161985A (en) | 2020-02-21 | 2022-08-15 | Three-dimensional measurement device and light source device |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2022007906A true JP2022007906A (en) | 2022-01-13 |
JP7498651B2 JP7498651B2 (en) | 2024-06-12 |
Family
ID=80109939
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020198371A Active JP7498651B2 (en) | 2020-02-21 | 2020-11-30 | 3D measuring device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7498651B2 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023021803A1 (en) * | 2021-08-18 | 2023-02-23 | 浜松ホトニクス株式会社 | Method for designing phase modulation layer and method for producing light-emitting element |
WO2024053314A1 (en) * | 2022-09-05 | 2024-03-14 | 浜松ホトニクス株式会社 | Semiconductor light emitting device |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004251699A (en) | 2003-02-19 | 2004-09-09 | Ricoh Co Ltd | Shape measuring device and three-dimensional position sensor |
JP4831703B2 (en) | 2008-04-23 | 2011-12-07 | 国立大学法人 和歌山大学 | Object displacement measurement method |
JP6104662B2 (en) | 2013-03-25 | 2017-03-29 | 株式会社東芝 | Measuring device, method and program |
JP6788574B2 (en) | 2015-03-13 | 2020-11-25 | 浜松ホトニクス株式会社 | Semiconductor light emitting device |
JP7005175B2 (en) | 2017-05-30 | 2022-01-21 | キヤノン株式会社 | Distance measuring device, distance measuring method and imaging device |
JP7149707B2 (en) | 2018-01-15 | 2022-10-07 | キヤノン株式会社 | Information processing device, its control method and program, and operation control system |
JP6994165B2 (en) | 2018-03-16 | 2022-01-14 | 日本電気株式会社 | 3D shape measuring device, 3D shape measuring method and program |
JP7090447B2 (en) | 2018-03-22 | 2022-06-24 | 株式会社キーエンス | Image processing equipment |
-
2020
- 2020-11-30 JP JP2020198371A patent/JP7498651B2/en active Active
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023021803A1 (en) * | 2021-08-18 | 2023-02-23 | 浜松ホトニクス株式会社 | Method for designing phase modulation layer and method for producing light-emitting element |
WO2024053314A1 (en) * | 2022-09-05 | 2024-03-14 | 浜松ホトニクス株式会社 | Semiconductor light emitting device |
Also Published As
Publication number | Publication date |
---|---|
JP7498651B2 (en) | 2024-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20230392920A1 (en) | Multiple channel locating | |
US11166004B2 (en) | Three-dimensional computational imaging method and apparatus based on single-pixel sensor, and non-transitory computer-readable storage medium | |
BR112020025649A2 (en) | apparatus for intraoral scanning and method for generating a three-dimensional digital model of an intraoral surface | |
US11327297B2 (en) | Generation method for programmable analog fringe pattern with extended depth of field | |
WO2020264035A1 (en) | Intraoral 3d scanner employing multiple miniature cameras and multiple miniature pattern projectors | |
US20070115484A1 (en) | 3d shape measurement system and method including fast three-step phase shifting, error compensation and calibration | |
JP2020531856A (en) | Flight time depth camera with low resolution pixel imaging | |
JP2022007906A (en) | Three-dimensional measuring device | |
JP2011064482A (en) | Device and method of high-speed three-dimensional measurement | |
JP2004191092A (en) | Three-dimensional information acquisition system | |
CN105953749A (en) | Optimal three-dimensional morphology measurement method | |
CN113218336B (en) | Grating plate, projection device for three-dimensional reconstruction of structured light and measurement device | |
JP2022161985A (en) | Three-dimensional measurement device and light source device | |
US20230375747A1 (en) | Single nanostructure-integrated metalens | |
JP5667891B2 (en) | Shape measurement method | |
CN110375640A (en) | A kind of rotary phase shifting interferometer and measurement method measuring transparent substance | |
JP5853284B2 (en) | Shape measuring apparatus and shape measuring method | |
US20240162679A1 (en) | Light emission device and light source device | |
US11085761B2 (en) | Determining surface structures of objects | |
Lei et al. | A phase-angle inspired calibration strategy based on MEMS projector for 3D reconstruction with markedly reduced calibration images and parameters | |
JP2020118533A (en) | Method of measuring outer-of-plane displacement distribution and three-dimensional shape of measurement object and device using the same | |
US20230204345A1 (en) | Three-dimensional measurement device | |
KR20130028594A (en) | Method for creating structured light, apparatus thereof and aparatus for motion detection by using apparatus thereof | |
WO2024053314A1 (en) | Semiconductor light emitting device | |
Cheng et al. | Optimization of bit-pairing codification with learning for 3D reconstruction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230516 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20231227 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240109 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20240305 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240423 |
|
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: 20240528 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20240531 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7498651 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |