CN103876764B - 一种血管显影方法及装置 - Google Patents

一种血管显影方法及装置 Download PDF

Info

Publication number
CN103876764B
CN103876764B CN201310594580.7A CN201310594580A CN103876764B CN 103876764 B CN103876764 B CN 103876764B CN 201310594580 A CN201310594580 A CN 201310594580A CN 103876764 B CN103876764 B CN 103876764B
Authority
CN
China
Prior art keywords
blood vessel
detected
image
interest
area
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
CN201310594580.7A
Other languages
English (en)
Other versions
CN103876764A (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.)
Beijing Neusoft Medical Equipment Co Ltd
Original Assignee
Neusoft Medical Systems Co 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 Neusoft Medical Systems Co Ltd filed Critical Neusoft Medical Systems Co Ltd
Priority to CN201310594580.7A priority Critical patent/CN103876764B/zh
Publication of CN103876764A publication Critical patent/CN103876764A/zh
Priority to US14/527,760 priority patent/US20150141818A1/en
Application granted granted Critical
Publication of CN103876764B publication Critical patent/CN103876764B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50—Clinical applications
    • A61B6/504—Clinical applications involving diagnosis of blood vessels, e.g. by angiography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/48—Other medical applications
    • A61B5/4887—Locating particular structures in or on the body
    • A61B5/489—Blood vessels
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means
    • A61B5/742—Details of notification to user or communication with user or patient ; user input means using visual displays
    • A61B5/7425—Displaying combinations of multiple images regardless of image source, e.g. displaying a reference anatomical image with a live image
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means
    • A61B5/7475—User input or interface means, e.g. keyboard, pointing device, joystick
    • A61B5/748—Selection of a region of interest, e.g. using a graphics tablet
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computerised tomographs
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computerised tomographs
    • A61B6/037—Emission tomography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/46—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient
    • A61B6/461—Displaying means of special interest
    • A61B6/466—Displaying means of special interest adapted to display 3D data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/46—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient
    • A61B6/467—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B6/469—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/48—Diagnostic techniques
    • A61B6/481—Diagnostic techniques involving the use of contrast agents
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50—Clinical applications
    • A61B6/501—Clinical applications involving diagnosis of head, e.g. neuroimaging, craniography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
    • A61B6/5241—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT combining overlapping images of the same imaging modality, e.g. by stitching
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5252—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data removing objects from field of view, e.g. removing patient table from a CT image
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging

Landscapes

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

Abstract

本发明实施例公开了一种血管显影方法及装置，其中方法为：扫描待检测血管所在区域得到平扫图像和增强图像，依据所述平扫图像中骨骼区域对所述增强图像进行减影处理得到减影图像；然后，在所述增强图像中检测所述待检测血管的血管区域；最后在所述减影图像中融合所述待检测血管的血管区域得到所述待检测血管的显影图像。可见本发明能够使得减影图像中保留了穿越骨骼的血管。

Description

