CN105030265B - 图像显示设备、方法以及存储图像显示程序的存储介质 - Google Patents

图像显示设备、方法以及存储图像显示程序的存储介质 Download PDF

Info

Publication number
CN105030265B
CN105030265B CN201510111794.3A CN201510111794A CN105030265B CN 105030265 B CN105030265 B CN 105030265B CN 201510111794 A CN201510111794 A CN 201510111794A CN 105030265 B CN105030265 B CN 105030265B
Authority
CN
China
Prior art keywords
image
images
point
display
peripheral edge
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
Application number
CN201510111794.3A
Other languages
English (en)
Other versions
CN105030265A (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.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
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 Fujitsu Ltd filed Critical Fujitsu Ltd
Publication of CN105030265A publication Critical patent/CN105030265A/zh
Application granted granted Critical
Publication of CN105030265B publication Critical patent/CN105030265B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T1/00—General purpose image data processing
    • G06T1/20—Processor architectures; Processor configuration, e.g. pipelining
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/50—Information retrieval; Database structures therefor; File system structures therefor of still image data
    • G06F16/54—Browsing; Visualisation therefor
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/22—Matching criteria, e.g. proximity measures
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74—Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/75—Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
    • G06V10/752—Contour matching
    • 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/10081—Computed x-ray tomography [CT]
    • 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/20—Special algorithmic details
    • G06T2207/20092—Interactive image processing based on input by user
    • G06T2207/20104—Interactive definition of region of interest [ROI]
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00—Indexing scheme relating to image or video recognition or understanding
    • G06V2201/03—Recognition of patterns in medical or anatomical images
    • G06V2201/031—Recognition of patterns in medical or anatomical images of internal organs

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Multimedia (AREA)
  • Software Systems (AREA)
  • Computing Systems (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Controls And Circuits For Display Device (AREA)
  • Image Analysis (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)

Abstract

本发明涉及图像显示设备、方法以及存储图像显示程序的存储介质，一种图像显示设备，包括：显示装置以及计算机，该计算机被配置成执行使计算机用作为显示控制部件的程序。显示部件同时显示第一CT图像和第二CT图像。第一CT图像和第二CT图像表示同一对象但是在不同时段中获得。显示控制部件进行：根据第一显示模式将第一CT图像和第二CT图像同时显示在显示装置上；确定在第一CT图像中指定的目标部分；以及将第二CT图像的对应部分的显示模式从第一显示模式改变成第二显示模式，该对应部分对应于在第一CT图像中指定的目标部分。

Description

图像显示设备、方法以及存储图像显示程序的存储介质
技术领域
本文中所讨论的实施例涉及图像显示技术。
背景技术
通常，在医疗场所存在下述情况：医生通过使用x射线计算机断层扫描图像(在下文中，称为“CT图像”)来进行诊断。例如，将在不同拍摄时段中获得的CT图像并排显示在显示装置上，并且将CT图像彼此比较以检查病变区域如何改变。
以下专利文献公开了背景技术。
专利文献1：日本专利申请特许公开No.2008-510247
由于患者的位置不同以及患者的呼吸时间不同，导致在不同时段中获得的CT图像中的病变区域的位置和角度可能存在变化。因此，执业医师可能难以准确地识别集中部分的位置。从而，执业医师可能难以对集中病变区域的改变进行比较。
发明内容
根据实施例的一方面，提供了一种图像显示设备，包括：显示装置，该显示装置被配置成同时显示第一CT图像和第二CT图像，第一CT图像和第二CT图像表示同一对象但是在不同时段中获得；以及计算机，该计算机被配置成执行使计算机用作为显示控制部件的程序，该显示控制部件进行：根据第一显示模式将第一CT图像和第二CT图像同时显示在显示装置上；确定在第一CT图像中指定的目标部分；以及将第二CT图像的对应部分的显示模式从第一显示模式改变成第二显示模式，该对应部分对应于在第一CT图像中指定的目标部分。
根据实施例的另一方面，提供了一种图像显示设备，该图像显示设备包括计算机以执行使该图像显示设备进行图像显示处理的程序，该图像显示处理包括：识别在不同时段中获得的两个CT图像中的每个CT图像中的预定部分的外围边缘；在CT图像中的一个CT图像中的外围边缘上确定相对的两个点，所述两个点被定位在穿过在CT图像中的该一个CT图像中指定的指定点的直线上；基于在CT图像中的外围边缘上的每个点的特征量来从另一CT图像中的外围边缘上的点中识别与CT图像中的该一个CT图像的外围边缘上的所述两个点对应的在外围边缘上的两个点；识别在CT图像中的该一个CT图像中的指定点距在CT图像中的该一个CT图像中确定的所述两个点中的每个点的距离；以及识别在该另一CT图像中满足所识别的距离的比率的点。
根据实施例的再一方面，提供了一种同时显示第一CT图像和第二CT图像的图像显示方法，第一CT图像和第二CT图像表示同一对象但是在不同时段中获得，该图像显示方法包括：根据第一显示模式将第一CT图像和第二CT图像同时显示在显示装置上；确定在第一CT图像中指定的目标部分；以及将第二CT图像的对应部分的显示模式从第一显示模式改变成第二显示模式，该对应部分对应于在第一CT图像中指定的目标部分。
附图说明
图1是CT图像拍摄系统的概要结构的图示；
图2是图像显示设备的硬件结构的框图；
图3(a)至图3(d)是在诊断支持程序的操作下的图像显示设备上的显示内容的图示；
图4(a)至图4(e)是用于说明对在不同时段中获得的CT图像的大区域定位的CT图像的图示；
图5是图像显示设备的用于实现局部区域定位和显示控制的功能结构的框图；
图6A和图6B是用于说明特征点/量提取部件的功能的概要的图示；
图7是由特征点/量提取部件进行的特征点/量提取处理的流程图；
图8A和图8B是存储在特征量数据库中的特征量信息的图示；
图9是特征点匹配部件的功能的概要的图示；
图10是由特征点匹配部件进行的特征点匹配处理的流程图；
图11是存储在对应矢量DB中的对应点信息的图示；
图12A和图12B是用于说明对应点估计部件的功能的图示；
图13是由对应点估计部件进行的对应点估计处理的流程图；
图14是由定位处理部件进行的定位处理的流程图；
图15A至图15D是在不同时段中获得的CT图像的第一显示示例的图示；
图16是由显示控制部件进行的第一显示处理的流程图；
图17A至图17D是在不同时段中获得的CT图像的第二显示示例的图示；
图18A是由显示控制部件进行的第二显示处理的流程图的一部分；以及
图18B是由显示控制部件进行的第二显示处理的流程图的另一部分。
具体实施方式
现在参照附图给出对实施例的描述。在描述和附图中，对具有基本相同的功能结构的结构元件给予相同的附图标记，并且省略对其的重复描述。
<第一实施例>
首先参照图1给出对根据第一实施例的包括图像显示设备120的CT(计算机断层扫描)图像拍摄系统100的描述。图1是CT图像拍摄系统100的示例的图示。
CT图像拍摄系统100包括CT设备110和图像显示设备120。CT设备110和图像显示设备120彼此电连接使得其可以与彼此交换数据。
CT设备110通过对通过使用放射性射线扫描患者身体内部而获得的图像数据进行处理来创建患者身体内部的切片图像(剖面图像)。也就是说，CT设备100获得患者身体内部的切片图像。CT设备100将所创建的切片图像(在下文中，称为CT图像)发送至图像显示设备120。
图像显示设备120存储由CT设备110获得的CT图像，并且通过对CT图像进行定位或对准使得CT图像中的图片位置彼此对应而将在不同时段中获得的CT图像同时并排显示。
图像显示设备120包括图像数据库130(在下文中，将该数据库缩写为DB)。此外，图像显示设备120包括在执业医师例如医生基于由CT设备110获得的CT图像来诊断患者时所使用的诊断支持程序140。
诊断支持程序140包括大区域定位程序141，在显示在不同时段中获得的CT图像时，该大区域定位程序141根据在显示不同时段中获得的CT图像时的线性变形来校正CT图像之间的未定位(未对准和未配准)。诊断支持程序140还包括局部区域定位程序142，该局部区域定位程序142通过对已经应用了大区域定位的CT图像中剩余的定位误差进行校正来进行CT图像之间的定位。诊断支持程序140还包括同步显示程序143，该同步显示程序143用于指定已经应用了大区域定位的CT图像中的目标点以及扩大并显示该目标点。
此外，图像显示设备120包括定位数据库(DB)150，该定位数据库(DB)150存储在进行局部区域定位程序142时创建的数据。
以下给出对图像显示设备120的硬件结构的描述。图2是示出了图像显示设备120的硬件结构的框图。图像显示设备120包括：CPU 201、ROM(只读存储器)202以及RAM(随机存取存储器)203。图像显示设备120还包括：存储部件204、连接部件205、显示部件(显示装置)206以及操作部件207。图像显示设备120的这些部件通过总线208互相连接。存储部件204可以是介质驱动单元例如CD-ROM驱动单元，介质驱动单元读取计算机可读存储介质例如CR-ROM存储程序如稍后提及的图像显示程序。
CPU 201通过执行存储在存储部件204中的各种程序的计算机来配置。
ROM 202通过非易失性存储器来配置。ROM 202对CPU 201执行存储在存储部件204中的各种程序所必需的各种程序和数据进行存储。具体地，ROM 202存储引导程序例如BIOS(基本输入/输出系统)、EFI(可扩展固件接口)等。
RAM 203是主存储装置例如DRAM(动态随机存取存储器)、SRAM(静态随机存取存储器)等。RAM 203用作为在CPU 201执行存储部件204中的各种程序时开发的工作区域。
存储部件204对安装在图像显示设备120中的各种程序以及通过执行程序创建的数据进行存储。连接部件205连接至CT设备110以与CT设备110进行数据交换。显示部件206显示由CT设备110获得的并且存储在图像DB 130中的CT图像。操作部件207接收由执业医师例如医生对图像显示设备120进行的各种操作。
接下来，给出对在由执业医师例如医生在图像显示设备120上进行的操作的内容与在图像显示设备120中进行诊断支持程序140的情况下在显示部件206上的显示内容之间的关系的描述。图3(a)至图3(d)是示出了在图像显示设备120中的操作内容和显示内容与诊断支持程序140的操作之间的关系的图示。
当在图像显示设备120中执行诊断支持程序140时，如图3(a)所示，要并排显示不同时段中获得的CT图像的平行显示画面被显示在显示部件206上。在平行显示画面被显示的情况下，当执业医师例如医生选择要用于诊断的CT图像时，诊断支持程序140从图像DB130读取所选择的CT图像。
然后，诊断支持程序140显示所读取的CT图像作为例如药物治疗之前的CT图像(参照图3(b))。
在药物治疗之前的CT图像被显示的情况下，当执业医师例如医生选择药物治疗之后的CT图像时，诊断支持程序140从图像DB 130读取所选择的CT图像。
此时，诊断支持程序140启动大区域定位程序141以将大区域定位应用至所读取的CT图像。此外，诊断支持程序140启动局部区域定位程序142以将局部区域定位应用至已经应用了大区域定位的CT图像。
然后，诊断支持程序140显示已经应用了大区域定位的CT图像作为例如药物治疗之后的CT图像(参照图3(c))。
在药物治疗之前的CT图像和药物治疗之后的CT图像被并排显示的情况下，诊断支持程序140使同步显示程序143被执行。由此，如图3(d)中所示，如果执业医师例如医生将CT图像中的一个CT图像中的任意位置指定为目标点，则另一CT图像中与所指定的点对应的点也被指定为目标点。此外，如果针对CT图像中的一个CT图像中的指定的目标点输入放大显示的指令，则在另一CT图像中也放大并显示指定的目标点。
如上所述，在图像显示设备120中，当并排显示不同时段中获得的CT图像时，CT图像预先经受大区域定位和局部区域定位。此外，在图像显示设备120中，如果针对不同时段中获得并且被并排显示的CT图像中的一个CT图像进行操作例如指定目标点，则在另一CT图像中的对应位置中自动地同步进行同一操作(指定目标点)。从而，例如，执业医师例如医生可以容易地掌握CT图像中的对应位置，由此容易地比较对应位置以检查病变区域如何改变。
以下详细给出对通过包括在诊断支持程序140中的大区域定位程序141、局部区域定位程序142以及同步显示程序143而实现的各种功能的描述。
首先，给出对通过执行大区域定位程序141而实现的针对在不同时段中获得的CT图像的大区域定位的概要的描述。图4(a)至图4(e)是表示针对在不同时段中获得的CT图像的大区域定位的概要的图示。在图4(a)至图4(e)中示出的CT图像中的每个CT图像表示从患者的头侧观看时患者的胸腔部的剖面图。在CT图像中的每个CT图像中，患者的左肺400a、右肺400b、心脏400c以及身体脂肪400d(称为CT图像401)被成像。
CT图像401是在不同时段中获得的CT图像中的比另一CT图像更早获得的一个CT图像(例如，医疗操作或药物治疗之前获得的CT图像)，并且在下文中称为“比较源CT图像”。
另一方面，CT图像411是在不同时段中获得的CT图像中的比另一CT图像更迟获得的一个CT图像(例如，医疗操作或药物治疗之后获得的CT图像)，并且在下文中称为“比较目标CT图像”。
根据大区域定位程序141的执行，对比较目标CT图像411施加线性变形。具体地，旋转比较目标CT图像411以匹配比较源CT图像401从而创建比较目标CT图像412。此外，将经旋转的比较目标CT图像412平行移动(移位)以创建比较目标CT图像413。尽管在图4(a)至图4(e)的示例中未示出，然而在必要时，可以对CT图像应用放大或缩小处理。
如上所述，通过对比较目标CT图像411施加线性变形，比较目标CT图像411可以被粗略地匹配目标源CT图像401。图像420示出了其中已经经受了大区域定位的比较目标CT图像413与比较源CT图像401交叠的图像。
根据图像420，可以解释的是，通过应用大区域定位，比较目标CT图像413可以被粗略地匹配比较源CT图像401。然而，只要应用大区域定位，则如由虚线421和422所示局部位置上的位移就会保留。这是因为在获得比较源CT图像401的时刻与获得比较目标CT图像411的时刻之间患者的呼吸和心跳的时间存在差异。
从而，在图像显示设备120中，通过执行局部区域定位程序142将局部区域定位应用至已经经受了大区域定位的比较目标CT图像413。此外，通过使用局部区域定位的结果而执行同步显示程序143来进行显示控制以使得在比较源CT图像401与已经应用了局部区域定位的比较目标CT图像413之间容易比较。以下给出对由图像显示设备120进行的局部区域定位和显示控制的描述。
首先，给出对由图像显示设备120进行的局部区域定位和显示控制的功能结构的描述。图5是示出了由图像显示设备120进行的局部区域定位和显示控制的框图。
图像显示设备通过CPU 201执行局部区域定位程序142和同步显示程序143来实现以下提及的每个部件的功能。在本实施例中，图像显示设备120的定位DB 150包括特征量DB521和对应矢量DB 522，其中的每个都设置在存储部件204中。
如图5中所示，图像显示设备120包括：显示图像输入部件511、特征点/量提取部件512、特征点匹配部件513、对应点估计部件514以及定位处理部件515。图像显示设备120还包括显示控制部件516。
显示图像输入部件511根据由执业医师例如医生发出的指令来输入比较源CT图像401和已经应用了大区域定位处理的比较目标CT图像413。
特征点/量提取部件512针对包含在所输入的比较源CT图像401中的多个特征点中的每个特征点计算特征量。特征点/量提取部件512还针对包含在所输入的比较目标CT图像413中的多个特征点中的每个特征点计算特征量。
特征点匹配部件513基于比较源CT图像401的特征点中的每个特征点的特征量以及比较目标CT图像413的特征点中的每个特征点的特征量来确定比较目标CT图像413的特征点中的每个特征点中与比较源CT图像401的特征点中的每个特征点对应的特征点。特征点匹配部件513将所确定的对应点作为对应点信息存储在对应矢量DB 522中。
对应点估计部件514针对比较源CT图像401中的在特征点匹配部件513中未确定对应点的每个特征点来估计比较目标CT图像413中的对应点。此外，对应点估计部件514在将所估计的特征点添加至对应点信息之后将该对应点信息存储在对应矢量DB 522中。
定位处理部件515基于存储在对应矢量DB 522中的对应点信息来对应用了大区域定位处理的比较目标CT图像413进行局部区域定位处理。
在比较源CT图像401和已经应用了局部区域定位处理的比较目标CT图像413并排显示的情况下，如果指定了CT图像中的一个CT图像的一部分，则显示控制部件516控制显示器使得在另一CT图像中的对应位置被自动地指定。此外，如果执业医师例如医生针对并排显示的CT图像中的一个CT图像的一部分输入放大显示的指令，则显示控制部件516控制显示器使得在另一CT图像中的对应位置被自动地放大并显示(放大显示)。
接下来，给出对用于由图像显示设备120进行的局部区域定位的各种功能的细节的描述。图6A和图6B是示出了特征点/量提取部件512的功能的概要的图示。
如图6A中所示，特征点/量提取部件512通过沿纵向方向延伸的多个网格线和沿横向方向延伸的多个网格线来划分比较源CT图像401，并且将网格线的每个交点确定为特征点。在图6A的示例中，将由此确定的特征点表示为α1、α2、α3、……。
此外，特征点/量提取部件512基于包含在特征点的周围的预定数量像素的区域(例如，半径为10个像素的区域)中的像素的亮度值来计算比较源CT图像401的每个特征点的特征量。
具体地，特征点/量提取部件512提取包含在特征点α1周围的区域601中的两个任意像素，并且计算这两个像素之间的亮度值的差。然后，如果所计算的亮度差大于预定阈值，则特征点/量提取部件512将特征点α1的特征量设定为“1”。如果所计算的差小于或等于预定阈值，则特征点/量提取部件512将特征点α1的特征量设定为“0”。
在图6A的示例中，提取点A1和点B1作为两个任意像素，并且将特征点α1的特征量设定为“1”，这是因为这两个点之间的亮度差大于阈值。此外，提取另外的点A2和点B2作为两个任意像素，并且将特征点α1的特征量设定为“0”，这是因为这两个点之间的亮度差小于阈值。将从一个特征点周围的区域提取两个像素的处理重复多次(例如，128次)。也就是说，在一个特征点的特征量中包含与预定数量的次数对应的位数(例如，128位)。
类似地，如图6B中所示，特征点/量提取部件512通过沿纵向方向延伸的多个网格线和沿横向方向延伸的多个网格线来划分比较目标CT图像413，并且将网格线的每个交点确定为特征点。在图6B的示例中，将由此确定的特征点表示为β1、β2、β3、……。
此外，特征点/量提取部件512基于包含在特征点的周围的预定数量像素的区域(例如，半径为10个像素的区域)中的像素的亮度值来计算比较目标CT图像413的每个特征点的特征量。
具体地，特征点/量提取部件512提取包含在特征点β1周围的区域611中的两个任意像素，并且计算这两个像素之间的亮度值的差。然后，如果所计算的亮度差大于预定阈值，则特征点/量提取部件512将特征点β1的特征量设定为“1”。如果所计算的差小于或等于预定阈值，则特征点/量提取部件512将特征点β1的特征量设定为“0”。
在图6B的示例中，提取点A1和点B1作为两个任意像素，并且将特征点β1的特征量设定为“1”，这是因为这两个点之间的亮度差大于阈值。此外，提取另外的点A2和点B2作为两个任意像素，并且将特征点β1的特征量设定为“0”，这是因为这两个点之间的亮度差小于阈值。将从一个特征点周围的区域提取两个像素的处理重复多次(例如，128次)。也就是说，在一个特征点的特征量中包含与预定数量的次数对应的位数(例如，128位)。
接下来，给出对由特征点/量提取部件512进行的特征点/量提取处理的流程的描述。图7是特征点/量提取处理的流程图。尽管为了简化说明起见以下给出了对比较源CT图像401进行特征点/量提取处理的情况的描述，但是同样也适用于比较目标CT图像413。
如图7中所示，在步骤S701中，特征点/量提取部件512使用多个纵向延伸的网格线和多个横向延伸的网格线来划分比较源CT图像401。此外，特征点/量提取部件512将网格线的每个交点确定为特征点。
在步骤S702中，特征点/量提取部件512用1代替计数器值m，计数器值m用于对特征点的数量进行计数。
在步骤S703中，特征点/量提取部件512用1代替计数器值n，计数器值n用于对特征量的数量进行计数。
在步骤S704中，特征点/量提取部件512提取特征点αm周围的区域601中的点An和点Bn，并且计算所提取的点An与点Bn之间的亮度值的差(亮度差)。
然后，在步骤S705中，特征点/量提取部件512确定在步骤S704中所计算的亮度差是否大于预定值。如果在步骤S705中确定亮度差大于预定值，则处理前进至步骤S706，在步骤S706中向特征点αm的特征量加1。另一方面，如果在步骤S705中确定亮度差小于或等于预定值，则处理前进至步骤S707，在步骤S707中向特征点αm的特征量加1。
在步骤S708中，特征点/量提取部件512确定亮度差是否已被计算了预定数量次数(例如，128次)。如果在步骤S708中确定亮度差尚未被计算预定数量次数，则处理前进至步骤S709。然后，在步骤S709中，特征点/量提取部件512增大计数器值n，计数器值n用于对特征量的数量进行计数，并且此后处理返回至步骤S704。
另一方面，如果在步骤S708中确定亮度差已被计算了预定数量次数，则处理前进至步骤S710。然后，在步骤S710中，特征点/量提取部件512计算包含在针对特征点αm所计算的特征量中的为“1”的片的数量。如果计算结果表明包含在特征量中的为“1”的数量大于预定阈值，则确定特征点αm具有大亮度差，并且处理前进至步骤S711。
此处，CT图像中具有大亮度差的特征点例如被定位在患者的身体表面、左肺和右肺的边缘部分或者患者的心脏的边缘部分上。
在步骤S711中，特征点/量提取部件512将特征点αm的特征量存储在特征量DB 521中。
另一方面，如果包含在所计算的特征量中的为“1”的片的数量小于预定阈值，则特征点/量提取部件512确定特征点αm具有小亮度差，并且处理前进至步骤S712。
注意，CT图像中具有小亮度差的特征点通常被定位在患者的左肺和右肺以及心脏内。
然后，特征点/量提取部件512确定不存在特征点αm的特征量。在这种情况下，不将特征点αm的特征量存储在特征量DB 521中。
在步骤S713中，特征点/量提取部件512确定是否针对包含在比较源CT图像401中的所有特征点计算了特征量。如果在步骤S713中确定存在尚未对其计算特征量的特征点，则处理前进至步骤S714。然后，在步骤S714中，特征点/量提取部件512增大计数器值m，并且此后处理返回至步骤S703。
另一方面，如果在步骤S713中确定对包含在比较源CT图像401中的所有特征点计算了特征量，则特征点/量提取处理结束。
接下来，给出对存储在特征量DB 521中的特征量信息的描述。图8A和图8B是存储在特征量DB 521中的特征量信息的图示。
如图8A和图8B中所示，特征量信息通过被划分成关于比较源CT图像401的信息(图8A)和关于比较目标CT图像413的信息(图8B)来管理。比较源CT图像401的特征量信息800包含“特征点”和“特征量”作为信息项目。
从比较源CT图像401提取的特征点存储在图8A中所示的表中的“特征点”的列中。从比较源CT图像401提取的特征点的特征量存储在图8A中所示的表中的“特征量”的列中。
类似地，比较目标CT图像413的特征量信息810包含“特征点”和“特征量”作为信息项目。从比较目标CT图像413提取的特征点存储在图8B中所示的表中的“特征点”的列中。从比较目标CT图像413提取的特征点的特征量存储在图8B中所示的表中的“特征量”的列中。
在特征点/量提取处理中，对于被确定为不具有特征量的特征点(图8A中的特征点α3和图8B中的特征点βM-1)不存储特征量。
接下来，给出对用于通过图像显示设备120进行局部区域定位的功能中的特征点匹配部件513的功能的描述。图9是表示特征点匹配部件513的功能的概要的图示。
如图9中所示，特征点匹配部件513针对特征点的特征量计算异或以便提取出比较目标CT图像413的与比较源CT图像401的特征点对应的特征点。
具体地，特征点匹配部件513计算比较源CT图像401的特征点α1的特征量与比较目标CT图像413的各个特征点β1、β2、β3、……的特征量中的每个特征量的异或。尽管在图9中示出了仅比较源CT图像401的特征点α1的特征量与比较目标CT图像413的特征点的特征量的异或，但是特征点匹配部件513还针对包含在比较源CT图像401中的其他特征点计算异或。
此外，特征点匹配部件513基于异或的计算结果来计算特征点之间的相似度。例如，特征点匹配部件513根据包含在异或的计算结果中的“0”的数量来计算相似度。
此外，特征点匹配部件513将所计算的特征量之间的相似度中具有最大相似度的特征点确定为“对应点”。例如，如果比较源CT图像401的特征点α1的特征量与比较目标CT图像413的特征点β1的特征量之间的相似度最大，则将比较源CT图像401的特征点α1的对应点确定为比较目标CT图像413的特征点β1。
接下来，给出对由特征点匹配部件513进行的特征点匹配处理的流程的描述。图10是由特征点匹配部件513进行的特征点匹配处理的流程图。
在步骤S1001中，特征点匹配部件513用1代替计数器值ma，计数器值ma用于对比较源CT图像401的特征点进行计数。
在步骤S1002中，特征点匹配部件513用1代替计数器值mb，计数器值mb用于对比较目标CT图像413的特征点进行计数。
在步骤S1003中，特征点匹配部件513从特征量DB 521读取特征点αma的特征量。在步骤S1004中，特征点匹配部件513确定是否存在特征点αma的特征量。
如果在步骤S1004中确定不存在特征点αma的特征量，则处理前进至步骤S1012以确定不存在针对特征点αma的对应点。另一方面，如果在步骤S1004中确定存在特征点αma的特征量，则处理前进至步骤S1005。
在步骤S1005中，特征点匹配部件513从特征量DB 521读取特征点βmb的特征量。在步骤S1006中，特征点匹配部件513确定是否存在特征点βmb的特征量。
如果在步骤S1006中确定不存在特征点βmb的特征量，则处理前进至步骤S1009。另一方面，如果在步骤S1006中确定存在特征点βmb的特征量，则处理前进至步骤S1007。
在步骤S1007中，特征点匹配部件513计算特征点αma的特征量与特征点βmb的特征量的异或。在步骤S1008中，特征点匹配部件513计算特征点αma的特征量与特征点βmb的特征量之间的相似度。
然后，在步骤S1009中，特征点匹配部件513计算是否针对比较目标CT图像413的所有特征点都进行了与特征点αma的异或的计算。如果在步骤S1009中确定在比较目标CT图像413的特征点中存在尚未计算与特征点αma的异或的特征点，则处理前进至步骤S1010。在步骤S1010中，特征点匹配部件513增大计数器值ma，并且处理返回至步骤S1005。
另一方面，如果在步骤S1009中确定针对比较目标CT图像413的所有特征点都计算了与特征点αma的异或，则处理前进至步骤S1011。
在步骤S1011中，特征点匹配部件513确定所计算的相似度中的所有相似度是否都小于阈值。如果在步骤S1011中确定所计算的相似度中的所有相似度都小于阈值，则处理前进至步骤S1012。在步骤S1012中，特征点匹配部件513确定不存在与比较源CT图像401的特征点αma对应的比较目标CT图像413的特征点，并且确定不存在针对特征点αma的对应点。
另一方面，如果在步骤S1011中确定在所计算的相似度中存在大于或等于阈值的相似度，则处理前进至步骤S1013。在步骤S1013中，特征点匹配部件513将比较目标CT图像413的特征点中具有最大相似度的特征点确定为比较源CT图像401的特征点αma的对应点。
在步骤S1014中，特征点匹配部件513将在步骤S1013中确定的对应点存储在对应矢量DB 522中。如果在步骤S1012中确定不存在对应点，则特征点匹配部件513将表示不存在与特征点αma对应的对应点的信息存储在对应矢量DB 522中。
然后，在步骤S1015中，特征点匹配部件513确定是否对比较源CT图像401的所有特征点都进行了特征点匹配处理。如果在步骤S1015中确定在比较源CT图像的特征点中存在尚未对其进行特征点匹配处理的特征点，则处理前进至步骤S1016。在步骤S1016中，特征点匹配部件513增大计数器值ma，并且处理返回至步骤S1002。
另一方面，如果在步骤S1015中确定对比较源CT图像401的所有特征点都进行了特征点匹配处理，则特征点匹配处理结束。
如上所述，在根据本实施例的图像显示设备中，如果在特征点匹配处理中确定不存在特征点的特征量，则不进行用于确定对应点的处理并且确定不存在对应点。此外，如果不存在特征点βmb的特征量，则不进行对相似度的计算。
因此，图像显示设备120在执行局部区域定位程序142时可以用较小的计算量来进行特征点匹配处理。
接下来，给出对通过特征点匹配处理存储在对应矢量DB 522中的对应点信息的描述。图11是示出了存储在对应矢量DB 522中的对应点信息的图示。
如图11中所示，对应点信息1100包括“比较源CT图像的特征点”和“比较目标CT图像的特征点”作为信息项目。比较源CT图像401的特征点存储在“比较源CT图像的特征点”中。比较目标CT图像413的特征点中被确定为比较源CT图像401的特征点的对应点的特征点存储在图11的表中的“比较目标CT图像的特征点”的列中。注意，在图11的示例中，比较源CT图像401的特征点α3和αM-1表明在比较目标CT图像413中不存在对应特征点，并且确定不存在对应点。
以下给出对用于通过图像显示设备120进行局部区域定位的功能中的对应点估计部件514的功能的描述。图12是对应点估计部件514的功能的概要的图示。
对应点估计部件514从比较源CT图像401的特征点中估计与在特征点匹配部件513中确定不存在对应点的特征点对应的特征点。
在图12A和图12B中，放大的显示部分1201和1211为出于容易理解说明起见对比较源CT图像401和比较目标CT图像413的部分的放大的图示。在放大的显示部分1201中，点αx为表示患者的肺内部的点，该点可以是被特征点匹配部件513确定为不具有对应点的特征点。可以通过使用周围的特征点例如特征点α11和α12的对应矢量来估计点αx。
边界线1202和1212是表示患者的肺的内部与外部之间的边界的外围边缘。特征点α11是边界线1202上的点，边界线1212上的特征点β11(指代放大的显示部分1211)被确定为针对特征点α11的对应点。
类似地，边界线1203和1213是表示患者肺的内部与外部之间的边界的外围边缘。特征点α12是边界线1203上的点，边界线1213上的特征点β12(指代放大的显示部分1211)被确定为针对特征点α12的对应点。
此处，特征点α11和α12被指定为穿过点αx的直线1204与边界线1202和1203中的每个边界线的交点，并且相对于介于特征点α11和α12之间的点αx彼此相对。假定，点αx与特征点α11之间的距离为r1，并且点αx与特征点α12之间的距离为r2。
对应点估计部件514基于特征点α11和α12以及特征点α11与α12之间的关系来估计作为点αx的对应点的点βx。
具体地，对应点估计部件514获取穿过特征点β11和β12的直线1214，并且估计下述点βx作为点αx的对应点，所述点βx为直线1214上的点，并且在该点处距特征点β11的距离与距点β12的距离的比率为r1：r2。
接下来，给出对由对应点估计部件514进行的对应点估计处理的流程的描述。图13是由对应点估计部件514进行的对应点估计处理的流程图。
在步骤S1301中，对应点估计部件514从对应矢量DB 522中顺序地读取比较源CT图像401的被确定为不具有对应点的特征点作为要处理的特征点。
在步骤S1302中，对应点估计部件514从对应矢量DB 522中读取最近的特征点作为基准点，所述特征点在穿过要处理的点(αx)的直线上并且与介于所述特征点之间的要处理的特征点彼此相对。
在步骤S1303中，对应点估计部件514计算要处理的特征点距在步骤S1302中读取的基准点(α11，α12)的距离的比率(r1：r2)。
在步骤S1304中，对应点估计部件514搜索对应矢量DB 522并且指定比较源CT图像413的与基准点(α11，α12)中的每个基准点对应的特征点。
在步骤S1305中，对应点估计部件514获取穿过在步骤S1304中指定的与基准点(α11，α12)对应的特征点(β11，β12)的直线。此外，对应点估计部件514计算在所获取的直线上的下述点，针对该点的距与基准点对应的特征点的距离的比率等于在步骤S1303中所计算的距离的比率(r1：r2)。此外，对应点估计部件514将所计算的点(βx)确定为要处理的特征点的对应点。
在步骤S1306中，对应点估计部件514将在步骤S1305中确定的对应点作为要处理的特征点的对应点存储在对应矢量DB 522中。
在步骤S1307中，对应点估计部件514确定是否针对比较源CT图像401的其中的每个特征点被确定为不具有对应点的所有特征点都进行了对应点估计处理。如果在步骤S1307中确定存在未对其进行对应点估计处理的特征点，则处理前进至步骤S1301以对被确定为不具有对应点的特征点之后的特征点进行步骤S1301至步骤S1306的处理。
另一方面，如果在步骤S1307中确定对被确定为不具有对应点的所有特征点都进行了对应点估计处理，则对应点估计处理结束。
如上所述，根据本实施例的图像显示设备120基于距两个基准点的距离的比率来估计对应点。从而，可以用比在通过进行特征点匹配处理而计算特征点的特征量之间的相似度之后确定对应点的情况更小的计算量来估计对应点。
具体地，根据本实施例的图像显示设备120通过针对具有大亮度差的特征点例如患者的身体表面和器官(肺和心脏)进行特征点匹配处理来高准确度地确定对应点。此外，根据本实施例的图像显示设备120通过下述方式来减小计算量：基于使用高准确度地确定的对应点而获得的距离的比率来估计对应点，其中，针对具有小亮度差的特征点例如患者的肺或心脏的内部进行特征点匹配处理来高准确度地确定对应点。
从而，根据本实施例的图像显示设备120可以用较小计算量实现高准确度定位。
接下来，给出对用于通过图像显示设备120进行局部定位的功能中的定位处理部件515的功能的描述。图14是由定位处理部件515进行的定位处理的流程图。
在步骤S1401中，定位处理部件515读取存储在对应矢量DB 522中的比较源CT图像401的特征点以及比较目标CT图像413的与比较源CT图像401的特征点对应的对应点。
在步骤S1402中，定位处理部件515计算由在步骤S1401中读取的特征点和对应点限定的对应矢量。此外，定位处理部件515基于所计算的对应矢量来确定比较目标CT图像413的每个特征点的终点。此时，定位处理部件515进行校正处理使得在移动特征点之后的比较目标CT图像413变成光滑图像(也就是说，校正对应矢量)。
在步骤S1403中，定位处理部件515基于经校正的对应矢量来对比较目标CT图像413的每个像素进行定位。在步骤S1404中，定位处理部件515将所定位的比较目标CT图像显示在显示部件206上。
接下来，给出对作为用于通过图像显示设备120进行显示控制的功能的显示控制部件516的功能的描述。显示控制部件516进行显示控制用于使得并排(横向地)显示的比较源CT图像401与经受了由定位处理部件515进行的定位的比较目标CT图像413之间容易比较。
图15A至图15D是在不同时段中获得的CT图像的第一显示示例的图示。图15A和图15C示出了比较源CT图像401。图15B和图15D示出了由定位控制部件515对其已经进行了定位的比较目标CT图像1501。
此处，图15A示出了由执业医师例如医生在所显示的比较源CT图像401上将坐标(x1，y1)上的点1511指定为目标点的情况。当显示控制部件516检测到在所显示的比较源CT图像401上对点1511的指定时，显示控制部件516识别所指定的点1511的坐标。此外，显示控制部件516自动地显示在比较目标CT图像1501上的与所指定的点1511对应的点1512作为比较目标CT图像1501上的目标点。
如上所述，在由定位处理部件515进行的定位处理中进行对对应矢量的校正。因此，在比较目标CT图像1501上的与比较源CT图像401上的点1511对应的点1512与比较源CT图像401上的点1511的坐标不能为相同坐标。从而，显示控制部件516在对对应矢量的校正已反映到其上的比较目标CT图像1501上指定坐标(x1'，y1')。
图15C示出了执业医师例如医生输入放大并显示(放大显示)包含指定点1511的区域1521的指令的情况。当显示控制部件516检测到放大显示区域1521的指令时，显示控制部件516在比较目标CT图像1501上自动地放大并显示(放大显示)与区域1521对应的区域1522。
如上所述，在不同时段中获得的比较源CT图像401和比较目标CT图像1501被并排同时显示的情况下，如果将CT图像中的一个CT图像上的点指定为目标点，则显示控制部件516自动地将另一CT图像上的点指定为目标点。此外，如果针对包含CT图像中的一个CT图像上的点的区域输入放大显示的指令，则在另一CT图像上的包含对应点的区域也被自动地放大并显示(放大显示)。从而，执业医师例如医生可以容易地掌握对应于目标点的位置，并且可以容易地比较CT图像以获知病变如何改变。
接下来，给出对由显示控制部件516进行的显示处理的流程的描述。图16是由显示控制部件516进行的第一显示处理的流程图。当根据图16中所示的流程图进行处理时，假定比较源CT图像401和比较目标CT图像1501被并排显示在显示部件206上。
在步骤S1601中，显示控制部件516确定执业医师例如医生是否指定了目标点。如果在步骤S1601中确定不存在对目标点的指定，则显示控制部件516等待对目标点的指定的确定。
另一方面，如果在步骤S1601中确定指定了目标点，则处理前进至步骤S1602。在步骤S1602中，显示控制部件516确定是否对比较源CT图像进行目标点的指定。
如果在步骤S1602中确定对比较源CT图像进行了目标点的指定，则处理前进至步骤S1603。在步骤S1603中，显示控制部件516识别所指定的目标点的坐标。
在步骤S1604中，显示控制部件516改变所识别的坐标在比较源CT图像401中的显示模式。
在步骤S1605中，显示控制部件516从对应矢量DB 522中读取所识别的坐标的对应点，并且计算所读取的对应点在由定位处理部件515对所读取的对应点进行定位处理之后的坐标。此外，显示控制部件516改变所计算的坐标在比较目标CT图像1501中的显示模式。
在步骤S1606中，显示控制部件516确定执业医师例如医生是否输入了对目标点放大显示的指令。如果在步骤S1606中确定未输入对目标点放大显示的指令，则处理返回至步骤S1601。
另一方面，如果在步骤S1606中确定输入了对目标点放大显示的指令，则处理前进至步骤S1607。在步骤S1607中，显示控制部件516在比较源CT图像401中放大并显示(放大显示)包含所识别的坐标的预定区域。
然后，在步骤S1608中，显示控制部件516在比较目标CT图像1501中放大显示包含所计算的坐标的预定区域。
另一方面，如果在步骤S1602中确定对比较目标CT图像1501进行了对目标点的指定，则处理前进至步骤S1609。在步骤S1609中，显示控制部件516识别所指定的目标点的坐标。
在步骤S1610中，显示控制部件516改变所识别的坐标在比较目标CT图像1501中的显示模式。
在步骤S1611中，显示控制部件516计算所识别的坐标在定位处理部件515进行定位处理之前的坐标，并且针对所识别的坐标提取在比较源CT图像401上的对应点的坐标。此外，显示控制部件516改变所提取的坐标在比较源CT图像401中的显示模式。
然后，在步骤S1612中，显示控制部件516确定执业医师例如医生是否输入了对目标点放大显示的指令。如果在步骤S1612中确定未输入对目标点放大显示的指令，则处理返回至步骤S1601。
另一方面，如果在步骤S1612中确定输入了对目标点放大显示的指令，则处理前进至步骤S1613。在步骤S1613中，显示控制部件516在比较目标CT图像1501中放大显示包含所识别的坐标的预定区域。
然后，在步骤S1614中，显示控制部件516在比较源CT图像401中放大显示包含所提取的坐标的预定区域。
在步骤S1615中，显示控制部件516确定是否输入了结束显示处理的指令。如果在步骤S1615中确定未输入结束显示处理的指令，则处理返回至步骤S1601。另一方面，如果确定输入了结束显示处理的指令，则显示处理结束。
如上所述，根据本实施例的图像显示设备120在并排显示在不同时段中获得的CT图像作为比较源CT图像和比较目标CT图像时针对比较源CT图像进行大区域定位和局部区域定位。此外，在进行局部区域定位时，根据本实施例的图像显示设备120计算比较目标CT图像的与比较源CT图像的特征点对应的特征点作为对应点。此外，当不存在比较源CT图像的特征点的对应点时，根据本实施例的图像显示设备120通过使用在比较源CT图像的特征点附近的外围边缘中并且位于相对位置处的特征点来估计对应点。此外，根据本实施例的图像显示设备120通过使用所计算的对应点和所估计的对应点来进行对比较目标CT图像的局部区域定位。此外，如果在并排显示的比较源CT图像和比较目标CT图像中的一个CT图像中指定了目标点，则根据本实施例的图像显示设备120在另一CT图像的对应位置处自动地指定目标点。此外，如果发出在并排显示的比较源CT图像和比较目标CT图像中的一个CT图像中放大显示包含目标点的区域的指令，则根据本实施例的图像显示设备120在另一CT图像中自动地放大显示包含对应点的区域。
由此，执业医师例如医生可以容易地比较CT图像以获知在CT图像中的一个CT图像中指定的病变随时间流逝如何改变。此外，可以用更小计算量实现局部定位。
<第二实施例>
以下给出对第二实施例的描述。在第二实施例中，当并排显示比较源CT图像和比较目标CT图像时，显示控制部件516使用尚未进行大区域定位和局部区域定位的CT图像作为比较目标CT图像。
图17A至图17D是在不同时段中获得的CT图像的第二显示示例的图示。图17A和图17C示出了比较源CT图像401。图17B和图17D示出了在进行大区域定位和局部区域定位之前的比较目标CT图像401。
此处，图17A示出了由执业医师例如医生在所显示的比较源CT图像401上指定坐标(x1，y1)上的点1511作为目标点的情况。当显示控制部件516检测到在所显示的比较源CT图像401上对点1511的指定时，显示控制部件516识别所指定的点1511的坐标。此外，显示控制部件516在比较目标CT图像411上自动地显示在比较目标CT图像411上的与所指定的点1511对应的点1712。
因为比较目标CT图像411为尚未进行大区域定位和局部区域定位的图像，所以显示控制部件516通过参照对应矢量DB 522得到对应于点1511的对应点。通过参照对应矢量DB 522得到的对应点是进行了大区域定位的比较目标CT图像413上的点。从而，显示控制部件516还将所得到的对应点转换成在进行大区域定位之前的点。在比较目标CT图像411上对由此获取到的点1712(坐标(x2，y2))进行显示。也就是说，显示控制部件516在比较目标CT图像411上自动地指定坐标(x2，y2)。
图17C示出了执业医师例如医生输入放大并显示(放大显示)包含指定点1511的区域1521的指令的情况。当显示控制部件516检测到放大显示区域1521的指令时，显示控制部件516在比较目标CT图像411上自动地放大并显示(放大显示)与区域1521对应的区域1722。
如上所述，在不同时段获得的比较源CT图像401和比较目标CT图像411被横向地并排显示的情况下，如果指定了CT图像中的一个CT图像上的点，则显示控制部件516自动地指定另一CT图像上的点。此外，如果针对包含CT图像中的一个CT图像上的点的区域输入放大显示的指令，则在另一CT图像上的包含对应点的区域也被自动地放大并显示(放大显示)。从而，执业医师例如医生可以容易地掌握对应于指定点的位置，并且可以容易地比较CT图像以获知病变如何改变。
接下来，给出对由显示控制部件516进行的显示处理的流程的描述。图18A和图18B是由显示控制部件516进行的第二显示处理的流程图的一部分。在根据图18A和图18B中示出的流程图进行处理时，假定比较源CT图像401和比较目标CT图像411横向地并排显示在显示部件206上。此外，假定针对比较目标CT图像411执行了大区域定位程序141和局部区域定位程序142，并且假定对应点信息1100存储在对应矢量DB 522中。
在步骤S1801中，显示控制部件516确定执业医师例如医生是否指定了目标点。如果在步骤S1801中确定不存在对目标点的指定，则显示控制部件516等待对目标点的指定的确定。
另一方面，如果在步骤S1801中确定指定了目标点，则处理前进至步骤S1802。在步骤S1802中，显示控制部件516确定是否对比较源CT图像进行了目标点的指定。
如果在步骤S1802中确定对比较源CT图像进行了目标点的指定，则处理前进至步骤S1803。在步骤S1803中，显示控制部件516识别所指定的目标点的坐标。
在步骤S1804中，显示控制部件516改变所识别的坐标在比较源CT图像401中的显示模式。
在步骤S1805中，显示控制部件516通过参照对应矢量DB 255来提取在比较目标CT图像411上的与所识别的坐标对应的对应点的坐标。此外，显示控制部件516将所提取的对应点的坐标转换成在进行大区域定位之前的坐标。
在步骤S1806中，显示控制部件516改变在步骤S1805中转换的坐标在比较目标CT图像411上的显示模式。
然后，在步骤S1807中，显示控制部件516确定执业医师例如医生是否输入了对目标点放大显示的指令。如果在步骤S1807中确定未输入对目标点放大显示的指令，则处理返回至步骤S1810。
另一方面，如果在步骤S1807中确定输入了对目标点放大显示的指令，则处理前进至步骤S1808。在步骤S1808中，显示控制部件516在比较源CT图像401中放大并显示(放大显示)包含所识别的坐标的预定区域。
然后，在步骤S1809中，显示控制部件516在比较目标CT图像411中放大显示包含所计算的坐标的预定区域。
在步骤S1810中，显示控制部件516确定是否输入了结束显示处理的指令。如果确定未输入结束显示处理的指令，则处理返回至步骤S1801。
另一方面，如果在步骤S1802中确定对比较目标CT图像411进行了对目标点的指定，则处理前进至步骤S1901。
在步骤S1901中，显示控制部件516识别所指定的目标点的坐标。在步骤S1902中，显示控制部件516改变所识别的坐标在比较目标CT图像411中的显示模式。
在步骤S1903中，显示控制部件516计算所识别的坐标在进行大区域定位之后的坐标。此外，显示控制部件516通过参照对应矢量DB 522提取在比较源CT图像401上的对应点的坐标，该坐标对应于进行大区域定位之后的坐标。
在步骤S1904中，显示控制部件516改变在步骤S1903中提取的坐标在比较源CT图像401上的显示模式。
然后，在步骤S1905中，显示控制部件516确定执业医师例如医生是否输入了对目标点放大显示的指令。如果在步骤S1905中确定未输入对目标点放大显示的指令，则处理返回至步骤S1810。
另一方面，如果在步骤S1905中确定输入了对目标点放大显示的指令，则处理前进至步骤S1906。在步骤S1906中，显示控制部件516在指定的比较目标CT图像411中放大显示包含所识别的坐标的预定区域。
然后，在步骤S1907中，显示控制部件516在比较源CT图像401中放大显示包含所提取的坐标的预定区域。
在步骤S1810中，显示控制部件516确定是否输入了结束显示处理的指令。如果在步骤S1810中确定输入了结束显示处理的指令，则显示控制处理结束。
如上所述，根据本实施例，即使在比较源CT图像401和比较目标CT图像411横向地并排显示的情况下，当在CT图像中的一个CT图像中指定目标点时，该目标点可以在另一CT图像中的对应位置处被自动地指定。此外，当指示在CT图像中的一个CT图像中放大显示包含目标点的区域时，在另一CT图像中包含对应点的区域可以被自动地放大显示。
因此，使用根据本实施例的图像显示设备的执业医师例如医生可以容易地比较CT图像以检查在CT图像中的一个CT图像中指定的病变随时间流逝如何改变。
<第三实施例>
在上文提及的实施例中的每个实施例中，任何方法可以适用于匹配沿患者的身体轴线获得的多个CT图像中的CT图像。例如，如果沿从患者的头部朝向腿的身体轴线获得多个CT图像，则距头部相同位置(头部图像到要比较的CT图像存在相同数量的图片胶卷)的CT图像可以被匹配用于比较。
此外，尽管在上文提及的实施例中的每个实施例中图像显示设备120包括大区域定位程序141、局部区域定位程序142以及同步显示程序143，然而图像显示设备120可以包括仅同步显示程序143。在这样的情况下，可以将其他程序安装在与图像显示设备120分开设置的服务器设备中。
此外，尽管在上文提及的实施例中的每个实施例中使用x射线CT图像，然而图像显示设备可适用于显示除x射线CT图像以外的医疗CT图像例如MRI(核磁共振成像)图像的情况。
此外，尽管在上文提及的实施例中同时并排显示CT图像(横向地彼此靠近)用于比较，然而可以通过竖直地或对角地彼此靠近地布置来同时显示CT图像用于比较。
本文中详述的所有示例和传统语言意在出于教学法目的辅助读者理解发明的原理以及由发明者贡献的深化本领域的构思，并且认为其并非限制于这样具体详述的示例和情况，在说明书中的这样的示例的组织也并非涉及示出发明的优势和劣势。尽管详细地描述了本发明的一个或多个实施例，应当理解，在不偏离本发明的精神和范围的情况下可以对其做出各种变化、替代以及变型。

Claims (9)

1.一种图像显示设备，包括：
显示装置，所述显示装置被配置成同时显示第一图像和第二图像，所述第一图像和所述第二图像表示同一对象的在不同时段中获得的图像；以及
处理器，所述处理器被配置成执行用于进行以下处理的程序，所述处理包括：
根据第一显示模式将所述第一图像和所述第二图像同时显示在所述显示装置上；
确定在所述第一图像中指定的目标部分；以及
将所述第二图像的对应部分的显示模式从所述第一显示模式改变成第二显示模式，所述第二图像的对应部分对应于在所述第一图像中指定的所述目标部分；
识别在所述第一图像和所述第二图像中的每个图像中的预定部分的外围边缘；
在所述第一图像中的外围边缘上确定相对的两个点，所述两个点被定位在穿过在所述第一图像中指定的指定点的直线上；
计算在所述第一图像的外围边缘和所述第二图像的外围边缘中的每个外围边缘上的每个点的特征量；
基于在所述第一图像和所述第二图像中的外围边缘上的每个点的特征量，从所述第二图像中的外围边缘上的点中识别与所述第一图像的外围边缘上的所述两个点对应的在外围边缘上的两个点；
识别所述第一图像中的所述指定点距在所述第一图像中确定的所述两个点中的每个点的距离；以及
识别在所述第二图像中满足所识别的距离的比率的点，
其中，所述第一图像和所述第二图像中的每个是计算机断层扫描图像或核磁共振成像图像。
2.根据权利要求1所述的图像显示设备，其中，所述处理器执行的处理还包括：
对于包含在所述第一图像中的多个点与包含在所述第二图像中的多个点之间的对应关系，基于根据在包含相应的点的预定区域中的像素的亮度值所计算的特征量来确定所述对应关系。
3.根据权利要求2所述的图像显示设备，其中，所述处理器执行的处理还包括：
在通过确定所述对应关系未确定包含在所述第二图像中的对应点时，对包含在所述第二图像中的所述对应点进行估计，所述对应点对应于包含在所述第一图像中的目标点。
4.根据权利要求3所述的图像显示设备，其中，所述处理器执行的处理还包括：
基于通过确定所述对应关系所确定的对应点或者通过所述估计所估计的对应点来计算与在所述第一图像中指定的所述指定点对应的点；以及根据不同于所述第一显示模式的所述第二显示模式将所计算的点显示在所述第二图像中。
5.根据权利要求1所述的图像显示设备，其中，所述处理器执行的处理还包括：
在对所述第二图像进行定位以匹配所述第一图像之后将所述第一图像和所述第二图像彼此靠近地同时进行显示。
6.根据权利要求1所述的图像显示设备，其中，所述第二显示模式为放大显示模式，所述第二图像的所述对应部分根据所述放大显示模式被放大并且被显示在所述第二图像中。
7.一种图像显示设备，包括：
计算机，执行用于进行图像显示处理的程序，所述图像显示处理包括：
识别在不同时段中获得的两个图像中的每个图像中的预定部分的外围边缘；
在所述两个图像中的一个图像中的外围边缘上确定相对的两个点，所述两个点被定位在穿过在所述两个图像中的所述一个图像中指定的指定点的直线上；
基于在所述两个图像中的外围边缘上的每个点的特征量，从所述两个图像中的另一图像中的外围边缘上的点中识别与所述两个图像中的所述一个图像的外围边缘上的所述两个点对应的在外围边缘上的两个点；
识别在所述两个图像中的所述一个图像中的所述指定点距在所述两个图像中的所述一个图像中确定的所述两个点中的每个点的距离；以及
识别在所述两个图像中的所述另一图像中满足所识别的距离的比率的点，
其中，所述两个图像中的每个是计算机断层扫描图像或核磁共振成像图像。
8.一种同时显示第一图像和第二图像的图像显示方法，所述第一图像和所述第二图像表示同一对象的在不同时段中获得的图像，所述图像显示方法包括：
通过计算机根据第一显示模式将所述第一图像和所述第二图像同时显示在显示装置上；
通过所述计算机确定在所述第一图像中指定的目标部分；以及
通过所述计算机将所述第二图像的对应部分的显示模式从所述第一显示模式改变成第二显示模式，所述第二图像的对应部分对应于在所述第一图像中指定的所述目标部分；
通过所述计算机识别在所述第一图像和所述第二图像中的每个图像中的预定部分的外围边缘；
通过所述计算机在所述第一图像中的外围边缘上确定相对的两个点，所述两个点被定位在穿过在所述第一图像中指定的指定点的直线上；
通过所述计算机计算在所述第一图像的外围边缘和所述第二图像的外围边缘中的每个外围边缘上的每个点的特征量；
基于在所述第一图像和所述第二图像中的外围边缘上的每个点的特征量，通过所述计算机从所述第二图像中的外围边缘上的点中识别与所述第一图像的外围边缘上的所述两个点对应的在外围边缘上的两个点；
通过所述计算机识别所述第一图像中的所述指定点距在所述第一图像中确定的所述两个点中的每个点的距离；以及
通过所述计算机识别在所述第二图像中满足所识别的距离的比率的点，
其中，所述第一图像和所述第二图像中的每个是计算机断层扫描图像或核磁共振成像图像。
9.一种非暂态计算机可读存储介质，其中存储有程序，当所述程序由计算机执行时，使得所述计算机进行根据权利要求8所述的图像显示方法。
CN201510111794.3A 2014-04-28 2015-03-13 图像显示设备、方法以及存储图像显示程序的存储介质 Active CN105030265B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2014092758A JP6344039B2 (ja) 2014-04-28 2014-04-28 画像表示装置、画像表示方法及びプログラム
JP2014-092758 2014-04-28

Publications (2)

Publication Number Publication Date
CN105030265A CN105030265A (zh) 2015-11-11
CN105030265B true CN105030265B (zh) 2018-04-13

Family

ID=54335251

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510111794.3A Active CN105030265B (zh) 2014-04-28 2015-03-13 图像显示设备、方法以及存储图像显示程序的存储介质

Country Status (3)

Country Link
US (1) US9514384B2 (zh)
JP (1) JP6344039B2 (zh)
CN (1) CN105030265B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7066315B2 (ja) 2016-09-29 2022-05-13 日本製鉄株式会社 バーリング加工装置、バーリング加工方法、金属部品の製造方法、バーリング加工品及び金属部品
CN114219779A (zh) * 2016-11-24 2022-03-22 株式会社尼康 图像处理装置
KR101923183B1 (ko) * 2016-12-14 2018-11-28 삼성전자주식회사 의료 영상 표시 방법 및 의료 영상 표시 장치
WO2018198327A1 (ja) * 2017-04-28 2018-11-01 オリンパス株式会社 内視鏡診断支援システム、内視鏡診断支援プログラム及び内視鏡診断支援方法
US10467723B2 (en) * 2017-12-21 2019-11-05 Qualcomm Incorporated Tile-based check values for data content integrity in a GPU subsystem
JP7394959B2 (ja) 2020-03-05 2023-12-08 富士フイルム株式会社 医用画像処理装置、医用画像処理方法及びプログラム、医用画像表示システム

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101006446A (zh) * 2004-08-18 2007-07-25 皇家飞利浦电子股份有限公司 用于乳房造影的评估的显示系统
CN101040801A (zh) * 2006-03-23 2007-09-26 株式会社东芝 读图报告显示装置及其方法
CN103239253A (zh) * 2012-02-14 2013-08-14 株式会社东芝 医用图像诊断装置

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60159972A (ja) 1984-01-30 1985-08-21 Toshiba Corp 画像の精密位置合わせ装置
JPH08103439A (ja) 1994-10-04 1996-04-23 Konica Corp 画像の位置合わせ処理装置及び画像間処理装置
US6868172B2 (en) * 2001-10-03 2005-03-15 Eastman Kodak Company Method for registering images in a radiography application
US7296239B2 (en) * 2002-03-04 2007-11-13 Siemens Corporate Research, Inc. System GUI for identification and synchronized display of object-correspondence in CT volume image sets
JP2008000536A (ja) * 2006-06-26 2008-01-10 Fujifilm Corp 画像表示装置
JP4935693B2 (ja) * 2008-01-25 2012-05-23 コニカミノルタホールディングス株式会社 画像生成装置、プログラム、画像生成方法
JP2012095791A (ja) 2010-11-01 2012-05-24 Toshiba Corp 医用画像生成装置、医用画像生成方法及び医用画像生成プログラム

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101006446A (zh) * 2004-08-18 2007-07-25 皇家飞利浦电子股份有限公司 用于乳房造影的评估的显示系统
CN101040801A (zh) * 2006-03-23 2007-09-26 株式会社东芝 读图报告显示装置及其方法
CN103239253A (zh) * 2012-02-14 2013-08-14 株式会社东芝 医用图像诊断装置

Also Published As

Publication number Publication date
JP2015208539A (ja) 2015-11-24
US9514384B2 (en) 2016-12-06
JP6344039B2 (ja) 2018-06-20
US20150310610A1 (en) 2015-10-29
CN105030265A (zh) 2015-11-11

Similar Documents

Publication Publication Date Title
CN105030265B (zh) 图像显示设备、方法以及存储图像显示程序的存储介质
EP2961322B1 (en) Segmentation of large objects from multiple three-dimensional views
US11723722B2 (en) Methods and systems for guidance in cardiac resynchronization therapy
US7778686B2 (en) Method and apparatus for medical intervention procedure planning and location and navigation of an intervention tool
US7996060B2 (en) Apparatus, method, and computer software product for registration of images of an organ using anatomical features outside the organ
US9179890B2 (en) Model-based positioning for intracardiac echocardiography volume stitching
JP5584006B2 (ja) 投影画像生成装置、投影画像生成プログラムおよび投影画像生成方法
US10460441B2 (en) Trachea marking
EP3291735B1 (en) System and method for motion compensation in medical procedures
JP2020032198A (ja) 身体部内の超音波プローブを追跡するためのシステム
US10052032B2 (en) Stenosis therapy planning
JP2019500146A (ja) 体部の３次元モデル
US9135696B2 (en) Implant pose determination in medical imaging
JP2007319689A (ja) 心臓の表示方法および磁気共鳴装置
US11327773B2 (en) Anatomy-aware adaptation of graphical user interface
US20120224752A1 (en) Image-based diagnosis assistance apparatus, its operation method and program
EP3456248A1 (en) Hemodynamic parameters for co-registration
US9224188B2 (en) Image processing device, method and program
CN109350059B (zh) 用于肘部自动对准的组合的转向引擎和界标引擎
JP6876200B2 (ja) 動的術中セグメンテーションデータに対する静的術前計画データの位置合わせ
JP5923067B2 (ja) 診断支援装置および診断支援方法並びに診断支援プログラム
US12062198B2 (en) Method and system for multi-modality joint analysis of vascular images
Perissinotto et al. Robust temporal alignment of multimodal cardiac sequences
US20240206907A1 (en) System and Method for Device Tracking in Magnetic Resonance Imaging Guided Inerventions
EP4060363A1 (en) Automated detection of critical stations in multi-station magnetic resonance imaging

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant