JP7218122B2 - 眼科撮影装置及びその制御方法 - Google Patents
眼科撮影装置及びその制御方法 Download PDFInfo
- Publication number
- JP7218122B2 JP7218122B2 JP2018147777A JP2018147777A JP7218122B2 JP 7218122 B2 JP7218122 B2 JP 7218122B2 JP 2018147777 A JP2018147777 A JP 2018147777A JP 2018147777 A JP2018147777 A JP 2018147777A JP 7218122 B2 JP7218122 B2 JP 7218122B2
- Authority
- JP
- Japan
- Prior art keywords
- optical path
- components
- light
- unit
- path length
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims description 21
- 238000000034 method Methods 0.000 title claims description 18
- 230000003287 optical effect Effects 0.000 claims description 51
- 238000012545 processing Methods 0.000 claims description 32
- 238000006243 chemical reaction Methods 0.000 claims description 17
- 238000005259 measurement Methods 0.000 claims description 17
- 230000001131 transforming effect Effects 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000002238 attenuated effect Effects 0.000 claims description 4
- 230000002194 synthesizing effect Effects 0.000 claims description 2
- 238000012014 optical coherence tomography Methods 0.000 description 12
- 238000002834 transmittance Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 7
- 238000010408 sweeping Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 210000001525 retina Anatomy 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000034162 Primary hypomagnesemia with secondary hypocalcemia Diseases 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000004424 eye movement Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 201000002062 intestinal hypomagnesemia 1 Diseases 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/12—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
- A61B3/1225—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes using coherent radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/102—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for optical coherence tomography [OCT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/0016—Operational features thereof
- A61B3/0025—Operational features thereof characterised by electronic signal processing, e.g. eye models
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02001—Interferometers characterised by controlling or generating intrinsic radiation properties
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02001—Interferometers characterised by controlling or generating intrinsic radiation properties
- G01B9/02002—Interferometers characterised by controlling or generating intrinsic radiation properties using two or more frequencies
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02015—Interferometers characterised by the beam path configuration
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/0209—Low-coherence interferometers
- G01B9/02091—Tomographic interferometers, e.g. based on optical coherence
-
- G—PHYSICS
- G02—OPTICS
- G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
- G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
-
- G—PHYSICS
- G02—OPTICS
- G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
- G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
- G02F1/0009—Materials therefor
- G02F1/0018—Electro-optical materials
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10101—Optical tomography; Optical coherence tomography [OCT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30041—Eye; Retina; Ophthalmic
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/456—Optical coherence tomography [OCT]
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Ophthalmology & Optometry (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Theoretical Computer Science (AREA)
- Nonlinear Science (AREA)
- Signal Processing (AREA)
- Optics & Photonics (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Eye Examination Apparatus (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Description
本発明は、眼科撮影装置及びその制御方法に関する。
現在、眼を観察や撮影等するための眼科装置としては、例えば、前眼部撮影装置、眼底カメラ、共焦点レーザー走査検眼鏡（Ｓｃａｎｎｉｎｇ Ｌａｓｅｒ Ｏｐｈｔｈａｌｍｏｓｃｏｐｅ：ＳＬＯ）装置、多波長光波干渉を利用した光干渉断層撮影（Ｏｐｔｉｃａｌ Ｃｏｈｅｒｅｎｃｅ Ｔｏｍｏｇｒａｐｈｙ：ＯＣＴ）装置等がある。中でもＯＣＴ装置は被写体の断層画像を高解像度に得ることができるため、特に眼科装置として網膜の専門外来では必要不可欠な装置になりつつある。
眼科用ＯＣＴ装置では、被検眼の眼底を広範囲で撮影する場合や、眼底の湾曲が大きい場合、あるいは撮影したい部位が深さ方向に長い場合に、参照光の光路長と照射光の光路長の差が大きくなることがある。ここで、照射光の光路長とは照射光の光路とその反射光の光路を合わせた光路の光路長を意味する。この光路長差が大きいとＯＣＴ装置では断層像の折り返しや欠けが発生し、眼底の形状を正しく表した断層像を得ることができない。
そこで、折り返しが発生した断層像を反転させて元の断層像と合成し、画像解析により網膜を検出する技術が特許文献１に開示されている。この技術により深さ方向に広い範囲で眼底の適切な形状解析を行うことができる。
ここで、従来の技術では、折り返しが発生した後の断層像を反転して接続し形状解析を行う技術である。このため、断層像から折り返しそのものを除去するものではない。すなわち、折り返しが発生した部位の詳細な構造の観察に適した断層像を得ることができない。
本発明の目的の一つは、深さ方向により広い断層像であって、折り返しが低減された断層像を取得することである。
本発明に係る眼科撮影装置の一つは、
測定光を照射した被検体からの戻り光と前記測定光に対応する参照光とによる干渉光を干渉信号として検出する検出部と、
前記検出された干渉信号をアナログ信号からディジタル信号に変換する変換部と、
前記変換して得た複数の信号に関する複数の成分を用いて、前記被検体の断層像を生成する演算処理部と、を有し、
前記演算処理部は、前記変換部のナイキスト周波数より高い周波数の成分と、前記変換部のナイキスト周波数より低い周波数の成分とを含む前記複数の成分を用いて、前記断層像を生成する。
測定光を照射した被検体からの戻り光と前記測定光に対応する参照光とによる干渉光を干渉信号として検出する検出部と、
前記検出された干渉信号をアナログ信号からディジタル信号に変換する変換部と、
前記変換して得た複数の信号に関する複数の成分を用いて、前記被検体の断層像を生成する演算処理部と、を有し、
前記演算処理部は、前記変換部のナイキスト周波数より高い周波数の成分と、前記変換部のナイキスト周波数より低い周波数の成分とを含む前記複数の成分を用いて、前記断層像を生成する。
本発明の一つによれば、深さ方向により広い断層像であって、折り返しが低減された断層像を取得することができる。
本実施形態に係る眼科撮影装置の一つは、変換部により変換された干渉信号から得られる複数の成分（例えば、強度と位相とを含む情報）であって、変換部のナイキスト周波数より高い周波数の成分と、変換部のナイキスト周波数より低い周波数の成分とを用いて、被検体の断層像を生成する。これにより、深さ方向により広い断層像であって、折り返しが低減された断層像を取得することができる。以下では本実施形態について説明する。
（光干渉断層計）
図１は、本実施形態に係るＯＣＴ（眼科撮影装置の一例）の構成を示す図である。波長掃引光源１０１は、出射する光の波長が短波長から長波長へ、または長波長から短波長へ時間的に変化する光源である。波長掃引光源１０１から出射された光Ｌが、分割手段の一例である分岐ファイバカップラ１０２で照射光ＬＡ（測定光）と参照光ＬＢに分波される。また、照射光ＬＡはコリメータレンズ１０３で平行光となり、走査ミラー１０８で反射され、集光レンズ１０５を経て、被検体１０７に照射される。また、被検体１０７に照射されて反射した光ＬＡ’（戻り光）は、カップラ１０２を経由して、カップラ１０４に入射する。一方、参照光ＬＢは、コリメータレンズ１０９、ミラー１１０、１１１、集光レンズ１１２を経て光ＬＢ’としてカップラ１０４に入射する。光ＬＡ’と光ＬＢ’とはカップラ１０４で合波すると同時に分波され、検出部の一例である差動検出器１２０に干渉光として入射する。
図１は、本実施形態に係るＯＣＴ（眼科撮影装置の一例）の構成を示す図である。波長掃引光源１０１は、出射する光の波長が短波長から長波長へ、または長波長から短波長へ時間的に変化する光源である。波長掃引光源１０１から出射された光Ｌが、分割手段の一例である分岐ファイバカップラ１０２で照射光ＬＡ（測定光）と参照光ＬＢに分波される。また、照射光ＬＡはコリメータレンズ１０３で平行光となり、走査ミラー１０８で反射され、集光レンズ１０５を経て、被検体１０７に照射される。また、被検体１０７に照射されて反射した光ＬＡ’（戻り光）は、カップラ１０２を経由して、カップラ１０４に入射する。一方、参照光ＬＢは、コリメータレンズ１０９、ミラー１１０、１１１、集光レンズ１１２を経て光ＬＢ’としてカップラ１０４に入射する。光ＬＡ’と光ＬＢ’とはカップラ１０４で合波すると同時に分波され、検出部の一例である差動検出器１２０に干渉光として入射する。
また、差動検出器１２０は、干渉光の強度をアナログ信号に変換する。差動検出器１２０から出力されたアナログ信号は分岐され、一方はローパスフィルタ１２１によって低周波成分のみとなり、変換部の一例であるＡＤ変換部１３１によりディジタル信号に変換される。他方はハイパスフィルタ１２２によって高周波成分のみとなり、別の変換部の一例であるＡＤ変換部１３２によってディジタル信号に変換される。ローパスフィルタ１２１とハイパスフィルタ１２２の帯域は電気的に調整可能である。なお、ローパスフィルタ１２１とハイパスフィルタ１２２は、検出された干渉信号を所定の周波数特性に従って減衰させるフィルタ部の一例である。これらのディジタル信号は、演算処理部の一例である演算処理装置１１４によって後述する方法で処理される。演算処理装置１１４は、これらのディジタル信号に対してフーリエ変換を始めとする処理を行うことによって被検体１０７の情報を取得する。表示部１１５は演算処理装置１１４によって得られた被検体の断層像を表示する。なお、演算処理装置は、眼科撮影装置に通信可能に接続されていれば良い。また、演算処理装置は、眼科撮影装置の内部に組み込まれていても良い。
以上は、被検体１０７のある１点における断層像の取得プロセスであり、このように被検体１０７の奥行き方向の断層に関する情報を取得することをＡ－ｓｃａｎと呼ぶ。また、Ａ－ｓｃａｎと直交する方向で被検体の断層に関する情報、すなわち２次元画像を取得するための走査をＢ－ｓｃａｎと呼ぶ。本実施形態においては走査ミラー１０８によって、Ｂ－ｓｃａｎを行う。
（断層像の取得）
断層像の取得について図２を用いて説明する。ここでは、図１におけるローパスフィルタ１２１の帯域が充分広く調整されている場合について説明する。この場合、差動検出器１２０から出力されるアナログ信号は全てローパスフィルタ１２１を通過するため、ローパスフィルタ１２１は無視してよいものとする。
断層像の取得について図２を用いて説明する。ここでは、図１におけるローパスフィルタ１２１の帯域が充分広く調整されている場合について説明する。この場合、差動検出器１２０から出力されるアナログ信号は全てローパスフィルタ１２１を通過するため、ローパスフィルタ１２１は無視してよいものとする。
ここで、図２（ａ）は、被検体１０７の眼底のある断層の様子を表したものである。照射光が深さ位置２４０で反射した場合に照射光の光路長と参照光の光路長とが一致するように、後述の方法によってミラー１１０、１１１の位置が調整されている。
また、図２（ｂ）は、図２（ａ）の位置２０１をＡ－ｓｃａｎすることによって得られた信号を演算処理装置１１４がフーリエ変換し、周波数成分ごとに信号強度を示したものである。横軸は周波数を、縦軸は信号強度の対数を表している。直流成分２４１は周波数０を表している。照射光の光路長と参照光の光路長とが完全に一致すると、掃引された全ての波長で照射光と参照光が強め合うため、直流成分の信号が得られる。照射光の光路長と参照光の光路長が異なる場合はその差に応じて干渉が発生し、信号の周波数が変化する。
また、位置２０１では、被検体１０７が深さ位置２４０よりも深い位置にあるため光路長差が生じ、実像２０２が得られる。ナイキスト周波数２４２はＡＤ変換部１３１のサンプリング周波数の半分の周波数であり、ＡＤ変換部１３１によって取得できる最大の周波数である。ＡＤ変換部１３１でアナログ／ディジタル変換を行うとナイキスト周波数を軸とした折り返し（エイリアシング）が発生する。鏡像２０３はＡＤ変換部１３１でサンプリングしたときに実像２０２と区別がつかない周波数を表した像である。
また、図２（ｃ）は、演算処理装置１１４によって生成されたある段階の断層像を示している。画像の上端２４３は直流成分２４１を表し、画像の下端２４４はナイキスト周波数２４２を表している。直線２０４は図２（ｂ）の信号を画像化したものである。実像２０２が直流成分からナイキスト周波数までの範囲に全て含まれているため、直線２０４上では画像の上端から下端の間に断層画像が含まれている。
また、図２（ｄ）は、位置２１１の信号を図２（ｂ）と同様に示したものである。位置２１１では位置２０１よりも被検体１０７が深い位置にあるため、実像２１２は実像２０２よりも高い周波数を持つ。実像２１２はナイキスト周波数を一部越えているため、折り返した鏡像２１３と一部重なっている。そのため、直線２１４上では本来の断層画像と折り返した断層画像が一部重なった画像になっている。
また、図２（ｅ）は、位置２２１の信号を図２（ｂ）と同様に示したものである。位置２２１では位置２１１よりも被検体１０７がさらに深い位置にあるため、実像２２２はさらに高い周波数となり、ナイキスト周波数を越えている。そのため直線２２４上では鏡像２２３すなわち折り返した部分のみが画像化されている。このようにナイキスト周波数を越えた周波数の信号をサンプリングすることを一般にアンダーサンプリングという。
以上から分かる通り、被検体１０７の深さ方向の範囲に対してＡＤ変換部１３１のサンプリング周波数が充分ではない場合、断層像に折り返しが発生し、被検体１０７の構造を正しく表した断層像を取得することができない。
なお、前述したミラー１１０、１１１の位置調整は以下のように行う。撮影者は表示部１１５に表示された図２（ｃ）またはそれに類する画像を見ながら、不図示の入力手段によって光路長変更部の一例であるミラー１１０、１１１の位置を調整する。ミラー１１０、１１１の位置を適切に調整すると、図２（ｂ）のように画像の下端のみで折り返しが発生する断層像を得ることができる。なお、ミラー１１０、１１１の位置の調整は、例えば、ミラー１１０、１１１を共通のステージで光軸に沿って移動することで行うことにより、参照光の光路長を変更することができる。なお、光路長変更部は、参照光の光路長の変更だけでなく、測定光の光路長を変更するものであっても良いし、参照光の光路長と測定光の光路長との両方を変更するものであっても良い。すなわち、光長長変更部は、参照光の光路長と測定光の光路長との差を変更できるものであれば、何でも良い。
なお、ミラー１１０、１１１の位置が適切でないと、画像の上端で折り返しが発生することがある。これは、深さ位置２４０と被検体１０７が重なっている場合である。このとき信号の周波数成分は直流成分２４１の両側に分布し、周波数が負の領域が鏡像として正の領域に折り返すことで断層像の折り返しが発生する。従来のＯＣＴ装置では、断層像の折り返しという表現はこの折り返しを指す場合が多い。一方で本発明では、このような画像の上端２４３における折り返しはミラー１１０、１１１の位置調整によって回避することを前提とする。その結果として画像の下端２４４でナイキスト周波数による折り返しが発生するが、以下ではこのナイキスト周波数による折り返しを抑制し、被検体１０７の構造を深さ方向に広い範囲で画像化するために演算処理装置１１４が追加で行う処理について述べる。
これ以降の説明では、説明を簡単にするためフィルタの位相シフトはないものとする。同様に、眼底をＡ－ｓｃａｎすることによって得られた信号は全ての周波数で位相が揃っているものとする。このとき、この信号をフーリエ変換すると実数のみからなる周波数分布が得られるため、以降の説明では実数のみを用いる。
（折り返しの低減）
また、図３は、被検体１０７の構造を深さ方向に広い範囲で画像化するための追加の処理について示した図である。ここで、図３（ａ）は、フィルタの周波数特性を示した図である。特性３０１はローパスフィルタ１２１の周波数特性であり、特性３１１はハイパスフィルタ１２２の周波数特性である。横軸は周波数を、縦軸はフィルタ透過率の対数を表している。いずれもカットオフ周波数はナイキスト周波数２４２付近になるように調整されている。なお、フィルタの周波数特性に関する情報として、図３（ａ）のようなグラフやグラフを示すテーブル等を記憶部に記憶させておくことが好ましい。
また、図３は、被検体１０７の構造を深さ方向に広い範囲で画像化するための追加の処理について示した図である。ここで、図３（ａ）は、フィルタの周波数特性を示した図である。特性３０１はローパスフィルタ１２１の周波数特性であり、特性３１１はハイパスフィルタ１２２の周波数特性である。横軸は周波数を、縦軸はフィルタ透過率の対数を表している。いずれもカットオフ周波数はナイキスト周波数２４２付近になるように調整されている。なお、フィルタの周波数特性に関する情報として、図３（ａ）のようなグラフやグラフを示すテーブル等を記憶部に記憶させておくことが好ましい。
また、図３（ｂ）は、ローパスフィルタ１２１の周波数特性と、ＡＤ変換部１３１によって得られたディジタル信号から取得した断層像３１２を示したものである。なお、断層像３１２は、変換部のナイキスト周波数より低い周波数の成分を用いて生成される第１の部分断層像の一例である。ここで、信号ＳＬと信号ＳＨはそれぞれの周波数において差動検出器１２０から出力された信号を示している。信号ＳＬと信号ＳＨの周波数はナイキスト周波数２４２を軸に対称になっている。信号ＳＬと信号ＳＨはローパスフィルタ１２１によって減衰する。透過率ＰＬ，ＰＨはそれぞれ信号ＳＬと信号ＳＨの周波数におけるローパスフィルタ１２１の透過率を示した値である。これらは装置の調整の際に決まる値である。ローパスフィルタ１２１によって信号ＳＬと信号ＳＨが減衰すると、ＡＤ変換部１３１に入力される信号はそれぞれＰＬＳＬとＰＨＳＨである。これらはＡＤ変換部１３１でアナログ／ディジタル変換されると折り返しによって重なり加算された結果１つの信号ＳＰとなる。演算処理装置１１４はフーリエ変換によって得られたこの信号ＳＰを記憶し、同時に信号ＳＰの強度を輝度とした断層像３１２を生成する。一方で、信号ＳＰは上記の説明より、以下の式で求めたものと一致する。
また、図３（ｃ）は、ハイパスフィルタ１２２の周波数特性と、ＡＤ変換部１３２によって得られたディジタル信号から取得した断層像３２２を図３（ｂ）と同様に示したものである。なお、断層像３２２は、変換部のナイキスト周波数より高い周波数の成分を用いて生成される第２の部分断層像の一例である。図３（ｂ）と同様に、信号ＳＱは以下の式で求めたものと一致する。
これらの式を連立方程式として解くと、以下のように信号ＳＬと信号ＳＨを得ることができる。
したがって、演算処理装置１１４は、断層像３１２と断層像３２２の生成に用いた信号およびフィルタの透過率を用いて、信号ＳＬと信号ＳＨを得ることができる。すなわち、折り返しによって重なっていた、ナイキスト周波数未満の信号ＳＬとナイキスト周波数以上の信号ＳＨを上記の方法で分離することができる。この計算をナイキスト周波数２４２以外の全ての周波数について行うことで、それぞれの周波数の信号を得ることができる。ナイキスト周波数２４２ではＰＬとＰＨが一致するため上記の計算は行うことができないが、アナログ／ディジタル変換によって得られた信号をそのまま用いることで、信号ＳＬ（この場合信号ＳＨと同じ信号である）を求めることができる。
また、図３（ｄ）は、以上により得られた信号の強度を画像化して被検体１０７を深い範囲で画像化した断層像を示したものである。図２（ｃ）や図３（ｂ）（ｃ）と比較して深さ方向に広い範囲の断層像が得られていることが分かる。
以上の説明ではフィルタの周波数シフトはないものとし、眼底をＡ－ｓｃａｎすることによって得られた信号は全ての周波数で位相が揃っているものとした。しかし実際の周波数フィルタでは位相シフトが発生する。また、実際に眼底をＡ－ｓｃａｎすることによって得られた信号には様々な位相の信号が含まれている。したがって、これらを考慮した計算を行うことが望ましい。具体的には、以上の説明で用いた信号ＳＬ、ＳＨ、透過率ＰＬ、ＰＨ、ＱＬ、ＱＨ、輝度ＳＰ、ＳＱを全て実数から複素数に置き換えればよい。透過率ＰＬ、ＰＨ、ＱＬ、ＱＨはフィルタのゲイン特性と位相特性から複素数として得ることができる。信号ＳＰ、ＳＱは眼底をＡ－ｓｃａｎすることによって得られた信号をフーリエ変換することで得ることができる。信号ＳＬ、ＳＨは以上の計算から複素数として得られるため、その強度を画像の輝度とすればよい。なお、本実施形態における成分とは、強度だけでなく、位相を含む情報であり、また、実装だけでなく、複素数を含む情報である。
図４は以上の説明に基づいて断層像を生成し表示するまでの流れを示したものである。ステップＳ４０１において、演算処理装置１１４はＡＤ変換部１３１とＡＤ変換部１３２によって得られた干渉信号を取得する。ステップＳ４０２において、演算処理装置１１４はそれらの干渉信号のフーリエ変換を行い、信号ＳＰ、ＳＱを得る。ステップＳ４０３において、演算処理装置１１４は上記で説明した方法で信号をナイキスト周波数で分離する計算を行い、信号ＳＬ、ＳＨを得る。ステップＳ４０４において、演算処理装置１１４は信号ＳＬ、ＳＨの強度を求めて図３（ｄ）に示す断層像を生成する。ステップＳ４０５において、表示部１１５が断層像を表示する。
本実施形態では短時間の撮影が可能である。なぜなら、折り返し強度の異なる画像を複数取得するために、ミラー１１０、１１１を駆動して参照光の光路長を変更する必要がなく、図３（ｄ）の断層像を得るためにミラー１０８が行う走査は１回でよいためである。
また、ローパスフィルタ１２１とハイパスフィルタ１２２の透過率の値は環境依存などを考慮して装置の調整後に更新してもよい。例えば、検査前や検査中に照射光を装置内部の不図示のミラーで反射させて得られた信号からこれらの透過率を計算してもよい。また、得られた図３（ｄ）の断層像の輝度分布に基づいて透過率の値を推定し、図３（ｄ）の断層像を再生成してもよい。
また、信号ＳＬと信号ＳＨを求める式は、上記以外のものでもよい。例えば、フィルタが理想フィルタに近く、カットオフ周波数がナイキスト周波数に近い場合は断層像３１２には折り返し部分がほぼ発生せず、断層像３２２にはほぼ折り返し部分のみが発生する。したがって、断層像３２２を上下反転させて断層像３１２に単純に接続することで、所望の断層像に近い断層像を得ることができる。これは上記の式でＰＨとＱＬを０とし、ＰＬとＱＨを１とした場合の計算と同じである。
また、使用するフィルタはローパスフィルタとハイパスフィルタの組み合わせでなくてもよい。例えば、ローパスフィルタのみでカットオフ周波数が異なる組み合わせを用いてもよい。ローパスフィルタのみであっても上記のＰＬとＱＬの値が異なるかもしくはＰＨとＱＨの値が異なれば同様の計算を行うことができる。また、ナイキスト周波数の２倍よりもさらに高い周波数を考慮してバンドパスフィルタを用いてもよい。
また、フィルタは３つ以上でもよい。周波数帯域の異なる３つ以上のフィルタを組み合わせて同様の計算を行うことで、さらに深い範囲の断層像を取得可能にしてもよい。また同じ帯域のフィルタを複数用いて得られる信号強度を平均化して輝度の算出精度を高くしてもよい。
また、フィルタは１つでもよい。フィルタを通した信号と通さない信号から同様の計算を行ってもよい。その場合、フィルタを通さない信号では減衰がほぼないため、ＰＬとＰＨの値はほぼ１とすればよい。さらに唯一のフィルタが理想ローパスフィルタに近ければ、ＱＬを１、ＱＨを０とすればよい。この場合、上記のＳＨの計算式は、フィルタを通さない信号からフィルタを通した信号の単純な減算に対応する。すなわちこのような構成では、信号の減算処理によって深さ方向に広い範囲の断層像を得ることが出来る。
また、１つのフィルタで時間的に帯域を切り替えてもよい。帯域を切り替えながら複数の断層像を連続で取得することにより、これら複数の断層像からより深い範囲の断層像を生成することができる。この場合、被検眼の動きの影響を小さくするため、短時間で複数の断層像を取得することが好ましい。そのためには、本実施形態のようにフィルタの帯域は電気的に変更できる構成が好ましい。なお、フィルタの帯域を切り替える構成としては、例えば、フィルタに設けられ、フィルタの周波数特性を変更する特性変更部（例えば、可変抵抗）と、特性変更部を制御する制御部とから成る構成が好ましい。
また、本実施形態ではフィルタ以外による信号の減衰もフィルタの透過率の値に含めている。すなわち、デフォーカスによる信号強度の低下やＡＤ変換部内部のフィルタ特性による減衰もローパスフィルタ１２１およびハイパスフィルタ１２２による減衰と見なして計算を行っている。しかし、これらは別のパラメータとして扱ってもよい。
また、ＡＤ変換部１３１とＡＤ変換部１３２のサンプリング周波数はお互いに異なっていてもよい。演算処理装置１１４が上記の処理を行うかどうかを切り替えられるようにしてもよい。
また、断層像３１２と断層像３２２を生成せず、得られた信号から直接上記の式で最終的な断層像図３（ｄ）を生成してもよい。その場合、断層像３１２や３２２の画像化が不要なため、より高速な処理が可能である。
また、本実施形態では、波長掃引光源１０１による波長掃引には安定性があり、波数に対して一定速度である場合を想定しているが、一定速度ではない光源であっても良い。このような場合には、光源または装置内部で等波数間隔のサンプリングを行うためのｋ－ｃｌｏｃｋを生成し、ＡＤ変換部に入力するような構成が好ましい。なお、ｋ－ｃｌｏｃｋは、変換部がアナログ信号をサンプリングするクロックを生成するクロック生成部の一例である。また、ｋ－ｃｌｏｃｋは、波長掃引光源１０１からの光のうち一部の光が通る光路が第一光路と第一光路に対して光路長差を有する第二光路とに分岐された干渉計として構成されることが好ましい。ここで、波長掃引光源１０１による波長掃引が一定速度ではない場合であっても、変換部のナイキスト周波数より高い周波数の成分は、変換部による干渉信号（アナログ信号）のサンプリングが２回行われる時間より短い時間で振動する成分とみなせるものとする。また、波長掃引光源１０１による波長掃引が一定速度ではない場合であっても、変換部のナイキスト周波数より低い周波数の成分は、変換部による干渉信号（アナログ信号）のサンプリングが２回行われる時間より長い時間で振動する成分とみなせるものとする。
また、本実施形態は、ＳＳ－ＯＣＴであるが、その他のＯＣＴに適用してもよい。例えば、ＳＤ－ＯＣＴではラインセンサの前に光学ローパスフィルタを配置して本実施形態のローパスフィルタ１２１と同等の役割をさせてもよい。また、ラインセンサの画素サイズや回折格子、レンズ等の設計に基づく分解能の低下をローパスフィルタと同等の効果と考え、本発明を適用してもよい。
以上の説明のように、本実施形態によれば、ナイキスト周波数より高い周波数の信号と低い周波数の信号をそれぞれ計算することができる。それにより、ナイキスト周波数より高い周波数の領域すなわち眼底のより深い領域を画像化することができ、深さ方向に広い断層像を得ることが出来る。さらに、光路長を変更する方法と比較すると信号取得を短時間で行うことができ、被検眼の固視微動や疲労の影響を抑えることが出来る。
（その他の実施形態）
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。
Claims (12)
- 測定光を照射した被検体からの戻り光と前記測定光に対応する参照光とによる干渉光を干渉信号として検出する検出部と、
前記検出された干渉信号をアナログ信号からディジタル信号に変換する変換部と、
前記変換して得た複数の信号に関する複数の成分を用いて、前記被検体の断層像を生成する演算処理部と、を有し、
前記演算処理部は、前記変換部のナイキスト周波数より高い周波数の成分と、前記変換部のナイキスト周波数より低い周波数の成分とを含む前記複数の成分を用いて、前記断層像を生成することを特徴とする眼科撮影装置。 - 前記演算処理部は、前記ナイキスト周波数より低い周波数の成分を用いて第１の部分断層像を生成し、前記ナイキスト周波数より高い周波数の成分を用いて第２の部分断層像を生成し、前記第１の部分断層像と前記第２の部分断層像とを合成することによって前記断層像を生成することを特徴とする請求項１に記載の眼科撮影装置。
- 前記検出された干渉信号を所定の周波数特性に従って減衰させるフィルタ部をさらに有し、
前記変換部は、前記減衰した干渉信号を変換することを特徴とする請求項１または２に記載の眼科撮影装置。 - 前記変換部とは異なり、前記検出された干渉信号を変換する別の変換部をさらに有し、
前記演算処理部は、前記減衰した干渉信号と前記別の変換部によって変換された干渉信号とから得られる前記複数の成分を用いて、前記断層像を生成することを特徴とする請求項３に記載の眼科撮影装置。 - 前記フィルタ部に設けられ、前記フィルタ部の周波数特性を変更する特性変更部と、
前記特性変更部を制御する制御部と、をさらに有し、
前記演算処理部は、前記周波数特性が変更される前に変換された干渉信号と、前記周波数特性が変更された後に変換された干渉信号とから得られる前記複数の成分を用いて、前記断層像を生成することを特徴とする請求項３に記載の眼科撮影装置。 - 前記演算処理部は、前記フィルタ部の周波数特性に関する情報と前記変換された干渉信号とを用いて前記複数の成分を取得し、前記取得された複数の成分を用いて前記断層像を生成することを特徴とする請求項３乃至５のいずれか１項に記載の眼科撮影装置。
- 前記参照光の光路長と前記測定光の光路長との少なくとも１つを変更する光路長変更部をさらに有し、
前記演算処理部は、前記光路長変更部による光路長の変更が行われずに取得された前記変換された干渉信号を用いて折り返し成分が分離された前記複数の成分を取得し、前記取得された複数の成分を用いて前記断層像を生成することを特徴とする請求項１乃至６のいずれか１項に記載の眼科撮影装置。 - 測定光を照射した被検体からの戻り光と前記測定光に対応する参照光とによる干渉光を干渉信号として検出する検出部と、
前記検出された干渉信号をアナログ信号からディジタル信号に変換する変換部と、
前記変換された干渉信号を用いて、前記被検体の断層像を生成する演算処理部と、
前記参照光の光路長と前記測定光の光路長との少なくとも１つを変更する光路長変更部と、を有し、
前記演算処理部は、前記光路長変更部による光路長の変更が行われずに取得された折り返し成分が異なる複数の前記変換された干渉信号を用いて前記折り返し成分が分離された複数の成分を取得し、前記取得された複数の成分を用いて前記断層像を生成することを特徴とする眼科撮影装置。 - 波長掃引光源と、
波長掃引光源からの光を前記測定光と前記参照光とに分割する分割手段と、
前記波長掃引光源からの光のうち一部の光が通る光路が第一光路と前記第一光路に対して光路長差を有する第二光路とに分岐された干渉計として構成され、前記変換部によるサンプリングのクロックを生成するクロック生成部と、
をさらに有することを特徴とする請求項１乃至８のいずれか１項に記載の眼科撮影装置。 - 測定光を照射した被検体からの戻り光と前記測定光に対応する参照光とによる干渉光を干渉信号として検出する検出部と、
前記検出された干渉信号をアナログ信号からディジタル信号に変換する変換部と、を有する眼科撮影装置の制御方法であって、
前記変換して得た複数の信号に関する複数の成分であって、前記変換部のナイキスト周波数より高い周波数の成分と、前記変換部のナイキスト周波数より低い周波数の成分とを含む前記複数の成分を用いて、前記被検体の断層像を生成することを特徴とする眼科撮影装置の制御方法。 - 測定光を照射した被検体からの戻り光と前記測定光に対応する参照光とによる干渉光を干渉信号として検出する検出部と、
前記検出された干渉信号をアナログ信号からディジタル信号に変換する変換部と、
前記参照光の光路長と前記測定光の光路長との少なくとも１つを変更する光路長変更部と、を有する眼科撮影装置の制御方法であって、
前記光路長変更部による光路長の変更が行われずに取得された折り返し成分が異なる複数の前記変換された干渉信号を用いて前記折り返し成分が分離された複数の成分を取得し、前記取得された複数の成分を用いて断層像を生成することを特徴とする眼科撮影装の制御方法。 - 請求項１０または１１に記載の眼科撮影装置の制御方法をコンピュータに実行させるプログラム。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/JP2018/039174 WO2019082840A1 (ja) | 2017-10-27 | 2018-10-22 | 眼科撮影装置及びその制御方法 |
| US16/857,034 US20200281462A1 (en) | 2017-10-27 | 2020-04-23 | Ophthalmic imaging apparatus, control method for ophthalmic imaging apparatus, and computer-readable medium |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017208428 | 2017-10-27 | ||
| JP2017208428 | 2017-10-27 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2019080904A JP2019080904A (ja) | 2019-05-30 |
| JP7218122B2 true JP7218122B2 (ja) | 2023-02-06 |
Family
ID=66670642
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2018147777A Active JP7218122B2 (ja) | 2017-10-27 | 2018-08-06 | 眼科撮影装置及びその制御方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US20200281462A1 (ja) |
| JP (1) | JP7218122B2 (ja) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010014514A (ja) | 2008-07-03 | 2010-01-21 | Fujifilm Corp | 光断層画像化装置及び光断層画像化装置における干渉信号の処理方法 |
| JP2014176566A (ja) | 2013-03-15 | 2014-09-25 | Canon Inc | 画像処理システム、画像処理方法及びプログラム。 |
| JP2015114284A (ja) | 2013-12-13 | 2015-06-22 | キヤノン株式会社 | 光干渉断層計 |
| JP2015208574A (ja) | 2014-04-28 | 2015-11-24 | キヤノン株式会社 | 眼科撮影装置、その制御方法、およびプログラム |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8624964B2 (en) * | 2005-12-02 | 2014-01-07 | Koninklijke Philips N.V. | Depth dependent filtering of image signal |
| US20140340634A1 (en) * | 2013-04-05 | 2014-11-20 | Wasatch Photonics, Inc. | Optical coherence tomography systems and methods |
| US10184782B2 (en) * | 2013-05-17 | 2019-01-22 | Ninepoint Medical, Inc. | Frequency-domain optical coherence tomography with extended field-of-view and reduction of aliasing artifacts |
| JP2015102537A (ja) * | 2013-11-28 | 2015-06-04 | キヤノン株式会社 | 光干渉断層計 |
| JP2015226579A (ja) * | 2014-05-30 | 2015-12-17 | キヤノン株式会社 | 光干渉断層撮影装置及び光干渉断層撮影装置の制御方法 |
| JP2016075585A (ja) * | 2014-10-07 | 2016-05-12 | キヤノン株式会社 | 撮像装置、断層画像のノイズ低減方法、及びプログラム |
| US9906726B2 (en) * | 2015-04-22 | 2018-02-27 | Canon Kabushiki Kaisha | Image stabilization control apparatus, optical apparatus and storage medium storing image stabilizing control program |
| JP6765786B2 (ja) * | 2015-05-01 | 2020-10-07 | キヤノン株式会社 | 撮像装置、撮像装置の作動方法、情報処理装置、及び情報処理装置の作動方法 |
| JP2017003832A (ja) * | 2015-06-12 | 2017-01-05 | キヤノン株式会社 | 自動焦点調節装置および光学機器 |
| JP6370280B2 (ja) * | 2015-09-16 | 2018-08-08 | 富士フイルム株式会社 | 断層画像生成装置、方法およびプログラム |
-
2018
- 2018-08-06 JP JP2018147777A patent/JP7218122B2/ja active Active
-
2020
- 2020-04-23 US US16/857,034 patent/US20200281462A1/en not_active Abandoned
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010014514A (ja) | 2008-07-03 | 2010-01-21 | Fujifilm Corp | 光断層画像化装置及び光断層画像化装置における干渉信号の処理方法 |
| JP2014176566A (ja) | 2013-03-15 | 2014-09-25 | Canon Inc | 画像処理システム、画像処理方法及びプログラム。 |
| JP2015114284A (ja) | 2013-12-13 | 2015-06-22 | キヤノン株式会社 | 光干渉断層計 |
| JP2015208574A (ja) | 2014-04-28 | 2015-11-24 | キヤノン株式会社 | 眼科撮影装置、その制御方法、およびプログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| US20200281462A1 (en) | 2020-09-10 |
| JP2019080904A (ja) | 2019-05-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6805539B2 (ja) | 眼科撮像装置 | |
| JP6765786B2 (ja) | 撮像装置、撮像装置の作動方法、情報処理装置、及び情報処理装置の作動方法 | |
| JP4902721B2 (ja) | 光断層画像生成装置及び光断層画像生成方法 | |
| JP5752955B2 (ja) | 光断層像撮影装置 | |
| JP6426974B2 (ja) | データ処理方法及びｏｃｔ装置 | |
| JP6632266B2 (ja) | 撮像装置 | |
| JP6195334B2 (ja) | 撮像装置、撮像方法およびプログラム | |
| US9615736B2 (en) | Optical interference tomographic apparatus, and method for controlling optical interference tomographic apparatus | |
| WO2016120933A1 (en) | Tomographic imaging apparatus, tomographic imaging method, image processing apparatus, image processing method, and program | |
| US10524663B2 (en) | Phase measurement, analysis, and correction methods for coherent imaging systems | |
| JP2015198757A (ja) | 断層撮影装置 | |
| JP2017046975A (ja) | 眼科撮影装置及び眼科撮影プログラム | |
| JP2022176282A (ja) | 眼科装置、及びその制御方法 | |
| JP6452990B2 (ja) | データ処理方法及びｏｃｔ装置 | |
| JP6180877B2 (ja) | 光画像計測装置 | |
| JP7218122B2 (ja) | 眼科撮影装置及びその制御方法 | |
| JP5784100B2 (ja) | 画像形成装置及び画像形成方法 | |
| JP2019063044A (ja) | Ｏｃｔ装置、およびｏｃｔ装置制御プログラム | |
| JP2019154764A (ja) | 涙液層厚み測定装置及び方法 | |
| JP6779674B2 (ja) | Ｏｃｔ装置 | |
| WO2019082840A1 (ja) | 眼科撮影装置及びその制御方法 | |
| WO2019203091A1 (ja) | 画像処理装置、画像処理方法及びプログラム | |
| US11974806B2 (en) | Ophthalmologic information processing apparatus, ophthalmologic apparatus, ophthalmologic information processing method, and recording medium | |
| US11925411B2 (en) | Ophthalmologic information processing apparatus, ophthalmologic apparatus, ophthalmologic information processing method, and recording medium | |
| WO2023282339A1 (ja) | 画像処理方法、画像処理プログラム、画像処理装置及び眼科装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210728 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220712 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220906 |
|
| 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: 20221227 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230125 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 7218122 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |