医用图像处理装置、X射线CT装置及医用图像处理方法与流程

本申请享受2017年2月2日提出申请的日本国专利申请号2017－017849的优先权的利益，该日本国专利申请的全部内容援引于本申请。此外，2018年1月31日提出申请的日本国专利申请号2018－015731的全部内容援引于本申请。

实施方式涉及医用图像处理装置、x射线ct装置及医用图像处理方法。

背景技术

以往，在肺疾患的诊断中，进行评价呼吸功能的处理。例如，使用肺活量计(肺活量测定法(spirometry))观察肺整体的换气量、换气曲线(肺活量曲线)，从而进行慢性堵塞性肺疾患(chronicobstructivepulmonarydisease：copd)等的诊断。但是，为了早期发现copd等的肺疾患，仅用肺活量测定法有时无法充分地诊断。另外，通过肺活量测定法，难以确认肺的哪个部分的功能降低。

技术实现要素：

本发明要解决的课题在于，提供能够高精度地检测呼吸功能有异常的肺的区域的医用图像处理装置、x射线ct装置及医用图像处理方法。

实施方式的医用图像处理装置具有提取部、计算部、检测部、输出控制部。提取部从沿着时间序列拍摄到的三维医用图像数据(data)中，提取多个与形成肺的肺叶及亚区域中的至少一方对应的区域。计算部对于所提取到的多个区域的每个区域，计算与呼吸功能有关的物理指标值。检测部通过将所述多个区域各自的所述物理指标值的随时间变化互相比较，检测所述多个区域中的与所述呼吸功能有关的异常区域。输出控制部输出表示所述异常区域的信息。

效果

根据实施方式的医用图像处理装置，能够高精度地检测呼吸功能有异常的肺的区域。

附图说明

图1是表示第1实施方式的x射线ct装置的构成的一例的图。

图2a～图2e是用于说明第1实施方式的提取功能的处理的图。

图3a及图3b是用于说明第1实施方式的检测功能的处理的图。

图4是用于说明第1实施方式的输出控制功能的处理的图。

图5是表示第1实施方式的x射线ct装置进行的处理步骤的流程图(flowchart)。

图6是用于说明第1实施方式的变形例的检测功能的处理的图。

图7是用于说明第2实施方式的检测功能的处理的图。

图8a～图8d是用于说明第2实施方式的检测功能的处理的图。

图9是用于说明第2实施方式的输出控制功能的处理的图。

图10是表示第2实施方式的x射线ct装置进行的处理步骤的流程图。

图11是用于说明其他的实施方式的x射线ct装置进行的处理的图。

图12是表示其他的实施方式的医用图像处理装置的构成例的框图。

图13是表示其他的实施方式的提供信息处理服务(service)的服务器装置的构成例的框图(block)图。

具体实施方式

本发明要解决的课题在于，提供能够高精度地检测呼吸功能有异常的肺的区域的医用图像处理装置、x射线ct装置及医用图像处理方法。

实施方式的医用图像处理装置具备处理电路。处理电路从沿着时间序列拍摄到的三维医用图像数据中，提取多个与形成肺的肺叶及亚区域中的至少一方对应的区域。处理电路对于所提取到的多个区域的每个区域，计算与呼吸功能有关的物理指标值。处理电路通过将所述多个区域各自的所述物理指标值的随时间变化互相比较，检测所述多个区域中的与所述呼吸功能有关的异常区域。处理电路输出表示所述异常区域的信息。

以下，参照附图，对医用图像处理装置及x射线ct(computedtomography)装置的实施方式详细地进行说明。在以下的实施方式中，将拍摄被检体的x射线ct图像数据的x射线ct装置作为例子进行说明。但是，实施方式并不限定于此，例如，也能够广泛应用于能够拍摄三维的x射线图像数据的x射线诊断装置、或者能够对三维的医用图像数据进行处理的医用图像处理装置(计算机(computer))。

(第1实施方式)

图1是表示第1实施方式的x射线ct装置1的构成的一例的图。如图1所示，第1实施方式的x射线ct装置1具有架台10、诊视床装置20及控制台(console)30。

架台10是对被检体p(患者)照射x射线，检测透射被检体p后的x射线，并输出至控制台30的装置，具有x射线照射控制电路11、x射线产生装置12、检测器13、数据收集电路(das：dataacquisitionsystem)14、旋转框架(frame)15及架台驱动电路16。

旋转框架15是将x射线产生装置12和检测器13以夹着被检体p相对的方式支撑，并通过后述的架台驱动电路16在以被检体p为中心的圆轨道上高速地旋转的圆环状的框架。

x射线照射控制电路11是作为高电压产生部，对x射线管12a供给高电压的装置，x射线管12a使用从x射线照射控制电路11供给的高电压，产生x射线。x射线照射控制电路11通过后述的扫描控制电路33的控制，调整对x射线管12a供给的管电压、管电流，从而调整对被检体p照射的x射线量。

另外，x射线照射控制电路11进行楔子(wedge)12b的切换。另外，x射线照射控制电路11通过调整准直器(collimator)12c的开口量，来调整x射线的照射范围(扇形(fan)角、锥(corn)角)。另外，本实施方式也可以是，操作人员手动切换多种楔子的情况。

x射线产生装置12是产生x射线，并将所产生的x射线向被检体p照射的装置，具有x射线管12a、楔子12b及准直器12c。

x射线管12a是通过由未图示的高电压产生部供给的高电压对被检体p照射x射线束的真空管，随着旋转框架15的旋转，对被检体p照射x射线束。x射线管12a产生具有扇形角及锥角并扩展的x射线束。例如，在x射线照射控制电路11的控制下，x射线管12a能够在全(full)重构用时在被检体p的整个周围连续照射x射线、或在半(half)重构用时在能够进行半重构的照射范围(180度+扇形角)内连续照射x射线。另外，在x射线照射控制电路11的控制下，x射线管12a能够在预先设定的位置(管球位置)间歇地照射x射线(脉冲(pulse)x射线)。另外，x射线照射控制电路11也能够对从x射线管12a照射的x射线的强度进行调制。例如，x射线照射控制电路11在特定的管球位置，使从x射线管12a照射的x射线的强度变强，在特定的管球位置以外的范围，使从x射线管12a照射的x射线的强度变弱。

楔子12b是用于调节从x射线管12a照射的x射线的x射线量的x射线过滤器。具体而言，楔子12b是使从x射线管12a照射的x射线透射并衰减以使从x射线管12a向被检体p照射的x射线成为预先决定的分布的过滤器(filter)。例如，楔子12b是将铝(aluminium)加工成规定的目标角度、规定的厚度的过滤器。另外，楔子12b也被称为楔形过滤器(wedgefilter)、领结式过滤器(bow-tiefilter)。

