CN117529268A - Oct信号处理方法、oct装置及程序 - Google Patents
Oct信号处理方法、oct装置及程序 Download PDFInfo
- Publication number
- CN117529268A CN117529268A CN202180099708.5A CN202180099708A CN117529268A CN 117529268 A CN117529268 A CN 117529268A CN 202180099708 A CN202180099708 A CN 202180099708A CN 117529268 A CN117529268 A CN 117529268A
- Authority
- CN
- China
- Prior art keywords
- light
- oct
- corrected
- signal
- optical system
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 13
- 230000003287 optical effect Effects 0.000 claims abstract description 117
- 238000005259 measurement Methods 0.000 claims abstract description 64
- 238000001514 detection method Methods 0.000 claims abstract description 43
- 230000001131 transforming effect Effects 0.000 claims abstract description 10
- 239000006185 dispersion Substances 0.000 claims description 105
- 238000000034 method Methods 0.000 claims description 25
- 230000008569 process Effects 0.000 claims description 17
- 238000012545 processing Methods 0.000 abstract description 56
- 238000012014 optical coherence tomography Methods 0.000 description 62
- 210000001508 eye Anatomy 0.000 description 48
- 238000010586 diagram Methods 0.000 description 13
- 238000003384 imaging method Methods 0.000 description 12
- 208000035405 autosomal recessive with axonal neuropathy spinocerebellar ataxia Diseases 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 10
- 238000010608 single channel array normalization Methods 0.000 description 10
- 238000012937 correction Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 9
- 230000004048 modification Effects 0.000 description 8
- 238000012986 modification Methods 0.000 description 8
- 239000005357 flat glass Substances 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 210000004220 fundus oculi Anatomy 0.000 description 3
- 210000001747 pupil Anatomy 0.000 description 3
- 210000001525 retina Anatomy 0.000 description 3
- 238000013473 artificial intelligence Methods 0.000 description 2
- 210000005252 bulbus oculi Anatomy 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 206010020675 Hypermetropia Diseases 0.000 description 1
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- BJQHLKABXJIVAM-UHFFFAOYSA-N bis(2-ethylhexyl) phthalate Chemical compound CCCCC(CC)COC(=O)C1=CC=CC=C1C(=O)OCC(CC)CCCC BJQHLKABXJIVAM-UHFFFAOYSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 210000004087 cornea Anatomy 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000004305 hyperopia Effects 0.000 description 1
- 201000006318 hyperopia Diseases 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 208000001491 myopia Diseases 0.000 description 1
- 230000004379 myopia Effects 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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/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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Ophthalmology & Optometry (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Eye Examination Apparatus (AREA)
Abstract
一种OCT信号处理方法,将从光源射出的光分离为测量光和参照光,使用物镜光学系统将上述测量光向被检眼引导,在将上述参照光引导到参照光学系统后,处理器对通过接收从上述被检眼反射的测量光和上述参照光发生干涉的干涉光而得到的检测信号进行处理,在OCT信号处理方法中包括:通过对上述检测信号进行傅里叶变换而生成OCT信号的步骤;检测上述OCT信号中的高峰部分的步骤;和通过从上述OCT信号除去上述高峰部分而生成第1修正OCT信号的步骤。
Description
技术领域
本公开的技术涉及OCT信号处理方法、OCT装置及程序。
背景技术
在美国专利第9109870号说明书中,公开了一种生成被检眼的断层图像的装置。期望得到抑制了伪影的断层图像。
发明内容
本公开的技术的第1方案的OCT信号处理方法将从光源射出的光分离为测量光和参照光,使用物镜光学系统将上述测量光向被检眼引导,在将上述参照光引导到参照光学系统后,处理器对通过接收从上述被检眼反射的测量光和上述参照光发生干涉的干涉光而得到的检测信号进行处理,该OCT信号处理方法包括:通过对上述检测信号进行傅里叶变换而生成OCT信号的步骤;检测上述OCT信号中的高峰部分的步骤;和通从上述OCT信号除去上述高峰部分而生成第1修正OCT信号的步骤。
本公开的技术的第2方案的OCT装置具备处理器,将从光源射出的光分离为测量光和参照光,使用物镜光学系统将上述测量光向被检眼引导,在将上述参照光引导到参照光学系统后,上述处理器对通过接收从上述被检眼反射的测量光和上述参照光发生干涉的干涉光而得到的检测信号进行处理,上述处理器执行以下步骤:通过对上述检测信号进行傅里叶变换而生成OCT信号的步骤;从上述OCT信号检测高峰部分的步骤;和通过除去上述高峰部分而生成第1修正OCT信号的步骤。
本公开的技术的第3方案的程序将从光源射出的光分离为测量光和参照光,使用物镜光学系统将上述测量光向被检眼引导,在将上述参照光引导到参照光学系统后,使计算机对通过接收从上述被检眼反射的测量光和上述参照光发生干涉的干涉光而得到的检测信号进行处理,上述程序使上述计算机执行以下步骤:通过对上述检测信号进行傅里叶变换而生成OCT信号的步骤;从上述OCT信号检测高峰部分的步骤;和通过除去上述高峰部分而生成第1修正OCT信号的步骤。
附图说明
图1是表示第1实施方式的眼科系统100的结构的框图。
图2是表示眼科装置110的结构的框图。
图3A是表示OCT单元20P的电气系统的结构的框图。
图3B是表示摄影光学系统19的结构的图。
图4是控制装置16的CPU16A的功能框图。
图5A是表示由控制装置16进行的处理的流程图。
图5B是表示由控制装置16进行的图5A的步骤258的高峰除去处理的详情的流程图。
图6是表示不使参照光透射过色散产生部25P的情况下的OCT信号的图。
图7A是表示使参照光透射过色散产生部25P的情况下的OCT信号的图。
图7B是表示使参照光透射过色散产生部25P的情况下的、从OCT信号除去了高峰部分的情况下的信号的图。
图7C是表示使参照光透射过色散产生部25P的情况下的、从OCT信号除去高峰部分并进行了除去范围的补全的情况下的信号的图。
图8是表示第1变形例的OCT单元20Q的电气系统的结构的框图。
图9是表示第2变形例的OCT单元20R的电气系统的结构的框图。
图10是表示第3变形例的OCT单元20P的电气系统的结构的框图。
具体实施方式
以下,参照附图详细说明本公开的技术的实施方式。
参照图1,说明第1实施方式的眼科系统100的结构。如图1所示,眼科系统100具备眼科装置110、眼轴长度测量装置120、管理服务器装置(以下,称为“服务器”)140和图像显示装置(以下,称为“观看器”)150。眼科装置110获取眼底图像。眼轴长度测量装置120测量患者的眼轴长度。服务器140将通过由眼科装置110拍摄患者的眼底而得到的眼底图像与患者的ID对应地存储。观看器150显示从服务器140获取到的眼底图像等医疗信息。
眼科装置110、眼轴长度测量装置120、服务器140及观看器150经由网络130而相互连接。网络130是LAN、WAN、因特网或广域以太网等任意网络。例如,在眼科系统100构筑于一个医院的情况下,能够对网络130采用LAN。
观看器150是客户端服务器系统中的客户端,经由网络连接多台。另外,为了担保系统的冗余性,服务器140也可以经由网络连接多台。或者,若眼科装置110具备图像处理功能及观看器150的图像阅览功能,则眼科装置110能够在独立(stand-alone)状态下实现眼底图像的获取、图像处理及图像阅览。该情况下,可以省略观看器150。另外,若服务器140具备观看器150的图像阅览功能,则在眼科装置110和服务器140的结构下,能够实现眼底图像的获取、图像处理及图像阅览。该情况下,可以省略观看器150。
此外,使用其他眼科设备(视野测量、眼压测量等检查设备)或AI(ArtificialIntelligence)进行图像解析的诊断支援装置可以经由网络130与眼科装置110、服务器140及观看器150连接。
接下来,参照图2说明眼科装置110的结构。
为了便于说明,将激光扫描检眼镜(Scanning Laser Ophthalmoscope)称为“SLO”。另外,将光学相干断层扫描仪(Optical Coherence Tomography)称为“OCT”。
此外,在眼科装置110设置于水平面的情况下将水平方向设为“X方向”,将相对于水平面垂直的方向设为“Y方向”,将连结被检眼12的前眼部的瞳孔的中心与眼球的中心的方向设为“Z方向”。因此,X方向、Y方向及Z方向相互垂直。
眼科装置110具备SLO单元18、OCT单元20P、摄影光学系统19及控制装置16,通过对被检眼12的眼底进行拍摄而获取眼底图像。以下,将由SLO单元18获取到的二维眼底图像称为SLO图像。另外,将基于由OCT单元20P获取到的OCT数据创建的视网膜的断层图像或正面图像(en-face图像)等称为OCT图像。
控制装置16具备具有CPU(Central Processing Unit(中央处理装置))16A、RAM(Random Access Memory)16B、ROM(Read-Only memory)16C及输入输出(I/O)端口16D的计算机。在ROM16C中存储有后述的信号处理程序。
CPU16A是本公开的技术的“处理器”的一例。ROM16C中存储的信号处理程序相当于本公开的技术的“程序”的一例。
控制装置16具备经由I/O端口16D与CPU16A连接的输入/显示装置16E。输入/显示装置16E具有显示被检眼12的图像、从用户受理各种指示的图形用户界面。作为图形用户界面,可列举触摸面板-显示器。
控制装置16基于通过上述拍摄得到的数据生成被检眼12的图像。控制装置16具备与I/O端口16D连接的通信接口(I/F)16F。眼科装置110经由通信接口(I/F)16F及网络130与眼轴长度测量装置120、服务器140及观看器150连接。
如上述那样,在图2中,眼科装置110的控制装置16具备输入/显示装置16E,但本公开的技术并不限定于此。例如,眼科装置110的控制装置16也可以不具备输入/显示装置16E,而具备物理上独立于眼科装置110的单独的输入/显示装置。该情况下,该显示装置具备在控制装置16的CPU16A的控制下动作的图像处理器单元。图像处理器单元可以基于CPU16A输出指示的图像信号,显示SLO图像等。
如图3A所示,OCT单元20P包括光源20A、检测器20B、第1光耦合器20C、参照光学系统20D1、准直透镜20E及第2光耦合器20F,执行如后述那样指定的进行OCT拍摄的部分的OCT拍摄。
如图3B所示,摄影光学系统19具备OCT单元20P用的聚焦机构19A、使来自OCT单元20P的光在X方向上扫描的第1扫描器19B、SLO单元18用的聚焦机构19C、以及使来自SLO单元18的光在X方向上扫描的第2扫描器19D。摄影光学系统19具备分色镜19E,将来自OCT单元20P的光朝向被检眼12反射且将由被检眼12的眼底反射的光朝向OCT单元20P反射,将来自SLO单元18的光向被检眼12透射且将由被检眼12的眼底反射的光向SLO单元18透射。摄影光学系统19具备公共扫描器19SY及物镜19T,公共扫描器使由分色镜19E反射的来自OCT单元20P的光和透射过分色镜19E的来自SLO单元18的光在Y方向上扫描。
物镜19T是本公开的技术的“物镜光学系统”的一例。
物镜19T具备通过将多个透镜以粘结剂粘结而形成的未图示的接合透镜(也称为贴合透镜)。第1扫描器19B、第2扫描器19D及公共扫描器19SY是能够使光束偏转的光学元件。例如,作为第1扫描器19B及公共扫描器19SY,能够使用镜式检流计(galvanometermirror)等。作为第2扫描器19D,能够使用多面镜(polygon mirror)等。
如图3A所示,OCT单元20P的光源20A射出多个不同波长的光。从光源20A射出的光在第1光耦合器20C分支。分支出的一方的光作为测量光,在准直透镜20E成为平行光后,入射到摄影光学系统19。测量光通过聚焦机构19A进行聚焦调整,通过第1扫描器19B和公共扫描器19SY,分别在X方向及Y方向上扫描。扫描光经由瞳孔,照射到眼底。由眼底反射的测量光经由摄影光学系统19向OCT单元20P入射,并经由准直透镜20E及第1光耦合器20C向第2光耦合器20F入射。
从光源20A射出并在第1光耦合器20C分支出的另一方的光作为参照光,向参照光学系统20D1入射,经由参照光学系统20D1入射到第2光耦合器20F。参照光学系统20D1除了多个透镜以外,具备使参照光产生波长色散的色散产生部25P。色散产生部25P是本公开的技术的“色散部”的一例。
色散产生部25P由以对参照光赋予规定量的色散的方式确定介质及厚度的平行平板玻璃构成。通过色散产生部25P在参照光产生的波长色散是与色散产生部25P的介质及厚度相应的固有现象。此外,在本实施方式中,通过规定介质在参照光的光路形成了色散产生部25P。对透射过该色散产生部25P的参照光赋予了基于色散产生部25P的色散。
根据色散产生部25P的色散的大小,干涉信号的深度方向上的分辨率的降低情况发生变化。干涉信号的分辨率越降低则后述的伪影越容易识别。因此,在以不赋予色散的情况下的分辨率设为基准的情况下,期望利用色散产生部25P对参照光赋予使分辨率降低二分之一以上的色散。但是,若过于增大色散产生部25P的色散(例如使析像度降低到百分之一的色散),则干涉信号的分辨率过低,因此难以通过色散修正进行修正。
由被检眼反射的测量光和被赋予了基于色散产生部25P的色散的参照光入射到第2光耦合器20F。然后参照光和测量光在第2光耦合器20F发生干涉而生成干涉光。干涉光由检测器20B接收。在控制装置16使用由检测器20B检测出的干涉信号生成OCT图像。
在此,由被检眼反射的测量光与被赋予了规定量的色散的参照光发生干涉而得到的干涉光受到基于色散产生部25P的色散的影响,而深度方向上的分辨率降低。但是,由摄影光学系统19内的物镜19T等光学元件产生的自干涉光由于成为测量光的干涉,所以不会受到基于色散产生部25P的色散的影响。在产生了该自干涉光的情况下,由摄影光学系统19内的物镜19T等光学元件产生的自干涉光与由被检眼反射的测量光和被赋予了规定量的色散的参照光发生干涉得到的干涉光一起,包含于检测器20B接收的光。
本公开的技术通过将被赋予了基于色散产生部25P的色散的参照光与干涉信号的色散修正组合,能够除去因由摄影光学系统19内的光学元件产生的自干涉光导致的伪影。除去伪影的方法的详情将后述。
图2所示的SLO单元18具备多个光源,例如,B光(蓝色光)的光源、G光(绿色光)的光源、R光(红色光)的光源、以及IR光(红外线(例如,近红外光))的光源、和将来自这些光源的光反射或透射并引导到通向摄影光学系统19的一个光路的多个光学系统。
从SLO单元18入射到摄影光学系统19的光由聚焦机构19C进行聚焦调整,通过第2扫描器19D和公共扫描器19SY分别在X方向及Y方向上扫描。扫描光经由瞳孔,照射到被检眼12的后眼部。由眼底反射的反射光经由摄影光学系统19向SLO单元18入射。
SLO单元18具备与上述多个光源对应的多个光检测元件、和将来自各光源的光的上述反射光向与各光源对应的光检测元件引导的多个分光器。然后,使用由各光检测元件检测出的信号在控制装置16生成SLO图像。
图1所示的眼轴长度测量装置120测量被检眼12的眼轴方向的长度即眼轴长度。测量出的眼轴长度被发送到服务器140。服务器140将患者的眼轴长度与患者名ID对应地存储。
图1所示的服务器140及观看器150分别由具有未图示的CPU、RAM、ROM及输入输出(I/O)端口等的计算机构成。在输入输出(I/O)端口连接有存储装置、显示器、鼠标、键盘及通信接口(I/F)。能够经由通信接口(I/F)与连接于网络130的其他装置进行通信。
接下来,参照图4,说明通过控制装置16的CPU16A执行信号处理程序而实现的各种功能。信号处理程序具备获取应处理的干涉信号的获取功能、获取到的干涉信号的信号处理功能、生成OCT图像的图像处理功能。通过CPU16A执行具有该各功能的信号处理程序,CPU16A如图4所示,作为获取部52、信号处理部54及图像处理部56发挥功能。
接下来,使用图5A、图5B,详细说明由控制装置16进行的信号处理。通过控制装置16的CPU16A执行信号处理程序,实现了图5A、图5B的流程图所示的信号处理。该信号处理在经由输入/显示装置16E指定了在患者的被检眼的眼底的视网膜中进行OCT拍摄的部分且操作了显示于输入/显示装置16E的未图示的开始按钮的情况下开始。
在图5A的步骤252中,获取部52扫描被检眼12。具体地说,获取部52控制OCT单元20P和摄影光学系统19,执行由用户指定的被检眼12的拍摄位置的OCT拍摄。然后,在步骤254中,获取部52获取通过OCT拍摄得到的干涉信号CS。在被检眼12的眼底为该拍摄位置的情况下,以测量光在眼底形成焦点的方式控制摄影光学系统19。由眼底反射的测量光和被色散产生部25P赋予了色散的参照光发生干涉而得到的干涉光由检测器20B接收(检测)。检测器20B将干涉信号CS(检测信号)输出到控制装置16并存储于RAM16B。然后,获取部52从RAM16B获取干涉信号CS。如上述那样,在通过摄影光学系统19的光学元件产生自干涉光的情况下,与由眼底反射的测量光和被色散产生部25P赋予了色散的参照光发生干涉而得到的干涉光一起,该自干涉光也入射到检测器20B。由此,该情况下,在干涉信号CS包含基于来自被检眼12的反射光产生的干涉光和基于摄影光学系统产生的自干涉光。
在步骤256中,图像处理部56通过对干涉信号CS执行FFT处理(快速傅里叶变换Fast Fourier Transform),获取OCT信号TS。
在步骤258中,图像处理部56从OCT信号TS中检测亮度急剧变化的尖锐的信号的部分即高峰部分,执行将该高峰部分除去的高峰除去处理。该高峰除去处理使用图5B详细进行说明。然后,图像处理部56生成高峰部分被除去的第1修正OCT信号TSA。
在步骤260中,图像处理部56通过对第1修正OCT信号TSA执行逆FFT处理(InversFast Fourier Transform处理),生成第1修正干涉信号CSA(噪声被除去的修正检测信号)。在此所说的色散为波长色散,因此为了进行基于频率范围的信号修正,对第1修正OCT信号TSA进行逆FFT,生成第1修正干涉信号CSA。
在步骤262中,图像处理部56执行色散修正处理。具体地说,图像处理部56从第1修正干涉信号CSA修正由色散产生部25P赋予的色散的影响。能够对该色散修正处理任意适用公知的色散修正算法。通过进行色散修正处理,图像处理部56生成修正了由色散产生部25P赋予的色散的影响、并且除去了伪影的第2修正检测信号CSB。
在步骤264中,图像处理部56通过对修正了由色散产生部25P赋予的色散的影响且除去了伪影的第2修正检测信号CSB执行FFT处理,获取修正了色散的影响且除去了伪影的第2修正OCT信号TSB。
在步骤266中,图像处理部56基于第2修正OCT信号TSB,进行眼球形状修正及全范围处理等,生成被检眼12的OCT图像(断层图像)。图像处理部56可以进一步进行将得到的视网膜和角膜等的断层图像按层可视化的分割处理、以计算求出分割出的各层的层厚的处理、生成Enface图像等的处理等。此外,图像处理部56在OCT信号为以C扫描得到的数据为三维数据(OCT体数据)的情况下,也可以从OCT体数据生成en-face图像。
在步骤266中,图像处理部56将生成的断层图像和第2修正OCT信号TSB输出到RAM16B或未图示的外部的存储装置。通过步骤266的处理结束,信号处理(即,信号处理程序的执行)结束。
像这样,通过从OCT信号除去因由摄影光学系统19的光学元件产生的自干涉光形成的高峰部分,生成除去了高峰部分的第1修正OCT信号,对该第1修正OCT信号进行逆傅里叶变换而生成修正检测信号,生成对修正检测信号进行了色散修正的修正干涉信号,对修正干涉信号进行傅里叶变换,得到第2修正OCT信号。由此,能够得到抑制了伪影的OCT图像(断层图像)。
接下来,使用图5B的流程图说明图5A的步骤258的高峰除去处理的详情。
在步骤2580中,图像处理部56将参数n设为1。在将构成OCT信号TS的A-SCAN的数量设为N时,参数n取1~N的值。在此N为1以上的自然数。
接着,在步骤2581中,图像处理部56从OCT信号TS提取第n个A-SCAN的信号ASn_#n。
然后,在步骤2582中,图像处理部56判断在信号ASn中是否存在高峰部分。在此,信号ASn所含的高峰部分是指例如如图7A所示,信号ASn的强度值的变化急剧的特征性的高峰部分111DA。特征性的高峰部分111DA具体地说,是强度值的半值宽度小于规定值的部分,更具体地说,例如,是半值全宽(Full Width at Half Maximum,FWHM)例如小于8um的部分。该情况下,图像处理部56从信号ASn中判断信号ASn中是否存在半值全宽小于8um的部分。或者能够使用公知的各种各样的高峰检测方法。
图7A是表示存在高峰部分111DA的A-SCAN的信号ASn的波形的图,纵轴示出信号ASn的强度值,横轴示出深度。高峰部分111DA在横轴上位于深度DP的位置。该高峰部分111DA是因由摄影光学系统19内的光学元件产生的测量光的自干涉光(该自干涉光是不受到色散影响的自干涉光)得到的检测信号的强度值。另一方面,图7A的信号部部分112DA是由被检眼12反射的测量光和被色散产生部25P赋予了色散的参照光发生干涉得到的干涉光的检测信号的强度值。信号部部分112DA因色散的影响而深度方向的析像度显著降低。
另一方面,参照光路中不存在色散产生部25P的情况下的信号ASn如图6所示,包含基于由被检眼12反射的测量光和未被赋予色散的参照光发生干涉得到的干涉光形成的检测信号的强度值的部分112A、以及基于由摄影光学系统19内的光学元件产生的测量光的自干涉光形成的检测信号的强度值即高峰部分111A。高峰部分111A位于与图7A的高峰部分111DA相同的横轴上的位置即深度DP。这是因为,图7A和图6所示的信号Asn的条件不同仅为参照光路中的参照光学系统20D1是否设有色散产生部25P这一不同之处,对于通过测量光路的测量光而言为相同条件。
假设在横轴的深度DP的位置存在与被检眼12的构造物对应的强度值,也会因赋予色散而在OCT信号中模糊,能够明确区分出因由摄影光学系统19内的光学元件产生的测量光的自干涉光形成的检测信号的强度值。也就是说,通过在参照光路中配置色散产生部25P并对参照光赋予规定的色散,从而在OCT信号中,能够在横轴的深度DP的位置区分与被检眼12的构造物对应的强度值和因由摄影光学系统19内的光学元件产生的测量光的自干涉光形成的强度值。
这样的高峰部分对于被检眼12的OCT信号而言为噪声。成为因由眼科装置110的摄影光学系统19的物镜19T和其他光学元件等产生的自干涉光导致的与原本的被检眼12的构造不对应的不需要像(以下,称为伪影(artifact))的原因。
像这样在信号ASn中,作为噪声,包含尖锐的高峰部分的理由如下所述。物镜19T具备通过将多个透镜以粘结剂粘结而形成的接合透镜。测量光当入射到物镜19T时,不仅透射过多个透镜,还在接合透镜中的多个透镜的透镜面反射。由于测量光在透镜面反射,所以在反射的光与测量光之间产生自干涉(标准具(etalon))。因该自干涉而导致在测量光包含噪声。此外,在被检眼12近视的情况下,与远视的情况相比,噪声更多产生,断层图像中产生更大的伪影。
另外,若因眼科装置110的摄影光学系统19的物镜19T等光学元件导致的伪影的存在是既知的,则在OCT信号中预先判明高峰部分111A的横轴上的深度DP。由此,在因摄影光学系统19的物镜19T等光学元件导致的伪影的存在是既知的情况下,通过具有表示被检眼12的拍摄位置与高峰部分111A的横轴上的深度DP的关系的查找表,而无需进行步骤2582的判断。在步骤2583中通过参照该查找表而可知是否为存在高峰部分111A的A扫描,在存在高峰部分111A的A扫描的情况下,判明横轴上的深度DP,能够除去该深度DP的强度值。
在步骤2582中,在由图像处理部56判断成信号ASn中不存在尖锐的高峰部分的情况下(步骤2582:否),前进到步骤2585。然后,在步骤2585中图像处理部56将信号ASn作为第n个A-SCAN的信号保存到RAM64。
在步骤2582中判断成信号ASn中存在尖锐的高峰部分的情况下(步骤2584:是),进入步骤2583。然后,在步骤2583中图像处理部56如图7B所示将高峰部分111DA的范围中的强度值设为零。由此,图像处理部56从信号ASn除去高峰部分111DA。
另外,也可以不是将高峰部分111DA的强度值设为零,而是以成为高峰部分的周边的强度值的平均值的方式将强度值减小至该平均值。
对于高峰部分111DA,由于为因由摄影光学系统19内的光学元件产生的测量光的自干涉光形成的检测信号的强度值,所以即使从信号ASn去除,也能够无视对因色散而模糊的与被检眼12的构造物对应的强度值的影响。
接着,在步骤2584中,图像处理部56进行补全存在过高峰部分111DA的范围的强度值111DC的处理,如图7C所示,生成将被除去了的高峰部分插补的信号CASn。然后,图像处理部56将信号CASn设为修正后的第n个A-SCAN的信号。
作为进行补全的处理,使用线性补全、样条补全、利用了高斯过程回归的补全等公知的补全手法,基于高峰部分111DA的周边的信号部部分112DA的强度值,能够求出存在过高峰部分111DA的范围的强度值111CD。
然后,在步骤2585中图像处理部56将信号CASn设为第n个A-SCAN的信号并保存到RAM64。
接着,在步骤2586中图像处理部56判断n是否为N,若n不为N(步骤2586:否),则进入步骤2587。在步骤2587中,图像处理部56设为n=n+1,重复步骤2581~2585的处理。
另一方面,若n为N(步骤2586:是),则对构成OCT信号TS的N个A-SCAN全部的高峰除去的处理结束,得到修正OCT信号TSA,图5B的流程图结束。然后,图5A的步骤260开始。
如以上说明那样在本实施方式中,由于从由被检眼反射的测量光和被赋予了基于色散产生部25P的色散的参照光发生干涉的干涉光得到OCT信号,所以能够有效除去因由眼科装置110的摄影光学系统19的物镜19T和其他光学元件等产生的测量光的自干涉形成的A-SCAN信号的噪声。然后,对除去了因测量光的自干涉形成的A-SCAN信号的噪声的OCT信号进行逆傅里叶变换而得到修正干涉信号。然后通过对该修正干涉信号进行色散修正而修正由色散产生部25P赋予的色散的影响。然后,能够再次对色散修正后的修正干涉信号进行傅里叶变换而得到除去了伪影的OCT信号。
接下来,说明本公开的技术的变形例。此外,对与上述实施方式相同的部分标注相同的附图标记并省略其说明,以不同部分为主进行说明。
在上述实施方式中,OCT单元20P的参照光学系统20D1具备色散产生部25P。在图3A所示的实施例中,将色散产生部25P由平行平板玻璃构成,但本公开的技术并不限定于此。
图8中示出了第1变形例的摄影光学系统19和OCT单元20Q的电气系统的结构。如图8所示,可以取代图3A所示的参照光学系统20D1,将以对参照光赋予上述色散的方式确定了长度的光纤25Q用作色散产生部。光纤25Q配置在参照光的光路上。
图9中示出了第2变形例的摄影光学系统19和OCT单元20R的电气系统的结构。如图9所示,可以取代图3A所示的第1光耦合器20C及第2光耦合器20F而具备光耦合器20G,取代图3A所示的参照光学系统20D1而使用具备多个透镜、色散产生部25R及反射镜的参照光学系统20D2。从光源20A射出的光在光耦合器20G分支。分支出的一方的光成为测量光。由眼底反射的测量光经由摄影光学系统19向OCT单元20R入射,并入射到光耦合器20G。从光源20A射出并在光耦合器20G分支出的另一方的光作为参照光,向参照光学系统20D2入射,从多个透镜及色散产生部25R透射,在反射镜反射,从多个透镜及色散产生部25R透射,并入射到光耦合器20G。色散产生部25R是以对参照光赋予上述色散的方式确定了介质及厚度的平行平板玻璃。入射到光耦合器20G的这些光,即由眼底反射的测量光和参照光在光耦合器20G发生干涉而生成干涉光。干涉光由检测器20B接收。
色散产生部25R配置在参照光的光路上。
另外,在上述实施方式中,将色散产生部配置在参照光的光路上,对参照光赋予色散。本公开的技术并不限定于此。也可以将色散产生部仅配置在测量光的光路,从而取代参照光仅对测量光赋予色散,或将色散产生部配置在参照光的光路和测量光的光路各自上,从而对参照光和测量光赋予色散。
图10中示出了第3变形例的摄影光学系统19和OCT单元20P的电气系统的结构。图10中示出了对参照光和测量光赋予色散的例子。如图10所示,OCT单元20P具备图3A所示的参照光学系统20D1,并且在摄影光学系统19,以供测量光透射的方式具备色散产生部19P。通过色散产生部25P和色散产生部19P,对参照光和测量光赋予使各光路的色散特性互不相同的色散。色散产生部25P和色散产生部19P是以对参照光和测量光赋予上位色散的方式确定了介质及厚度的平行平板玻璃。此外,在图8所示的例子中也是,可以在摄影光学系统19以供测量光透射的方式具备色散产生部19P。在第3变形例中也是,在OCT信号中,能够区分与被检眼的构造物对应的强度值和因由摄影光学系统内的光学元件产生的测量光的自干涉光得到的强度值。
如以上说明那样,在本公开的技术的实施方式中,由于除去因在物镜光学系统中产生的测量光的自干涉形成的OCT信号中的噪声,所以能够得到抑制了伪影的断层图像。
在以上说明的实施方式中,作为产生上述自干涉而在断层图像产生伪影的至少一个光学系统,说明了在物镜19T具备接合透镜的情况,但本公开的技术并不限定于此。
例如,第1,可以是物镜19T包括隔着规定厚度的极薄空气层而接近配置的多枚(例如,2枚)透镜的情况。该情况下的多枚透镜各自以大致相等的曲率面(曲率半径的面)接近配置。
第2,可以是物镜19T包括无粘结剂而贴紧的多枚(例如,2枚)透镜的情况。该情况下的多枚透镜是相等的曲率面(曲率半径的面)接触且光学上贴紧的透镜组。
第3,可以是物镜19T包括接近配置的多枚(例如,2枚)凸透镜的情况。该情况下的多枚透镜的彼此顶点不接触。
第4,可以是物镜19T包括彼此顶点接触的多枚(例如,2枚)凸透镜的情况。该情况下的多枚透镜的彼此顶点接触。
第5,可以是物镜19T包括平行平板玻璃的情况。
另外,在以上说明的各例中,说明了在物镜19T设有产生上述自干涉而在断层图像产生伪影的至少一个光学系统的情况,但本公开的技术并不限定于此。本公开的技术也能够适用于在从被检眼12到检测器20B之间,例如,在摄影光学系统19的聚焦机构19A、19C及/或到OCT单元20P内的检测器20B的光路上存在产生上述自干涉而在断层图像产生伪影的至少一个光学系统的情况。
在本公开中,各结构要素(装置等)只要不产生矛盾,则可以仅存在一个或存在两个以上。
在以上说明的各例中,例示了通过利用计算机的软件结构实现图像处理的情况,但本公开的技术并不限定于此。例如,也可以取代利用计算机的软件结构,仅通过FPGA(Field-Programmable Gate Array)或ASIC(Application Specific IntegratedCircuit)等硬件结构执行图像处理。也可以利用软件结构执行图像处理中的一部分处理,利用硬件结构执行剩余的处理。
以上说明的各信号处理只不过为一例。因此,当然可以在不脱离主旨的范围内删除不需要的步骤、追加新的步骤或替换处理顺序。
本说明书中记载的全部文献、专利申请及技术标准与通过参照具体引入且单独记载各个文献、专利申请及技术标准的情况同样地,通过参照引入于本说明书中。
Claims (11)
1.一种OCT信号处理方法,将从光源射出的光分离为测量光和参照光,使用物镜光学系统将所述测量光向被检眼引导,在将所述参照光引导到参照光学系统后,处理器对通过接收从所述被检眼反射的测量光和所述参照光发生干涉的干涉光而得到的检测信号进行处理,
所述OCT信号处理方法中包括:
通过对所述检测信号进行傅里叶变换而生成OCT信号的步骤;
检测所述OCT信号中的高峰部分的步骤;和
通过从所述OCT信号除去所述高峰部分而生成第1修正OCT信号的步骤。
2.根据权利要求1所述的OCT信号处理方法,其中,
通过设于所述参照光学系统的色散部而对所述参照光赋予了规定量的色散,
还包括:
通过对所述第1修正OCT信号进行逆傅里叶变换而生成第1修正检测信号的步骤;
对所述第1修正检测信号所含的所述色散进行修正并输出第2修正检测信号的步骤;和
通过对所述第2修正检测信号进行傅里叶变换而生成修正了所述色散的第2修正OCT信号的步骤。
3.根据权利要求2所述的OCT信号处理方法,其中,
包括基于所述第2修正OCT信号生成断层图像的步骤。
4.根据权利要求1至3中任一项所述的OCT信号处理方法,其中,
所述高峰部分是强度值的半值宽度小于规定值的部分。
5.根据权利要求1至4中任一项所述的OCT信号处理方法,其中,
所述高峰部分因由所述物镜光学系统产生的自干涉光而产生。
6.一种OCT装置,具备处理器,所述OCT装置将从光源射出的光分离为测量光和参照光,使用物镜光学系统将所述测量光向被检眼引导,在将所述参照光引导到参照光学系统后,所述处理器对通过接收从所述被检眼反射的测量光和所述参照光发生干涉的干涉光而得到的检测信号进行处理,其中,
所述处理器执行以下步骤:
通过对所述检测信号进行傅里叶变换而生成OCT信号的步骤;
从所述OCT信号检测高峰部分的步骤;和
通过除去所述高峰部分而生成第1修正OCT信号的步骤。
7.根据权利要求6所述的OCT装置,其中,
通过设于所述参照光学系统的色散部而对所述参照光赋予了规定量的色散,
所述处理器还执行以下步骤:
通过对所述第1修正OCT信号进行逆傅里叶变换而生成第1修正检测信号的步骤;
对所述第1修正检测信号所含的所述色散进行修正并输出第2修正检测信号的步骤;和
通过对所述第2修正检测信号进行傅里叶变换而生成修正了所述色散的第2修正OCT信号的步骤。
8.根据权利要求7所述的OCT装置,其中,
所述处理器还执行基于所述第2修正OCT信号生成断层图像的步骤。
9.一种程序,将从光源射出的光分离为测量光和参照光,使用物镜光学系统将所述测量光向被检眼引导,在将所述参照光引导到参照光学系统后,使计算机对通过接收从所述被检眼反射的测量光和所述参照光发生干涉的干涉光得到的检测信号进行处理,
所述程序使所述计算机执行以下步骤:
通过对所述检测信号进行傅里叶变换而生成OCT信号的步骤;
从所述OCT信号检测高峰部分的步骤;和
通过除去所述高峰部分而生成第1修正OCT信号的步骤。
10.根据权利要求9所述的程序,其中,
通过设于所述参照光学系统的色散部而对所述参照光赋予了规定量的色散,
使所述计算机还执行以下步骤:
通过对所述第1修正OCT信号进行逆傅里叶变换而生成第1修正检测信号的步骤;
对所述第1修正检测信号所含的所述色散进行修正并输出第2修正检测信号的步骤;和
通过对所述第2修正检测信号进行傅里叶变换而生成修正了所述色散的第2修正OCT信号的步骤。
11.根据权利要求10所述的程序,其中,
使所述计算机还执行基于所述第2修正OCT信号生成断层图像的步骤。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2021/024001 WO2022269876A1 (ja) | 2021-06-24 | 2021-06-24 | Oct信号処理方法、oct装置、およびプログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117529268A true CN117529268A (zh) | 2024-02-06 |
Family
ID=84544276
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202180099708.5A Pending CN117529268A (zh) | 2021-06-24 | 2021-06-24 | Oct信号处理方法、oct装置及程序 |
Country Status (4)
Country | Link |
---|---|
EP (1) | EP4344613A1 (zh) |
JP (1) | JPWO2022269876A1 (zh) |
CN (1) | CN117529268A (zh) |
WO (1) | WO2022269876A1 (zh) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012024551A (ja) * | 2010-06-25 | 2012-02-09 | Panasonic Corp | 光断層画像取得装置 |
JP6195334B2 (ja) | 2012-08-30 | 2017-09-13 | キヤノン株式会社 | 撮像装置、撮像方法およびプログラム |
JP6779690B2 (ja) * | 2016-07-27 | 2020-11-04 | 株式会社トプコン | 眼科画像処理装置及び眼科撮影装置 |
JP6776076B2 (ja) * | 2016-09-23 | 2020-10-28 | 株式会社トプコン | Oct装置 |
-
2021
- 2021-06-24 WO PCT/JP2021/024001 patent/WO2022269876A1/ja active Application Filing
- 2021-06-24 EP EP21947161.2A patent/EP4344613A1/en active Pending
- 2021-06-24 JP JP2023529381A patent/JPWO2022269876A1/ja active Pending
- 2021-06-24 CN CN202180099708.5A patent/CN117529268A/zh active Pending
Also Published As
Publication number | Publication date |
---|---|
JPWO2022269876A1 (zh) | 2022-12-29 |
EP4344613A1 (en) | 2024-04-03 |
WO2022269876A1 (ja) | 2022-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5511437B2 (ja) | 光断層像撮影装置 | |
US11941788B2 (en) | Image processing method, program, opthalmic device, and choroidal blood vessel image generation method | |
JP6195334B2 (ja) | 撮像装置、撮像方法およびプログラム | |
JP2023009530A (ja) | 画像処理方法、画像処理装置、及びプログラム | |
US10660514B2 (en) | Image processing apparatus and image processing method with generating motion contrast image using items of three-dimensional tomographic data | |
US10561311B2 (en) | Ophthalmic imaging apparatus and ophthalmic information processing apparatus | |
JP6909109B2 (ja) | 情報処理装置、情報処理方法、及びプログラム | |
JP6411792B2 (ja) | 調節機能評価装置 | |
JP2022040372A (ja) | 眼科装置 | |
JP7368581B2 (ja) | 眼科装置、及び眼科情報処理装置 | |
JP2019177032A (ja) | 眼科画像処理装置、および眼科画像処理プログラム | |
CN117529268A (zh) | Oct信号处理方法、oct装置及程序 | |
US20210049742A1 (en) | Image processing apparatus, image processing method, and non-transitory computer-readable storage medium | |
JP7264254B2 (ja) | 画像処理方法、画像処理装置、およびプログラム | |
EP3375349B1 (en) | Information processing apparatus, image generation method, and computer-readable medium | |
JP2019033956A (ja) | 計測装置 | |
JP2022089086A (ja) | 画像処理方法、画像処理装置、及び画像処理プログラム | |
WO2023182011A1 (ja) | 画像処理方法、画像処理装置、眼科装置、及びプログラム | |
WO2023282339A1 (ja) | 画像処理方法、画像処理プログラム、画像処理装置及び眼科装置 | |
WO2022177028A1 (ja) | 画像処理方法、画像処理装置、及びプログラム | |
WO2022250048A1 (ja) | 画像処理方法、画像処理装置、及びプログラム | |
US20230380680A1 (en) | Ophthalmic apparatus and method of controlling the same | |
WO2023199847A1 (ja) | 画像処理方法、画像処理装置、及びプログラム | |
WO2021256132A1 (ja) | 眼科装置、眼科装置の制御方法、及びプログラム | |
JP6833455B2 (ja) | 眼科撮影装置及びその制御方法、眼科撮影システム、並びに、プログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |