CN101203180A - 对血管进行成像的方法和设备 - Google Patents

对血管进行成像的方法和设备 Download PDF

Info

Publication number
CN101203180A
CN101203180A CNA2006800224720A CN200680022472A CN101203180A CN 101203180 A CN101203180 A CN 101203180A CN A2006800224720 A CNA2006800224720 A CN A2006800224720A CN 200680022472 A CN200680022472 A CN 200680022472A CN 101203180 A CN101203180 A CN 101203180A
Authority
CN
China
Prior art keywords
interest
area
length
elementary path
user
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
CNA2006800224720A
Other languages
English (en)
Inventor
C·P·维瑟
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
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 Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN101203180A publication Critical patent/CN101203180A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/02007Evaluating blood vessel condition, e.g. elasticity, compliance
    • 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/467Arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B6/469Arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10116X-ray image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20092Interactive image processing based on input by user
    • G06T2207/20104Interactive definition of region of interest [ROI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20112Image segmentation details
    • G06T2207/20132Image cropping
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular

Landscapes

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

Abstract

一种对血管(2)进行成像的方法,该方法包括步骤:提供管的二维视图;计算感兴趣区域(40);限定管的表征基本路径(13);在感兴趣区域(40)内部,计算垂直于基本路径(2a、2b)的基线(22);沿着基线(22),确定基线与管的侧边(2a、2b)交叉之处的转变;为许多基线寻找这种转变;考虑转变的集合来限定管的侧边(27、28);允许用户输入修改命令,以及响应于接收用户输入的修改命令,修改感兴趣区域的至少一部分并对修改的感兴趣区域重复这些计算。

Description

对血管进行成像的方法和设备
技术领域
本发明总体来说涉及对超长结构进行成像的方法。更准确地说,本发明涉及对(人或动物)体内血管进行成像的方法,而且,将在血管的环境中具体说明本发明。然而,应该注意,这个例子不应视为对本发明范围的限制,因为本发明的要点还可以被用于对其它结构进行成像。无论如何,如血管之类的超长结构的特征在于:可以沿一条路径划出中心线,该路径可以是直的或弯的,而沿着该中心线的长度来看,宽度是不变的,但是宽度也可以沿中心线的长度变化。
背景技术
成像过程典型地包括步骤:显示来源于三维体的二维图像。通过穿过所讨论的躯体而传送辐射(典型地,X-射线辐射),可以获得这个二维图像。辐射敏感接收器(典型地,感光板)接收穿过躯体的辐射,并且还接收依赖于该躯体内部结构的吸收和投射特性而形成的较暗和较亮形状的图案。在这个图案内部,可以识别血管。为了提高血管和周围组织之间的反差,通常在血管中注入造影剂。然而,对人眼来说,精确地识别血管的侧边有时是困难的。为了帮助用户,成像装置能够计算血管的侧边的装置,还能够通过在投射影象中画线而突出这些侧边。该装置还能够例如识别狭窄,并以毫米或与原状管直径的百分比或两者一起来计算该狭窄的尺寸。
用于计算血管的轮廓的方法在实践中和文献中都是已知的;举例来说,参考US-6 829 379。如图1A-1F所示,这些方法通常按下述方式运行。
首先,如图1A所示,获得关于血管2及其邻近结构3的图像,并且,在计算机系统的屏幕上显示这个图像1。典型地,图像1是作为使用X-射线辐射的投影图像而被获得的,但是,它也可能从驻留在计算机系统的存储器中的数据集处获得。
典型地,通过限定第一段端点11和第二段端点12,使用户可以确定感兴趣段,两个端点11和12都位于管2的图像内。为限定该段的端点,计算机系统可以具有鼠标之类的输入设备,其允许用户在屏幕上移动指针和在期望的位置上进行″点击″。由于在屏幕上移动指针和在期望的位置上进行″点击″是为人所熟知的,因而这里不再对其赘述。
其次,在端点11和12之间计算表征基本路径13,该基本路径跟随管2的形状。可以通过计算机使用最小成本算法来计算这个基本路径13,但是,也可以由用户确定该基本路径,方式是利用作为画图工具的指针图标在管2的图像内画直线或曲线。应该注意,这个基本路径并不一定等于血管的中心线,但是,典型地,该基本路径是中心线的良好近似。
如图1B中的较大比例所示,在基线13上定义基点21。为这个基点21计算基线22,其贯穿基点21并与基本路径13垂直。基线22具有有限的长度:在23和24处标明基线22的端点。
图1B示出基线22分别在边缘点25和26穿越管2的侧边2a和2b。通过扫描基线22以及检测从第一亮度到第二亮度的转变,计算机程序能够确定这些边缘点25和26。在管与其周围之间反差良好的情况下,侧边2a和2b将对应于从暗到亮的转变,反之亦然。
对沿着基本路径13长度的多个基点重复上述步骤。典型地,这些基点彼此间的位置相当接近。因此,获得第一组边缘点25,其表明了第一侧边2a的位置,并且,获得第二组边缘点26,其表明了第二侧边2b的位置。如图1C所示,在每组边缘点中,分别通过边缘线27和28来连接邻近的边缘点。边缘线27和28可以被定义为从边缘点到边缘点的一系列直线段,但是,它也可以是平滑的线。这两个边缘线表示计算的血管边缘,并显示血管2的轮廓。监测者可以观察血管的轮廓,并且可以使用这个信息来做出诊断。图1D示出一个例子,其中的血管具有狭窄,即,该管的流路变窄的地方。典型地,计算机能够识别该狭窄,并能够计算狭窄31处的流路直径D2,流路直径D2或者被表示为以毫米计的绝对值,或者被表示为相对于原状血管直径D1的比率。由于已经知道用于识别和计算狭窄的算法,而本发明的目的也不是旨在提供用于识别和计算狭窄的改进算法,因此在此不必再对这种算法赘述。
在确定基线22上的边缘点25和26的步骤中,遇到一个问题。在较暗和较亮图像图案之间的反差的基础上,计算边缘点,即,基线与管侧边的交叉点。根据血管在躯体内的位置和观测方向,躯体的其它部分(如,骨骼)可能会产生干扰,从而使得计算机算法采用错误的转变作为边缘点的位置。图1E中示出了这个问题,其中,图像1包括接近于管2的图像的某些躯体部分5的图像。管2与干扰躯体部分5之间的反差相对较小,并且,这个转变可能被计算机算法所忽略,而计算机算法采用干扰躯体部分5和其周围之间相对较大的反差作为管的侧边的指示。在此情况下,计算机算法将假定错误的边缘线29,如图1E所示。血管2的这个错误的表示可以导致错误的诊断,例如,计算错误的狭窄位置和/或错误的狭窄大小。
在图1E的例子中,如果有经验的用户留心且注意到算法采用了躯体部分5的边缘而不是管2的实际边缘,该有经验的用户可以在视觉上辨别出所述错误。然而,即使由用户辨别出错误,现有技术水平也没有给用户提供改正错误的工具。
本发明的目的是克服这个问题。
发明内容
如上所述,基线22具有端点23和24。只有在寻找位于所述端点23和24之间的边缘点25和26时才扫描该基线。所有基线的端点共同限定了感兴趣区域,基本路径的形状和基线的长度限定该感兴趣区域的轮廓。在现有技术水平下,所有的基线具有同样的长度,而且这个长度是固定的。根据本发明,提供一种工具,其允许用户改变感兴趣区域。用户可以在视觉上分析被计算和被呈现在屏幕上的管轮廓,他可以看出哪部分轮廓看起来是正确的,而哪部分轮廓看起来是错误的,并且,他可以改变感兴趣区域以促使计算机算法重新计算轮廓。
根据本发明,感兴趣区域也显示在显示屏上。这总是可以做到的,但是,优选地,只有在接收相应的用户请求后才作此显示。可以通过按压一个键或多个键的组合来提供用户请求,但也可能的是,当计算的管的轮廓显示在显示屏上时,计算机询问用户是否对该结果满意;如果该用户的反应表明该用户不满意,那么可以自动显示感兴趣区域的轮廓。
附图说明
在下面的描述中,将参照附图进一步说明本发明的这些及其它的方面、特征和优点,同样的附图标记代表同样的或相似的部分。
图1A-E示意性地示出在2D图像中寻找血管的侧壁的过程;
图2是示意性地示出成像系统的框图;
图3示意性地示出正显示在显示屏上的感兴趣区域;
图4A-B示出感兴趣区域的总体移动;
图5A-B示出修改感兴趣区域的基本路径的形状;
图6A-D示出修改感兴趣区域的宽度;
图7A-C示出按比例修改感兴趣区域的宽度;
图7D示出正由计算机产生的操控点;
图8A-E示出修改一段感兴趣区域的宽度的几个变化;
图9A示出增加固定长度增量的修改;
图9B示出乘以恒定因数的修改;
图10A-C以3D形式示出修改感兴趣区域的宽度的几个变化;
图11A-11B示出本发明在防止计算误差方面的效果。
具体实施方式
图2是示意性地示出根据本发明成像系统30的框图。成像系统30包括控制设备31(典型地,适当编程的计算机)、显示屏32、数据源33和用户输入设备34。数据源33可以是含有数据集的存储器,但也可以包括X-射线辐射装置等等。用户输入设备34可以包括键盘,但优选地,包括鼠标。在优选的实施例中,成像系统30具有图形用户界面,其包括显示在屏32上的指针50,可以通过用户的命令移动该指针。此外,用户可以输入选择命令,例如,通过点击鼠标键,来选择指针当时所指向的对象。进而,移动指针将移动所选择的对象(拖动)。由于点击和拖动的概念是为人所熟知的,这里不必再赘述。
图3示意性地示出正显示在显示屏32上的感兴趣区域40。图3示出:感兴趣区域40是由基本路径13、基线22的集合以及基线22的端点23和24所限定的。位于基本路径13的一侧的所有端点23的集合共同限定了感兴趣区域40的第一侧边41。同样地,位于基本路径13的另一侧的所有端点24的集合共同限定了感兴趣区域40的第二侧边42。作为参照,图3中也显示了管2的图像。
成像系统30响应用于改变感兴趣区域40的用户输入命令。如在下面的例子中所假定的那样,可以经由图形界面来提供用户输入命令。
对用户修改来说,可能有几个变化。
在一个实施例中,可以整体上移动感兴趣区域40,包括基本路径13。这在图4A和4B中示出。图4A示出初始的感兴趣区域40,还示出可以放置指针图标50的可能位置,以选取感兴趣区域40:如果这个位置位于基本路径13上则非常直观。图4B示出基本路径13以及相对于管图像2进行了移动的感兴趣区域40。
在第二实施例中,可以修改基本路径13的形状。这在图5A和5B中示出。图5A示出具有初始的基本路径13的初始的感兴趣区域40。通过点击基本路径13,用户放置两个固定图标(anchor icon)52和53,并在这两个固定图标52和53之间放置操作图标(handle icon)。利用指针图标50,用户可以点击和拖动操作图标51。作为响应,控制设备31计算新的基本路径13′。在超出固定图标52和53的区域,基本路径13′与初始的基本路径13相同。在固定图标52和53之间的区域,计算作为经过固定图标52和53以及操作图标51的平滑线的新的基本路径13′,其在固定图标52和53与初始的基本路径13平滑地连接。图5B示出新的基本路径13′和相应的感兴趣区域40;其中夸大了基本路径的修改。
在第三实施例中,可以在保持基本路径13的形状的同时来修改感兴趣区域40的宽度。修改感兴趣区域的宽度包括修改基线的长度,而保持基线的方向与基本路径垂直。
通常,所有的基线具有同样的长度,而且每个基线的位置相对于基本路径13是对称的,即,各个基线22的中点位于基本路径13上。在第三实施例的第一变化中,在增加或减小基线长度的同时保持这些特征。用户可以通过按压适当的键来进行输入,但是,也可以使用鼠标来进行输入。图6A示出具有初始的基本路径13的初始的感兴趣区域40。通过点击感兴趣区域40的侧边42,用户在这个侧边42上放置操作图标61。在62处表明穿过这个操作图标61的相应基线。利用指针图标50,用户可以沿着这个基线62点击并拖动操作图标61,从而移动这个基线62的端点。控制设备31计算基线62的修改的长度L62,并更改所有的基线使其具有同样的长度L62,同时保持中点的位置仍与基本路径重合。图6B中示出了可能结果的例子。
在第三实施例的第二变化中,可以只移动具有操作图标61的侧边42,而另一侧边41保持不变。图6C中示出了可能结果的例子。除了需要另外按压一个例如控制按钮或转换按钮的命令按钮之外,用于实现这个第二变化的用户操作可以等同于用于实现第一变化的用户操作。可能的是,使用该命令按钮的用户操作导致第一变化,而不使用该命令按钮的用户操作导致第二变化,或反之亦然。也可能的是,使用第一命令按钮的用户操作导致第一变化,使用第二命令按钮的用户操作导致第二变化,而忽略不使用命令键的用户操作。
操作图标61是否沿着相应的基线62移动并不重要。可替换地,该操作图标可能移动至任何期望的位置。作为响应,控制设备31计算从这个位置到基本路径13的距离L61,并使用这个距离,或者对称地(第一变化),或者只在基本路径的一侧(第二变化),来修改所有的基线。图6D中示出了第二变化的可能结果的例子。
在第三变化中,沿着基本路径的长度基线长度可能并不是不变的,而是沿着基本路径的长度增大或减小。这在图7A-C示出。图7A示出具有初始的基本路径13的初始的感兴趣区域40。通过点击感兴趣区域40的角43,用户在这个角43上放置操作图标71。在72处表明经过这个操作图标71的相应基线,其对应于感兴趣区域的端点12。利用指针图标50,用户可以沿着这个基线72点击并拖动操作图标71,从而移动这个基线72的端点。控制设备31计算基线72的修改长度L72,并改变所有基线的长度,这样,经过相反端点11的基线(末端边缘)的长度保持在其初始长度L0,而所有中间的基线的长度按与相应的基点到所述相反端点11的距离成比例地改变,同时保持中点的位置仍与基本路径重合。这在图7B中示出,其中,为了简化起见,基本路径13显示为直线,因此,所有的基线都假定为平行。
在第四变化中,可以只移动具有操作图标71的侧边42,而另一侧边41保持不变。图7C示出可能的结果的例子,其可与图7B进行比较。除了需要另外按压一个例如控制按钮或转换按钮的命令按钮之外,用于实现这个第四变化的用户操作可以等同于用于实现第三变化的用户操作。
如上所述,用户可以通过在期望的位置进行点击来放置操作图标。或者,如图7D所示,也可以由用户给出一般的修改命令,并且作为响应控制设备31在侧边41、42上和/或感兴趣区域40的角中放置操作图标61和71。如上所述,通过点击和拖动这些操作图标中所选择的那一个,用户可以进行上述相应的修改。
在上述变化中,沿基本路径的长度方向基线长度或者保持不变,或者从最小长度L0(端点12处)到最大长度L72(相对的端点11处)线性地改变。在第五变化中,在两个端点11和12之间的某个基点处,可以给基线长度的一个极值(最大值/最小值)。这在图8A-D示出。图8A示出具有初始的基本路径13的初始的感兴趣区域40。通过点击侧边42,用户放置两个固定图标82和83,并在这两个固定图标82和83之间放置操作图标81。利用指针图标50,用户可以点击和拖动操作图标81,或者是沿着相应的基线84(见第一变化),或者是以任意的方向。作为响应,控制设备31计算相应的基线84的修改长度,或是由操作图标81的新位置所限定的新基线84′的修改长度。控制设备31还计算在两个固定图标82和83之间的区域中的所有基线的修改长度,同时保持这个区域外的所有基线的长度。
在两个固定图标82和83之间的区域中的修改的基线可以都获得同样的长度,其等于所述相应基线84或84′的长度L84。图8B中示出了可能结果的例子。该图也示出:相对于基本路径13,修改的基线可以保持对称。
在两个固定图标82和83之间的区域中的修改的基线长度也可以从初始的长度L0(在固定图标82和83处)到极值L84按比例增大或减少。图8C中示出了可能结果的例子。该图也示出:可以只修改位于基本路径13的一侧的那一半基线的长度,而保持另一半基线的长度。
在两个固定图标82和83之间的区域中的修改的基线长度也可以从初始的长度L0(在固定图标82和83处)到极值L84按比例增大或减少,这样,感兴趣区域40的侧边42就不会有锐利的边缘。图8D中示出了可能结果的例子。
在上述所有例子中,用户通过按压或不按压某个键可以给出用于单侧修改或对称修改的命令。
应该注意,上述例子将初始的感兴趣区域40作为出发点,其中所有的基线都具有同样的长度。然而,也可以执行两个或更多的连续的修改步骤:在此情况下,用于后续修改步骤的出发点应是先前修改步骤的结果,而在这样情况下,起始位置可能具有某种不对称性;例如,感兴趣区域的起始轮廓可能如图8B中所示。在下一个修改步骤中,通过增加基线长度的绝对值或将基线增加一定的百分比来进行修改。图9A中示出了可能结果的例子。利用如上参照图6C所述的操作图标61,通过向所有的基线增加固定的长度增量Δ91,来单侧地增加如图8B中所示的感兴趣区域的宽度(即,[新的长度]减[旧的长度]所得的差对所有的基线都是相等的)。可替换地,如图9B所示,通过将所有的基线增加同样的百分比,来单侧地增加如图8B中所示的感兴趣区域的宽度(即,[新的长度]除以[旧的长度]所得的比率对所有的基线都是相等的)。
同样地,可以通过增加绝对长度增量ΔL75来实现图7B中所示的修改,其中,增量大小与沿着基本路径13所测量的离端点11的距离成比例,或者可以通过将每个基线长度乘以相乘因子R76=L76′/L76来实现图7B中所示的修改,该相乘因子与沿着基本路径13所测量的离端点11的距离成比例。
参照图7B和7C,已经描述了多个修改,其中指针50被用于选取放置在感兴趣区域40的角43的操作图标71,以将感兴趣区域40的整个侧边42作为围绕相对角44的转轴,为了保持感兴趣区域40的长度,相应地增加侧边42的长度。同样地,参照图8A,用户可以使用指针50来选取固定图标82或83。将所选择的固定点83远离或朝向基本路径13移动将与移动角图标71具有相似的效果,即,将位于固定图标82和固定图标83之间的部分侧边42作为围绕相对的固定图标82的转轴,如图8E中所示出的单侧修改。
应该注意,图1A的图像1是三维实体的二维可视化。可以以不同的方式获得这个二维的可视化。
首先,可以作为投影图像而获得该二维图像,例如,X-射线照片。促使成像辐射经过处于监视之下的躯体,从辐射源到辐射敏感表面。位于该源与接收器之间的所有的躯体部分都有助于该图像,即,该图像具有″深度″。
其次,可以作为横截面而获得该二维图像,例如,CT扫描。只有位于该横截面内的那些躯体部分才有助于图像,即,该图像没有″深度″。
第三,也可以根据自三维数据集来计算图像,该数据集可以通过例如MRI扫描而获得。在这样的三维数据集中,感兴趣的血管实际上是三维对象,而感兴趣区域实际上具有三维形状,类似于绕着基本路径的弯曲的圆筒(管子)。在此情况下,实际上可以在三维空间中执行用于在感兴趣区域内部寻找血管的边缘(轮廓)的计算,而结果则作为二维视图呈现在二维图形显示界面(显示屏)上。同时,如上所述,可以利用二维的图形界面工具(指针、操作图标;点击、拖动)来给出用户命令。在此情况下,用户操作的结果也可是三维的:即使是感兴趣区域40仅仅作为二维轮廓而显示在二维图像1中时(这使得修改看起来只影响二维图像1的平面中的感兴趣区域40),感兴趣区域40实际上在三维空间也受影响。
图10A根据垂直于基本路径13的平面示出感兴趣区域40的示意性的横截面;在本节中,假定感兴趣区域40的轮廓是圆的。显示了直角坐标系X、Y和R、;第三坐标Z垂直于所画平面,即,平行于基本路径13。假定图1A中的图像1位于XZ平面,以使得图1A中的图像1的观测方向对应于Y方向。参照图8A的例子,图10A中示出操作图标81。
当用户使用指针50来移动操作图标81远离基本路径13时,对应于同一基点21的所有基线可以被扩大同样的量(或者按照绝对值,或者按照百分比),而与它们的坐标无关,如图10B所示(也见图8B)。
该增加也可能在Y方向为零,而与从X方向到Y方向的坐标成比例,如图10C所示。该增加可以线性地依赖于坐标,或根据正弦曲线,或任何其它适当的曲线。
因此,本发明提供高度直观的、易于使用的工具,以供修改感兴趣区域40。因此,用户可以清除某些错误。参照图11A-B对此进行解释。图11A示意性地示出图1D-E的管2和躯体部分5,并且,还示出基本路径13以及相对较宽的具有侧边41和42的感兴趣区域40。管2和躯体部分5两者都位于感兴趣区域40内部。在此情况下,算法可能提供错误的计算结果,如图1E中所示。
图11B示出与图11A一样的图像,感兴趣区域40的宽度现在被单侧地减小了。虽然右侧边42并没有精确地沿着管2的侧边,而且,虽然感兴趣区域40仍包括躯体部分5的一部分,但是,算法却不再考虑从躯体5到周围的转变,这迫使算法使用从躯体部分5到管2的转变,从而获得正确的计算结果,如图1D所示。
因此,本发明的成功之处在于提供用于对血管进行成像的方法。该方法包括的步骤为:
提供该管的二维视图;
限定该管的表征基本路径13;
计算感兴趣区域40作为与基本路径13对齐的带(strip);
在感兴趣区域40内部,计算垂直于基本路径13的基线22;
沿着基线22,确定基线与管的侧边2a和2b交叉之处的转变;
对许多基线确定这种转变;
考虑转变的集合限定了管的侧边27和28;
允许用户输入修改命令,以及响应于接收用户输入的修改命令,修改感兴趣区域的至少一部分和重复对修改的感兴趣区域的计算。
对所属技术领域的专业人员来说可以理解:本发明不局限于上面讨论的示例性的实施例,而且,在权利要求所限定的本发明的保护范围内,还可以做出若干变更和修改。
如上所述,已经参照框图解释了本发明,该框图示出根据本发明的设备的功能块。可以理解,一个或多个这些功能块可以实现为硬件,其中,这种功能块的功能是由单独的硬件元件来执行的,但是,一个或多个这些功能块也可以实现为软件,这样,这种功能块的功能是由计算机程序或可编程设备例如微处理器、微控制器、数字信号处理器等等的一个或多个程序行来执行的。

Claims (18)

1.一种对超长结构(2)进行成像的方法,该方法包括步骤:
a)提供结构(2)和其周围的二维视图;
b)限定该结构的一段;
c)限定该结构段的端点(11、12)之间的表征基本路径(13);
d)限定感兴趣区域(40)为与基本路径(13)对齐的带;
e)在感兴趣区域(40)内,计算该结构段的侧边(2a、2b);
f)在该结构的二维视图上投影所计算的边缘的表示(27、28);
其中,步骤(e)包括步骤::
e1)在基本路径上限定至少一个基点(21);
e2)计算经过基点(21)且垂直于基本路径(13)的基线(22),该基线具有预定的长度和对应于感兴趣区域(40)的侧边(41、42)的两个相对的端点(23、24);
e3)沿着基线(22),在基本路径(13)的两侧,扫描二维视图并确定表示基线与该结构的侧边(2a、2b)交叉之处的交叉点(25、26)的转变;
其中,步骤(f)包括限定连接在多个基点处所获得的交叉点(25、26)的线(27、28)的步骤;
该方法的特征在于以下步骤:
g)在该结构的二维视图上投影感兴趣区域(40)的表示;
h)接收用户输入的修改命令,以供修改感兴趣区域;
i)响应于接收用户输入的修改命令,修改感兴趣区域的至少一部分;
j)就所修改的感兴趣区域,重复步骤(e)和(f)。
2.如权利要求1所述的方法,其中,步骤(i)包括在总体上移动感兴趣区域的步骤。
3.如权利要求1所述的方法,其中,步骤(i)包括修改感兴趣区域(40)的至少一部分的形状的步骤。
4.如权利要求3所述的方法,其中,步骤(i)包括修改基本路径(13)的至少一部分的形状的步骤。
5.如权利要求3所述的方法,其中,步骤(i)包括修改基本路径(13)的至少一部分的宽度的步骤。
6.如权利要求5所述的方法,其中,在基本路径(13)的所述部分内,修改基线(22)的长度,并且在基线(22)的端点(23、24)的已修改位置的基础上计算被修改的感兴趣区域(40)的被修改的侧边(41′、42′)。
7.如权利要求6所述的方法,其中,长度修改是以加上某一增量/减量(ΔL75;Δ91)而计算的,或者是以乘以某一因数(L76′/L76;L92′/L92)而计算的。
8.如权利要求7所述的方法,其中,在感兴趣区域的所述部分内,所有基线的长度修改都是相等的。
9.如权利要求7所述的方法,其中,在感兴趣区域的所述部分内,所有基线的长度修改相对于基本路径都是对称的。
10.如权利要求7所述的方法,其中,在感兴趣区域的所述部分内,所有基线的长度修改沿着基本路径的长度连续地增加或减小。
11.如权利要求7所述的方法,其中,在感兴趣区域的所述部分内,在所述部分内的位置处所有基线的长度修改都具有极值(最大值;最小值),而在所述部分的末端处基线的长度修改为零。
12.如权利要求1所述的方法,其中,步骤(h)包括接收图形输入命令的步骤。
13.如权利要求12所述的方法,其中,步骤(h)包括以下步骤:在感兴趣区域上投影至少一个操作图标(51;61;71),投影至少一个指针图标(50),接收用于移动该指针图标的用户输入命令,接收用于在视觉上将该指针图标附着在该操作图标上的用户输入命令(″点击″),和接收用于同指针图标(50)一起移动操作图标(51;61;71)的用户输入命令(″拖动″)。
14.如权利要求13所述的方法,其中,操作图标(51;61;71)由控制设备(31)限定。
15.如权利要求13所述的方法,其中,操作图标(51;61;71)由用户放置。
16.如权利要求13所述的方法,进一步包括步骤:投影两个固定点(82、83)以限定感兴趣区域的所述部分,以及在该两个固定点(82、83)之间投影操作图标(81)。
17.一种用于对血管(2)中的狭窄(31)进行计算的方法,该方法包括步骤:
使用权利要求1-16中的任一方法来对血管(2)成像;
基于计算的侧边(27、28)来计算血管(2)的平均直径(D1);
基于计算的侧边(27、28)来确定血管(2)的较窄部分(31);
基于计算的侧边(27、28)来计算较窄部分(31)的通道的直径(D2);
提供表示狭窄(31)的直径(D2)的输出信号,或者表示狭窄(31)的直径(D2)与血管(2)的平均直径(D1)之间的比率(D2/D1)的输出信号。
18.一种成像装置(30),被设计为执行根据权利要求1-17中的任一方法。
CNA2006800224720A 2005-06-21 2006-06-20 对血管进行成像的方法和设备 Pending CN101203180A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP05105496 2005-06-21
EP05105496.3 2005-06-21