准直器12c是在x射线照射控制电路11的控制下，将通过楔子12b调节了x射线量的x射线的照射范围缩小所用的缝隙。

架台驱动电路16通过旋转驱动旋转框架15，使x射线产生装置12和检测器13在以被检体p为中心的圆轨道上转动。

检测器(x射线检测器)13是检测透射被检体p后的x射线的2维阵列(array)型检测器(面检测器)，配置多通道(channel)量的x射线检测元件而成的检测元件列沿着z轴方向排列有多列。具体而言，第1实施方式中的检测器13，具有沿着z轴方向排列有320列等多列的x射线检测元件，例如，能够检测在包括被检体p的肺、心脏的范围等大范围内透射被检体p后的x射线。另外，z轴方向与架台10为非倾斜(tilt)时的状态下的旋转框架15的旋转中心轴方向对应。

数据收集电路14是das，从由检测器13检测到的x射线的检测数据中，收集投影数据。例如，数据收集电路14对于通过检测器13检测到的x射线强度分布数据，进行放大处理、a/d转换处理、通道间的灵敏度校正处理等而生成投影数据，并将所生成的投影数据对后述的控制台30发送。例如，在旋转框架15的旋转中从x射线管12a连续照射x射线的情况下，数据收集电路14收集整个周围量(360度量)的投影数据群。另外，数据收集电路14将所收集到的各投影数据与管球位置建立对应，并发送至后述的控制台30。管球位置为表示投影数据的投影方向的信息。另外，通道间的灵敏度校正处理也可以由后述的前处理电路34进行。

诊视床装置20是承载被检体p的装置，如图1所示，具有诊视床驱动装置21和顶板22。诊视床驱动装置21使顶板22向z轴方向移动，使被检体p在旋转框架15内移动。顶板22是载置被检体p的板。

另外，架台10例如一边使顶板22移动一边使旋转框架15旋转而执行对被检体p进行螺旋状地扫描(scan)的螺旋式(helical)扫描。或者，架台10在使顶板22移动后将被检体p的位置固定不动的状态下使旋转框架15旋转而执行按圆轨道对被检体p进行扫描的常规(conventional)扫描。或者，架台10使顶板22的位置以一定间隔移动而执行用多个扫描区域(area)进行常规扫描的调强(stepandshoot)方式。

控制台30是受理操作人员对x射线ct装置1的操作，并且使用通过架台10所收集到的投影数据重构x射线ct图像数据的装置。控制台30如图1所示，具有输入电路31、显示器(display)32、扫描控制电路33、前处理电路34、存储电路35、图像重构电路36及处理电路37。输入电路31、显示器32、扫描控制电路33、前处理电路34、存储电路35、图像重构电路36及处理电路37连接为能够相互通信。

输入电路31具有x射线ct装置1的操作人员进行各种指示、各种设定的输入所用的鼠标(mouse)、键盘(keyboard)、跟踪球(trackball)、开关(switch)、按钮(button)、操纵杆(joystick)等，将从操作人员受理的指示、设定的信息转送给处理电路37。例如，输入电路31从操作人员受理x射线ct图像数据的摄影条件、重构x射线ct图像数据时的重构条件、针对x射线ct图像数据的图像处理条件等。另外，输入电路31受理用于选择对被检体p的检查的操作。另外，输入电路31受理用于指定图像上的部位的指定操作。

显示器32是供操作人员参照的监视器，在处理电路37的控制下，将从x射线ct图像数据生成的图像数据显示给操作人员、或者显示经由输入电路31从操作人员受理各种指示、各种设定等所用的gui(graphicaluserinterface)。另外，显示器32显示扫描计划的计划画面、扫描中的画面等。

扫描控制电路33在处理电路37的控制之下，控制x射线照射控制电路11、架台驱动电路16、数据收集电路14及诊视床驱动装置21的动作，从而控制架台10中的投影数据的收集处理。具体而言，扫描控制电路33对收集定位图像(扫描图像)的定位摄影及收集诊断中使用的图像的主摄影(主扫描)中的投影数据的收集处理分别进行控制。

例如，扫描控制电路33将x射线管12a固定在0度的位置(相对于被检体为正视方向的位置)，使顶板22定速移动同时连续地进行摄影，从而摄影到2维的扫描图像。或者，扫描控制电路33将x射线管12a固定在0度的位置，使顶板22断续地移动同时与顶板移动同步地反复进行断续地摄影，从而摄影到2维的扫描图像。这里，扫描控制电路33不仅能够从相对于被检体p为正视方向摄影到定位图像，而且能够从任意的方向(例如，侧视方向等)摄影到定位图像。

另外，扫描控制电路33通过收集针对被检体的整周量的投影数据，进行三维的x射线ct图像数据(体数据)的摄影。例如，扫描控制电路33通过螺旋式扫描或者非螺旋式扫描来收集针对被检体p的整周量的投影数据。另外，扫描控制电路33通过以比主摄影低的线量收集整周量的投影数据，也能够拍摄到三维的扫描(scano)图像。

另外，扫描控制电路33通过使体(volume)数据的拍摄持续进行规定期间，从而能够进行拍摄到沿着时间序列的多个体数据的动态(dynamic)体扫描(也称为“动态扫描”)。例如，通过在被检体p正在进行某关节的运动期间持续收集整周量的投影数据达规定期间，从而能够拍摄到用规定的帧(frame)率(体率；volumerate)重构的多个体数据。另外，通过动态扫描拍摄到的时间序列的体数据，被称为4维x射线ct图像数据或者4dct图像数据。

前处理电路34对通过数据收集电路14所生成的投影数据，进行对数转换处理、偏移校正、灵敏度校正及射束硬化校正等的校正处理，生成校正后的投影数据。具体而言，前处理电路34对于通过数据收集电路14所生成的定位图像的投影数据及通过主摄影所收集到的投影数据，分别生成校正后的投影数据，并保存在存储电路35中。

存储电路35存储通过前处理电路34所生成的投影数据。具体而言，存储电路35存储通过前处理电路34所生成的定位图像的投影数据及通过主摄影所收集的诊断用的投影数据。另外，存储电路35存储通过后述的图像重构电路36所生成的x射线ct图像数据等。另外，存储电路35适当存储后述的处理电路37的处理结果。

图像重构电路36使用存储电路35存储的投影数据，重构x射线ct图像数据。具体而言，图像重构电路36根据定位图像的投影数据及诊断中所用的图像的投影数据，分别重构x射线ct图像数据。这里，作为重构方法，有各种方法，例如列举出反投影处理。另外，作为反投影处理，列举出例如基于fbp(filteredbackprojection)法的反投影处理。或者，图像重构电路36使用逐次近似法，也能够重构x射线ct图像数据。另外，图像重构电路36通过针对x射线ct图像数据进行各种图像处理，生成图像数据。并且，图像重构电路36，将重构出的x射线ct图像数据、通过各种图像处理生成的图像数据保存于存储电路35。另外，图像重构电路36是图像重构部的一例。