一种血管显影方法及装置
技术领域
本发明涉及医学图像技术领域，特别是涉及一种血管显影方法及装置。
背景技术
血管疾病是严重影响人类健康的疾病之一，因此对于血管显影技术的要求越来越高，由于人体中血管形态结构复杂，使得部分血管被骨骼遮挡，在三维体绘制结果中很难直观的显示血管的全部形态和结构。比如：头颈部血管从心脏出发穿过颅骨底部骨骼区域进入大脑，颈动脉穿过颅骨、椎动脉穿过六块脊椎骨进入颅骨内部等。
目前，通常利用减影技术减除图像中的骨骼区域，以保留全部的血管区域。现有减影技术的实现是需要对患者进行两次扫描，并且在第二次扫描时，需要对患者注射显影剂，这样就导致第一步扫描与第二部扫描之间的时间间隔较长。由于患者很难长时间保持相同体位，两次扫描过程中患者的被检测区域会发生运动位移，这样会引起血管穿越骨骼的部分被误认为是骨骼区域而减除问题，因此，利用现有减影技术生成的减影图像无法准确完整地显示血管区域的结构。
发明内容
为了解决上述技术问题，本发明中提供了一种血管显影方法及装置，通过检测增强图中的待检测血管区域，将待检测血管区域与减影图像对齐融合得到待检测血管的显影图像。这样处理能够保证穿过骨骼被误认为是骨骼部分而减去的血管图像完整保留下来，提高显影成像的准确性和可靠性。
本发明公开了如下技术方案：
本发明提供了一种血管显影方法，包括：
扫描待检测血管所在区域得到平扫图像和增强图像，依据所述平扫图像中骨骼区域对所述增强图像进行减影处理得到减影图像；
在所述增强图像中检测所述待检测血管的血管区域；
在所述减影图像中融合所述待检测血管的血管区域得到所述待检测血管的显影图像。
优选的，在所述增强图像中检测所述待检测血管的血管区域，包括：
在所述增强图像中根据血管灰度分布情况估计待检测血管的血管路径；
利用灰度平滑度大小关系确定所述血管路径中每一个位置点的血管半径；
按照所述血管半径和血管路径分割所述增强图像得到待检测血管的血管区域。
优选的，所述在所述增强图像中根据血管灰度分布情况估计待检测血管的血管路径，包括：
将所述增强图像与所述待检测血管穿越的骨骼区域模型匹配，确定在所述增强图像中骨骼的位置；
依据骨骼的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管的起止位置；
计算所述增强图像中各位置与所述待检测血管的起止位置之间的灰度相似度，选择灰度相似度差异最小的位置点组成起止位置间的一条血管路径。
优选的，所述利用灰度平滑度大小关系确定所述血管路径中每一个位置点的血管半径，包括：
计算所述血管路径上各位置点的不同尺度半径范围内的灰度平滑度；
选择满足平滑度阈值条件的最大半径作为该路径位置点处的血管半径。
优选的，所述依据骨骼的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管的起止位置，包括：
预先建立至少一种尺度下的血管截面灰度分布模板；
依据所述待检测骨骼的位置估计待检测血管所在的感兴趣区域，在所述感兴趣区域内逐层与所述血管截面灰度分布模板进行匹配得到待检测血管所在的感兴趣区域中的潜在位置；
采用聚类算法对所述待检测血管所在的感兴趣区域的潜在位置进行定位分类，根据待检测血管与所处骨骼的位置起止关系选择最大类的定位点作为待检测血管的起止位置。
优选的，若待检测血管为头颈部的血管时，
所述依据骨骼的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管的起止位置，包括：
预先建立至少一种尺度下的血管截面灰度分布模板；
依据颅骨的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管在感兴趣区域的终止潜在位置；
依据颅骨的位置估计颈部的位置，利用边缘检测算法和圆形检测算子对所述颈部位置的增强图像逐层检测得到待检测血管在感兴趣区域的起始潜在位置；
采用聚类算法对所述起始潜在位置和终止潜在位置分别进行定位分类，分别选择最大类的定位点作为待检测血管的起止位置。
优选的，还包括：
利用三维体绘制显示所述待检测血管的显影图像。
本发明还提供了一种血管显影装置，包括：
扫描单元，用于扫描待检测血管所在区域得到平扫图像和增强图像，依据所述平扫图像中骨骼区域对所述增强图像进行减影处理得到减影图像；
检测单元，用于在所述增强图像中检测所述待检测血管的血管区域；
融合单元，用于在所述减影图像中融合所述待检测血管的血管区域得到所述待检测血管的显影图像。
优选的，所述检测单元，包括：
估计子单元，用于在所述增强图像中根据血管灰度分布情况估计待检测血管的血管路径；
确定子单元，用于利用灰度平滑度大小关系确定所述血管路径中每一个位置点的血管半径；
分割子单元，用于按照所述血管半径和血管路径分割所述增强图像得到待检测血管的血管区域。
优选的，所述估计子单元，包括：
第一匹配模块，用于将所述增强图像与所述待检测血管穿越的骨骼区域模型匹配，确定在所述增强图像中骨骼的位置；
第二匹配模块，用于依据骨骼的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管的起止位置；
第一选择模块，用于计算所述增强图像中各位置与所述待检测血管的起止位置之间的灰度相似度，选择灰度相似度差异最小的位置点组成起止位置间的一条血管路径。
优选的，所述确定子单元，包括：
第一计算模块，用于计算所述血管路径上各位置点的不同尺度半径范围内的灰度平滑度；
第二选择模块，用于选择满足平滑度阈值条件的最大半径作为该路径位置点处的血管半径。
优选的，所述第二匹配模块，包括：
建立子模块，用于预先建立至少一种尺度下的血管截面灰度分布模板；
匹配子模块，用于依据所述待检测骨骼的位置估计待检测血管所在的感兴趣区域，在所述感兴趣区域内逐层与所述血管截面灰度分布模板进行匹配得到待检测血管所在的感兴趣区域中的潜在位置；
选择子模块，用于采用聚类算法对所述待检测血管所在的感兴趣区域的潜在位置进行定位分类，根据待检测血管与所处骨骼的位置起止关系选择最大类的定位点作为待检测血管的起止位置。
优选的，若待检测血管为头颈部的血管时，
所述第二匹配模块，具体包括：
建立子模块，用于预先建立至少一种尺度下的血管截面灰度分布模板；
匹配子模块，用于依据颅骨的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管在感兴趣区域的终止潜在位置；
检测子模块，用于依据颅骨的位置估计颈部的位置，利用边缘检测算法和圆形检测算子对所述颈部位置的增强图像逐层检测得到待检测血管在感兴趣区域的起始潜在位置；
选择子模块，用于采用聚类算法对所述起始潜在位置和终止潜在位置分别进行定位分类，分别选择最大类的定位点作为待检测血管的起止位置。
优选的，还包括：绘制单元，用于利用三维体绘制显示所述待检测血管的显影图像。
本发明的血管显影方法及装置，为了确保血管显影过程中待检测血管穿越骨骼的血管区域被正确显示，首先扫描待检测血管所在区域得到平扫图像和增强图像，依据所述平扫图像中骨骼区域对所述增强图像进行减影处理得到减影图像；然后，在所述增强图像中检测所述待检测血管的血管区域；这样处理能够确定出待检测血管在增强图像中的血管区域此血管区域就是待检测血管穿越骨骼或者临近骨骼部分被误减除的血管区域，最后，在所述减影图像中融合所述待检测血管的血管区域得到所述待检测血管的显影图像。可见本发明的有益效果在于：在增强图像中检测出待检测血管的血管区域，然后在减影图像的基础上融合待检测血管的血管区域图像，保留穿越骨骼或者临近骨骼的血管不被减除，从而确保待检测血管区域结构完整的被显示。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例血管显影方法实施例1的流程图；
图2为本发明实施例血管显影方法实施例2的流程图；
图3为本发明实施例血管显影装置实施例1的示意图；
图4为本发明实施例血管显影装置实施例2的示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明实施例进行详细描述。
需要说明的是本发明的方法可以用于对血管扫描显示的系统，比如：CT（ComputedTomography）电子计算机X射线扫描系统、PET-CT(PositronEmissionTomography)正电子发射断层成像扫描系统、MRI(MagneticResonanceImaging）磁共振成像系统等系统。
参见图1，示出了本发明实施例血管显影方法实施例1的流程图，所述方法包括：步骤101￣步骤103。
步骤101，扫描待检测血管所在区域得到平扫图像和增强图像，依据所述平扫图像中骨骼区域对所述增强图像进行减影处理得到减影图像。
在实际扫描过程中，针对不同的待检测血管确定对应的扫描区域，比如当需要检测心血管时则扫描患者的胸腔区域、需要检测头颈部血管则扫描患者的头颈部区域等。当确定待检测血管所在区域后，首先，扫描待检测区域的断层图像得到平扫图像，可以通过阈值分割平扫图像确定平扫图像中的骨骼区域；然后，对患者注射显影剂再扫描待检测血管所在区域得到增强图像；然后将增强图像与平扫图像按照配准方法对齐，并在增强图像中减除对应的骨骼区域得到减影图像。
步骤102，在所述增强图像中检测所述待检测血管的血管区域。
针对上述步骤102本发明实施例提供了一种实现方式，可以包括：步骤1021￣步骤1023。
步骤1021，在所述增强图像中根据血管灰度分布情况估计待检测血管的血管路径。
步骤1022，利用灰度平滑度大小关系确定所述血管路径中每一个位置点的血管半径。
步骤1023，按照所述血管半径和血管路径分割所述增强图像得到待检测血管的血管区域。
针对上述步骤1021本发明实施例提供了一种确定血管路径的实现方式，可以包括：步骤1021A￣步骤1021C，
步骤1021A，将所述增强图像与所述待检测血管穿越的骨骼区域模型匹配，确定在所述增强图像中骨骼的位置。
本步骤的骨骼区域模型可以使用主动轮廓模型（ActuveShapeModels，ASM）或者主动表现模型（ActiveAppearanceModel，AAM）。本步骤的匹配操作可以采用寻优匹配比如：共轭梯度法、powell寻优法等。
步骤1021B，依据骨骼的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管的起止位置。
步骤1021C，计算所述增强图像中各位置与所述待检测血管的起止位置之间的灰度相似度，选择灰度相似度差异最小的位置点组成起止位置间的一条血管路径。
针对上述步骤1022本发明实施例提供了一种确定血管半径的实现方式，可以包括：步骤1022A和步骤1022B。
步骤1022A，计算所述血管路径上各位置点的不同尺度半径范围内的灰度平滑度；
步骤1022B，选择满足平滑度阈值条件的最大半径作为该路径位置点处的血管半径。
步骤103，在所述减影图像中融合所述待检测血管的血管区域得到所述待检测血管的显影图像。
本步骤的融合可以采用将待检测血管的血管区域按照坐标位置对应的填充到减影图像中。
本发明血管显影过程实质上通过扫描待检测血管所在区域得到增强图像，对所述增强图像进行减影处理得到减影图像；然后在所述增强图像中检测所述待检测血管的血管区域；最后，在所述减影图像中融合所述待检测血管的血管区域得到所述待检测血管的显影图像。可见：本发明将被误删减的血管区域融合在减影图像中，使得在减影图像基础上保留穿越骨骼的血管区域，确保血管区域结构完整。
针对上述步骤1021B本发明实施例还给出了的一种实现方式，可以包括：
预先建立至少一种尺度下的血管截面灰度分布模板；
依据所述待检测骨骼的位置估计待检测血管所在的感兴趣区域，在所述感兴趣区域内逐层与所述血管截面灰度分布模板进行匹配得到待检测血管所在的感兴趣区域中的潜在位置；
采用聚类算法对所述待检测血管所在的感兴趣区域的潜在位置进行定位分类，根据待检测血管与所处骨骼的位置起止关系选择最大类的定位点作为待检测血管的起止位置。
由于人体不同部位的骨骼和血管结构特点复杂程度不同，比如：头颈部、胸部、腿部等部位的血管与骨骼之间的结构完全不同。其中，头颈部的结构最为复杂，血管从心脏出发穿过颅骨底部骨骼区域进入大脑，头颈部区域主要包括四条血管：左颈动脉、右颈动脉、左椎动脉、右椎动脉。为了更快速更准确的检测出这四根血管，最重要的步骤是快速准确的确定出待检测血管的起止位置。因此，为了保证确定头颈部的血管显影的精确性，本发明实施例提供了一种优选方案。
参见图2，示出了本发明实施例血管显影方法实施例2的示意图，所述方法包括：
步骤201，扫描待检测血管所在区域得到平扫图像和增强图像，依据所述平扫图像中骨骼区域对所述增强图像进行减影处理得到减影图像。
步骤202，预先建立至少一种尺度下的血管截面灰度分布模板。
不同患者的血管粗细程度不同，为了满足不同患者的扫描需求，为了使得后续的匹配处理得到潜在中心位置精准度较高，需要获取多种不同尺度下的血管截面灰度分布模板，当然只采用一种常规尺度下的血管截面灰度分布模板也是可以的。
步骤203，依据颅骨的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管在感兴趣区域的终止潜在位置。
在步骤202中确定的颅骨的位置的基础上，估计出这四根血管的终止潜在位置。所谓终止潜在位置是考虑到血管从心脏出发进入脑部，所以将颅骨内部的血管位置称之为终止位置，同时将颈部的血管位置称为起始位置。
根据医学临床经验中人体的颅骨中动脉位置关系，可以依据颅骨的位置估计出待检测血管的潜在位置，比如：左颈动脉在颅骨中心的左侧、左椎动脉在颅骨中心的右侧、两椎动脉在颅骨底部枕骨大孔内，左椎动脉在左颈动脉的右侧，右椎动脉在左椎动脉的右侧。
在实际应用中通常使用颅骨内威利斯环感兴趣区域和基底动脉感兴趣区域，所谓威利斯环感兴趣区域包括：左前交通动脉感兴趣区域和右前交通感兴趣区域。左颈动脉在左前交通动脉感兴趣区域左侧方向延伸，右颈动脉在右前交通感兴趣区域右侧方向延伸。所谓基底动脉感兴趣区域是由左椎动脉和右椎动脉交汇合并成的一个动脉。因此，常常根据颅骨内脑组织的定位确定出威利斯环感兴趣区域和基底动脉感兴趣区域；然后，在这两个感兴趣区域内分别估计这四根血管的潜在位置。
根据颅骨的位置估计出这几根血管的潜在位置，还只是一个大概的区域范围，然后利用血管截面灰度分布模板与这些大概区域范围进行匹配能够得出待检测血管的终止潜在中心位置。这样找出的待检测血管的终止位置较为准确且缩小了查找范围。
步骤204，依据颅骨的位置估计颈部的位置，利用边缘检测算法和圆形检测算子对所述颈部位置的增强图像逐层检测得到待检测血管在感兴趣区域的起始潜在位置。
步骤205，采用聚类算法对所述起始潜在位置和终止潜在位置分别进行定位分类，分别选择最大类的定位点作为待检测血管的起止位置。
上述步骤204和步骤205能够得出四根血管的终止位置、起始位置，终止位置和起始位置分别是位置区域，每个位置区域中都有多个定位点。利用聚类算法将定位点进行分类，剔除错误定位点，将具有最大类的定位点作为待检测血管的起止位置，由于在头颈部区域存在四根待检测血管，可以通过血管解剖结构以及相互位置关系，区分出左颈动脉、右颈动脉、左椎动脉、右椎动脉这四根血管各自的起止位置。
步骤206，计算所述增强图像中各位置与所述待检测血管的起止位置之间的灰度相似度，选择灰度相似度差异最小的位置点组成起止位置间的一条血管路径。
由于颅骨的血管结构和颈部的血管结构不同，采用步骤204的方法也能够确定待检测血管的起始潜在中心位置。为了更为精确的确定待检测血管的起始潜在中心位置，通过步骤205的方法进行检测。
灰度相似度就是指各位置点与起止点的灰度差异大小。在实际应用中可以采用在指定全部点与起止点的最小灰度差异并进行归一化处理，再取倒数变到0和1之间的值作为灰度相似度。
以起始位置为基准位置，首先选择起始位置相邻区域中与起始位置灰度相似度差异最小的位置点，然后再选择该位置点的邻域中与该位置点灰度相似度差异最小的位置点，依次选择出相似度差异最小的位置点，直到选择至终止位置为止，从而根据选择出的所有相似度差异最小的位置点确定一条从起始位置至终止位置的血管路径。
可以利用优先队列检测连接待检测血管起始和终止点间的具有高相似度的路径，比如：先将血管起始点压入优先队列，迭代将待检测点邻域各点放入优先队列，优先队列弹出与起止点相似度最高的位置点，并记录与前一弹出点的邻接关系，若该点为终止点则停止迭代，否则将弹出点邻域内未检测点压入优先队列，继续迭代过程到优先队列弹出点为血管终止点为止。根据优先弹出点的先后顺序及与相邻接关系，回溯获得从起始点到终止点的最大相似度测度的路径，该路径即为连接血管起止点间的血管路径。
当然，同样的也可以先以终止位置为基准位置，依次确定灰度相似度差异最小的位置点，直到确定出起始位置为止。
当然，还可以是从起始位置和终止位置的邻域范围内，选择出多个与其相比相似度值差异最小的位置点，利用这些选择的位置点，确定出一条从起始位置至终止位置的血管路径。
步骤207，利用灰度平滑度大小关系确定所述血管路径中每一个位置点的血管半径。
灰度平滑度是用于描述血管半径范围内的全部点的灰度差异程度，灰度平滑度可以采用不同尺度半径范围内的灰度方差来衡量大小，本步骤通过计算各个位置点在不同尺度半径范围内的灰度平滑度，然后分别将这些灰度平衡度与阈值进行比较并选择出大于阈值的灰度平滑度，再从中选择出最大半径作为该位置点处的血管半径。
步骤208，按照所述血管半径和血管路径分割所述增强图像得到待检测血管的血管区域。
步骤209，在所述减影图像中融合所述待检测血管的血管区域得到所述待检测血管的显影图像。
针对头颈部区域的四根待检测血管，将这四根待检测血管的血管区域处的增强图像与减影图像进行融合，使得减影图像中待检测血管穿越骨骼被减去的部分得到保留。
对应本发明实施例血管显影方法实施例1本发明实施例提供了一种血管显影装置，具体通过以下实施例进行解释说明。
参见图3，示出了本发明实施例血管显影装置实施例1的示意图，所述装置包括：
扫描单元301、检测单元302和融合单元303。下面结合该装置的内部结构和工作原理对各个单元进行解释说明。
扫描单元301，用于扫描待检测血管所在区域得到平扫图像和增强图像，依据所述平扫图像中骨骼区域对所述增强图像进行减影处理得到减影图像；
检测单元302，用于在所述增强图像中检测所述待检测血管的血管区域；
融合单元303，用于在所述减影图像中融合所述待检测血管的血管区域得到所述待检测血管的显影图像。
可选的，所述检测单元，包括：
估计子单元，用于在所述增强图像中根据血管灰度分布情况估计待检测血管的血管路径；
确定子单元，用于利用灰度平滑度大小关系确定所述血管路径中每一个位置点的血管半径；
分割子单元，用于按照所述血管半径和血管路径分割所述增强图像得到待检测血管的血管区域。
在上述检测单元的基础上，可选的，所述估计子单元，包括：
第一匹配模块，用于将所述增强图像与所述待检测血管穿越的骨骼区域模型匹配，确定在所述增强图像中骨骼的位置；
第二匹配模块，用于依据骨骼的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管的起止位置；
第一选择模块，用于计算所述增强图像中各位置与所述待检测血管的起止位置之间的灰度相似度，选择灰度相似度差异最小的位置点组成起止位置间的一条血管路径。
在上述检测单元的基础上所述确定子单元，包括：
第一计算模块，用于计算所述血管路径上各位置点的不同尺度半径范围内的灰度平滑度；
第二选择模块，用于选择满足平滑度阈值条件的最大半径作为该路径位置点处的血管半径。
在上述估计子单元的基础上，可选的，所述第二匹配模块，包括：
建立子模块，用于预先建立至少一种尺度下的血管截面灰度分布模板；
匹配子模块，用于依据所述待检测骨骼的位置估计待检测血管所在的感兴趣区域，在所述感兴趣区域内逐层与所述血管截面灰度分布模板进行匹配得到待检测血管所在的感兴趣区域中的潜在位置；
选择子模块，用于采用聚类算法对所述待检测血管所在的感兴趣区域的潜在位置进行定位分类，根据待检测血管与所处骨骼的位置起止关系选择最大类的定位点作为待检测血管的起止位置。
本发明首先扫描待检测血管所在区域的增强图像，在所述增强图像中去除所述骨骼区域得到减影图像；然后，在所述增强图像中检测所述待检测血管的血管区域；这样处理能够确定出待检测血管在增强图像中的血管区域此血管区域就是待检测血管穿越骨骼或者临近骨骼部分被误减除的血管区域，最后，在所述减影图像中融合所述待检测血管的血管区域得到所述待检测血管的显影图像。
针对头颈部区域的血管进行显影时，本发明还提供了一种血管显影装置，具体是在上述血管显影装置实施1的基础上所述检测单元中的估计子模块中的第二匹配模块与其不同。其他单元和模块都与装置实施例1中的单元相同。具体参见图4，示出了本发明实施例血管显影装置实施例2的示意图，所述装置包括：
该装置中包含扫描单元401、检测单元402和融合单元403。除了检测单元中的第二匹配模块之外，其他单元和模块都与装置实施例1中的单元相同，在此不再赘述。下面对该装置的内部结构和工作原理对各个单元进行解释说明。
检测单元包括：估计子单元、确定子单元和分割子单元；
其中，估计子单元包括：第一匹配模块，第二匹配模块和第一选择模块；若待检测血管为头颈部的血管时，
所述第二匹配模块，具体包括：
建立子模块，用于预先建立至少一种尺度下的血管截面灰度分布模板；
匹配子模块，用于依据颅骨的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管在感兴趣区域的终止潜在位置；
检测子模块，用于依据颅骨的位置估计颈部的位置，利用边缘检测算法和圆形检测算子对所述颈部位置的增强图像逐层检测得到待检测血管在感兴趣区域的起始潜在位置；
选择子模块，用于采用聚类算法对所述起始潜在位置和终止潜在位置分别进行定位分类，分别选择最大类的定位点作为待检测血管的起止位置。
本发明实施例是依据头颈部骨骼和血管的结构特征，利用血管截面灰度分布模板匹配、边缘检测算法和圆形检测算子检测分别确定出待检测血管的潜在终止位置和潜在起始位置，这种方式能够更准确的确定潜在位置，减少计算复杂度，提高显影处理速率。
需要说明的是，在本文中诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
需要说明的是，本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体（Read-OnlyMemory，ROM）或随机存储记忆体（RandomAccessMemory，RAM）等。
以上对本发明所提供的一种血管显影方法和装置进行了详细介绍，本文中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (12)

1.一种血管显影方法，其特征在于，包括：
扫描待检测血管所在区域得到平扫图像和增强图像，依据所述平扫图像中骨骼区域对所述增强图像进行减影处理得到减影图像；
在所述增强图像中检测所述待检测血管的血管区域；
在所述减影图像中融合所述待检测血管的血管区域得到所述待检测血管的显影图像；
其中，在所述增强图像中检测所述待检测血管的血管区域，包括：
在所述增强图像中根据血管灰度分布情况估计待检测血管的血管路径；
利用灰度平滑度大小关系确定所述血管路径中每一个位置点的血管半径；
按照所述血管半径和血管路径分割所述增强图像得到待检测血管的血管区域。
2.根据权利要求1所述的方法，其特征在于，所述在所述增强图像中根据血管灰度分布情况估计待检测血管的血管路径，包括：
将所述增强图像与所述待检测血管穿越的骨骼区域模型匹配，确定在所述增强图像中骨骼的位置；
依据骨骼的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管的起止位置；
计算所述增强图像中各位置与所述待检测血管的起止位置之间的灰度相似度，选择灰度相似度差异最小的位置点组成起止位置间的一条血管路径。
3.根据权利要求1所述的方法，其特征在于，所述利用灰度平滑度大小关系确定所述血管路径中每一个位置点的血管半径，包括：
计算所述血管路径上各位置点的不同尺度半径范围内的灰度平滑度；
选择满足平滑度阈值条件的最大半径作为该路径位置点处的血管半径。
4.根据权利要求2所述的方法，其特征在于，所述依据骨骼的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管的起止位置，包括：
预先建立至少一种尺度下的血管截面灰度分布模板；
依据所述待检测骨骼的位置估计待检测血管所在的感兴趣区域，在所述感兴趣区域内逐层与所述血管截面灰度分布模板进行匹配得到待检测血管所在的感兴趣区域中的潜在位置；
采用聚类算法对所述待检测血管所在的感兴趣区域的潜在位置进行定位分类，根据待检测血管与所处骨骼的位置起止关系选择最大类的定位点作为待检测血管的起止位置。
5.根据权利要求2所述的方法，其特征在于，若待检测血管为头颈部的血管时，
所述依据骨骼的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管的起止位置，包括：
预先建立至少一种尺度下的血管截面灰度分布模板；
依据颅骨的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管在感兴趣区域的终止潜在位置；
依据颅骨的位置估计颈部的位置，利用边缘检测算法和圆形检测算子对所述颈部位置的增强图像逐层检测得到待检测血管在感兴趣区域的起始潜在位置；
采用聚类算法对所述起始潜在位置和终止潜在位置分别进行定位分类，分别选择最大类的定位点作为待检测血管的起止位置。
6.根据权利要求1至5中任一项所述的方法，其特征在于，还包括：
利用三维体绘制显示所述待检测血管的显影图像。
7.一种血管显影装置，其特征在于，包括：
扫描单元，用于扫描待检测血管所在区域得到平扫图像和增强图像，依据所述平扫图像中骨骼区域对所述增强图像进行减影处理得到减影图像；
检测单元，用于在所述增强图像中检测所述待检测血管的血管区域；
融合单元，用于在所述减影图像中融合所述待检测血管的血管区域得到所述待检测血管的显影图像；
其中，所述检测单元，包括：
估计子单元，用于在所述增强图像中根据血管灰度分布情况估计待检测血管的血管路径；
确定子单元，用于利用灰度平滑度大小关系确定所述血管路径中每一个位置点的血管半径；
分割子单元，用于按照所述血管半径和血管路径分割所述增强图像得到待检测血管的血管区域。
8.根据权利要求7所述的装置，其特征在于，所述估计子单元，包括：
第一匹配模块，用于将所述增强图像与所述待检测血管穿越的骨骼区域模型匹配，确定在所述增强图像中骨骼的位置；
第二匹配模块，用于依据骨骼的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管的起止位置；
第一选择模块，用于计算所述增强图像中各位置与所述待检测血管的起止位置之间的灰度相似度，选择灰度相似度差异最小的位置点组成起止位置间的一条血管路径。
9.根据权利要求7所述的装置，其特征在于，所述确定子单元，包括：
第一计算模块，用于计算所述血管路径上各位置点的不同尺度半径范围内的灰度平滑度；
第二选择模块，用于选择满足平滑度阈值条件的最大半径作为该路径位置点处的血管半径。
10.根据权利要求8所述的装置，其特征在于，所述第二匹配模块，包括：
建立子模块，用于预先建立至少一种尺度下的血管截面灰度分布模板；
匹配子模块，用于依据所述待检测骨骼的位置估计待检测血管所在的感兴趣区域，在所述感兴趣区域内逐层与所述血管截面灰度分布模板进行匹配得到待检测血管所在的感兴趣区域中的潜在位置；
选择子模块，用于采用聚类算法对所述待检测血管所在的感兴趣区域的潜在位置进行定位分类，根据待检测血管与所处骨骼的位置起止关系选择最大类的定位点作为待检测血管的起止位置。
11.根据权利要求8所述的装置，其特征在于，若待检测血管为头颈部的血管时，
所述第二匹配模块，具体包括：
建立子模块，用于预先建立至少一种尺度下的血管截面灰度分布模板；
匹配子模块，用于依据颅骨的位置估计待检测血管所在的感兴趣区域，将所述感兴趣区域的增强图像逐层与血管截面灰度分布模板匹配确定待检测血管在感兴趣区域的终止潜在位置；
检测子模块，用于依据颅骨的位置估计颈部的位置，利用边缘检测算法和圆形检测算子对所述颈部位置的增强图像逐层检测得到待检测血管在感兴趣区域的起始潜在位置；
选择子模块，用于采用聚类算法对所述起始潜在位置和终止潜在位置分别进行定位分类，分别选择最大类的定位点作为待检测血管的起止位置。
12.根据权利要求7至11中任一项所述的装置，其特征在于，还包括：
绘制单元，用于利用三维体绘制显示所述待检测血管的显影图像。
CN201310594580.7A 2013-11-21 2013-11-21 一种血管显影方法及装置 Active CN103876764B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201310594580.7A CN103876764B (zh) 2013-11-21 2013-11-21 一种血管显影方法及装置
US14/527,760 US20150141818A1 (en) 2013-11-21 2014-10-29 Vascular imaging method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310594580.7A CN103876764B (zh) 2013-11-21 2013-11-21 一种血管显影方法及装置

Publications (2)

Publication Number Publication Date
CN103876764A CN103876764A (zh) 2014-06-25
CN103876764B true CN103876764B (zh) 2016-03-30

Family

ID=50946124

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310594580.7A Active CN103876764B (zh) 2013-11-21 2013-11-21 一种血管显影方法及装置

Country Status (2)

Country Link
US (1) US20150141818A1 (zh)
CN (1) CN103876764B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017057316A1 (ja) * 2015-09-30 2017-04-06 富士フイルム株式会社 光学フィルム、輝度向上フィルム、輝度向上フィルム付きバックライトユニット及び液晶表示装置
DE102016203857B4 (de) 2016-03-09 2021-05-12 Siemens Healthcare Gmbh Verfahren zur Erfassung und Verarbeitung von Bilddaten eines Untersuchungsobjekts
CN105640583A (zh) * 2016-03-31 2016-06-08 上海联影医疗科技有限公司 血管造影方法
WO2018001099A1 (zh) * 2016-06-30 2018-01-04 上海联影医疗科技有限公司 一种血管提取方法与系统
CN107545584B (zh) * 2017-04-28 2021-05-18 上海联影医疗科技股份有限公司 医学图像中定位感兴趣区域的方法、装置及其系统
CN108514425B (zh) * 2018-05-10 2021-07-30 东软医疗系统股份有限公司 一种造影剂跟踪扫描的方法和装置
CN110755097A (zh) * 2018-07-25 2020-02-07 通用电气公司 医学成像方法和系统、计算机可读存储介质
CN109646035B (zh) * 2019-01-04 2022-04-22 北京永新医疗设备有限公司 骨断层图像重建方法及系统
CN110070083A (zh) * 2019-04-24 2019-07-30 深圳市微埃智能科技有限公司 图像处理方法、装置、电子设备和计算机可读存储介质
CN110448267B (zh) * 2019-09-06 2021-05-25 重庆贝奥新视野医疗设备有限公司 一种多模眼底动态成像分析系统及其方法
CN111062997B (zh) * 2019-12-09 2023-09-12 上海联影医疗科技股份有限公司 一种血管造影成像方法、系统、设备和存储介质
CN111415335B (zh) * 2020-03-11 2023-12-22 北京深睿博联科技有限责任公司 血管标注方法和装置
CN111915538B (zh) * 2020-08-19 2024-03-19 南京普爱医疗设备股份有限公司 一种用于数字血管减影的图像增强方法及系统
CN114533096B (zh) * 2022-02-21 2023-05-19 郑州市中心医院 一种脑血管造影中伪影去除方法及系统
CN117274130B (zh) * 2023-11-22 2024-02-09 中南大学 一种口腔癌医疗影像自适应增强方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1591477A (zh) * 2003-08-28 2005-03-09 株式会社东芝 3d图象处理设备
CN1620990A (zh) * 2003-11-25 2005-06-01 通用电气公司 在ct血管造影术中分割结构的方法及设备
CN1983333A (zh) * 2005-11-17 2007-06-20 通用电气公司 用于选择和/或标记管分枝的方法和装置
TW201225922A (en) * 2010-12-31 2012-07-01 Univ Nat Central Method of CT angiography to visualize trans-osseous blood vessels

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4542459A (en) * 1982-11-26 1985-09-17 General Electric Company Matched filter for x-ray hybrid subtraction
JP2536024B2 (ja) * 1988-02-29 1996-09-18 株式会社島津製作所 Ｘ線画像処理装置
FR2829268A1 (fr) * 2001-09-04 2003-03-07 Koninkl Philips Electronics Nv Procede de traitement d'images pour angiographie soustractiv e numerisee
EP1938271A2 (en) * 2005-10-21 2008-07-02 The General Hospital Corporation Methods and apparatus for segmentation and reconstruction for endovascular and endoluminal anatomical structures
US8090177B2 (en) * 2008-08-01 2012-01-03 Sti Medical Systems, Llc Methods for detection and characterization of atypical vessels in cervical imagery
DE102009006636B4 (de) * 2008-12-30 2016-02-18 Siemens Aktiengesellschaft Verfahren zur Ermittlung einer 2D-Kontur einer in 3D-Bilddaten abgebildeten Gefäßstruktur
US10217238B2 (en) * 2009-05-20 2019-02-26 Algotec Systems Ltd. Soft tissue segmentation using a bone atlas
US20120022366A1 (en) * 2010-07-21 2012-01-26 Marcus Pfister Registration of aorta to patient via two 2d images for placement of a stent
JP2012040204A (ja) * 2010-08-19 2012-03-01 Toshiba Corp Ｘ線画像撮影装置
US9595088B2 (en) * 2013-11-20 2017-03-14 Toshiba Medical Systems Corporation Method of, and apparatus for, visualizing medical image data

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1591477A (zh) * 2003-08-28 2005-03-09 株式会社东芝 3d图象处理设备
CN1620990A (zh) * 2003-11-25 2005-06-01 通用电气公司 在ct血管造影术中分割结构的方法及设备
CN1983333A (zh) * 2005-11-17 2007-06-20 通用电气公司 用于选择和/或标记管分枝的方法和装置
TW201225922A (en) * 2010-12-31 2012-07-01 Univ Nat Central Method of CT angiography to visualize trans-osseous blood vessels

Also Published As

Publication number Publication date
US20150141818A1 (en) 2015-05-21
CN103876764A (zh) 2014-06-25

Similar Documents

Publication Publication Date Title
CN103876764B (zh) 一种血管显影方法及装置
JP6703082B2 (ja) 医用画像における構造物をセグメンテーションする医用画像処理装置、医用画像をセグメンテーションするための方法及び医用画像をセグメンテーションするコンピュータプログラムを記憶する記憶媒体
CN103339653B (zh) 勾画组织损伤的方法
Sundaresan et al. Automated characterization of the fetal heart in ultrasound images using fully convolutional neural networks
JP6059261B2 (ja) マルチモーダル画像統合におけるレジストレーション精度を改善するインテリジェントなランドマーク選択
EP2149121B1 (en) Image registration method
CN103732134B (zh) 用于自动左心室功能评价的系统、装置、设备和方法
US8218839B2 (en) Automatic localization of the left ventricle in cardiac cine magnetic resonance imaging
US9002078B2 (en) Method and system for shape-constrained aortic valve landmark detection
US20080247622A1 (en) Methods, Systems, and Computer Program Products For Hierarchical Registration Between a Blood Vessel and Tissue Surface Model For a Subject and a Blood Vessel and Tissue Surface Image For the Subject
CN114170134B (zh) 基于颅内dsa影像的狭窄评估方法及装置
JP2004535874A (ja) 磁気共鳴血管造影法及びその装置
Yang et al. Renal compartment segmentation in DCE-MRI images
Neubert et al. Automated 3D segmentation of vertebral bodies and intervertebral discs from MRI
Suri et al. Angiography and plaque imaging: advanced segmentation techniques
US20210035290A1 (en) Flow analysis in 4d mr image data
CN107507212B (zh) 数字脑可视化方法、装置、计算设备及存储介质
WO2014017853A2 (ko) Ｔｏｆ-ｍｒａ를 이용한 혈류특성 및 ｍｒ-신호강도구배(전단율) 유도방법
CN109886966A (zh) 一种从超声波图像中提取目标特征的图像处理方法
JP2013223620A (ja) 医用画像のレジストレーション処理装置、方法およびプログラム
US20190038150A1 (en) Method for analyzing blood flow by using medical image
KR20150115438A (ko) 심장 부위를 포함하는 의료 영상 오버레이 방법 및 그 장치
Likhite et al. Deformable and rigid model-based image registration for quantitative cardiac perfusion
Kuijf et al. Assessment of methods to extract the mid-sagittal plane from brain MR images
CN111612867B (zh) 运动伪影校正方法、装置、计算机设备和可读存储介质

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20170307

Address after: Room 6, building 3, building No. 8, Wang Xi Road, Haidian District, Beijing, 100193, China

Patentee after: Beijing Neusoft Medical Equipment Co., Ltd.

Address before: Hunnan New Century Road 110179 Shenyang city of Liaoning Province, No. 16

Patentee before: Dongruan Medical Systems Co., Ltd., Shenyang

TR01 Transfer of patent right