CN102458221B - 图像处理设备及其控制方法 - Google Patents
图像处理设备及其控制方法 Download PDFInfo
- Publication number
- CN102458221B CN102458221B CN201080024458.0A CN201080024458A CN102458221B CN 102458221 B CN102458221 B CN 102458221B CN 201080024458 A CN201080024458 A CN 201080024458A CN 102458221 B CN102458221 B CN 102458221B
- Authority
- CN
- China
- Prior art keywords
- image
- layer
- region
- artifact region
- brightness
- 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.)
- Expired - Fee Related
Links
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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0073—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by tomography, i.e. reconstruction of 3D images from 2D projections
-
- 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/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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0062—Arrangements for scanning
- A61B5/0066—Optical coherence imaging
-
- 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/02083—Interferometers characterised by particular signal processing and presentation
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T1/00—General purpose image data processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- General Physics & Mathematics (AREA)
- Ophthalmology & Optometry (AREA)
- Pathology (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Eye Examination Apparatus (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明提供一种图像处理设备,所述图像处理设备处理通过断层图像摄像设备拍摄被检眼的图像所获得的断层图像的图像,包括:层候选检测部件,用于从所述断层图像检测所述被检眼的视网膜的层候选;伪影区域判断部件,用于基于使用所述层候选所获得的图像特征,判断所述断层图像中的伪影区域;以及图像校正部件,用于基于所述伪影区域判断部件的判断结果和所述区域中的图像特征,校正所述伪影区域中的亮度。
Description
技术领域
本发明涉及一种图像处理设备及其控制方法和计算机程序。
背景技术
为了早期诊断作为生活方式病和致盲的前几位原因的各种疾病,普遍进行眼科检查。由于诸如OCT(光学相干断层成像)等的针对眼部的断层图像摄像设备允许三维观察视网膜层内部的状态,所以期望其对于疾病有效给出更恰当的诊断。通过从断层图像测量例如神经纤维层或视网膜的层厚度的变化、以及诸如视网膜色素上皮的不平坦性等的层几何形状的变化,可以定量诊断诸如青光眼、黄斑水肿和年龄相关性黄斑变性等的疾病的进展程度、以及治疗之后的康复水平。为了定量测量这些层的厚度,提出了一种用于使用计算机从断层图像检测视网膜的各个层、并测量这些层的厚度的技术(参考日本特开2008-073099号公报)。
另一方面,在OCT断层图像中,当物体强烈反射或吸收测量光时,通常在物体后方生成由信号的衰减或缺失所导致的伪影。注意,物体包括例如诸如血管等的组织和诸如白斑和出血等的病变部位。如图1A所示,当与视网膜的深度方向(以下称为z轴方向或A扫描方向)相关联地正常显示亮度时,最大亮度出现在视网膜色素上皮2附近。然而,如图1B所示,当在视网膜血管4的z轴的正向侧生成伪影区域5时,伪影区域5中在视网膜色素上皮6附近的亮度衰减或缺失。因此,在生成伪影的区域中,根据层的亮度的衰减程度来提取层、以及测量层厚度和层几何 形状通常变得困难。为了解决这一问题,提出了这样一种技术,该技术从眼底的表面图像提取血管区域,将血管区域逆投影至OCT断层图像,并且在逆投影区域附近对层进行插值,以估计由血管导致的伪影区域中的层位置(参考日本特开2007-325831号公报)。
发明内容
然而,日本特开2008-073099号公报所述的技术没有讨论用于计算生成伪影的区域中的层位置的任何方法。
另外,日本特开2007-325831号公报对作为很可能生成伪影的区域逆投影至断层图像上的区域附近的层进行简单插值,但是未通过检测逆投影区域中衰减的层的亮度信号来计算本来的层位置。此外,即使在为了层检测而执行图像校正时,由于伪影区域和其它区域具有不同的诸如直方图和对比度水平等的与亮度有关的性质,所以需要判断伪影区域以适应性地应用图像校正。
因此,本发明判断亮度由于诸如血管等的组织或者诸如白斑或出血等的病变部位的影响而衰减的区域,并且应用图像校正以便于检测该区域中的层。
本发明的实施例的一个方面,涉及一种图像处理设备,其对通过断层图像摄像设备拍摄被检眼的图像所获得的断层图像进行处理,所述图像处理设备包括:层候选检测部件,用于从所述断层图像检测所述被检眼的视网膜的层候选;伪影区域判断部件,用于基于使用所述层候选所获得的图像特征,判断所述断层图像中的伪影区域;以及图像校正部件,用于基于所述伪影区域判断部件的判断结果和所述伪影区域中的图像特征,校正所述伪影区域中的亮度。
本发明的实施例的另一方面,涉及一种图像处理设备,其对通过断层图像摄像设备拍摄被检眼的图像所获得的断层图像进行处理,所述图像处理设备包括:层候选检测部件,用于从所述断层图像检测所述被检眼的视网膜的层候选;伪影区域判断部件,用于基于使用所述层候选所获得的图像特征,判断所述断层图像中的伪影区域;以及层确定部件,用于基于所述伪影区域判断部件的判断结果,确定所述伪影区域中的所述视网膜的层的位置;其中,所述层确定部件使用用于指定所述伪影区域中所包括的层几何形状的并且通过多个控制点所定义的几何模型,并且基于与所述几何模型的几何形状相关联的评价函数和与所述控制点附近的亮度相关联的评价函数来确定所述层的位置。
通过以下(参考附图)对典型实施例的说明,本发明的其它特征将显而易见。
附图说明
图1A和1B是用于说明由于存在视网膜血管而导致视网膜色素上皮上的图像信号衰减的例子的图;
图2是示出根据第一实施例与图像处理设备10连接的设备的结构的框图;
图3是示出根据第一实施例的图像处理设备10的功能结构的框图;
图4是示出根据第一实施例的图像处理设备10的处理序列的流程图;
图5是用于说明根据第一实施例的伪影区域的判断方法的图;
图6是示出根据第一实施例判断伪影区域所需的图像处理 序列的流程图;
图7是示出根据第一实施例对被判断为伪影区域的区域的图像处理序列的流程图;
图8是示出根据第二实施例的图像处理设备10的功能结构的框图;
图9是示出根据第二实施例对被判断为伪影区域的区域的图像处理序列的流程图;
图10是示出根据第三实施例的图像处理设备10的功能结构的框图;
图11是示出根据第三实施例的图像处理设备10的处理序列的流程图;
图12是用于说明根据第三实施例的伪影区域的判断方法的图;
图13是示出根据第四实施例的图像处理设备10的功能结构的框图;
图14是示出根据第五实施例与图像处理设备10连接的设备的结构的框图;
图15是示出根据第五实施例的图像处理设备10的功能结构的框图;
图16是示出根据第五实施例的图像处理设备10的处理序列的流程图;
图17A和17B是用于说明根据第五实施例的伪影区域的判断方法的图;
图18是示出根据第六实施例的图像处理设备10的功能结构的框图;
图19是示出根据第七实施例的图像处理设备10的功能结构的框图;
图20是示出根据第七实施例的图像处理设备10的处理序列的流程图;
图21是示出根据第七实施例对被判断为伪影区域的区域的图像处理序列的流程图;
图22是用于说明根据第七实施例的用于测量层的不平坦几何形状的方法的图;
图23是示出根据第八实施例的图像处理设备10的功能结构的框图;
图24是示出根据第八实施例对被判断为伪影区域的区域的图像处理序列的流程图;
图25是用于说明根据第八实施例的图像处理方法确定方法和层确定方法的图;
图26是示出根据第九实施例的图像处理设备10的功能结构的框图;
图27是示出根据第九实施例的图像处理设备10的处理序列的流程图;
图28是用于说明根据第九实施例的伪影区域判断方法的图;
图29是示出根据第十实施例的图像处理设备10的功能结构的框图;
图30是示出根据第十一实施例的图像处理设备10的功能结构的框图;
图31是示出根据第十一实施例的图像处理设备10的处理序列的流程图;
图32是示出根据第十二实施例的图像处理设备10的功能结构的框图;以及
图33是示出可以通过软件实现图像处理设备10的各个单元 的计算机的基本配置的框图。
具体实施方式
下面将参考附图详细说明根据本发明的图像处理设备和方法的优选实施例。然而,本发明的范围不局限于所示出的例子。
第一实施例
当要利用OCT设备拍摄被检眼(作为检查对象的眼)的图像时,如果视网膜上存在血管或白斑,则由于测量光的强度降低,所以在所获得的图像上视网膜色素上皮减弱,并且难以检测层。因此,在本实施例中,从被检眼的断层图像判断伪影区域,并且根据该区域中的统计量对该区域进行图像校正。
图2是示出根据本实施例与图像处理设备10连接的设备的结构的框图。如图2所示,经由光纤和例如USB或IEEE1394的接口将图像处理设备10连接到断层图像摄像设备20。断层图像摄像设备20经由基于例如 的局域网(LAN)30与数据服务器40连接。注意,图像处理设备10可以经由诸如因特网等的外部网络与这些设备连接。断层图像摄像设备20获得眼部的断层图像,并且包括例如时域OCT或傅立叶域OCT。断层图像摄像设备20响应于操作者(未示出)的操作拍摄被检眼(未示出)的断层图像。设备20将所获得的断层图像发送给图像处理设备10。数据服务器40保持被检眼的断层图像和图像特征量等。数据服务器40存储从断层图像摄像设备20输出的被检眼的断层图像,并且分析从图像处理设备10输出的结果。数据服务器40响应于来自图像处理设备10的请求,将与被检眼相关联的先前的数据发送给图像处理设备10。
注意,本实施例将说明要获取视网膜色素上皮候选的情况。然而,要获取的候选不局限于视网膜色素上皮的外边界2。可以 检测其它层边界(视网膜色素上皮的内边界(未示出)、光感受器的内节和外节之间的边界3、或者外界膜(未示出))。当代替黄斑部,针对视神经乳头进行层候选的获取、伪影区域的判断和图像校正时,可以使用已知部位识别方法从要处理的区域预先排除诸如乳头中心部(凹陷部)等的不存在层的区域。
下面将参考图3说明根据本实施例的图像处理设备10的功能结构。图3是图像处理设备10的功能框图。如图3所示,图像处理设备10包括断层图像获取单元310、存储单元320、图像处理单元330、显示单元340、结果输出单元350和指示获取单元360。此外,图像处理单元330包括层候选检测单元331、伪影区域判断单元332和图像校正单元333。下面将参考与本实施例的图像处理设备10所执行的实际处理序列相关联的、图4所示的流程图,说明构成图像处理设备10的各个块的功能。
在步骤S410,断层图像获取单元310请求断层图像摄像设备20发送断层图像,并且获取从断层图像摄像设备20所发送的断层图像。然后,单元310将所获取的信息发送给存储单元320。存储单元320存储断层图像。在步骤S420,层候选检测单元331从存储单元320获取断层图像,并且从断层图像检测内界膜1和视网膜色素上皮候选点序列{Pi}。然后单元331将这些结果输出给存储单元320。
由于视网膜层对于各个层具有不同亮度,所以在两个相邻层之间的边界处生成亮度值的对比(边缘)。因此,通过关注该对比提取层边界。可以使用用于提取包括这类对比的区域的各种方法。例如,将对比认为是边缘,并且可以通过检测该边缘来提取层位置。更具体地,通过对断层图像应用边缘检测滤波器来检测边缘成分,并且从玻璃体侧开始在眼底的深度方向上搜索边缘。然后,将第一个峰值位置检测为玻璃体和视网膜层 之间的边界,并且将最大峰值位置检测为视网膜色素上皮边界。
可以通过应用诸如Snakes或水平集方法等的可变形模型来检测层边界。在水平集方法的情况下,定义比要检测的区域的维数高一个维数的水平集函数,并且将要检测的层边界认为是零水平线。通过更新水平集函数来控制轮廓,从而检测层边界。另外,可以使用诸如GraphCut等的图形定理检测层边界。在这种情况下,设置与图像的各个像素相对应的节点和被称为“汇”(sink)和“源”(source)的端子,并且设置在节点之间连接的边缘(n-link)和在端子之间连接的边缘(t-link)。通过基于向这些边缘赋予权重所创建的图计算最小切割,检测层边界。
上述层位置提取方法可以三维应用于作为要处理的对象的三维(3D)断层图像,或者在将输入的3D断层图像当作2D断层图像集时,可以独立应用于各二维(2D)断层图像。注意,用于检测层边界的方法不局限于这些方法,并且可以使用任何其它方法,只要它们可以从眼部的断层图像检测层边界即可。
在步骤S430,伪影区域判断单元332基于在步骤S420所检测到的视网膜色素上皮的候选点序列{Pi}的连续性,判断是否在各层候选点附近生成了伪影(是否生成了伪影区域)。如果已判断了伪影区域,则单元332计算与伪影区域中的亮度相关联的统计量。此外,单元332将判断结果输出给存储单元320。后面将使用图6示出的流程图详细说明该步骤的伪影区域判断处理。
在步骤S440,伪影区域判断单元332根据在步骤S430所获得的判断结果分支处理。也就是说,对于判断为生成了伪影的层候选点,单元332发送信号以指示图像校正单元333执行预定处理(处理序列进入步骤S450)。另一方面,如果单元332判断为未生成伪影的除伪影区域以外的区域(以下称为真图像区域),则发送信号以指示显示单元340执行预定处理(处理序列进入步 骤S455)。
此外,在步骤S450,当在预定层的候选点附近生成了伪影时,图像处理单元330执行分析处理。后面将使用图7示出的流程图详细说明该步骤的处理。另一方面,在步骤S455,作为在预定层的候选点附近未生成伪影时的处理,显示单元340执行用于与真图像区域相关联地显示断层图像的正常图像显示处理。
在步骤S460,指示获取单元360从外部获取关于是否将与被检眼相关联的当前处理结果存储在数据服务器40中的指示。操作者经由例如键盘和鼠标(未示出)输入该指示。如果操作者指示存储当前结果,则处理进入步骤S470;否则,处理跳至步骤S480。此后,在步骤S470,结果输出单元350将检查日期和时间、识别被检眼所使用的信息、被检眼的断层图像和通过图像处理单元330所获得的分析结果相互关联、作为要储存的信息,并且将该信息发送给数据服务器40。
在步骤S480,指示获取单元360从外部获取关于是否结束图像处理设备10的断层图像分析处理的指示。操作者经由例如键盘和鼠标(未示出)输入该指示。如果获取用于结束该处理的指示,则图像处理设备10结束其处理。另一方面,如果获取用于继续该处理的指示,则处理返回到步骤S410以执行用于下一被检眼的处理(或者用于相同被检眼的再处理)。以此方式执行图像处理设备10的处理。
下面将使用图5说明在步骤S430所执行的伪影区域判断处理的序列。图5示出包括伪影区域的断层图像的例子,并且示出以虚线围绕的区域5作为伪影区域。作为这些生成伪影的区域的特征,已知有下面两个特征。注意,在图5示出的断层图像中,将与视网膜的深度方向相对应的断层图像的纵向定义为z轴,并且将与深度方向垂直的横向定义为x轴。z轴方向对应于A扫描 方向。
(1)这些区域5中的亮度的平均值、方差和最大值都小于真图像区域中的亮度的平均值、方差和最大值。
(2)当在包括这些区域5的区域中计算视网膜色素上皮候选点时,由于对除视网膜色素上皮以外的高亮度区域、例如视网膜血管区域的错误提取,容易生成如点Pi和Pi+1一样的不连续部分。
因此,本实施例使用这些特征如下判断各伪影区域。
(i)检测不连续的层候选点对。
(ii)检查被判断为不连续点的各层候选点对的伪影生成侧。
(iii)跟踪伪影生成侧的层候选点序列,直到发现下一不连续点为止。
(iv)计算在穿过所跟踪的层候选点的A扫描线上、该层候选点的深层侧的亮度的统计量(平均值、方差、最大值等)。注意,深层侧表示图5中的两个点中在z轴方向上具有较大坐标值的一个点。
(v)将以不连续的候选点作为端点、并且包括低亮度统计量的区域判断为伪影区域。此外,由于认为亮度统计量是反映各伪影区域中的亮度衰减的程度的统计量,所以,在步骤S450中确定伪影区域中的图像处理方法时使用该亮度统计量。
下面将参考图6示出的流程图说明用于基于视网膜色素上皮的候选点序列{Pi}进行伪影判断的实际方法。在这种情况下,将作为处理对象的3D断层图像当作2D断层图像集,并且向各2D断层图像应用下面的2D图像处理。
在步骤S610,针对所有层候选点计算相邻层候选点之间的连续性C。如下给出连续性C:
其中,i为层候选点数量,并且Si为通过在层候选点Pi和Pi+1之间进行插值所获得的曲线上的像素的亮度的统计量。另外,S是由整个层候选点序列{Pi}所定义的曲线上的像素的亮度的统计量,并且Ts是阈值。在这种情况下,使用所定义的曲线上的像素的亮度的平均值作为亮度的统计量Si和S。在公式(1)中,如果|Si-S|等于或小于阈值Ts,则判断为相邻点连接。
注意,亮度统计量不局限于上述量,并且可以使用其它统计量,例如最大值、方差、众数值和中值。可选地,可以基于这些统计量的组合来判断连续性。另外,作为S,可以使用针对各摄像设备或被检者预先计算出的、与预定层的亮度相关联的统计量,或者可以使用预先所设置的标准值。注意,连续性的判断中所使用的指标利用与连接层候选点的曲线上的亮度相关联的统计量。然而,本发明不局限于这些特定指标。例如,作为血管或白斑的错误检测的结果,在不连续的候选点附近生成与z轴方向平行的边缘(如图5中的附图标记7所示),并且作为捕获这类特征的指标,可以使用与A扫描方向相关联的多个亮度剖面图之间的差异程度。注意,亮度剖面图表示示出A扫描方向上的空间位置和这些位置处的亮度之间的关系的图形,如图1A或1B的右侧的图所示,并且相邻亮度剖面图之间的差异通常小。
例如,将与亮度剖面图之间的差异程度相关联的指标的实际计算方法描述为:
(1)对于由图5示出的四个点Pi、Pi+1′、Pi″和Pi+1″所包围的局部区域中的像素,计算x轴方向上的亮度的差D。注意,Pi″是具有与Pi相同的x坐标值、并且具有最大z坐标值的点。另外,Pi+1″是具有与Pi+1相同的x坐标值、并且具有最大z坐标值的点。
(2)通过将等于或大于阈值Td1的差D相加,计算总值∑D。随着边缘在更宽的范围清晰出现,值∑D取更大的值。
(3)当总值∑D大于阈值Td2时,判断为不连续点。
还可以通过组合这多个指标来判断连续性。
在步骤S620,检查被判断为不连续点的层候选点对的伪影生成侧,以指定伪影区域的边缘部分。对各个不连续层候选点对执行伪影区域边缘部分指定处理。更具体地,在图5中的一对Pi和Pi+1的情况下:
(i)计算在穿过Pi的A扫描线上Pi的z轴的正向侧的亮度的统计量;
(ii)计算在穿过Pi+i的A扫描线上Pi+1′的z轴的正向侧的亮度的统计量。注意,Pi+1′是具有与Pi+1相同的x坐标和与Pi相同的z坐标的点;以及
(iii)比较这两个统计量,并且判断为在具有较小统计量的候选点侧生成伪影。将所确定的具有较小统计量的候选点定义为伪影区域的边缘部分。
在该例子中,由于Pi+1侧的亮度的统计量成为较小,所以判断为在Pi+1侧生成伪影。注意,用于选择计算亮度的统计量时所使用的点对的方法不局限于上述方法。例如,如图5所示,可以使用Pi′(具有与Pi相同的x坐标和与Pi+1相同的z坐标的点)和Pi+1 来进行判断。
在步骤S630,跟踪伪影区域侧的层候选点,直到发现下一个不连续点为止,以计算伪影区域的范围。例如,对于图5中的Pi+1,将直到Pi+3为止的区域判断为伪影区域。
在步骤S640,计算被判断为伪影区域的区域中的各个候选点的z轴正向侧的区域中亮度的平均值、方差或最大值。然而,计算与亮度相关联的统计量的空间范围不局限于这类特定范围。例如,可以将伪影区域分成任意的局部区域,并且可以针对各个局部区域计算统计量。在这种情况下,不必一定作为独立步骤来执行层候选点跟踪处理和与亮度信号相关联的统计量的计算处理,并且可以在每当针对任意范围跟踪层候选点时计算亮度信号的统计量。
利用上述处理,计算以不连续部分作为边的伪影区域的范围和伪影区域中亮度的统计量。注意,伪影区域的判断方法不局限于上述方法。例如,不仅可以对B扫描图像(与y轴垂直的断层图像)、而且还可以对与x轴垂直的断层图像执行伪影区域判断处理,并且这两个判断处理中所确定的伪影区域均可以被判断为伪影区域。可选地,可以将该判断处理三维应用于3D断层图像。
下面将参考图7说明在步骤S450所执行的处理的序列。
在步骤S710,图像校正单元333基于通过伪影区域判断单元332所计算出的与伪影区域中的亮度相关联的统计量,校正各伪影区域中的亮度。作为亮度校正方法,将说明基于伪影区域中的直方图转换的方法。更具体地,将伪影区域中的亮度平均值和方差调整成与真图像区域中的相同。也就是说,假定x是校正之前的信号,y是校正之后的信号,Sf和Af分别是伪影区域中的亮度的标准差和平均值,并且St和At分别是除伪影区域以外 的整个图像中的亮度的标准差和平均值,如下进行校正:
y=(St/Sf)*(x-Af)+At
注意,图像校正方法不局限于此。例如,可以通过下面的(i)~(iii)所述的方法实现图像校正方法。
可选地,可以使用任意图像校正,只要其可以使伪影区域中校正之前的亮度和校正之后的亮度之间的递增函数关系成立。
(i)例如,对亮度进行线性转换以使得伪影区域中的最大亮度与真图像区域中的最大亮度相一致。在这种情况下,假定y是图像校正之后伪影区域的亮度,x为图像校正之前伪影区域的亮度,ImaxF是伪影区域中的最大亮度,IminF为该区域中的最小亮度,并且ImaxT为真图像区域中的最大亮度,如下进行校正:
y=ImaxT*(x-IminF)/(ImaxF-IminF)
(ii)例如,执行诸如Sobel滤波或Laplacian滤波等的边缘强调处理或者仅允许高频成分通过的空间频率滤波处理。
(iii)执行强调层结构的基于Hessian矩阵的特征值的层结构强调滤波处理。该滤波基于Hessian矩阵的三个特征值(λ1、λ2、λ3)之间的关系来强调3D浓度分布的二次局部结构。Hessian矩阵是由图像的浓度值I的二阶偏导数所形成的矩形矩阵,如下给出:
Hessian矩阵的特征值之间的关系描述如下:
λ3≤λ2≤λ1 ...(4)
用于强调层结构所需的特征值的条件公式描述如下:
λ3<<λ2=λ1=0 ...(5)
可以通过根据基于这些公式所计算出的这三个特征值进行计算来强调视网膜的层结构:
其中,ω(λs;λt)是如下给出的权重函数:
其中,γ和α是权重。
注意,不必一定独立执行上述图像处理方法,而且还可以组合执行它们。当在步骤S640中将伪影区域分成多个局部区域并且对各个区域计算与亮度相关联的统计量时,还对各个局部区域应用图像校正。
在步骤S720,显示单元340将步骤S710中图像的校正结果重叠在断层图像上。当用线表示各伪影区域的边界时,可以针对各个边界使用预定颜色的线,或者可以在不明确表示边界的情况下用半透明颜色表示层。对于使用例如GUI所指定的区域,可以选择性显示校正前后的图像,并且可以显示在步骤S640所计算出的与各伪影区域中的亮度相关联的统计量的信息。如上所述,执行步骤S450中的处理。
根据上述结构,图像处理设备10指定伪影区域,并且基于例如与该区域中的亮度相关联的统计量执行图像校正,从而获得可以容易地检测到伪影区域中所包括的层的图像。
第二实施例
在本实施例中,不仅如第一实施例中一样在伪影区域中进行图像校正,而且还从校正后的图像检测预定层。本实施例甚至对于生成伪影和亮度衰减的区域应用图像校正,以使得便于检测其余边缘信息,并且允许通过检测该边缘信息计算更精确的层位置。
由于与根据本实施例的图像处理设备10连接的设备与第一实施例中的相同,因此不再重复对其的说明。图8是根据本实施例的图像处理设备10的功能框图。参考图8,本实施例的图像处理单元801与第一实施例的图像处理设备10的图像处理单元330的结构的不同在于添加了层确定单元334。
下面将说明本实施例的图像处理序列。本实施例的处理序列除步骤S450和S455外与图4的流程图所示的序列相同。因此,在本实施例中,仅说明步骤S450和S455,并且不再重复其它步骤的说明。
在步骤S450,执行伪影区域中的图像处理。下面将使用图9详细说明该步骤的处理。
在步骤S1010,图像校正单元333基于与通过伪影区域判断单元332所计算出的伪影区域的亮度相关联的统计量,校正该区域中的亮度。注意,该处理与步骤S710中的图像校正处理相同,并且不给出对其的详细说明。在步骤S1020,层确定单元334基于经过了图像校正单元333的图像校正的区域的亮度信息,获取要提取的层的图像特征,并且连接这些特征点作为层位置。
例如,视网膜色素上皮本来是各A扫描线上的最高亮度区 域,并且即使在伪影区域中,也趋向具有较高亮度。因此,通过在x轴方向上连接如下像素以确定层位置,其中,该像素是在图像校正后的区域的各个A扫描线上位于层候选点的z轴正向侧的、亮度最大的像素。然而,用于根据图像校正结果检测层的方法不局限于这一特定方法。
例如,可以通过在x轴方向上连接如下的像素来提取层,其中,该像素是在各个A扫描线上位于层候选点的z轴正向侧的、具有等于或大于预定值的亮度并且具有最大z坐标的像素。可选地,可以通过计算伪影区域的各个A扫描线上各像素在图像校正前后亮度的线性和、并且在x轴方向上连接与最大和相对应的像素来提取层。
可选地,在伪影区域的各个A扫描线上选择视网膜色素上皮的多个候选点,并且将所有层候选点定义为伪影区域中的层候选点集。对于从层候选点集所获得的层候选点序列的各个组合,可以设置与以下相关联的评价函数:
(i)层候选点序列的亮度和的大小;以及
(ii)层候选点序列几何形状的平滑度,
并且可以将使这两个评价函数值的线性和最大的层候选点组合确定为层位置。注意,在选择层候选点的组合时,可以使用不仅包括伪影区域中的层候选点序列、而且还包括该区域附近的层候选点的层候选点序列来计算评价函数值。
在步骤S1030,通过计算所计算出的与视网膜色素上皮相对应的层候选点序列与步骤S420所计算出的内界膜1之间的各个x坐标和y坐标的距离,测量视网膜层厚度。然而,测量内容不局限于此。例如,可以计算层候选点之间的角度分布,以检查层几何形状的不平坦性。要测量的层厚度不局限于视网膜层厚度。例如,可以分析诸如光感受器细胞层等的其它层几何形 状。将包括所计算出的层厚度和层几何形状的信息输出给存储单元320。
在步骤S1040,显示单元340在断层图像上重叠:
(i)层确定结果;
(ii)层几何形状测量结果;以及
(iii)伪影区域的范围和该区域中的图像校正结果。
对于(i),显示单元340将步骤S1020的层确定结果重叠在断层图像上。当以线表示层的边界时,可以对于各边界使用预定颜色的线,或者可以在不明确表示边界的情况下以半透明颜色表示层的区域。注意,在进行这类显示时,希望采用允许使用例如GUI来选择感兴趣的断面的配置。另外,可以使用已知的体绘制技术三维显示这些结果。
对于(ii),显示单元340显示整个3D断层图像(x-y平面)的层厚度的分布映射图,作为层几何形状测量结果。然而,本发明不局限于这一特定显示方法,并且可以显示感兴趣的断面中的各个层的面积,或者可以显示整个预定层的体积。可选地,可以计算并显示操作者在x-y平面上所指定的区域的体积。
对于(iii),当在步骤S1020所检测的层的可靠性低时(例如,当所检测到的层的信号弱时),未检测到层,并且将图像校正结果重叠在断层图像上。然而,本发明不局限于这一图像校正结果的特定显示方法,并且即使在层检测结果满意时,也可以将图像校正结果显示在显示单元340上。
如上所述,执行步骤S450的伪影区域图像处理。下面将说明步骤S455中的真图像区域图像处理。在步骤S455,作为真图像区域中的处理,
(i)将在步骤S420所获取的层候选点序列定义为层位置,并且根据层位置测量层几何形状;以及
(ii)将层位置和层几何形状的测量结果显示在显示单元340上。
因而,下面将详细说明该步骤的处理。
以与图9的步骤S1030相同的方式测量层几何形状。然而,在真图像区域图像处理的情况下,分析对象是不同于步骤S1030的未校正的图像。然后,如图9的步骤S1040一样,显示单元340显示结果。然而,在真图像区域图像处理的情况下,不同于步骤S1040,不显示图像校正结果。
根据上述结果,图像处理设备10对指定伪影区域执行图像校正,并且从校正结果检测与层位置相对应的图像特征,从而更精确地计算伪影区域中的层位置。
第三实施例
在本实施例中,代替第一和第二实施例中仅使用断层图像对伪影区域的判断,根据被检眼的断层图像生成投影图像,并且将从投影图像所提取的组织或病变部位的位置信息逆投影至断层图像,以预先缩窄伪影候选区域。通常,与仅根据断层图像相比,根据投影图像来计算由例如血管(或出血)所导致的伪影区域的位置信息更为容易。因此,本实施例将说明下面的情况:从投影图像提取血管(出血)区域,将该位置信息映射在断层图像上,并且围绕映射后的区域搜索并指定伪影区域的边缘部分,以更高精度计算伪影区域的范围。
由于与根据本实施例的图像处理设备10连接的设备的结构与第一实施例中的相同,所有不再重复对其的说明。图10是根据本实施例的图像处理设备10的功能框图。参考图10,不同于第一实施例的图像处理设备10的图像处理单元330,图像处理单元1001包括投影图像生成单元335和特征提取单元336。由于其余单元与图3中的相同,所有不再重复对其的说明。
下面将参考图11所示的流程图说明根据本实施例的图像处理的内容。注意,该流程图中的大多数步骤与图4的流程图中的相同,并且相同步骤编号表示相同处理。不再重复对相同处理的说明。在本实施例中,不同于第一实施例,在步骤S420的层候选检测处理和步骤S430的伪影区域判断处理之间执行投影图像生成处理S1110和特征提取处理S1120。
在步骤S1110,投影图像生成单元335通过将断层图像投影生成图像。更具体地,利用作为简单在z轴的正向上将断层图像上的像素的亮度相加所获得的值的像素值来定义投影图像。然而,投影图像的各像素值不局限于这类值,并且可以将亮度和除以相加的像素的数量。可选地,可以使用各深度位置处的亮度的最大值或最小值作为投影图像的各像素值。另外,不必将z轴方向上的所有像素的亮度相加,并且可以仅将任意范围中的或特定层之间的像素的亮度相加。
在步骤S1120,从通过投影图像生成单元335所生成的投影图像来提取特征区域,其中,在该特征区域中,存在诸如被检眼的视网膜血管等的生物组织或病变部位。由于视网膜血管具有细的线状结构,所以使用强调线状结构的滤波器来进行提取。在这种情况下,使用基于对比度的线段强调滤波器,例如,计算被定义为结构要素的线段中的图像浓度值的平均值和围绕该结构元素的局部区域中的平均值之差的滤波器。注意,可以使用作为滤波器的处理结果所获得的多值区域作为血管提取结果,或者可以使用利用特定阈值二值化的区域作为提取结果。
然而,用于强调线状结构的方法不局限于此。例如,可以使用诸如Sobel滤波器或Laplacian滤波器等的差分滤波器。可以针对浓度值图像的各个像素计算Hessian矩阵的特征值,并且可以从作为结果所获得的两个特征值的组合来提取线状区域。此 外,可以使用诸如仅用线段作为结构元素的顶帽(top hat)运算等的任意已知的血管提取方法。
当将这样在步骤S1120所计算出的投影图像上的特征区域(x,y)逆投影至断层图像时,获得逆投影区域,如图12中的虚线区域8所示。通常,在视网膜血管的z轴的正向侧容易发生亮度衰减。因此,当将所提取的特征(在x-y方向上)的位置逆投影至断层图像时,逆投影得到的虚线区域8更可能包括伪影区域5。然而,当对错误提取的区域进行逆投影时,在逆投影区域中不发生亮度衰减。即使当对正确提取的视网膜血管进行逆投影时,逆投影区域下方的亮度衰减轻微,并且在一些情况下,对层提取几乎没有影响。
因此,判断在逆投影区域中和该区域的边界附近是否生成了伪影。如果生成了伪影区域,则计算与该区域中的亮度相关联的统计量。因此,步骤S430中的伪影区域判断方法除作为连续性的计算对象的层候选点的范围以外、与第一实施例的步骤S610~S640中的方法基本相同。更具体地,不是针对所有层候选点、而是针对逆投影区域内部及x-y方向上的该区域附近执行连续性计算处理。
根据上述结构,本实施例的图像处理设备10根据断层图像和投影图像指定伪影区域,并且基于例如该区域中的亮度统计量执行图像校正,从而获得可以更容易地检测到伪影区域中所包括的层区域的图像。
第四实施例
不同于第三实施例,本实施例在判断伪影区域之后不仅执行伪影区域中的图像校正,而且还从校正后的图像检测预定层。本实施例应对下面几点。
(i)当由于血管(出血)而生成伪影时,将根据投影图像所计 算出的血管(出血)区域的位置信息映射至断层图像,并且从映射得到的区域的周围区域搜索并指定伪影区域的边缘部分,从而允许以更高精度计算伪影区域的范围。
(ii)即使在生成伪影并且亮度衰减的区域中,也使该区域经过图像校正以允许容易地检测其余边缘信息,并且从校正后的区域检测边缘信息,从而允许计算更精确的层位置。
由于与根据本实施例的图像处理设备10连接的设备的结构与第三实施例中的结构相同,所以不再重复对其的说明。图13是根据本实施例的图像处理设备10的功能框图。不同于第三实施例的图像处理单元1001,本实施例的图像处理单元1301包括层确定单元334。本实施例的图像处理的内容除步骤S450和S455的处理以外与图13中的相同。因此,下面仅说明本实施例的步骤S450和S455的处理,并且不给出其它步骤的说明。
在步骤S450,作为伪影区域中的图像处理,执行图像校正、层确定、层几何形状测量和结果显示处理。该步骤的处理与第二实施例中的步骤S1010~S1040中的相同,并且不再对其进行详细说明。在步骤S455,作为未生成伪影时所执行的处理,根据在步骤S420所获取的层位置测量层几何形状,并且将层位置和层几何形状测量结果重叠在断层图像上。该重叠方法的细节与第二实施例中的步骤S1110~S1120中的相同,并且不再对其进行详细说明。
根据上述结构,本实施例的图像处理设备10根据断层图像和投影图像指定伪影区域,并且在该区域中执行图像校正。由于该设备根据校正结果检测层的图像特征,所以可以更精确地计算该区域中的层位置。
第五实施例
在本实施例中,在第三实施例以外,将从被检眼的表面图 像和投影图像中的至少一者所提取的组织或病变部位的位置信息逆投影至断层图像上,以预先缩窄伪影候选区域。这是因为,当生成诸如白斑等的特别容易从表面图像提取的病变部位时,使用表面图像计算白斑区域,并且从所提取的病变部位的周围区域搜索并指定伪影区域的边缘部分,从而允许以更高精度计算伪影区域的范围。
图14示出与根据本实施例的图像处理设备10连接的设备的结构。在本实施例中,不同于第三实施例,除断层图像摄像设备20以外,该结构还包括表面图像摄像设备50。表面图像摄像设备50拍摄眼部的表面图像,并且包括例如眼底照相机或SLO(激光扫描检眼镜)。图15是本实施例的图像处理设备10的功能框图。本实施例的图像处理设备10包括表面图像获取单元315,并且不同于第三实施例的图像处理设备10的结构,图像处理单元1501包括对准单元337。
下面将参考图16所示的流程图说明本实施例的图像处理单元1501的图像处理。本实施例的图像处理序列除步骤S1610~S1650的处理外与图11所示的几乎相同。因此,下面将说明这些步骤的处理。
在步骤S1610,除通过断层图像获取单元310获取断层图像以外,表面图像获取单元315请求表面图像摄像设备50发送表面图像,并且获取从表面图像摄像设备50所发送的表面图像。假定输入眼底照相机图像作为表面图像。单元315将所获取的信息发送给存储单元320。
在步骤S1110的投影图像生成处理之后的步骤S1620,特征提取单元336从通过表面图像获取单元315所获取的表面图像提取诸如血管等的组织或者诸如白斑等的病变区域。由于视网膜血管具有线状结构,所以使用强调线状结构的滤波器进行提取。 由于线状结构提取方法与步骤S1120中的相同,所以不给出对其的说明。由于白斑作为颗粒状高亮度区域存在,所以通过诸如顶帽变换等的形态学运算来计算白斑。在这种情况下,通过形态学运算获得作为高亮度多值区域的白斑区域,并且可以使用该多值区域本身作为提取结果,或者可以使用利用特定阈值二值化的区域作为提取结果。然而,白斑提取方法不局限于此,并且可以使用表面图像的亮度和强调对比度的已知滤波器的输出结果作为特征量,通过诸如支持向量机等的识别器或者诸如Ada Boost等的识别器组来识别白斑。
在步骤S1630,对准单元337在投影图像和表面图像之间进行对准,以将投影图像的坐标与表面图像的坐标相关联。在对准时,预先定义表示两个图像之间的相似度的评价函数,并且将图像变形以获得最佳评价值。作为相似度评价方法,使用用于利用相互信息内容基于像素值来评价相似度的方法。然而,本发明不局限于这一特定方法,并且可以使用通过特征提取单元336从表面图像和投影图像所提取的均方差、相关系数、或者血管区域之间的重叠面积或血管分支部之间的距离等。在仿射变换的前提下,通过变换或转动图像、或者改变放大因数来实现图像变形。
在步骤S1640,使用在步骤S1630所计算出的对准参数,将在步骤S1620计算出的来自表面图像的血管或白斑提取结果逆投影至断层图像。作为逆投影区域,获得图17A中的虚线区域8所示的区域。通常,由于亮度衰减常常发生在视网膜血管的z轴的正向侧,所以当将各特征提取区域(在x-y方向上)的坐标和亮度逆投影至断层图像时,各逆投影虚线区域8更可能包括伪影。然而,当如本实施例中一样,从表面图像提取白斑时,在白斑候选区域中可能错误地提取诸如脉络膜小疣等的颗粒状高 亮度区域。在这种情况下不发生层区域中的亮度衰减,如图17B所示。即使当对正确提取的视网膜血管区域进行逆投影时,逆投影区域下方的亮度衰减轻微,并且在一些情况下,对层提取几乎没有影响。因此,判断在逆投影区域中和该区域边界附近是否生成了伪影。如果生成了伪影,则计算该区域中的亮度的统计量。
该伪影区域判断方法与第一实施例的步骤S610~S640的情况下的判断方法大体相同,但是作为连续性的计算对象的层候选点的范围不同于第一实施例。更具体地,不是对所有层候选点、而是对逆投影区域内部和x-y方向上的该区域附近执行连续性计算处理。注意,除仅从断层图像所获得的信息以外,还可以参考从投影图像和眼底图像所获得的信息来判断伪影区域。例如,当从投影图像所获得的视网膜血管区域与从眼底图像所获得的视网膜血管区域重叠时,可以认为该区域很可能是血管,并且可以判断为该区域的边缘部分不连续。可选地,可以计算根据断层图像所计算出的连续性值与血管区域重叠程度的值的线性和,并且可以使用阈值将该线性和二值化,以判断连续性。
由于步骤S1650的伪影区域图像处理显示图像校正之后的校正结果,也就是说,采用与第一实施例中的大体相同的序列,所以不给出对其的详细说明。然而,在本实施例中,在伪影区域的图像校正时,还可以参考从眼底图像所获得的信息。例如,当在眼底图像上白斑的亮度非常高时,由于即使在断层图像上、亮度在层候选点的z轴的正向侧也很可能衰减,所以与白斑区域的亮度信号值成比例地放大或强调亮度。在这种情况下,作为白斑区域的亮度,直接参考眼底图像上该区域的像素值。然而,白斑区域的亮度不局限于这些值,并且可以参考通过例如形态 学运算等所获得的处理结果的值(多值数据)。
根据上述结构,本实施例的图像处理设备10基于使用表面图像和投影图像所指定的伪影区域中的亮度的统计量执行图像校正,从而获得可以更容易地检测该区域中所存在的层区域的图像。
第六实施例
本实施例不仅执行第五实施例中的伪影区域的图像校正,而且还从校正后的图像检测预定层。尤其在由于白斑而生成伪影时,本实施例利用以下几点。
(i)将根据表面图像所计算出的白斑区域的位置信息映射至断层图像,并且从该区域的周围区域搜索并指定伪影区域的边缘部分,从而以更高精度计算伪影区域的范围。
(ii)即使在生成伪影并且亮度衰减的区域中,也使该区域经过图像校正以允许容易地检测其余边缘信息,并且从校正后的区域检测边缘信息,从而允许计算更精确的层位置。
与根据本实施例的图像处理设备10连接的设备的结构与第五实施例中的相同。图18是根据本实施例的图像处理设备10的功能框图。不同于第五实施例的图像处理单元1501,本实施例的图像处理单元1801包括层确定单元334。除步骤S1650和S455的处理以外,本实施例的图像处理序列与图16所示的相同。因此,仅说明步骤S1650和S455的处理,并且不再重复对其它步骤的说明。
在步骤S1650,作为伪影区域中的图像处理,执行图像校正、层确定、层几何形状测量和结果显示处理。该步骤的处理与第二实施例的步骤S1010~S1040的处理相同,并且不再重复对其的详细说明。然而,在本实施例中,还可以使用在步骤S1010从眼底图像所获得的信息来执行图像校正。用于在图像 校正时参考从眼底图像所获得的信息的实际序列与第五实施例步骤S1650的情况下的序列相同,并且不再重复对其的说明。
在步骤S455,作为在没有生成伪影时所执行的处理,根据在步骤S420所获取的层位置测量层几何形状,并且将层位置和层几何形状测量结果重叠在断层图像上。这种情况下的重叠方法与第二实施例中的相同,并且不再重复对其的详细说明。
根据上述结构,图像处理设备10在根据表面图像和投影图像所指定的伪影区域中执行图像校正。通过从校正结果检测层的图像特征,可以更精确地计算该区域中的层位置。
第七实施例
在本实施例中,根据被检眼的断层图像判断伪影区域,并且考虑到伪影区域中的亮度和该区域周围的层几何形状,使用这两者的信息来计算伪影区域中的层位置。
由于与根据本实施例的图像处理设备10连接的设备的结构与第一实施例的图2所示的相同,所以没有给出对其的说明。图19是根据本实施例的图像处理设备10的功能框图。参考图19,代替图像校正单元333,本实施例的图像处理单元1901包括图像处理方法确定单元1910,并且不同于第一实施例的图像处理设备10的图像处理单元330的结构,还包括层确定单元334。图像处理方法确定单元1910包括亮度使用判断单元1911和评价函数设置单元1912。下面将与通过本实施例的图像处理设备10所执行的实际处理序列相关联,参考图20所示的流程图说明构成图像处理设备10的各个块的功能。
下面将说明本实施例的图像处理序列。注意,除步骤S2010~S2030的处理以外,本实施例的处理序列与图4所示的流程图中的相同。因此,下面仅说明步骤S2010~S2030,并且不再重复对其它步骤的说明。
在步骤S2010,伪影区域判断单元332根据在步骤S430所获得的判断结果分支处理。也就是说,对于被判断为生成了伪影的层候选点,单元332发送信号以指示图像处理方法确定单元1910执行预定处理。另一方面,如果单元332判断为未生成伪影的真图像区域,则发送信号以指示显示单元340执行预定处理。
在步骤S2020,当在预定层的候选点附近生成伪影时,图像处理单元1901执行分析处理。后面将使用图21所示的流程图详细说明该步骤的处理。在步骤S2030,显示单元340将层判断结果重叠在断层图像上。当以线表示层的边界时,可以对各个边界使用预定颜色的线,或者在不明确表示边界的情况下,以半透明颜色表示层的区域。注意,在进行这一显示时,希望采用允许使用例如GUI来选择感兴趣的断面的结构。另外,可以使用已知的体绘制技术三维显示这些结果。
此外,可以通过针对各个坐标(x,y)计算所计算出的与视网膜色素上皮相对应的层候选点序列与步骤S420中所计算出的内界膜1之间的距离,测量视网膜层厚度。在这种情况下,显示单元340呈现与所测量的层几何形状相关联的信息作为对于整个3D断层图像(x-y平面)的层厚度的分布映射图。然而,本发明不局限于这一特定显示方法。例如,显示单元340可以与检测结果的显示处理同步地显示感兴趣的断面中各个层的面积。可选地,显示单元340可以显示整个体积,或者可以计算并显示操作者在x-y平面上所指定的区域的体积。
下面将说明本实施例的步骤S2020中的伪影区域图像处理。在本实施例中,向层位置应用可变形模型,从而使得即使在包括噪声的断层图像中也可以计算层位置。下面将说明使用Snakes作为可变形模型的例子。在这种情况下,通过使与模型几何形状相关联的评价函数值以及与构成该模型的控制点附近 的亮度相关联的评价函数值的线性和最小化,确定层位置。
作为与几何形状相关联的评价函数,使用构成模型的控制点位置的差与二阶微分值的线性和。模型几何形状随着线性和的减小而变得更平滑。作为与亮度相关联的评价函数,使用通过向构成模型的控制点附近的亮度梯度分配负号所获得的值。这是为了随着与边缘的距离的减小而减小评价函数值。通常将为了使可变形模型变形所使用的评价函数的权重设置成固定值,而不管构成模型的控制点是否包括在伪影区域中。在这种情况下,由于伪影区域中的亮度衰减,并且该区域中的亮度变化小,所以实际上基于与模型几何形状相关联的评价函数值的大小来确定层位置。当在伪影区域中残留例如边缘的信息时,在重视与亮度相关联的信息的情况下确定层位置(优先于基于模型几何形状的平滑度计算层位置的情况),从而更精确地检测层几何形状。为此,根据伪影区域中的控制点处的亮度的衰减程度,与真图像区域相比,将与亮度相关联的评价函数的权重设置得更大。
然而,当伪影区域中的亮度小,并且几乎没有残留例如边缘的信息时,在确定层位置时不能使用例如边缘的亮度信息。因此,不增大与亮度相关联的评价函数的权重,并且基于与模型几何形状相关联的函数评价值的大小来确定层位置。
下面将参考图21所示流程图说明用于实际设置可变形模型的各个评价函数的权重的方法。
在步骤S2110,图像处理方法确定单元1910从存储单元320读出在图6的流程图的步骤S640所计算出的、与伪影区域中的亮度相关联的统计量,其中,图6详细示出步骤S430中的伪影区域判断处理。单元1910基于该统计量判断:由于在该区域中残留(衰减的)边缘所以在层检测时使用亮度信息、还是由于亮 度不足而不使用例如边缘的亮度信息。
更具体地,计算指标E(i),并且如果E(i)=1,则判断为在确定层位置时使用亮度信息;如果E(i)=0,在判断为不使用亮度信息。如下给出指标E(i):
其中,i是可变形模型的控制点编号,B是与背景区域(例如,在内界膜1的z轴的负方向侧的区域)中的亮度相关联的统计量,并且Fi是与控制点i所属的伪影区域中的亮度相关联的统计量。对于与亮度相关联的统计量,使用最大值。另外,Ts是阈值。
注意,与亮度相关联的统计量不局限于此。例如,可以使用平均值、方差或标准差。另外,可以将伪影区域分成任意局部区域,并且可以使用各个控制点所属的局部区域中的亮度的统计量作为Fi。
在步骤S2120,图像处理方法确定单元1910获取与伪影区域周围的层候选点序列的不平坦性相关联的信息。这是因为,如果围绕伪影区域的层几何形状具有不平坦性,则在该区域中也很可能生成不平坦性,并且不得不通过降低与几何形状的平滑度的评价相关联的权重来计算层位置。作为用于计算层几何形状的实际不平坦性的方法,计算与伪影区域周围的层候选点之间的角度相关联的统计量。在这种情况下,使用最大值作为统计量。
如图22所示,计算层候选点i处的层候选点之间的角度,作为通过将线段Qi-1-Qi延长至Qi+1侧所获得的线段与线段Qi-Qi+1 之间的角度θi。不平坦性随着θi增大而变大。对于伪影区域周围的层候选点执行这一角度计算,并且使用所计算出的角度中的最大值作为表示层候选点序列的不平坦程度的指标。注意,与不平坦程度相关联的指标不局限于层候选点之间的角度,并且可以计算层候选点位置处的二阶微分值(Qi-1-2Qi+Qi+1)的统计量(平均值、方差、最大值等)。可选地,可以计算在将伪影区域周围的层候选当作曲线时的极值或拐点的数量。与层候选点的不平坦程度相关联的统计量不局限于最大值。例如,可以使用平均值、方差或标准差。
在步骤S2130,图像处理方法确定单元1910使用步骤S2110的判断结果和在步骤S2120所计算出的与层几何形状的不平坦程度相关联的指标,设置可变形模型的评价函数的权重。作为与几何形状相关联的评价函数的权重,设置与在步骤S2120所计算出的、表示层几何模型的不平坦程度的指标成反比的值。根据与在步骤S2110所计算出的伪影区域中的亮度信息的使用相关联的判断结果,如下设置与亮度相关联的评价函数的权重。
(i)当在确定层位置时使用伪影区域中的亮度信息时
根据伪影区域中的亮度的衰减的程度,增大与亮度相关联的评价函数的权重。然后,将与在步骤S640所计算出的区域中的亮度统计量Fs和真图像区域中的亮度统计量Ts之间的比Ts/Fs成比例的值设置为与亮度相关联的评价函数的权重。然而,与亮度相关联的评价函数的权重的设置方法不局限于此。例如,可以设置任意权重函数,只要在Fs和与亮度相关联的评价函数的权重之间建立递减函数的关系即可。
(ii)当在确定层位置时不使用伪影区域中的亮度信息时
在伪影区域中,将与几何模型的亮度相关联的评价函数的权重设置成与真图像区域相同的值。注意,(ii)中与亮度相关联 的评价函数的权重的设置方法不局限于此。例如,可以降低与亮度相关联的评价函数的权重,或者将其设置成0。
在步骤S2140,层确定单元334根据在步骤S2130所设置的评价函数的权重计算评价值,并且使用诸如Greedy算法等的优化方法进行迭代计算,从而最小化评价函数值。当评价值的变化量小于预定值、或者当迭代计算次数超过预定计数时,单元334结束层几何模型的变形,并且将在结束时的层几何模型的位置确定为层位置。
注意,可以作为2D曲线或3D曲面模型来计算这些几何模型。本实施例说明了使用Snakes作为可变形模型的例子,但是可以使用水平集。此外,可以使用任意模型,只要它们在模型的变形时参考亮度的基于模型的分割方法中设置与亮度相关联的评价函数的权重即可。
根据上述结构,图像处理设备10指定伪影区域,并且考虑该区域周围的层几何形状的不平坦性和该区域中的边缘信息来执行图像处理,因此以比传统方法更高的精度来计算层位置。
第八实施例
在本实施例中,在计算伪影区域中的层位置时,判断是否使用该区域中的例如边缘的亮度信息。当使用亮度信息时,校正该区域中的亮度,然后计算层位置。当不使用亮度信息时,通过插值计算层位置。本实施例覆盖以下几点。
(i)在亮度衰减但并非不足的血管或小白斑区域中,转换亮度以便于使用例如残留的边缘的信息,然后检测层位置,从而计算更精确的层位置。
(ii)在亮度不足并且不能使用边缘信息的大白斑或严重出血区域中,考虑伪影区域的生成位置和周围层几何形状,对周围层位置进行插值,从而计算更精确的层位置。
由于与根据本实施例的图像处理设备10连接的设备的结构与第七实施例中的相同,所以不再重复对其的说明。图23是根据本实施例的图像处理设备10的功能框图。在本实施例中,不同于第七实施例,图像处理单元2301包括图像校正单元333,并且图像处理方法确定单元1910包括图像校正方法设置单元1913和插值函数设置单元1914。
下面将参考图20示出的图像处理序列说明本实施例中的图像处理的内容。注意,除步骤S2020的处理以外,本实施例的图像处理序列与第七实施例的相同。因此,下面将参考图24说明步骤S2020的处理的改变点,并且不再重复对其它步骤的说明。
在步骤S2410进行判断:由于在伪影区域中残留弱的边缘所以在计算层位置时使用亮度信息、还是由于亮度不足而不使用例如边缘的亮度信息。由于实际判断序列与第七实施例的步骤S2110的相同,所以不再重复对其的说明。
在步骤S2420,亮度使用判断单元1911根据步骤S2410的判断结果分支处理。也就是说,如果判断为在确定层位置时使用伪影区域中的亮度信息,则单元1911发送指示图像校正方法设置单元1913执行预定处理的信号。另一方面,如果判断为不使用该区域中的亮度信息,则单元1911发送指示插值函数设置单元1914执行预定处理的信号。
在步骤S2430,设置在伪影区域中执行亮度的转换(图像校正)所需的参数。可以使用各种图像校正方法。在本实施例中,将说明下面的图像校正方法中的参数的设置序列。
(i)基于线性转换的方法
(ii)基于直方图转换的方法
(iii)强调层结构的方法
在执行(i)中基于线性转换的图像校正时,使用伪影区域中的最大亮度ImaxF、该区域中的最小亮度IminF、以及真图像区域中的最大亮度ImaxT,如下设置线性转换参数:
y=ImaxT*(x-IminF)/(ImaxF-IminF)
其中,y是图像校正之后的伪影区域中的亮度,并且x是图像校正之前的伪影区域中的亮度。在这种情况下,执行将伪影区域中的最大亮度调整成真图像区域中的最大亮度的图像校正。
在(ii)的方法中,执行各种直方图转换以使伪影区域中的直方图特征近似于真图像区域的直方图特征。在该方法中,将伪影区域中的亮度平均值和方差调整成与真图像区域中的相同。在这种情况下,使用伪影区域中的亮度的标准差Sf和平均值Af、以及除伪影区域以外的整个图像中的亮度的标准差St和平均值At,可以如下设置图像校正函数的参数:
y=(St/Sf)*(x-Af)+At
其中,x是校正之前的信号,并且y是校正之后的信号。
此外,在(iii)的方法中,当使用Hessian矩阵的特征值执行层结构强调处理时,通过上述公式(5)描述与特征值λ1、λ2和λ3(λ1≥λ2≥λ3)相关联的条件公式。Hessian矩阵是由多变量函数的全部二阶偏导数所形成的矩形矩阵,如公式(3)所示,并且I是图像的浓度值。根据这三个特征值,可以使用上述公式(6)强调层结构。注意,ω(λs;λt)是对于层结构强调处理所使用的权重,并且如上述公式(7)一样对其进行设置。然而,在公式(7)中,值s和t的实际组合为(s,t)=(1,3)或(2,3),并且将γ和α分别设置成固定值。
为了根据要检测的层的厚度实现层结构强调处理,如下设置基于作为层结构强调处理的预处理所执行的Gaussian函数的平滑参数s:
(i)根据伪影区域周围的层位置处的A扫描线上的亮度剖面图,计算层位置处的层厚度。在这种情况下,计算与层位置处的亮度的差等于或小于预定值的线的范围,并且使用其长度作为层厚度。
(ii)与伪影区域周围的层厚度的值成比例地设置在进行平滑时所使用的Gaussian滤波器的分辨率s的值。
注意,本实施例的图像校正方法不局限于这些方法。例如,可以使用任意图像校正方法,只要在伪影区域中的校正前后的亮度之间建立递增函数的关系、并且该方法包括可调整的参数即可。
在步骤S2440,图像校正单元333基于在步骤S2430所设置的图像校正方法,对伪影区域中的亮度执行转换(图像校正),从而使得便于检测层位置。此外,在步骤S2450,层确定单元334基于经过了图像校正单元333的图像校正的区域的亮度信息,获取要提取的层的图像特征,并且通过连接这些特征点来定义层位置。例如,视网膜色素上皮本来为各个A扫描线上的最高亮度区域,并且即使在伪影区域中也趋于具有较高亮度。因此,通过在x轴方向上连接在图像校正之后的区域中的各个A扫描线上位于层候选点的z轴的正向侧的、亮度最大的像素,确定层位置。
利用上述序列,确定在亮度使用判断单元1911判断为使用伪影区域中的亮度信息时的层位置。下面将说明在亮度使用判断单元1911判断为不使用伪影区域中的亮度信息时图像处理方法确定单元1910的处理内容。
在步骤S2460,图像处理方法确定单元1910获取通过伪影区域判断单元332所计算出的、与伪影区域的范围相关联的信息。更具体地,假定i为图25中的伪影区域的标志,单元1910 获取与以下相关联的信息:伪影的生成位置(xi,yi)和宽度Wi、以及属于该区域附近的真图像区域ni的层候选点的数量di。在步骤S2470,单元1910计算表示伪影区域周围所存在的层候选点序列的不平坦程度的指标。由于这一指标与步骤S2120所计算的相同,所以不再重复对其的详细说明。
此外,在步骤S2480,图像处理方法确定单元1910根据在步骤S2460和S2470所获取的信息,选择对伪影区域之间的层候选点序列插值时所使用的插值函数的类型或阶次、以及插值中所使用的层候选点。
(i)首先,在各伪影区域中,根据在步骤S2470所计算出的伪影区域的宽度Wi和与该区域附近的层几何形状相关联的指标的值(例如,层候选点之间的角度的统计量θi),选择插值函数的类型或阶次。例如,如果Wi小于预定值,则选择线性插值;如果Wi等于或大于预定值,则选择B样条插值。可选地,可以使用下面的选择方法。也就是说,如果层几何形状具有大的不平坦性(如果θi的平均值或最大值等于或大于预定值),则使用自然样条插值,其中,自然样条插值的插值曲线穿过控制点,并且与B样条插值相比,自然样条插值允许计算更精确的层位置。另外,可以使用下面的选择方法。也就是说,即使对于相同类型的插值函数,与伪影区域周围的层候选点之间的角度θi的统计量(平均值、最大值等)的大小成比例地设置插值函数的阶次。
(ii)然后,选择伪影区域中的层位置的插值所使用的层候选点。更具体地,检查属于伪影区域i附近的真图像区域ni-1和ni的层候选点的数量di-1和di是否满足为通过所选择的插值函数进行插值所需的最小必需数量。如果这些数量满足最小必需数量,则根据插值所需的数量来选择属于这些真图像区域的层候选点中在伪影区域附近的层候选点。另一方面,如果这些数量不满 足最小必需数量,则选择属于其它真图像区域的层候选点。例如,当真图像区域(例如,图25中的真图像区域ni+1)中存在的层候选点的数量不足以在插值中使用时,还选择相邻区域ni+2中的层候选点。
当在当前要插值的方向上不存在具有足够数量的层候选点的真图像区域时(例如,当伪影区域存在于图像区域的边缘时),如下选择层候选点。也就是说,将要插值的方向改变成存在足够的能用于插值的层候选点的方向,并且选择属于该方向附近的真图像区域的层候选点。例如,当伪影区域存在于图像的边缘,如图25所示,并且与x方向相关联的插值所使用的层候选点的数量短缺时,可以选择存在于通过伪影区域的y-z平面上的伪影区域附近的真图像区域中的层候选点。注意,要插值的方向的类型不一定局限于与x轴或y轴平行的方向,并且可以将其改变成存在足够的能用于插值的层候选点的任意方向。例如,考虑到下面的情况,可以选择存在于如图25一样通过圆形扫描所生成的平面上的伪影区域附近的真图像区域中的层候选点:例如,在黄斑部或视神经乳头中,层几何形状趋于具有类似于同心形状的图像特征。将与在该步骤所确定的图像处理方法相关联的信息发送给层确定单元334。
在步骤S2450,层确定单元334利用在该步骤中所选择的类型的插值函数,通过在步骤S2480所选择的层候选点之间进行插值,确定伪影区域中的层位置。然后,将所计算出的层位置的信息输出给存储单元320。
利用上述序列,执行与本实施例相对应的伪影区域图像处理。注意,在本实施例中,在确定图像处理方法之后执行图像校正。然而,执行图像校正的时刻不局限于此。例如,在伪影区域判断单元332的处理之后,图像校正单元333可以根据该区 域中亮度衰减的程度来执行图像校正。在这种情况下,图像处理方法确定单元接收到该图像校正结果,并且响应于亮度使用判断单元的判断结果,进行与插值处理相关联的设置。另外,图像校正方法设置单元1913包括在图像校正单元333中。
根据上述结构,图像处理设备10指定伪影区域,并且判断是否使用伪影区域中的例如边缘的亮度信息。如果使用例如边缘的亮度信息,则在校正该区域中的亮度之后执行层确定处理。当不使用该信息时,根据伪影区域的范围和该区域附近的层几何形状执行插值处理,从而以高精度计算层位置。
第九实施例
在本实施例中,代替第七实施例中仅使用断层图像来计算伪影区域中的层位置,根据被检眼的断层图像生成投影图像,并且将从投影图像所提取的组织或病变部位的位置信息逆投影至断层图像,以缩窄伪影候选区域。
本实施例覆盖以下几点。
(i)将根据投影图像所计算出的血管(出血)区域的位置信息映射至断层图像,并且从所映射的区域的周围区域搜索并指定伪影区域的边缘部分,从而允许以更高精度计算伪影区域的范围。
(ii)由于考虑伪影区域中残留的边缘信息和伪影区域周围的层几何形状的不平坦性,通过对评价函数进行加权来应用层模型,所以可以计算更精确的层位置。
由于与根据本实施例的图像处理设备10连接的设备的结构与第七实施例的相同,所以不给出对其的说明。图26是根据本实施例的图像处理设备10的功能框图。不同于第七实施例的图19所示的图像处理单元1901,本实施例的图像处理单元2601包括投影图像生成单元335和特征提取单元336。由于其它单元与 图19中相同,所以不再重复对其的说明。
下面将参考图27所示的图像处理序列说明本实施例的图像处理的内容。注意,除步骤S2710、S2720和S2730的处理以外,本实施例的处理与第七实施例的图20所示的相同。因此,下面说明步骤S2710、S2720和S2730的处理。
在步骤S2710,投影图像生成单元335通过对断层图像进行投影生成图像。由于实际生成方法与第三实施例的图11的步骤S1110中所述的相同,所以不再重复对其的说明。在步骤S2720,从由投影图像生成单元335所生成的投影图像提取诸如视网膜血管等的组织或病变部位的特征。由于实际生成方法与第三实施例的图11的步骤S1120中所述的相同,所以不再重复对其的说明。
在步骤S2730,当将在步骤S2720所计算出的投影图像上的血管区域(x,y)逆投影至断层图像时,获得图28中以虚线区域2801所表示的区域(以下称为逆投影区域)。通常,在视网膜血管的z轴的正向侧,容易发生亮度衰减。因此,当将所提取的特征(在x-y方向上)的位置逆投影至断层图像时,逆投影得到的虚线区域2801很可能包括伪影区域5。然而,当对错误提取的区域进行逆投影时,在逆投影区域中不发生亮度衰减。即使在对正确提取的视网膜血管进行逆投影时,逆投影得到的区域下方的亮度衰减也轻微,并且在一些情况下,对层提取几乎没有影响。因此,判断在逆投影区域中和该区域的边界附近是否生成了伪影。如果生成了伪影,则计算与该区域中的亮度相关联的统计量。
除作为连续性的计算对象的层候选点的范围以外,伪影区域判断方法与第一实施例的步骤S610~S640中的大体相同。更具体地,不是对于所有层候选点,而是对逆投影区域内部和x-y 方向上该区域附近执行连续性计算处理。
根据上述结构,使用断层图像和投影图像指定伪影区域,并且不仅考虑该区域周围的层几何形状,而且还考虑该区域中的边缘信息,通过对评价函数进行加权来应用层模型,从而以高精度计算层位置。
第十实施例
本实施例在第九实施例中在亮度使用判断单元的判断之后,在使用例如边缘的亮度信息时对伪影区域执行图像校正然后计算层位置,并且在不使用该信息时通过插值处理计算层位置。
本实施例覆盖以下几点:
(i)将根据投影图像所计算出的血管(出血)区域的位置信息映射至断层图像,并且从映射得到的区域的周围区域搜索并指定伪影区域的边缘部分,从而允许以更高精度计算伪影区域的范围。
(ii)在例如残留有衰减的边缘的信息的血管区域中,在转换亮度之后检测层位置,从而更精确地计算层位置。
(iii)在亮度不足并且不能使用例如边缘的亮度信息的区域中,考虑伪影区域的生成位置和周围层几何形状对周围层位置进行插值,从而更精确地计算层位置。
由于与根据本实施例的图像处理设备10连接的设备的结构与第九实施例的相同,所以不再重复对其的说明。图29是根据本实施例的图像处理设备10的功能框图。不同于第九实施例,本实施例的图像处理单元2901还包括对准单元337,并且代替评价函数设置单元1912,图像处理方法确定单元1910包括图像校正方法设置单元1913和插值函数设置单元1914。下面将参考图24和27所示的图像处理序列说明本实施例的图像处理的内容。 注意,除步骤S2020以外的处理与第九实施例的相同。因此,下面将说明步骤S2020的处理,并且不再重复对其它步骤的说明。
在步骤S2020,执行伪影区域中的图像处理。由于该步骤的处理与第八实施例的图24的步骤S2410~S2480中的相同,所以不再重复对其的详细说明。
根据上述结构,图像处理设备10使用断层图像和投影图像指定伪影区域,并且判断是否使用该区域中的例如边缘的亮度信息。当使用该信息时,在校正该区域中的亮度之后执行层确定处理。当不使用该信息时,根据伪影区域的范围和该区域周围的层几何形状执行插值处理,从而以高精度计算层位置。
第十一实施例
本实施例向第九实施例添加了下面的处理:将从被检眼的表面图像和投影图像中至少一个所提取的组织或病变部位的位置信息逆投影至断层图像,以预先缩窄伪影候选区域。
本实施例覆盖以下几点:
(i)当生成从表面图像可容易地提取的诸如白斑等的病变部位时,从表面图像提取白斑区域,将其位置信息映射至断层图像,并且从该区域的周围区域搜索并指定伪影区域的边缘部分,从而以更高精度计算伪影区域的范围。
(ii)由于考虑伪影区域中残留的边缘信息和伪影区域周围的层几何形状的不平坦性,通过对评价函数进行加权来应用层模型,所以可以计算更精确的层位置。
不同于第九实施例,与根据本实施例的图像处理设备10连接的设备的结构还包括表面图像摄像设备50,如图14所示。另外,在根据本实施例的图像处理设备10中,不同于第九实施例,图像处理单元3001包括表面图像获取单元315,如图30所示。
下面将参考图31所示的图像处理序列说明本实施例的图像处理的内容。注意,由于除步骤S3110~S3150以外的处理与第九实施例的相同,所以不再重复对其的说明。另外,步骤S3110~S3140的处理与上述第五实施例的步骤S1610~S1640的相同,并且不再重复对其的说明。
此外,在步骤S3150,执行伪影区域中的图像处理。该步骤的处理与第七实施例的步骤S2110~S2140的大体相同。然而,在本实施例中,在步骤S2130,还可以使用从眼底图像所获得的信息来设置图像处理参数。例如,当眼底图像上的白斑的亮度非常高时,由于亮度在断层图像上的白斑区域的z轴的正向侧很可能衰减,所以与白斑区域的亮度信号值成比例地增大与几何形状相关联的评价函数的权重。作为白斑区域的亮度信号值,可以直接参考眼底图像上的该区域的像素值,或者可以参考作为形态学运算等的处理结果所获得的区域的值(多值数据)。
根据上述结构,在使用表面图像和投影图像所指定的伪影区域中,由于不仅考虑该区域周围的层几何形状,而且还考虑该区域中的亮度信息,通过对评价函数进行加权来应用层模型,因而以高精度来计算层位置。
第十二实施例
本实施例在第十一实施例中在亮度使用判断单元的判断之后,在使用例如边缘的亮度信息时对伪影区域执行图像校正然后确定层位置,并且在不使用该信息时通过插值处理来计算层位置。
尤其当由于白斑而生成伪影时,本实施例覆盖以下几点。
(i)将根据表面图像所计算出的白斑区域的位置信息映射至断层图像,并且从所映射的区域的周围区域搜索并指定伪影区 域的边缘部分,从而允许以更高精度计算伪影区域的范围。
(ii)在例如残留有衰减的边缘的信息的血管区域中,在转换亮度之后检测层位置,从而更精确地计算层位置。
(iii)在亮度缺失并且不能使用例如边缘的亮度信息的区域中,考虑伪影区域的生成位置和周围层几何形状对周围层位置进行插值,从而更精确地计算层位置。
由于与根据本实施例的图像处理设备10连接的设备的结构与第十一实施例的相同,所以不再重复对其的说明。图32是根据本实施例的图像处理设备10的功能框图。在本实施例中,不同于第五实施例,代替评价函数设置单元1912,图像处理方法确定单元1910包括图像校正方法设置单元1913和插值函数设置单元1914。另外,本实施例中的图像处理序列与第十一实施例中的大体相同。然而,如下执行步骤S3150的处理。
然而,在步骤S3150所执行的处理具有与第十实施例中的相同的序列。也就是说,如图24所示,判断是否使用伪影区域中的例如边缘的亮度信息。当使用该信息时,在校正该区域中的亮度之后执行层确定处理。当不使用该信息时,根据该区域的范围和该区域周围的层几何形状设置插值函数的类型和参数,然后执行插值处理。
根据上述结构,在使用表面图像和投影图像所指定的伪影区域中,判断是否使用该区域中的例如边缘的亮度信息。当使用该信息时,在校正该区域中的亮度之后执行层确定处理。当不使用该信息时,根据伪影区域的范围和该区域周围的层几何形状执行插值处理。这样,可以以高精度计算层位置。
其它实施例
上述实施例作为图像处理设备来实现本发明。然而,本发明的实施例不仅局限于图像处理设备,并且可以作为在通过计 算机的CPU执行时实现功能的软件来实现。
图33是示出作为软件实现图像处理设备10的各个单元的功能所使用的计算机的基本结构的框图。CPU 3301使用存储在RAM 3302和ROM 3303中的程序和数据来控制整个计算机。另外,CPU 3301通过控制对与图像处理设备10的各个单元相对应的软件程序的执行来实现各个单元的功能。RAM 3302包括用于临时存储从外部存储装置3304所装载的计算机程序和数据的区域,并且还包括CPU 3301执行各种处理所需的工作区。通过RAM 3302实现存储单元320的功能。
ROM 3303通常存储计算机的BIOS和设置数据等。外部存储装置3304用作诸如硬盘驱动器等的大容量信息存储装置,并且存储由CPU 3301所执行的操作系统和程序。外部存储装置3304存储本实施例的说明中所给出的信息,并且根据需要将这些信息装载进RAM 3302。由例如液晶显示器构成监视器3305。例如,监视器3305可以显示从显示单元340输出的内容。
键盘3306和鼠标3307是输入装置。操作者可以使用这些输入装置向图像处理设备10输入各种指示。接口3308用于在图像处理设备10和外部设备之间交换各种数据,并且可以由例如IEEE 1394、USB或 端口等构成。将经由接口3308所获取的数据取入RAM 3302。经由接口3308实现断层图像获取单元310和结果输出单元350的功能。上述组件经由总线3309相互连接。
还可以利用读出并执行记录在存储器装置上的程序以进行上述实施例的功能的系统或设备的计算机(或者CPU或MPU等装置)和通过下面的方法实现本发明的各方面,其中,利用系统或设备的计算机通过例如读出并执行记录在存储器装置上的程序以进行上述实施例的功能来进行上述方法的各步骤。为此, 例如,通过网络或者通过用作存储器装置的各种类型的记录介质(例如,计算机可读介质)将该程序提供给计算机。
尽管参考典型实施例说明了本发明,但是应该理解,本发明不局限于所公开的典型实施例。所附权利要求书的范围符合最宽的解释,以包含所有这类修改、等同结构和功能。
本申请要求2009年6月2日提交的日本2009-133455号专利申请的优先权,其全部内容通过引用包含于此。
Claims (21)
1.一种图像处理设备,其对通过断层图像摄像设备拍摄被检眼的图像所获得的断层图像进行处理,所述图像处理设备包括:
层候选检测部件,用于从所述断层图像检测所述被检眼的视网膜的层候选;
伪影区域判断部件,用于基于使用所述层候选所获得的图像特征,判断所述断层图像中的伪影区域;以及
图像校正部件,用于基于所述伪影区域判断部件的判断结果和所述伪影区域中的图像特征,校正所述伪影区域中的亮度。
2.根据权利要求1所述的图像处理设备,其特征在于,还包括层确定部件,所述层确定部件用于基于通过所述图像校正部件校正后的所述伪影区域中的亮度,确定所述伪影区域中的视网膜的层的位置。
3.根据权利要求2所述的图像处理设备,其特征在于,还包括亮度使用判断部件,所述亮度使用判断部件用于基于所述伪影区域中的亮度的大小,判断在所述层的位置的确定中是否使用所述伪影区域中的亮度,
其中,当所述亮度使用判断部件判断为使用所述伪影区域中的亮度时,所述层确定部件使用所述伪影区域中的亮度执行所述层的位置的确定处理,以及
当所述亮度使用判断部件判断为不使用所述伪影区域中的亮度时,所述层确定部件通过使用所述伪影区域附近的层候选进行插值处理,确定所述层的位置。
4.一种图像处理设备,其对通过断层图像摄像设备拍摄被检眼的图像所获得的断层图像进行处理,所述图像处理设备包括:
层候选检测部件,用于从所述断层图像检测所述被检眼的视网膜的层候选;
伪影区域判断部件,用于基于使用所述层候选所获得的图像特征,判断所述断层图像中的伪影区域;以及
层确定部件,用于基于所述伪影区域判断部件的判断结果,确定所述伪影区域中的所述视网膜的层的位置,
其中,所述层确定部件使用用于指定所述伪影区域中所包括的层几何形状并且由多个控制点所定义的几何模型,并且基于与所述几何模型的几何形状相关联的评价函数和与所述控制点附近的亮度相关联的评价函数来确定所述层的位置。
5.根据权利要求4所述的图像处理设备,其特征在于,所述层确定部件以随着所述伪影区域中像素的亮度的降低而增大与所述控制点附近的亮度相关联的评价函数的权重的方式来确定与所述控制点附近的亮度相关联的评价函数的权重,或者根据所述伪影区域中像素的亮度和所述断层图像中的除所述伪影区域以外的像素的亮度之间的比来确定与所述控制点附近的亮度相关联的评价函数的权重,以及
所述层确定部件根据由所述伪影区域附近的所述层候选所指定的层的几何形状的不平坦程度,确定与所述几何模型的几何形状相关联的评价函数的权重,以使得随着不平坦程度的增大而降低该权重。
6.根据权利要求1~5中任一项所述的图像处理设备,其特征在于,使用所述层候选所获得的图像特征包括所述层候选之间的连续性,以及
所述伪影区域判断部件判断为:所述断层图像中的所述伪影区域具有连续性低且被判断为不连续的层候选作为所述伪影区域的端点。
7.根据权利要求6所述的图像处理设备,其特征在于,所述伪影区域判断部件使用如下的层候选作为所述伪影区域的端点:所述层候选是被判断为不连续的两个层候选中、比各个层候选深的区域中的亮度较低的层候选。
8.根据权利要求1~5中任一项所述的图像处理设备,其特征在于,使用所述层候选所获得的图像特征包括由两个层候选所指定的区域中像素的亮度在与所述断层图像的A扫描方向垂直的方向上的差,以及
所述伪影区域判断部件基于A扫描方向上的多个亮度剖面图之间的差异程度,判断所述断层图像中的所述伪影区域,其中,所述A扫描方向为与所述视网膜的深度方向相对应的断层图像的纵向。
9.根据权利要求1~5中任一项所述的图像处理设备,其特征在于,还包括:
投影图像生成部件,用于通过在与所述视网膜的深度方向相对应的所述断层图像的A扫描方向上将所述断层图像的各个像素相加,生成投影图像,其中,所述A扫描方向为与所述视网膜的深度方向相对应的断层图像的纵向;以及
特征提取部件,用于从所述投影图像提取包括所述被检眼中的生物组织和病变部位中的至少一个的特征区域,
其中,所述伪影区域判断部件在所述特征区域中和所述特征区域附近判断所述伪影区域。
10.根据权利要求9所述的图像处理设备,其特征在于,所述特征提取部件还从通过拍摄所述被检眼的图像所获得的表面图像提取所述特征区域,以及
所述伪影区域判断部件通过还使用从所述表面图像所提取的特征区域判断所述伪影区域。
11.一种图像处理设备的控制方法,所述图像处理设备对通过断层图像摄像设备拍摄被检眼的图像所获得的断层图像进行处理,所述控制方法包括以下步骤:
层候选检测步骤,用于控制层候选检测部件,以从所述断层图像检测所述被检眼的视网膜的层候选;
伪影区域判断步骤,用于控制伪影区域判断部件,以基于使用所述层候选所获得的图像特征,判断所述断层图像中的伪影区域;以及
图像校正步骤,用于控制图像校正部件,以基于所述伪影区域判断步骤中的判断结果和所述伪影区域中的图像特征,校正所述伪影区域中的亮度。
12.一种图像处理设备的控制方法,所述图像处理设备对通过断层图像摄像设备拍摄被检眼的图像所获得的断层图像进行处理,所述控制方法包括以下步骤:
层候选检测步骤,用于控制层候选检测部件,以从所述断层图像检测所述被检眼的视网膜的层候选;
伪影区域判断步骤,用于控制伪影区域判断部件,以基于使用所述层候选所获得的图像特征,判断所述断层图像中的伪影区域;以及
层确定步骤,用于控制层确定部件,以基于所述伪影区域判断步骤中的判断结果,确定所述伪影区域中的所述视网膜的层的位置,
其中,在所述层确定步骤中,使用用于指定所述伪影区域中所包括的层几何形状并且由多个控制点所定义的几何模型,并且基于与所述几何模型的几何形状相关联的评价函数和与所述控制点附近的亮度相关联的评价函数,确定所述层的位置。
13.一种用于对通过断层图像摄像设备拍摄被检眼的图像所获得的断层图像进行处理的方法,所述方法包括以下步骤:
识别步骤,用于基于在所述断层图像中所检测到的预定层边界的候选点的连续性来识别所述断层图像中的伪影区域;
第一检测步骤,用于在除所述伪影区域以外的区域中,检测构成所述被检眼的视网膜的层的边界;
第二检测步骤,用于在所述断层图像的所述伪影区域中,基于所述伪影区域的亮度检测构成所述被检眼的视网膜的层的边界;以及
显示步骤,用于通过将表示所检测到的层的边界的线重叠在所述断层图像上来显示所述线。
14.根据权利要求13所述的方法,其特征在于,在所述第二检测步骤中,基于所述伪影区域的亮度和利用所述第一检测步骤所检测到的层的边界的位置,检测所述伪影区域中的所述层的边界。
15.根据权利要求13所述的方法,其特征在于,在所述识别步骤中,使用所述断层图像和/或所述被检眼的眼底图像提取血管区域和/或白斑区域,并且将所述断层图像中与所提取的区域相对应的区域识别为所述伪影区域。
16.根据权利要求13所述的方法,其特征在于,在所述识别步骤中,判断在所述断层图像中与血管区域或白斑区域相对应的区域中是否存在具有衰减的亮度的所述伪影区域。
17.根据权利要求13所述的方法,其特征在于,在所述识别步骤中,将在从所述被检眼的前方看时位于所述视网膜的血管区域的后方的区域识别为所述伪影区域。
18.一种图像处理设备,其进行根据权利要求13~17中任一项所述的方法,所述图像处理设备包括:
识别部件,用于进行所述识别步骤;
检测部件,用于进行所述第一检测步骤和所述第二检测步骤;以及
显示部件,用于进行所述显示步骤。
19.一种断层图像摄像系统,包括:
根据权利要求18所述的图像处理设备;以及
断层图像摄像设备。
20.一种图像处理设备,其对通过断层图像摄像设备拍摄被检眼的图像所获得的断层图像进行处理,所述图像处理设备包括:
识别部件,用于基于在所述断层图像中所检测到的预定层边界的候选点的连续性来识别所述断层图像中的伪影区域;
检测部件,用于在除所述伪影区域以外的区域中,检测构成所述被检眼的视网膜的层的边界,并且用于在所述断层图像的所述伪影区域中,基于所述伪影区域的亮度检测构成所述被检眼的视网膜的层的边界;以及
显示部件,用于通过将表示所检测到的层的边界的线重叠在所述断层图像上来显示所述线。
21.一种用于对通过断层图像摄像设备拍摄被检眼的图像所获得的断层图像进行处理的方法,所述方法包括以下步骤:
第一检测步骤,用于基于构成所述被检眼的视网膜的层的各自的亮度,在所述断层图像中检测所述层的边界;
识别步骤,用于识别在所述第一检测步骤中所检测到的边界不连续的区域;
第二检测步骤,用于基于所识别出的区域的亮度检测所识别出的区域中的层的边界;以及
显示步骤,用于通过将表示所检测到的层的边界的线重叠在所述断层图像上来显示所述线。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009133455A JP4909378B2 (ja) | 2009-06-02 | 2009-06-02 | 画像処理装置及びその制御方法、コンピュータプログラム |
JP2009-133455 | 2009-06-02 | ||
PCT/JP2010/059302 WO2010140601A1 (en) | 2009-06-02 | 2010-05-26 | Image processing apparatus, control method thereof, and computer program |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102458221A CN102458221A (zh) | 2012-05-16 |
CN102458221B true CN102458221B (zh) | 2014-08-27 |
Family
ID=43297734
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201080024458.0A Expired - Fee Related CN102458221B (zh) | 2009-06-02 | 2010-05-26 | 图像处理设备及其控制方法 |
Country Status (6)
Country | Link |
---|---|
US (1) | US8861817B2 (zh) |
EP (1) | EP2437648B1 (zh) |
JP (1) | JP4909378B2 (zh) |
KR (1) | KR101318685B1 (zh) |
CN (1) | CN102458221B (zh) |
WO (1) | WO2010140601A1 (zh) |
Families Citing this family (68)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4850927B2 (ja) * | 2009-06-02 | 2012-01-11 | キヤノン株式会社 | 画像処理装置、画像処理方法及びコンピュータプログラム |
JP5436076B2 (ja) | 2009-07-14 | 2014-03-05 | キヤノン株式会社 | 画像処理装置、画像処理方法およびプログラム |
JP5645432B2 (ja) * | 2010-03-19 | 2014-12-24 | キヤノン株式会社 | 画像処理装置、画像処理システム、画像処理方法、及び画像処理をコンピュータに実行させるためのプログラム |
JP5657941B2 (ja) * | 2010-07-30 | 2015-01-21 | 株式会社トプコン | 光断層画像化装置及びその作動方法 |
EP2609853A4 (en) * | 2010-08-27 | 2016-03-09 | Sony Corp | APPARATUS AND METHOD FOR IMAGE PROCESSING |
JP5702991B2 (ja) | 2010-11-19 | 2015-04-15 | キヤノン株式会社 | 画像処理装置及び画像処理方法 |
JP5701024B2 (ja) * | 2010-11-26 | 2015-04-15 | キヤノン株式会社 | 画像処理装置及び方法 |
JP5242667B2 (ja) | 2010-12-22 | 2013-07-24 | 株式会社東芝 | マップ変換方法、マップ変換装置及びマップ変換プログラム |
JP5904711B2 (ja) | 2011-02-01 | 2016-04-20 | キヤノン株式会社 | 画像処理装置、画像処理方法およびプログラム |
CN102136135B (zh) * | 2011-03-16 | 2012-12-19 | 清华大学 | 眼前节光学相干层析图像的角膜及眼前房内轮廓提取方法 |
JP5836634B2 (ja) | 2011-05-10 | 2015-12-24 | キヤノン株式会社 | 画像処理装置及び方法 |
JP6025311B2 (ja) | 2011-08-01 | 2016-11-16 | キヤノン株式会社 | 眼科診断支援装置および方法 |
JP5955163B2 (ja) | 2011-09-06 | 2016-07-20 | キヤノン株式会社 | 画像処理装置および画像処理方法 |
JP5912358B2 (ja) * | 2011-09-14 | 2016-04-27 | 株式会社トプコン | 眼底観察装置 |
US8864309B2 (en) | 2011-12-05 | 2014-10-21 | Bioptigen, Inc. | Optical imaging systems having input beam shape control and path length control |
JP6101048B2 (ja) | 2012-02-20 | 2017-03-22 | キヤノン株式会社 | 画像処理装置及び画像処理方法 |
JP5947561B2 (ja) | 2012-02-20 | 2016-07-06 | キヤノン株式会社 | 画像形成装置及び画像形成方法 |
JP6114495B2 (ja) | 2012-02-20 | 2017-04-12 | キヤノン株式会社 | 画像表示装置及び画像表示方法、撮影システム |
JP5979904B2 (ja) * | 2012-02-20 | 2016-08-31 | キヤノン株式会社 | 画像処理装置、眼科撮影システム、及び画像処理方法 |
US8777412B2 (en) | 2012-04-05 | 2014-07-15 | Bioptigen, Inc. | Surgical microscopes using optical coherence tomography and related methods |
JP2014045869A (ja) * | 2012-08-30 | 2014-03-17 | Canon Inc | 撮影装置、画像処理装置、及び画像処理方法 |
JP6099903B2 (ja) * | 2012-08-30 | 2017-03-22 | キヤノン株式会社 | 画像処理装置および画像処理方法 |
JP6075844B2 (ja) * | 2012-10-18 | 2017-02-08 | キヤノン株式会社 | 眼科装置および眼科方法並びに記憶媒体 |
JP2014110883A (ja) * | 2012-10-30 | 2014-06-19 | Canon Inc | 画像処理装置及び画像処理方法 |
GB201302887D0 (en) * | 2013-02-19 | 2013-04-03 | Optos Plc | Improvements in or relating to image processing |
JP6092659B2 (ja) | 2013-02-28 | 2017-03-08 | キヤノン株式会社 | 画像処理装置及び画像処理方法 |
JP6198410B2 (ja) | 2013-02-28 | 2017-09-20 | キヤノン株式会社 | 画像処理装置及び画像処理方法 |
JP6200168B2 (ja) | 2013-02-28 | 2017-09-20 | キヤノン株式会社 | 画像処理装置及び画像処理方法 |
CN104036472A (zh) * | 2013-03-06 | 2014-09-10 | 北京三星通信技术研究有限公司 | 用于3d图像质量增强的方法和设备 |
US9401009B2 (en) | 2013-03-06 | 2016-07-26 | Samsung Electronics Co., Ltd. | Method and apparatus for enhancing quality of 3D image |
US10497124B2 (en) * | 2013-03-15 | 2019-12-03 | Kabushiki Kaisha Topcon | Optic disc image segmentation method and apparatus |
EP3003123B1 (en) | 2013-06-04 | 2020-08-19 | Bioptigen, Inc. | Optical coherence tomography imaging system and optical laser scanning system comprising a beam shaping optical system with a +-+ lens triplet, where second and third lens are movable, and method |
US10456030B2 (en) | 2013-07-29 | 2019-10-29 | Bioptigen, Inc. | Procedural optical coherence tomography (OCT) for surgery and related methods |
WO2015031589A1 (en) | 2013-08-28 | 2015-03-05 | Bioptigen, Inc. | Heads up displays for optical coherence tomography integrated surgical microscopes |
JP6512738B2 (ja) * | 2014-02-28 | 2019-05-15 | キヤノン株式会社 | 画像処理装置、画像処理方法及びプログラム |
DE102014113901A1 (de) | 2014-09-25 | 2016-03-31 | Carl Zeiss Meditec Ag | Verfahren zur Korrektur eines OCT-Bildes und Kombinationsmikroskop |
US9675247B2 (en) * | 2014-12-05 | 2017-06-13 | Ricoh Co., Ltd. | Alpha-matting based retinal vessel extraction |
JP6470604B2 (ja) | 2015-03-25 | 2019-02-13 | キヤノン株式会社 | 画像処理装置及び画像処理方法及びプログラム |
US9639917B2 (en) * | 2015-05-19 | 2017-05-02 | Novartis Ag | OCT image modification |
US9508127B1 (en) * | 2015-06-30 | 2016-11-29 | Konica Minolta Laboratory U.S.A., Inc. | Processing for creating a transmission image without artificial noises |
JP6736270B2 (ja) | 2015-07-13 | 2020-08-05 | キヤノン株式会社 | 画像処理装置及び画像処理装置の作動方法 |
WO2017030057A1 (ja) * | 2015-08-19 | 2017-02-23 | 興和株式会社 | 画像処理装置、画像処理方法、及び画像処理プログラム |
PL3362987T3 (pl) | 2015-10-14 | 2022-01-10 | Shanghai United Imaging Healthcare Co., Ltd. | Układ i sposób korekcji obrazu |
CN106651986B (zh) * | 2016-01-21 | 2021-05-18 | 上海联影医疗科技股份有限公司 | 计算机断层成像伪影校正方法 |
CN105809188B (zh) * | 2016-02-26 | 2017-08-08 | 山东大学 | 一种基于ambp改进算法的真菌性角膜炎图像识别方法 |
JP6156855B2 (ja) * | 2016-03-02 | 2017-07-05 | 富士フイルム株式会社 | 分類装置、分類プログラムおよび分類装置の動作方法 |
JP6971559B2 (ja) * | 2016-10-31 | 2021-11-24 | キヤノン株式会社 | 情報処理装置、情報処理方法およびプログラム |
US10896490B2 (en) * | 2016-12-23 | 2021-01-19 | Oregon Health & Science University | Systems and methods for reflectance-based projection-resolved optical coherence tomography angiography |
CN108959215B (zh) * | 2017-05-24 | 2022-06-10 | 阿里健康信息技术有限公司 | 一种数据处理方法和设备、显示方法和设备 |
KR20200005406A (ko) * | 2018-07-06 | 2020-01-15 | 연세대학교 산학협력단 | 진단 보조 시스템 |
CN111712186A (zh) | 2017-12-20 | 2020-09-25 | 医鲸股份有限公司 | 用于辅助心血管疾病的诊断的方法和装置 |
JP7108227B2 (ja) * | 2018-03-02 | 2022-07-28 | 株式会社ニデック | 自覚式検眼装置及び自覚式検眼プログラム |
JP7130989B2 (ja) * | 2018-03-05 | 2022-09-06 | 株式会社ニデック | 眼科画像処理装置、および眼科画像処理プログラム |
US10973406B2 (en) | 2018-03-06 | 2021-04-13 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, and non-transitory computer readable medium |
IT201800004940A1 (it) * | 2018-04-27 | 2019-10-27 | Metodo per determinare una curva di trabecolato basato sulla applicazione di soglia a una immagine a canale singolo derivata da una immagine a colori | |
IT201800004934A1 (it) * | 2018-04-27 | 2019-10-27 | Metodo per determinare una curva di trabecolato basato su minimi locali negativi multipli in una immagine a toni-di-grigio | |
WO2019207481A1 (en) * | 2018-04-27 | 2019-10-31 | Nidek Technologies Srl | Methods of determining a trabecular meshwork curve |
IT201800004938A1 (it) * | 2018-04-27 | 2019-10-27 | Metodo per determinare una curva di trabecolato basato sulla applicazione di soglie ai canali di una immagine a colori | |
JP7166182B2 (ja) * | 2019-01-28 | 2022-11-07 | 株式会社トプコン | 眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラム |
US11417144B2 (en) | 2019-02-22 | 2022-08-16 | Nec Corporation | Processing apparatus, fingerprint image extraction processing apparatus, system, processing method, and computer readable medium |
US10832074B2 (en) | 2019-03-08 | 2020-11-10 | International Business Machines Corporation | Uncertainty region based image enhancement |
JP6888643B2 (ja) * | 2019-03-22 | 2021-06-16 | 株式会社ニデック | Oct解析処理装置、及びoctデータ処理プログラム |
JP2020171664A (ja) * | 2019-04-09 | 2020-10-22 | キヤノン株式会社 | 画像解析方法及び画像解析装置 |
WO2020209012A1 (ja) * | 2019-04-09 | 2020-10-15 | キヤノン株式会社 | 画像解析方法及び画像解析装置 |
CN110310235B (zh) * | 2019-05-21 | 2021-07-27 | 北京至真互联网技术有限公司 | 眼底图像处理方法、装置及设备和存储介质 |
EP4055553A1 (en) * | 2020-01-08 | 2022-09-14 | Haag-Streit Ag | Optical coherence tomography system for ophthalmology |
US11875480B1 (en) * | 2020-05-15 | 2024-01-16 | Verily Life Sciences Llc | Distinguishing artifacts from pathological features in digital images |
CN111861917A (zh) * | 2020-07-10 | 2020-10-30 | 温州医科大学 | 一种基于信号逆向补偿的脉络膜oct图像增强方法及装置 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4527251B2 (ja) * | 2000-07-14 | 2010-08-18 | 興和株式会社 | 眼科測定装置 |
WO2004098396A2 (en) | 2003-05-01 | 2004-11-18 | The Cleveland Clinic Foundation | Method and apparatus for measuring a retinal sublayer characteristic |
JP4461258B2 (ja) * | 2005-10-31 | 2010-05-12 | 国立大学法人 筑波大学 | 光断層画像化法における補正方法 |
EP1976424B1 (en) * | 2006-01-19 | 2011-09-21 | Optovue, Inc. | A method of eye examination by optical coherence tomography |
JP4855150B2 (ja) | 2006-06-09 | 2012-01-18 | 株式会社トプコン | 眼底観察装置、眼科画像処理装置及び眼科画像処理プログラム |
JP5095167B2 (ja) | 2006-09-19 | 2012-12-12 | 株式会社トプコン | 眼底観察装置、眼底画像表示装置及び眼底観察プログラム |
JP5058627B2 (ja) * | 2007-02-26 | 2012-10-24 | 株式会社トプコン | 眼底観察装置 |
JP5101969B2 (ja) * | 2007-09-27 | 2012-12-19 | 株式会社トプコン | 眼科装置 |
JP5101975B2 (ja) * | 2007-10-04 | 2012-12-19 | 株式会社トプコン | 眼底観察装置及び眼底画像処理装置 |
JP2010142428A (ja) | 2008-12-18 | 2010-07-01 | Canon Inc | 撮影装置及び撮影方法、プログラム、記録媒体 |
JP4909377B2 (ja) * | 2009-06-02 | 2012-04-04 | キヤノン株式会社 | 画像処理装置及びその制御方法、コンピュータプログラム |
JP5436076B2 (ja) | 2009-07-14 | 2014-03-05 | キヤノン株式会社 | 画像処理装置、画像処理方法およびプログラム |
JP5025715B2 (ja) | 2009-12-08 | 2012-09-12 | キヤノン株式会社 | 断層画像撮影装置、画像処理装置、画像処理システム、画像処理装置の制御方法及びプログラム |
JP5582772B2 (ja) | 2009-12-08 | 2014-09-03 | キヤノン株式会社 | 画像処理装置及び画像処理方法 |
JP5702991B2 (ja) | 2010-11-19 | 2015-04-15 | キヤノン株式会社 | 画像処理装置及び画像処理方法 |
JP5904711B2 (ja) | 2011-02-01 | 2016-04-20 | キヤノン株式会社 | 画像処理装置、画像処理方法およびプログラム |
-
2009
- 2009-06-02 JP JP2009133455A patent/JP4909378B2/ja active Active
-
2010
- 2010-05-26 CN CN201080024458.0A patent/CN102458221B/zh not_active Expired - Fee Related
- 2010-05-26 US US13/318,848 patent/US8861817B2/en not_active Expired - Fee Related
- 2010-05-26 EP EP10783383.2A patent/EP2437648B1/en active Active
- 2010-05-26 KR KR1020117030891A patent/KR101318685B1/ko active IP Right Grant
- 2010-05-26 WO PCT/JP2010/059302 patent/WO2010140601A1/en active Application Filing
Non-Patent Citations (3)
Title |
---|
JP特开2002-28134A 2002.01.29 |
JP特开2007-127425A 2007.05.24 |
JP特开2009-78088A 2009.04.16 |
Also Published As
Publication number | Publication date |
---|---|
KR20120024842A (ko) | 2012-03-14 |
US8861817B2 (en) | 2014-10-14 |
EP2437648A4 (en) | 2014-09-17 |
KR101318685B1 (ko) | 2013-10-16 |
CN102458221A (zh) | 2012-05-16 |
EP2437648B1 (en) | 2019-09-11 |
EP2437648A1 (en) | 2012-04-11 |
US20120063660A1 (en) | 2012-03-15 |
WO2010140601A1 (en) | 2010-12-09 |
JP2010279440A (ja) | 2010-12-16 |
JP4909378B2 (ja) | 2012-04-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102458221B (zh) | 图像处理设备及其控制方法 | |
CN102458225B (zh) | 图像处理设备及其控制方法 | |
CN102469935B (zh) | 图像处理设备和图像处理方法 | |
JP5955163B2 (ja) | 画像処理装置および画像処理方法 | |
CN103717122B (zh) | 眼科诊断支持设备和眼科诊断支持方法 | |
US8699774B2 (en) | Image processing apparatus, control method thereof, and program | |
CN102525405B (zh) | 图像处理设备和图像处理方法 | |
Lupascu et al. | Automated detection of optic disc location in retinal images | |
Gustavsson et al. | A dynamic programming procedure for automated ultrasonic measurement of the carotid artery | |
JP5174232B2 (ja) | 画像処理方法、画像処理装置、断層像撮像システム及び、コンピュータプログラム | |
EP2772185A1 (en) | Image processing apparatus and image processing method | |
US20140307933A1 (en) | Methods to reduce variance in oct analysis of the macula | |
CN103565404B (zh) | 识别关注对象的特征及处理断层图像的方法和设备 | |
JP2013153844A (ja) | 画像処理装置、画像処理方法及びプログラム。 | |
Li et al. | Retinal vessel detection and measurement for computer-aided medical diagnosis | |
CN112164043A (zh) | 一种用于多眼底图像拼接的方法及系统 | |
JP6864450B2 (ja) | 眼科撮影装置 | |
JP6892234B2 (ja) | 眼科撮影装置 | |
CN110211021A (zh) | 图像处理装置、图像处理方法和存储介质 | |
US6134341A (en) | Method for analyzing a stereoscopic image of a fundus, and an apparatus for executing that method | |
Niu et al. | Automated detection of foveal center in SD‐OCT images using the saliency of retinal thickness maps | |
CA2731385C (en) | Image analysis system and method | |
US20180240240A1 (en) | Image processing apparatus, image processing method, and image processing program | |
Corona et al. | Digital stereo-optic disc image analyzer for monitoring progression of glaucoma | |
Wang | Estimation of papilledema severity using spectral-domain optical coherence tomography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140827 Termination date: 20190526 |
|
CF01 | Termination of patent right due to non-payment of annual fee |