另外，图像重构电路36重构通过动态扫描所拍摄到的时间序列的三维医用图像数据(4dct图像数据)。例如，图像重构电路36通过以规定的帧率重构持续规定期间收集到的整周量的投影数据，从而重构出沿着时间序列的多个体数据。由此，例如，对某关节的运动的样子进行表示的连续的多个帧(时相/相位(phase))的体数据(4dct图像数据)得以重构。图像重构电路36将重构出的4dct图像数据保存于存储电路35。

处理电路37通过控制架台10、诊视床装置20及控制台30的动作，进行x射线ct装置1的整体控制。具体而言，处理电路37控制扫描控制电路33，从而控制通过架台10进行的ct扫描。另外，处理电路37通过控制图像重构电路36，从而控制控制台30中的图像重构处理、图像生成处理。另外，处理电路37进行控制以将由存储电路35存储的各种图像数据显示于显示器32。

另外，处理电路37如图1所示，执行提取功能371、计算功能372、检测功能373及输出控制功能374。这里，例如，图1所示的处理电路37的构成要素即提取功能371、计算功能372、检测功能373及输出控制功能374执行的各处理功能，以通过计算机可执行的程序(program)的方式记录于存储电路35中。处理电路37是通过从存储电路35读出并执行各程序来实现与各程序对应的功能的处理器(processor)。换言之，读出了各程序的状态的处理电路37，具有图1的处理电路37内所示的各功能。另外，关于提取功能371、计算功能372、检测功能373及输出控制功能374执行的各处理功能，在后面叙述。

另外，在本实施方式中，对通过单一的处理电路37实现以下说明的各处理功能的情况进行了说明，但也可以是将多个独立的处理器组合而构成处理电路37，并由各处理器执行程序来实现功能的情况。

在上述说明中使用的“处理器”这一说法，意味着例如cpu(centralprocessingunit)、gpu(graphicsprocessingunit)、或者面向特定用途的集成电路(applicationspecificintegratedcircuit：asic)、可编程(programmable)逻辑器件(device)(例如，简单可编程逻辑器件(simpleprogrammablelogicdevice：spld)、复杂可编程逻辑器件(complexprogrammablelogicdevice：cpld)及现场可编程门阵列(fieldprogrammablegatearray：fpga))等的电路。处理器通过读出并执行在存储电路35中所保存的程序来实现功能。另外，也可以代替在存储电路35中保存程序，而在处理器的电路内中直接装入程序而构成。在此情况下，处理器通过读出并执行装入到电路内的程序来实现功能。另外，本实施方式的各处理器不限于每个处理器作为单一的电路而构成的情况，也可以将多个独立的电路组合而构成为1个处理器并实现其功能。进而，也可以将各图中的多个构成要素统合到1个处理器中实现其功能。

以上，对第1实施方式的x射线ct装置1的构成进行了说明。在该构成的基础上，第1实施方式的x射线ct装置1为了高精度地检测呼吸功能有异常的肺的区域，而执行以下的各处理功能。

提取功能371，从沿着时间序列拍摄到的三维医用图像数据中，提取多个与形成肺的肺叶及亚区域中的至少一方对应的区域。例如，提取功能371从存储电路35读出在存储电路35中所存储的4dct图像数据。然后，提取功能371基于ct值，从所读出的4dct图像数据中提取与肺整体对应的区域。然后，提取功能371通过区域划分处理(分割)从肺整体的区域中提取多个与亚区域对应的区域。另外，提取功能371是提取部的一例。

图2a至图2e是用于说明第1实施方式的提取功能371的处理的图。图2a中例示了从正面观看到的左右的肺的示意图。另外，在图2b中例示了右外侧面(右肺的外侧面)的示意图。另外，在图2c中，例示了右内侧面(右肺的内侧面)的示意图。另外，在图2d中，例示了左外侧面(左肺的外侧面)的示意图。另外，在图2e中，例示了左内侧面(左肺的内侧面)的示意图。

如图2a至图2e所示，例如提取功能371从肺整体的区域中，提取与多个亚区域分别对应的多个区域。这里，所谓的亚区域，是形成肺的肺叶的区域。若举具体例，提取功能371，通过使对肺中的多个亚区域的位置进行表示的模板(template)图像变形为被检体的肺形状，从而提取多个与被检体的亚区域对应的区域。另外，模板图像，预先与肺的解剖学的特征点的三维的位置关系建立对应。

另外，上述的提取功能371的处理只不过是一例，并不限定于上述的例子。例如，在上述的例子中，对提取功能371提取与肺的亚区域对应的区域的情况进行了说明，但也能够提取与肺叶对应的区域。另外，在上述的说明中，对使用模板图像的情况进行了说明，但不限于此，例如，也能够对4dct图像数据中的肺的活动进行解析并提取活动不同的部位作为肺叶、亚区域的边界。另外，提取肺叶、亚区域的方法，可以应用现有的任何技术。

计算功能372对于所提取到的多个区域的每个区域，计算与呼吸功能有关的物理指标值(参数(parameter))。例如，计算功能372计算通过提取功能371提取到的多个区域各自的体积。举出具体例时，计算功能372对于4dct图像数据中包括的各个时相的体数据，计算多个区域各自的体积。由此，计算功能372针对多个区域的每个区域，计算各时相的体积。另外，计算功能372是计算部的一例。

另外，上述的计算功能372的处理只不过是一例，并不限定于上述的例子。例如，在上述的例子中，对计算功能372计算体积作为物理指标值的情况进行了说明，但不限于此，例如也能够计算表面积、比表面积(表面积除以体积而得到的值)、或者ct值。另外，各区域的ct值例如是各区域中包括的全部像素的ct值的平均值，ct值相应于各区域中包括的空气的量而变化，因此可被作为表示呼吸功能的指标来利用。即，计算功能372能够计算体积、表面积、比表面积及ct值中的至少一个。

检测功能373基于物理指标值的随时间变化，检测多个区域中的、呼吸功能有异常的异常区域。例如，检测功能373通过将多个区域各自的物理指标值的随时间变化互相进行比较，从而检测异常区域。另外，检测功能373是检测部的一例。换言之，检测功能373通过将多个区域各自的物理指标值的随时间变化互相进行比较，从而检测多个区域中的与呼吸功能有关的异常区域。

例如，检测功能373通过对多个区域的各时相的参数的值描绘拟合(fitting)曲线(curve)，从而生成表示各区域的参数的随时间变化的曲线(曲线)。并且，检测功能373从所生成的各曲线检测极大点及极小点。这里，曲线的极大点与“最大吸气”对应，极小点与“最大呼气”对应。并且，检测功能373检测多个区域各自的参数的随时间变化的曲线的倾向不同的区域，作为异常区域。

