CN106473744B - 一种基于平面切割的产轴线计算方法 - Google Patents
一种基于平面切割的产轴线计算方法 Download PDFInfo
- Publication number
- CN106473744B CN106473744B CN201610831020.2A CN201610831020A CN106473744B CN 106473744 B CN106473744 B CN 106473744B CN 201610831020 A CN201610831020 A CN 201610831020A CN 106473744 B CN106473744 B CN 106473744B
- Authority
- CN
- China
- Prior art keywords
- pelvis
- axis
- plane
- cutting
- planes
- 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
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 17
- 238000004364 calculation method Methods 0.000 title claims abstract description 10
- 210000004197 pelvis Anatomy 0.000 claims abstract description 71
- 230000002146 bilateral effect Effects 0.000 claims abstract description 10
- 239000000203 mixture Substances 0.000 claims abstract description 8
- 238000000034 method Methods 0.000 claims description 15
- 238000010586 diagram Methods 0.000 claims description 12
- 210000003754 fetus Anatomy 0.000 claims description 5
- 230000033001 locomotion Effects 0.000 claims description 4
- 230000002269 spontaneous effect Effects 0.000 claims description 4
- 230000001154 acute effect Effects 0.000 claims description 3
- 210000000988 bone and bone Anatomy 0.000 claims description 3
- 230000001605 fetal effect Effects 0.000 abstract description 17
- 238000005259 measurement Methods 0.000 abstract description 9
- 238000005516 engineering process Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 230000035606 childbirth Effects 0.000 description 3
- 210000001215 vagina Anatomy 0.000 description 3
- 238000000691 measurement method Methods 0.000 description 2
- 208000036029 Uterine contractions during pregnancy Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 238000012730 cephalometry Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000004700 fetal blood Anatomy 0.000 description 1
- 210000002458 fetal heart Anatomy 0.000 description 1
- 238000003706 image smoothing Methods 0.000 description 1
- 208000015181 infectious disease Diseases 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 210000003689 pubic bone Anatomy 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 210000001693 sacrococcygeal joint Anatomy 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
- A61B5/1075—Measuring physical dimensions, e.g. size of the entire body or parts thereof for measuring dimensions by non-invasive methods, e.g. for determining thickness of tissue layer
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/43—Detecting, measuring or recording for evaluating the reproductive systems
- A61B5/4306—Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
- A61B5/4343—Pregnancy and labour monitoring, e.g. for labour onset detection
- A61B5/4362—Assessing foetal parameters
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Gynecology & Obstetrics (AREA)
- Pediatric Medicine (AREA)
- Pregnancy & Childbirth (AREA)
- Reproductive Health (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
本发明公开了一种基于平面切割的产轴线计算方法,包括以下步骤:1)数据预处理;以骨盆的医学影像数据作为输入,进行预处理提取网格化的骨盆数据,获得骨盆的三角网格模型,调整模型位于标准解剖学位置,即盆骨左右对称部位处于同一水平面且双边髂前上棘与耻骨结节组成的平面平行于z轴;2)使用一系列平面截取处于标准解剖学位置的骨盆的三角网格模型,这些平面与骨盆轴垂直;3)在得到的一系列截平面上求最大空圆,得到一系列圆心坐标,其连线及延长线组成的产轴线。本发明为胎头位置数字化和精确化测量提供了一种简便、准确的表示手段,能够更好地与现代产程监护电子设备结合,避免了目前医院中普遍采用的手动测量的低效和低精度的缺点。
Description
技术领域
本发明涉及骨盆测量和产程监护技术领域,尤其是指一种基于平面切割的产轴线计算方法。
背景技术
随着电子技术和信息技术的发展,各种胎儿监护设备越来越普及,如多普勒胎心仪、胎儿监护仪、脐带血流检测仪等,而监护参数包括胎心率、宫缩压力、胎动、母亲血压、血氧饱和度、脉博、心电等。但是,目前所有的电子胎儿监护技术主要针对产前监护,很少有针对分娩过程(产程)的监护,而后者则是影响出生人口质量的关键因素,对新生儿的智力和健康有关键作用。另外,通过实时的产程监护,能够及时发现胎头位置异常情况,为医生对患者实行剖宫产还是自然分娩的决策提供依据,避免因决策不及时造成的各类母儿近远期并发症甚至死亡。产程图是分娩学中描述产程的直观依据,是决定分娩方式和控制分娩进程的主要依据。胎头位置是产程图的重要参数。胎头位置的测量是以骨盆轴为参考的。骨盆轴是指连接骨盆各假想平面中点的曲线。此轴上段向下向后,中段向下,下段向下向前。分娩时,胎儿沿此轴娩出。在解剖学中,骨盆轴通常根据骨盆的特征点得到,包括骨盆入口横径、坐骨棘平面位置、耻骨联合上缘、骶尾关节等,如图1所示。但由于胎头的特殊解剖位置以及测量时的特殊要求,长期以来的临床实践中,一般将骨盆轴简化为一条至上而下的直线,如图2所示。由于没有切实可行的宫颈扩张和胎头位置自动测量方法,产程进展判断仍主要依赖传统阴道指检,这种方式存在以下问题:
1、高主观性,客观性差,高度依赖医生经验。
2、胎头位置估测不准确:间接估测胎头位置,胎头变形和形成可能造成误导。
3、胎位估测不准:30%-60%阴道检查不能准确测定胎位。
4、反复阴道检查,给产妇带来不适、增加感染的几率。
发明内容
本发明的目的在于克服现有技术的缺点和不足,提供一种简便、准确的基于平面切割的产轴线计算方法,避免了目前医院中普遍采用的手动测量的低效和低精度的缺点,有利于骨盆测量和产程监护的数字化、精确化。
为实现上述目的,本发明所提供的技术方案为:一种基于平面切割的产轴线计算方法,该产轴线定义为自然分娩过程中胎儿头部中心的运动轨迹,该方法包括以下步骤:
1)数据预处理
以骨盆的医学影像数据作为输入,进行预处理提取网格化的骨盆数据,获得骨盆的三角网格模型,调整模型位于标准解剖学位置,即盆骨左右对称部位处于同一水平面且双边髂前上棘与耻骨结节组成的平面平行于z轴;
2)使用一系列平面截取处于标准解剖学位置的骨盆的三角网格模型,这些平面与骨盆轴垂直;
3)在得到的一系列截平面上求最大空圆,得到一系列圆心坐标,其连线及延长线组成的产轴线;其中,产轴线在骨盆以外的部分被称作延长线,该延长线由尾椎的弯曲趋势决定。
在步骤1)中,骨盆的基准位置定义如下:
对于骨盆点云数据,调整骨盆方向,使得骨盆从初始的标准解剖学位置绕x轴沿顺时针方向旋转90°,使得髂前上棘与耻骨结节处于同一水平面;其中,所述标准解剖学位置是指盆骨左右对称部位处于同一水平面且双边髂前上棘与耻骨结节组成的平面平行于z轴,其冠状面和矢状面视图旋转后得到的初始骨盆位置。
在步骤2)中,采用保持平面不变,而不断旋转骨盆的方法实现从垂直于骨盆轴的方向对三维骨盆数据进行切割,其具体实现方式如下:
2.1)对于网格化骨盆数据,设切割平面α:zα=k,其与骨盆相交得到的截平面内的最大空圆圆心为C1;
2.2)另设切割平面β:zβ=zα+j,其与骨盆相交得到的截平面内的最大空圆圆心为C2;
2.3)求向量C1C2与z轴的所夹的锐角θ,若直线C1C2的斜率为正,则使骨盆绕过点C2的x轴平行线逆时针旋转θ;若直线C1C2的斜率为负,则使骨盆绕过点C2的x轴平行线顺时针旋转θ,旋转骨盆后,向量C1C2平行于z轴;
2.4)转到步骤2.1),在新的位置上对骨盆进行切割,直到切割平面遍历整个骨盆;
在上述步骤中,k、j为常数,k取值根据不同骨盆模型对象,通常选取相对靠底部的位置。
在步骤3)中,求截平面的最大空圆是使用基于Voronoi图的算法求解。
本发明与现有技术相比,具有如下优点与有益效果:
本发明方法采用一系列垂直于骨盆轴平面对骨盆各段进行切割并在相应的截平面上求最大空圆圆心,这些圆心的连线及其延长线即是所求的产轴线。这种产轴线自动计算方法,为胎头位置数字化和精确化测量提供了一种简便、准确的表示手段,能够更好地与现代产程监护电子设备结合,避免了目前医院中普遍采用的手动测量的低效和低精度的缺点。
附图说明
图1为基于特征点的骨盆轴解剖示意图。
图2为传统的胎头位置测量方法示意图。
图3为本发明方法的流程示意图。
图4a为标准解剖学位置冠状面视图。
图4b为标准解剖学位置矢状面视图。
图4c为初始骨盆位置调整示意图。
图5a为直线C1C2斜率为正的结构示意图。
图5b为C2的x轴平行线逆时针旋转θ示意图。
图5c为直线C1C2斜率为负的结构示意图。
图5d为C2的x轴平行线顺时针旋转θ示意图。
具体实施方式
下面结合具体实施例对本发明做进一步的说明。
本实施例所述的基于平面切割的产轴线计算方法,产轴线定义为自然分娩过程中胎儿头部中心的运动轨迹。该方法为胎头测量提供了一种新的概念,为胎头位置数字化和精确化测量提供了一种简便、准确的表示手段,弥补了传统胎头位置测量方法无法与现代产程监护电子设备结合的不足。该自动计算方法使用一系列垂直于骨盆轴平面对骨盆各段进行切割并在相应的截平面上求最大空圆圆心,这些圆心的连线及其延长线即是所求的产轴线。如图3所示,包括以下步骤:
1)数据预处理
对某患者的CT扫描进行图像平滑、图像增强处理后用阈值分割的方法提取轮廓线并进行三维重构和三角剖分,得到骨盆的三角网格模型。调整模型方向,使得骨盆从初始的标准解剖学位置绕x轴沿顺时针方向旋转90°,得到的初始骨盆位置。其中,骨盆的基准位置定义为:对于骨盆点云数据,调整骨盆方向,使得骨盆从初始的标准解剖学位置绕x轴沿顺时针方向旋转90°,使得髂前上棘与耻骨结节处于同一水平面;其中,所述标准解剖学位置是指盆骨左右对称部位处于同一水平面且双边髂前上棘与耻骨结节组成的平面平行于z轴,其冠状面和矢状面视图旋转后得到的初始骨盆位置,如图4a至4c所示。
2)采用保持平面不变,通过旋转骨盆的方法实现从垂直于骨盆轴的方向对三维骨盆数据进行切割。每次骨盆旋转角度两次平面切割所得两个截平面内的最大空圆圆心决定,其具体方法如下:
2.1)取平面α:zα=k。则平面与骨盆数据相交后截平面的所有点坐标集合为P1={p11,p12.....p1n}。此截平面最大内接圆圆心为C1。
2.2)保持骨盆不动,取平面β:zβ=zα+0.5。则平面与骨盆数据相交后截平面的所有点坐标集合为P2={p2l,p22.....p2n}。此截平面最大内接圆圆心为C2。
2.3)求向量C1C2与z轴的所夹的锐角θ。若直线C1C2的斜率为正(如图5a所示),则使骨盆绕过点C2的x轴平行线逆时针旋转θ(如图5b所示);若直线C1C2的斜率为负(如图5c所示),则使骨盆绕过点C2的x轴平行线顺时针旋转θ(如图5d所示)。旋转骨盆后,向量C1C2平行于z轴。
2.4)k=k+0.5。转到步骤2.1),在新的位置上对骨盆进行切割,直到切割平面遍历整个骨盆。
考虑到骨盆模型在z轴上坐标范围是[-2.8,3],k的初始值取-2.5。
3)对于截平面点集Pn={P1,P2},通过基于Voronoi图求截平面内的最大空圆圆心。算法具体描述如下:
3.1)计算点集Pn的Voronoi图,记为V(Pn)。
3.2)计算点集Pn的凸壳,记为CH(Pn);令rmax=0。
3.3)对于V(Pn)的每个Voronoi顶点v,如果v在CH(Pn)内部,则计算以v为圆心的空圆的半径并修改rmax。
3.4)对于每条Voronoi边,计算CH(Pn)的边与Voronoi边的交点s,计算圆心在s的空圆的半径并修改rmax。
3.5)输出rmax和对应的圆心。
以上所述实施例只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。
Claims (3)
1.一种基于平面切割的产轴线计算方法,该产轴线定义为自然分娩过程中胎儿头部中心的运动轨迹,其特征在于,该方法包括以下步骤:
1)数据预处理
以骨盆的医学影像数据作为输入,进行预处理提取网格化的骨盆数据,获得骨盆的三角网格模型,调整模型位于标准解剖学位置,即盆骨左右对称部位处于同一水平面且双边髂前上棘与耻骨结节组成的平面平行于z轴;
2)使用一系列平面截取处于标准解剖学位置的骨盆的三角网格模型,这些平面与骨盆轴垂直;
3)在得到的一系列截平面上求最大空圆,得到一系列圆心坐标,其连线及延长线组成的产轴线;其中,产轴线在骨盆以外的部分被称作延长线,该延长线由尾椎的弯曲趋势决定。
2.根据权利要求1所述的一种基于平面切割的产轴线计算方法,其特征在于,在步骤2)中,采用保持平面不变,而不断旋转骨盆的方法实现从垂直于骨盆轴的方向对三维骨盆数据进行切割,其具体实现方式如下:
2.1)对于网格化骨盆数据,设切割平面α:zα=k,其与骨盆相交得到的截平面内的最大空圆圆心为C1;
2.2)另设切割平面β:zβ=zα+j,其与骨盆相交得到的截平面内的最大空圆圆心为C2;
2.3)求向量C1C2与z轴的所夹的锐角θ,若直线C1C2的斜率为正,则使骨盆绕过点C2的x轴平行线逆时针旋转θ;若直线C1C2的斜率为负,则使骨盆绕过点C2的x轴平行线顺时针旋转θ,旋转骨盆后,向量C1C2平行于z轴;
2.4)转到步骤2.1),在新的位置上对骨盆进行切割,直到切割平面遍历整个骨盆;
在上述步骤中,k、j为常数,k取值根据不同骨盆模型对象,选取相对靠底部的位置。
3.根据权利要求1所述的一种基于平面切割的产轴线计算方法,其特征在于:在步骤3)中,求截平面的最大空圆是使用基于Voronoi图的算法求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610831020.2A CN106473744B (zh) | 2016-09-19 | 2016-09-19 | 一种基于平面切割的产轴线计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610831020.2A CN106473744B (zh) | 2016-09-19 | 2016-09-19 | 一种基于平面切割的产轴线计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106473744A CN106473744A (zh) | 2017-03-08 |
CN106473744B true CN106473744B (zh) | 2019-08-20 |
Family
ID=58267484
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610831020.2A Expired - Fee Related CN106473744B (zh) | 2016-09-19 | 2016-09-19 | 一种基于平面切割的产轴线计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106473744B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110151222A (zh) * | 2018-02-14 | 2019-08-23 | 王幼萍 | 产程导航曲线生成方法、胎头位置测算方法和装置 |
CN110507362B (zh) * | 2019-09-09 | 2023-06-23 | 邹金沙 | 一种产轴线构建方法、胎头位置测量方法及装置 |
CN112348862B (zh) * | 2020-11-03 | 2022-11-22 | 广州三瑞医疗器械有限公司 | 一种基于骨盆模型的骨盆平面中心测算方法及其应用 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7850625B2 (en) * | 2003-08-06 | 2010-12-14 | Trig Medical Ltd. | Method and apparatus for monitoring labor parameter |
CN202342032U (zh) * | 2011-09-02 | 2012-07-25 | 暨南大学 | 基于磁场跟踪定位技术的电子产程图系统 |
CN102314712B (zh) * | 2011-09-23 | 2013-11-06 | 暨南大学 | 基于三维成像技术的骨盆建模与骨盆参数测算方法 |
WO2014069553A1 (ja) * | 2012-10-31 | 2014-05-08 | 株式会社上島電興社 | 骨盤回旋角度算出装置及び骨盤回旋角度算出方法、臼蓋被覆3次元形状推定装置及び臼蓋被覆3次元形状推定方法、骨盤回旋角度算出プログラム、並びに臼蓋被覆3次元形状推定プログラム |
US9247998B2 (en) * | 2013-03-15 | 2016-02-02 | Intellijoint Surgical Inc. | System and method for intra-operative leg position measurement |
US10478151B2 (en) * | 2014-08-12 | 2019-11-19 | General Electric Company | System and method for automated monitoring of fetal head descent during labor |
-
2016
- 2016-09-19 CN CN201610831020.2A patent/CN106473744B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN106473744A (zh) | 2017-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110464380B (zh) | 一种对中晚孕期胎儿的超声切面图像进行质量控制的方法 | |
CN106473744B (zh) | 一种基于平面切割的产轴线计算方法 | |
CN106963378A (zh) | 一种基于电磁定位与超声成像的胎头方位测量方法 | |
CN106999157A (zh) | 用于通过计算和显示超声测量和图形模型进行胎儿显像的方法与系统 | |
JP2008534082A (ja) | 胎児、新生児及び成体の器官のオペレータに依存しない標準超音波画像を取得及び生成するシステム、方法及び媒体 | |
CN110279433A (zh) | 一种基于卷积神经网络的胎儿头围自动精确测量方法 | |
CN106570856A (zh) | 水平集动态规划颈总动脉内中膜厚度测量装置和方法 | |
CN115222937A (zh) | 一种脊柱侧弯检测方法及装置 | |
WO2022089051A1 (zh) | 颅骨矫正方案生成系统、构建方法、获取方法及装置 | |
CN108836394B (zh) | 一种胎头下降角度自动测量方法 | |
CN106875377A (zh) | 一种基于高斯曲率流的脊椎特征点自动识别方法 | |
Zhen et al. | Locating multiple standard planes in first-trimester ultrasound videos via the detection and scoring of key anatomical structures | |
CN202342032U (zh) | 基于磁场跟踪定位技术的电子产程图系统 | |
Santoso et al. | Health monitoring of fetal ultrasound image using active contour models | |
WO2020103098A1 (zh) | 超声成像方法、设备、存储介质,处理器及计算机设备 | |
Smith et al. | Automated torso contour extraction from clinical cardiac MR slices for 3D torso reconstruction | |
Chervenak et al. | Current status of fetal age and growth assessment | |
CN112348862B (zh) | 一种基于骨盆模型的骨盆平面中心测算方法及其应用 | |
CN113143306B (zh) | 一种用于颅内动脉瘤cfd诊断可视化设备 | |
Lu et al. | Computerised interpretation systems for cardiotocography for both home and hospital uses | |
CN112638267B (zh) | 超声成像方法及系统、存储介质、处理器和计算机设备 | |
CN116211353B (zh) | 可穿戴超声膀胱容量测定与多模态影像形貌评估系统 | |
KR20200117896A (ko) | 태아 신경계의 질환을 결정하기 위한 시스템 및 방법 | |
Guo et al. | Artificial Intelligence in Prenatal Ultrasound: Clinical Application and Prospect. | |
CN110648739B (zh) | 一种基于空间解剖位置的阑尾ct图像报告方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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: 20190820 |