CN111246792B - 解析装置、解析方法、以及程序 - Google Patents
解析装置、解析方法、以及程序 Download PDFInfo
- Publication number
- CN111246792B CN111246792B CN201880067903.8A CN201880067903A CN111246792B CN 111246792 B CN111246792 B CN 111246792B CN 201880067903 A CN201880067903 A CN 201880067903A CN 111246792 B CN111246792 B CN 111246792B
- Authority
- CN
- China
- Prior art keywords
- region
- unit
- space
- detection region
- detection
- 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
Images
Classifications
-
- 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
-
- 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
-
- 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/1025—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for confocal scanning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0002—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
- A61B5/0004—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network characterised by the type of physiological signal transmitted
- A61B5/0013—Medical image data
-
- 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
-
- 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/0071—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by measuring fluorescence emission
-
- 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/10076—4D tomography; Time-sequential 3D tomography
-
- 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
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Ophthalmology & Optometry (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Quality & Reliability (AREA)
- Geometry (AREA)
- Signal Processing (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Eye Examination Apparatus (AREA)
- Image Analysis (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
解析装置具备:获取部,获取将物理量与划分三维空间、三维时空间、或者四维时空间而成的每个单位区域建立对应关系的数据;设定部,在上述空间或者时空间中,设定作为三维以上的时空间的检测区域,在与上述检测区域的维相同的空间中,设置包围上述检测区域的间隙,并在包围上述间隙的位置设定对照区域;以及判定部,基于相互比较由上述设定部设定的上述检测区域以及上述对照区域各自所包含的一个以上的上述单位区域的物理量,判定上述检测区域所包含的一个以上的上述单位区域是否为规定区域。
Description
技术领域
本发明涉及解析装置、解析方法、以及程序。
背景技术
利用二维或者三维地拍摄物体的装置。该装置例如包括激光扫描显微镜、共聚焦激光显微镜或者荧光显微镜等光学显微镜、TEM(Transmission Electron Microscope:透射电子显微镜)或SEM(Scanning Electron Microscope:扫描电子显微镜)等电子显微镜、利用OCT(Optical Coherence Tomography:光学相干断层扫描)、CT(ComputedTomography:计算机断层扫描)、PET(Positron Emission Tomography:正电子发射断层扫描)、MRI(Magnetic Resonance Imaging:磁共振成像)、NMR(Nuclear MagneticResonance:核磁共振)或者超声波等拍摄物体的断层图像的断层图像拍摄装置等各种装置。这样的装置例如在医疗现场被利用于拍摄患者的患部等。
专利文献1:日本特开平9-84763号公报
非专利文献1:Yabusaki,Katsumi,Joshua D.Hutcheson,Payal Vyas,SergioBertazzo,Simon C.Body,Masanori Aikawa,and Elena Aikawa.2016.“Quantificationof Calcified Particles in Human Valve Tissue Reveals Asymmetry of CalcificAortic Valve Disease Development.”Frontiers in Cardiovascular Medicine 3(1):44.doi:10.3389/fcvm.2016.00044.http://dx.doi.org/10.3389/fcvm.2016.00044.
然而,在现有技术中,在想要检测的对象物较微小的情况下,有不能精度良好地检测该对象物的情况。
发明内容
本发明是考虑这样的情况而完成的,其目的之一在于提供能够提高对象物的检测精度的解析装置、解析方法、以及程序。
解决上述问题的本发明是对非专利文献1所示的方法的思想进行改进以使其也能够应用于三维时空间(三维空间或者二维空间和时间轴)、或者对三维空间加上时间轴的四维时空间的发明,其一方式是解析装置,具备:获取部,获取将物理量与每个单位区域建立对应关系的数据,其中,上述单位区域是划分三维空间、三维时空间、或者四维时空间而成的;设定部,在上述空间或者时空间中,设定作为三维以上的时空间的检测区域,在与上述检测区域的维相同的空间中,设置包围上述检测区域的间隙,并在包围上述间隙的位置设定对照区域;以及判定部,基于相互比较由上述设定部设定的上述检测区域以及上述对照区域各自所包含的一个以上的上述单位区域的物理量,判定上述检测区域所包含的一个以上的上述单位区域是否为规定区域。
根据本发明,能够提高对象物的检测精度。
附图说明
图1是表示包含第一实施方式的解析装置200的检查系统1的构成的一个例子的图。
图2是表示测量装置100的构成的一个例子的图。
图3是表示解析装置200的构成的一个例子的图。
图4是示意地表示体数据212的图。
图5是表示扫描位置信息214的一个例子的图。
图6是表示控制部220所进行的一系列的处理的一个例子的流程图。
图7是表示控制部220所进行的一系列的处理的一个例子的流程图。
图8A是表示检测区域Ra的设定位置的一个例子的图。
图8B是表示检测区域Ra的设定位置的一个例子的图。
图8C是表示检测区域Ra的设定位置的一个例子的图。
图8D是表示检测区域Ra的设定位置的一个例子的图。
图8E是表示检测区域Ra的设定位置的一个例子的图。
图8F是表示检测区域Ra的设定位置的一个例子的图。
图9A是表示检测区域Ra的设定位置的其它的例子的图。
图9B是表示检测区域Ra的设定位置的其它的例子的图。
图9C是表示检测区域Ra的设定位置的其它的例子的图。
图9D是表示检测区域Ra的设定位置的其它的例子的图。
图9E是表示检测区域Ra的设定位置的其它的例子的图。
图10是表示检测区域Ra的设定位置的其它的例子的图。
图11是表示设置与多个体素相对应的间隙而设定的对照区域Rb的一个例子的图。
图12是表示设置与多个体素相对应的间隙而设定的对照区域Rb的一个例子的图。
图13是表示配置为圆环状的对照区域Rb的一个例子的图。
图14是表示配置为圆环状的对照区域Rb的一个例子的图。
图15是表示被单个对照区域Rb包围的检测区域Ra的一个例子的图。
图16是表示对照区域Rb的设定方法的其它的例子的图。
图17是用于说明在以包围检测区域Ra的方式设定对照区域Rb时,在设置间隙时和不设置间隙时判定结果不同的图。
图18是用于说明在以包围检测区域Ra的方式设定对照区域Rb时,在设置间隙时和不设置间隙时判定结果不同的图。
图19是表示体素聚类CL的生成方法的一个例子的图。
图20是表示体素聚类CL的生成方法的一个例子的图。
图21是表示体素聚类CL的生成方法的一个例子的图。
图22是表示显示了基于评价判定部232的评价结果的图像的显示部204的一个例子的图。
图23是示意地表示眼球E的剖面的图。
图24A是表示观测时刻不同的两个体数据212的一剖面的图。
图24B是表示观测时刻不同的两个体数据212的一剖面的图。
图25是表示第三实施方式的检查系统2的构成的一个例子的图。
图26是表示由OCT拍摄装置300生成的体数据212的一个例子的图。
图27是表示实施方式的解析装置200的硬件构成的一个例子的图。
具体实施方式
以下,参照附图,对本发明的解析装置、解析方法、以及程序的实施方式进行说明。
[概要]
实施方式的解析装置获取将物理量与划分三维空间、三维时空间、或者四维时空间而成的每个单位区域建立对应关系的体数据,在该体数据内,设定作为三维以上的时空间的检测区域,并在进一步包围了包围该检测区域的间隙的位置,设定一个或者多个对照区域。
体数据例如是按某单位区域(以下,称为体素)划分将宽度、高度、进深作为维的三维空间,并且将物理量与该划分后的每个体素建立对应关系的数据。物理量是标量值或者向量值,例如,是激光的反射光强度、电磁波(电波或者电子等)的反射波强度、声波的反射波强度、图像的像素值(亮度值等)、压力、流速、温度、吸收光谱等。另外,体数据例如也可以是按体素划分将宽度、高度、进深三个中的任意两个、和时间作为维的三维时空间得到的数据,也可以是按体素(或者动态体素)划分将宽度、高度、进深、时间作为维的四维时空间而得到的数据,也可以是按体素划分将其它的要素作为维的多维空间而得到的数据。
然后,解析装置基于相互比较检测区域以及对照区域各自所包含的一个以上的体素的物理量,判定检测区域所包含的一个以上的体素是否是规定区域,并基于该判定结果,评价在生成体数据时观测到的观测对象。
(第一实施方式)
[系统构成]
图1是表示包含第一实施方式的解析装置200的检查系统1的构成的一个例子的图。检查系统1例如是检查人或者其它的动物的眼球E的系统。例如,检查系统1具备测量装置100和解析装置200。测量装置100例如向作为被检体的眼球E的前房内照射激光,并接收该照射的激光的散射光(反射光),并将表示该受光结果的信息提供给解析装置200。前房是指虹膜与角膜的最内层的内皮细胞之间的区域,充满被称为眼房水的液体。解析装置200基于由测量装置100提供的信息,求出在眼球E的前房内存在的生物体细胞(cell)的数目,对眼球E的状态进行评价。以下,将激光的照射方向(眼球E的深度方向)设为z方向(轴),将与该z方向正交的某一个方向设为x方向(轴),并将与z方向以及x方向双方正交的方向设为y方向(轴)来进行说明。即,z方向表示三维空间的进深方向,x方向表示三维空间的宽度方向以及高度方向中一方,y方向表示三维空间的宽度方向以及高度方向中另一方。此外,也可以x方向、y方向、以及z方向中任意一个是时间的维(t方向)。
[测量装置的构成]
图2是表示测量装置100的构成的一个例子的图。测量装置100例如具备第一光源101、透镜102、103、107、108、110、以及114~117、电流镜104以及105、棱镜106、半透明反射镜109、受光掩模112、光电倍增管113、以及第二光源118。例如,第一光源101是照射激光的激光源,第二光源118是照射红色等可见光的LED(Light Emitting Diode:发光二极管)。
由第一光源101照射的激光在透镜102以及103放大为平行光等,并射入电流镜104。射入电流镜104的激光在电流镜104的镜面被反射,并射入电流镜105。射入电流镜105的激光在电流镜105的镜面被反射,并射入棱镜106。射入棱镜106的激光经由透镜107聚光到眼球E的前房内。此时,通过未图示的驱动机构驱动电流镜104以及105,变更相对于激光的射入角度的镜面的角度。例如,电流镜104有助于前房内的激光的水平扫描(x方向),电流镜105有助于垂直扫描(y方向)。由此,使聚光到眼球E的前房内的激光以x方向以及y方向的二维方式扫描。
聚光到眼球E的前房内的激光的一部分被反射为散射光。例如,在前房内漂浮蛋白质分子或者细胞等生物体物质的情况下,激光的一部分被生物体物质反射。从前房反射的散射光在由透镜108聚光之后转换为平行光,并射入半透明反射镜109。射入半透明反射镜109的光的一部分透过半透明反射镜109,被导光至透镜110侧,射入半透明反射镜109的剩余的光被半透明反射镜109反射,被导光至透镜114侧。透过了半透明反射镜109的散射光在由透镜110聚光之后通过用于限定视场的受光掩模112,该通过的光射入光电倍增管13。受光掩模112例如是具有圆形的开口的针孔。由此,仅焦点所在的位置的光射入光电倍增管13。光电倍增管13通过光子计数法将射入的光转换为离散的电信号。然后,光电倍增管13将转换后的电信号输出给解析装置200。
被半透明反射镜109反射的散射光由透镜114进行漫射或者聚光,并经由透镜115以及116射入进行检查的利用者的眼睛。由第二光源118照射的光被透镜117漫射或者聚光之后通过半透明反射镜109射入利用者的眼睛。透镜117设置在使由第二光源118照射的光的像成像在与受光掩模112共轭的位置的位置。由此,例如,通过第二光源118的光将框状的窗口的像成像,从而能够使利用者视觉确认相对于窗口的反射散射光的生物物体的相对位置。
测量装置100通过使电流镜104以及105驱动使激光二维(x-y方向)地扫描,之后使受光掩模112的位置向通过了透镜110的光的光轴方向移动,来使激光的扫描位置向前房的进深方向亦即z方向移动。由此,能够三维地测量眼球E的前房。
另外,也可以测量装置100还具备使上述的光学系统的位置沿z方向移动的工作台。该情况下,测量装置100也可以在使电流镜104以及105驱动使激光二维(x-y方向)地扫描时,代替使受光掩模112的位置移动,而通过使工作台沿z方向移动,来变更激光的扫描位置。
[解析装置的构成]
图3是表示解析装置200的构成的一个例子的图。如图示的例子那样,解析装置200具备通信接口202、显示部204、存储部210、以及控制部220。
通信接口202例如通过有线或者无线与测量装置100进行通信。另外,通信接口202也可以与测量装置100以外的其它装置进行通信。
显示部204例如是LCD(Liquid Crystal Display:液晶显示器)、有机EL(Electroluminescence:电致发光)显示器等显示装置。
存储部210例如通过HDD(Hard Disc Drive:硬盘驱动器)、闪存、EEPROM(Electrically Erasable Programmable Read Only Memory:电可擦可编程只读存储器)、ROM(Read Only Memory:只读存储器)、RAM(Random Access Memory:随机存储器)等实现。存储部210例如储存在控制部220中参照的程序、体数据212、扫描位置信息214等。
图4是示意地表示体数据212的图。如图示那样,第一实施方式中的体数据212是将从散射光转换来的电信号的信号强度值即散射光的强度与将宽度方向x、高度方向y、进深方向z作为维的三维空间所包含的各体素建立对应关系的数据。以下,将与体素建立对应关系的物理量称为体素值进行说明。此外,在图示的例子中,将体素表现为具有某体积的三维空间,但也可以表现为无维(零维)的点。在由点表现体素的情况下,体数据212可以是具有储存了物理量作为阵列元素的多维的阵列结构的数据。
图5是表示扫描位置信息214的一个例子的图。扫描位置信息214是与后述的检测区域Ra的设定位置相关的信息,是将规定区域标志与体数据212示出的空间的各坐标建立对应关系的信息。规定区域标志是指利用二值表示各坐标的体素是否为规定区域的标志。在本实施方式中,将在眼球E的前房充满的眼房水中漂浮的生物体细胞作为观测对象,所以规定区域定义为生物体细胞存在的区域。
例如,在体数据212示出的空间设定了某检测区域Ra时,在判定为与该检测区域Ra重合的空间区域为规定区域的情况下,将标志“1”与该空间区域所包含的一个以上的体素的各坐标建立对应关系,否则,将标志“0”建立对应关系。此时,对于各体素的坐标的规定区域标志来说,在是否为规定区域的判定前预先与“0”建立对应关系。每当改变检测区域Ra的位置(使位置偏移)就进行该处理,最终被赋予一次标志“1”的体素(标志从“0”改写为“1”的体素)是作为规定区域的体素,图5所示的规定区域标志成为“1”。即,一旦规定区域标志从“0”变更为“1”,该标志保持为“1”的状态。
控制部220例如具备获取部222、检测区域设定部224、指标值运算部226、规定区域判定部228、聚类生成部230、评价判定部232、以及显示控制部234。这些构成要素的一部分或者全部通过CPU(Central Processing Unit:中央处理器)、GPU(Graphics ProcessingUnit:图形处理器)等处理器执行储存于存储部210的程序来实现。另外,控制部220的构成要素的一部分或者全部也可以通过LSI(Large Scale Integration:大规模集成电路)、ASIC(Application Specific Integrated Circuit:专用集成电路)、或者FPGA(Field-Programmable Gate Array:现场可编程门阵列)等硬件实现,也可以通过软件与硬件的配合实现。
以下,使用流程图对控制部220所进行的一系列的处理进行说明。图6以及图7是表示控制部220所进行的一系列的处理的一个例子的流程图。例如若经由通信接口202从测量装置100或者其它的装置获取电信号则执行本流程图的处理。另外,可以按规定的周期反复进行本流程图的处理。
首先,若获取部222经由通信接口202从测量装置100获取电信号,则基于该信号生成体数据212(步骤S100),并使其存储于存储部210。另外,也可以在测量装置100侧进行生成体数据212的处理,该情况下,获取部222也可以从测量装置100获取体数据212。另外,在云服务器等外部存储装置存储已经生成的体数据212的情况下,获取部222也可以经由通信接口202与外部存储装置进行通信,从该装置获取体数据212。
接下来,检测区域设定部224根据体数据212的分辨率(体素相对于体数据212的密度),决定检测区域Ra的尺寸。“决定尺寸”例如是指决定各维的区域的长度、区域的面积或者体积等。
例如,在体数据212示出的空间是由x-y-z表示的三维空间的情况下,对于检测区域设定部224来说,体数据212的分辨率越大越减小检测区域Ra的宽度、高度、进深的各边的长度,体数据212的分辨率越小越增大检测区域Ra的宽度、高度、进深的各边的长度。另外,检测区域设定部224也可以以成为体数据212的单位区域亦即体素的整数倍的体积的方式决定检测区域Ra的尺寸。另外,检测区域设定部224也可以不依赖于体数据212的分辨率地决定宽度、高度、进深中一部分的边的长度。此外,检测区域Ra的形状(各边的长度的比率)被预先决定。
接下来,检测区域设定部224在体数据212示出的空间的任意的位置设定决定的尺寸的检测区域Ra(步骤S104)。
接下来,检测区域设定部224以设置某恒定的间隙包围设定的检测区域Ra的方式,设定一个以上的对照区域Rb(步骤S106)。“包围”是指至少在设定了检测区域Ra的空间的部分空间中,检测区域Ra落入由一个或者多个对照区域Rb形成的封闭空间中。例如,在检测区域Ra为立方体等六面体的情况下,在六面体的某一个剖面落入由一个或者多个对照区域Rb形成的封闭空间中的情况下,对照区域Rb包围检测区域Ra。
图8A~图8F是表示检测区域Ra的设定位置的一个例子的图。图8A~图8E表示作为体数据212示出的x-y-z空间的部分空间的x-y平面。图8A表示z=1的x-y平面,图8B表示z=2的x-y平面,图8C表示z=3的x-y平面,图8D表示z=4的x-y平面,图8E表示z=5的x-y平面。另外,在图示的例子中,各体素表现为x、y、z的各边的长度为1的立方体。此外,这里称为“平面”的情况下,与一个体素相对应的z方向的长度可以是单位长度(例如1)。
例如,在检测区域设定部224将与一个体素相同的体积的检测区域Ra设定为(x,y,z)=(3,3,3)的坐标的情况下,在z=3的x-y平面,将与检测区域Ra相邻的八个体素作为间隙,并将与这八个体素相邻的十六个体素设定为对照区域Rb。另外,检测区域设定部224在z方向上与z=3的x-y平面相邻的z=2的x-y平面、和z=4的x-y平面双方,将(x,y)=(3,3)的坐标的体素、和与该体素相邻的八个体素作为间隙,并将与这八个体素相邻的十六个体素设定为对照区域Rb。另外,检测区域设定部224在z方向上与z=2的x-y平面相邻的z=1的x-y平面、和在z方向上与z=4的x-y平面相邻的z=5的x-y平面双方,将全部的体素设定为对照区域Rb。由此,如图8F所示,在三维空间中,至少隔开检测区域Ra的周围的一个体素设定对照区域Rb。
另外,检测区域设定部224也可以也在对照区域Rb彼此之间设置间隙。即,检测区域设定部224也可以使多个对照区域Rb相互分离地设定。图9A~图9E是表示检测区域Ra的设定位置的其它的例子的图。例如,检测区域设定部224在将与一个体素相同的体积的检测区域Ra设定为(x,y,z)=(3,3,3)的坐标的情况下,在z=3的x-y平面,将与检测区域Ra相邻的八个体素作为间隙,并将与这八个体素相邻的十六个体素中的(x,y)=(1,1)、(1,3)、(1,5)、(3,1)、(3,5)、(5,1)、(5,3)、(5,5)八个体素设定为对照区域Rb。另外,检测区域设定部224在z方向上与z=3的x-y平面相邻的z=2的x-y平面、和z=4的x-y平面双方不设定对照区域Rb,并在z方向上与z=2的x-y平面相邻的z=1的x-y平面、和在z方向上与z=4的x-y平面相邻的z=5的x-y平面双方,将(x,y)=(1,1)、(1,3)、(1,5)、(3,1)、(3,3)、(3,5)、(5,1)、(5,3)、(5,5)九个体素设定为对照区域Rb。这样,在三维空间中,可以隔开检测区域Ra的周围的一个体素设定对照区域Rb,并且在对照区域Rb彼此之间也设置间隙。
另外,检测区域设定部224也可以设定与多个体素相同的体积的检测区域Ra。图10是表示检测区域Ra的设定位置的其它的例子的图。如图示的例子那样,检测区域设定部224也可以设定与九个体素相同的体积的检测区域Ra。该情况下,检测区域设定部224可以也使在与检测区域Ra之间设置间隙而设定的对照区域Rb为将九个体素加在一起的体积。
另外,检测区域设定部224并不限定于设置与一个体素相对应的间隙,也可以设置两个、三个或者其以上的间隙。图11以及图12是表示设置与多个体素相对应的间隙而设定的对照区域Rb的一个例子的图。在图11的例子中,设置两个体素的间隙,在图12的例子中,设置三个体素的间隙。
另外,虽然以检测区域设定部224在检测区域Ra的周围设定对照区域Rb时,将对照区域Rb配置为矩形形状为例进行了说明但是并不限定于此,也可以配置为圆环状。图13以及图14是表示配置为圆环状的对照区域Rb的一个例子的图。检测区域设定部224可以如图13的例子那样,在对检测区域Ra设置间隙的同时将八个对照区域Rb配置为圆环状。另外,检测区域设定部224也可以如图14的例子那样,以对照区域Rb彼此相互重合一部分的方式设定多个对照区域Rb。由此,在检测区域Ra的周围设定了外观上其区域连续的一个对照区域Rb。
另外,虽然以检测区域设定部224在检测区域Ra的周围设定多个对照区域Rb为例进行了说明但是并不限定于此,也可以设定包围检测区域Ra的一个对照区域Rb。
另外虽然在图10~图12、以及图14所示的例子中,检测区域设定部224在对照区域Rb的内侧设置了间隙,但也可以如图13所例示的那样,以向对照区域Rb的外侧突出一定程度的方式设定间隙。
图15是表示被单个对照区域Rb包围的检测区域Ra的一个例子的图。如图示那样,例如,检测区域设定部224可以设定球形的检测区域Ra,并设定在其周围设置了间隙的相同的球形的对照区域Rb。例如,检测区域设定部224可以使对照区域Rb的内周为对检测区域Ra的外周考虑了间隙的厚度的长度。由此,能够以从360°的全部方向包住检测区域Ra的方式设定对照区域Rb。
另外,检测区域设定部224也可以在设定球形的检测区域Ra的情况下,设定图16所示那样的对照区域Rb。图16是表示对照区域Rb的设定方法的其它的例子的图。如图示那样,在检测区域Ra为球形的情况下,可以在球形的检测区域Ra的z方向上的中心轴与检测区域Ra相交的各交点P1以及P2分别设定以三片弧形的区域为一组的对照区域Rb。即,检测区域设定部224可以设定将球形的上侧划分成3部分,并将下侧划分成3部分的对照区域Rb。此时,在从球形的检测区域Ra的中心到间隙的外周面为止的距离为A[m]的情况下,各对照区域Rb的内周面的弧的长度可以至少为π/3[rad]×A[m]。
接下来,指标值运算部226计算由检测区域设定部224设定的检测区域Ra的最大体素值A(步骤S108)。最大体素值A是指与检测区域Ra所包含的一个以上的体素建立对应关系的体素值中最大的体素值。
接下来,规定区域判定部228判定最大体素值A是否在预先决定的阈值亦即最大体素阈值THa以上(步骤S110)。规定区域判定部228在判定为最大体素值A小于最大体素阈值THa的情况下,判定为检测区域Ra所包含的一个以上的体素全部不为规定区域(步骤S112)。最大体素阈值THa是“第二阈值”的一个例子。
另一方面,在判定为最大体素值A在最大体素阈值THa以上的情况下,指标值运算部226计算检测区域Ra的平均体素值B(步骤S114)。平均体素值B是与检测区域Ra所包含的一个以上的体素建立对应关系的体素值的平均。平均体素值B是“第一平均值”的一个例子。
接下来,规定区域判定部228判定平均体素值B是否在预先决定的阈值亦即平均体素阈值THb以上(步骤S116)。平均体素阈值THb例如是比最大体素阈值THa小的值。规定区域判定部228在判定为小于平均体素阈值THb的情况下,移至S112的处理,判定为检测区域Ra所包含的一个以上的体素全部不为规定区域。平均体素阈值THb是“第三阈值”的一个例子。
另一方面,在判定为平均体素值B在平均体素阈值THb以上的情况下,指标值运算部226计算通过S106的处理设定的一个以上的对照区域Rb各自的平均体素值C(步骤S118)。平均体素值C是与对照区域Rb所包含的一个以上的体素建立对应关系的体素值的平均。平均体素值C是“第二平均值”的一个例子。
接下来,规定区域判定部228判定从各平均体素值B减去平均体素值C后的值(B-C)是否在差分体素阈值THc以上(步骤S120)。差分体素阈值THc例如是比平均体素阈值THb小的值。差分体素阈值THc是“第一阈值”的一个例子。
例如,规定区域判定部228在判定为即使从各平均体素值B减去平均体素值C后的各差值中的任意一个小于差分体素阈值THc的情况下,也移至S112的处理,判定为检测区域Ra所包含的一个以上的体素全部不为规定区域。
另一方面,规定区域判定部228在判定为从各平均体素值B减去平均体素值C后的各差值中的任意一个均在差分体素阈值THc以上的情况下,判定为检测区域Ra所包含的一个以上的体素全部为规定区域(步骤S122)。
接下来,规定区域判定部228更新扫描位置信息214中,检测区域Ra所包含的一个以上的体素各自的坐标所对应的规定区域标志(步骤S124)。
例如,规定区域判定部228在S112的处理中判定为检测区域Ra的体素不为规定区域的情况下,在扫描位置信息214中,使该检测区域Ra的体素的规定区域标志为“0”。另外,规定区域判定部228在S122的处理中判定为检测区域Ra的体素为规定区域的情况下,在扫描位置信息214中,使该检测区域Ra的体素的规定区域标志为“1”。
图17以及图18是用于说明在以包围检测区域Ra的方式设定对照区域Rb时,在设置间隙时和不设置间隙时判定结果不同的图。例如,有作为检测对象的生物体细胞的x方向的大小与检测区域Ra以及对照区域Rb双方的x方向的大小同等程度的情况。这样的情况下,在如图17所示,不在检测区域Ra与对照区域Rb之间设置间隙的情况下,有作为检测对象的生物体细胞与检测区域Ra和对照区域Rb双方重叠的可能性,而从检测区域Ra的平均体素值B减去对照区域Rb的平均体素值C后的值(B-C)容易小于差分体素阈值THc,而有可能本来应该判定为检测区域Ra是表示生物体细胞的规定区域,但判定为检测区域Ra不是表示生物体细胞的规定区域。与此相对,在如图18所示,在检测区域Ra与对照区域Rb之间设置了间隙的情况下,即使是生物体细胞较大,也与检测区域Ra相邻的区域重叠这样的情况,由于该重叠的区域是作为间隙设置的区域,所以从检测区域Ra的平均体素值B减去对照区域Rb的平均体素值C后的值(B-C)也容易在差分体素阈值THc以上,能够精度良好地检测表示生物体细胞的规定区域。
接下来,检测区域设定部224判定是否对体数据212示出的空间的全部区域设定了检测区域Ra(步骤S126)。检测区域设定部224在判定为未对全部区域设定检测区域Ra的情况下,返回到S104的处理,变更检测区域Ra的设定位置。由此,重复进行新变更了位置的检测区域Ra的体素为是否规定区域的判定。
另一方面,在判定为对全部区域设定了检测区域Ra的情况下,聚类生成部230参照扫描位置信息214,生成将规定区域标志为“1”的体素中平均体素值B同等程度的体素聚类后的体素聚类CL(步骤S128)。“体素值同等程度”例如是指在对作为比较对象的体素值允许了几[%]~十几[%]左右的误差的范围内相同。因此,“体素值同等程度”也包含有体素值相同的情况。
图19~图21是表示体素聚类CL的生成方法的一个例子的图。例如,如图19以及图20所示的例子那样,若设定多个检测区域Ra,则聚类生成部230通过标记处理,将多个检测区域Ra中平均体素值B同等程度的区域彼此结合,生成一个体素聚类CL。
例如,聚类生成部230选择体数据212所包含的多个体素中还未赋予标签的体素作为关注的体素(以下,称为关注体素),并对该关注体素赋予某标签。聚类生成部230对赋予了标签的关注体素的周围的体素(例如,在x、y、z方向上与关注体素相邻的体素),判定是否已经赋予了标签。若未对周围的体素赋予标签,并且为与关注体素同等程度的平均体素值B,则聚类生成部230对周围的体素赋予与对关注体素赋予的标签相同的标签。此时,属于相同的体素聚类CL的体素作为相同的体素值(例如平均体素值B)进行处理,所以赋予相同的标签。聚类生成部230将赋予了标签的周围的体素作为新的关注体素进行处理,进一步对周围的体素进行标签的赋予的有无和体素值的确认并赋予标签。这样一来,聚类生成部230对体数据212的全部体素赋予标签,并将赋予了相同的标签并且相互处于相邻关系的体素的集合生成为一个体素聚类CL。
接下来,聚类生成部230判定体素聚类CL的大小是否在允许尺寸内(步骤S130)。允许尺寸是指根据检查系统1的检查对象设定的阈值,例如,在检查对象为眼球E内的生物体细胞等的情况下,是根据测量装置100的分辨率对假定存在于眼球E内的最大的生物体细胞的实际尺寸进行缩放而成的预计的尺寸。
例如,如上述的图21所示,聚类生成部230对体素聚类CL的x方向的最大尺寸ΔLX与x方向上的允许尺寸THΔLX进行比较,对体素聚类CL的z方向的最大尺寸ΔLZ与z方向上的允许尺寸THΔLZ进行比较。同样地,聚类生成部230对体素聚类CL的y方向的最大尺寸ΔLY与y方向上的允许尺寸THΔLY进行比较。聚类生成部230在最大尺寸ΔLX在THΔLX以下、最大尺寸ΔLZ在THΔLZ以下、并且最大尺寸ΔLY在THΔLY以下的情况下,判定为体素聚类CL的大小在允许尺寸内。
另一方面,聚类生成部230在最大尺寸ΔLX超过THΔLX的情况下、最大尺寸ΔLZ超过THΔLZ的情况下、或者最大尺寸ΔLY超过THΔLY的情况下,判定为体素聚类CL的大小不在允许尺寸内。此外,也可以即使最大尺寸ΔLX、ΔLZ、ΔLY中的任意一个在允许尺寸以下,聚类生成部230也判定为体素聚类CL的大小在允许尺寸内。
若体素聚类CL的大小不在允许尺寸内,则聚类生成部230判断为该体素聚类CL示出的规定区域是与提取对象物(例如生物体细胞)不同的对象物(例如虹膜等),而将与成为体素聚类CL的基础的多个体素建立对应关系的规定区域标志从“1”变更为“0”(步骤S132)。由此,放弃对成为体素聚类CL的基础的多个体素做出的是规定区域这样的判定结果,将该多个体素作为不是规定区域的体素进行处理。另一方面,若体素聚类CL的大小在允许尺寸内,则将该体素聚类CL检测为生物体细胞。
接下来,聚类生成部230计算未放弃是规定区域这样的判定结果的体素聚类CL的数目,即检测出的生物体细胞的数目(步骤S134)。
接下来,控制部220判定是否对体数据212的全部区域进行了S128~S134的一系列的聚类处理(步骤S136)。在判定为未对体数据212的全部区域进行上述的一个例子的聚类处理的情况下,控制部220返回到S128的处理。
另一方面,在判定为对体数据212的全部区域进行了上述的一个例子的聚类处理的情况下,评价判定部232基于由聚类生成部230计算出的生物体细胞的数目,进行被检体亦即眼球E的评价(步骤S138)。
例如,评价判定部232参照各体数据212的生物体细胞的数目,在生物体细胞的数目比基准值大的情况下,可以给出“有受检者患患有特定的疾病(例如角膜炎症)的可能性”这样的评价结果。例如,基准值可以根据观测到的生物体细胞的数目与特定的疾病的发症的有无的相关结果等适当地决定。
接下来,显示控制部234使基于评价判定部232的评价结果的图像显示于显示部204(步骤S140)。由此,本流程图的处理结束。
图22是表示显示了基于评价判定部232的评价结果的图像的显示部204的一个例子的图。例如,显示控制部234可以控制显示部204,使其显示包含生物体细胞的数目、生物体细胞的数目的基准值、以及特定的疾病的有无(或者其概率)的评价结果。另外,显示控制部234也可以在体数据212的一个剖面(例如,生物体细胞的数目最多的剖面图像),使提取出的生物体细胞与其提取位置重叠显示。通过使这样的图像显示于显示部204,能够减少依赖于各阅读者(例如医师等)的经验、技能等的影响,即个体偏差的影响,能够辅助医师等基于定量并且客观的数值给出稳定的诊断结果。此外,显示控制部234也可以仅使将由聚类生成部230计算出的生物体细胞的数目与作为其计算源的体数据212的一个剖面建立相关关系的图像显示于显示部204。由此,至少能够客观地对被检体评价生物体细胞的数目。
根据以上说明的第一实施方式,通过具备获取将物理量与划分三维空间、三维时空间、或者四维时空间而成的各体素建立对应关系的体数据212的获取部222、在体数据212示出的空间或者时空间中,设定作为三维以上的时空间的检测区域Ra,并在与检测区域Ra的维相同的空间中,设置包围检测区域Ra的间隙,并在包围该间隙的位置设定对照区域Rb的检测区域设定部224、以及基于相互比较检测区域Ra以及对照区域Rb各自所包含的一个以上的体素的物理量,判定检测区域Ra所包含的一个以上的体素是否是规定区域的规定区域判定部228,能够提高对象物的检测精度。
(第二实施方式)
以下,对第二实施方式进行说明。在第二实施方式中,在基于扫描时刻(观测时刻)相互不同的多个体数据212,进行作为被检体的眼球E的评价这一点,与第一实施方式不同。以下,以与第一实施方式的不同点为中心进行说明,并省略与第一实施方式相同的部分的说明。
在第二实施方式的检查系统1中,例如,通过每经过规定的时间观测注入了对身体无害的微小粒子(例如,作为标记物质的微珠等)的眼球E的眼房水,获取按时间序列连续的多个体数据212。即,在第二实施方式中,获取进一步对x-y-z的三维添加了时间t的维后的四维的体数据212。
例如,在第二实施方式中,将微小粒子存在的区域作为规定区域进行处理。该情况下,评价判定部232基于通过聚类生成部230的聚类处理生成的体素聚类CL中,大小在允许尺寸内的体素聚类CL(即微小粒子)的与时间经过对应的移动量,评价作为被检体的眼球E。
一般而言,可以说前眼房水流动,若该眼房水的流动性较低则诊断为青光眼的嫌疑较高。因此,评价判定部232根据注入前房内的微小粒子的与时间经过对应的移动,判定青光眼等病状的有无。
图23是示意地表示眼球E的剖面的图。一般而言,眼房水在图中的箭头所示的流路中流动,例如从流动睫状体上皮向后眼房分泌,在润湿晶状体之后,通过瞳孔流至前眼房,从虹膜附着部与角膜所成的角的部分(前房角)通过小梁网的间隙(丰塔纳腔),集中于巩膜静脉窦,并流出到前睫状体静脉。另外,也有眼房水的一部分从前眼房经由虹膜的基部被吸收的情况。这样,眼房水通过在眼球E内流动,参与晶状体、虹膜、角膜等的营养,或者将眼压保持恒定来保持眼球E的形状。然而,若眼房水的排出功能等发生异常,则有眼房水过度地积蓄在前眼房内的情况。该情况下,成为眼压异常高的状态,容易引发青光眼。
这样,由于在眼房水的流动性与青光眼之间存在因果关系,所以评价判定部232根据四维的体数据212观测微小粒子的移动,评价眼房水的流动性。
图24A以及图24B是表示观测时刻不同的两个体数据212的一个剖面的图。图24A示出在某一时刻ti,在x-y平面上观测到的微小粒子,图24B示出在从某一时刻ti开始进一步经过了规定时间ΔT后的时刻ti+ΔT,在x-y平面上观测到的微小粒子。在对图24A与图24B进行比较的情况下,在时刻ti+ΔT,微小粒子从时刻ti的位置沿x轴方向移动Δx的位移量。例如,评价判定部232判定每规定时间ΔT的微小粒子的位移量是否小于阈值,在微小粒子的位移量小于阈值的情况下,判定为该微小粒子留在相同的位置,在微小粒子的位移量在阈值以上的情况下,判定为该微小粒子未留在相同的位置而流动。然后,显示控制部234使基于评价判定部232的评价结果的图像显示于显示部204。由此,阅读者能够根据眼房水的流动的有无、其流动的大小,定量并且客观地给出青光眼等病状的诊断结果。
根据以上说明的第二实施方式,基于四维以上的体数据212,求出作为规定区域检测出的对象物的与时间经过对应的移动的有无、其移动量,所以能够更定量并且客观地评价作为被检体的眼球E。
(第三实施方式)
以下,对第三实施方式进行说明。在第一以及第二实施方式中,以作为生成体数据212的装置,具备向前眼房照射激光并测量该前眼房内的生物体细胞的测量装置100为例进行了说明。与此相对,在第三实施方式中,在作为生成体数据212的装置,代替测量装置100,具备光干涉断层图像拍摄装置(以下,称为OCT拍摄装置)这一点,与第一以及第二实施方式不同。以下,以与第一以及第二实施方式的不同点为中心进行说明,省略与第一以及第二实施方式相同的部分的说明。
图25是表示第三实施方式的检查系统2的构成的一个例子的图。第三实施方式的检查系统(眼底拍摄系统)2例如具备解析装置200和OCT拍摄装置300。OCT拍摄装置300通过对眼球E照射光,并测量该光的反射光与照射的一部分的光相互干涉而成的干涉光,来测定眼球E内部的位移。由此,例如,获取拍摄到包含视网膜的眼底Er的OCT图像IM作为体数据212。以下,OCT拍摄装置300作为光谱域OCT(Spectral-domain OCT;SD-OCT)、波长扫描OCT(Swept-source OCT;SS-OCT)等傅里叶域OCT(Fourier-domain OCT;FD-OCT)进行说明,但并不限定于此。OCT拍摄装置300例如也可以采用时域OCT(Time-domain OCT;TD-OCT)或者其它的方式。
例如,第三实施方式的解析装置200从由OCT拍摄装置300生成的OCT图像IM,提取HRF(Hyperreflective foci:高反射信号)等高亮度部位。例如,可以说HRF表示脂蛋白质或者巨噬细胞,近年来,进行了根据患有糖尿病视网膜症的患者的眼球的OCT图像观测到的HRF的数目与该患者的视功能相关这样的学术报告。在第三实施方式中,规定区域被定义为HRF等高亮度部位存在的区域。
照明光学系统320例如具备卤素灯或者氙气灯等照明用的光源(未图示),通过将从该光源照射的光导光至眼底Er来照亮眼底Er。
拍摄光学系统330经由光纤Fa将在眼底Er反射的反射光导光至OCT单元310侧。另外,拍摄光学系统330使从OCT单元310经由光纤Fa导光的照射光扫描并且导光至眼球E。例如,拍摄光学系统330具备准直仪、电流镜(均未图示)等,将对眼球E照射的照射光的照射方向(图中z方向)在与该照射方向正交的水平方向(图中x方向或者y方向)上变更。即,拍摄光学系统330通过光栅扫描方式使照射光扫描。由此,照射到眼球E的照射光在x方向以及y方向上扫描。
解析装置200的获取部222经由通信接口202从OCT拍摄装置300获取OCT图像IM,并使获取到的OCT图像IM作为体数据212存储于存储部210。
图26是表示通过OCT拍摄装置300生成的体数据212的一个例子的图。如图示那样,第三实施方式的体数据212由多个OCT图像IMn(n是任意的自然数)构成。各OCT图像IMj(1≤j≤n)的z方向表示沿着照射光的照射方向(照射光的光轴)的方向,表示眼球E的深度方向。沿着照射光的照射方向的方向例如是指相对于照射光的光轴允许几度~十几度左右的误差(角度宽度)的方向。另外,各OCT图像IMj的x方向表示与z方向正交的平面的任意一个方向。这些OCT图像IMn在与z方向以及x方向双方正交的y方向上并排。y方向与各OCT图像IM的拍摄时刻t对应。即,各OCT图像IM按拍摄时刻依次排列。
另外,获取部222也可以从OCT拍摄装置300一个一个获取OCT图像IM。该情况下,获取部222可以通过参照各OCT图像IM的拍摄时刻t、其拍摄时的y方向上的位置信息,在y方向上排列依次获取的OCT图像IM。此时,获取部222例如可以以成为时间序列顺序的方式适当地在y方向上对OCT图像IM进行排序。
第三实施方式的解析装置200基于由这样的多个OCT图像IM构成的体数据212,判定检测区域Ra所包含的体素是否为规定区域,在存在作为规定区域的体素的情况下,判定将该体素聚类后的体素聚类CL是否在允许尺寸内。解析装置200在体素聚类CL在允许尺寸内的情况下,判定为体素聚类CL是HRF等高亮度部位。即,解析装置200判定为在眼底Er存在HRF等高亮度部位。解析装置200计算在体数据212示出的空间(三维图像)中存在几个高亮度部位,并基于计算出的高亮度部位的数目评价作为被检体的眼球E。
根据以上说明的第三实施方式,在由多个OCT图像IM构成的体数据212示出的空间中,设定检测区域Ra、和包围该区域的一个以上的对照区域Rb,并基于相互比较检测区域Ra以及对照区域Rb各自所包含的一个以上的体素的物理量(例如像素值),判定检测区域Ra所包含的一个以上的体素是否为规定区域(HRF等高亮度部位),所以与第一或者第二实施方式相同,能够提高对象物的检测精度。
(其它的实施方式)
以下,对其它的实施方式进行说明。在上述的第一以及第二实施方式中,以检测存在于眼球E的前房内的生物体细胞为例进行说明,在第三实施方式中,以检测存在于眼球E的眼底Er的HRF等高亮度部位为例进行了说明但并不限定于此。例如,也可以使用上述的方法检测对内脏等进行荧光标记的微小物体、血液中的血球等。该情况下,可以利用荧光显微镜等光学显微镜、TEM或SEM等电子显微镜、CT、PET、MRI、NMR、或者超声波等其它的剖面图像拍摄装置等各种测量装置获取体数据212。
例如,在生物成像的领域,通过使内脏等的细胞与为了检测特定的抗体而进行标记的物质(以下,称为标记物质)接触,使标记物质与特定的抗体结合,来进行标记细胞。该情况下,利用荧光显微镜等生成测量了内脏等的细胞的三维的体数据212。解析装置200可以参照该体数据212,检测标记物质存在的区域作为规定区域。标记物质例如是通过FITC(Fluorescein Isothiocyanate:异硫氰酸荧光素)、Alexa Fluor色素、Cy色素等荧光色素、PE(Phycoerythrin藻红蛋白)、APC(Allophycocyanin:别藻蓝蛋白)等荧光蛋白质使其染色的抗体等。另外,标记物质也可以是使酵素结合的抗体,也可以是与磁珠、琼脂糖等载体结合的抗体。通过从体数据212检测这样的标记物质存在的区域作为规定区域,例如能够辨别标记物质集中于内脏的哪个位置(将内脏的哪个位置染色)。
另外,在体数据212是利用CT血管拍摄法(CT血管造影;CTA)、MR血管拍摄法(MR血管造影;MRA)等三维血管拍摄法观测血管内部的三维数据的情况下,解析装置200既可以从血管的体数据212检测特定的血球存在的区域作为规定区域,也可以检测利用荧光等进行了标记的特定的血球存在的区域作为规定区域。
这样在想要检测的对象为微小物体的情况下,能够应用本方法。
[硬件构成]
例如,通过图27所示那样的硬件的构成实现上述的实施方式的解析装置200。图27是表示实施方式的解析装置200的硬件构成的一个例子的图。
解析装置200构成为通过内部总线或者专用通信线相互连接NIC(NetworkInterface Card:网卡)等通信接口200-1、CPU200-2、RAM200-3、ROM200-4、闪存或者HDD等二次存储装置200-5、以及驱动装置200-6。在驱动装置200-6安装有光盘等移动存储介质。储存于二次存储装置200-5的程序200-5a通过DMA控制器(未图示)等展开于RAM200-3,并被CPU200-2执行,从而实现控制部220。另外,CPU200-2参照的程序既可以储存于安装于驱动装置200-6的移动存储介质,也可以经由网络从其它的装置下载。
上述实施方式能够如以下那样表现。
一种解析装置,构成为具备:
存储器,储存程序;以及
处理器,
上述处理器通过执行上述程序进行以下处理:
获取将物理量与每个单位区域建立对应关系的数据,其中,上述单位区域是划分三维空间、三维时空间、或者四维时空间而成的,
在述空间或者时空间中,设定作为三维以上的时空间的检测区域,
在与上述检测区域的维相同的空间中,设置包围上述检测区域的间隙,并在包围上述间隙的位置设定对照区域,
基于相互比较上述设定的上述检测区域以及上述对照区域各自所包含的一个以上的上述单位区域的物理量,判定上述检测区域所包含的一个以上的上述单位区域是否为规定区域。
以上,使用实施方式对用于实施本发明的方式进行了说明,但本发明并不限定于这样的实施方式,能够在不脱离本发明的主旨的范围内施加各种变形以及置换。
Claims (8)
1.一种存储介质,存储用于使计算机执行的程序,该程序包括:
获取将物理量与每个单位区域建立对应关系的数据,其中,上述单位区域是划分三维空间、或者三维时空间、或者四维时空间而成的;
在上述空间或者时空间中,设定作为三维以上的时空间的检测区域;
在与上述检测区域的维相同的空间中,设置包围上述检测区域的间隙,并在包围上述间隙的位置设定对照区域;
基于相互比较上述设定的上述检测区域以及上述对照区域各自所包含的一个以上的上述单位区域的物理量,判定上述检测区域所包含的一个以上的上述单位区域是否为规定区域;
判定上述检测区域所包含的一个以上的上述单位区域的物理量的最大值是否在第二阈值以上;以及
在判定为上述最大值小于上述第二阈值的情况下,判定为上述检测区域所包含的一个以上的上述单位区域不为上述规定区域。
2.根据权利要求1所述的存储介质,其中,上述程序包括:
在包围上述间隙的位置设定多个上述对照区域。
3.根据权利要求2所述的存储介质,其中,上述程序包括:
在设定多个上述对照区域的情况下,使多个上述对照区域相互分离地设定。
4.根据权利要求2或者3所述的存储介质,其中,上述程序包括:
在设定多个上述对照区域的情况下,使多个上述对照区域相互相邻地设定。
5.根据权利要求2或者3所述的存储介质,其中,上述程序包括:
在设定多个上述对照区域的情况下,使多个上述对照区域相互重叠一部分地设定。
6.根据权利要求1或者2所述的存储介质,其中,上述程序包括:
判定上述检测区域所包含的一个以上的上述单位区域的物理量的第一平均值与一个以上的上述对照区域各自所包含的一个以上的上述单位区域的物理量的第二平均值的差分是否在第一阈值以上;以及
在判定为上述第一平均值与一个以上的上述第二平均值的一部分或者全部的差分在上述第一阈值以上的情况下,判定为上述检测区域所包含的一个以上的上述单位区域为上述规定区域。
7.一种存储介质,存储用于使计算机执行的程序,该程序包括:
获取将物理量与每个单位区域建立对应关系的数据,其中,上述单位区域是划分三维空间、或者三维时空间、或者四维时空间而成的;
在上述空间或者时空间中,设定作为三维以上的时空间的检测区域;
在与上述检测区域的维相同的空间中,设置包围上述检测区域的间隙,并在包围上述间隙的位置设定对照区域;
基于相互比较上述设定的上述检测区域以及上述对照区域各自所包含的一个以上的上述单位区域的物理量,判定上述检测区域所包含的一个以上的上述单位区域是否为规定区域;
判定上述检测区域所包含的一个以上的上述单位区域的物理量的第一平均值是否在第三阈值以上;以及
在判定为上述第一平均值不在上述第三阈值以上的情况下,判定为上述检测区域所包含的一个以上的上述单位区域不为上述规定区域。
8.一种存储介质,存储用于使计算机执行的程序,该程序包括:
获取将物理量与每个单位区域建立对应关系的数据,其中,上述单位区域是划分三维空间、或者三维时空间、或者四维时空间而成的;
在上述空间或者时空间中,设定作为三维以上的时空间的检测区域;
在与上述检测区域的维相同的空间中,设置包围上述检测区域的间隙,并在包围上述间隙的位置设定对照区域;
基于相互比较上述设定的上述检测区域以及上述对照区域各自所包含的一个以上的上述单位区域的物理量,判定上述检测区域所包含的一个以上的上述单位区域是否为规定区域;
生成使判定为是上述规定区域的一个以上的上述单位区域成为一个的聚类;以及
使上述生成的聚类的数目与上述获取的数据相关联地显示于显示部。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017209546 | 2017-10-30 | ||
JP2017-209546 | 2017-10-30 | ||
PCT/JP2018/040107 WO2019088020A1 (ja) | 2017-10-30 | 2018-10-29 | 解析装置、解析方法、およびプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111246792A CN111246792A (zh) | 2020-06-05 |
CN111246792B true CN111246792B (zh) | 2022-09-13 |
Family
ID=66333123
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201880067903.8A Active CN111246792B (zh) | 2017-10-30 | 2018-10-29 | 解析装置、解析方法、以及程序 |
Country Status (5)
Country | Link |
---|---|
US (1) | US11328417B2 (zh) |
EP (1) | EP3705027A4 (zh) |
JP (1) | JP7234126B2 (zh) |
CN (1) | CN111246792B (zh) |
WO (1) | WO2019088020A1 (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010514530A (ja) * | 2007-01-04 | 2010-05-06 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 獲得投影データから関心領域の補正画像を生成する装置、方法、及びコンピュータ・プログラム |
JP2015516192A (ja) * | 2012-03-20 | 2015-06-11 | ポール リード スミス ギターズ、リミテッド パートナーシップ | 画像内のデータを明らかにするための近傍関係に基づく画素の得点化及び調整 |
WO2016161198A1 (en) * | 2015-04-02 | 2016-10-06 | Hedronx Inc. | Virtual three-dimensional model generation based on virtual hexahedron models |
JP2016200444A (ja) * | 2015-04-08 | 2016-12-01 | 一般財団法人電力中央研究所 | ボイドの計測方法、計測装置、及び計測プログラム |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3608854B2 (ja) | 1995-09-22 | 2005-01-12 | 興和株式会社 | 眼科測定装置 |
US20160081663A1 (en) * | 2014-09-18 | 2016-03-24 | General Electric Company | Method and system for automated detection and measurement of a target structure |
-
2018
- 2018-10-29 JP JP2019550372A patent/JP7234126B2/ja active Active
- 2018-10-29 WO PCT/JP2018/040107 patent/WO2019088020A1/ja unknown
- 2018-10-29 EP EP18873441.2A patent/EP3705027A4/en not_active Withdrawn
- 2018-10-29 US US16/756,928 patent/US11328417B2/en active Active
- 2018-10-29 CN CN201880067903.8A patent/CN111246792B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010514530A (ja) * | 2007-01-04 | 2010-05-06 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 獲得投影データから関心領域の補正画像を生成する装置、方法、及びコンピュータ・プログラム |
JP2015516192A (ja) * | 2012-03-20 | 2015-06-11 | ポール リード スミス ギターズ、リミテッド パートナーシップ | 画像内のデータを明らかにするための近傍関係に基づく画素の得点化及び調整 |
WO2016161198A1 (en) * | 2015-04-02 | 2016-10-06 | Hedronx Inc. | Virtual three-dimensional model generation based on virtual hexahedron models |
JP2016200444A (ja) * | 2015-04-08 | 2016-12-01 | 一般財団法人電力中央研究所 | ボイドの計測方法、計測装置、及び計測プログラム |
Non-Patent Citations (1)
Title |
---|
Quantification of Calcified Particles in Human Valve Tissue Reveals Asymmetry of Calcific Aortic Valve Disease Development;Katsumi Yabusaki et al.;《Frontiers in Cardiovascular Medicine》;20161104;第3卷(第1期);第1-10页 * |
Also Published As
Publication number | Publication date |
---|---|
WO2019088020A1 (ja) | 2019-05-09 |
EP3705027A1 (en) | 2020-09-09 |
CN111246792A (zh) | 2020-06-05 |
EP3705027A4 (en) | 2021-09-15 |
US11328417B2 (en) | 2022-05-10 |
JP7234126B2 (ja) | 2023-03-07 |
US20200250818A1 (en) | 2020-08-06 |
JPWO2019088020A1 (ja) | 2020-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11935241B2 (en) | Image processing apparatus, image processing method and computer-readable medium for improving image quality | |
Sampson et al. | Towards standardizing retinal optical coherence tomography angiography: a review | |
JP7193343B2 (ja) | 機械学習技法を用いたoctアンギオグラフィにおけるアーチファクトを減少させるための方法及び装置 | |
US8684528B2 (en) | Fundus analyzing appartus and fundus analyzing method | |
JP6702764B2 (ja) | 光干渉断層データの処理方法、該方法を実行するためのプログラム、及び処理装置 | |
JP2016002380A (ja) | 画像処理装置、その作動方法及びプログラム | |
JP5628636B2 (ja) | 眼底画像処理装置及び眼底観察装置 | |
CN106028912B (zh) | 具有改进的运动对比的光学相干断层扫描(oct)系统 | |
US9060711B2 (en) | Automated detection of uveitis using optical coherence tomography | |
JP2013176497A (ja) | 眼底観察装置及び眼底画像解析装置 | |
CN115443480A (zh) | 使用通道编码板坯进行oct正面病理学分割 | |
Chu et al. | Complex signal-based optical coherence tomography angiography enables in vivo visualization of choriocapillaris in human choroid | |
JP2020103579A (ja) | 画像処理装置、画像処理方法及びプログラム | |
CN111246792B (zh) | 解析装置、解析方法、以及程序 | |
JP2020054812A (ja) | 画像処理装置、画像処理方法及びプログラム | |
JP7111874B2 (ja) | 眼科撮影装置 | |
US20190200859A1 (en) | Patterned beam analysis of iridocorneal angle | |
WO2018079765A1 (ja) | 画像処理装置、眼底撮像システム、画像処理方法、およびプログラム | |
JP7384656B2 (ja) | 眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラム | |
EP4289338A1 (en) | Ophthalmic information processing device, ophthalmic device, ophthalmic information processing method, and program | |
JP7387812B2 (ja) | 画像処理装置、画像処理方法及びプログラム | |
JP2024003982A (ja) | 画像処理方法、画像処理装置、プログラム、及び記録媒体 | |
JP2023070264A (ja) | 走査型イメージング装置、その制御方法、プログラム、及び記録媒体 | |
JP2024003558A (ja) | 眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラム | |
CN117425424A (zh) | 对眼睛上的手术进行诊断和规划和/或过程监控的装置和方法 |
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 | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 40026693 Country of ref document: HK |
|
GR01 | Patent grant | ||
GR01 | Patent grant |