图3a及图3b是用于说明第1实施方式的检测功能373的处理的图。图3a中例示对区域a～区域d中的体积的随时间变化进行表示的曲线(graph)(曲线)。另外，在图3b中，例示对区域e～区域h中的体积的随时间变化进行表示的曲线。在图3a及图3b中，纵轴表示各区域的体积(volume)，横轴表示时间(time)。另外，区域a～h分别与亚区域对应。另外，在图3a中，区域b～d的最大吸气的时刻(timing)是“t1”，最大呼气的时刻是“t2”。另外，区域a的最大呼气的时刻是“t3”，最大吸气的时刻是“t4”。

如图3a所示，检测功能373基于最大吸气及最大呼气的时刻的偏移，检测异常区域。例如，检测功能373计算各区域的最大呼气的时刻的差，并互相进行比较，从而检测最大呼气的时刻的偏移。在图3a中，区域a的最大呼气的时刻t3(极小点50)，与其他的区域b～d各自的最大呼气的时刻t2不同。具体而言，即使使得区域a的最大吸气的时刻与各区域b～d的最大吸气的时刻t1大致一致，下一个的最大呼气的时刻也有偏移。同样地，区域a的最大吸气的时刻t4(极大点51)，与其他的区域b～d各自的最大吸气的时刻不同。在此情况下，检测功能373检测区域a作为异常区域。换言之，检测功能373在最大呼气(或者最大吸气)的时刻的偏移为阈值以上的情况下，检测为异常区域。

如图3b所示，检测功能373基于最大吸气时与最大呼气时之间的体积差(peaktopeak：p－p)，检测异常区域。例如，检测功能373计算各区域的p－p，并互相比较，从而检测异常区域。在图3b中，区域f的p－p与其他的区域e、g、h各自的p－p相比较较小，暗示未将呼气吐净(极小点52)。在此情况下，检测功能373将区域f检测为异常区域。

这样，检测功能373将对多个区域各自的参数的随时间变化进行表示的曲线的倾向互相进行比较。并且，检测功能373将与其他的区域的倾向不同的区域，检测为异常区域。

另外，上述的检测功能373的处理只不过是一例，并不限定于上述的例子。例如，纵轴也可以用各区域的相对于最大体积的比例[％]表示。另外，在上述的例子中，对检测功能373基于最大吸气及最大呼气的时刻的偏移检测异常区域的情况进行了说明，但实施方式并不限定于此。例如，检测功能373也能够使用最大吸气及最大呼气中的、任一方的时刻的偏移，来检测异常区域。另外，例如，检测功能373能够将无法与正弦曲线(sin曲线)近似的区域或者曲线没有周期性的区域检测为异常区域。

另外，在上述的例子中，对将曲线的倾向不同的区域检测为异常区域的情况进行了说明，但实施方式并不限定于此。关于其他的检测方法，作为变形例而在后面叙述。

输出控制功能374输出表示异常区域的信息。例如，输出控制功能374在检测到异常区域的时相的三维医用图像数据上，使异常区域强调显示。另外，输出控制功能374是输出控制部的一例。

图4是用于说明第1实施方式的输出控制功能374的处理的图。在图4中，例示显示于显示器32的被检体的肺的图像。另外，在图4中，例示了作为异常区域而检测到“区域a”的情况。

如图4所示，例如，输出控制功能374使基于三维医用图像数据的肺的图像显示于显示器32。这里，该图像，是基于检测到区域a的时相(t3或者t4)的三维医用图像数据的体绘制图像(或者mpr(multiplanarreconstructions)图像等)。换言之，输出控制功能374显示基于检测到异常区域的时相的三维医用图像数据的显示用图像。并且，输出控制功能374将图像中的、与通过检测功能373检测到的异常区域对应的位置强调显示(以与其他的区域不同的颜色显示等)。

这样，输出控制功能374将表示异常区域的信息显示于显示器32。另外，图4只不过是一例，并不限定于图示的例子。例如，输出控制功能374也可以不必须是基于检测到异常区域的时相的三维医用图像数据的图像，而能够在任意的时相的图像上显示异常区域。

另外，在图4中，对肺的图像被显示为静止图像的情况进行了说明，但实施方式并不限定于此。例如，输出控制功能374也可以将肺的图像显示为动画。在动画显示中，输出控制功能374，也可以以随着接近检测到异常区域的时相而强调显示异常区域，并随着远离检测到异常区域的时相而使强调显示消失(不强调)的方式，显示异常区域。

另外，在图4中，对表示异常区域的信息被显示于显示器32的情况进行了说明，但实施方式并不限定于此。例如，输出控制功能374，也能够作为“区域a”等的文本数据来显示，还能够通过文本的阅读功能而进行声音输出。具体而言，输出控制功能374也可以将对检测到异常区域的时相(t3)进行表示的信息(时刻)、对通过检测功能373检测到的偏移的大小进行表示的信息(例如，t2与t3之差)进行表示(或者声音输出)。另外，表示异常区域的信息被输出的输出目的地，不限于显示器32、声音输出用器件，可以是任意的存储介质、其他的装置(报表(report)制作用的应用程序(application)等)。

图5是表示第1实施方式的x射线ct装置1进行的处理步骤的流程图。图5所示的处理步骤，例如通过由操作人员输入意在使检测异常区域的处理开始的指示而开始。

如图5所示，在步骤(step)s101中，处理电路37判定是否是处理时刻。例如，处理电路37在由操作人员输入意在使检测异常区域的处理开始的指示时，判定为是处理时刻，并开始步骤s102以后的处理。另外，在步骤s101为否定的情况下，处理电路37不使步骤s102以后的处理开始，而为待机状态。

在步骤s101为肯定时，在步骤s102中，提取功能371从4dct图像数据中，提取形成肺的多个区域。例如，提取功能371基于ct值，从4dct图像数据中提取与肺整体对应的区域。并且，提取功能371通过区域划分处理(分割(segmentation))从肺整体的区域中提取多个与亚区域对应的区域。

在步骤s103中，计算功能372针对多个区域分别计算与呼吸功能有关的参数(物理指标值)。例如，计算功能372针对4dct图像数据中包括的各个时相的体数据，计算多个区域各自的体积。由此，计算功能372针对多个区域分别计算各时相的体积。

在步骤s104中，检测功能373基于参数的随时间变化，从多个区域中检测异常区域。例如，检测功能373将对多个区域各自的体积的随时间变化进行表示的曲线的倾向互相进行比较。并且，检测功能373基于最大吸气及最大呼气的时刻的偏移，检测异常区域。