Publications (1)

Publication Number Publication Date
CN101203180A true CN101203180A (zh) 2008-06-18

Family

ID=37055972

Family Applications (1)

Application Number Title Priority Date Filing Date
CNA2006800224720A Pending CN101203180A (zh) 2005-06-21 2006-06-20 对血管进行成像的方法和设备

Country Status (5)

Country Link
US (1) US20100067753A1 (zh)
EP (1) EP1903943A1 (zh)
JP (1) JP2008543482A (zh)
CN (1) CN101203180A (zh)
WO (1) WO2006137016A1 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103109308A (zh) * 2010-09-20 2013-05-15 皇家飞利浦电子股份有限公司 管状结构的管腔特性的量化
CN103505288A (zh) * 2012-06-29 2014-01-15 通用电气公司 超声成像方法和超声成像设备
CN103109308B (zh) * 2010-09-20 2016-11-30 皇家飞利浦电子股份有限公司 管状结构的管腔特性的量化

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9788729B2 (en) * 2008-11-21 2017-10-17 Toshiba Medical Systems Corporation Image processing apparatus and image processing method
US8290231B2 (en) * 2009-01-23 2012-10-16 Naveen Garg Method and apparatus for providing measurement data of an anomaly in a medical image
US8659603B2 (en) * 2009-03-31 2014-02-25 General Electric Company System and method for center point trajectory mapping
US10799209B2 (en) 2012-12-26 2020-10-13 Philips Image Guided Therapy Corporation Measurement navigation in a multi-modality medical imaging system
EP3229672B1 (en) * 2014-12-08 2021-11-17 Koninklijke Philips N.V. Bedside interface for percutaneous coronary intervention planning
EP3503026A1 (en) * 2017-12-20 2019-06-26 Koninklijke Philips N.V. Device, system and method for interacting with vessel images
CN109993743B (zh) * 2019-04-09 2023-06-06 飞依诺科技(苏州)有限公司 血管图像处理方法、装置、设备及存储介质

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4709333A (en) * 1986-01-03 1987-11-24 General Electric Company Method and apparatus for imaging in the presence of multiple high density objects
US5383463A (en) * 1993-08-02 1995-01-24 Friedman; Zvi Mapping of flow parameters
JPH0765154A (ja) * 1993-08-31 1995-03-10 Toshiba Corp 血管像の定量解析装置及びその定量解析方法
US5704028A (en) * 1994-12-21 1997-12-30 Micrografx, Inc. Graphics systems and method having data fields and shape placement control
US6047080A (en) * 1996-06-19 2000-04-04 Arch Development Corporation Method and apparatus for three-dimensional reconstruction of coronary vessels from angiographic images
US6331116B1 (en) * 1996-09-16 2001-12-18 The Research Foundation Of State University Of New York System and method for performing a three-dimensional virtual segmentation and examination
US5971767A (en) * 1996-09-16 1999-10-26 The Research Foundation Of State University Of New York System and method for performing a three-dimensional virtual examination
US7194117B2 (en) * 1999-06-29 2007-03-20 The Research Foundation Of State University Of New York System and method for performing a three-dimensional virtual examination of objects, such as internal organs
AU5101898A (en) * 1996-11-08 1998-05-29 E.I. Du Pont De Nemours And Company Method of quantifying visual uniformity of planar objects
US6077226A (en) * 1999-03-30 2000-06-20 General Electric Company Method and apparatus for positioning region of interest in image
US6546271B1 (en) * 1999-10-01 2003-04-08 Bioscience, Inc. Vascular reconstruction
US6829379B1 (en) * 2000-11-27 2004-12-07 Ge Medical Systems Global Technology Company, Llc Methods and apparatus to assist and facilitate vessel analysis
US6643533B2 (en) * 2000-11-28 2003-11-04 Ge Medical Systems Global Technology Company, Llc Method and apparatus for displaying images of tubular structures
US7570802B2 (en) * 2001-12-27 2009-08-04 The United States Of America As Represented By The Secretary Of The Department Of Health And Human Services Automated centerline detection algorithm for colon-like 3D surfaces
JP4060615B2 (ja) * 2002-03-05 2008-03-12 株式会社東芝 画像処理装置及び超音波診断装置
US7471814B2 (en) * 2002-11-27 2008-12-30 The Board Of Trustees Of The Leland Stanford Junior University Curved-slab maximum intensity projections
US7170517B2 (en) * 2002-11-27 2007-01-30 The Board Of Trustees Of The Leland Stanford Junior University Curved-slab maximum intensity projections
JP4564233B2 (ja) * 2003-01-31 2010-10-20 株式会社東芝 画像処理装置、画像データ処理方法、及びプログラム
JP4421203B2 (ja) * 2003-03-20 2010-02-24 株式会社東芝 管腔状構造体の解析処理装置
CA2521857C (en) * 2003-04-07 2012-06-12 E.I. Du Pont De Nemours And Company Method and apparatus for quantifying visual showthrough of printed images on the reverse of planar objects
US20050119564A1 (en) * 2003-11-28 2005-06-02 Anders Rosholm Pre-operative planning of implantations
JP5129480B2 (ja) * 2003-09-25 2013-01-30 パイエオン インコーポレイテッド 管状臓器の3次元再構成を行うシステム及び血管撮像装置の作動方法
US20050113680A1 (en) * 2003-10-29 2005-05-26 Yoshihiro Ikeda Cerebral ischemia diagnosis assisting apparatus, X-ray computer tomography apparatus, and apparatus for aiding diagnosis and treatment of acute cerebral infarct

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103109308A (zh) * 2010-09-20 2013-05-15 皇家飞利浦电子股份有限公司 管状结构的管腔特性的量化
CN103109308B (zh) * 2010-09-20 2016-11-30 皇家飞利浦电子股份有限公司 管状结构的管腔特性的量化
CN103505288A (zh) * 2012-06-29 2014-01-15 通用电气公司 超声成像方法和超声成像设备
CN103505288B (zh) * 2012-06-29 2017-11-17 通用电气公司 超声成像方法和超声成像设备

