CN102812380A - 医用图像处理装置 - Google Patents

医用图像处理装置 Download PDF

Info

Publication number
CN102812380A
CN102812380A CN2012800002792A CN201280000279A CN102812380A CN 102812380 A CN102812380 A CN 102812380A CN 2012800002792 A CN2012800002792 A CN 2012800002792A CN 201280000279 A CN201280000279 A CN 201280000279A CN 102812380 A CN102812380 A CN 102812380A
Authority
CN
China
Prior art keywords
mentioned
bulls
eye chart
image
section
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012800002792A
Other languages
English (en)
Other versions
CN102812380B (zh
Inventor
中野史树
川崎友宽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Medical Systems Corp
Original Assignee
Toshiba Corp
Toshiba Medical Systems Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Toshiba Corp, Toshiba Medical Systems Corp filed Critical Toshiba Corp
Publication of CN102812380A publication Critical patent/CN102812380A/zh
Application granted granted Critical
Publication of CN102812380B publication Critical patent/CN102812380B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/20Drawing from basic elements, e.g. lines or circles
    • G06T11/206Drawing of charts or graphs
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/037Emission tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/503Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5235Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5247Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from an ionising-radiation diagnostic technique and a non-ionising radiation diagnostic technique, e.g. X-ray and ultrasound
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0883Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/483Diagnostic techniques involving the acquisition of a 3D volume of data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A61B8/5261Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from different diagnostic modalities, e.g. ultrasound and X-ray

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Optics & Photonics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Dentistry (AREA)
  • Pulmonology (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • General Physics & Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Nuclear Medicine (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

提高靶心图的诊断能力。本实施方式涉及的医用图像处理装置具备:存储部、提取部、归一化部、生成部、以及显示部。存储部存储与心脏的功能指标相关的三维的功能图像数据。提取部从上述三维的功能图像数据中提取心肌区域。归一化部使上述心肌区域的内壁与外壁之间的距离在既定的数值范围内进行归一化。生成部生成由二维的极坐标来表现与上述既定的数值范围内的规定值对应的上述心肌区域上的位置的像素值的空间分布的靶心图。显示部显示上述靶心图。

Description

医用图像处理装置
技术领域
本实施方式涉及医用图像处理装置。
背景技术
在缺血性心脏病的诊断中,利用了基于SPECT装置的心肌SPETCT。在心肌SPECT中,活用有靶心图。靶心图(bullseye map)的生成方法如以下那样。首先,根据通过SPECT装置收集到的体数据产生多个短轴断面图像。在各短轴断面图像中,设定从心肌区域的中心点放射状地延伸的多条直线。从各直线上的心肌区域的内膜与外膜之间(心肌内部)的多个像素的像素值中确定最大值。对极坐标上的对应位置的像素分配所确定的最大值。由此生成靶心图。这样,在靶心图中,分配有贯穿心肌内部的直线上的最大值。从而,不能通过靶心图把握心肌内部的像素值分布。因此,难以进行关注于心肌的局部部位的缺血区域的危险度的判定。
现有技术文献
专利文献1:日本特开2009-18005号公报
发明内容
实施方式的目的在于提供一种能够提高靶心图的诊断能力的医用图像处理装置。
本实施方式涉及的医用图像处理装置具备:存储部,存储与心脏的功能指标相关的三维的功能图像数据;提取部,从上述三维的功能图像数据中提取心肌区域;归一化部,将上述心肌区域的内壁与外壁之间的距离在既定的数值范围内进行归一化;生成部,生成用二维的极坐标来表现与上述既定的数值范围内的规定值对应的上述心肌区域上的位置的像素值的空间分布的靶心图;显示部,显示上述靶心图。
实现提高靶心图的诊断能力。
附图说明
图1是表示第1实施方式涉及的图像处理装置的结构的图。
图2是表示在图1的控制部的控制下进行的、第1实施方式涉及的心肌SPECT浏览显示处理的典型的流程的图。
图3是用于说明图2的步骤SA1的图。
图4是用于说明图2的步骤SA3的图。
图5是用于说明图2的步骤SA5的图。
图6是表示图2的步骤SA7中的靶心图的显示例的图。
图7是用于说明图2的步骤SA8的图。
图8是表示第2实施方式涉及的医用图像处理装置的结构的图。
图9是表示在图8的控制部的控制下进行的、第2实施方式涉及的心肌SPECT浏览显示处理的典型的流程的图。
图10是用于说明图9的步骤SB3的图。
图11是用于说明通过图8的显示/非显示判定部与血管区域配置部进行的、冠状动脉区域的局部部位的显示/非显示处理的图。
图12是表示第3实施方式涉及的医用图像处理装置的结构的图。
图13是表示在图12的控制部的控制下进行的、第3实施方式涉及的心肌SPECT浏览显示处理的典型的流程的图。
图14是表示图13的步骤SC3中的、配置了表示显示曲断面的空间分布的曲线的CT断面图像的显示例的图。
图15是用于说明通过图12的血管配置而进行的、伴随着指定值的变更的CT断面图像上的曲线的移动的图。
图16是表示基于变形例1涉及的医用图像处理装置的多帧显示处理的流程的示意图。
图17是表示变形例2涉及的医用图像处理装置的结构的图。
图18是表示通过图17的三维模型产生部进行的三维靶心图的产生处理的典型的流程的示意图。
图19是表示通过图17的显示部显示的三维模型图像的显示的应用例的图。
图20是表示变形例3涉及的医用图像处理装置的结构的图。
图21是用于说明基于图20的图像处理部的圆环图像的产生处理的图。
图22是表示通过图20的图像处理部产生的圆环图像的一个例子的图。
图23是表示通过图20的图像处理部产生的三维圆环图像的一个例子的图。
图24是表示通过图20的图像处理部产生的杯子模型图像的一个例子的图。
符号说明
1:医用图像处理装置;11:存储部;13:断面图像产生部;15:心肌区域提取部;17:归一化部;19:指定位置确定部;21:靶心图生成部;23:显示部;25:操作部;27:控制部。
具体实施方式
以下,参照附图说明本实施方式涉及的医用图像处理装置。
(第1实施方式)
图1是表示第1实施方式涉及的医用图像处理装置1的结构的图。如图1所示,医用图像处理装置1具有:存储部11、断面图像产生部13、心肌区域提取部15、归一化部17、指定位置确定部19、靶心图生成部21、显示部23、操作部25、以及控制部27。
存储部11存储与心脏的功能指标相关的三维的功能图像数据。作为三维的功能图像数据,能够应用通过X射线计算机断层摄影装置或磁共振成像装置、超音波诊断装置、SPECT装置、PET装置等医用图像诊断装置产生的功能图像数据。在通过SPECT装置或PET装置产生的体数据中,例如,作为像素值配置有从集聚于心肌的放射性同位素放射出的γ射线的累计值(或者基于累计值的心脏的功能指标的指标值)。另外,在通过X射线计算机断层摄影装置或磁共振成像装置、超音波诊断装置产生的体数据中,配置有与心脏的运动功能相关的功能指标的指标值。例如,作为X射线计算机断层摄影装置或磁共振成像装置中的功能指标,可知有BP(心肌组织内的每单位体积以及单位时间的血流量)、BV(心肌组织内的每单位体积的血流量)、以及MTT(平均通过时间)等。作为超音波诊断装置中的功能指标,例如,可知有心肌的壁运动信息。作为该壁运动信息,例如,可知有与向壁压方向的变化相关的信息(radial-strain)、与向心轴方向的变化相关的信息(longitudinal-strain)、与向圆周方向的变化相关的信息(circumferential-strain)、短轴断面内的与面积重心相关的信息(rotation)、以及作为不同的短轴断面间的旋转的差分的信息(torsion)。以下,假设三维的功能图像是通过心肌SPETCT的技术由SPECT装置产生的体数据(以下,称为心肌SPECT数据)。将心肌SPECT数据的像素的像素值称为心肌SPECT值。另外,存储部11存储用于后述的心肌SPECT浏览显示处理的图像处理程序。
断面图像产生部13根据心肌SPECT数据,产生分别与沿着心脏区域或心脏区域的特定部位的心轴而排列的多个短轴断面对应的多个短轴断面图像。短轴断面意味着与心轴正交的断面。
心肌区域提取部15从心肌SPECT数据中提取心脏区域的特定部位。特定部位可以是心脏区域的任何的解剖学部位。以下,假设特定部位是与临床上有用的左心室的心肌相关的像素区域(以下,称为心肌区域)。典型地,心肌区域提取部15分别从多个短轴断面图像的各个中提取心肌区域。
归一化部17将心肌区域的内壁与外壁之间的距离在既定的数值范围内进行归一化。典型地,归一化部17分别针对多个短轴断面图像,使心肌区域的内壁与外壁之间的距离在既定的数值范围内进行归一化。数值范围能够由用户经由操作部25任意地设定。
指定位置确定部19确定与既定的数值范围内的用户指定的指定值对应的心肌区域上的位置。分别针对多个短轴断面图像的各个,确定与用户指定的指定值对应的心肌区域上的位置。指定值由用户经由操作部25在既定的数值范围内任意地设定。
靶心图生成部21生成通过由左心室区域的心轴周围的旋转角与来自心轴上的基准点的距离而规定的二维的极坐标,来表现与指定值对应的心肌区域上的位置的心肌SPECT值的空间分布的靶心图。可以说指定值与靶心图的显示断面位置相对应。另外,靶心图的显示断面在通过直角坐标系规定的心肌SPECT数据中,形成曲面。从而,将靶心图的显示断面称为显示曲断面。
显示部23将生成的靶心图显示在显示装置上。作为显示装置,例如,能够适当地利用CRT显示器、液晶显示器、有机EL显示器、等离子显示器等。
操作部25按照来自用户的经由输入装置的指示来输入各种指示或信息。例如,操作部25按照来自用户的经由输入装置的指示,指定既定的数值范围内的指定值。作为输入装置,能够适当地利用鼠标或轨迹球等定位设备、模式切换开关等选择设备、或者键盘等输入设备。
控制部27作为医用图像处理装置1的中枢而发挥作用。控制部27通过读出存储于存储部11的图像处理程序,并按照程序控制各部,从而进行心肌SPECT浏览显示处理。
以下,详细说明在控制部27的控制下进行的心肌SPECT浏览显示处理。图2是表示在控制部27的控制下进行的、心肌SPECT浏览显示处理的典型的流程的图。
如图2所示,控制部27以由用户经由操作部25形成心肌SPECT浏览显示处理的开始指示为契机,使断面图像产生部13进行短轴断面图像的产生处理(步骤SA1)。在步骤SA1中,断面图像产生部13根据心肌SPECT数据,产生分别与沿着左心室区域的心轴而排列的多个短轴断面对应的多个短轴断面图像。
图3是用于说明步骤SA1的图。具体而言,首先,断面图像产生部13为了提取左心室区域而在心肌SPECT数据MSD中应用区域扩张法(region growing)。左心室区域是心肌所能具有的心肌SPECT值范围所包含的像素的集合。首先,断面图像产生部13通过用户经由操作部25,或者通过图像处理,在左心室区域设定种子点。接着,断面图像产生部13从种子点中探索位于26附近的体素,并综合满足综合条件的附近体素。综合条件被设定为左心室所能具有的心肌SPECT值范围。由此,提取满足综合条件且与种子点连结的像素区域,即,左心室区域。如果提取左心室区域,则断面图像产生部13按照被提取出的左心室区域的形状来确定左心室区域的心尖部PA与心基部PB。断面图像产生部13将连接心尖部PA与心基部PB的轴设定为心轴(短轴)AC。如果设定心轴AC,则断面图像产生部13设定沿着被设定的心轴AC而排列的多个短轴断面SP。并且,断面图像产生部13根据心肌SPECT数据MSD来产生分别与所设定的多个短轴断面SP对应的多个断面图像IS。各短轴断面图像IS通过直角坐标系来表现出短轴断面SP中的心肌SPECT值的空间分布。各短轴断面图像IS包含有左心室区域。左心室区域包含有心肌区域。在此,假设短轴断面SP的切片编号例如从心尖部PA朝向心基部PB规定为从1到N(整数)。N被规定为短轴断面SP的数量。N能够由用户经由操作部25任意地进行设定。
如果进行步骤SA1,则控制部27使心肌区域提取部15进行心肌区域的提取处理(步骤SA2)。在步骤SA2中,心肌区域提取部15从各短轴断面图像中提取心肌区域。具体而言,心肌区域提取部15在各短轴断面图像中应用区域扩张法(region growing)。心肌区域是心肌所能具有的心肌SPECT值范围所包含的像素的集合。首先,心肌区域提取部15通过用户经由操作部25,或者通过图像处理,在心肌区域设定种子点。接着,心肌区域提取部15从种子点中探索位于26附近的体素,并综合满足综合条件的附近体素。综合条件设定为心肌所能具有的心肌SPECT值范围。由此,提取满足综合条件且与种子点连结的像素区域,即,心肌区域。
如果进行步骤SA2,则控制部27使归一化部17进行归一化处理(步骤SA3)。在步骤SA3中,归一化部17使心肌区域的内壁与外壁之间的间隔在既定的数值范围内进行归一化。
图4是用于说明步骤SA3的图。图4表示各短轴断面图像所包含的心肌区域RM。具体而言,首先,归一化部17从各短轴断面图像中提取心肌区域RM的内壁WI与外壁WO。更具体而言,归一化部17通过描绘心肌区域RM的内侧的边界来提取心肌区域的内侧的闭合曲线(内壁),并通过描绘心肌区域的外侧的边界来提取心肌区域的外侧的闭合曲线(外壁)。内侧的边界是心肌区域RM的两壁中的距离心轴近的边界,外侧的边界是心肌区域RM的两壁中的距离心轴远的边界。如果提取内侧闭合曲线(内壁WI),则归一化部17计算内侧闭合曲线的重心,如果提取外侧闭合曲线,则归一化部17计算外侧闭合曲线(外壁WO)的重心。接着,归一化部17计算内侧闭合曲线的重心与外侧闭合曲线的重心的中心点PC。将所计算出的中心点PC作为各短轴断面图像中的心肌区域RM的中心点进行处理。归一化部17设定从中心点PC放射状地延伸的多个直线LI。在此,为了识别直线LI,在各直线LI上设定编号。归一化部17确定直线LI与内壁WI的交点PI、直线LI与外壁WO的交点PO。归一化部17确定被确定的交点PI以及交点PO的短轴断面图像上的坐标。短轴断面图像上的坐标通过直角坐标系来规定。归一化部17使各直线LI上的交点PI与交点PO之间的距离在既定的数值范围内进行归一化。例如,将既定的数值范围设定为0~1。此时,归一化部17将交点PI的位置设定为0,将交点PO的位置设定为1。交点PI与交点PO之间的直线LI上的各点的位置设定为随着从交点PI向交点PO前进,从交点PI到各点的位置的距离从0到1线性地变化。由此,进行归一化。将通过归一化而规定的各点的位置称为归一化位置。
如果进行步骤SA3,则控制部27待机由用户经由操作部25指定既定的数值范围内的数值(步骤SA4)。用户经由操作部25将数值范围的下限值以上、上限值以下的所希望的数值指定为指定值。当数值范围为从0到1时,将0以上1以下的所希望的数值指定为指定值。指定值也可以经由键盘或触摸屏等输入装置直接地输入数值。另外,也可以利用滑杆以及其调节器来指定指定值。滑杆与数值范围相对应。调节器与指定值相对应。用户经由输入装置,使调节器在滑杆上滑动。控制部27保存有使滑杆上的调节器的位置与指定值对应的列表。控制部27利用该列表来计算与调节器的位置对应的指定值。指定值的指定与靶心图的显示曲断面位置的指定同义。
另外,以往,将直线上的最大心肌SPECT值映射到靶心图上。在这样的以往的处理中,不存在靶心图的显示曲断面这样的概念。
如果由用户指定指定值,则控制部27使指定位置确定部19进行指定位置的确定处理(步骤SA5)。在步骤SA5中,指定位置确定部19确定与在步骤SA4中指定的指定值对应的心肌区域上的位置。具体而言,指定位置确定部19针对各短轴断面图像的各直线,确定与指定值对应的像素的位置(坐标)。所确定的像素的坐标通过直角坐标系来规定。更具体而言,指定位置确定部19确定与指定值对应的像素的坐标和像素所属于的短轴断面的切片编号与像素所属于的直线的编号。另外,指定值也可以预先决定。
如果进行步骤SA5,则控制部27使靶心图生成部21进行靶心图的生成处理(步骤SA6)。在步骤SA6中,靶心图生成部21生成由通过心轴周围的旋转角与来自心轴上的基准点的距离而规定的二维的极坐标,来表现在步骤SA5中确定的位置的心肌SPECT值的空间分布的靶心图。作为基准点,典型的情况是利用心尖部。
图5是用于说明步骤SA6的图。如图5所示,假设在各短轴断面图像IS的各直线LI中,确定与指定值PD对应的像素的坐标。首先,靶心图生成部21确定所确定的坐标的心肌SPECT值。并且,靶心图生成部21通过将所确定的心肌SPECT值分配(映射)到模板上的对应的位置来生成靶心图。模板的中心与左心室区域的心尖部对应。模板的外缘与左心室区域的心基部对应。从模板的中心点到各像素的距离与短轴断面的切片编号(断面位置)对应。模板的中心点周围的角度与心轴周围的角度对应。即,模板的各像素与各直线LI一一对应。模板的各像素根据短轴断面的切片编号与直线LI的编号唯一地确定。在模板的各像素中,分配与该像素对应的直线LI上的像素中的与指定值PD对应的像素的心肌SPECT值。
即,靶心图生成部21确定与在步骤SA5中确定的切片编号与直线LI的编号对应的模板上的像素,对所确定的模板上的像素分配在步骤SA5中确定的像素的心肌SPECT值。通过对模板上的所有的像素进行该处理,从而生成靶心图。
如果进行步骤SA6,则控制部27使显示部23进行显示处理(步骤SA7)。在步骤SA7中,显示部23将在步骤SA6中显示出的靶心图显示在显示装置上。
图6是表示靶心图BM的显示例的图。如图6所示,靶心图BM的各像素由根据心肌SPECT值的颜色来显示。具体而言,显示部23保存有使心肌SPECT值与色值对应的色表。显示部23利用色表来对在步骤SA6中生成的靶心图的每个像素确定色值。并且,显示部23由与所确定的色值对应的颜色来显示靶心图BM的各像素。
另外,在显示画面中,为了对用户提示指定值,配置显示区域RN或显示区域RS即可。在显示区域RN中,显示指定值的数值。在显示区域RS中显示滑杆。在与滑杆上的指定值对应的位置上,重叠有调节器TU。用户通过观察指定值的数值,或者通过观察滑杆上的调节器TU的位置,能够把握指定值。另外,也可以只设置显示区域RN与显示区域RS的任一个。
如果进行步骤SA7,则控制部27待机由用户经由操作部25来变更指定值(步骤SA8)。即,医用图像处理装置1具有能够任意地变更靶心图的显示曲断面的功能。
图7是用于说明步骤SA8的图。如图7所示,例如,指定值PD与输入装置(鼠标)的移动联动地变更即可。具体而言,当在显示画面(或者应用程序窗口)上的特定区域中,将鼠标向右方向拖动时,指定值PD增加与向该右方向的拖动量对应的数值。另外,当在显示画面(或者应用程序窗口)上的特定区域中,将鼠标向左方向拖动时,指定值PD减少与向该右方向的拖动量对应的数值。例如,从基于用户的指定值的变更操作的操作性的观点来看,特定区域设定为靶心图BM的显示区域即可。控制部27按照拖动量与拖动方向,实时计算变更后的指定值。
另外,如上述那样,指定值的变更方法并不只限定于与输入装置的移动连动的方法。例如,也可以通过经由键盘或触摸屏等输入装置直接地输入数值,或移动滑杆上的调节器的位置,来变更指定值。
如果变更指定值(靶心图的显示曲断面),则控制部27为了生成与变更后的指定值相关的靶心图,重复步骤SA5、SA6、以及SA7。具体而言,在步骤SA5中,指定位置确定部19确定与在步骤SA8中设定的变更后的指定值对应的心肌区域上的位置。在步骤SA6中,靶心图生成部21按照所确定的位置的心肌SPECT值,生成与变更后的指定值相关的靶心图。在步骤SA7中,显示部23显示与变更后的指定值相关的靶心图。显示部23随着指定值的变更,实时显示与变更后的指定值相关的靶心图。
这样,控制部27能够任意地变更指定值,即,能够任意地变更靶心图的显示曲断面,并通过实时显示与变更后的指定值对应的靶心图,能够进行靶心图的浏览显示。
不变更指定值,例如,以由用户经由操作部25形成心肌SPECT浏览显示处理的结束指示为契机(步骤SA8:否定),控制部27结束心肌SPECT浏览显示处理。
医师等用户不只关心心肌区域的厚度方向的最大值。例如,当诊断有无缺血区域时,需要沿着心肌区域的厚度方向观察心肌SPECT值。心肌区域的厚度(各直线上的内壁与外壁之间的距离)不相同,根据心肌区域中的部位的不同而不同。从而,如果根据绝对距离(归一化前的距离)来指定心肌区域的厚度方向的位置,则在某个部位中会靠近内壁,或者在其他的部位中会靠近外壁。
通过上述构成,医用图像处理装置1使心肌区域的厚度在既定的数值范围内进行归一化,并按照在该数值范围内的与任意的指定值对应的位置的心肌SPECT值来生成靶心图。这样,通过进行归一化,医用图像处理装置1能够通过单一的指定值,来统括地指定内壁附近或外壁附近、沿着心肌内部等心肌区域的厚度方向的解剖学上同质的位置。即,通过归一化,能够将指定值作为靶心图的显示曲断面的位置进行处理。另外,医用图像处理装置1能够通过按照来自用户的指示任意地变更指定值,从而能够任意地变更靶心图的显示曲断面。由此,医用图像处理装置1能够将沿着心肌区域的厚度方向的任意的位置的心肌SPECT值的空间分布通过靶心图进行浏览显示。通过浏览显示,能够在短时间内调查大范围的部位,例如,实现缺血区域的有无等迅速的诊断或缺血区域的早期发现。
这样,第1实施方式涉及的医用图像处理装置能够实现提高靶心图的诊断能力。
(第2实施方式)
第2实施方式涉及的医用图像处理装置在靶心图中重叠显示血管区域。以下,针对第2实施方式涉及的医用图像处理装置进行说明。另外,在以下的说明中,针对具有与第1实施方式大致相同的功能的构成要素,添加同一符号,重复说明只在必要时进行。
图8是表示第2实施方式涉及的医用图像处理装置2的结构的图。如图8所示,医用图像处理装置2具有:存储部11、断面图像产生部13、心肌区域提取部15、归一化部17、指定位置确定部19、靶心图生成部21、显示部23、操作部25、控制部27、位置对准部29、血管区域提取部31、以及血管区域配置部33。
第2实施方式涉及的存储部11还存储有与心脏的形态指标相关的三维的形态图像数据。作为三维的形态图像数据,能够适当地使用通过X射线计算机断层摄影装置或磁共振成像装置、超音波诊断装置等医用图像诊断装置产生的形态图像数据。以下,假设第2实施方式涉及的三维的形态图像数据是通过由多切片CT,将被注入了造影剂的心脏进行扫描而产生的体数据(以下,称为CT体数据)。在CT体数据中,分配有形态指标之一的CT值。
位置对准部29使心肌SPECT数据与CT体数据进行位置对准。通过位置对准,计算用于使心肌SPECT数据与CT体数据位置对准的位置对准信息。
血管区域提取部31从CT体数据中提取出与血管相关的像素区域(以下,称为血管区域)。血管区域典型的情况是与冠状动脉相关的像素区域(以下,称为冠状动脉区域)。
血管区域配置部33将用与靶心图相同的二维的极坐标系,来表现被提取出的血管区域的空间分布的血管区域配置在靶心图上。具体而言,血管区域配置部33计算通过与靶心图相同的二维的极坐标系而规定的冠状动脉区域的位置。血管区域配置部33在计算出的靶心图上的冠状动脉区域的位置的像素区域绘制点。由此,在靶心图上配置冠状动脉区域。
第2实施方式涉及的显示部23将被配置了冠状动脉区域的靶心图显示在显示装置上。
以下,详细说明在控制部27的控制下进行的第2实施方式涉及的心肌浏览显示处理。图9是表示在控制部27的控制下进行的第2实施方式涉及的心肌浏览显示处理的典型的流程的图。另外,假设在第2实施方式涉及的心肌浏览显示处理的开始时刻,已经利用第1实施方式生成了靶心图。
如图9所示,控制部27首先使位置对准部29进行位置对准处理(步骤SB1)。在步骤SB1中,位置对准部29使心肌SPECT数据与CT体数据位置对准。心肌SPECT数据与CT体数据摄像时的心相位(舒张期或收缩期等)有时不同。从而,心肌SPECT数据所包含的左心室区域与CT体数据所包含的左心室区域典型的情况是大小不同。从而,位置对准部29提取各数据的解剖学特征点,并利用该解剖学特征点执行位置对准。例如,作为解剖学特征点,提取左心室区域的心尖部与心基部。此时,作为位置对准信息,位置对准部29例如计算从CT体数据的心尖部以及心基部向心肌SPECT数据的心尖部以及心基部的坐标转换式。并且,位置对准部29通过将该坐标转换式乘以CT体数据,从而使心肌SPECT数据与CT体数据位置对准。
如果进行步骤SB1,则控制部27使血管区域提取部31进行血管区域提取处理(步骤SB2)。在步骤SB2中,血管区域提取部31从CT体数据中提取冠状动脉区域。具体而言,血管区域提取部31为了提取冠状动脉区域而在CT体数据中应用区域扩张法。首先,血管区域提取部31通过用户经由操作部25,或者通过图像处理,来在各短轴断面图像的冠状动脉区域设定种子点。接着,血管区域提取部31从种子点中探索位于26附近的体素,并综合满足综合条件的附近体素。综合条件设定为冠状动脉所能具有的心肌SPECT值范围。由此,提取满足综合条件且与种子点连结的像素区域,即,提取冠状动脉区域。
如果进行步骤SB2,则控制部27使血管区域配置部33进行配置处理(步骤SB3)。在步骤SB3中,血管区域配置部33在靶心图中配置冠状动脉区域。
图10是用于说明步骤SB3的图。具体而言,首先,血管区域配置部33利用坐标转换式,将CT体数据的冠状动脉区域RCA投影到心肌SPECT数据的各短轴断面图像IS上。并且,血管区域配置部33针对各短轴断面图像IS,计算与靶心图相同的极坐标系中的冠状动脉的位置。并且,血管区域配置部33在靶心图上的被计算出的位置的像素中绘制表示冠状动脉区域的点。该绘制的点的集合是靶心图上的冠状动脉区域RCA’。更具体而言,血管区域配置部33针对各短轴断面图像IS的各直线LI,判定在直线LI上是否存在冠状动脉区域RCA。当判定为在直线LI上没有存在冠状动脉区域RCA时,血管区域配置部33不对与该直线LI对应的靶心图的像素分配表示冠状动脉区域RCA’的点。当判定为在直线LI上存在冠状动脉区域RCA时,血管区域配置部33对与该直线LI对应的靶心图的像素分配表示冠状动脉区域RCA’的点。由此,对靶心图配置冠状动脉区域RCA’。
如果进行步骤SB3,则控制部27使显示部23进行显示处理(步骤SB4)。在步骤SB4中,显示部23在显示装置上显示在步骤SB3中生成的、被配置了冠状动脉区域的靶心图。通过在靶心图上配置冠状动脉区域,从而,用户能够容易地把握靶心图中的解剖学的位置关系。
如果进行步骤SB4,则控制部27结束第2实施方式涉及的心肌SPECT浏览显示处理。
另外,也可以将靶心图上的冠状动脉区域的局部部位的显示/非显示随着指定值的变更进行切换。以下,说明冠状动脉区域的局部部位的显示/非显示处理。
如图8所示,为了进行该冠状动脉区域的局部部位的显示/非显示处理,医用图像处理装置2还具备有显示/非显示判定部35。图11是用于说明通过显示/非显示判定部35与血管区域配置部33进行的、冠状动脉区域的局部部位的显示/非显示处理的图。图11的(a)是表示短轴断面图像(直角坐标系)中的心肌区域RM以及其周边的模式图。虚线L1表示与第1指定值对应的靶心图的显示曲断面。与第1指定值对应的显示曲断面被设置在内壁附近。虚线L2表示与第1指定值对应的靶心图的显示曲断面。与第2指定值对应的显示曲断面被设定在外壁附近。假设在心肌区域RM中,包含有与缺血部位相关的像素区域RI(以下,称为缺血区域)。缺血部位是与正常部位相比血液没有浸润的心肌部位。缺血部位与正常部位相比放射性同位素的集聚量少。从而,缺血区域与和正常部位相关的心肌区域相比心肌SPECT值低。在部位A中,缺血区域RI涵盖从内壁到外壁附近,与部位B相比较,在临床上危险度高。图11的(b)表示与显示曲断面L1相关的靶心图。图11的(b)表示与显示曲断面L2相关的靶心图。
显示/非显示判定部35分别针对与指定值所对应的显示曲断面相关的靶心图所包含的多个像素,按照该像素的心肌SPECT值判定是否将血管区域显示在该像素上。具体而言,显示/非显示判定部35分别针对靶心图的多个像素,判定该像素的心肌SPECT值是否在阈值以上。阈值被设定为能够区别与正常部位相关的心肌区域与缺血区域的心肌SPECT值。
当判定为心肌SPECT值在阈值以下时,血管区域配置部33对靶心图的像素配置表示冠状动脉区域RCA’的点。当判定为心肌SPECT值在阈值以上时,血管区域配置部33不对靶心图的像素配置表示冠状动脉区域RCA’的点。
每当显示/非显示判定部35与血管区域配置部33由用户经由操作部25来变更指定值(显示曲断面),则执行上述的处理。这样,每当通过显示/非显示判定部35与血管区域配置部33,来变更指定值(靶心图的显示曲断面位置),则能够将靶心图上的冠状动脉区域,对每个局部部位进行显示/非显示的切换。即,存在于比阈值低的心肌SPECT值的像素区域附近的冠状动脉区域被配置在靶心图中,存在于比阈值高的心肌SPECT值的像素区域附近的冠状动脉区域被从靶心图中除去。这样,通过将存在于比阈值低的心肌SPECT值的像素区域附近的冠状动脉区域配置在靶心图中,从而用户能够评价各部位的临床上的危险度。
这样,第2实施方式涉及的医用图像处理装置2能够实现提高靶心图的诊断能力。
(第3实施方式)
第3实施方式涉及的医用图像处理装置将表示靶心图的显示曲断面的空间分布的曲线重叠显示于通过直角坐标系规定的形态图像。以下,针对第3实施方式涉及的医用图像处理装置进行说明。另外,在以下的说明中,针对具有与第1以及第2实施方式大致相同的功能的构成要素,添加同一符号,重复说明只在必要时进行。
图12是表示第3实施方式涉及的医用图像处理装置3的结构的图。如图12所示,第3实施方式涉及的医用图像处理装置3具有:存储部11、断面图像产生部13、心肌区域提取部15、归一化部17、指定位置确定部19、靶心图生成部21、显示部23、操作部25、控制部27、位置对准部29、血管区域提取部31、血管区域配置部33、显示/非显示判定部35、曲线配置部37。
第3实施方式涉及的断面图像产生部13根据CT体数据,产生与心肌区域交差的断面(以下,称为CT断面)相关的断面图像(以下,称为CT断面图像)。如果CT断面与心肌区域交差,则可以具有任何位置以及朝向。另外,CT断面的数量可以是一个也可以是多个。
曲线配置部37计算与指定值对应的靶心图的显示曲断面的CT断面图像上的位置。并且,曲线配置部37在所计算出的位置上描绘表示曲线的点。由此,在各CT断面图像中配置表示靶心图的显示曲断面的空间分布的曲线。
图13是表示在第3实施方式涉及的控制部27的控制下进行的心肌SPECT浏览显示处理的典型的流程的图。另外,假设在第3实施方式涉及的心肌浏览显示处理的开始时刻,已经利用第1实施方式生成了靶心图。
如图13所示,控制部27首先使断面图像产生部13进行断面图像的产生处理(步骤SC1)。在步骤SC1中,断面图像产生部13根据CT体数据,产生与正交的3个断面相关的3个CT断面图像。具体而言,首先,断面图像产生部13与步骤SA1相同,在CT体数据中设定心轴。如果设定心轴,则断面图像产生部13产生以心轴为基准的断面图像。在心肌SPECT中,由于临床上的有用性,频繁采用了以心轴为基准的正交的3个断面。即,该3个断面是正交于心轴的短轴断面、与心轴平行相交的垂直长轴断面、以及与心轴平行相交、且与垂直长轴断面正交的水平长轴断面。此时,断面图像产生部13对CT体数据实施MPR(multiplanar reconstruction)处理,产生与短轴断面相关的CT断面图像、与垂直长轴断面相关的CT断面图像、与水平长轴断面相关的CT断面图像。
如果进行步骤SC1,则控制部27使曲线配置部37进行曲线位置的计算处理(步骤SC2)。在步骤SC2中,曲线配置部37对各CT断面图像,配置表示与指定值对应的靶心图的显示曲断面的空间分布的曲线。具体而言,首先,曲线配置部37根据在步骤SB1中计算出的心肌SPECT数据与CT体数据的位置对准信息,计算与指定值对应的靶心图的显示曲断面的各CT断面图像中的位置。更具体而言,首先,曲线配置部37根据指定值,确定构成显示曲断面的心肌SPECT数据上的多个像素。曲线配置部37根据所确定的多个像素,利用位置对准信息,确定构成显示曲断面的各CT断面图像的多个像素。如果确定多个像素,则曲线配置部37对所确定的多个像素,绘制用于绘制曲线的点。由此,在各CT断面图像上配置表示显示曲断面的空间分布的曲线。另外,在点之间不连续,离间的情况下,也可以将点之间用直线或曲线连结。
如果进行步骤SC2,则控制部27使显示部23进行显示处理(步骤SC3)。在步骤SC3中,显示部23显示配置了表示显示曲断面的空间分布的曲线的CT断面图像。
图14是表示配置了表示显示曲断面的空间分布的曲线的CT断面图像IC的显示例。如图14所示,在显示画面上,与靶心图BM一起,配置在步骤SC1中产生的、与水平长轴断面相关的CT断面图像IC1、与垂直长轴断面相关的CT断面图像IC2、以及与短轴断面相关的CT断面图像IC3。在各断面图像IC1、IC2、IC3上,分别重叠有表示在步骤S2中计算出的显示曲断面的空间分布的曲线LC1、LC2、LC3。这样,将表示靶心图的显示曲断面的空间分布的曲线配置于通过直角坐标系规定的形态图像之一的CT断面图像。由此,用户能够容易地把握靶心图的显示曲断面的位置。
当由用户经由操作部25变更了指定值时,曲线配置部37按照变更后的指定值使曲线在CT断面图像上移动。
图15是用于说明伴随着指定值的变更的CT断面图像IC1、IC2、IC3上的曲线LC1、LC2、LC3的移动的图。如图15所示,例如,指定值与输入装置(鼠标)的移动连动地变更即可。血管配置部37按照拖动量与拖动方向实时计算变更后的指定值。并且,血管配置部37对各CT断面图像实时配置表示与变更后的指定值对应的靶心图的显示曲断面的空间分布的曲线。
具体而言,当在显示画面(或者应用程序窗口)上的特定区域中,将鼠标向右方向拖动时,曲线LC1、LC2、LC3在CT断面图像IC1、IC2、IC3上向外壁侧移动与向其右方向的拖动量对应的数值。另外,当在显示画面(或者应用程序窗口)上的特定区域中将鼠标向左方向拖动时,曲线LC1、LC2、LC3在CT断面图像IC1、IC2、IC3上向内壁侧移动与向其左方向的拖动量对应的数值。例如,从基于用户的指定值的变更操作的操作性的观点来看,特定区域设定在靶心图BM的显示区域即可。
另外,如上述那样,指定值的变更方法并不只限定于与输入装置的移动连动的方法。例如,也可以通过经由键盘或触摸屏等输入装置直接地输入数值,或者移动滑杆上的调节器的位置,从而变更指定值。
这样,与经由操作部25的指定值,即,靶心图的显示曲断面的变更连动,能够使CT断面图像上的曲线移动。
这样,第3实施方式涉及的医用图像处理装置能够实现提高靶心图的诊断能力。
接着,针对本实施方式涉及的各种的变形例进行说明。
(变形例1)
变形例1涉及的医用图像处理装置以多帧形式来显示与多个显示曲断面相关的多个靶心图。以下,针对变形例1涉及的医用图像处理装置进行说明。另外,在以下的说明中,针对具有与本实施方式大致相同的功能的构成要素,添加同一符号,重复说明只在必要时进行。
图16是表示基于变形例1涉及的医用图像处理装置的多帧显示处理的流程的示意图。以通过用户经由操作部25形成多帧显示指示为契机,变形例1涉及的控制部27开始多帧显示处理。与本实施方式相同,通过心肌区域提取部15来从心肌SPECT数据中提取心肌区域,并通过归一化部17将心肌区域的内壁与外壁之间的距离在既定的数值范围内进行归一化。
并且,如图16(a)所示,例如,通过指定位置确定部19,在心肌区域MR内,设定厚度方向的位置不同的多个显示曲断面。显示曲断面CP在通过归一化部17归一化后的心肌区域MR的内壁与外壁之间以一定的间隔来设定。如图16(b)所示,靶心图生成部21根据心肌SPECT数据,生成分别与多个显示曲断面CP对应的多个靶心图BM。如图16(c)所示,显示部23以多帧形式在画面上显示分别与多个显示曲断面CP对应的多个靶心图BM。具体而言,多个靶心图BM从内壁侧向内壁侧依次排列并显示于一画面。换而言之,与不同的显示曲断面相关的多个靶心图BM排列显示于一画面。
这样,通过以多帧形式显示多个靶心图,从而,用户能够一览地观察多个靶心图。
这样,变形例1涉及的医用图像处理装置能够实现提高靶心图的诊断能力。
(变形例2)
变形例2涉及的医用图像处理装置显示三维的靶心图。以下,针对变形例2涉及的医用图像处理装置进行说明。另外,在以下的说明中,针对具有与本实施方式大致相同的功能的构成要素,添加同一符号,重复说明只在必要时进行。
图17是表示变形例2涉及的医用图像处理装置4的结构的图。如图17所示,医用图像处理装置4具有:存储部11、断面图像产生部13、心肌区域提取部15、归一化部17、指定位置确定部19、靶心图生成部21、显示部23、操作部25、控制部27、以及三维模型产生部41。
三维模型产生部41根据心肌SPECT数据,产生由变形例2中特有的三维坐标来表现心肌区域内的心肌SPECT值的空间分布的三维模型图像。变形例2中特有的三维坐标通过二维极坐标、归一化后的内壁与外壁之间的高度来规定。二维极坐标被规定为与靶心图相同的极坐标系。高度例如通过靶心图的显示曲断面的归一化数值、或来自内壁或者外壁的靶心图的个数来规定。以下,将通过该三维坐标来表现的图像处理空间称为三维模型空间。
接着,针对通过变形例2涉及的三维模型产生部41进行的三维模型图像的产生处理的流程进行说明。
图18是表示通过三维模型产生部41进行的三维靶心图的产生处理的典型的流程的示意图。如图18(a)所示,假设通过靶心图生成部21,生成厚度方向的位置不同的多个靶心图BM。三维模型产生部41通过阈值处理等,从多个靶心图BM中提取可视化对象区域。可视化对象区域可以是临床上关心的任何区域,但为了具体地进行以下的说明,假设为缺血区域。此时,三维模型产生部41通过阈值处理等从多个靶心图BM中提取缺血区域IR。
典型地,缺血区域IR从内壁朝朝向外壁而形成。从而,在三维模型数据中,缺血区域IR的高度表现Transmural Extent(与壁厚方向的缺血区域的宽度相关的指标)。Transmural Extent在临床上,被用于心肌梗塞的严重度的指标。
如图18(b)所示,三维模型产生部41确定缺血区域IR所包含的各像素的极坐标与高度坐标,对与所确定的极坐标与高度坐标对应的三维模型空间内的体素分配体素值。分配的体素可以是心肌SPECT值,也可以是表示为缺血区域IR的既定值(例如,“1”等)。这样,通过对三维模型空间分配缺血区域IR内的所有的像素的像素值,从而产生三维模型数据MD。Transmural Extent在三维模型数据MD中,由从内壁向外壁延伸的柱状图来表示。
另外,如图18(b)所示,为了容易地把握位置关系,模式性的外壁MAO与模式性的内壁MAI进行位置匹配,并配置在三维模型数据MD中即可。模式性的外壁MAO与模式性的内壁MAI也可以分别是与外壁的显示曲断面相关的靶心图BMO和与内壁的显示曲断面相关的靶心图BMI。另外,也可以不配置模式性的外壁MAO与模式性的内壁MAI。
三维模型产生部41对三维模型数据实施现有的绘制处理,产生显示用三维模型图像。绘制中的视点以及视线能够设定为任意的位置以及方向。然而,为了使医师等用户容易地把握位置关系,将视点以及视线初期设定成隔壁部等配置在画面的附近即可。所产生的三维模型图像通过显示部23来显示。
用户通过观察根据三维模型数据的三维模型图像,能够在图像上容易地把握Transmural Extent。
另外,在上述的说明中,假设三维模型图像根据心肌SPECT数据经由多个靶心图产生。然而,变形例2涉及的三维模型图像也可以不经由多个靶心图,一起根据心肌SPECT数据直接地产生三维模型图像。
接着,参照图19针对三维模型图像的显示的应用例进行说明。如图19所示,显示部23将目前显示出的靶心图的位置重叠显示于三维模型图像MI。具体而言,如变形例1那样,显示部23以多帧形式显示厚度方向的位置不同的多个靶心图。用户经由操作部25选择想要观察的靶心图。所选择出的靶心图BM以原始的尺寸通过显示部23来显示。另外,所选择出的靶心图BM’的位置通过显示部23显示在三维模型图像MI上。具体而言,显示部23确定所选择出的靶心图BM’的高度坐标,并将靶心图BM’的标记MA1重叠显示于三维模型图像MI上的所确定的高度坐标。
为了容易地把握三维模型图像MI与靶心图BM的位置关系,在三维模型图像MI中显示标记MA2,在靶心图BM中显示标记MA3即可。标记MA2表示三维模型图像MI的基准位置。典型地,基准位置是三维模型图像MI的正面。标记MA3表示与标记MA2对应的位置。例如,当标记MA2表示三维模型图像MI的正面时,标记MA3表示与三维模型图像MI的正面对应的靶心图BM上的位置。另外,基准位置并不限定于三维模型图像MI的正面,在靶心图的观察中,也可以表示一般所利用的解剖学基准点。
另外,为了提高三维模型图像的识别性,也可以限定于关心区域来显示缺血区域。例如,由用户经由操作部25将关心区域设定在三维模型图像MI的模式性的外壁或者模式性的内壁上即可。显示部23将不属于关心区域内的极坐标的缺血区域从三维模型图像中除去,只将属于关心区域内的极坐标的缺血区域描绘在三维模型图像中。由此,能够由用户在三维模型图像中只描绘出关心的某个缺血区域。
另外,与第2实施方式相同,也可以将血管区域重叠于三维模型图像的模式性的外壁或者模式性的内壁。例如,显示部23显示Transmural Extent值比阈值大的区域的血管区域(例如,冠状动脉的血管区域)即可。另外,血管区域的显示的开与关能够根据来自操作部25的指示进行切换即可。由此,用户例如能够容易地识别冠状动脉与梗塞心肌的关联。
这样,变形例2涉及的医用图像处理装置能够实现提高靶心图的诊断能力。
(变形例三维)
变形例3涉及的医用图像处理装置根据与靶心图不同的显示样态,显示使心壁厚方向的信息可视化的圆环图像。以下,针对变形例3涉及的医用图像处理装置进行说明。另外,在以下的说明中,针对具有与本实施方式大致相同的功能的构成要素,添加同一符号,重复说明只在必要时进行。
图20是表示变形例3涉及的医用图像处理装置5的结构的图。如图20所示,医用图像处理装置4具有存储部11、断面图像产生部13、心肌区域提取部15、归一化部17、指定位置确定部19、靶心图生成部21、显示部23、操作部25、控制部27、以及图像处理部51、。
图像处理部51根据心肌SPECT数据,产生由模式性的圆环形状来表现心肌区域的像素值分布的图像(圆环图像)。圆环图像通过显示部23来显示。
图21是用于说明基于图像处理部51的圆环图像的产生处理的图。如图21所示,图像处理部51将与心轴AC正交的多个断面CSn设定为心肌SPECT数据。在图21中,为了列举一个例子,示出了5个断面CS1、CS2、CS3、CS4、以及CS5。在各断面CSn中,心肌区域MR具有大致圆环形状。图像处理部51将各断面CSn中的心肌区域MR的形状,基准化(标度)为预先设定的圆环形状。换而言之,图像处理部51对圆环模型上的对应的像素分配各断面CSn中的心肌区域MR的各像素的像素值。圆环模型随着距离心尖部的距离变远,而具有大的半径。通过使各圆环模型中心一致地配置,产生圆环图像。
图22是表示通过图像处理部51产生的圆环图像RI的一个例子的图。如图22所示,圆环图像RI包含有分别与多个断面对应的多个圆环区域RRN。在图22中,如与图21对应的那样,示出5个圆环区域RR1、RR2、RR3、RR4、以及RR5。多个圆环区域RRn是共有中心CP的同心圆环。各圆环区域RRn以不空间性地重叠的方式进行基准化。各圆环区域RRn通过二维直角坐标系来规定。例如,各圆环区域RRn的内缘与心肌区域的内壁对应,各圆环区域RRn的外缘与心肌区域的外壁对应。
这样,圆环图像能够由基准化后的直角坐标系来表现沿着心轴排列的多个断面上的像素值的空间分布。从而,能够容易地把握心壁厚方向的心肌SPECT值的分布。另外,由于在圆环图像中,心肌区域的形状以标准的圆环模型进行基准化,因此,不被个体差或检查状况等影响,能够容易地把握解剖学的位置关系。
另外,图像处理部51也能够根据心肌SPECT数据,产生由模式性的圆环形状来三维地表现心肌区域的像素值分布的图像(三维圆环图像)。三维圆环图像通过显示部23来显示。
图23是表示通过图像处理部51产生的三维圆环图像RI’的一个例子的图。在三维圆环图像RI’中,多个圆环区域RRn按照与心轴AC方向相关的位置,沿着心轴AC进行排列。三维圆环图像RI’的视点以及视线按照经由基于用户的操作部25的指示任意地设定即可。例如,如图23所示,视点以及视线设定为能够从斜上方俯瞰多个圆环区域RRn即可。各圆环区域RRn也可以通过与心轴方向相关的带有厚度的MIP来产生。由此,显示部23能够同时显示与心轴方向相关的心肌SPECT值的空间分布。用户通过经由操作部25使视点在心轴AC周围旋转,能够观察全方位的心肌SPECT值的分布。
另外,图像处理部51能够根据心肌SPECT数据产生杯子模型图像。图24是表示通过图像处理部51产生的杯子模型图像CMI的一个例子的图。如图24所示,杯子模型图像CMI通过对具有模式性的心肌形状的杯子形状的三维模型分配内壁以及外壁的心肌SPECT值而产生。杯子模型图像CMI的视点以及视线,按照经由基于用户的操作部25的指示任意地设定即可。由此,通过经由操作部25使视点在心轴AC周围旋转,用户能够观察全方位的心肌SPECT值的分布。
虽然说明了本发明的几个实施方式,但这些实施方式是作为例子而提示的,并不意图限定本发明的范围。这些实施方式能够以其他的各种形态进行实施,在不脱离发明的要旨的范围内,能够进行各种的省略、置换、变更。这些实施方式或其变形与包含于发明的范围或要旨中一样,包含于权利要求书记载的发明及其等同的范围中。

Claims (9)

1.一种医用图像处理装置,其特征在于,具备:
存储部,存储与心脏的功能指标相关的三维的功能图像数据;
提取部,从上述三维的功能图像数据中提取心肌区域;
归一化部,使上述心肌区域的内壁与外壁之间的距离在既定的数值范围内进行归一化;
生成部,用二维的极坐标来表现与上述既定的数值范围内的规定值对应的上述心肌区域上的位置的像素值的空间分布;以及
显示部,显示上述靶心图。
2.根据权利要求1所述的医用图像处理装置,其特征在于,
还具备确定部,确定与上述既定的数值范围内的上述规定值对应的上述心肌区域上的位置,
上述生成部按照上述确定的位置的像素值生成上述靶心图。
3.根据权利要求2所述的医用图像处理装置,其特征在于,
还具备产生部,根据上述三维的功能图像数据,产生沿着连接心尖部与心基部的心轴而排列的多个断面图像,
上述提取部从上述多个断面图像中分别提取多个心肌区域,
上述归一化部使上述多个心肌区域的各个的内壁与外壁之间的距离在既定的数值范围内进行归一化,
上述确定部确定与上述既定的数值范围内的用户指定的指定值对应的上述多个心肌区域的各自的位置。
4.根据权利要求1所述的医用图像处理装置,其特征在于,
还具备指定部,按照来自用户的指示在上述既定的数值范围内指定上述规定值。
5.根据权利要求1所述的医用图像处理装置,其特征在于,
还具备:位置对准部,按照解剖学特征点,将与形态指标相关的三维的形态图像数据与上述三维的功能图像数据进行位置对准;
提取部,从上述三维的形态图像数据中提取血管区域;以及
配置部,将用二维的极坐标系来表现上述提取出的血管区域的空间分布的血管区域配置在上述靶心图上,其中
上述显示部显示配置了上述血管区域的靶心图。
6.根据权利要5所述的医用图像处理装置,其特征在于,
还具备判定部,分别针对上述靶心图所包含的多个像素的各个,按照上述像素的像素值来判定是否将上述血管区域显示在上述像素上,
当判定为显示上述血管区域时,上述配置部将上述血管区域配置于上述靶心图上的上述像素,当判定为没有显示上述血管区域时,不将上述血管区域配置于上述靶心图上的上述像素。
7.根据权利要求1所述的医用图像处理装置,其特征在于,
还具备产生部,产生根据与形态指标相关的三维的形态图像数据的断面图像;以及
配置部,在上述断面图像上,配置表示与上述指定值对应的上述靶心图的显示曲断面的空间分布的曲线,其中,
上述显示部显示被配置了上述曲线的靶心图。
8.根据权利要求1所述的医用图像处理装置,其特征在于,
上述生成部生成分别与上述既定的数值范围内的多个规定值对应的多个靶心图,
上述显示部排列显示上述多个靶心图。
9.根据权利要求1所述的医用图像处理装置,其特征在于,
还具备三维模型图像产生部,产生通过用二维极坐标与从归一化的内壁到外壁的高度来规定的三维坐标,来表现上述心肌区域内的像素值的空间分布的三维模型图像,
上述显示部显示上述三维模型图像。
CN201280000279.2A 2011-03-09 2012-03-01 医用图像处理装置 Active CN102812380B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2011051396 2011-03-09
JP2011-051396 2011-03-09
PCT/JP2012/055240 WO2012121105A1 (ja) 2011-03-09 2012-03-01 医用画像処理装置

Publications (2)

Publication Number Publication Date
CN102812380A true CN102812380A (zh) 2012-12-05
CN102812380B CN102812380B (zh) 2016-01-20

Family

ID=46798072

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201280000279.2A Active CN102812380B (zh) 2011-03-09 2012-03-01 医用图像处理装置

Country Status (4)

Country Link
US (1) US8942460B2 (zh)
JP (1) JP6129476B2 (zh)
CN (1) CN102812380B (zh)
WO (1) WO2012121105A1 (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106102583A (zh) * 2014-01-10 2016-11-09 泰勒顿国际公司 瘢痕和纤维心脏区的检测
CN107249427A (zh) * 2015-07-06 2017-10-13 奥林巴斯株式会社 医疗装置、医疗图像生成方法以及医疗图像生成程序
CN109859165A (zh) * 2018-12-24 2019-06-07 新绎健康科技有限公司 一种取脉点的定位方法及装置
CN111816313A (zh) * 2020-07-01 2020-10-23 杭州脉流科技有限公司 基于左心室三维模型分析心肌运动的模拟显示方法、装置、计算机设备和存储介质
CN111816282A (zh) * 2020-07-01 2020-10-23 杭州脉流科技有限公司 基于ct影像的左心室功能学分析方法、装置、计算机设备和存储介质
CN113643285A (zh) * 2021-09-16 2021-11-12 人工智能与数字经济广东省实验室(广州) 基于人工神经网络的心肌靶心图分析方法、***和介质

Families Citing this family (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012121105A1 (ja) * 2011-03-09 2012-09-13 株式会社 東芝 医用画像処理装置
WO2014084398A1 (ja) * 2012-11-30 2014-06-05 株式会社 東芝 医用画像診断装置
WO2014192751A1 (ja) * 2013-05-27 2014-12-04 株式会社東芝 画像処理装置及び画像処理方法
JP5669977B1 (ja) * 2013-08-29 2015-02-18 日本メジフィジックス株式会社 心筋輪郭抽出技術
JP5968928B2 (ja) * 2014-02-14 2016-08-10 富士フイルム株式会社 医用画像表示制御装置およびその作動方法並びに医用画像表示制御プログラム
JP2015150325A (ja) * 2014-02-18 2015-08-24 株式会社東芝 医用画像処理装置及び医用画像処理方法
FR3020700B1 (fr) * 2014-04-30 2016-05-13 Univ Bordeaux Methode pour la quantification de la presence de graisses dans une region du cœur
JP6285304B2 (ja) * 2014-07-15 2018-02-28 富士フイルムRiファーマ株式会社 コンピュータプログラム、画像処理装置及び方法
JP6452974B2 (ja) * 2014-07-15 2019-01-16 キヤノンメディカルシステムズ株式会社 医用画像処理装置
JP6357049B2 (ja) * 2014-08-11 2018-07-11 キヤノン株式会社 画像処理装置、医用画像診断システム、画像処理方法およびプログラム
JP6098677B2 (ja) * 2014-10-10 2017-03-22 キヤノンマーケティングジャパン株式会社 医用画像処理装置、医用画像処理装置に搭載可能なプログラム、及び医用画像処理方法
JP6734028B2 (ja) * 2014-11-14 2020-08-05 キヤノンメディカルシステムズ株式会社 医用画像診断装置、画像処理装置及び画像生成方法
EP3250128B1 (en) * 2015-01-26 2020-10-21 Spectrum Dynamics Medical Limited Systems and methods for medical image registration
KR101690655B1 (ko) * 2015-01-29 2016-12-28 삼성전자주식회사 의료 영상 처리 장치 및 의료 영상 처리 방법
KR101690654B1 (ko) 2015-01-29 2016-12-28 삼성전자주식회사 의료 영상 장치 및 의료 영상을 처리하는 방법
CN105279796A (zh) * 2015-10-10 2016-01-27 苏州工业园区凌志软件股份有限公司 一种3d柱状图的绘制方法及装置
US10346993B2 (en) * 2015-11-17 2019-07-09 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for image processing in magnetic resonance imaging
JP6639935B2 (ja) * 2016-02-09 2020-02-05 キヤノンメディカルシステムズ株式会社 画像処理装置およびmri装置
KR102380216B1 (ko) 2016-11-29 2022-03-28 블랙모어 센서스 앤드 애널리틱스, 엘엘씨 포인트 클라우드 데이터 세트에서 객체의 분류를 위한 방법 및 시스템
EP3548841A4 (en) 2016-11-30 2020-06-10 Blackmore Sensors And Analytics Inc. METHOD AND SYSTEM FOR DOPPLER DETECTION AND DOPPLER CORRECTION OF THE OPTICAL DETECTION OF THE CHIRPED AREA
JP6811862B2 (ja) 2016-11-30 2021-01-13 ブラックモア センサーズ アンド アナリティクス エルエルシー 光学測距システムによる適応走査の方法およびシステム
US10422880B2 (en) 2017-02-03 2019-09-24 Blackmore Sensors and Analytics Inc. Method and system for doppler detection and doppler correction of optical phase-encoded range detection
WO2019209727A1 (en) 2018-04-23 2019-10-31 Blackmore Sensors and Analytics Inc. Method and system for controlling autonomous vehicle using coherent range doppler optical sensors
US10813612B2 (en) 2019-01-25 2020-10-27 Cleerly, Inc. Systems and method of characterizing high risk plaques
US11308610B2 (en) * 2019-12-11 2022-04-19 Shanghai United Imaging Intelligence Co., Ltd. Systems and methods for machine learning based automatic bullseye plot generation
KR20220124217A (ko) 2020-01-07 2022-09-13 클리어리, 인크. 의료 영상 분석, 진단, 위험도 층화, 의사 결정 및/또는 질환 추적을 위한 시스템, 방법 및 디바이스
US20220392065A1 (en) 2020-01-07 2022-12-08 Cleerly, Inc. Systems, methods, and devices for medical image analysis, diagnosis, risk stratification, decision making and/or disease tracking
US11969280B2 (en) 2020-01-07 2024-04-30 Cleerly, Inc. Systems, methods, and devices for medical image analysis, diagnosis, risk stratification, decision making and/or disease tracking
US20230289963A1 (en) 2022-03-10 2023-09-14 Cleerly, Inc. Systems, devices, and methods for non-invasive image-based plaque analysis and risk determination

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006075216A (ja) * 2004-09-07 2006-03-23 Toshiba Corp 医用画像処理システム、医用画像処理プログラムおよび医用画像処理方法
CN1914617A (zh) * 2004-02-03 2007-02-14 美国西门子医疗解决公司 与心脏有关的疾病和病症的自动诊断和决策支持***和方法
JP2008173236A (ja) * 2007-01-17 2008-07-31 Toshiba Corp 医用画像表示装置
US20090005672A1 (en) * 2007-06-28 2009-01-01 Satoshi Sugiura Image processing apparatus, image diagnostic apparatus and image processing method
CN101601071A (zh) * 2006-07-26 2009-12-09 皇家飞利浦电子股份有限公司 应力水平心脏功能分析结果的可视化
EP2130490A1 (en) * 2007-03-14 2009-12-09 Fujifilm Corporation Heart function display device and program therefor
CN101796544A (zh) * 2007-09-03 2010-08-04 皇家飞利浦电子股份有限公司 体素数据的可视化

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102005002949A1 (de) * 2005-01-21 2006-08-03 Siemens Ag Verfahren zur Visualisierung von Schädigungen im Myokard
JP5405046B2 (ja) * 2007-06-28 2014-02-05 株式会社東芝 画像処理装置、画像診断装置および画像処理方法
JP2009018005A (ja) 2007-07-11 2009-01-29 Toshiba Corp ブルズアイマップ作成装置、ブルズアイマップ作成方法、及び画像処理装置
WO2010098444A1 (ja) * 2009-02-26 2010-09-02 株式会社 日立メディコ 医用画像処理装置及び方法
JP4926199B2 (ja) * 2009-04-16 2012-05-09 富士フイルム株式会社 診断支援装置、診断支援プログラムおよび診断支援方法
JP5502583B2 (ja) * 2009-11-12 2014-05-28 学校法人日本医科大学 画像表示装置および画像表示プログラム並びに画像診断装置
WO2012121105A1 (ja) * 2011-03-09 2012-09-13 株式会社 東芝 医用画像処理装置
WO2012174495A2 (en) * 2011-06-17 2012-12-20 Carnegie Mellon University Physics based image processing and evaluation process of perfusion images from radiology imaging

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1914617A (zh) * 2004-02-03 2007-02-14 美国西门子医疗解决公司 与心脏有关的疾病和病症的自动诊断和决策支持***和方法
JP2006075216A (ja) * 2004-09-07 2006-03-23 Toshiba Corp 医用画像処理システム、医用画像処理プログラムおよび医用画像処理方法
CN101601071A (zh) * 2006-07-26 2009-12-09 皇家飞利浦电子股份有限公司 应力水平心脏功能分析结果的可视化
JP2008173236A (ja) * 2007-01-17 2008-07-31 Toshiba Corp 医用画像表示装置
EP2130490A1 (en) * 2007-03-14 2009-12-09 Fujifilm Corporation Heart function display device and program therefor
US20090005672A1 (en) * 2007-06-28 2009-01-01 Satoshi Sugiura Image processing apparatus, image diagnostic apparatus and image processing method
CN101796544A (zh) * 2007-09-03 2010-08-04 皇家飞利浦电子股份有限公司 体素数据的可视化

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
赵书俊等: "心肌灌注SPECT靶心图定量分析方法的现状和发展", 《国外医学 放射医学核医学分册》 *
赵书俊等: "断层核素心肌显像定量分析方法的研究", 《中原工学院学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106102583A (zh) * 2014-01-10 2016-11-09 泰勒顿国际公司 瘢痕和纤维心脏区的检测
CN107249427A (zh) * 2015-07-06 2017-10-13 奥林巴斯株式会社 医疗装置、医疗图像生成方法以及医疗图像生成程序
CN107249427B (zh) * 2015-07-06 2019-05-07 奥林巴斯株式会社 医疗装置、医疗图像生成方法以及医疗图像生成程序
CN109859165A (zh) * 2018-12-24 2019-06-07 新绎健康科技有限公司 一种取脉点的定位方法及装置
CN109859165B (zh) * 2018-12-24 2023-06-09 新绎健康科技有限公司 一种取脉点的定位方法及装置
CN111816313A (zh) * 2020-07-01 2020-10-23 杭州脉流科技有限公司 基于左心室三维模型分析心肌运动的模拟显示方法、装置、计算机设备和存储介质
CN111816282A (zh) * 2020-07-01 2020-10-23 杭州脉流科技有限公司 基于ct影像的左心室功能学分析方法、装置、计算机设备和存储介质
CN113643285A (zh) * 2021-09-16 2021-11-12 人工智能与数字经济广东省实验室(广州) 基于人工神经网络的心肌靶心图分析方法、***和介质

Also Published As

Publication number Publication date
JP6129476B2 (ja) 2017-05-17
US20120263368A1 (en) 2012-10-18
JP2012196436A (ja) 2012-10-18
US8942460B2 (en) 2015-01-27
CN102812380B (zh) 2016-01-20
WO2012121105A1 (ja) 2012-09-13

Similar Documents

Publication Publication Date Title
CN102812380B (zh) 医用图像处理装置
US10937544B2 (en) Medical imaging apparatus for displaying x-ray images of different types
CN102301393B (zh) 将生物结构中的时变参数可视化
CN101336844B (zh) 医用图像处理装置以及医用图像诊断装置
CN102208114B (zh) 投影图像生成技术
JP2024063008A (ja) 血管内の特徴を検出し且つ表示するためのシステム、方法及びコンピュータ可読媒体
JP5309109B2 (ja) 医用画像処理装置および方法、並びにプログラム
US10524765B2 (en) Refinement of an anatomical model using ultrasound
JP6200249B2 (ja) 情報処理装置、情報処理方法
JP5901930B2 (ja) 画像処理装置
CN106456125A (zh) 用于将医学图像中的特征链接到解剖结构模型的***以及其操作方法
CN101809628A (zh) 解剖数据的可视化
CN103491877B (zh) 医用图像显示装置、医用图像显示方法
CN101961245A (zh) 在采集的放射线照相术图像中补偿呼吸运动的***和方法
CN103732151A (zh) 医用图像装置以及医用图像构成方法
US20140055448A1 (en) 3D Image Navigation Method
CN106097422A (zh) 基于大数据的肝脏三维图像动态演示***
JP2001087228A (ja) 読影支援装置
JP5159195B2 (ja) 医用画像処理装置
JP6387166B2 (ja) 情報処理装置、情報処理方法
Pieper Visualization and Display for Image-Guided Therapy

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
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20160713

Address after: Japan Tochigi

Patentee after: Toshiba Medical System Co., Ltd.

Address before: Tokyo, Japan, Japan

Patentee before: Toshiba Corp

Patentee before: Toshiba Medical System Co., Ltd.