在步骤s105中，输出控制功能374显示异常区域。例如，输出控制功能374，在检测到异常区域的时相的三维医用图像数据上，使异常区域强调显示。

这样，x射线ct装置1，受理意在使检测异常区域的处理开始的操作人员的指示，执行步骤s102～s105的各处理，并使异常区域显示。另外，图5的内容只不过是一例，并不限定于此。

如上述那样，在第1实施方式的x射线ct装置1中，提取功能371，从沿着时间序列拍摄到的三维医用图像数据中，提取多个与形成肺的肺叶及亚区域中的至少一方对应的区域。计算功能372，对于所提取到的多个区域的每个区域，计算与呼吸功能有关的物理指标值(参数)。检测功能373基于物理指标值的随时间变化，检测多个区域中的、呼吸功能有异常的异常区域。即，检测功能373，通过将多个区域各自的物理指标值的随时间变化互相进行比较，从而检测异常区域。输出控制功能374输出表示异常区域的信息。由此，x射线ct装置1能够高精度地检测呼吸功能有异常的肺的区域。

例如，x射线ct装置1并不对被检体p的肺的状态整体进行评价，而以肺叶、亚区域的单位进行评价。由此，x射线ct装置1不仅能够检测出被检体p的肺是否有异常并提示给操作人员，还能够检测肺的哪部分(区域)有异常并提示给操作人员。

(第1实施方式的变形例1)

检测功能373的处理并不限定于上述的实施方式，也可以通过其他的实施方式来实现。例如，检测功能373也能够将基于参数的随时间变化的最大吸气时与最大呼气时的差分的评价值(index)比其他的区域低的区域，检测为异常区域。

例如，检测功能373使用下述的式(1)，计算评价值f[％]。另外，在式(1)中，vi与最大吸气时的体积对应。另外，ve与最大呼气时的体积对应。

【数式1】

即，检测功能373通过将多个区域各自的最大吸气时的体积与最大呼气时的体积的差分除以最大吸气时的体积并作为百分率，从而计算评价值f[％]。并且，检测功能373将对于多个区域分别计算出的评价值f[％]互相进行比较。比较的结果，检测功能373将与其他的区域相比较评价值f[％]较低的区域检测为异常区域。

这样，检测功能373将基于参数的随时间变化的最大吸气时与最大呼气时的差分的评价值比其他的区域低的区域，检测为异常区域。另外，上述的评价值f[％]不限于体积，也可以使用表面积、比表面积或者ct值来计算。

另外，在变形例1中，如果有最大吸气时及最大呼气时的体数据的话，就能够检测异常区域。换言之，检测功能373通过使用在最大吸气时(或者最大呼气时)让被检体p屏息期间拍摄到的体数据，不使用4dct图像数据就能够检测异常区域。

(第1实施方式的变形例2)

另外，例如，检测功能373能够将参数的随时间变化的包括最大呼气时的期间的微分系数比其他的区域小的区域，检测为异常区域。

例如，作为肺疾患的重要的观察结果之一，列举出呼气是否被吐净。在呼气未被吐净的情况下，认为最大呼气时的周边的曲线变得平缓。

因此，检测功能373对多个区域的每个区域，从表示参数的随时间变化的曲线确定极小点(向下凸的拐点)。并且，检测功能373计算包括所特定的极小点的时相的规定期间的曲线的微分系数。检测功能373将对于多个区域分别计算出的微分系数互相进行比较。比较的结果，检测功能373将与其他的区域相比较而言微分系数较小的区域检测为异常区域。

这样，检测功能373将参数的随时间变化的包括最大呼气时的期间的微分系数比其他的区域小的区域检测为异常区域。另外，用于计算微分系数的期间能够任意地设定。另外，检测功能373也能够计算包括最大吸气时的规定期间的微分系数。

(第1实施方式的变形例3)

另外，例如，检测功能373能够通过对左右的肺中成对的区域彼此进行比较，检测异常区域。

例如，关于左右的肺中的哪个有异常，存在会作为被检体p的自觉症状而获得的情况。在这种情况下，通过对左右的肺中成对的区域彼此进行比较，能够进行异常区域的检测。

例如，检测功能373将各区域的参数的随时间变化与各区域的成对的区域的参数的随时间变化进行比较。这里，形成左右的肺的肺叶、亚区域，不一定位于左右对称的位置。因此，检测功能373是与距作为对象的区域的左右对称的位置最近的区域进行比较。

作为一例，在被检体p自觉左肺有异常而右肺正常的情况下，检测功能373将左肺的各区域与右肺中距对称的位置最近的区域进行比较。并且，检测功能373将与右肺的区域的参数的随时间变化(曲线)的倾向大幅背离的左肺的区域，检测为异常区域。

这样，检测功能373将左右的肺中的与正常一方相比较有大幅背离的区域，检测为异常区域。另外，在上述的说明中，说明了将亚区域彼此进行比较的情况，但不限于此，也能够以肺叶的单位进行比较。

(第1实施方式的变形例4)

另外，例如，检测功能373通过将多个区域中的成为基准的基准区域与各区域进行比较，能够检测异常区域。

例如，检测功能373将多个区域中的、ct值具有标准的值的区域设定为基准区域(正常区域)。并且，检测功能373将多个区域各自参数除以基准区域的参数来计算相对值。这里，如果各区域是正常的，则各区域的相对值为大致相同程度的值。因此，检测功能373将多个区域各自的相对值互相比较，将与其他的区域相比较相对值脱离的区域检测为异常区域。另外，正常的亚区域的ct值预先设定。

这样，检测功能373通过将多个区域中的成为基准的基准区域与各区域进行比较，从而能够检测异常区域。另外，在上述的说明中，对将ct值具有标准的值的区域作为基准区域的情况进行了说明，但并不限定于此。例如，检测功能373也可以将最大吸气与最大呼气间的体积变化最大的区域、或者由操作人员所指定的区域设定为基准区域。即，检测功能373将被认为正常的亚区域的区域设定为基准区域。

(第1实施方式的变形例5)

另外，例如，检测功能373能够将参数的随时间变化的曲线的规定时相的值未到达按正弦曲线变化的情况下的规定时相的值的区域，检测为异常区域。

图6是用于说明第1实施方式的变形例的检测功能373的处理的图。图6中例示对某区域中的体积的随时间变化予以表示的曲线。在图6中，纵轴表示各区域的体积(volume)，横轴表示时间(time)。

如图6所示，例如，检测功能373为，假定为某区域的体积的随时间变化的曲线成为正弦曲线，并计算某时相的体积的值作为阈值。在图6的例子中，检测功能373计算最大吸气时与最大呼气时的中间点(50％)的值作为阈值。并且，检测功能373基于某区域的体积的随时间变化的曲线是否到达阈值，来检测异常区域。在图6的例子中，某区域的曲线，在中间点53处未到达50％。因此，检测功能373将该区域检测为异常区域。

