CN102764124B - 一种基于磁共振成像对穿支皮瓣血管定位与测量的方法 - Google Patents
一种基于磁共振成像对穿支皮瓣血管定位与测量的方法 Download PDFInfo
- Publication number
- CN102764124B CN102764124B CN 201210235419 CN201210235419A CN102764124B CN 102764124 B CN102764124 B CN 102764124B CN 201210235419 CN201210235419 CN 201210235419 CN 201210235419 A CN201210235419 A CN 201210235419A CN 102764124 B CN102764124 B CN 102764124B
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- perforator flap
- length
- flap blood
- artery perforator
- 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.)
- Expired - Fee Related
Links
- 210000004204 blood vessel Anatomy 0.000 title claims abstract description 88
- 238000002595 magnetic resonance imaging Methods 0.000 title abstract 2
- 238000000691 measurement method Methods 0.000 title abstract 2
- 238000000034 method Methods 0.000 claims abstract description 21
- 238000003384 imaging method Methods 0.000 claims abstract description 8
- 206010033675 panniculitis Diseases 0.000 claims abstract description 7
- 238000006073 displacement reaction Methods 0.000 claims abstract description 5
- 210000003205 muscle Anatomy 0.000 claims abstract description 4
- 210000001367 artery Anatomy 0.000 claims description 59
- 238000005481 NMR spectroscopy Methods 0.000 claims description 9
- 239000002872 contrast media Substances 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 3
- 238000003709 image segmentation Methods 0.000 claims description 3
- 210000004304 subcutaneous tissue Anatomy 0.000 abstract 1
- 238000002372 labelling Methods 0.000 description 8
- 230000003187 abdominal effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- 210000001015 abdomen Anatomy 0.000 description 2
- 230000036770 blood supply Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- 210000000577 adipose tissue Anatomy 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000002316 cosmetic surgery Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003387 muscular Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 210000001113 umbilicus Anatomy 0.000 description 1
- 210000005166 vasculature Anatomy 0.000 description 1
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种基于磁共振成像对穿支皮瓣血管定位与测量的方法,包括以下步骤:采集得到人体皮下组织的磁共振图像,在图像中选取参考点;寻找图像中穿支皮瓣血管穿入肌肉处的位置,设置标记点;读取参考点与标记点的成像坐标,并计算两者之间的位移;对血管进行追踪和分割,并计算出穿支皮瓣血管的长度和直径;并将所有标记点投影至皮肤表面,得到皮肤表面的穿支皮瓣血管标记图像。本发明方法能够精确定位,并精确测量穿支血管的长度和直径大小,所需人工操作的时间少且效率高。
Description
技术领域
本发明方法属于磁共振成像的技术领域,尤其涉及一种基于磁共振成像对穿支皮瓣血管定位与测量的方法。
背景技术
穿支血管是指茎状细小、连接并穿入肌肉组织中的血管,对皮瓣组织供血及再生皮下组织的功能。穿支皮瓣是指以穿支血管供血、由皮肤和脂肪组织构成的大型皮瓣组织。在现今临床的整形外科中,穿支皮瓣手术有着越来越广泛的应用,如植皮,乳房再造手术等。准确测量穿支皮瓣血管位置和直径大小能给外科医生提供有效的解剖信息,有利于穿支皮瓣手术的进行。
传统的穿支皮瓣血管定位方法是利用脉冲多普勒超声仪器在对人体皮下组织进行探测,通过声音信号强弱来判定穿支皮瓣血管位置。该方法的缺点是定位不精确,而且无法反映穿支血管的长度和直径大小的信息。
发明内容
本发明针对现有的定位穿支皮瓣血管的技术方法的不足,提出一种借助于磁共振成像手段以及磁共振图像显示和后处理技术,对穿支皮瓣血管进行精确定位,测量血管长度和直径并对穿支皮瓣血管显示标记的方法。
本发明公开了一种基于磁共振成像对穿支皮瓣血管定位与测量的方法,包括以下步骤:
步骤一:利用磁共振成像技术采集得到人体皮下组织的造影剂增强T1加权磁共振图像,在所述图像中通过感兴趣区工具选取参考点;
步骤二:寻找所述图像中穿支皮瓣血管穿入肌肉处的位置,并用通过所述感兴趣区工具设置标记点;
步骤三:利用计算机程序读取所述参考点与标记点的成像坐标,并计算所述标记点到所述参考点的位移;
步骤四:利用图像分割算法对所述被标记的穿支皮瓣血管进行追踪和分割,根据追踪得到的血管边界长度计算出所述被标记的穿支皮瓣血管的长度,根据分割得到的血管边界长度和血管面积计算出所述被标记的穿支皮瓣血管和直径;
步骤五:利用3D表面皮肤重建技术对磁共振图像进行皮肤重建显示,并将所有所述标记点投影至皮肤表面,得到皮肤表面的穿支皮瓣血管标记图像。
其中,所述步骤三进一步包括:对所述标记点进行命名。
其中,依据所述标记点相对于所述参考点的方位命名所述标记点。
其中,所述步骤四中,所述被标记的穿支皮瓣血管的长度为所述追踪得到的血管边界长度的一半,以如下公式(I)表示:
L=B1/2 (I)
其中,B1为所述追踪得到的血管边界长度,L为所述被标记的穿支皮瓣血管的长度。
其中,所述步骤四中,所述被标记的穿支皮瓣血管的直径为所述分割得到的血管边界长度和血管面积的商,以如下公式(II)表示:
D=A/(B2/2) (II)
其中,D为所述被标记的穿支皮瓣血管的直径,A为所述分割得到的血管面积,B2为所述分割得到的血管边界长度。
相较于传统的利用多普勒超声仪器对穿支皮瓣血管进行定位的方法,本发明利用磁共振成像技术采集得到人体皮下组织的造影剂增强高对比度T1加权图像,并通过磁共振图像显示和后处理技术,不仅定位精确,而且能精确测量穿支血管的长度和直径大小。此外,利用本发明方法定位及测量穿支皮瓣血管大部分过程为计算机自动测量过程,只需借助少量人工手动介入,所需时间少且效率高。
附图说明
图1是本发明实施例中对穿支皮瓣血管的边界检测的示意图。
图2是本发明实施例中对穿支皮瓣血管进行正方形范围分段截取的示意图。
图3是本发明实施例中对穿支皮瓣血管标记点的皮肤表面投影的示意图。
图4是本发明实施例中穿支皮瓣血管在3D重建得到的皮肤表面图上的显示和标记的示意图。
具体实施方式
下面结合附图和实施例对本发明的具体实施方式做进一步详细的说明,但不应以此限制本发明的保护范围。
本发明基于磁共振成像对穿支皮瓣血管定位与测量的方法,包括以下步骤:
步骤一:利用磁共振成像技术采集得到人体皮下组织的造影剂增强T1加权磁共振图像,在图像中通过感兴趣区工具选取一个参考点。
步骤二:寻找图像中穿支皮瓣血管穿入肌肉处的位置,并用通过感兴趣区工具设置标记点。
步骤三:利用计算机程序读取参考点与标记点的成像坐标,并计算标记点到参考点的位移。对所有穿支血管进行命名标记,命名规则依据相对于参考点左右位置,分别命名为L和R,依据坐标位置从上到下的大小顺序将穿支血管标记为IL1、L2、L3,R1、R2、R3等。
步骤四:利用图像分割算法对被标记的穿支皮瓣血管进行追踪和分割,分别根据追踪得到的血管边界长度,以及分割得到的血管边界长度和血管面积计算出被标记的穿支皮瓣血管的长度和直径。
其中,被标记穿支皮瓣血管的长度为追踪得到的血管边界长度的一半,以如下公式(I)表示:
L=B1/2 (I)
式(I)中,B1为追踪得到的血管边界长度,L为被标记穿支皮瓣血管的长度。
被标记的穿支皮瓣血管的直径为分割得到的血管边界长度和血管面积的商,以如下公式(II)表示:
D=A/(B2/2) (II)
式(II)中,D为被标记的穿支皮瓣血管的直径,A为分割得到的血管面积,B2为分割得到的血管边界长度。
步骤五:利用3D表面皮肤重建技术对磁共振图像进行皮肤重建显示,并将所有标记点投影至皮肤表面,得到皮肤表面的穿支皮瓣血管标记图像。
实施例:
本实施例采集的磁共振图像数据为人体腹部的造影剂增强的T1加权图像,数据来源于GE1.5T磁共振成像仪。
步骤一:利用磁共振成像技术采集得到人体皮下组织的造影剂增强T1加权磁共振图像。所使用的成像序列为脂肪信号受抑制的梯度回波序列,其成像参数如下:恢复时间为6.064ms,回波时间为2.92ms,倾倒角为15度,观察野为480*480mm,层厚为3mm,扫描矩阵为512*512。所使用的图像显示工具为0siriX软件。
步骤二:选取参考点并标记穿支皮瓣血管。
本实施例选取人体腹部中肚脐为参考点,并用借助磁共振图像软件系统中的感兴趣区工具(ROI)在磁共振图像中的腹部左侧位置找到并标记了三处穿支皮瓣血管,在腹部右侧位置处找到并标记了两处穿支皮瓣血管。
步骤三:对标记处进行命名。
利用编写好的0siriX的插件程序自动读取出包括参考点的所有标记点的成像坐标,并计算各穿支血管标记点到参考点的位移。对所有穿支血管进行命名标记,命名规则依据相对于参考点左右位置,分别命名为L和R,依据坐标位置从上到下的大小顺序将穿支血管标记为L1、L2、L3以及R1、R2、R3等。本实施例在腹部左侧位置有三处穿支皮瓣血管,分别标记为L1、L2、L3,在腹部右侧位置处有两处穿支皮瓣血管,分别标记为R1、R2,如图4所示。
步骤四:自动测量穿支血管的长度。
如图1,对所有标记点处的邻近区域利用图像梯度算子进行边界检测,检测出血管的边界并做追踪,以此检测并计算出穿支皮瓣血管的边界长度B,将边界长度B的一半值定义为穿支皮瓣血管的长度L=B1/2,如公式(I)所示。
步骤五:自动测量穿支血管的直径。
将标记点周围11×11像素正方形范围内的穿支血管进行分断截取。如图2所示,利用对步骤三中已检测出的边界做截取,保留在正方形范围内的血管边界,并计算该部分边界长度B2。此时该部分的血管长度为B2/2。计算边界所包围的血管部分的面积A。最终计算血管直径D=A/(B2/2)。步骤四中五处穿支血管的长度与直径计算值分别如表1所示:
表1图4中标记的穿支血管的长度与直径
步骤五:利用3D皮肤重建对穿支皮瓣血管的显示标记。
如图3,首先将所有标记的穿支皮瓣血管标记直线投影至皮肤表面。然后利用3D皮肤重建技术对图像进行皮肤重建,最后穿支皮瓣血管的位置在皮肤表面图像上得到显示和标记。穿支皮瓣血管在3D重建得到的皮肤表面图上的显示和标记的实例见图4。
以上所述仅为本发明的较佳实施例,并非用来限定本发明的实施范围。任何所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作各种变动与润饰,本发明保护范围应以权利要求书所界定的保护范围为准。
Claims (3)
1.一种基于磁共振成像对穿支皮瓣血管定位与测量的方法,其特征在于,包括以下步骤:
步骤一:利用磁共振成像技术采集得到人体皮下组织的造影剂增强T1加权磁共振图像,在所述图像中通过感兴趣区工具选取参考点;
步骤二:寻找所述图像中穿支皮瓣血管穿入肌肉处的位置,并用通过所述感兴趣区工具设置标记点;
步骤三:利用计算机程序读取所述参考点与标记点的成像坐标,并计算所述标记点到所述参考点的位移;
步骤四:利用图像分割算法对被标记的穿支皮瓣血管进行追踪和分割,根据追踪得到的血管边界长度计算出所述被标记的穿支皮瓣血管的长度,根据分割得到的血管边界长度和血管面积计算出所述被标记的穿支皮瓣血管的直径;
步骤五:利用3D表面皮肤重建技术对磁共振图像进行皮肤重建显示,并将所有所述标记点投影至皮肤表面,得到皮肤表面的穿支皮瓣血管标记图像;
其中,步骤四中,
所述被标记的穿支皮瓣血管的长度为所述追踪得到的血管边界长度的一半,以如下公式(I)表示:
L=B1/2 (I)
B1为所述追踪得到的血管边界长度,L为所述被标记的穿支皮瓣血管的长度;
所述被标记的穿支皮瓣血管的直径为所述分割得到的血管面积和血管边界长度的商,以如下公式(II)表示:
D=A/(B2/2) (II)
其中,D为所述被标记的穿支皮瓣血管的直径,A为所述分割得到的血管面积,B2为所述分割得到的血管边界长度。
2.如权利要求1所述基于磁共振成像对穿支皮瓣血管定位与测量的方法,其特征在于,所述步骤三进一步包括:对所述标记点进行命名。
3.如权利要求2所述基于磁共振成像对穿支皮瓣血管定位与测量的方法,其特征在于,依据所述标记点相对于所述参考点的方位命名所述标记点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201210235419 CN102764124B (zh) | 2012-07-09 | 2012-07-09 | 一种基于磁共振成像对穿支皮瓣血管定位与测量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201210235419 CN102764124B (zh) | 2012-07-09 | 2012-07-09 | 一种基于磁共振成像对穿支皮瓣血管定位与测量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102764124A CN102764124A (zh) | 2012-11-07 |
CN102764124B true CN102764124B (zh) | 2013-12-25 |
Family
ID=47091869
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201210235419 Expired - Fee Related CN102764124B (zh) | 2012-07-09 | 2012-07-09 | 一种基于磁共振成像对穿支皮瓣血管定位与测量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102764124B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106345062B (zh) * | 2016-09-20 | 2018-01-16 | 华东师范大学 | 一种基于磁共振成像的经颅磁刺激线圈定位方法 |
US10893849B2 (en) * | 2017-02-28 | 2021-01-19 | Canon Medical Systems Corporation | Ultrasound image diagnosis apparatus, medical image processing apparatus, and computer program product |
CN107616839A (zh) * | 2017-11-10 | 2018-01-23 | 岳春华 | 一种基于3d打印重建的皮瓣穿支动脉定位方法 |
CN111358553B (zh) * | 2020-03-11 | 2021-08-10 | 宁波市第一医院 | 一种穿支皮瓣构建方法、装置、计算机设备及可存储介质 |
CN114463557B (zh) * | 2022-01-25 | 2023-03-24 | 腾讯科技(深圳)有限公司 | 轮廓标记方法、装置、设备、存储介质及系统 |
CN115456950B (zh) * | 2022-08-15 | 2023-08-11 | 上海交通大学医学院附属第九人民医院 | 一种穿支血管自动识别和定位方法、设备及存储介质 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040133129A1 (en) * | 2003-01-03 | 2004-07-08 | Mindguard Ltd. | Measurements in a body lumen using guidewire with spaced markers |
JP2005198708A (ja) * | 2004-01-13 | 2005-07-28 | Toshiba Corp | 血管狭窄率解析装置及び血管狭窄率解析方法 |
US7970187B2 (en) * | 2005-09-06 | 2011-06-28 | Pie Medical Imaging B.V. | Method, apparatus and computer program for contour detection of vessels using x-ray densitometry |
CN101140659B (zh) * | 2007-09-29 | 2010-08-04 | 华中科技大学 | 在数字血管造影图像中分割血管数据的方法 |
BRPI0907272A2 (pt) * | 2008-04-14 | 2015-07-21 | Novadaq Technologies Inc | "método de avaliação e aparelho de avaliação de perfusão de vasos perfurantes" |
CN100571637C (zh) * | 2008-05-28 | 2009-12-23 | 华中科技大学 | 动态模型指导下的血管造影三维重建方法 |
CN101984916B (zh) * | 2010-11-17 | 2012-10-31 | 哈尔滨工程大学 | 基于数字图像处理技术的血管管径测量方法 |
-
2012
- 2012-07-09 CN CN 201210235419 patent/CN102764124B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102764124A (zh) | 2012-11-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102764124B (zh) | 一种基于磁共振成像对穿支皮瓣血管定位与测量的方法 | |
CN106344153B (zh) | 一种柔性穿刺针针尖自动跟踪装置及方法 | |
CN103997971B (zh) | 用于超声心动描记的自动成像平面选择 | |
CN103735268B (zh) | 一种体位检测方法及系统 | |
US8358818B2 (en) | Apparatus and methods of compensating for organ deformation, registration of internal structures to images, and applications of same | |
JP5543444B2 (ja) | 生検実施のための方法及びシステム | |
CN204909663U (zh) | 一种基于图像的外科手术导航系统 | |
CN103040525A (zh) | 一种多模医学影像手术导航方法及系统 | |
RU2014102345A (ru) | Анатомическая пометка обнаруженных изменений в данных изображений, полученных в ходе длительных наблюдений | |
Hill et al. | Correcting scaling errors in tomographic images using a nine degree of freedom registration algorithm | |
Orczyk et al. | Preliminary experience with a novel method of three-dimensional co-registration of prostate cancer digital histology and in vivo multiparametric MRI | |
CN112168346A (zh) | 三维医学影像与患者实时重合的方法及手术辅助系统 | |
CN106456084A (zh) | 超声成像装置 | |
CN107106128A (zh) | 用于分割解剖目标的超声成像装置和方法 | |
Li et al. | Augmenting intraoperative ultrasound with preoperative magnetic resonance planning models for percutaneous renal access | |
Wen et al. | A novel ultrasound probe spatial calibration method using a combined phantom and stylus | |
CN114404039B (zh) | 三维模型的组织漂移校正方法、装置、电子设备及存储介质 | |
CN205964158U (zh) | 一种与肺部手术导航系统配用的胸部约束带 | |
CN103892918B (zh) | 定位装置、影像叠合方法及影像投影系统 | |
CN116236280B (zh) | 一种基于多模态图像融合的介入治疗引导方法及系统 | |
JP6944492B2 (ja) | 画像取得方法、関連装置及び読み取り可能な記憶媒体 | |
CN214157490U (zh) | 一种应用三维医学影像与患者实时重合方法的手术辅助系统 | |
Chaoui et al. | Virtual movements-based calibration method of ultrasound probe for computer assisted surgery | |
Chen et al. | Three dimensional ultrasound guided percutaneous renal puncture: A phantom study | |
Li et al. | Comparison of 2d and 3d ultrasound guided percutaneous renal puncture |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131225 |
|
CF01 | Termination of patent right due to non-payment of annual fee |