CN108133512A - 基于磁共振扫描的胎儿体表结构可视化三维成像方法 - Google Patents

基于磁共振扫描的胎儿体表结构可视化三维成像方法 Download PDF

Info

Publication number
CN108133512A
CN108133512A CN201711380521.4A CN201711380521A CN108133512A CN 108133512 A CN108133512 A CN 108133512A CN 201711380521 A CN201711380521 A CN 201711380521A CN 108133512 A CN108133512 A CN 108133512A
Authority
CN
China
Prior art keywords
fetus
data
magnetic resonance
imaging
pixel
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.)
Granted
Application number
CN201711380521.4A
Other languages
English (en)
Other versions
CN108133512B (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201711380521.4A priority Critical patent/CN108133512B/zh
Publication of CN108133512A publication Critical patent/CN108133512A/zh
Application granted granted Critical
Publication of CN108133512B publication Critical patent/CN108133512B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/00—3D [Three Dimensional] image rendering
    • G06T15/08—Volume rendering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00—Indexing scheme for image data processing or generation, in general
    • G06T2200/08—Indexing scheme for image data processing or generation, in general involving all processing steps from image acquisition to 3D model generation

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Public Health (AREA)
  • General Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Computer Graphics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开了一种基于磁共振扫描的胎儿体表结构可视化三维成像方法，包括如下步骤：S1胎儿磁共振数据采集，S2数据转换，S3三维后处理，S4胎儿体表结构三维成像，应用光线投射法完成体绘制显示，从而计算出屏幕上的该像素点处的颜色值。本发明利用真实稳态自由进动梯度回波序列及现有的磁共振成像硬件设施进行胎儿感兴趣区组织结构三维数据采集，通过数据处理，生成包含胎儿组织学特征的新数据，利用三维后处理软件对此新数据进行表面重建，获得胎儿体表结构的三维图像。

Description

基于磁共振扫描的胎儿体表结构可视化三维成像方法
技术领域
本发明涉及医学影像学技术领域，尤其涉及一种基于磁共振扫描的胎儿体表结构可视化三维成像方法。
背景技术
目前胎儿体表结构的三维成像只能依靠三维或四维超声成像来获取，图1为目前三维超声获得的胎儿面部图像；技术手段单一，超声成像条件受到胎儿体位、结构遮挡和羊水多少等诸多因素的影响，有时不能获得满意的胎儿体表三维图像。另外，超声检查的阳性发现也缺少其它影像检查的印证与支持。
磁共振扫描安全、无辐射，具有更高的组织分辨率，目前已广泛应用于胎儿产前疾病诊断，但由于其成像信号特点，使得其三维后处理过程中难以对胎儿体表结构进行表面遮蔽成像及最大强度投影的方法，因而难以通过常规方法获取基于磁共振检查的胎儿体表结构三维图像。
发明内容
本发明的目的是为了解决现有技术中存在的缺点，而提出的一种基于磁共振扫描的胎儿体表结构可视化三维成像方法。
为了实现上述目的，本发明采用了如下技术方案：
基于磁共振扫描的胎儿体表结构可视化三维成像方法，包括如下步骤：
S1胎儿磁共振数据采集，真实稳态自由进动梯度回波序列，增加扫描技术和扫描方位；孕妇取仰卧位、脚先进或头先进，采用腹部柔线圈覆盖中、下腹，采用单次激发快速自旋回波T2加权像序列进行胎儿头部横断位、矢状位、冠状位3个常规磁共振成像，用于观察胎儿面部结构及三维扫描定位；采用真实稳态自由进动梯度回波序列进行胎儿颌面部冠状位扫描，利用冠状位真实稳态自由进动梯度回波序列图像数据进行后处理和重建；
S2数据转换，依次读取MRI扫描序列中每一幅DICOM图像，获取每一个像素灰度值g,统计获得最大的灰度值gmax，设转换后的图像像素灰度为g'，则g'＝gmax-g,用g'组成的像素矩阵与原有DICOM文件的“头信息”组合，重新生成新的DICOM文件序列；
S3三维后处理，采用数据转换后的新图像序列作为可视化的数据来源，将其堆叠成三维体数据；
S4胎儿体表结构三维成像，应用光线投射法完成体绘制显示，从屏幕上的每一像素点处发，根据设定的视点方向和位置，发出一条光线，这条光线穿过体数据空间，在光线与体空间的交点之间进行等距采样，对每个采样点通过由距离某一采样点最近的8个体素的颜色值及不透明度做三线性插值计算求出该采样点的不透明度值及颜色值；在求出该条射线上所有采样点的颜色值及不透明度值以后，再计算在给定的光学模型中的光强，并与前面的采样点结果合成，从而计算出屏幕上的该像素点处的颜色值。
优选的，S1中MRI扫描序列的序列参数：TR/TE：556.60ms-612.00ms/1.80ms-1.93ms，FOV：310×380，Metrix：174×304，Flipangle：79，Bandwidth：484HZ/Px,层厚4mm，间隔-50％～-80％：扫描时间：12s-18s。本发明利用真实稳态自由进动梯度回波序列及现有的磁共振成像硬件设施进行胎儿感兴趣区组织结构三维数据采集，通过数据处理，生成包含胎儿组织学特征的新数据，利用三维后处理软件对此新数据进行表面重建，获得胎儿体表结构的三维图像。
本发明提供一种除超声以外的新的胎儿体表三维成像方法；是现有超声三维成像的必要补充和磁共振成像数据处理的创新；对于胎儿产前诊断(如面部畸形、唇腭裂等)有重要意义。目前我们已经应用此方法成功获取了胎儿面部三维图像，图像质量优良，可以满足胎儿颜面畸形的诊断，且图像直观、显示清晰，应用前景广阔。
附图说明
图1为目前三维超声获得的胎儿面部图像；
图2为常规磁共振图像；
图3为磁共振常规三维后处理得到的图像；
图4是本发明方法得到的孕26周正常胎儿颌面部图像；
图5是本发明方法得到的孕22周单侧唇腭裂胎儿的颌面部图像，箭头示胎儿唇裂位置及大小；
图6是本发明采用光线投射法示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。
参照图2-5，本发明提出的一种基于磁共振扫描的胎儿体表结构可视化三维成像方法，包括如下步骤：
S1胎儿磁共振数据采集，胎儿磁共振数据采集可以使用现有的各种类型磁共振扫描仪进行，技术员经过简单的培训即可掌握扫描的技术要领。
MRI扫描与数据采集要求：
MRI扫描序列：真实稳态自由进动梯度回波序列，目前主流磁共振生产厂家的各机型均配置有该扫描序列，增加扫描技术和扫描方位；磁共振扫描技术(以颌面部成像为例)：孕妇取仰卧位、脚先进或头先进，采用腹部柔线圈覆盖中、下腹(根据感兴趣区位置调整)，采用单次激发快速自旋回波T2加权像序列进行胎儿头部横断位、矢状位、冠状位3个常规磁共振成像，用于观察胎儿面部结构及三维扫描定位，采用真实稳态自由进动梯度回波T2加权像序列进行胎儿颌面部冠状位扫描。
序列参数：TR/TE：556.60ms-612.00ms/1.80ms-1.93ms，FOV：310×380，Metrix：174×304，Flipangle：79，Bandwidth：484HZ/Px,层厚4mm，间隔-50％～-80％：扫描时间：12s-18s。利用冠状位真实稳态自由进动梯度回波T2加权像序列图像数据进行后处理和重建。
S2数据转换：依次读取MRI扫描序列中每一幅DICOM图像，获取每一个像素灰度值g,统计获得最大的灰度值gmax，设转换后的图像像素灰度为g'，则g'＝gmax-g。用g'组成的像素矩阵与原有DICOM文件的“头信息”组合，重新生成新的DICOM文件序列。
S3三维后处理，采用数据转换后的新图像序列作为可视化的数据来源，将其堆叠成三维体数据。
S4胎儿体表结构三维成像，应用光线投射法完成体绘制(volumeRendering,VR)显示。原理如图6所示，从屏幕上的每一像素点处发，根据设定的视点方向和位置，发出一条光线，这条光线穿过体数据空间，在光线与体空间的交点之间进行等距采样，对每个采样点通过由距离某一采样点最近的8个体素的颜色值及不透明度做三线性插值计算求出该采样点的不透明度值及颜色值。在求出该条射线上所有采样点的颜色值及不透明度值以后，再计算在给定的光学模型中的光强，并与前面的采样点结果合成，从而计算出屏幕上的该像素点处的颜色值。
根据光线组合方式，组合分为从后到前和从前到后两种组合方式。设光线穿过体素前的光强和不透明值为cin,αin；穿出体素为cout,αout,采样点为s1,s2......sn。
光线从后到前的组合为
cout＝c(si)*α(si)+cin(1-αin)
其中C(s0)＝Cbackground,α(s0)＝1。
光线由前到后的组合为
αout＝αin+α(si)(1-αin)
其中，在光线从前到后累加的过程中，当发现αout＞1时，停止累加，避免无意义的计算。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。

Claims (2)

1.基于磁共振扫描的胎儿体表结构可视化三维成像方法，其特征在于，包括如下步骤：
S1胎儿磁共振数据采集，真实稳态自由进动梯度回波序列，增加扫描技术和扫描方位；孕妇取仰卧位、脚先进或头先进，采用腹部柔线圈覆盖中、下腹，采用单次激发快速自旋回波T2加权像序列进行胎儿头部横断位、矢状位、冠状位3个正交方位常规磁共振成像，用于观察胎儿面部结构及三维扫描定位；采用真实稳态自由进动梯度回波序列进行胎儿颌面部冠状位扫描，利用冠状位真实稳态自由进动梯度回波序列图像数据进行后处理和重建；
S2数据转换，依次读取MRI扫描序列中每一幅DICOM图像，获取每一个像素灰度值g,统计获得最大的灰度值gmax，设转换后的图像像素灰度为g'，则g'＝gmax-g,用g'组成的像素矩阵与原有DICOM文件的“头信息”组合，重新生成新的DICOM文件序列；
S3三维后处理，采用数据转换后的新图像序列作为可视化的数据来源，将其堆叠成三维体数据；
S4胎儿体表结构三维成像，应用光线投射法完成体绘制显示，从屏幕上的每一像素点处发，根据设定的视点方向和位置，发出一条光线，这条光线穿过体数据空间，在光线与体空间的交点之间进行等距采样，对每个采样点通过由距离某一采样点最近的8个体素的颜色值及不透明度做三线性插值计算求出该采样点的不透明度值及颜色值；在求出该条射线上所有采样点的颜色值及不透明度值以后，再计算在给定的光学模型中的光强，并与前面的采样点结果合成，从而计算出屏幕上的该像素点处的颜色值。
2.根据权利要求1所述的基于磁共振扫描的胎儿体表结构可视化三维成像方法，其特征在于，S1中MRI扫描序列的序列参数：TR/TE：556.60ms-612.00ms/1.80ms-1.93ms，FOV：310×380，Metrix：174×304，Flipangle：79，Bandwidth：484HZ/Px,层厚4mm，间隔-50％～-80％：扫描时间：12s-18s。
CN201711380521.4A 2017-12-20 2017-12-20 基于磁共振扫描的胎儿体表结构可视化三维成像方法 Expired - Fee Related CN108133512B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711380521.4A CN108133512B (zh) 2017-12-20 2017-12-20 基于磁共振扫描的胎儿体表结构可视化三维成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711380521.4A CN108133512B (zh) 2017-12-20 2017-12-20 基于磁共振扫描的胎儿体表结构可视化三维成像方法

Publications (2)

Publication Number Publication Date
CN108133512A true CN108133512A (zh) 2018-06-08
CN108133512B CN108133512B (zh) 2021-08-06

Family

ID=62390782

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711380521.4A Expired - Fee Related CN108133512B (zh) 2017-12-20 2017-12-20 基于磁共振扫描的胎儿体表结构可视化三维成像方法

Country Status (1)

Country Link
CN (1) CN108133512B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109727240A (zh) * 2018-12-27 2019-05-07 深圳开立生物医疗科技股份有限公司 一种三维超声图像的遮挡组织剥离方法及相关装置

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101794460A (zh) * 2010-03-09 2010-08-04 哈尔滨工业大学 基于光线投射体绘制算法的人体心脏三维解剖组织结构模型可视化方法
CN101925382A (zh) * 2008-01-28 2010-12-22 耶达研究及发展有限公司 内窥镜成像光动力学治疗系统及其使用方法
CN102631196A (zh) * 2011-12-08 2012-08-15 中国科学院深圳先进技术研究院 磁共振弹性成像三维可视化方法及系统
CN103021019A (zh) * 2013-01-10 2013-04-03 广东工业大学 一种基于膝关节ct图像的高逼真度模型体绘制方法
CN103068299A (zh) * 2009-09-28 2013-04-24 伊露米纳雷控股有限公司 阴道内监视设备和网络
CN103222879A (zh) * 2012-01-25 2013-07-31 通用电气公司 用于识别超声成像的最佳图像帧的系统和方法
CN103793905A (zh) * 2012-11-02 2014-05-14 上海联影医疗科技有限公司 确定脑部图像冠状位方向的方法及装置
CN103839285A (zh) * 2012-11-20 2014-06-04 三星麦迪森株式会社 用于显示医学图像的方法和设备
CN104123752A (zh) * 2014-07-28 2014-10-29 魏程 一种能够进行3d打印的人体器官的三维建模方法
JP2015013109A (ja) * 2013-06-05 2015-01-22 株式会社東芝 超音波診断装置およびプローブ加減圧情報表示プログラム
CN104887316A (zh) * 2015-04-24 2015-09-09 长春理工大学 基于主动立体显示技术的虚拟立体内窥镜显示方法
US20150297756A1 (en) * 2014-04-21 2015-10-22 Brown University Protein-Based Nanobubble X-Ray Imaging Method and Agent
DE102014208650A1 (de) * 2014-05-08 2015-11-12 Siemens Aktiengesellschaft Verfahren zur Ermittlung eines Synchronisationssignals für eine Bildgebungsvorrichtung
CN105744895A (zh) * 2013-11-21 2016-07-06 三星麦迪森株式会社 显示超声图像的方法和设备

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101925382A (zh) * 2008-01-28 2010-12-22 耶达研究及发展有限公司 内窥镜成像光动力学治疗系统及其使用方法
CN103068299A (zh) * 2009-09-28 2013-04-24 伊露米纳雷控股有限公司 阴道内监视设备和网络
CN101794460A (zh) * 2010-03-09 2010-08-04 哈尔滨工业大学 基于光线投射体绘制算法的人体心脏三维解剖组织结构模型可视化方法
CN102631196A (zh) * 2011-12-08 2012-08-15 中国科学院深圳先进技术研究院 磁共振弹性成像三维可视化方法及系统
CN103222879A (zh) * 2012-01-25 2013-07-31 通用电气公司 用于识别超声成像的最佳图像帧的系统和方法
CN103793905A (zh) * 2012-11-02 2014-05-14 上海联影医疗科技有限公司 确定脑部图像冠状位方向的方法及装置
CN103839285A (zh) * 2012-11-20 2014-06-04 三星麦迪森株式会社 用于显示医学图像的方法和设备
CN103021019A (zh) * 2013-01-10 2013-04-03 广东工业大学 一种基于膝关节ct图像的高逼真度模型体绘制方法
JP2015013109A (ja) * 2013-06-05 2015-01-22 株式会社東芝 超音波診断装置およびプローブ加減圧情報表示プログラム
CN105744895A (zh) * 2013-11-21 2016-07-06 三星麦迪森株式会社 显示超声图像的方法和设备
US20150297756A1 (en) * 2014-04-21 2015-10-22 Brown University Protein-Based Nanobubble X-Ray Imaging Method and Agent
DE102014208650A1 (de) * 2014-05-08 2015-11-12 Siemens Aktiengesellschaft Verfahren zur Ermittlung eines Synchronisationssignals für eine Bildgebungsvorrichtung
CN104123752A (zh) * 2014-07-28 2014-10-29 魏程 一种能够进行3d打印的人体器官的三维建模方法
CN104887316A (zh) * 2015-04-24 2015-09-09 长春理工大学 基于主动立体显示技术的虚拟立体内窥镜显示方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘萍 等: "基于MRI的在体足月胎儿体表数字化三维模型的构建", 《中国医学影像学杂志》 *
李晓清: "医学三维重建图像光线投射并行化算法研究及GPU实现", 《中国硕士学位论文全文数据库信息科技辑》 *
陈小松: "DSA医学图像的增强技术及应用研究", 《中国博硕士学位论文全文数据库（硕士）信息科技辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109727240A (zh) * 2018-12-27 2019-05-07 深圳开立生物医疗科技股份有限公司 一种三维超声图像的遮挡组织剥离方法及相关装置
CN109727240B (zh) * 2018-12-27 2021-01-19 深圳开立生物医疗科技股份有限公司 一种三维超声图像的遮挡组织剥离方法及相关装置

Also Published As

Publication number Publication date
CN108133512B (zh) 2021-08-06

Similar Documents

Publication Publication Date Title
Nelson et al. Visualization of 3D ultrasound data
Nelson et al. Three‐dimensional echocardiographic evaluation of fetal heart anatomy and function: acquisition, analysis, and display.
US6252599B1 (en) Image display method and image display apparatus
Pretorius et al. 3‐dimensional sonographic analysis based on color flow Doppler and gray scale image data: a preliminary report
US20110125016A1 (en) Fetal rendering in medical diagnostic ultrasound
US20090091323A1 (en) Magnetic resonance imaging apparatus and method
Kelly et al. Three-dimensional US of the fetus. Work in progress.
JP2007537770A (ja) 内視鏡画像内の管腔状構造のディスプレイ最適化のための動的なクロップボックス決定方法
CN102106741A (zh) 一种二维超声图像的三维重建方法
Bega et al. Application of three‐dimensional ultrasonography in the evaluation of the fetal heart.
CN102300505A (zh) 超声波诊断装置以及图像数据显示用控制程序
CN103792502A (zh) 用于成像中时序的自动系统
CN104246528A (zh) 具有自动选择记录序列的磁共振成像
Morgan et al. Versatile low-cost volumetric 3-D ultrasound platform for existing clinical 2-D systems
JP3878456B2 (ja) 画像診断支援システム
US7319328B1 (en) System, method, and computer-readable medium for magnetic resonance diffusion anisotropy image processing
CN108226833A (zh) 用于俘获牙科对象的方法
EP2893369A1 (en) Method and apparatus for ultrasound image acquisition
Sepulveda et al. Fetal magnetic resonance imaging and three-dimensional ultrasound in clinical practice: general aspects
JP4406139B2 (ja) Ｍｒｉ装置
JP2007289569A (ja) 医用画像処理装置及び医用画像処理方法
CN108133512A (zh) 基于磁共振扫描的胎儿体表结构可视化三维成像方法
Herickhoff et al. Low-cost sensor-enabled freehand 3d ultrasound
Brookes et al. Three-dimensional imaging of the postmortem fetus by MRI: early experience
JP5575695B2 (ja) Ｍｒｉ装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210806

Termination date: 20211220

CF01 Termination of patent right due to non-payment of annual fee