这样，检测功能373将参数的随时间变化的曲线的规定时相的值未到达以正弦曲线变化的情况下的规定时相的值的区域，检测为异常区域。另外，图6只不过是一例，并不限定于上述的说明。例如，在图6中，对最大吸气时与最大呼气时的中间的时相被设定为阈值的情况进行了说明，但不限于此，能够设定为任意的时相。另外，在图6中，对将一个时间点设定为阈值的情况进行了说明，但不限于此，能够将多个时相设定为阈值。

(第2实施方式)

在第1实施方式中，对x射线ct装置1检测与异常所在的肺叶、亚区域对应的区域的情况进行了说明，但实施方式并不限定于此。例如，x射线ct装置1对于对肺叶、亚区域供给空气的支气管也能够进行检测异常所在的区域的处理。

第2实施方式的x射线ct装置1，具有与图1所例示的x射线ct装置1同样的构成，处理电路37的处理的一部分不同。因此，在第2实施方式中，以与第1实施方式不同的点为中心进行说明，对于具有与第1实施方式中已说明的构成同样的功能的点，将说明省略。

例如，提取功能371进一步提取多个与对多个区域分别供给空气的支气管对应的支气管区域。这里，通常，人体的气管，被分支为对左右的肺供给空气的主支气管(左主支气管，右主支气管)，进而分支为对肺叶供给空气的肺叶支气管、对亚区域供给空气的亚区域支气管(末端部)。例如，提取功能371将与对亚区域供给空气的亚区域支气管对应的区域，作为支气管区域而提取。另外，提取方法可以应用使用了模板图像的方法等、现有的任意的技术。

另外，提取功能371作为支气管区域而提取的范围，并不仅限定于亚区域支气管。例如，提取功能371除了亚区域支气管以外，还可以将包括肺叶支气管、左右的主支气管的范围作为支气管区域而提取。但是，为了取与亚区域的对应关系，优选提取至少包括亚区域支气管的区域。

例如，计算功能372对于所提取到的多个支气管区域，分别计算物理指标值。举出具体例，计算功能372对于支气管区域，与亚区域的区域同样地、计算表面积、比表面积或者ct值作为物理指标值。另外，计算功能372也可以计算各支气管区域的截面积，作为支气管区域的物理指标值。该截面积例如是与支气管区域的长轴正交的截面的面积。

例如，检测功能373通过将有对应关系的区域及支气管区域组合并进行比较，从而检测异常区域。这里，所谓的有对应关系的区域及支气管区域，表示亚区域的区域、及包括对该亚区域供给空气的亚区域支气管的支气管区域。

图7及图8a～图8d是用于说明第2实施方式的检测功能373的处理的图。在图7及图8a～图8d中，纵轴表示各区域的体积(volume)，横轴表示时间(time)。

在图7所示的例子中，对检测功能373使用有对应关系的区域及支气管区域的峰值(peak)的偏移来检测异常区域的情况进行说明。这里，在有对应关系的亚区域与亚区域支气管之间，存在着亚区域支气管先膨胀然后亚区域膨胀的解剖学的关系。该膨胀的时间差，被认为在各个亚区域大致一定。因此，检测功能373针对有对应关系的亚区域与亚区域支气管的每个对(pair)(组合)，计算参数的随时间变化的曲线的峰值的偏移。在图7所示的例子中，检测功能373计算支气管区域的曲线(实线)的最大吸气的时刻与亚区域的区域的曲线(虚线)的最大吸气的时刻的差分。并且，检测功能373将每个对的峰值的偏移互相比较，并将与其他的对相比较而言峰值的偏移较大的对的亚区域的区域，检测为异常区域。

在图8a～图8d所示的例子中，进一步对使用支气管扩张剂检测异常区域的情况进行说明。图8a中，例示了支气管扩张剂的非投放时的亚区域的区域的体积的随时间变化(曲线)。另外，图8b中，例示了支气管扩张剂的非投放时的支气管区域的体积的随时间变化(曲线)。另外，图8c中，例示了支气管扩张剂的投放时的亚区域的区域的体积的随时间变化(曲线)。另外，图8d中，例示了支气管扩张剂的投放时的支气管区域的体积的随时间变化(曲线)。另外，区域a～d表示亚区域的区域。另外，区域a’～d’表示对各区域a～d供给空气的支气管区域。

如图8a及图8c所示，在支气管扩张剂的投放时/非投放时，亚区域的体积的随时间变化看不到显著的变化。另一方面，如图8b及图8d所示，支气管扩张剂非投放时的区域a’的p－p较小(图8b)，但支气管扩张剂投放时的区域a’的p－p较大(图8d)。在这种情况下，暗示着肺气肿性的copd。

因此，检测功能373在支气管扩张剂的投放时/非投放时，计算亚区域的区域及支气管区域的p－p，并对计算出的p－p彼此进行比较。并且，检测功能373，在通过支气管扩张剂的投放，亚区域的区域的p－p没有变化而支气管区域的p－p变大的情况下，将该亚区域的区域检测为异常区域。

这样，检测功能373通过对有对应关系的亚区域的区域及支气管区域组合并进行比较，从而检测异常区域。另外，图7及图8a～图8d只不过是一例，并不限定于上述的说明。例如，在暗示着支气管的异常的情况下，检测功能373也可以检测呼吸功能有异常的异常支气管区域。

例如，输出控制功能374将亚区域的区域的图像、及对该亚区域的区域供给空气的支气管区域的图像同时显示。

图9是用于说明第2实施方式的输出控制功能374的处理的图。图9中例示了显示器32上显示的被检体的肺的图像。另外，在图9中，例示了检测为异常区域到“区域a”的情况。

如图9所示，例如，输出控制功能374使基于三维医用图像数据的肺的图像显示在显示器32上。并且，输出控制功能374将图像中的、通过检测功能373检测到的区域a和与对区域a供给空气的支气管对应的区域a’(支气管区域)的位置强调显示(以与其他的区域不同的颜色显示等)。

这样，输出控制功能374将表示异常区域的信息显示于显示器32。另外，图9只不过是一例，并不限定于图示的例子。例如，输出控制功能374也能够将亚区域的区域的随时间变化的曲线和对该亚区域的区域供给空气的支气管区域的随时间变化的曲线同时显示。

图10是表示第2实施方式的x射线ct装置1进行的处理步骤的流程图。图10所示的处理步骤中的、步骤s201、步骤s202a及步骤s203a的各处理，与图5所示的步骤s101、步骤s102及步骤s103的各处理是同样的，所以将说明省略。

如图10所示，步骤s201中为肯定时，在步骤s202b中，提取功能371从4dct图像数据提取多个支气管区域。例如，提取功能371使用模板图像，提取多个与支气管对应的支气管区域。

