CN103565455A - 医用图像处理装置、医用图像诊断装置以及医用图像处理方法 - Google Patents

医用图像处理装置、医用图像诊断装置以及医用图像处理方法 Download PDF

Info

Publication number
CN103565455A
CN103565455A CN201310312557.4A CN201310312557A CN103565455A CN 103565455 A CN103565455 A CN 103565455A CN 201310312557 A CN201310312557 A CN 201310312557A CN 103565455 A CN103565455 A CN 103565455A
Authority
CN
China
Prior art keywords
bifurcation point
sections
profile
vascular pattern
downside
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
CN201310312557.4A
Other languages
English (en)
Other versions
CN103565455B (zh
Inventor
A·格雷厄姆
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 CN103565455A publication Critical patent/CN103565455A/zh
Application granted granted Critical
Publication of CN103565455B publication Critical patent/CN103565455B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/46—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient
    • A61B6/461—Displaying means of special interest
    • A61B6/466—Displaying means of special interest adapted to display 3D data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50—Clinical applications
    • A61B6/504—Clinical applications involving diagnosis of blood vessels, e.g. by angiography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computerised tomographs
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computerised tomographs
    • A61B6/037—Emission tomography

Abstract

提供一种对于在分岔点或者分岔点附近自动地决定血管树构造的血管的轮廓有益的装置等。具备血管模型生成单元，关于至少包括分岔点、相对该分岔点位于血流的上游侧的上侧节段、相对该分岔点位于血流的下游侧的第1下侧节段、相对该分岔点位于血流的下游侧且不同于所述第1下侧节段的第2下侧节段在内的被检体的血管，至少生成第1血管模型和第2血管模型，该第1血管模型包括所述上侧节段、所述分岔点以所述第1下侧节段，该第2血管模型包括所述上侧节段、所述分岔点以及所述第2下侧节段；以及轮廓分岔点确定单元，至少使用所述第1血管模型和所述第2血管模型，确定所述血管的轮廓分岔点。

Description

医用图像处理装置、医用图像诊断装置以及医用图像处理方法
技术领域
本说明书记载的实施方式，总的说来，例如涉及医用图像用途中的分岔构造的模型化的方法以及装置。实施方式特别地应用于特别地在分岔点或者分岔点附近自动地决定血管树构造的血管的轮廓。
背景技术
众所周知利用三维医用图像技术来获得表现患者的身体或其他的对象的内部的构造的图像数据。
例如，周知使用计算机断层摄影(CT)扫描仪来获得表现血管或者其他的构造的三维图像数据。通常，造影剂被注入到被检者的血管，使得通过关心血管。如果存在造影剂，则使用CT扫描而获得的表现表现血管的图像数据具有比较高的亮度，易于从周围的组织及其他的构造体中区别出来。
另外，例如，能够使用包括MRI、PET、或者X射线血管造影在内的其它的各种三维图像技术对血管进行图像化。
已知非常多的为了对血管以及其它的分岔构造进行识别并模型化，而自动地处理三维图像数据的技术。在众多的既知技术中，利用沿着该构造的一串点，例如解析该构造的强度分布来决定血管或者其它的分岔构造的中心线的位置。对这些点彼此进行接合或者拟合，能够生成沿着血管或者其他的构造体的长度延伸的中心线。
这对决定表现血管或其他的构造体的边界的轮廓也是有用的。例如，在血管的情况下，关于轮廓，因为根据表现轮廓内的点的图像数据的亮度等级以及亮度变化自动地计算血管内的狭窄的程度，所以能够在后续的步骤中使用。因此，迅速且正确地决定轮廓是很重要的。
在一部分既知技术中，将针对沿着血管的一串点的轮廓生成在大致垂直于这些点处的血管的中心线的平面上。关于沿着血管的一串点的各个，决定轮廓面，分析该轮廓面上的数据的强度分布，决定血管的边界。沿着血管的该点处的边界接下来用轮廓表现。
但是，中心线也许是连续的，但中心线的弯曲一般不是平滑的(曲线的一次导函数不是连续的)。例如能够通过利用最佳适应算法而将中心线的曲线设为平滑，在此情况下，有修正后的中心线通过钙、即根据场所而处于血管的外侧的担心。关于这一点，参照图1a、图1b以及图1c进行说明。
图1a示出血管4的通过计算求出的中心线2。可知，中心线是连续的，但沿着长度方向急剧地变化。血管包括钙沉着或者其它的闭塞部6。可知，(例如如图1b所示那样的)如果在一串的点，将轮廓面7a、7b、7c、7d、7e、7f、7g设为相对中心线垂直，则其结果得到的轮廓过于急剧地改变方向。根据情况，在方向变化中的一个的紧前面和紧后面决定的轮廓有时重叠，使得将通过轮廓表现的血管返回原形时看起来是错误的。
图1c示出图1a的血管4，但是此处的中心线处于平滑状态。可知如果使中心线8平滑，则也许能够得到一连串的平滑的轮廓，但却由于设为平滑而导致中心线8不再沿着血管4的中心而是通过钙沉着区域的结果。
操作者观察重叠有计算而得到的中心线和轮廓的血管或者其他的构造的图像，能够手动地修正中心线或者轮廓的位置。但是，这对于操作者来说需要花费时间，存在效率差的问题。而且，在操作者在计算而得到的中心线或者轮廓中发现重大的错误的情况下，操作者对计算出中心线或者轮廓的图像处理系统的信赖度降低。
多数情况下，体内的血管为图2中示意地示出的树状构造的方式。图2的血管构造具备根节点10，根节点10经由分岔点20、22、24到达各种各样的叶节点12、14、16、18。所谓血管，能够视为在从根节点10至叶节点12的路径中可通过一个或者多个分岔点的路径、或者这样的路径的一部分。所谓节段能够视为从一个分岔点至下一分岔点的(例如，分岔点20至分岔点22的)、或者根10至下一分岔点20的构造的一部分。
在一部分既知的系统中，有时关于不同的血管独立地计算轮廓，这意味着关于形成不同的血管的一部分的构造的部分多次地再次计算轮廓。例如，分岔点20与22之间的节段形成根节点10与叶节点12之间的血管的一部分，而且还形成根节点与叶节点14之间的血管的一部分。如果在操作者必须针对某个血管的节段或者节段的一部分手动调整轮廓的情况下，对于该相同节段或者节段的一部分，在考虑另一血管而已经手动调整过轮廓时，对于操作者来说，将变得特别浪费时间。
发明内容
实施方式提供一种医用画像处理装置，其特征在于，具备：血管模型生成单元，关于至少包括分岔点、相对该分岔点位于血流的上游侧的上侧节段、相对该分岔点位于血流的下游侧的第1下侧节段、相对该分岔点位于血流的下游侧且不同于所述第1下侧节段的第2下侧节段在内的被检体的血管，至少生成第1血管模型和第2血管模型，该第1血管模型包括所述上侧节段、所述分岔点、所述第1下侧节段，该第2血管模型包括所述上侧节段、所述分岔点、所述第2下侧节段；以及轮廓分岔点确定单元，至少使用所述第1血管模型和所述第2血管模型来确定所述血管的轮廓分岔点。
另外，实施方式还提供一种医用画像诊断装置，其特征在于，具备：血管模型生成单元，关于至少包括分岔点、相对该分岔点位于血流的上游侧的上侧节段、相对该分岔点位于血流的下游侧的第1下侧节段、相对该分岔点位于血流的下游侧且不同于所述第1下侧节段的第2下侧节段在内的被检体的血管，至少生成第1血管模型和第2血管模型，该第1血管模型包括所述上侧节段、所述分岔点以及所述第1下侧节段，该第2血管模型包括所述上侧节段、所述分岔点以及所述第2下侧节段；以及轮廓分岔点确定单元，至少使用所述第1血管模型和所述第2血管模型，确定所述血管的轮廓分岔点。
进而，实施方式提供一种医用图像处理方法，其特征在于包括：血管模型生成步骤，关于至少包括分岔点、相对该分岔点位于血流的上游侧的上侧节段、相对该分岔点位于血流的下游侧的第1下侧节段、相对该分岔点位于血流的下游侧且不同于所述第1下侧节段的第2下侧节段在内的被检体的血管，至少生成第1血管模型和第2血管模型，该第1血管模型包括所述上侧节段、所述分岔点以及所述第1下侧节段，该第2血管模型包括所述上侧节段、所述分岔点以及所述第2下侧节段；以及轮廓分岔点确定步骤，至少使用所述第1血管模型和所述第2血管模型，确定所述血管的轮廓分岔点。
接下来，作为非限定性的例子说明实施方式，并且表示在以下的附图中。
附图说明
图1a是表示血管与中心线的示意图。
图1b是表示血管与中心线的示意图。
图1c是表示血管与中心线的示意图。
图2是分岔构造的概图。
图3是一实施方式的装置的概图。
图4是表示按照一实施方式进行的处理的概要的流程图。
图5是血管的另一概图。
图6是分岔构造的概图。
图7是分岔构造的概图。
图8是分岔构造的概图。
具体实施方式
一实施方式涉及的医用画像处理装置具备：血管模型生成单元，关于至少包括分岔点、相对该分岔点位于血流的上游侧的上侧节段、相对该分岔点位于血流的下游侧的第1下侧节段、相对该分岔点位于血流的下游侧且不同于所述第1下侧节段的第2下侧节段在内的被检体的血管，至少生成包括所述上侧节段、所述分岔点、所述第1下侧节段在内的第1血管模型和包括所述上侧节段、所述分岔点、所述第2下侧节段在内的第2血管模型；以及轮廓分岔点确定单元，至少使用所述第1血管模型和所述第2血管模型，确定所述血管的轮廓分岔点。
如图3中示意地示出的那样，一实施方式的图像处理装置构成为实施前一段落记载的方法。在这种情况下，该装置具备作为个人计算机(PC)或者工作站的处理装置32，处理装置32与显示设备34和CT扫描仪36连接，在该情况下，与作为计算机键盘和鼠标的用户输入设备或者设备38连接。处理装置还与数据保存部40连接。
关于患者或者其它的对象也可以使用能够进行三维CT测量的任意的适当的类型的CT扫描仪。图1的实施方式关于CT扫描数据进行了记载，但也可以在替代实施方式中使用生成任意的适当的类型的图像数据的、其它的任意适当的类型的扫描仪，例如，可以使用适当的形式的MR数据，如果进行适当的预处理，则也可以使用PET数据或者数字减影X射线血管造影数据、超声波图像数据。进而，执行本实施方式的医用画像处理的装置能够通过医用工作站等以单体来实现，但也可以将执行该医用画像处理的功能安装到医用画像诊断装置（Ｘ射线ＣＴ装置、ＭＲＩ装置、ＰＥＴ、Ｘ射线诊断装置、超音波诊断装置等）。另外，本实施方式的医用画像处理也可以由计算机依据专用程序来实现。
处理装置32是提供自动或者半自动地处理图像数据的处理用资源的装置，并且具备中央处理装置(CPU)10，中央处理装置(CPU)10能够以加载各种软件模块或者其它的软件组件并执行的方式进行动作，该各种软件模块或者其它的软件组件构成为边参照图4边执行以下详述的方法。
处理装置包括:图像绘制单元42、用户交互单元44、轮廓计算单元46、图像数据处理单元48。
图像数据处理单元构成为决定表现上侧节段的至少一部分和第1下侧节段的至少一部分的第1血管模型的路径，并且决定表现上侧节段的至少一部分和第2下侧节段的至少一部分的第2血管模型的路径，上侧节段在分岔点处被分开为第1下侧节段和第2下侧节段。
轮廓计算单元46构成为决定轮廓分岔点的位置，并且关于轮廓分岔点以上的至少一个位置或者全部的位置，保证第1血管模型的轮廓和第2血管模型的轮廓实质相等。
处理装置32中还包括个人计算机的其它的组件，该个人计算机的其它的组件包括硬盘驱动器、RAM、ROM、数据总线、包括各种设备驱动的操作系统、包括显卡的硬件设备。为了易于观察附图而未在图1中示出这样的组件。
在图3的实施方式中，处理装置32具备2.0GHz的芯片组、2×4芯处理器，提供使用9线程的多线程环境，设置有6.0GB的RAM。但是，也可以应用于任意的适当的CPU或其它的组件中。
在图3的实施方式中，在通过CT扫描仪36扫描后，由处理装置32从CT扫描仪36接收图像数据组37，并保存到数据保存部40，之后由处理装置进行处理。图3的实施方式示出的扫描仪36是CT扫描仪，但在替代实施方式中，为了得到图像数据组，能够使用其它的任意的适当的类型的扫描仪。
在图3的实施方式的变形例中，处理装置32不是从扫描仪36而是从远程的数据保存部(未图示)接收图像数据组。远程的数据保存部在某个期间，将从多个不同的扫描仪得到的大量的不同的数据组与关联的患者的数据一起保存。数据保存部既可以是保存大量的患者的数据的服务器，也可以形成医用图像保存通信系统(PACS)、例如东芝医疗系统株式会社的Rapideye(注册商标)系统的一部分。
如图4的流程图中示出概要那样，图3的系统构成为进行具有一连串的阶段的处理。
在处理的第1阶段，图像数据处理单元48从数据保存部40接收体图像数据组37。在此情况下，数据组37具有作为注入了造影剂的患者的CT图像数据而使用CT扫描仪36获得的数据。
图像数据处理单元48执行用于自动地识别表示血管构造、例如冠状动脉构造的数据的既知的手法。
接下来，图像数据处理单元48在阶段52中，执行生成沿着血管的中心的中心线点，并连结这些点或者进行拟合而得到中心线的既知的手法。中心线形成树状构造，该树状构造例如如图2中示意的那样，具备血管、节段、节点以及分岔点。
在接下来的阶段，轮廓计算单元46关于构造的血管模型的每一个计算轮廓平面。各血管模型表示从一个节点至另一节点的路径，例如表示从根节点至叶节点的路径，或者表示这样的路径的一部分。各血管模型有时包括几个节段和分岔点。各血管的始点和终点，例如以包括关心点或者关心区域的方式，既可以由用户手动地选择，也可以自动地选择。能够理解某个血管模型也可以与其它的血管模型共有一个或者多个节段以及分岔点的情况。
所记载的实施方式的特征在于，关于沿着血管模型的中心线的某个特定的位置处的轮廓平面，根据比该特定位置靠上的中心线点和靠下的中心线点的双方的位置来决定，而替代仅根据该特定位置处的中心线来决定。该特征在图5中示意地示出。图5示出图1a至1c示出的血管构造4、计算而得到的中心线2、石灰化特征部(calcificationfeature：钙化特征部)6。
在图5的情况下，轮廓计算单元46在比要计算轮廓的中心线点靠上以及靠下的位置，选择沿着中心线而处于规定的距离处的中心线点70、72。接着，轮廓计算单元46将轮廓平面74决定为在连结中心线点70、72的直线76的中间点垂直的平面。由此，某个特定位置处的轮廓平面不仅由该特定位置处的直线的方向，而且还根据血管的区域中的中心线的方向来决定。这样，不会对中心线的位置造成影响，并且能够获得将轮廓的朝向的变化设为平滑的效果。
轮廓计算单元46为了在图6所示那样的分岔点或者其附近的位置决定构造的轮廓面，能够利用与图5关联地记载的处理。图6示出在分岔点86处，被分为第1下侧节段82和第2下侧节段84的上侧节段80。第1血管模型用于表现上侧节段80和第1下侧节段82。第2血管模型用于表现上侧节段80和第2下侧节段84。
轮廓计算单元46，最初针对第1和第2血管模型的双方，计算与分岔点86或者其附近的位置有关的轮廓面。
例如，参照图6，轮廓计算单元46针对第1血管模型和第2血管模型的双方决定分岔点附近的位置104的轮廓面100。关于位置104处的第1血管模型的轮廓面，选择位置比104靠上的点106和靠下的点108，用直线112连结点106和点108，关于轮廓面100，通过决定为在直线112的中间点垂直于直线112的面来确定。同样地，关于第2血管模型，在计算在相同的位置104的轮廓面102的情况下，选择比该位置104靠上的点106和靠下的点110，用直线114连结点106和110，将轮廓面102决定为在直线114的中间点垂直于直线114的面。
但是，下侧的点108、110比分岔点86靠下，所以点108和110在第1血管模型和第2血管模型中不同。因此，如图6所示意的那样，轮廓面100、102针对第1血管模型以及第2血管模型决定为朝向不同的方向。
在接下来的阶段56，轮廓计算单元决定针对轮廓面100、102的每一个的垂直方向，将这些轮廓面的法线之间的角度124与预先决定的阈值进行比较。在该情况下，阈值是阈值角。图7中示出法线之间的角度124。如果角度124大于预先决定的阈值，则该位置处的轮廓面100、102被视为不同的轮廓面。垂直于轮廓面100的方向能够视为第1轮廓参数，垂直轮廓面102的方向能够视为第2轮廓参数。
接着，关于处于比中心线分岔点86靠上的下一位置，针对第1血管模型和第2血管模型的双方重复进行阶段54中的轮廓面决定手续。针对轮廓面的法线之间的角度再次与阈值比较，如果该角度大于预先决定的阈值，则关于比中心线分岔点86靠上的下一位置重复进行手续。
如果行进到比中心线分岔点86靠上，则针对两个血管的轮廓面的法线之间的角度逐渐消失而变为零。假设，在阶段56中，如果决定为针对轮廓面的法线之间的角度124为预先决定的阈值以下，则该位置被决定为是轮廓分岔点120。因此，轮廓分岔点被配置在针对轮廓面的法线之间的角度被认为以能够允许的程度靠近零的位置。在图8中示出轮廓分岔点120，可知轮廓分岔点120与血管分岔点86不同。
如果决定了轮廓分岔点，则手续转移到阶段58，关于剩下的全部的位置决定轮廓面。所记载的实施方式的特征在于，关于比轮廓分岔点120靠上的位置，针对一方的第1血管模型以及第2血管模型决定一次轮廓面，而对另一方的第1血管模型以及第2血管模型同样地设定针对这些位置的轮廓面。因此，轮廓面的一部分有时不是计算2次以上，而是仅计算1次，由此，能够节约时间，减轻处理负担。
在图6的例子中，针对从上侧节段上的从分岔点86远离的位置的轮廓面88、90对于双方的血管模型相同，仅计算1次，而接下来在双方的血管模型中使用。比轮廓分岔点靠下的轮廓面92、94、96、98、100、102在不同的血管模型应用中各自不同。
接着，关于各个位置，在这些位置根据轮廓面上出现的图像数据决定轮廓。在此，所记载的实施方式的特征在于:关于比轮廓分岔点120靠上的位置，针对一方的第1血管模型以及第2血管模型中的一方决定一次的轮廓，这些位置的轮廓针对另一方的第1血管模型以及第2血管模型也同样地分配。因此，轮廓的一部分不是计算2次以上，而是仅计算1次，由此，在此也能够节约时间，并减轻处理负担。
在接下来的阶段60，将表现构造体的图像与表现构造体的边界的轮廓一起，能够在图像绘制单元42的管理下根据需要而在图形用户接口(GUI)上绘制并显示。轮廓也可以在图像内接合，表现构造体的连续的边界。
血管模型的一方或者双方的轮廓以及中心线的位置，在特定的实施方式中，之后也可以由用户手动地编辑。例如，在用于冠状动脉解析时，用户观察表现血管的图像数据，也许想要移动血管模型的中心线或者轮廓的位置，以更好地表现血管的路径或者边界。虽然基于GUI的详细的内容，但用户也许通过操作鼠标，例如使用托拽操作来变更中心线或者轮廓的位置。用户交互单元44接收用户的输入，针对轮廓以及/或者中心线进行所请求的变更。
关于图4而记载的实施方式的特征在于，能够将轮廓计算单元构成为:关于比轮廓分岔点靠上的位置，通过自动地链接某位置处的第1血管模型的轮廓和相同位置处的第2血管模型的轮廓，从而第2血管模型的轮廓对应于第1轮廓模型中的变更而自动地改变。因此，如果用户在比轮廓分岔点靠上的位置手动地变更第1血管模型的轮廓，则对第2血管模型的对应的轮廓自动地进行相同的变更。对照地，关于比轮廓分岔点靠下的位置，用户能够相互独立地变更第1以及第2血管的每一个的轮廓。
关于图4而记载的实施方式设置为能够将沿着树状构造的轮廓中的变化平滑化。该实施方式对于高效地编辑树状构造也是有用的。利用图像解析技术抽出的原始的血管中心线不需要修正，就连在分岔点，轮廓看起来似乎也平滑地连结。用户能够简单地编辑共有的轮廓，并且针对靠近分岔点的每个血管用户也能够分别地编辑轮廓。通过避免不必要地重复生成轮廓，从而能够提高性能。
在所记载的实施方式中，将沿着中心线而选择的位置处的轮廓面决定为垂直于连结比该所选择的位置靠上以及靠下的中心线点的直线，但在替代实施方式中，也可以将比所选择的位置靠上以及靠下的多个中心线点用于决定轮廓面。例如，也可以利用比所选择的位置靠上的两个以上的中心线点、比所选择的位置靠下的两个以上的中心线点，另外，也可以使用任意的适当的曲线拟合技术而拟合为曲线。这样的实施方式中的轮廓面被视为垂直于所拟合的曲线。
在关联图6至8而记载的实施方式中，通过比较针对轮廓面的法线之间的角度与预先决定的阈值来决定轮廓分岔点的位置。在替代实施方式中，使用另一轮廓参数决定轮廓或者轮廓面的类似性，并由此决定轮廓分岔点的位置。例如，在一实施方式中，比较轮廓的形状，决定沿着中心线的一串点处的、各轮廓位置之间的平均间隔。接着，轮廓分岔点被决定为平均间隔为阈值以下的位置。在其它的实施方式中也能够使用其它的任意的适当的轮廓参数。
进而，在另一替代实施方式中，轮廓分岔点被决定在比分岔点靠上预定的距离的位置。因此，在这样的实施方式中，不需要比较不同的血管模型的轮廓参数也能够决定轮廓分岔点的位置。
参照将单一的节段在分岔点处分开为两个下侧节段的例而说明了实施方式。然而，实施方式也能够用于在更复杂的树构造中将包括多个分岔点以及/或者来自各分岔点的多个节段的结构模型化。
另外，通过关联CT测量而记载了实施方式，但所记载的实施方式，例如也能够用于对从包括MRI、X射线血管造影术、或者PET摄影法在内的其它的医用图像摄影装置获得的图像数据进行模型化。实施方式能够使用于对包括至少一个分岔点的任意的适当的构造进行模型化，实施方式不仅限于对血管构造进行模型化。例如，实施方式也可以使用于对构成淋巴系统或者神经系统的一部分的血管进行模型化。
在另一实施方式中，提供一种获得平滑的血管中心线的方法，在该方法中，将具备能够(不是连续且平滑地)不连续地改变方向的、由被连结的节段构成的血管中心线的血管树设为输入。在被连结的节段的顶点(分岔点)，在顶点的两侧和从顶点离开的部位，能够将位置配置到中心线上。将这样的位置连结起来的矢量能够视作顶点处的轮廓的法线，这样的轮廓的集合能够描绘血管壁。关于其它的轮廓，这些轮廓与它们的法线一起描绘原始数据的有用且可使用的抽象化了的血管的平滑的路径。
在另一实施方式中，提供一种利用血管中心线来识别分岔点处的血管节段的方法，该方法通过以下步骤来进行：决定比上述分岔点靠上且位于平滑的中心线上的位置，关于比所决定的位置靠下的各血管节段，检查轮廓的法线和这些法线的方向的差异，在该方向的差异小到能够允许的程度的部位，识别为各个节段在此相连。因此，就连在分岔点也能够清清楚楚地识别血管。
为了进行本方法，通过使用具有可执行的计算机可读命令的计算机程序，实施方式能够实施特定的功能，这是本领域技术人员非常清楚的。计算机程序的功能能够利用硬件(例如，使用CPU)来实施。但是，实施方式还可以通过一个或者多个ASIC(面向特定用途的集成电路)、或者通过组合硬件和软件的结构来实现。
本说明书记载了特定的单元，但在替代实施方式中，能够通过单一的单元、处理用资源、或者其它的计算机来提供一个或者多个这些单元的功能，或者，能够组合两个以上的单元来提供在单一的单元中提供的功能。单元能够单独地设置，或者也能够组合设置。
在本说明书中使用的情况下，上侧和下侧、以及靠上和靠下这些用语在分岔点仅意图表示一侧或者另一侧。上侧以及下侧、以及靠上和靠下这些用语并不用于表示特定的任意的坐标系中的相对高度或者位置。例如，虽然基于患者或者其它的对象的朝向，但上侧节段或者血管也可以位于下侧节段或者血管之上、之下或者旁边，不作限定。
另外，在上述实施方式中，也可以区别地显示轮廓分岔点和中心线的分岔点。进而，在上述实施方式中，以将对象设为血管来确定轮廓分岔点的情况作为了例子。但是，不限于该例子，关于具有分岔点的管状构造物（例如，肺等）都能够应用本实施方式的手法。
记载了特定的实施方式，但这些实施方式仅作为例子而提示，不用于限定本发明的范围。实际上，本说明书中记载的新的方法以及系统能够通过各种各样的其它的方式来实现。进而，在本说明书记载的方法以及系统的方式中，能够不脱离本发明的宗旨地进行各种省略、置换以及变更。所附的权利要求书及其均等范围意图包括能够进入发明的范围的这样的方式以及修正。

Claims (12)

1.一种医用画像处理装置，其特征在于，具备：
血管模型生成单元，关于至少包括分岔点、相对该分岔点位于血流的上游侧的上侧节段、相对该分岔点位于血流的下游侧的第1下侧节段、相对该分岔点位于血流的下游侧且不同于所述第1下侧节段的第2下侧节段在内的被检体的血管，至少生成第1血管模型和第2血管模型，该第1血管模型包括所述上侧节段、所述分岔点、所述第1下侧节段，该第2血管模型包括所述上侧节段、所述分岔点、所述第2下侧节段；以及
轮廓分岔点确定单元，至少使用所述第1血管模型和所述第2血管模型，确定所述血管的轮廓分岔点。
2.根据权利要求1所述的医用画像处理装置，其特征在于，
所述轮廓分岔点确定单元使用位于比第1分岔点还上游侧的位置的所述上侧节段的中心线来确定所述血管的轮廓分岔点，在该第1分岔点，所述上侧节段的中心线被分岔为所述第1下侧节段的中心线和所述第2下侧节段的中心线。
3.根据权利要求1或2所述的医用画像处理装置，其特征在于，
所述轮廓分岔点确定单元使用所述血管的第1轮廓面和所述血管的第2轮廓面来确定所述血管的轮廓分岔点，所述血管的第1轮廓面是使用所述第1血管模型来决定的，所述血管的第2轮廓面是使用所述第2血管模型来决定的。
4.根据权利要求3所述的医用画像处理装置，其特征在于，
所述轮廓分岔点确定单元使用在所述上侧节段中选择出的位置和在所述第1下侧节段中选择出的位置来决定所述第1轮廓面，使用在所述上侧节段中选择出的位置和在所述第2下侧节段中选择出的位置来决定所述第2轮廓面。
5.根据权利要求3所述的医用画像处理装置，其特征在于，
所述轮廓分岔点确定单元使用所述上侧节段的规定位置处的中心线和所述第1下侧节段的规定位置处的中心线来决定所述第1轮廓面，使用所述上侧节段的规定位置处的中心线和所述第2下侧节段的规定位置处的中心线来决定所述第2轮廓面。
6.根据权利要求3至5中任意一项所述的医用画像处理装置，其特征在于，
所述轮廓分岔点确定单元根据比较所述第1轮廓面的法线与所述第2轮廓面的法线所成的角和规定的阈值得到的结果来确定所述血管的轮廓分岔点。
7.根据权利要求3至5中任意一项所述的医用画像处理装置，其特征在于，
所述轮廓分岔点确定单元根据所述第1轮廓面和所述第2轮廓面的类似性来确定所述血管的轮廓分岔点。
8.根据权利要求1至7中任意一项所述的医用画像处理装置，其特征在于，
在所述轮廓分岔点确定单元中，在变更了相对所述轮廓分岔点位于血流的上流侧的所述第1血管模型的轮廓面、以及相对所述轮廓分岔点位于血流的上流侧的所述第2血管模型的轮廓面中的一方的情况下，也按照该变更来变更另一方的所述轮廓面。
9.根据权利要求1至8中任意一项所述的医用画像处理装置，其特征在于，
所述轮廓分岔点确定单元相互独立地变更相对所述轮廓分岔点位于血流的下游侧的所述第1血管模型的轮廓面、和相对所述轮廓分岔点位于血流的下游侧的所述第2血管模型的轮廓面。
10.根据权利要求1至9中任意一项所述的医用画像处理装置，其特征在于，
还具备显示单元，该显示单元区别显示所述轮廓分岔点和所述第1分岔点。
11.一种医用画像诊断装置，其特征在于，具备：
血管模型生成单元，关于至少包括分岔点、相对该分岔点位于血流的上游侧的上侧节段、相对该分岔点位于血流的下游侧的第1下侧节段、相对该分岔点位于血流的下游侧且不同于所述第1下侧节段的第2下侧节段在内的被检体的血管，至少生成第1血管模型和第2血管模型，该第1血管模型包括所述上侧节段、所述分岔点以及所述第1下侧节段，该第2血管模型包括所述上侧节段、所述分岔点以及所述第2下侧节段；以及
轮廓分岔点确定单元，至少使用所述第1血管模型和所述第2血管模型，确定所述血管的轮廓分岔点。
12.一种医用图像处理方法，其特征在于包括：
血管模型生成步骤，关于至少包括分岔点、相对该分岔点位于血流的上游侧的上侧节段、相对该分岔点位于血流的下游侧的第1下侧节段、相对该分岔点位于血流的下游侧且不同于所述第1下侧节段的第2下侧节段在内的被检体的血管，至少生成第1血管模型和第2血管模型，该第1血管模型包括所述上侧节段、所述分岔点以及所述第1下侧节段，该第2血管模型包括所述上侧节段、所述分岔点以及所述第2下侧节段；以及
轮廓分岔点确定步骤，至少使用所述第1血管模型和所述第2血管模型，确定所述血管的轮廓分岔点。
CN201310312557.4A 2012-07-24 2013-07-24 医用图像处理装置、医用图像诊断装置以及医用图像处理方法 Active CN103565455B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US13/556,705 US9020232B2 (en) 2012-07-24 2012-07-24 Per vessel, vessel tree modelling with shared contours
US13/556,705 2012-07-24

Publications (2)

Publication Number Publication Date
CN103565455A true CN103565455A (zh) 2014-02-12
CN103565455B CN103565455B (zh) 2016-04-20

Family

ID=49994934

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310312557.4A Active CN103565455B (zh) 2012-07-24 2013-07-24 医用图像处理装置、医用图像诊断装置以及医用图像处理方法

Country Status (3)

Country Link
US (1) US9020232B2 (zh)
JP (1) JP6169912B2 (zh)
CN (1) CN103565455B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015136853A1 (ja) * 2014-03-14 2015-09-17 テルモ株式会社 画像処理装置、画像処理方法及びプログラム
JP6780936B2 (ja) * 2016-01-13 2020-11-04 ザイオソフト株式会社 医用画像処理装置、医用画像処理方法、及び医用画像処理プログラム
EP3270355B1 (en) * 2017-01-27 2019-07-31 Siemens Healthcare GmbH Determining a complexity value of a stenosis or a section of a vessel
CA3100495A1 (en) 2018-05-16 2019-11-21 Benevis Informatics, Llc Systems and methods for review of computer-aided detection of pathology in images
CN111429404B (zh) * 2020-03-04 2020-12-15 国科赛赋(深圳)新药研发科技有限公司 一种用于心脑血管检测的成像系统和方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002237996A (ja) * 2001-02-13 2002-08-23 Hitachi Medical Corp 医療用ｘ線透視撮影装置
US20020118869A1 (en) * 2000-11-28 2002-08-29 Knoplioch Jerome F. Method and apparatus for displaying images of tubular structures
US20040066958A1 (en) * 2002-10-08 2004-04-08 Chen Shiuh-Yung James Methods and systems for display and analysis of moving arterial tree structures
CN1954340A (zh) * 2004-05-18 2007-04-25 皇家飞利浦电子股份有限公司 使用三维可变形网格模型对物体的三维树状管形表面自动分段的图像处理系统
US20100021025A1 (en) * 2008-07-22 2010-01-28 Jacobus Adrianus Engelbertus Johannes Hof Method, Apparatus and Computer Program for Quantitative Bifurcation Analysis in 3D Using Multiple 2D Angiographic Images
CN102208105A (zh) * 2010-03-31 2011-10-05 富士胶片株式会社 医学图像处理技术

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6842638B1 (en) * 2001-11-13 2005-01-11 Koninklijke Philips Electronics N.V. Angiography method and apparatus
JP4421203B2 (ja) 2003-03-20 2010-02-24 株式会社東芝 管腔状構造体の解析処理装置
US8086000B2 (en) * 2007-08-27 2011-12-27 Pie Medical Imaging B.V. Method, apparatus and computer program for quantitative bifurcation analysis on angiographic images
JP5305821B2 (ja) 2008-10-10 2013-10-02 株式会社東芝 医用画像処理装置及び医用画像診断装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020118869A1 (en) * 2000-11-28 2002-08-29 Knoplioch Jerome F. Method and apparatus for displaying images of tubular structures
JP2002237996A (ja) * 2001-02-13 2002-08-23 Hitachi Medical Corp 医療用ｘ線透視撮影装置
US20040066958A1 (en) * 2002-10-08 2004-04-08 Chen Shiuh-Yung James Methods and systems for display and analysis of moving arterial tree structures
CN1954340A (zh) * 2004-05-18 2007-04-25 皇家飞利浦电子股份有限公司 使用三维可变形网格模型对物体的三维树状管形表面自动分段的图像处理系统
US20100021025A1 (en) * 2008-07-22 2010-01-28 Jacobus Adrianus Engelbertus Johannes Hof Method, Apparatus and Computer Program for Quantitative Bifurcation Analysis in 3D Using Multiple 2D Angiographic Images
CN102208105A (zh) * 2010-03-31 2011-10-05 富士胶片株式会社 医学图像处理技术

Also Published As

Publication number Publication date
JP2014026653A (ja) 2014-02-06
US20140029816A1 (en) 2014-01-30
CN103565455B (zh) 2016-04-20
US9020232B2 (en) 2015-04-28
JP6169912B2 (ja) 2017-07-26

Similar Documents

Publication Publication Date Title
US11615894B2 (en) Diagnostically useful results in real time
US10993687B2 (en) Synthetic data-driven hemodynamic determination in medical imaging
US11138733B2 (en) Vascular characteristic determination with correspondence modeling of a vascular tree
US10470730B2 (en) Creating a vascular tree model
JP6553675B2 (ja) 脈管構造を数値的に評価するためのシステムおよび方法
US10134129B2 (en) Method and system for hemodynamic computation in coronary arteries
US10595807B2 (en) Calculating a fractional flow reserve
EP2946321B1 (en) Diagnostically useful results in real time
CN106659399B (zh) 使用患病和假想正常解剖学模型中的流计算的冠状动脉狭窄的非侵入功能评价的方法和系统
WO2018133118A1 (zh) 血流状态分析系统及方法
CN105184086A (zh) 用于冠状动脉中改进的血液动力学计算的方法及系统
CN103565455B (zh) 医用图像处理装置、医用图像诊断装置以及医用图像处理方法
CN117115150B (zh) 用于确定分支血管的方法、计算设备和介质
WO2015150320A1 (en) Segmentation of tubular organ structures
KR101479577B1 (ko) 심근 및 심혈관 정보의 통합 분석 방법
CN109564685B (zh) 鲁棒的肺叶分割
CN117036530B (zh) 基于跨模态数据的冠状动脉血流储备分数预测方法及装置
Uus et al. Design of a 0D image-based coronary blood flow model
CN109199418A (zh) 至少一个特征变量的可视化
Mukai et al. Construction Method of 3D Aorta Model

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for 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: 20160708

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.