CN100586368C - 基于人体解剖结构对称性的磁共振成像扫描自动定位方法 - Google Patents
基于人体解剖结构对称性的磁共振成像扫描自动定位方法 Download PDFInfo
- Publication number
- CN100586368C CN100586368C CN200610046121A CN200610046121A CN100586368C CN 100586368 C CN100586368 C CN 100586368C CN 200610046121 A CN200610046121 A CN 200610046121A CN 200610046121 A CN200610046121 A CN 200610046121A CN 100586368 C CN100586368 C CN 100586368C
- Authority
- CN
- China
- Prior art keywords
- symmetry
- image
- target area
- axis
- mentioned
- 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
Links
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开一种基于人体解剖结构对称性的磁共振成像扫描自动定位方法,包含如下步骤:a)获取MRI定位图像;b)利用人体解剖结构的对称性,在MRI定位图像上自动计算不同方位的定位信息;c)输出MRI定位图像的定位信息;其中步骤b)还包括:1)在整幅图像中分割出目标区域;2)初始化上述目标区域的对称轴;3)计算上述目标区域中关于上述对称轴对称的点集;4)基于上述点集定义代价函数,通过优化处理得到上述代价函数最优值;5)得到与上述代价函数最优值对应的作为定位信息的对称轴。本发明实现了MRI扫描的自动准确定位,保证断层扫描图像的一致性,提高了定位精度和工作效率。
Description
技术领域
本发明涉及一种医学影像处理定位方法,具体地说是基于人体解剖结构对称性的磁共振成像扫描自动定位方法。
背景技术
在MRI(Magnetic Resonance Imaging,磁共振成像)扫描中,需要预先扫描获得定位图像,在定位图像中设定参考定位线,由此来确定MRI断层扫描的位置和方向。
现有技术中,一般会在定位图像中由系统自动生成以定位图像的几何中心为中心、与x轴或者y轴方向平行的参考定位线,参考定位线的位置和方向决定了MRI断层扫描的位置和方向。但是在实际应用中,临床要求参考定位线的位置和方向与扫描的部位相对应,如以脑横断面的图像作为定位图像,在此定位图像上设定扫描矢状面断层的参考定位线时,参考定位线的中心位于大脑裂中心位置、方向与大脑裂平行;而以腰椎冠状面的图像作为定位图像,在此定位图像上设定扫描矢状面断层的参考定位线时,参考定位线的中心位于脊髓中心位置、方向与脊髓平行。因此,仅以定位图像的几何中心为中心、与x或者y轴方向平行的参考定位线不能满足临床的需求。
为解决这个问题,现有技术中又提供了手工调整参考定位线的位置和方向的方法,以满足对不同部位进行MRI断层扫描的需要。但是,采用手工调整参考定位线的位置和方向存在以下不足:
1、由于不可避免的人为误差,无法保证参考定位线的一致性;
2、人工操作定位精度低;
3、增大了医生的工作量,浪费扫描时间,降低了医生的工作效率。
发明内容
为了克服上述不足,本发明的目的在于提供一种方便医生操作、缩短病人的总扫描时间的基于人体解剖结构对称性的磁共振成像扫描自动定位方法。
为实现上述目的,本发明采用的技术方案是:
包含如下步骤:a)获取MRI定位图像;b)利用人体解剖结构的对称性,在MRI定位图像上自动计算不同方位的定位信息;c)输出MRI定位图像的定位信息。
上述步骤b)具体为:
1)在整幅图像中分割出目标区域;2)初始化上述目标区域的对称轴;3)计算上述目标区域中关于上述对称轴对称的点集;4)基于上述点集定义代价函数,通过优化处理得到上述代价函数最优值;5)得到与上述代价函数最优值对应的作为定位信息的对称轴;
上述步骤2)初始化上述目标区域的对称轴是基于计算上述目标区域的主轴进行的;上述步骤4)代价函数定义为上述目标区域中对称点灰度差的绝对值之和。
本发明具有以下有益效果及优点:
1.实现了MRI扫描的自动准确定位。由于本发明的定位线是基于人体解剖结构对称性特征自动生成,可一次性满足临床定位需求,不需要再进行手工调整;
2.保证断层扫描图像的一致性。本发明可以自动地、准确地计算出目标区域的定位线,保证同一名医生在不同的时间内或者不同医生在同一时间内对同一部位定位时的一致性,从而确保断层扫描图像的一致性;
3.提高定位精度。本发明在临床应用方面可降低手工定位时的人为误差,准确地计算出目标区域的定位线,尤其是在定位图像成像质量不佳、目标区域特征难以辨认时,定位效果依然很好;
4.提高了工作效率:本发明由于可以一次性自动确定满足需要的定位线,不用手工反复调整,因而可以减少病人的扫描时间,提高医生的工作效率;
5.应用范围广:本发明除可应用于脑横断面、颈椎冠状面、腰椎冠状面的定位外,还可以扩展推广到其他具备解剖结构对称性的部位。
附图说明
图1为本发明方法流程图;
图2为本发明图像对称性原理说明图;
图3~5为本发明的效果图(一)~(三)。
具体实施方式
如图1所示,本发明具体包括如下步骤:
S10:在获取的MRI整幅图像中分割出目标区域,即使用图像分割法,如双峰法、大津法或迭代法等(本实施例采用大津法),计算图像中目标区域和背景区域分割的阈值,使用此阈值分割出图像中的目标区域;
S20:计算目标区域的主轴,把主轴作为初始对称轴,对称轴方程可表示为:y=ax+b;
本实施例中,主轴作为对称轴的初始值,计算主轴的部分可以分为两步,一步是先计算目标区域的重心(x,y),作为主轴上的点;第二步是计算主轴的角度α0,计算主轴的角度可以根据公式 计算,其中:Ix和Iy是惯性距,Ixy是惯性积,它们的计算公式是Ix=∫Ay2dA,Iy=∫Ax2dA和Ixy=∫AxydA,A为目标区域;
S30:计算目标区域中关于对称轴对称的点集;在对称轴垂直的方向上找与对称轴的距离相等的点,这些点就是对称点,所有对称点就组成了对称点集;
S40:计算代价函数f(u,v),即计算目标区域中对称点灰度差的绝对值之和:
其中:G(xi0,yi0,z)、G(xi1,yi1,z)分别是对称点的灰度值,n为目标区域内对称点组数,z为对称轴;
S50:判断步骤S40中代价函数值是否最优(本实施例采用最小值为最优值);
S60:当代价函数值最小时,搜索完成,输出与上述代价函数最优值对应的作为定位信息的对称轴;
S70:否则修正对称轴,返回S30,重新计算代价函数进行搜索。
上述步骤S30至步骤S50是以目标区域的主轴方向和重心计算出的主轴作为对称轴的初始值,使用最优化搜索算法搜索准确的对称轴位置和方向,如Powell搜索算法。
考虑人体解剖结构具有对称性的特点,这种对称不仅体现在几何结构上,而且在图像的灰度上也具有一定的对称性,如图2所示,直线为对称轴,圆形、三角形、正方形分别为关于对称轴对称的三组对称点,假设灰度值分别为1、2、3,如果图像完全对称,那么每组对称点灰度值之差等于零,整幅图像对称点集的灰度差的绝对值之和等于零。实际图像中对称的两个像素点的灰度值并不一定相等,而是相近的,因此可以在图像中移动或者旋转图像的对称轴,每移动一个位置,计算整幅图像关于对称轴对称的所有像素点的灰度差的总和,当这个总和达到最小值时,此位置的对称轴即为目标区域准确的作为定位信息的对称轴。
本发明计算出目标区域的主轴的角度α0和重心(x,y)作为对称轴初始值,把目标区域对称轴两侧对称点集灰度差之和的表达式作为代价函数;然后在一定范围内改变重心点(x,y)坐标和对称轴方向三个参数,使用最优化搜索算法搜索代价函数的最小值,当搜索完成时得到的即是所求取的对称轴,此所求取的对称轴即为MRI定位图像的定位线。
此方法同时应用于脑横断面、颈椎冠状面、腰椎冠状面以及其他具备解剖结构对称性部位的定位,效果如图3~图5所示,分别为脑部横断面定位效果图、颈椎冠状面定位效果图以及腰椎冠状面定位效果图。
Claims (3)
1.一种基于人体解剖结构对称性的磁共振成像扫描自动定位方法,其特征在于包含如下步骤:
a)获取MRI定位图像;
b)利用人体解剖结构的对称性,在MRI定位图像上自动计算不同方位的定位信息;
c)输出MRI定位图像的定位信息;
所述步骤b)具体为:
1)在整幅图像中分割出目标区域;
2)初始化上述目标区域的对称轴;
3)计算上述目标区域中关于上述对称轴对称的点集;
4)基于上述点集定义代价函数,通过优化处理得到上述代价函数最优值;
5)得到与上述代价函数最优值对应的作为定位信息的对称轴。
2.按权利要求1所述的基于人体解剖结构对称性的磁共振成像扫描自动定位方法,其特征在于:所述步骤2)初始化上述目标区域的对称轴是基于计算上述目标区域的主轴进行的。
3.按权利要求1所述的基于人体解剖结构对称性的磁共振成像扫描自动定位方法,其特征在于:所述步骤4)代价函数定义为上述目标区域中对称点灰度差的绝对值之和。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200610046121A CN100586368C (zh) | 2006-03-20 | 2006-03-20 | 基于人体解剖结构对称性的磁共振成像扫描自动定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200610046121A CN100586368C (zh) | 2006-03-20 | 2006-03-20 | 基于人体解剖结构对称性的磁共振成像扫描自动定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1843295A CN1843295A (zh) | 2006-10-11 |
CN100586368C true CN100586368C (zh) | 2010-02-03 |
Family
ID=37062285
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200610046121A Active CN100586368C (zh) | 2006-03-20 | 2006-03-20 | 基于人体解剖结构对称性的磁共振成像扫描自动定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100586368C (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011160309A1 (zh) * | 2010-06-25 | 2011-12-29 | 中国科学院自动化研究所 | 基于鲁棒统计信息传播的多模态三维磁共振图像脑肿瘤分割方法 |
US11666241B2 (en) | 2019-02-02 | 2023-06-06 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for medical imaging |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103142230B (zh) * | 2013-03-27 | 2015-05-27 | 深圳安科高技术股份有限公司 | 一种磁共振定位扫描的方法及装置 |
DE102013210855A1 (de) | 2013-06-11 | 2014-12-11 | Siemens Aktiengesellschaft | Verfahren zu einer Anpassung einer Schichtpositionierung innerhalb eines Schichtprotokolls für eine Magnetresonanzuntersuchung sowie eine Magnetresonanzvorrichtung zur Ausführung des Verfahrens |
-
2006
- 2006-03-20 CN CN200610046121A patent/CN100586368C/zh active Active
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011160309A1 (zh) * | 2010-06-25 | 2011-12-29 | 中国科学院自动化研究所 | 基于鲁棒统计信息传播的多模态三维磁共振图像脑肿瘤分割方法 |
US9129382B2 (en) | 2010-06-25 | 2015-09-08 | Institute Of Automation, Chinese Academy Of Sciences | Method and system for brain tumor segmentation in multi-parameter 3D MR images via robust statistic information propagation |
US11666241B2 (en) | 2019-02-02 | 2023-06-06 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for medical imaging |
Also Published As
Publication number | Publication date |
---|---|
CN1843295A (zh) | 2006-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110946654B (zh) | 一种基于多模影像融合的骨科手术导航系统 | |
AU2019432052B2 (en) | Three-dimensional image measurement method, electronic device, storage medium, and program product | |
CN110264504B (zh) | 一种用于增强现实的三维配准方法和系统 | |
EP1892668B1 (en) | Registration of imaging data | |
CN110101452A (zh) | 一种用于外科手术的光磁一体化定位导航方法 | |
CN100418489C (zh) | 手术导航中基于基准面膜的多模式医学图像配准系统 | |
CN111260786A (zh) | 一种智能超声多模态导航系统及方法 | |
CN100586368C (zh) | 基于人体解剖结构对称性的磁共振成像扫描自动定位方法 | |
CN1326092C (zh) | 手术导航中基于基准面膜的多模式医学图像配准方法 | |
CN114305685B (zh) | 一种用于髋关节置换手术中的髋骨配准方法 | |
CN106137395B (zh) | 应用于无标记点光学手术导航系统的全自动病人注册方法 | |
US20130156281A1 (en) | Image registration method | |
CN109934861B (zh) | 一种头颈部多模态医学图像自动配准方法 | |
CN109771850B (zh) | 一种自适应放疗计划修正装置 | |
KR20210115010A (ko) | 이미지 처리 방법 및 장치, 전자 기기 | |
Bächler et al. | Restricted surface matching—numerical optimization and technical evaluation | |
US6757417B2 (en) | Method and apparatus for defining a three-dimensional imaging section | |
CN103729827A (zh) | 影像增强器c形臂x射线机三维重建重叠伪影校正方法 | |
CN1732848A (zh) | 对患者进行磁共振检查的方法 | |
CN113850810B (zh) | 用于转正股骨的方法及手术系统、存储介质以及电子设备 | |
US9254106B2 (en) | Method for completing a medical image data set | |
CN109003281A (zh) | 一种基于三维空间模型获得单块脊椎骨的装置及方法 | |
US20110163750A1 (en) | Method and apparatus for correcting the uniformity of a magnetic field | |
CA3088626C (en) | Method, system and computer program for determining position and/or orientation parameters of an anatomical structure | |
Song et al. | Computer-aided modeling and morphological analysis of hip joint |
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 | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Patentee after: DongSoft Medical System Co., Ltd. Address before: Hunnan New Century Road 110179 Shenyang city of Liaoning Province, No. 16 Patentee before: Dongruan Medical Systems Co., Ltd., Shenyang |