在步骤s203b中，计算功能372针对每个支气管，计算与呼吸功能有关的参数(物理指标值)。例如，计算功能372对于4dct图像数据中包括的各个时相的体数据，计算多个支气管区域各自的体积。

在步骤s204中，检测功能373将亚区域和支气管的参数组合并解析。例如，检测功能373将有对应关系的区域及支气管区域组合并比较，从而检测异常区域。

在步骤s205中，输出控制功能374显示异常区域及/或异常支气管区域。例如，输出控制功能374在检测到异常区域的时相的三维医用图像数据上，将异常区域和有对应关系的支气管区域强调显示。

这样，x射线ct装置1，受理意在使检测异常区域的处理开始的操作人员的指示，执行步骤s202～s205的各处理，显示异常区域。另外，图10的内容只不过是一例，并不限定于此。例如，在图10中，对步骤s202a、203a的处理和步骤s202b、203b的处理通过并行处理来执行的情况进行了说明，但实施方式并不限定于此。例如，步骤s202a、203a的处理和步骤s202b、203b的处理也可以不必须通过并行处理来执行。即，也能够将步骤s202a、203a的处理和步骤s202b、203b的处理中的任一方先处理。

如上述那样，第2实施方式的x射线ct装置1，对于对肺叶、亚区域供给空气的支气管也进行检测有异常的区域的处理。由此，例如，x射线ct装置1能够更详细地解析被检体p的肺。

另外，通过第1实施方式说明的内容，除了进行与支气管有关的处理这点以外，在第2实施方式中也能够应用。

(其他的实施方式)

除了上述的实施方式以外，可以以各种不同的方式实施。

(呼吸功能正常的情况下的曲线的显示)

除了上述的实施方式以外，也能够进一步生成并显示异常区域的呼吸功能正常的情况下的曲线。

例如，输出控制功能374，显示对异常区域的呼吸功能正常的情况下的物理指标值的随时间变化进行表示的正常曲线。具体而言，输出控制功能374基于异常区域的体积及多个区域中的与异常区域不同的区域的时相而使正常曲线的模板的波形变更，从而生成正常曲线。并且，输出控制功能374显示所生成的正常曲线。

图11是用于说明其他的实施方式的x射线ct装置进行的处理的图。图11中例示依次进行读出与异常波形(区域a)对应的模板波形的处理(s11)、调整模板波形的振幅和时相的处理(s12)及显示区域a的推定正常波形的处理(s13)的情况。另外，在图11中，对在从图3a所示的各区域a～d中检测区域a作为异常区域的情况下将区域a的正常曲线予以显示的处理进行说明。

如图11所示，在s11中，进行读出与异常波形(区域a)对应的模板波形的处理。在存储电路35中预先存储有模板波形。这里，所谓的模板波形，例如是对形成肺的多个区域(肺叶或者亚区域)各自的体积的随时间变化进行表示的代表性的曲线(曲线)。并且，输出控制功能374，在检测为异常区域到区域a的情况下，从存储电路35中读出与区域a对应的区域的模板波形。另外，对各区域的体积的随时间变化进行表示的代表性的曲线，例如通过多个健康人的各区域的体积的随时间变化的统计的处理来决定，但并不限定于此。例如，模板波形也可以单纯利用正弦波形。

在s12中，进行调整模板波形的振幅和时相的处理。例如，输出控制功能374基于被检体p的区域a的体积，将模板波形的振幅a1调整为振幅a2。具体而言，输出控制功能374基于代表性的健康人的区域a的体积与被检体p的区域a的体积的比较，将模板波形的振幅a1调整为振幅a2。更具体而言，输出控制功能374在代表性的健康人的区域a的平均体积为“v0”、被检体p的区域a的平均体积为“v1”的情况下，通过“a2＝a1×v1/v0”来计算振幅a2。

另外，输出控制功能374基于正常的区域b～d的曲线的时相，调整模板波形的时相。具体而言，输出控制功能374使模板波形在时间方向上变形，以使模板波形的最大吸气时的时相与“t1”一致，并使模板波形的最大呼气时的时相与“t2”一致。

这样，输出控制功能374调整模板波形的振幅和时相，从而生成区域a的推定正常波形的曲线。另外，推定正常波形的曲线是正常曲线的一例。

在s13中，进行将区域a的推定正常波形予以显示的处理。例如，输出控制功能374使所生成的区域a的推定正常波形的曲线与区域a的测定波形一起显示于显示器32。另外，区域a的测定波形，是图3a的区域a的曲线。

另外，输出控制功能374也可以不必须与区域a的测定波形一起显示。例如，输出控制功能374既可以仅显示区域a的推定正常波形的曲线，也可以将其与其他的区域b～d的测定波形一起显示。

(医用图像处理装置)

例如，在上述的实施方式中，对处理电路37的构成要素即提取功能371、计算功能372、检测功能373及输出控制功能374执行的各处理功能，在x射线ct装置1中执行的情况进行了说明，但实施方式并不限定于此。例如，上述的各处理功能也可以在工作站(workstation)等的医用图像处理装置中执行。

(x射线诊断装置)

另外，例如，上述的各处理功能也可以在能够拍摄时间序列的三维医用图像数据的x射线诊断装置中执行。例如，如果是从互不相同的3个方向对被检体p照射x射线，并能够拍摄到时间序列的体数据的x射线诊断装置，则能够使用拍摄到的时间序列的体数据执行上述的各处理功能。

图12是表示其他的实施方式的医用图像处理装置200的构成例的框图。医用图像处理装置200例如与个人计算机、工作站等信息处理装置、或者x射线ct装置中包括的控制台装置等的医用图像诊断装置的操作终端对应。

如图12所示，医用图像处理装置200具备输入接口(interface)201、显示器202、存储电路210及处理电路220。输入接口201、显示器202、存储电路210及处理电路220连接为能够相互通信。

输入接口201是鼠标、键盘、触摸面板等用于受理来自操作人员的各种的指示、设定请求的输入装置。显示器202是显示医用图像、或者显示操作人员使用输入接口201输入各种设定请求所用的gui的显示装置。

存储电路210例如是nand(notand)型闪存存储器、hdd(harddiskdrive)，存储用于显示医用图像数据、gui的各种程序、及被该程序所使用的信息。

处理电路220是控制医用图像处理装置200中的处理整体的电子设备(处理器)。处理电路220执行提取功能221、计算功能222、检测功能223及输出控制功能224。处理电路220执行的各处理功能例如以通过计算机可执行的程序的方式记录在存储电路210内。处理电路220读出并执行各程序，从而实现与所读出的各程序对应的功能。提取功能221、计算功能222、检测功能223及输出控制功能224，能够执行与图1所示的提取功能371、计算功能372、检测功能373及输出控制功能374基本上同样的处理。