Also Published As

Publication number Publication date
JP2008543482A (ja) 2008-12-04
US20100067753A1 (en) 2010-03-18
WO2006137016A1 (en) 2006-12-28
EP1903943A1 (en) 2008-04-02

Similar Documents

Publication Publication Date Title
CN101203180A (zh) 对血管进行成像的方法和设备
CN101243475B (zh) 具有根据临床任务工作流程的简单点击风格交互的特征的方法和设备
US6792071B2 (en) Method of performing geometric measurements on digital radiological images
CN105637559B (zh) 使用深度传感器的结构建模
US6734847B1 (en) Method and device for processing imaged objects
CN101075346B (zh) 数字图像中用于感兴趣点吸引的方法
EP1349098A1 (en) Method of performing geometric measurements on digital radiological images using graphical templates
US20150348272A1 (en) Image processing apparatus and method for obtaining position and orientation of imaging apparatus
US9542529B2 (en) Medical image data processing apparatus and method
Bornik et al. A hybrid user interface for manipulation of volumetric medical data
RU2706231C2 (ru) Визуализация объемного изображения анатомической структуры
RU2721078C2 (ru) Сегментация анатомической структуры на основе модели
CN104145237A (zh) 用于处理指针移动的系统和方法
CN106960431B (zh) 场景视图的可视化
JP6239694B1 (ja) 画像処理装置、画像処理方法、およびプログラム
US11460990B2 (en) Precise positioning of a marker on a display
KR101847577B1 (ko) 표시 장치용 3차원 인터랙션 방법
JP6293293B2 (ja) マルチセンサ計測装置のルーティンを確立する方法
Teichrieb et al. Virtual reality interfaces applied to correct elevation errors in digital elevation models
JPH0644358A (ja) 画像表示装置
John et al. A Manual Editor for Model Creation from Medical Images

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Open date: 20080618