CN104205172A - 图像生成装置、方法以及程序 - Google Patents

图像生成装置、方法以及程序 Download PDF

Info

Publication number
CN104205172A
CN104205172A CN201380014286.2A CN201380014286A CN104205172A CN 104205172 A CN104205172 A CN 104205172A CN 201380014286 A CN201380014286 A CN 201380014286A CN 104205172 A CN104205172 A CN 104205172A
Authority
CN
China
Prior art keywords
pixel
value
opacity curve
pixel value
chromatic graph
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.)
Pending
Application number
CN201380014286.2A
Other languages
English (en)
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.)
Fujifilm Corp
Original Assignee
Fujifilm 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 Fujifilm Corp filed Critical Fujifilm Corp
Publication of CN104205172A publication Critical patent/CN104205172A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/08Volume rendering
    • 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/46Arrangements for interfacing with the operator or the patient
    • A61B6/461Displaying means of special interest
    • A61B6/466Displaying means of special interest adapted to display 3D 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/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/504Apparatus 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 blood vessels, e.g. by angiography
    • 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/0891Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • A61B8/466Displaying means of special interest adapted to display 3D data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00Medical imaging apparatus involving image processing or analysis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4887Locating particular structures in or on the body
    • A61B5/489Blood vessels
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00Indexing scheme for image data processing or generation, in general
    • G06T2200/04Indexing scheme for image data processing or generation, in general involving 3D image data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00Indexing scheme for image generation or computer graphics
    • G06T2210/41Medical
    • 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
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/40ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Graphics (AREA)
  • Vascular Medicine (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Optics & Photonics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Generation (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

本发明在通过进行体绘制而生成伪三维图像时生成更准确地描绘了预定对象物的三维形态的图像。在三维图像中,使用定义像素值和不透明度关系的不透明度曲线进行体绘制从而生成伪三维图像时,从三维图像确定表示预定对象物的整个区域,针对所确定的整个区域设定基础的不透明度曲线,针对所确定的整个区域内的至少一部分像素的每一个,取得该像素的附近区域内的像素值的代表值,将使用该取得的代表值对所述设定的基础的不透明度曲线施加变更而成的不透明度曲线设定为在所述体绘制中应用于该像素的不透明度曲线。

Description

图像生成装置、方法以及程序
技术领域
本发明涉及在三维图像中通过进行体绘制而生成伪三维图像的图像生成装置、方法以及程序。
背景技术
以往,为了容易地掌握被摄体的三维结构等,进行生成/显示伪三维图像的处理,所述伪三维图像是使用计算机图形学的技术等将通过CT装置、MRI装置、超声波诊断装置等所取得的被摄体的三维图像数据在二维平面上立体地可视化后的图像。作为生成这种伪三维图像的方法,公知的有针对构成三维图像的各像素值(体素值)设定不透明度以及R、G、B的颜色信息,进行从观察的一侧向投影面的各像素的光线跟踪(ray-casting),从而进行可视化的体绘制法。
一般情况下,对该各像素值而设定的不透明度、颜色信息的集合(分别称作不透明度曲线、色图)对于一个三维图像而赋予一个组合。可是,该情况下,产生了不能够对具备相同信号值的不同的组织加以区别地显示的问题。对此,非专利文献1中,提出了如下方案:即利用体绘制生成以1张图像显示多个对象(例如,骨、血管、心脏、肝脏等)的图像时,能够在各自的对象范围应用不同的色图和/或不透明度曲线,从而以不同的色、不透明度表现具有相同信号值的不同的组织。
现有技术文献
专利文献
专利文献1:日本特开2011-212219号公报
非专利文献
非专利文献1:今井裕著,《なるほど!!医用3次元画像―考え方と処理法の虎の巻》,秀润社出版,2003年10月发行,第105页
发明内容
发明所要解决的课题
然而,例如血管的CT摄影中,有时采用如下方法:通过在血管中注入被称作造影剂的特殊液体来进行摄影,使得三维图像中血管的像素值成为与其他脏器不同的像素值,从而使得更加可视化且易于对该区域进行提取。可是,该情况下,由于摄影时机、造影剂的量等,即使在相同血管区域也会因场所不同而使像素值产生变动,因此存在如下问题:在生成以1张图像显示较广的范围的血管区域的图像的情况下,即使是被称为血管的一个的对象组织,若对其仅应用一个色图和/或不透明度曲线则不能够准确地显示整个血管的三维形态。
基于该像素值的变动的问题,有时也是因分体积(partial volume)的效果而引起。特别是在由几个像素(pixel)的直径所表示的较细的血管中,实际应该具有较高的像素值的血管在图像化时受到周围的像素值的影响,而作为比实际低的值被图像化。由此,例如在如图6所示的应用了与较粗的血管的像素值分布对应而设定的一个色图以及不透明度曲线而生成的图像中,较粗的血管部分清晰可见,但是较细的血管部分消失而不可见,在如图7所示的应用了使所述色图以及不透明度曲线分别与较细的血管的像素值分布对应而平衡移动而成的一个色图以及不透明度曲线而生成的图像中,较细的血管部分变得可见(参照被白线所围成的区域),但是较粗的血管部分看上去过度膨胀。
另外,在专利文献1中,提出了一种针对每1张图像而动态地对色图和/或不透明度曲线进行校正的方法,但是该方法是以基于每个图像的描绘范围小到能够忽略该范围内的像素值的变动的程度为前提并与如何决定在该1张图像的生成中所使用的一个不透明度曲线相关的方法,并没有解决生成对更广范围的对象区域进行显示的图像的情况下的上述问题。
本发明鉴于上述情形而提出,目的在于提供一种图像生成装置、方法以及程序,在通过进行体绘制而生成伪三维图像时,能够生成更加准确地描绘了预定对象物的三维形态的图像。
用于解决课题的手段
本发明的图像生成装置(第一图像生成装置)的特征在于,具备:图像生成部,在三维图像中,使用定义像素值和不透明度的关系的不透明度曲线进行体绘制,从而生成伪三维图像;区域确定部,从三维图像确定表示预定对象物的整个区域;以及不透明度曲线设定部,针对所确定的整个区域设定基础的不透明度曲线,针对所确定的整个区域内的至少一部分像素的每一个,取得该像素的附近区域内的像素值的代表值,将使用该取得的代表值对所设定的基础的不透明度曲线施加变更而成的不透明度曲线设定为在体绘制中应用于该像素的不透明度曲线。
优选的是,不透明度曲线设定部针对至少一部分像素的每一个,基于从该像素的附近区域内的像素值的代表值减去整个区域的像素值的代表值而得的减法值求得不透明度曲线的移位量,将以所求得的移位量使基础的不透明度曲线沿像素值方向移位而成的不透明度曲线设定为在体绘制中应用于该像素的不透明度曲线。
这里,在不透明度曲线是设定于在一方的轴(例如横轴)表示像素值、在另一方的轴(例如纵轴)表示不透明度的坐标系中的曲线的情况下,所谓“像素值方向”是指表示像素值的轴向。
另外,在基于上述减法值而求得移位量时,将具有与减法值的正负符号相同符号的值作为移位量而求得。此时,该移位量的绝对值可以与减法值的绝对值相同,也可以是对其乘以1以外的预定的系数而得到的值。由此,所谓“以移位量移位”是指在移位量为正的情况下沿表示像素值的轴的正方向移位,在移位量为负的情况下沿表示像素值的轴的负方向移位。
优选的是,不透明度曲线设定部针对以预定的间隔对整个区域内的像素进行采样而得的一部分像素的每一个,取得该像素的附近区域内的像素值的代表值,基于所取得的代表值和整个区域的像素值的代表值的差分来求得不透明度曲线的移位量,并将以所求得的移位量使所述基础的不透明度曲线沿像素值方向移位而成的不透明度曲线设定为在体绘制中应用于该像素的不透明度曲线,不透明度曲线设定部针对整个区域内的所采样的一部分像素以外的像素的每一个,利用位于该像素的附近的2个以上的所采样的像素各自的移位量通过插补求得该像素中的移位量的推定值,将以所求得的推定值使基础的不透明度曲线沿像素值方向移位而成的不透明度曲线设定为在体绘制中应用于该像素的不透明度曲线。
另外,附近区域内的像素值的代表值可以是附近区域内的像素值的众数、中位数或者平均值,也可以是附近区域内的像素值的整体中的相对于附近区域内的像素值的众数或者中位数处于预定的像素值宽度的范围内的像素值的平均值,也可以是附近区域内的像素值的整体中的相对于附近区域内的像素值的众数或者中位数处于预定的像素值宽度的范围内的像素值的众数或者中位数。
此时,整个区域的像素值的代表值可以是与附近区域内的像素值的代表值相同种类的值,也可以是不同的种类的值。例如,在附近区域内的像素值的代表值是附近区域内的像素值的平均值的情况下,整个区域的像素值的代表值可以是整个区域内的像素值的平均值,也可以是整个区域内的像素值的众数。
另外,优选的是,上述本发明的图像生成装置中,图像生成部在三维图像中使用定义像素值和显示色的关系的色图进行体绘制,从而生成伪三维图像,所述图像生成装置还具备色图设定部,该色图设定部针对所确定的整个区域设定基础的色图,针对所确定的整个区域内的至少一部分像素的每一个,计算该像素的附近区域内的像素值的代表值,并将使用所算出的代表值对所设定的基础的色图施加变更而成的色图设定为在体绘制中应用于该像素的色图。
本发明的图像生成装置(第二图像生成装置)的特征为,具备:图像生成部,在三维图像中,使用定义像素值和显示色的关系的色图进行体绘制,从而生成伪三维图像;区域确定部,从三维图像确定表示预定对象物的整个区域;以及色图设定部,针对所确定的整个区域设定基础的色图,针对所确定的整个区域内的至少一部分像素的每一个,计算该像素的附近区域内的像素值的代表值,将使用所算出的代表值对所设定的基础的色图施加变更而成的色图设定为在体绘制中应用于该像素的色图。
优选的是,上述第一以及第二图像生成装置中,预定对象物是血管。
本发明的第一以及第二图像处理方法中,由至少1台的计算机分别执行上述第一以及第二图像生成装置的各部所进行的处理。
本发明的第一以及第二图像处理程序是使至少1台的计算机分别执行上述第一以及第二图像显示装置的各部所进行的处理的程序。该程序记录在CD-ROM、DVD等的记录介质,或者以能够下载的状态记录于附属在服务器计算机的储存器、网络储存器并提供给用户。
发明效果
根据本发明的第一图像生成装置、方法以及程序,在三维图像中使用定义像素值和不透明度的关系的不透明度曲线进行体绘制从而生成伪三维图像时,从三维图像确定表示预定对象物的整个区域,并针对所确定的整个区域设定基础的不透明度曲线,针对所确定的整个区域内的至少一部分像素的每一个,取得该像素的附近区域内的像素值的代表值,将使用该取得的代表值对基础的不透明度曲线施加变更而成的不透明度曲线设定为在体绘制中应用于该像素的不透明度曲线,因此能够针对表示预定对象物的区域内所有像素的每一个,应用适合该描绘的不透明度曲线,并能够生成更准确地描绘了预定对象物的三维形态的图像。该效果在预定对象物是血管的情况下,进而在三维图像是在血管中注入了造影剂而进行摄影的图像的情况下,变得特别显著。
在上述本发明的图像生成装置、方法以及程序中,针对以预定的间隔对整个区域内的像素进行采样而得的一部分像素的每一个,取得该像素的附近区域内的像素值的代表值,基于所取得的代表值和整个区域的像素值的代表值的差分求得不透明度曲线的移位量,将以所求得的移位量使基础的不透明度曲线沿像素值方向移位而成的不透明度曲线设定为在体绘制中应用于该像素的不透明度曲线,针对整个区域内的所采样的一部分像素以外的像素的每一个,利用位于该像素的附近的2个以上的所采样的像素各自的移位量通过插补求得该像素中的移位量的推定值,将以该求得的推定值使基础的不透明度曲线沿像素值方向移位而成的不透明度曲线设定为在体绘制中应用于该像素的不透明度曲线。在该情况下,能够使运算处理高速化。
另外,在附近区域内的像素值的代表值是附近区域内的像素值的众数或者中位数的情况下,即使在附近区域内的像素值中存在由噪音等引起的极端的值,也能够使该值难以反映到代表值。
另外,在附近区域内的像素值的代表值是附近区域内的像素值的整体中的相对于附近区域内的像素值的众数或者中位数处于预定的像素值宽度的范围内的像素值的平均值、众数、或者中位数的情况下,能够使得附近区域内的像素值的中存在于最大侧或者最小侧的极端的值难以反映到代表值。
另外,上述本发明的图像生成装置,方法以及程序中,在针对所确定的整个区域设定基础的色图,针对所确定的整个区域内的至少一部分像素的每一个,计算该像素的附近区域内的像素值的代表值,将使用所算出的代表值对基础的色图施加变更而成的色图设定为在体绘制中应用于该像素的色图的情况下,能够针对表示预定对象物的区域内的所有像素的每一个应用适合该描绘的色图,并能够生成更准确地描绘了预定对象物的三维形态的图像。
根据本发明的第二图像生成装置、方法以及程序,在三维图像中使用定义像素值和不透明度的关系的色图进行体绘制从而生成伪三维图像时,从所述三维图像确定表示预定对象物的整个区域,针对所确定的整个区域设定基础的色图,针对所确定的整个区域内的至少一部分像素的每一个,取得该像素的附近区域内的像素值的代表值,将使用该取得的代表值对基础的色图施加变更而成的色图设定为在体绘制中应用于该像素的色图,因此能够针对表示预定对象物的区域内的所有像素的每一个应用适合该描绘的色图,并能够生成更准确地描绘了预定对象物的三维形态的图像。
附图说明
图1是第一实施方式中的图像生成装置的概略框图。
图2是表示第一实施方式中的图像生成装置的动作的流程图。
图3是第二实施方式中的图像生成装置的概略框图。
图4是表示第二实施方式中的图像生成装置的动作的流程图。
图5是表示由第二实施方式中的图像生成装置生成的图像的一例的图。
图6是表示利用以往的技术生成的图像的一例的图。
图7是表示利用以往的技术生成的图像的一例的图。
具体实施方式
以下,参照附图对本发明的第一实施方式进行说明。图1是表示图像生成装置1的概略结构的框图。另外,通过在计算机上执行读入到辅助存储装置的图像生成程序来实现图1那样的图像生成装置1的结构。此时,该图像生成程序存储在CD-ROM等存储介质或者经由互联网等网络散发,并安装于计算机。作为由计算机的CPU所执行的处理,图像处理程序对图像取得处理、图像生成处理、区域确定处理、不透明度曲线设定处理和显示控制处理进行规定,根据该规定,通过CPU执行上述各处理,使得计算机作为后述的图像取得部11、图像生成部12、区域确定部13、不透明度曲线设定部14、显示控制部16发挥功能。
另外,在图像生成装置1连接硬盘驱动器等存储装置2以及显示器等显示装置3。在存储装置2,存储有由CT、MRI、PET、SPECT、超声波图像等摄影装置对预定对象物进行摄影而得到的三维图像数据(三维图像)。三维图像是三维空间中的像素数据的集合,能够通过将由摄影装置得到的多张断层图像叠层而构成。图像取得部11取得存储于存储装置2的所述三维图像,并存储在内置于图像生成装置1或者与图像生成装置1连接的硬盘驱动器等存储装置。
区域确定部13根据由图像取得部11取得的三维图像对表示预定对象物的整个区域将那些确定。这里,所谓表示预定对象物的整个区域是指全部包含表示预定对象物的区域而不包含除此之外的区域的区域。以下,设预定对象物为血管地进行说明。区域确定部13从三维图像提取表示血管的区域(以下,也称作血管区域),将该提取的区域确定为整个区域。此时,能够利用阈值法、区域生长(Region Growing)法、水平集合(Level Set)法以及其他各种图像处理来进行血管区域的提取。
区域确定部13例如针对三维图像内的每一个局部区域计算3×3的海森(Hessian)矩阵的固有值,从而进行线状构造的探索。在包含线状构造的区域中,海森矩阵的三个固有值中的一个成为接近于0的值,而其他的两个成为相对较大的值。另外,与值接近于0的固有值对应的固有矢量表示线状构造的主轴方向。区域确定部13利用该关系,针对每个局部区域,基于海森矩阵的固有值判定线状构造相似度,针对识别为线状构造的局部区域,以其中心点为候补点进行检测。并且,基于预定的算法将利用探索所检测出的候补点连结。由此,构筑由将候补点以及候补点彼此连结的血管分支(边缘)构成的树构造。将所检测出的多个候补点的坐标信息、表示血管分支的方向的矢量信息与候补点、血管分支的识别符一起存储在存储器。接下来,针对每个所检测出的候补点,基于周边的像素值,在与血管路径垂直的截面中识别血管的轮廓(血管的外壁)。使用由图像分割(Graph-Cut)所代表的公知的分割方法进行形状的辨别。通过以上的处理,提取成为管状构造的血管区域,生成所提取的血管区域的确定所需的信息,并存储于存储器。
不透明度曲线设定部14对体绘制中应用于由区域确定部13所确定的血管区域内的各像素的不透明度曲线进行设定。不透明度曲线对像素值和不透明度的关系进行定义,并能够表示成以像素值为变量的函数。具体来说,首先,针对血管区域设定基础的不透明度曲线。例如,参照三维图像而对血管区域的像素值的分布以及其周边的像素值的分布进行研究,基于该分布求得成为对血管区域和其以外的区域进行区分的边界的像素值,设定在该像素值附近不透明度从“0”向“0”以外的值或者从“1”向“1”以外的值变化的不透明度曲线。不透明度曲线可以是不透明度从“0”向“1”阶梯状地变化的曲线,也可以是不透明度相对于像素值的增加或减少以预定的斜率变化的曲线。
接下来,不透明度曲线设定部14针对以预定的间隔对血管区域内的像素三维地进行采样而得的一部分像素的每一个,将从该像素的附近区域(例如,上下左右前后1cm范围的区域)内的像素值的平均值减去血液区域的像素值的平均值而得的减法值作为不透明度曲线的移位量而求得,并将以所决定的移位量使基础的不透明度曲线沿像素值方向移位而成的不透明度曲线设定为应用于该像素的不透明度曲线。即,在将像素值设为变量v、将基础的不透明度曲线设为OD(v)、将所决定的移动量设为m时,不透明度曲线设定部14针对所述采样的像素的每一个设定由O(v)=OD(v-m)所表示的不透明度曲线。
另外,不透明度曲线设定部14针对血管区域内的所述采样的一部分像素以外的像素的每一个,利用位于该像素附近的2个以上的所述采样的像素各自的移位量通过插补求得该像素中的移位量的推定值,并将以所决定的推定值使基础的不透明度曲线沿像素值方向移位而成的不透明度曲线设定为应用于该像素的不透明度曲线。
图像生成部12在三维图像中使用不透明度曲线设定部14所设定的不透明度曲线来进行体绘制,从而生成伪三维图像。显示控制部16在显示装置3显示由图像生成部12所生成的伪三维图像。
图2是表示图像生成装置1的动作的流程图。如图示的那样,首先,图像取得部11取得存储于存储装置2的三维图像(S1)。接下来,区域确定部13从由图像取得部11取得的三维图像提取表示血管的整个区域,将所提取的区域确定为整个区域(S2)。接下来,不透明度曲线设定部14针对血管区域设定基础的不透明度曲线(S3)。另外,不透明度曲线设定部14针对以预定的间隔对血管区域内的像素三维地进行采样而得的一部分像素的每一个,将从该像素的附近区域内的像素值的平均值减去血液区域的像素值的平均值而得的减法值作为不透明度曲线的移位量而求得,并针对所采样的一部分像素以外的像素的每一个,利用位于该像素附近的2个以上的所述采样的像素各自的移位量通过插补求得该像素的移位量的推定值(S4)。
此外,不透明度曲线设定部14针对所采样的一部分像素以外的像素的每一个,将以在步骤S3中决定的移位量使基础的不透明度曲线沿像素值方向移位而成的不透明度曲线设定为应用于该像素的不透明度曲线,并针对所采样的一部分像素以外的像素的每一个,将以在步骤S3中所决定的推定值使基础的不透明度曲线沿像素值方向移位而成的不透明度曲线设定为应用于该像素的不透明度曲线(S5)。
其后,图像生成部12使用不透明度曲线设定部14所设定的不透明度曲线针对三维图像进行体绘制,从而生成伪三维图像(S6)。并且,显示控制部16在显示装置3上显示所生成的伪三维图像(S7),并结束处理。
在本实施方式中,区域确定部13从三维图像中确定表示血管的整个区域,不透明度曲线设定部14针对所确定的整个区域设定基础的不透明度曲线,针对整个区域内的至少一部分像素的每一个,取得该像素的附近区域内的像素值的代表值,将使用所取得的代表值对基础的不透明度曲线施加变更而成的不透明度曲线设定为在体绘制中应用于该像素的不透明度曲线,图像生成部12使用所设定的不透明度曲线进行体绘制,从而生成伪三维图像。如此,针对表示血管的区域内的所有像素的每一个能够应用适合于该描绘的不透明度曲线,能够生成更准确地描绘了血管的三维形态的图像。
接下来,针对本发明的第二实施方式进行说明。图3表示本发明的第二实施方式的图像生成装置。本实施方式的图像生成装置100除了图1所示的第一实施方式的投影图像生成装置10的结构,还具备色图设定部15。其他方面与第一实施方式同样。
色图设定部15对在体绘制中由区域确定部13所确定的血管区域内的各像素中所应用的色图进行设定。色图对像素值和显示色的关系进行定义,能够表达成以像素值为变量的函数。具体来说,首先,针对血管区域设定基础的色图。接下来,色图设定部15,针对以预定的间隔对血管区域内的像素三维地进行采样而得的一部分像素的每一个,将从该像素的附近区域(例如,上下左右前后1cm范围的区域)内的像素值的平均值减去血液区域的像素值的平均值而得的减法值作为色图的移位量而求得,并将以所决定的移位量使基础的色图沿像素值方向移位而成的色图设定为应用于该像素的色图。即,在将像素值设为变量v、将基础的色图设为ClrMapD(v)、将所决定的移动量设为m时,色图设定部15针对所述采样的像素的每一个设定由ClrMap(v)=ClrMapD(v-m)所表达的色图。
另外,色图设定部15针对血管区域内的所述采样的一部分像素以外的像素的每一个,利用位于该像素的附近的2个以上的所述采样的像素各自的移位量通过插补求得该像素中的移位量的推定值,并将以所决定的推定值使基础的色图沿像素值方向移位而成的色图设定为应用于该像素的色图。
此时,在已经在不透明度曲线设定部14中求得了所述移位量以及推定值的情况下,色图设定部15能够不进行计算这些值的处理地使用由不透明度曲线设定部14所求出的移位量以及推定值来对应用于血管区域内的各像素的色图进行设定。当然,也可以是,不透明度曲线设定部14使用由色图设定部15所求得的移位量以及推定值来对应用于血管区域内的各像素的不透明度曲线进行设定。
图像生成部12在三维图像中使用由不透明度曲线设定部14所设定的不透明度曲线以及由色图设定部15所设定的色图进行体绘制,从而生成伪三维图像。图5表示由本实施方式的图像生成装置100所生成的伪三维图像的一例。图5中,与针对相同的三维图像利用以往的技术进行体绘制从而生成的图6以及图7的图像相比,能够更加准确地对整个血管的三维形态进行描绘。
图4是表示图像生成装置100的动作的流程图。如图所示,首先,图像取得部11取得在存储装置2中存储的三维图像(S11)。接下来,区域确定部13从由图像取得部11所取得的三维图像提取表示血管的整个区域,将所提取的区域确定为整个区域(S12)。接下来,不透明度曲线设定部14针对血管区域设定基础的不透明度曲线,色图设定部15针对血管区域设定基础的色图(S13)。
另外,不透明度曲线设定部14或者色图设定部15针对以预定的间隔对血管区域内的像素三维地进行采样而得的一部分像素的每一个,将从该像素的附近区域内的像素值的平均值减去血液区域的像素值的平均值而得的减法值作为移位量而求得,并针对所采样的一部分像素以外的像素的每一个,利用位于该像素的附近的2个以上的所述采样的像素各自的移位量通过插补求得该像素中的移位量的推定值(S14)。
接下来,不透明度曲线设定部14针对所采样的一部分像素的每一个,将以在步骤S14中所决定的移位量使基础的不透明度曲线沿像素值方向移位而成的不透明度曲线设定为应用于该像素的不透明度曲线,针对所采样的一部分像素以外的像素的每一个,将以在步骤S14中所决定的推定值使基础的不透明度曲线沿像素值方向移位而成的不透明度曲线设定为应用于该像素的不透明度曲线。另外,色图设定部15针对所采样的一部分像素的每一个,将以在步骤S14中所决定的移位量使基础的色图沿像素值方向移位而成的色图设定为应用于该像素的色图,针对所采样一部分像素以外的像素的每一个,将以在步骤S14中所决定的推定值使基础的色图沿像素值方向移位而成的色图设定为应用于该像素的色图(S15)。
其后,图像生成部12使用由不透明度曲线设定部14所设定的不透明度曲线以及由色图设定部15所设定的色图针对三维图像进行体绘制,从而生成伪三维图像(S16)。并且,显示控制部16将所生成的伪三维图像显示到显示装置3上(S17),结束处理。
本实施方式中,区域确定部13从三维图像确定表示血管的整个区域,不透明度曲线设定部14针对所确定的整个区域设定基础的不透明度曲线,针对整个区域内的至少一部分像素的每一个,取得该像素的附近区域内的像素值的代表值,将使用所取得的代表值对基础的不透明度曲线施加变更后而成的不透明度曲线设定为在体绘制中应用于该像素的不透明度曲线,色图设定部15针对所确定的整个区域设定基础的色图,针对整个区域内的至少一部分像素的每一个,取得该像素的附近区域内的像素值的代表值,将使用所取得的代表值对基础的色图施加变更而成的色图设定为在体绘制中应用于该像素的色图,图像生成部12使用这些设定的不透明度曲线以及色图进行体绘制,从而生成伪三维图像。由此,能够针对表示血管的区域内的所有像素的每一个应用适合于该描绘的不透明度曲线以及色图,能够生成更准确地描绘了血管的三维形态的图像。
另外,本实施方式中,对分别使基础的不透明度曲线或者基础的色图变形而求得应用于血管区域内的各像素的不透明度曲线以及色图双方的情况进行了说明,但是例如也可以是,对于不透明度曲线直接使用基础的不透明度曲线,仅对于色图使用对基础的色图进行变形而求得的色图。
另外,在上述各实施方式中,对仅针对以预定的间隔对血管区域内的像素三维地进行采样而得的一部分像素计算移位量并针对除此以外的像素的每一个通过插补而求得移位量的推定值的情况进行了说明,但是也可以是,针对血管区域内的所有像素的每一个,基于从该像素的附近区域内的像素值的代表值减去血液区域的像素值的代表值而得的减法值而求得移位量,使用所决定的移位量来设定应用于该像素的不透明度曲线和/或色图。
另外,上述各实施方式中,对针对每个像素求得移位量或者推定值的情况进行了说明,但是也可以是,对由2以上的像素构成的每一个小区域求得移位量或者推定值,并将基于该求得的移位量或者推定值而求得的不透明度曲线和/或色图应用于该小区域内的所有像素。即,可以是不透明度曲线设定部14针对构成由区域确定部13所确定的整个区域的预定大小的部分区域的每一个,取得该部分区域内的像素值的代表值,将使用该取得的代表值对基础的不透明度曲线施加变更而成的不透明度曲线设定为在体绘制中应用于该部分区域的不透明度曲线,也可以是不透明度曲线设定部14针对构成由区域确定部13所确定的整个区域的预定大小的部分区域的每一个,取得包含该部分区域的附近区域内的像素值的代表值,将使用该取得的代表值对所述设定的基础的不透明度曲线施加变更而成的不透明度曲线设定为在所述体绘制中应用于该部分区域的不透明度曲线。
同样地,可以是色图设定部15针对构成由区域确定部13所确定的整个区域的预定大小的部分区域的每一个,取得该部分区域内的像素值的代表值,将使用该取得的代表值对基础的色图施加变更而成的色图设定为在体绘制中应用于该部分区域的色图。也可以是色图设定部15针对构成由区域确定部13所确定的整个区域的预定大小的部分区域的每一个,取得包含该部分区域的附近区域内的像素值的代表值,将使用该取得的代表值对所述设定的基础的色图施加变更而成的色图设定为在所述体绘制中应用于该部分区域的色图。
另外,上述各实施方式中,对预定对象物是血管的情况进行了说明,但是预定对象物可以是肠、气管分支等具有管腔状结构的结构物,也可以是心脏、肝脏等各种脏器。
另外,上述各实施方式中,对附近区域的像素值的代表值和整个区域的像素值的代表值均是平均值的情况进行了说明,但是附近区域的像素值的代表值以及整个区域的像素值的代表值也可以分别是该区域内的像素值的众数或者中位数,也可以是该区域内的像素值的整体中的相对于该区域内的像素值的众数或者中位数处于预定的像素值宽度的范围内的像素值的平均值,也可以是该区域内的像素值的整体中的相对于该区域内的像素值的众数或者中位数处于预定的像素值宽度的范围内的像素值的众数或者中位数。此时,整个区域的像素值的代表值可以是与附近区域内的像素值的代表值相同种类的值,也可以是不同种类的值。
另外,在生成基于体绘制的图像时,在仅将整个三维图像中表示预定对象物的区域作为可视化对象区域的情况下,针对由区域确定机构13所确定的表示预定对象物的区域进行1体素程度的膨胀(Dilation)处理,在使用该膨胀处理后的区域进行体绘制的情况下,也能够利用体绘制描绘预定对象物的表面。

Claims (13)

1.一种图像生成装置,其特征在于,具备:
图像生成部,在三维图像中,使用定义像素值和不透明度的关系的不透明度曲线进行体绘制,从而生成伪三维图像;
区域确定部,从所述三维图像确定表示预定对象物的整个区域;以及
不透明度曲线设定部,针对所述确定的整个区域设定基础的不透明度曲线,针对所述确定的整个区域内的至少一部分像素的每一个,取得该像素的附近区域内的像素值的代表值,将使用该取得的代表值对所述设定的基础的不透明度曲线施加变更而成的不透明度曲线设定为在所述体绘制中应用于该像素的不透明度曲线。
2.根据权利要求1所述的图像生成装置,其特征在于,
所述不透明度曲线设定部针对所述至少一部分像素的每一个,基于从该像素的附近区域内的像素值的代表值减去所述整个区域的像素值的代表值而得的减法值求得不透明度曲线的移位量,将以该求得的移位量使所述基础的不透明度曲线沿像素值方向移位而成的不透明度曲线设定为在所述体绘制中应用于该像素的不透明度曲线。
3.根据权利要求1或2所述的图像生成装置,其特征在于,
所述不透明度曲线设定部针对以预定的间隔对所述整个区域内的像素进行采样而得的一部分像素的每一个,取得该像素的附近区域内的像素值的代表值,基于该取得的代表值和所述整个区域的像素值的代表值的差分来求得不透明度曲线的移位量,并将以该求得的移位量使所述基础的不透明度曲线沿像素值方向移位而成的不透明度曲线设定为在所述体绘制中应用于该像素的不透明度曲线,
所述不透明度曲线设定部针对所述整个区域内的所述采样的一部分像素以外的像素的每一个,利用位于该像素的附近的2个以上的所述采样的像素各自的移位量通过插补求得该像素中的移位量的推定值,将以该求得的推定值使所述基础的不透明度曲线沿像素值方向移位而成的不透明度曲线设定为在所述体绘制中应用于该像素的不透明度曲线。
4.根据权利要求1~3中任一项所述的图像生成装置,其特征在于,所述附近区域内的像素值的代表值是所述附近区域内的像素值的众数、中位数或者平均值。
5.根据权利要求1~3中任一项所述的图像生成装置,其特征在于,所述附近区域内的像素值的代表值是所述附近区域内的像素值的整体中的相对于所述附近区域内的像素值的众数或者中位数处于预定的像素值宽度的范围内的像素值的平均值。
6.根据权利要求1~3中任一项所述的图像生成装置,其特征在于,所述附近区域内的像素值的代表值是所述附近区域内的像素值的整体中的相对于所述附近区域内的像素值的众数或者中位数处于预定的像素值宽度的范围内的像素值的众数或者中位数。
7.根据权利要求1~6中任一项所述的图像生成装置,其特征在于,
所述图像生成部在所述三维图像中使用定义像素值和显示色的关系的色图进行体绘制,从而生成伪三维图像,
所述图像生成装置还具备色图设定部,该色图设定部针对所述确定的整个区域设定基础的色图,针对所述确定的整个区域内的至少一部分像素的每一个,计算该像素的附近区域内的像素值的代表值,并将使用该算出的代表值对所述设定的基础的色图施加变更而成的色图设定为在所述体绘制中应用于该像素的色图。
8.一种图像生成装置,其特征在于,具备:
图像生成部,在三维图像中,使用定义像素值和显示色的关系的色图进行体绘制,从而生成伪三维图像;
区域确定部,从所述三维图像确定表示预定对象物的整个区域;以及
色图设定部,针对所述确定的整个区域设定基础的色图,针对所述确定的整个区域内的至少一部分像素的每一个,计算该像素的附近区域内的像素值的代表值,将使用该算出的代表值对所述设定的基础的色图施加变更而成的色图设定为在所述体绘制中应用于该像素的色图。
9.根据权利要求1~8中任一项所述的图像生成装置,其特征在于,所述预定对象物是血管。
10.一种图像生成方法,在三维图像中,使用定义像素值和不透明度的关系的不透明度曲线进行体绘制,从而生成伪三维图像,其特征在于,
从所述三维图像确定表示预定对象物的整个区域,
针对所述确定的整个区域设定基础的不透明度曲线,针对所述确定的整个区域内的至少一部分像素的每一个,取得该像素的附近区域内的像素值的代表值,将使用该取得的代表值对所述设定的基础的不透明度曲线施加变更而成的不透明度曲线设定为在所述体绘制中应用于该像素的不透明度曲线。
11.一种图像生成方法,在三维图像中,使用定义像素值和显示色的关系的色图进行体绘制,从而生成伪三维图像,其特征在于,
从所述三维图像确定表示预定对象物的整个区域,
针对所述确定的整个区域设定基础的色图,针对所述确定的整个区域内的至少一部分像素的每一个,计算该像素的附近区域内的像素值的代表值,将使用该算出的代表值对所述设定的基础的色图施加变更而成的色图设定为在所述体绘制中应用于该像素的色图。
12.一种图像生成程序,其特征在于,
用于使计算机作为下述部件而发挥功能:
图像生成部,在三维图像中,使用定义像素值和不透明度的关系的不透明度曲线进行体绘制,从而生成伪三维图像;
区域确定部,从所述三维图像确定表示预定对象物的整个区域;以及
不透明度曲线设定部,针对所述确定的整个区域设定基础的不透明度曲线,并针对所述确定的整个区域内的至少一部分像素的每一个取得该像素的附近区域内的像素值的代表值,将使用该取得的代表值对所述设定的基础的不透明度曲线施加变更而成的不透明度曲线设定为在所述体绘制中应用于该像素的不透明度曲线。
13.一种图像生成程序,其特征在于,
用于使计算机作为下述部件而发挥功能:
图像生成部,在三维图像中,使用定义像素值和显示色的关系的色图进行体绘制,从而生成伪三维图像;
区域确定部,从所述三维图像确定表示预定对象物的整个区域;
色图设定部,针对所述确定的整个区域设定基础的色图,针对所述确定的整个区域内的至少一部分像素的每一个,计算该像素的附近区域内的像素值的代表值,将使用该算出的代表值对所述设定的基础的色图施加变更而成的色图设定为在所述体绘制中应用于该像素的色图。
CN201380014286.2A 2012-03-14 2013-03-13 图像生成装置、方法以及程序 Pending CN104205172A (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2012-056977 2012-03-14
JP2012056977A JP5814838B2 (ja) 2012-03-14 2012-03-14 画像生成装置、方法及びプログラム
PCT/JP2013/001634 WO2013136783A1 (ja) 2012-03-14 2013-03-13 画像生成装置、方法及びプログラム

Publications (1)

Publication Number Publication Date
CN104205172A true CN104205172A (zh) 2014-12-10

Family

ID=49160711

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201380014286.2A Pending CN104205172A (zh) 2012-03-14 2013-03-13 图像生成装置、方法以及程序

Country Status (6)

Country Link
US (1) US9563978B2 (zh)
EP (1) EP2827301B1 (zh)
JP (1) JP5814838B2 (zh)
CN (1) CN104205172A (zh)
IN (1) IN2014DN07725A (zh)
WO (1) WO2013136783A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110069310A (zh) * 2019-04-23 2019-07-30 北京小米移动软件有限公司 切换桌面壁纸的方法、装置及存储介质
CN110211216A (zh) * 2019-06-14 2019-09-06 北京理工大学 一种基于体绘制不透明度加权的三维图像空域融合方法

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6543099B2 (ja) * 2015-06-01 2019-07-10 シーフォーシー・メディカル・ソフトウェア,インコーポレーテッド 情報処理装置、情報処理方法、およびプログラム
JP6521250B2 (ja) * 2015-09-25 2019-05-29 京セラドキュメントソリューションズ株式会社 画像形成装置、色変換プログラムおよび色変換方法
JP6671946B2 (ja) * 2015-12-11 2020-03-25 キヤノン株式会社 情報取得装置、撮像装置及び情報取得方法
US10467798B2 (en) * 2016-12-19 2019-11-05 Canon Medical Systems Corporation Rendering a global illumination image from a volumetric medical imaging data set
JP7013849B2 (ja) * 2017-12-22 2022-02-01 大日本印刷株式会社 コンピュータプログラム、画像処理装置及び画像処理方法
JP7003635B2 (ja) * 2017-12-22 2022-01-20 大日本印刷株式会社 コンピュータプログラム、画像処理装置及び画像処理方法
JP2019114034A (ja) * 2017-12-22 2019-07-11 大日本印刷株式会社 コンピュータプログラム、画像処理装置、画像処理方法及びボクセルデータ
JP6436258B1 (ja) * 2018-03-27 2018-12-12 大日本印刷株式会社 コンピュータプログラム、画像処理装置及び画像処理方法
US11443476B2 (en) * 2020-05-27 2022-09-13 Canon Medical Systems Corporation Image data processing method and apparatus
JP2022090453A (ja) * 2020-12-07 2022-06-17 富士フイルムヘルスケア株式会社 磁気共鳴イメージング装置、画像処理装置および画像処理方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4421203B2 (ja) 2003-03-20 2010-02-24 株式会社東芝 管腔状構造体の解析処理装置
CN101563706B (zh) * 2006-07-31 2012-10-10 皇家飞利浦电子股份有限公司 用于为图像数据集的显像创建预设置图的方法、设备
JP2008086658A (ja) * 2006-10-04 2008-04-17 Fujifilm Corp 画像表示装置および画像表示プログラム
US7990378B2 (en) * 2007-05-07 2011-08-02 General Electric Company Methods and apparatus for volume rendering
JP5551955B2 (ja) * 2010-03-31 2014-07-16 富士フイルム株式会社 投影画像生成装置、方法、及びプログラム

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110069310A (zh) * 2019-04-23 2019-07-30 北京小米移动软件有限公司 切换桌面壁纸的方法、装置及存储介质
CN110069310B (zh) * 2019-04-23 2022-04-22 北京小米移动软件有限公司 切换桌面壁纸的方法、装置及存储介质
CN110211216A (zh) * 2019-06-14 2019-09-06 北京理工大学 一种基于体绘制不透明度加权的三维图像空域融合方法

Also Published As

Publication number Publication date
US20140375631A1 (en) 2014-12-25
EP2827301B1 (en) 2019-08-21
IN2014DN07725A (zh) 2015-06-26
JP2013191030A (ja) 2013-09-26
WO2013136783A1 (ja) 2013-09-19
EP2827301A1 (en) 2015-01-21
US9563978B2 (en) 2017-02-07
JP5814838B2 (ja) 2015-11-17
EP2827301A4 (en) 2015-11-25

Similar Documents

Publication Publication Date Title
CN104205172A (zh) 图像生成装置、方法以及程序
JP6871250B2 (ja) 医療撮像および医療撮像情報の効率的共有
JP5976431B2 (ja) 運動対象輪郭トラッキング装置及び運動対象輪郭トラッキング方法
Van Pelt et al. Exploration of 4D MRI blood flow using stylistic visualization
US20180174311A1 (en) Method and system for simultaneous scene parsing and model fusion for endoscopic and laparoscopic navigation
CN1806758A (zh) 在心脏的三维数据组中自动确定左心室位置和取向的方法
WO2017039664A1 (en) Visualization of surface-volume hybrid models in medical imaging
CN103356155A (zh) 虚拟内窥镜辅助的腔体病灶检查系统
CN103051915B (zh) 一种交互式立体视频关键帧的制作方法及制作装置
KR20020087946A (ko) 변형 가능한 3-d 대상의 이미지 시퀀스를 대상 벽 운동의표시들로 디스플레이하는 이미지 프로세싱 방법
RU2721078C2 (ru) Сегментация анатомической структуры на основе модели
CN110458834B (zh) 一种乳腺肿瘤图像处理系统、方法及装置
CN101996415B (zh) 眼球的三维建模方法
Golyanik et al. NRSfM-Flow: Recovering Non-Rigid Scene Flow from Monocular Image Sequences.
CN111105476A (zh) 一种基于Marching Cubes的CT图像三维重建方法
Sun et al. Multi-view stereo for large-scale scene reconstruction with MRF-based depth inference
Liu et al. Optimization of reconstruction of 2D medical images based on computer 3D reconstruction technology
Hentschke et al. Estimating blood flow velocity in angiographic image data
US20080144901A1 (en) Cartoon-like exaggeration of medical images to emphasize abnormalities
US10424062B2 (en) Representing an interior of a volume
Chebbi et al. 3D human heart anatomy: simulation and visualization based on MR images
US9269186B2 (en) Visualization method, and visualization apparatus
Xu et al. Animated 3d line drawings with temporal coherence
Schüle et al. Differentiable Rendering for Endoscopic Scene Reconstruction
Hong et al. Virtual angioscopy based on implicit vasculatures

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20141210