例如，提取功能221从沿着时间序列拍摄到的三维医用图像数据中，提取多个与形成肺的肺叶及亚区域中的至少一方对应的区域。另外，计算功能222对于所提取到的多个区域的每个区域，计算与呼吸功能有关的物理指标值。另外，检测功能223对多个区域各自的物理指标值的随时间变化进行互相比较，从而检测多个区域中的、与呼吸功能有关的异常区域。另外，输出控制功能224输出表示异常区域的信息。由此，医用图像处理装置200能够高精度地检测呼吸功能有异常的肺的区域。

(x射线诊断装置)

另外，例如，上述的各处理功能也可以在能够拍摄到时间序列的三维医用图像数据的x射线诊断装置中执行。例如，如果是从互不相同的3个方向对被检体p照射x射线，并能够拍摄到时间序列的体数据的x射线诊断装置，则能够使用拍摄到的时间序列的体数据，执行上述的各处理功能。

(作为网络(network)上的服务(service)的提供)

另外，例如，上述的实施方式的处理，能够作为经由网络的信息处理服务(云(crowd)服务)而提供。

图13是表示提供其他的实施方式的信息处理服务的服务器装置的构成例的框图。如图13所示，例如，在提供信息处理服务的服务中心(center)设置服务器装置300。服务器装置300与操作终端301连接。另外，服务器装置300经由网络302与多个客户(client)终端303a、303b、···、303n连接。另外，服务器装置300及操作终端301也可以经由网络302连接。另外，在不将多个客户终端303a、303b、···、303n予以区别而进行总称的情况下，记载为“客户终端303”。

操作终端301是操作服务器装置300的人员(操作人员)利用的信息处理终端。例如，操作终端301具备鼠标、键盘、触摸面板等用于受理来自操作人员的各种的指示、设定请求的输入装置。另外，操作终端301具备显示图像、或者显示操作人员使用输入装置输入各种设定请求所用的gui的显示装置。操作人员通过操作操作终端301，能够将各种的指示、设定请求对服务器装置300发送或者阅览服务器装置300内部的信息。另外，网络302是因特网、wan(wideareanetwork)、lan(localareanetwork)等任意的通信网。

客户终端303是利用信息处理服务的利用者所操作的信息处理终端。这里，利用者例如是从业于医疗机构的医师、技师等的医疗从业者。例如，客户终端303与个人计算机(personalcomputer)、工作站等信息处理装置、或者x射线ct装置中包括的控制台装置等的医用图像诊断装置的操作终端对应。客户终端303具有能够利用通过服务器装置300提供的信息处理服务的客户功能。另外，该客户功能以通过计算机可执行的程序的方式预先记录在客户终端303中。

服务器装置300具备通信接口310、存储电路320及处理电路330。通信接口310、存储电路320及处理电路330连接为能够相互通信。

通信接口310例如是网卡、网络适配器。通信接口310通过与网络302连接，进行服务器装置300与外部装置间的信息通信。

存储电路320例如是nand(notand)型闪存存储器、hdd(harddiskdrive)，存储用于显示医用图像数据、gui的各种程序、及被该程序所使用的信息。

处理电路330是对服务器装置300中的处理整体进行控制的电子设备(处理器)。处理电路330执行提取功能331、计算功能332、检测功能333及输出控制功能334。处理电路330执行的各处理功能例如以计算机可执行的程序的方式记录在存储电路320内。处理电路330读出并执行各程序，从而实现与读出的各程序对应的功能。提取功能331、计算功能332、检测功能333及输出控制功能334，能够执行与图1所示的提取功能371、计算功能372、检测功能373及输出控制功能374基本同样的处理。

例如，利用者操作客户终端303，输入旨在对位于服务中心的服务器装置300发送三维医用图像数据的指示(上载(upload))。被输入该指示时，客户终端303对服务器装置300发送三维医用图像数据。这里，该三维医用图像数据是沿着时间序列拍摄到包括被检体的肺的区域的体数据(4dct图像数据)。

并且，服务器装置300接收从客户终端303发送的三维医用图像数据。并且，在服务器装置300中，提取功能331从沿着时间序列拍摄到的三维医用图像数据中，提取多个与形成肺的肺叶及亚区域中的至少一方对应的区域。另外，计算功能332对于所提取到的多个区域的每个区域，计算与呼吸功能有关的物理指标值。另外，检测功能333对多个区域各自的物理指标值的随时间变化互相进行比较，从而检测多个区域中的与呼吸功能有关的异常区域。另外，输出控制功能334输出表示异常区域的信息。具体而言，输出控制功能334将表示异常区域的信息发送至客户终端303(下载(download))。由此，客户终端303的利用者能够阅览高精度地检测到例如呼吸功能有异常的肺的区域的信息。

即，上述的实施方式的处理，能够作为医用图像处理方法来提供。医用图像处理方法包括：服务器装置300从沿着时间序列拍摄到的三维医用图像数据中，提取多个与形成肺的肺叶及亚区域中的至少一方对应的区域的步骤。医用图像处理方法包括：服务器装置300对于所提取到的多个区域的每个区域，计算与呼吸功能有关的物理指标值的步骤。医用图像处理方法包括：服务器装置300通过将所述多个区域各自的所述物理指标值的随时间变化互相比较，检测所述多个区域中的与所述呼吸功能有关的异常区域的步骤。医用图像处理方法包括：服务器装置300输出表示所述异常区域的信息的步骤。

另外，作为在上述的实施方式及变形例中说明的各处理中的、自动的进行的处理而说明的处理的全部或者一部分也能够通过手动来进行，或者也能够通过公知的方法自动地进行作为手动进行的处理而说明的处理的全部或者一部分。除此之外，关于上述文中、图中所示的处理步骤、控制步骤、具体的名称、包括各种的数据、参数的信息，除了特殊记载的情况以外，能够任意地变更。

另外，在上述的实施方式及变形例中说明的医用图像处理方法，能够通过用个人计算机、工作站等的计算机执行预先准备的医用图像处理程序来实现。该医用图像处理方法能够经由因特网等网络来发布。另外，该医用图像处理方法也能够记录在硬盘、软盘(fd)、cd－rom、mo、dvd等计算机可读取的记录介质中，并通过计算机从记录介质中读出来执行。

根据以上说明的至少一个实施方式，能够高精度地检测呼吸功能有异常的肺的区域。

对本发明的几个实施方式进行了说明，但这些实施方式是作为例子提示的，意图不是限定发明的范围。这些实施方式，能够以其他各种各样的方式实施，在不脱离发明的宗旨的范围内，能够进行各种省略、置换、变更。这些实施方式及其变形，包含于发明的范围及宗旨中，同样地包含于权利要求书所记载的发明及其等同的范围中。