JP2009031229A - Measurement signal noise processing method - Google Patents
Measurement signal noise processing method Download PDFInfo
- Publication number
- JP2009031229A JP2009031229A JP2007198155A JP2007198155A JP2009031229A JP 2009031229 A JP2009031229 A JP 2009031229A JP 2007198155 A JP2007198155 A JP 2007198155A JP 2007198155 A JP2007198155 A JP 2007198155A JP 2009031229 A JP2009031229 A JP 2009031229A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- measurement
- value
- polarization
- noise
- 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
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
【課題】被計測物の構造、材質等に起因する計測上のばらつきを解決し、被計測物の厚さ等の計測誤差を小さくする。
【解決手段】通常の測定では、データを繰り返し計測し、その分布を求め平均を持って真の値とするが、ゼロ付近ではデータの値そのものが小さいので、元々ある測定ノイズなどに埋もれてしまい、ゼロ点では計測データの検出はできないので、分布がゼロになったり、データの分布は経験的に実際はなめらかになるので、真の値の周りに非対称に広がってしまい、単に移動平均では真の値を推定することができないが、未知のランダムノイズが重畳され対称でない観測信号の真の値を推定するために、観測信号の最大値をとることより、微小な(観測)信号領域での計測精度を改善する。
【選択図】図3An object of the present invention is to solve measurement variations caused by the structure and material of an object to be measured, and to reduce measurement errors such as the thickness of the object to be measured.
In normal measurement, data is measured repeatedly, and its distribution is obtained to obtain a true value with an average. However, since the data value itself is small near zero, it is buried in the original measurement noise. Since the measurement data cannot be detected at the zero point, the distribution becomes zero, or the data distribution is actually empirically smooth, so it spreads asymmetrically around the true value. The value cannot be estimated, but in order to estimate the true value of the observed signal with unknown random noise superimposed and not symmetric, the maximum value of the observed signal is taken, and measurement in the minute (observed) signal region is performed. Improve accuracy.
[Selection] Figure 3
Description
本発明は、計測機器(例.光コヒーレンストモグラフィー装置等)に利用される測定信号のノイズ処理方法に関する。 The present invention relates to a noise processing method for a measurement signal used in a measuring instrument (eg, an optical coherence tomography apparatus or the like).
従来、各種被計測物体の非破壊断層計測技術の1つとして、光断層画像装置「光コヒーレンストモグラフィー」(OCT)がある(特許文献1参照)。OCTは、光を計測プローブとして用いるため、被計測物体の屈折率分布、分光情報、偏光情報(複屈折率分布)等が計測できるという利点がある。 Conventionally, as one of non-destructive tomographic measurement techniques for various objects to be measured, there is an optical tomographic imaging apparatus “optical coherence tomography” (OCT) (see Patent Document 1). Since OCT uses light as a measurement probe, it has the advantage that it can measure the refractive index distribution, spectral information, polarization information (birefringence distribution), etc. of the measured object.
基本的なOCT43は、マイケルソン干渉計を基本としており、その原理を図5で説明する。光源44から射出された光は、コリメートレンズ45で平行化された後に、ビームスプリッター46により参照光と物体光に分割される。物体光は、物体アーム内の対物レンズ47によって被計測物体48に集光され、そこで散乱・反射された後に再び対物レンズ47、ビームスプリッター46に戻る。
The
一方、参照光は参照アーム内の対物レンズ49を通過した後に参照鏡50によって反射され、再び対物レンズ49を通してビームスプリッター46に戻る。このようにビームスプリッター46に戻った物体光と参照光は、物体光とともに集光レンズ51に入射し光検出器52(フォトダイオード等)に集光される。
On the other hand, the reference light passes through the
OCTの光源44は、時間的に低コヒーレンスな光(異なった時刻に光源から出た光同士は極めて干渉しにくい光)の光源を利用する。時間的低コヒーレンス光を光源としたマイケルソン型の干渉計では、参照アームと物体アームの距離がほぼ等しいときにのみ干渉信号が現れる。この結果、参照アームと物体アームの光路長差(τ)を変化させながら、光検出器52で干渉信号の強度を計測すると、光路長差に対する干渉信号(インターフェログラム)が得られる。
The
そのインターフェログラムの形状が、被計測物体48の奥行き方向の反射率分布を示しており、1次元の軸方向走査により被計測物体48の奥行き方向の構造を得ることができる。このように、OCT43では、光路長走査により、被計測物体48の奥行き方向の構造を計測できる。
The shape of the interferogram shows the reflectance distribution in the depth direction of the
このような軸方向の走査のほかに、横方向の機械的走査を加え、2次元の走査を行うことで被計測物体の2次元断面画像が得られる。この横方向の走査を行う走査装置としては、被計測物体を直接移動させる構成、物体は固定したままで対物レンズをシフトさせる構成、被計測物体も対物レンズも固定したままで、対物レンズの瞳面付近においたガルバノミラーの角度を回転させる構成等が用いられている。 In addition to the scanning in the axial direction, a two-dimensional cross-sectional image of the object to be measured can be obtained by performing a two-dimensional scanning by adding a horizontal mechanical scanning. The scanning device that performs the horizontal scanning includes a configuration in which the object to be measured is directly moved, a configuration in which the objective lens is shifted while the object is fixed, and a pupil of the objective lens while the object to be measured and the objective lens are fixed. The structure etc. which rotate the angle of the galvanometer mirror in the surface vicinity are used.
従来、OCT装置は、複屈折性を有する試料の構造を測定し、しかも試料の奥行方向に高精度な分解能が要求される各種の技術分野、例えば、半導体製品の製造等の工業分野、動物の生体観察や植物の構造観察のような動植物の研究・観察分野、各種文化財の解析・鑑定技術の分野、ロボット技術分野(植物、昆虫、動物等の各部器官を観察し、その構造・機能をロボット技術に応用するための技術)に有用である。 Conventionally, an OCT apparatus measures the structure of a sample having birefringence, and also requires various technical fields that require high-precision resolution in the depth direction of the sample, such as industrial fields such as semiconductor product manufacturing, animal Research and observation field of animals and plants such as biological observation and plant structure observation, field of analysis and appraisal technology of various cultural assets, robot technology field (observing each part organ such as plants, insects, animals, etc. This is useful for applications to robot technology).
しかしながら、被測定対象である物体における各部位における構造、材料、厚さの違い違いが計測上のバラツキとなり、厚さ等の計測誤差が増大してしまうという問題があった。 However, the difference in structure, material, and thickness at each part of the object to be measured has a variation in measurement, and there is a problem that measurement errors such as thickness increase.
本発明は、上記のようにOCT装置等の計測機器によって計測した際に生じる、測定信号、特に微小測定信号のノイズに起因する計測上のばらつきや誤差を少なくし、精密な測定を可能とする測定信号のノイズ処理方法を実現することを課題とする。 The present invention reduces measurement variations and errors caused by noise of measurement signals, particularly minute measurement signals, which occur when measurement is performed by a measuring instrument such as an OCT apparatus as described above, and enables precise measurement. It is an object to realize a noise processing method for a measurement signal.
本発明は上記課題を解決するために、計測機器による測定信号は、該測定信号に重畳されたノイズ信号が0付近であって対称ではない分布曲線を有する微小信号領域については、測定信号の最大値を真値として利用することを特徴とする測定信号のノイズ処理方法を提供する。 In order to solve the above-described problems, the present invention provides a measurement signal obtained by a measuring instrument, and a maximum signal of a measurement signal is measured in a minute signal region having a distribution curve in which a noise signal superimposed on the measurement signal is near 0 and is not symmetric. A noise processing method for a measurement signal, characterized in that a value is used as a true value.
本発明は上記課題を解決するために、計測機器による測定信号は、該測定信号の信号領域が限られている場合、真値が信号領域の端に近い値の場合、それに重畳されたノイズ信号によって微小測定信号が対称ではない分布曲線を有する信号領域については、測定信号の最大値を真値として利用することを特徴とする信号領域の限られた測定信号のノイズ処理方法を提供する。 In order to solve the above-described problems, the present invention provides a measurement signal from a measuring instrument in which the signal region of the measurement signal is limited, or when the true value is a value close to the end of the signal region, a noise signal superimposed thereon. For a signal region having a distribution curve in which a minute measurement signal is not symmetric, the maximum value of the measurement signal is used as a true value, and a noise processing method for a measurement signal having a limited signal region is provided.
本発明は上記課題を解決するために、計測機器による測定信号は、真値に重畳されたノイズ信号によって測定信号が対称ではない分布曲線を有する信号については、測定信号の最大値を真値として利用することを特徴とする測定信号のノイズ処理方法を提供する。 In order to solve the above-described problems, the present invention provides a measurement signal obtained by a measuring instrument that has a distribution curve in which the measurement signal is not symmetrical due to a noise signal superimposed on the true value. The present invention provides a noise processing method for a measurement signal, which is characterized by being utilized.
以上のような構成から成る本発明に係る測定信号のノイズ処理方法によれば、各種の計測機器による計測の際に生じる、測定信号、特に微小測定信号のノイズに起因する計測上のばらつきや誤差を少なくし、精密な測定を可能とする。 According to the measurement signal noise processing method according to the present invention having the above-described configuration, measurement variations and errors caused by noise of measurement signals, particularly minute measurement signals, that occur during measurement by various measuring devices. And precise measurement is possible.
例えば、光コヒーレンストモグラフィーの画像処理装置に適用すれば、光コヒーレンストモグラフィーによって得た動植物のいろいろな器官の構造解析のための計測において、位相ノイズを低減して明瞭なOCT画像を得ることが可能である。 For example, when applied to an image processing apparatus for optical coherence tomography, it is possible to obtain clear OCT images by reducing phase noise in measurement for structural analysis of various organs of animals and plants obtained by optical coherence tomography. is there.
本発明に係る測定信号のノイズ処理方法を実施するための最良の形態を、実施例に基づいて図面を参照して、以下に説明する。 The best mode for carrying out the noise processing method for measurement signals according to the present invention will be described below with reference to the drawings based on the embodiments.
この実施例1では、本発明に係る測定信号のノイズ処理方法を偏光感受型光コヒーレンストモグラフィー装置に適用した例で説明する。 In the first embodiment, the measurement signal noise processing method according to the present invention is applied to a polarization-sensitive optical coherence tomography apparatus.
「背景技術」の項で説明したとおり、光コヒーレンストモグラフィー装置(OCT)では、光源からのビームを参照アームと試料アームに分離して送り、試料アームでは試料(被検体)の深さ方向(A方向)に垂直な方向に走査(Bスキャン)して試料を照射し、この反射光と参照アームから反射される参照光との干渉スペクトルからA−B画像を得る(OCT計測を行う)ものであり、本発明は、FD−OCT(フーリエドメインOCT)等に適用される。 As described in the “Background Art” section, in the optical coherence tomography apparatus (OCT), the beam from the light source is sent separately to the reference arm and the sample arm, and in the sample arm, the depth direction (A The sample is irradiated in a direction perpendicular to (direction) (B-scan), and an AB image is obtained from the interference spectrum between the reflected light and the reference light reflected from the reference arm (performs OCT measurement). The present invention is applied to FD-OCT (Fourier domain OCT) and the like.
なお補足すると、試料の深さ(光軸)方向の走査(この走査を「A−スキャン」と言い、この方向を「A−方向」、「Aスキャン方向」とも言う。)は、一回の光照射により深さ方向の後方散乱データを取得することで可能であり、ガルバノ鏡によって横方向(A方向に垂直な方向)の走査(この走査を「B−スキャン」と言い、この方向を「B−方向」、「Bスキャン方向」とも言う。)を行うことで、2次元断層画像(偏光感受型のOCT画像)を得ることができる。 Supplementally, scanning in the depth (optical axis) direction of the sample (this scanning is referred to as “A-scan”, and this direction is also referred to as “A-direction”, “A scan direction”) is performed once. This is possible by acquiring backscattering data in the depth direction by light irradiation. Scanning in the lateral direction (direction perpendicular to the A direction) by a galvano mirror (this scanning is called “B-scan”). By performing “B-direction” and “B-scan direction”), a two-dimensional tomographic image (polarization-sensitive OCT image) can be obtained.
本発明に係る測定信号のノイズ処理方法を適用する偏光感受型光コヒーレンストモグラフィー装置は、Bスキャンと同時に(同期して)光源からの偏光ビーム(偏光子により直線的に偏光されたビーム)をEO変調器(偏光変調器、電気光学変調器)によって連続的に変調し、この連続的に偏光を変調した偏光ビームを分けて、一方を入射ビームとして走査して試料に照射し、その反射光(物体光)を得ると共に、他方を参照光として、両者のスペクトル干渉によりOCT計測を行うものである。 The polarization-sensitive optical coherence tomography apparatus to which the noise processing method of the measurement signal according to the present invention is applied is a EO-converted polarization beam (a beam linearly polarized by a polarizer) from a light source simultaneously (synchronously) with a B-scan. The light is continuously modulated by a modulator (polarization modulator, electro-optic modulator), the polarization beam is continuously modulated, and the sample is scanned as an incident beam to irradiate the sample, and the reflected light ( (Object light) and OCT measurement by spectral interference between the two using the other as reference light.
そして、このスペクトル干渉成分のうち、垂直偏光成分(H)と水平偏光成分(V)を同時に2つの光検出器で測定することにより、試料の偏光特性を表すジョーンズベクトルを得る(H画像とV画像)構成を特徴とするものである。 Of the spectral interference components, the vertical polarization component (H) and the horizontal polarization component (V) are simultaneously measured by two photodetectors, thereby obtaining a Jones vector representing the polarization characteristics of the sample (H image and V). Image) structure.
図1は、本発明に係る測定信号のノイズ処理方法を適用する偏光感受型光コヒーレンストモグラフィー装置の光学系の全体構成を示す図である。図1に示す偏光感受光画像計測装置1は、光源2、偏光子3、EO変調器4、ファイバーカプラー(光カプラー)5、参照アーム6、試料アーム7、分光器8等の光学要素を備えている。この偏光感受光画像計測装置1の光学系は、光学要素が互いにファイバー9で結合されているが、ファイバーで結合されていないタイプの構造(フリースペース型)であってもよい。 FIG. 1 is a diagram showing an overall configuration of an optical system of a polarization-sensitive optical coherence tomography apparatus to which a measurement signal noise processing method according to the present invention is applied. 1 includes optical elements such as a light source 2, a polarizer 3, an EO modulator 4, a fiber coupler (optical coupler) 5, a reference arm 6, a sample arm 7, a spectroscope 8, and the like. ing. The optical system of this polarization-sensitive received light image measuring apparatus 1 may have a structure of a type (free space type) in which optical elements are coupled to each other by a fiber 9 but not coupled to a fiber.
光源2は、広帯域スペクトルを有するスーパールミネッセントダイオード(SLD:Super Luminescent Diode)を使用する。なお、光源2は、パルスレーザでもよい。光源2には、コリメートレンズ11、光源2からの光を直線偏光にする偏光子3、進相軸を45°の方向にセットされたEO変調器4、集光レンズ13(場合によってはさらにレンズ13’が使用される。)及びファイバーカプラー5が、順次、接続されている。
The light source 2 uses a super luminescent diode (SLD) having a broadband spectrum. The light source 2 may be a pulse laser. The light source 2 includes a
EO変調器4は、進相軸を45°の方向に固定して、該EO変調器4にかける電圧を正弦的に変調することで、進相軸とそれに直交する遅相軸との間の位相差(リタデーション)を連続的に変えるもので、これにより、光源2から出て偏光子3で(縦)直線偏光となった光がEO変調器4に入射すると、上記変調の周期で、直線偏光→楕円偏光→直線偏光………などのように変調される。EO変調器4は、市販されているEO変調器を使用すればよい。 The EO modulator 4 fixes the fast axis in a 45 ° direction and modulates the voltage applied to the EO modulator 4 sinusoidally, so that the phase between the fast axis and the slow axis orthogonal thereto is obtained. The phase difference (retardation) is continuously changed. As a result, when light that has been emitted from the light source 2 and changed into (longitudinal) linearly polarized light by the polarizer 3 is incident on the EO modulator 4, the light is linearly changed at the above-described modulation period. Polarized light → elliptical polarized light → linearly polarized light, etc. As the EO modulator 4, a commercially available EO modulator may be used.
ファイバーカプラー5には分岐するファイバー9を介して、参照アーム6と試料アーム7が接続されている。参照アーム6には、偏波コントローラ(polarization controller)10、コリメートレンズ11、偏光子12、集光レンズ13及び参照鏡(固定鏡)14が、順次、設けられている。参照アーム6の偏光子12は、上記のとおり偏光状態を変調しても参照アーム6から戻ってくる光の強度が変化しないような方向を選択するために用いている。この偏光子12の方向(直線偏光の偏光方向)の調整は偏波コントローラ10とセットで行う。
A reference arm 6 and a sample arm 7 are connected to the fiber coupler 5 via a branching fiber 9. The reference arm 6 is sequentially provided with a polarization controller 10, a
試料アーム7では、偏波コントローラ15、コリメートレンズ11、固定鏡24、ガルバノ鏡16、集光レンズ13が、順次、設けられ、ファイバーカプラー5からの入射ビームが2軸のガルバノ鏡16により走査されて試料17に照射される。試料17からの反射光(後方散乱光)は物体光として再びファイバーカプラー5に戻り、参照光と重畳されて干渉ビームとして分光器8に送られる。
In the sample arm 7, a
分光器8は、順次接続される偏波コントローラ18、コリメートレンズ11、(偏光感受型体積位相ホログラフィック)回折格子19、フーリエ変換レンズ20、偏光ビームスプリッター21及び2つの光検出器22、23を備えている。この実施例では、光検出器22、23として、ラインCCDカメラ(1次元CCDカメラ)を利用する。ファイバーカプラー5から送られてくる干渉ビームは、コリメートレンズ11でコリメートされ、回折格子19によって干渉スペクトルに分光される。
The spectroscope 8 includes a polarization controller 18, a
回折格子19で分光された干渉スペクトルビームは、フーリエ変換レンズ20でフーリエ変換され偏光ビームスプリッター21で水平及び垂直成分に分けられ、それぞれ2つラインCCDカメラ(光検出器)22、23で検出される。この2つラインCCDカメラ22、23は、水平および垂直偏光信号両方の位相情報を検知するために使われるので、2つのラインCCDカメラ22、23は同一の分光器の形成に寄与するものでなくてはならない。
The interference spectrum beam split by the
なお、光源2、参照アーム6、試料アーム7及び分光器8には、それぞれ偏波コントローラ10、15、18が設けられているが、これらは、光源2から参照アーム6、試料アーム7、分光器8に送られるそれぞれのビームの初期偏光状態を調整して、EO変調器4で連続的に変調された偏光状態が、参照光と物体光においても互いに一定の振幅と一定の相対偏光状態の関係が維持され、さらにファイバーカプラー5に接続された分光器8において一定の振幅と一定の相対偏光状態を保たれるようにコントロールする。
The light source 2, the reference arm 6, the sample arm 7, and the spectrometer 8 are provided with
また、2つラインCCDカメラ22、23を含む分光器8を校正するときはEO変調器4を止める。参照光をブロックし、スライドガラスと反射鏡を試料アーム7におく。この配置は水平および垂直偏光成分のピークの位置が同じであることを保証する。そして、スライドガラスの後ろの面と反射鏡からのOCT信号は2つの分光器8で検知される。OCT信号のピークの位相差はモニターされる。 When the spectroscope 8 including the two line CCD cameras 22 and 23 is calibrated, the EO modulator 4 is stopped. The reference light is blocked, and the slide glass and the reflecting mirror are placed on the sample arm 7. This arrangement ensures that the positions of the peaks of the horizontal and vertical polarization components are the same. The OCT signals from the rear surface of the slide glass and the reflecting mirror are detected by the two spectrometers 8. The phase difference of the peak of the OCT signal is monitored.
この位相差はすべての光軸方向の深さでゼロであるべきである。次に、信号は2つラインCCDカメラ22、23を含む分光器8で複素スペクトルを得るために、ウィンドウされ逆フーリエ変換される。この位相差はすべての周波数でゼロであるべきなので、これらの値をモニターすることによって2つラインCCDカメラ22、23の物理的な位置は位相差が最小になるようにアライメントされる。 This phase difference should be zero at all optical axis depths. The signal is then windowed and inverse Fourier transformed to obtain a complex spectrum with a spectrometer 8 that includes two line CCD cameras 22,23. Since this phase difference should be zero at all frequencies, by monitoring these values, the physical positions of the two line CCD cameras 22, 23 are aligned so that the phase difference is minimized.
上記偏光感受光画像計測装置の動作は次のとおりである。光源2からの光を直線偏光し、この直線偏光されたビームをEO変調器4により連続的に偏光状態の変調を行う。即ち、EO変調器4は、進相軸を45°の方向に固定して、EO変調器4にかける電圧を正弦的に変調することで、進相軸とそれに直交する遅相軸との間の位相差(偏光角:リタデーション)を連続的に変えるもので、これにより、光源2から出て直線偏光子で(縦)直線偏光となった光がEO変調器4に入射すると、上記変調の周期で、直線偏光→楕円偏光→直線偏光………などのように変調される。 The operation of the polarized light-sensitive received image measuring device is as follows. The light from the light source 2 is linearly polarized, and the linearly polarized beam is continuously modulated in the polarization state by the EO modulator 4. That is, the EO modulator 4 fixes the fast axis in a 45 ° direction and modulates the voltage applied to the EO modulator 4 sinusoidally, so that the phase between the fast axis and the slow axis orthogonal thereto is obtained. The phase difference (polarization angle: retardation) of the light is continuously changed, so that when the light that has been emitted from the light source 2 and converted into (longitudinal) linearly polarized light by the linear polarizer enters the EO modulator 4, In a cycle, the light is modulated as follows: linearly polarized light → elliptical polarized light → linearly polarized light.
そして、直線偏光された偏光ビームをEO変調器4により連続的に偏光状態の変調を行うと同時に、Bスキャンを同期して行う。即ち、1回のBスキャンの間に、EO変調器4による偏光の連続的な変調を複数周期行う。ここで、1周期とは、偏光角(リターデーション)φが0〜2πと変化する期間である。要するに、この1周期の間に、偏光子からの光の偏光が、直線偏光(垂直偏光)→楕円偏光→直線偏光(水平偏光)………などのように連続的に変調する。 The linearly polarized polarized beam is continuously modulated in the polarization state by the EO modulator 4, and at the same time, the B scan is performed in synchronization. That is, during one B scan, the EO modulator 4 continuously modulates the polarization for a plurality of periods. Here, one period is a period during which the polarization angle (retardation) φ changes from 0 to 2π. In short, during this one cycle, the polarization of light from the polarizer is continuously modulated as linearly polarized light (vertical polarized light) → elliptical polarized light → linearly polarized light (horizontal polarized light).
このように偏光ビームの偏光を連続的に変調しながら、試料アーム7では、入射ビームをガルバノ鏡16により試料17に走査してBスキャンを行い、分光器8において、その反射光である物体光と参照光の干渉スペクトルについて、その水平偏光成分および垂直偏光成分を2つのラインCCDカメラ22、23で検出する。これにより、1回のBスキャンによって、それぞれ水平偏光成分及び垂直偏光成分に対応する2枚のA−Bスキャン画像が得られる。 In this way, while continuously modulating the polarization of the polarized beam, the sample arm 7 scans the incident beam on the sample 17 by the galvano mirror 16 to perform the B scan, and the spectroscope 8 performs the object light as the reflected light. As for the interference spectrum of the reference light, the horizontal polarization component and the vertical polarization component are detected by the two line CCD cameras 22 and 23. Thus, two AB scan images corresponding to the horizontal polarization component and the vertical polarization component are obtained by one B scan.
上記のとおり、1回のBスキャンの間に、偏光ビームの偏光の連続的な変調を複数周期行うが、各周期(1周期)の連続的な変調の間に2つのラインCCDカメラ22、23で検出した水平偏光成分および垂直偏光成分それぞれの偏光情報が1画素分の偏光情報となる。1周期の連続的な変調の間に2つのラインCCDカメラ22、23で偏光情報を検出タイミング信号に同期して行い、1周期に検出回数(取込回数)を、4回、8回等、適宜決めればよい。 As described above, during one B-scan, continuous modulation of the polarization of the polarized beam is performed for a plurality of periods, but the two line CCD cameras 22 and 23 are continuously modulated during each period (one period). The polarization information of each of the horizontal polarization component and the vertical polarization component detected in step 1 becomes polarization information for one pixel. Polarization information is synchronized with the detection timing signal by two line CCD cameras 22 and 23 during one period of continuous modulation, and the number of detections (number of acquisitions) is four times, eight times, etc. What is necessary is just to decide suitably.
このようにして1回のBスキャンの間に得た2枚のA−Bスキャン画像のデータを、Bスキャン方向に1次元フーリエ変換を行う。すると、0次、1次、−1次のピークが出る。ここで、0次のピークをそれぞれ抽出し、そのデータのみを用いて逆フーリエ変換すると、H0、V0画像が得られる。同様に、1次のピークをそれぞれ抽出し、そのデータのみを用いて逆フーリエ変換すると、H1、V1画像が得られる。 Two-dimensional AB scan image data obtained during one B scan in this way is subjected to a one-dimensional Fourier transform in the B scan direction. Then, 0th, 1st and −1st order peaks appear. Here, when 0th-order peaks are extracted and inverse Fourier transform is performed using only the data, H0 and V0 images are obtained. Similarly, H1 and V1 images are obtained by extracting primary peaks and performing inverse Fourier transform using only the data.
H0、H1画像から、試料17の偏光特性である式(後記する式(18)参照)のジョーンズマトリックスの成分のうち、J(1,1)およびJ(1,2)を求める事ができる。そして、V0、V1画像から、試料17の偏光特性である式(1)のジョーンズマトリックスの成分のうち、J(2,1)およびJ(2,2)を求める事ができる。 From the H0 and H1 images, J (1,1) and J (1,2) can be obtained from the Jones matrix components of the equation (see equation (18) below) which is the polarization characteristic of the sample 17. From the V0 and V1 images, J (2,1) and J (2,2) among the components of the Jones matrix of the expression (1), which is the polarization characteristic of the sample 17, can be obtained.
このようにして、1回のBスキャンにおついて4つの偏光特性を含む情報が得られる。そして、この4つの情報をそれぞれ、通常のFD−OCTと同様にAスキャン方向にフーリエ変換すると、1次のピークが試料17の深さ方向の情報を有し、しかもそれぞれ偏光特性に応じた4枚のA−B画像が得られる。 In this way, information including four polarization characteristics can be obtained for one B-scan. Then, when these four pieces of information are Fourier-transformed in the A-scan direction in the same manner as in ordinary FD-OCT, the primary peak has information in the depth direction of the sample 17, and each 4 corresponds to the polarization characteristic. A sheet of AB image is obtained.
ところで、偏光感受型光コヒーレンストモグラフィー装置では、入射光の偏光状態は 横方向のBスキャン中に偏光情報を埋め込むために直線偏光(垂直偏光)→楕円偏光→直線偏光(水平偏光)などのように連続的に、あるいは、3ステップ変調される。 By the way, in the polarization-sensitive optical coherence tomography apparatus, the polarization state of incident light is such as linearly polarized light (vertical polarized light) → elliptical polarized light → linearly polarized light (horizontal polarized light) in order to embed the polarization information during the horizontal B scan. Continuously or three-step modulated.
2つのCCDカメラ22、23によって検出されたそれぞれの光強度信号をフーリエ変換することにより、位相情報を含んだOCT信号を得る。このOCT信号の0次と1次の周波数成分を取り出し、それぞれそれを逆のフーリエ変換する。これらの値を使って、試料の偏光特性を表すジョーンズ行列のすべての要素が得られる。最終的に、試料の物理的な偏光依存性を表す値であるところの、位相遅延量(複屈折量)、複屈折の相対的な方位分布(方向)と二色性が計算される。 OCT signals including phase information are obtained by Fourier transforming the respective light intensity signals detected by the two CCD cameras 22 and 23. The 0th-order and 1st-order frequency components of this OCT signal are taken out and inverse Fourier transformed respectively. Using these values, all the elements of the Jones matrix representing the polarization properties of the sample are obtained. Finally, the phase retardation amount (birefringence amount), the relative orientation distribution (direction) of birefringence, and dichroism, which are values representing the physical polarization dependence of the sample, are calculated.
(測定信号のノイズ処理方法)
偏光感受型光コヒーレンストモグラフィーは試料の位相遅延量(位相情報)は、一般に、位相角度で表されその値(信号)は、たとえば、0〜180度のような正の値で計算される。このように分布する信号の場合、0度あるいは180度付近の微小信号領域では、それに重畳されるノイズの分布が対称(たとえばガウス分布)であっても測定信号は対象とならない。これは、たとえば偏光感受型光コヒーレンストモグラフィーで得られる位相情報のような観測信号が正の限られた値(0度より大180度より小)しか取り得ないからである。
(Measurement signal noise processing method)
In polarization-sensitive optical coherence tomography, the amount of phase delay (phase information) of a sample is generally represented by a phase angle, and the value (signal) is calculated as a positive value such as 0 to 180 degrees. In the case of a signal distributed in this way, in a minute signal region near 0 degrees or 180 degrees, the measurement signal is not a target even if the noise distribution superimposed on it is symmetric (for example, Gaussian distribution). This is because, for example, an observation signal such as phase information obtained by polarization-sensitive optical coherence tomography can take only a positive limited value (greater than 0 degrees and smaller than 180 degrees).
図2及び図3は、位相に対する信号強度の一般的な測定における、真の値(図2(a)、3(a)参照)、ノイズ分布(図2(b)、図3(b)参照)、観測値の関係(図2(c)、図3(c)参照)を説明する図(グラフ)である。図2(a)に示すように真の値が0度あるいは180度近傍の微小信号ではない値の場合、未知のランダムノイズ(図2(b)参照)が真の信号の値に対して加わる場合、観測される信号は真の信号に対して対称に広がる(図2(c)参照)。従って、単純な移動平均でランダムノイズは平均化されキャンセルされるため、信号の真の値を推定することができる。 2 and 3 show true values (see FIGS. 2 (a) and 3 (a)) and noise distribution (see FIGS. 2 (b) and 3 (b)) in a general measurement of signal intensity with respect to phase. ) And the relationship between the observed values (see FIGS. 2C and 3C) (graph). As shown in FIG. 2A, when the true value is not a minute signal near 0 degrees or 180 degrees, unknown random noise (see FIG. 2B) is added to the true signal value. In this case, the observed signal spreads symmetrically with respect to the true signal (see FIG. 2C). Therefore, since the random noise is averaged and canceled by a simple moving average, the true value of the signal can be estimated.
しかしながら、図3に示すように、真の値が0度あるいは180度近傍の微小信号の場合(図3(a)参照)は、ランダムノイズ(図3(b)参照)が重畳された観測信号の分布(図3(c)参照)は、真の信号に対して対称でなくなるため、このような移動平均では真の値を推定することができない。 However, as shown in FIG. 3, in the case of a minute signal whose true value is 0 degree or in the vicinity of 180 degrees (see FIG. 3A), an observation signal on which random noise (see FIG. 3B) is superimposed. Since the distribution (see FIG. 3C) is not symmetric with respect to the true signal, the true value cannot be estimated with such a moving average.
即ち、通常の測定では、データを繰り返し計測し、その分布を求め平均を持って真の値とする。しかし、ゼロ付近ではデータの値そのものが小さいので、元々ある測定ノイズなどに埋もれてしまい、ゼロ点では計測データの検出はできないので、分布がゼロになる。また、データの分布は経験的に実際はなめらかになるので、真の値の周りに非対称に広がってしまい、単に移動平均では真の値を推定することができない。 That is, in normal measurement, data is repeatedly measured, its distribution is obtained, and an average is obtained as a true value. However, since the data value itself is small near zero, it is buried in the original measurement noise and the measurement data cannot be detected at the zero point, so the distribution becomes zero. In addition, since the data distribution is actually empirically smooth, it spreads asymmetrically around the true value, and the true value cannot be simply estimated by the moving average.
本発明の測定信号のノイズ処理方法を適用することで、偏光感受型光コヒーレンストモグラフィー装置では、こうした未知のランダムノイズが重畳され対称でない観測信号の真の値を推定するために、図3(c)に示すように、観測信号の最大値(モード値)をとる構成とした。これにより、微小な(観測)信号領域での画像の解像度を上げることができる。 By applying the measurement signal noise processing method of the present invention, in the polarization-sensitive optical coherence tomography apparatus, in order to estimate the true value of the observation signal superimposed with such unknown random noise, FIG. ), The maximum value (mode value) of the observation signal is taken. Thereby, the resolution of the image in a minute (observation) signal region can be increased.
このような微小信号領域での微小測定信号のノイズ処理方法について、さらに以下説明する。本実施例において、偏光感受型光コヒーレンストモグラフィー装置で得られた画像信号から信号を位相遅延量を取得する信号処理装置(具体的には、コンピュータ)では、次に説明するような位相遅延量の処理を行うことで、前記「位相遅延量は実際の値より高く出ている。」という不適切さを回避することができる構成にある。 The noise processing method of the minute measurement signal in such a minute signal region will be further described below. In this embodiment, a signal processing device (specifically, a computer) that acquires a phase delay amount from an image signal obtained by a polarization-sensitive optical coherence tomography device has a phase delay amount as described below. By performing the processing, it is possible to avoid the inappropriateness that “the phase delay amount is higher than the actual value”.
ところで、偏光感受型光コヒーレンストモグラフィーは試料の位相遅延量を測定することができる。3次元の測定対象の特定の位置の位相遅延量は、2つのパラメータを持っており、1つは「位相遅延量そのものの大きさ」、もう一つは、「位相遅延量の軸の方向」である。 By the way, polarization sensitive optical coherence tomography can measure the amount of phase delay of a sample. The phase delay amount at a specific position of the three-dimensional measurement object has two parameters, one is “the size of the phase delay amount itself” and the other is “the direction of the axis of the phase delay amount”. It is.
「位相遅延量そのものの大きさ」は、特定の方向の位相遅延量、あるいは、位相遅延量の最大値、最小値などである。位相遅延量は、1周期を360度(2π)として、づれ量を角度の単位で表し、180度で折り返す性質があるので、範囲は0〜180度となる。 “The magnitude of the phase delay amount itself” is a phase delay amount in a specific direction, or a maximum value or a minimum value of the phase delay amount. The phase delay amount is set to 360 degrees (2π) for one period, and the amount of deviation is expressed in units of angles and has a property of turning back at 180 degrees, so the range is 0 to 180 degrees.
「位相遅延量の軸の方向」は、遅延量が、最大あるいは最小となる軸の方向であり、0〜180度の角度で表される。この軸を「速軸」あるいは「遅軸」という。 “Axis direction of phase delay amount” is an axis direction in which the delay amount is maximized or minimized, and is represented by an angle of 0 to 180 degrees. This axis is called “fast axis” or “slow axis”.
以下、偏光感受型光コヒーレンストモグラフィーは試料の位相延量を測定では、「位相遅延量そのものの大きさ」と「位相遅延量の軸の方向」の2つの測定測定信号のノイズ処理方法について、それぞれ、上記同様に図2及び図3を参照し、さらに図4を参照して説明する。 In the following, the polarization-sensitive optical coherence tomography measures the phase extension of the sample. For the noise processing methods of the two measurement measurement signals, “the size of the phase delay itself” and “the direction of the axis of the phase delay”, respectively. The above will be described with reference to FIGS. 2 and 3 as well as FIG.
(位相遅延量そのものの大きの測定信号のノイズ処理方法)
位相遅延量そのもの大きさは角度として計算されるため、0〜180度などの正の値をとる。一方、未知のノイズは一般的に対象な分布を持っている(図2(b)参照)。観測値は真の値にノイズが重畳したものである(図2(c)参照)。
(Noise processing method for measurement signals with large phase delay)
Since the phase delay amount itself is calculated as an angle, it takes a positive value such as 0 to 180 degrees. On the other hand, unknown noise generally has a target distribution (see FIG. 2B). The observed value is obtained by superimposing noise on the true value (see FIG. 2C).
真の値が0または180度から離れている時(図2(a)参照)は、観測値はノイズ(図2(b)参照)のため真の値の周りに対照的に広がって測定される(図2(c)参照)。そのため、観測値の単純移動平均をとることにより真の値を得ることができる(図2(c)参照)。 When the true value is away from 0 or 180 degrees (see Fig. 2 (a)), the observed value is measured in contrast to the true value due to noise (see Fig. 2 (b)). (See FIG. 2 (c)). Therefore, a true value can be obtained by taking a simple moving average of the observed values (see FIG. 2C).
真の値が小さく、0度付近(あるいは180度付近)の場合(図3(a)参照)は、ノイズ(図3(b)参照)により観測値は非対称に広がる(図3(c)参照)。この場合、単純移動平均では真の値より高くなる。この場合、観測値のピーク(最頻値、モード値)の位置をもって真の値とするのがよい(図3(c)参照)。 When the true value is small and near 0 degrees (or around 180 degrees) (see FIG. 3A), the observed value spreads asymmetrically due to noise (see FIG. 3B) (see FIG. 3C). ). In this case, the simple moving average is higher than the true value. In this case, the true value is preferably determined from the position of the peak (mode, mode value) of the observed value (see FIG. 3C).
具体的には、適当な離散化を行い、観測値をヒストグラム化する。もっとも高いヒストグラムの中心値などを真の値とする(ヒストグラム法)(図4参照)。 Specifically, appropriate discretization is performed, and the observed values are converted into a histogram. The center value of the highest histogram is set as a true value (histogram method) (see FIG. 4).
なお、上記信号処理方法は、OCT信号の処理に限るものではなく、一般に正の値しか取り得ないデータについて、そのデータの値が小さい場合の信号処理方法として適用されるものである。 The signal processing method is not limited to the processing of the OCT signal, and is generally applied as a signal processing method when the data value is small for data that can only take a positive value.
(位相遅延量の軸の方向の測定信号のノイズ処理方法)
偏光感受型光コヒーレンストモグラフィーは、前述のとおり、試料の位相遅延量の最大(あるいは最小)の方向(速軸あるいは遅軸方向)を測定することができる。速軸方向は一般に角度で表され、その値は、0〜180度のような正の値で計算される。
(Measurement signal noise processing method in the direction of the axis of phase delay)
As described above, the polarization-sensitive optical coherence tomography can measure the maximum (or minimum) direction (fast axis or slow axis direction) of the phase delay amount of the sample. The fast axis direction is generally expressed as an angle, and the value is calculated as a positive value such as 0 to 180 degrees.
このように分布する信号の場合、0度あるいは180度付近の信号領域では、それに重畳される未知のノイズの分布が対称(たとえばガウス分布)であっても測定信号は対称とはならない。これは、たとえば偏光感受型光コヒーレンストモグラフィーで得られる速軸方向のような角度情報が正の限られた値(0度より大180度より小)しか取り得ないからである。 In the case of a signal distributed in this way, in the signal region near 0 degree or 180 degrees, the measurement signal is not symmetric even if the distribution of unknown noise superimposed thereon is symmetric (for example, Gaussian distribution). This is because, for example, angle information such as the fast axis direction obtained by polarization-sensitive optical coherence tomography can take only a limited positive value (greater than 0 degrees and smaller than 180 degrees).
このように、偏光感受型光コヒーレンストモグラフィー装置における速軸方向(角度)の測定において、上記のとおり未知のランダムノイズが重畳されて生じた対称でない信号の真の値を推定する場合について、本発明の測定信号のノイズ処理方法を適用した例を、上記位相遅延量自体の測定信号のノイズ処理方法と同様に、図2〜4を参照して、以下説明する。 As described above, in the measurement in the fast axis direction (angle) in the polarization-sensitive optical coherence tomography apparatus, the present invention relates to the case of estimating the true value of a non-symmetric signal generated by superimposing unknown random noise as described above. An example of applying the measurement signal noise processing method will be described below with reference to FIGS. 2 to 4, similarly to the measurement signal noise processing method of the phase delay amount itself.
本実施例において、偏光感受型光コヒーレンストモグラフィー装置で得られた画像信号から信号を例えば速軸方向(角度)を取得する信号処理装置(具体的には、コンピュータ)では、次に説明するような処理を行うことで、前記「速軸方向角度は実際の値より高く出ている。」という不適切さを回避することができる構成にある。 In this embodiment, a signal processing apparatus (specifically, a computer) that acquires a signal, for example, a fast axis direction (angle) from an image signal obtained by a polarization-sensitive optical coherence tomography apparatus will be described below. By performing the processing, it is possible to avoid the inappropriateness that “the angle in the fast axis direction is higher than the actual value”.
例えば、速軸方向は角度として計算されるため、0〜180度などの正の値をとる。一方、未知のノイズは一般的に対象な分布を持っている(図2(b)参照)。観測値は真の値にノイズが重畳したものである(図2(c)参照)。 For example, since the fast axis direction is calculated as an angle, it takes a positive value such as 0 to 180 degrees. On the other hand, unknown noise generally has a target distribution (see FIG. 2B). The observed value is obtained by superimposing noise on the true value (see FIG. 2C).
真の値が0または180度から離れている時(図2(a)参照)は、観測値はノイズ(図2(b)参照)のため真の値の周りに対照的に広がって測定される(図2(c)参照)。そのため、観測値の単純移動平均をとることにより真の値を得ることができる(図2(c)参照)。 When the true value is away from 0 or 180 degrees (see Fig. 2 (a)), the observed value is measured in contrast to the true value due to noise (see Fig. 2 (b)). (See FIG. 2 (c)). Therefore, a true value can be obtained by taking a simple moving average of the observed values (see FIG. 2C).
真の値が信号領域の端に近く、0度付近(あるいは180度付近)の場合(図3(a)参照)は、ノイズ(図3(b)参照)により観測値は非対称に広がる(図3(c)参照)。この場合、単純移動平均では真の値より高くなる。この場合、観測値のピーク(最頻値、モード値)の位置をもって真の値とするのがよい(図3(c)参照)。この図3(c)に示すように、信号の最大値(モード値)をとる構成とすることで、信号領域の限られた信号領域での画像の解像度を上げることができる。 When the true value is close to the edge of the signal region and near 0 degrees (or near 180 degrees) (see FIG. 3A), the observed value spreads asymmetrically due to noise (see FIG. 3B) (see FIG. 3B). 3 (c)). In this case, the simple moving average is higher than the true value. In this case, the true value is preferably determined from the position of the peak (mode, mode value) of the observed value (see FIG. 3C). As shown in FIG. 3C, by adopting a configuration that takes the maximum value (mode value) of the signal, the resolution of the image in the signal region where the signal region is limited can be increased.
具体的には、適当な離散化を行い、観測値をヒストグラム化する。もっとも高いヒストグラムの中心値などを真の値とする(ヒストグラム法)(図4参照)。 Specifically, appropriate discretization is performed, and the observed values are histogrammed. The center value of the highest histogram is set as a true value (histogram method) (see FIG. 4).
なお、上記信号処理方法は、OCT信号の処理に限るものではなく、一般に限られた信号領域しか取り柄ないデータについて、そのデータの値が小さい場合の信号処理方法として適用されるものである。 The signal processing method is not limited to the OCT signal processing, and is generally applied as a signal processing method when the data value is small for data dealing with only a limited signal region.
また、上記信号処理方法は、信号領域が限られた微小信号の処理に限るものではなく、一般にノイズの分布が非対称になるようなデータについて、そのデータの値が小さい場合の信号処理方法として適用されるものである。 The signal processing method is not limited to the processing of a minute signal with a limited signal area, and is generally applied as a signal processing method for data whose noise distribution is asymmetric when the data value is small. It is what is done.
以上、本発明に係る測定信号のノイズ処理方法を実施するための最良の形態を実施例に基づいて説明したが、本発明はこのような実施例に限定されることなく、特許請求の範囲記載の技術的事項の範囲内で、いろいろな実施例があることは言うまでもない。 The best mode for carrying out the measurement signal noise processing method according to the present invention has been described above based on the embodiments. However, the present invention is not limited to such embodiments, and the scope of the claims is described. It goes without saying that there are various embodiments within the scope of the technical matters.
本発明に係る測定信号のノイズ処理方法は、以上のような構成であるから、光断層画像取得用OCTだけでなく、各種の計測機器の測定信号のノイズ処理方法として適用可能である。 Since the measurement signal noise processing method according to the present invention has the above-described configuration, it can be applied not only to optical tomographic image acquisition OCT but also as a measurement signal noise processing method for various measuring instruments.
1 偏光感受光画像計測装置
2 光源
3、12 偏光子
4 EO変調器(偏光変調器、電気光学変調器)
5 ファイバーカプラー(光カプラー)
6 参照アーム
7 試料アーム
8 分光器
9 ファイバー
10、15、18 偏波コントローラ
11 コリメートレンズ
13 集光レンズ
14 参照鏡(固定鏡)
16 ガルバノ鏡
17 試料
19 回折格子
20 フーリエ変換レンズ
21 偏光ビームスプリッター
22、23 光検出器(ラインCCDカメラ)
24 固定鏡
1 Polarized light receiving image measuring device
2 Light source
3, 12 Polarizer
4 EO modulator (polarization modulator, electro-optic modulator)
5 Fiber coupler (optical coupler)
6 Reference arm
7 Sample arm
8 Spectrometer
9 Fiber
10, 15, 18 Polarization controller
11 Collimating lens
13 Condensing lens
14 Reference mirror (fixed mirror)
16 Galvano mirror
17 samples
19 Diffraction grating
20 Fourier transform lens
21 Polarizing beam splitter
22, 23 Photodetector (Line CCD camera)
24 fixed mirror
Claims (3)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007198155A JP5173305B2 (en) | 2007-07-30 | 2007-07-30 | Measurement signal noise processing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007198155A JP5173305B2 (en) | 2007-07-30 | 2007-07-30 | Measurement signal noise processing method |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2009031229A true JP2009031229A (en) | 2009-02-12 |
JP2009031229A5 JP2009031229A5 (en) | 2010-09-02 |
JP5173305B2 JP5173305B2 (en) | 2013-04-03 |
Family
ID=40401873
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007198155A Expired - Fee Related JP5173305B2 (en) | 2007-07-30 | 2007-07-30 | Measurement signal noise processing method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5173305B2 (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102042849A (en) * | 2010-12-23 | 2011-05-04 | 北京大学 | Truth value estimation method based on multi-dimensional vector information and device thereof |
WO2013008516A1 (en) * | 2011-07-12 | 2013-01-17 | 国立大学法人筑波大学 | Program to correct data measured by ps-oct and ps-oct system that comes with the program |
JP2014232103A (en) * | 2013-05-01 | 2014-12-11 | キヤノン株式会社 | Method and device for processing polarization data of polarization oct |
WO2015186726A1 (en) * | 2014-06-06 | 2015-12-10 | 国立大学法人筑波大学 | Polarization sensitive optical image measurement system, and program loaded into said system |
KR20210130046A (en) * | 2020-04-21 | 2021-10-29 | 주식회사 필로포스 | Compact oct spectrometer with improved crosstalk, echo signal, and nonlinearity |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002503134A (en) * | 1997-06-02 | 2002-01-29 | イザット,ジョーゼフ,エイ. | Imaging of Doppler flow using optical coherence tomography |
JP2004286515A (en) * | 2003-03-20 | 2004-10-14 | Olympus Corp | Method for obtaining existence ratio of a plurality of natural spectrum constituents from spectral image |
JP2007127425A (en) * | 2005-10-31 | 2007-05-24 | Univ Of Tsukuba | Correction method in optical tomography |
-
2007
- 2007-07-30 JP JP2007198155A patent/JP5173305B2/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002503134A (en) * | 1997-06-02 | 2002-01-29 | イザット,ジョーゼフ,エイ. | Imaging of Doppler flow using optical coherence tomography |
JP2004286515A (en) * | 2003-03-20 | 2004-10-14 | Olympus Corp | Method for obtaining existence ratio of a plurality of natural spectrum constituents from spectral image |
JP2007127425A (en) * | 2005-10-31 | 2007-05-24 | Univ Of Tsukuba | Correction method in optical tomography |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102042849A (en) * | 2010-12-23 | 2011-05-04 | 北京大学 | Truth value estimation method based on multi-dimensional vector information and device thereof |
EP2733478A4 (en) * | 2011-07-12 | 2014-12-31 | Univ Tsukuba | PROGRAM TO CORRECT PS-OCT MEASURED DATA AND PS-OCT SYSTEM THAT ACCOMPANIES THE PROGRAM |
KR101593429B1 (en) | 2011-07-12 | 2016-02-12 | 고쿠리쯔 다이가쿠 호징 츠쿠바 다이가쿠 | A program for correcting the measurement data of the PSOCT and a PSOCT system |
KR20140006057A (en) * | 2011-07-12 | 2014-01-15 | 고쿠리쯔 다이가쿠 호징 츠쿠바 다이가쿠 | Program to correct data measured by ps-oct and ps-oct system that comes with the program |
CN103814286A (en) * | 2011-07-12 | 2014-05-21 | 国立大学法人筑波大学 | Program to correct data measured by PS-OCT and PS-OCT system that comes with the program |
WO2013008516A1 (en) * | 2011-07-12 | 2013-01-17 | 国立大学法人筑波大学 | Program to correct data measured by ps-oct and ps-oct system that comes with the program |
US9471277B2 (en) | 2011-07-12 | 2016-10-18 | University Of Tsukuba | Program for correcting data measured by PS-OCT and PS-OCT system equipped with the program |
JP2013019773A (en) * | 2011-07-12 | 2013-01-31 | Univ Of Tsukuba | Program to correct data measured by ps-oct and ps-oct system comprising the program |
JP2014232103A (en) * | 2013-05-01 | 2014-12-11 | キヤノン株式会社 | Method and device for processing polarization data of polarization oct |
JP2015230297A (en) * | 2014-06-06 | 2015-12-21 | 国立大学法人 筑波大学 | Polarization perceptive type optical image measurement system and program installed therein |
WO2015186726A1 (en) * | 2014-06-06 | 2015-12-10 | 国立大学法人筑波大学 | Polarization sensitive optical image measurement system, and program loaded into said system |
US9863869B2 (en) | 2014-06-06 | 2018-01-09 | University Of Tsukuba | Polarization-sensitive optical image measuring system and program installed in said system |
EP3153843A4 (en) * | 2014-06-06 | 2018-02-14 | University of Tsukuba | Polarization sensitive optical image measurement system, and program loaded into said system |
KR20210130046A (en) * | 2020-04-21 | 2021-10-29 | 주식회사 필로포스 | Compact oct spectrometer with improved crosstalk, echo signal, and nonlinearity |
KR102339811B1 (en) | 2020-04-21 | 2021-12-16 | 주식회사 필로포스 | Compact oct spectrometer with improved crosstalk, echo signal, and nonlinearity |
Also Published As
Publication number | Publication date |
---|---|
JP5173305B2 (en) | 2013-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4344829B2 (en) | Polarized light receiving image measuring device | |
US8879070B2 (en) | Two beams formed by Wollaston prism in sample arm in an optical coherence tomography apparatus | |
JP5787255B2 (en) | Program for correcting measurement data of PS-OCT and PS-OCT system equipped with the program | |
CN103344569B (en) | Polarization complex frequency domain optical coherence tomography method and system | |
CN104523239B (en) | Full-depth spectral domain optical coherent tomography device and method | |
CN110779464B (en) | A time-domain and frequency-domain joint analysis wide-spectrum coherence measurement method and system | |
JP4378533B2 (en) | Calibration method for components of optical coherence tomography | |
US7903256B2 (en) | Methods, systems, and computer program products for performing real-time quadrature projection based Fourier domain optical coherence tomography | |
JP6256879B2 (en) | Polarization-sensitive optical image measurement system and program installed in the system | |
JP2004028970A (en) | Polarization-sensitive optical spectrum interference coherence tomography apparatus and method for measuring polarization information inside sample using the apparatus | |
JP2009008393A (en) | Optical image measuring device | |
JP2009025245A (en) | Device for observing optical interference | |
CN103792192B (en) | Polarization Frequency Domain Optical Coherence Tomography System Based on Single Detector | |
JP5173305B2 (en) | Measurement signal noise processing method | |
JP4461258B2 (en) | Correction method in optical tomography | |
JP2010151684A (en) | Polarization sensitive optical image measuring instrument for extracting local double refraction information | |
JP2001272332A (en) | Angular dispersion optical spatial coherence tomographic imaging system | |
JP2019537726A (en) | Method and apparatus for optimizing optical performance of an interferometer | |
JP3934131B2 (en) | Coaxial spatial optical coherence tomographic image measurement device | |
KR100868439B1 (en) | Interference System for Polarization-sensitive Photocoherence Biological Imaging Equipment | |
JP2006250849A (en) | Optical image measurement method and apparatus using optical coherence tomography apparatus | |
JP2009031230A (en) | Displaying measurement data | |
JP2010151685A (en) | Device for analyzing, texture-discriminating and visualizing image information in optical image measuring instrument | |
JP2009258131A (en) | Method for calibrating constituent device of optical coherence tomography | |
Suzuki et al. | 3D thickness measurement using pulse-driven optical coherence tomography based on wavelet transform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100720 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20100720 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20120328 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120403 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120528 |
|
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: 20121225 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20121227 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5173305 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 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |