CN102933156A - 确定子宫准备分娩的情况的方法和系统 - Google Patents

确定子宫准备分娩的情况的方法和系统 Download PDF

Info

Publication number
CN102933156A
CN102933156A CN2011800287646A CN201180028764A CN102933156A CN 102933156 A CN102933156 A CN 102933156A CN 2011800287646 A CN2011800287646 A CN 2011800287646A CN 201180028764 A CN201180028764 A CN 201180028764A CN 102933156 A CN102933156 A CN 102933156A
Authority
CN
China
Prior art keywords
displacement
uterus
stress
parameter
alpha
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
CN2011800287646A
Other languages
English (en)
Other versions
CN102933156B (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.)
Cornell University
Original Assignee
Cornell University
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 Cornell University filed Critical Cornell University
Publication of CN102933156A publication Critical patent/CN102933156A/zh
Application granted granted Critical
Publication of CN102933156B publication Critical patent/CN102933156B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61FFILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F2/00Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
    • A61F2/82Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
    • A61F2/86Stents in a form characterised by the wire-like elements; Stents in the form characterised by a net-like or mesh-like structure
    • A61F2/90Stents in a form characterised by the wire-like elements; Stents in the form characterised by a net-like or mesh-like structure characterised by a net-like or mesh-like structure
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0833Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
    • A61B8/085Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/485Diagnostic techniques involving measuring strain or elastic properties
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61FFILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F2/00Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
    • A61F2/95Instruments specially adapted for placement or removal of stents or stent-grafts
    • A61F2/962Instruments specially adapted for placement or removal of stents or stent-grafts having an outer sleeve
    • DTEXTILES; PAPER
    • D03WEAVING
    • D03DWOVEN FABRICS; METHODS OF WEAVING; LOOMS
    • D03D1/00Woven fabrics designed to make specified articles
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • G06T7/0014Biomedical image inspection using an image reference approach
    • G06T7/0016Biomedical image inspection using an image reference approach involving temporal comparison
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Pathology (AREA)
  • Vascular Medicine (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Transplantation (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Cardiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Quality & Reliability (AREA)
  • Physiology (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Textile Engineering (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Media Introduction/Drainage Providing Device (AREA)

Abstract

一种确定子宫准备分娩的情况的方法,该方法包括使用至少部分子宫壁在收缩前后的影像获得组织位移信息,并使用组织位移信息获得的参量来确定子宫准备分娩的情况。在一个例子中,该方法还包括从组织位移信息获得造成组织位移的应力,在这个例子中,由组织位移信息获得的参量包括由应力获得的参量。本文还公开了实施该方法和制品的系统。

Description

确定子宫准备分娩的情况的方法和系统
背景技术
这些教导主要涉及用于检测波所诱发的组织变形的基于成像技术相关的成像方式。
在怀孕期间早产是母亲和胎儿发病及死亡的首要原因。妇产科可以使用药物预防早产,但这种治疗如若延迟则无效,并且治疗过度会对母子有害。怀孕期间的另一个问题是由于子宫收缩无效导致分娩没有进展。同样可以使用药物治疗以防需要剖腹产,但对无效宫缩延迟和过度治疗可能与治疗早产有同样的风险。因而,诊断早产和无效宫缩是产科医生面临的最重要的挑战。
目前诊断分娩和无效宫缩的方法有局限性。将宫腔压力导管插入阴道内。分娩力计、胎儿纤维连接蛋白检测和阴道检查是主观的。另一种技术是在腹部表面使用肌电图,该技术现在备受关注,是目前产科很有前途的一项技术。研究表明有几个由肌电图获得的试验变量,如电脉冲的频率,可用于诊断分娩情况。然而,对于分娩期和非分娩期的女性,这些变量值的范围不一致,投入临床使用前需要进一步研究。
已经提出在怀孕期间使用组织弹性成像来评估宫颈的成熟程度,组织弹性成像是一种评估组织硬度的无创技术。组织弹性成像技术由J.Ophir开发,已经用于医学,因为许多病症有机械性变化的表现。例如,癌变组织被按压时可能比正常组织僵硬,肝脏组织的剪切模量可能因纤维化而变大。总体上,当给定本构关系和边界条件后,可以用应变来计算应力和材料参数。弹性成像技术已被应用于肿瘤、肌肉骨骼、心血管疾病和其他应用,并且被证明是有效的。
需要一种诊断分娩情况并区分有效和无效子宫收缩的改进的方法和系统。
发明内容
在一个实施方式中,用于确定子宫准备分娩的情况的这些教导的方法包括从至少部分子宫壁在收缩前后的影像获得组织位移信息,以及利用由组织位移信息获得的参量来确定子宫准备分娩的情况。在一个实施方式中,该方法还包括由组织位移信息获得造成组织位移的应力,并且在这个例子中,由组织位移信息获得的参量包括由应力获得的参量。
还公开了用于实施这些教导的方法的系统的实施方式,以及其中包含计算机可读代码的计算机可用介质的实施方式,所述计算机可读代码的特征在于能驱动处理器执行这些教导的方法的各个步骤。
为更好地理解本文的教导以及其他进一步的目的,请参考附图和详细说明,其范围将由所附的权利要求确定。
附图简述
图1a、1b和1c是表示这些教导的方法的实施方式的示意性流程图;和
图2是表示这些教导的系统的实施方式的示意性方框图。
具体实施方式
下面详细描述的是实施这些教导的最佳设想模式。这些说明并无限制性意义,仅用于对这些教导的一般原理做解释说明,这些教导的范围由所附的权利要求做最佳的界定。
本文使用的“确定子宫准备分娩的情况”是指区分属于非分娩性子宫生理活动的收缩和导致分娩的子宫收缩。
本文所使用的“平衡关系”是指确保满足守恒率的情况。在一个例子中,这种关系要求应力张量发散基本为零,这确保当体力几乎为零时线性动量守恒。
本文所使用的“肌节”是由横纹肌的原纤维分割而成的其中一个部分。
本文所使用的“材料硬度”是指应力关于应变的导数。在一个实施方式中,用于确定子宫准备分娩的情况(区分子宫内的有效收缩和无效收缩)的这些教导的方法包括由至少部分子宫壁在收缩前后的影像获得组织位移信息,并利用由组织位移信息获得的参量来确定子宫准备分娩的情况。
在另一个实施方式中,这些教导的用于确定子宫准备分娩的情况的方法(区分子宫有效收缩和无效收缩)包括获得至少部分子宫壁在收缩前后的影像、组织位移信息,由组织移位信息获得造成组织位移的应力,并用所得到的应力参量确定子宫准备分娩的情况。
在一个例子中,应力参量包括子宫壁的压力。在一个实施方式中,获得应力的步骤包括使用平衡关系的步骤。
在一个例子中,至少部分子宫壁在收缩前后的影像是超声影像。
这些教导的系统的一个实施方式包括一个或多个处理器,以及其上具有计算机可读代码的计算机可用介质,该计算机可读代码驱动一个或多个处理器执行这些教导的方法。
图1a显示表示这些教导的方法的一个实施方式的流程图。参考图1a,获得至少部分子宫壁在收缩前后的影像(步骤110,图1a)。由收缩前后的影像获得组织位移信息(步骤120,图1a)。本文所使用的组织位移是指收缩过程中子宫壁本质上的各个点偏离收缩前位置的运动。本文所使用的组织位移应区别于涉及两个在不同点位上的传感器所测定的位移。两个在不同点位上的传感器所测定的位移不能提供有用的应变信息。用组织位移信息获得的参量确定子宫准备分娩的情况(或获得可以帮助医师做出判断的信息)(步骤130,图1a)。
图1b显示表示这些教导的方法的另一实施方式的流程图。参考图1b,获得至少部分子宫壁在收缩前后的影像(步骤110,图1b)。由收缩前后的影像获得组织位移信息(步骤120,图1b)。获得组织位移信息后,能获得造成组织位移的应力(步骤125,图1b)。用组织位移信息获得的参量确定子宫准备分娩的情况(或获得可以帮助医师做出判断的信息)(步骤130,图1b)。
图1c显示了表示这些教导的方法的实施方式的整合流程图。
图2是表示这些教导的系统的实施方式的方框图。图2所示的这些教导的系统的实施方式包括一个或多个处理器220,以及一个或多个具有计算机可读代码的计算机可用介质230,计算机可读代码使一个或多个处理器可由至少部分子宫壁在收缩前后的图像210获得组织位移信息,在一个例子中,由获得的组织位移信息获得应力,并使用由组织位移信息(在一个例子中,还有应力)获得的参量来确定子宫准备分娩的情况。输出接口240可输出这些教导的方法的结果。使用连接组件235(例如,但不限于计算机总线)可操作地连接图像接收接口210、一个或多个处理器220、输出接口240和计算机可用介质230。在一个例子中,还可操作地连接有数据库(未显示)。数据库可以包括确定子宫准备分娩的情况的数据或预报器。
通过下列示例性实施方式来解释这些实施方式和这些教导的原理。
为阐明这些教导,提供以下信息。
应变估计值
使用一般的弹性成像技术,通过在组织变形过程中捕获一系列超声数据,计算软组织承受变形时的应变。在各时间点和组织内感兴趣的各点,通过在变形的组织内定位与未变形组织中的感兴趣的点周围区域最超声相似的区域,估计组织位移。采用几个指标衡量各区域间的相似程度。最常用的指标例如有交叉相关性、方差和、协方差、绝对差总和、这些方法的各种归一化版本、混合信号相关性、极性重合相关性和零相位。Viola等得出结论,虽然他们的研究中不包括零相位,但归一化交叉相关性是最佳算法(F.Viola,W.F.Walker,A comparison of the performance of time-delay estimators in medicalultrasound,.Dept.of Biomed.Eng.,vol.50,no.4,pp.392.401,April 2003,其全部内容通过引用并入本文以供所有目的使用)。对于各个指标,通过在变形的组织中定位使该指标最大或最小的区域来估计所感兴趣的点的位移,或者在零相位的情况下,通过定位该指标为零的区域来估计所感兴趣的点的位移。
未变形区域与变形区域的超声结果几乎从不完全相似的原因有四个。首先,组织可能会移出传感器2-1)的运动平面。其次,应变变化导致一个区域的超声波外观改变。最后,第三条和第四条原因是抬升方向和横向上的超声波数据的采样低于奈奎斯特率,并且超声波有随机噪声。为了提高位移估计的准确度和分辨率,同时保持可接受的计算效率,研究人员一直在寻求改变弹性成像技术,以解决这四个问题。
最初的一种迭代弹性成像算法由Konofagou和Ophir引入(E.Konfagou,J.Ophir,″ANew Elastographic Method For Estimation and Imaging of Lateral Displacements,Lateral Strains,Corrected Axial Strains and Poisson’s Ratios in Tissues,″Ultrasound in Medicine and Biology,vol.24,no.8,pp.1183-1199,October 1998,其通过引用整体并入本文以供所有目的使用);还可见于公布的美国专利US20080097202,US20070049824,其通过引用整体并入本文以供所有目的使用,以及出版物Ophir J,AlamSK,Garra B,Kallel F,Konofagou E,Krouskop T,Varghese T,Elastography:ultrasonic estimation and imaging of the elastic properties of tissues,Proc Inst Mech Eng H.,1999;213(3):203-33;Lee W-N.and Konofagou E.E.,Angle-Independentand Multi-Dimensional Myocardial Elastography:From Theory to ClinicalValidation,Ultrasonics,48(6-7):563-7,2008(Invited),其通过引用整体并入本文供所有目的使用)。从本质上改变窗口大小和搜索区域,该算法首先沿轴向搜索最佳匹配,然后由轴向位移估计横向搜索,然后由最初的横向和轴向位移估计再次进行轴向搜索。他们使用播种式位移评估获得更精确的位移估计值,这正是我们的弹性成像算法的一个重要构思。
比较超声影像时,超声波图像——即超声斑点——去相关,造成位移估计错误。造成斑点去相关的原因有四个:平面外运动,局部应变,在垂直于波动方向的方向上的低分辨率和随机噪音。研究人员修改跟踪算法解决这四个问题,并改善位移估计的精度,同时保持可接受的计算效率。
在这些教导的方法的一个实施方式中,使用了混合法(L.Chen,G.M.Treece,J.E.Lindop,A.H.Gee,R.W.Prager,(2009),13(2),pp.286-296.B.Garra,E.Cespedes,J.Ophir,S.Spratt,R.Zuurbier和C.M.CM,″A quality-guided displacementtracking algorithm for ultrasonic elasticity imaging,″Medical Image Analysis,vol.13,no.2,pp.286-296,2009,其通过引用整体并入本文以供所有目的使用)。混合方法结合了四种跟踪算法:多重网格法(H.Chen,H.Shi,T.Varghese,.Improvementof displacement estimation using a two-step cross-correlation method,UltrasoundMed Biol,vol.33,no.1,pp.48.56,January 2007其通过引用整体并入本文供所有目的使用)、质量导引(L.Chen,G.M.Treece,J.E.Lindop,A.H.Gee,R.W.Prager,A quality-guided displacement tracking algorithm for ultrasonic elasticityimaging,.Medical Image Analysis,vol.13,no.2,pp.286.296,2009,其通过引用整体并入本文供所有目的使用)、零相位(A.Pesavento,C.Perrey,M.Krueger,H.Ermert.A time-efficient and accurate strain estimation concept for ultrasonicelastography using iterative phase zero estimation.Ultrasonics,Ferroelectricsand Frequency Control,IEEE Transactions on,vol.46,no.5,pp.1057-1067,1999,其通过引用整体并入本文以供所有目的使用)、和交叉相关性(J.Ophir,E.Cespedes,H.Ponnekanti,Y.Yazdi,and X.Li,.Elastography:A quantitative method for imagingthe elasticity of biological materials,Ultrasonic Imaging,vol.13,no.2,pp.111.134,April 1991,通过引用整体并入本文以供所有目的使用;还可参见公开的美国专利US20080097202,US20070049824,其通过引用整体并入本文以供所有目的使用),分为三级计算。在每一级,在不同的帧之间比较矩形窗口中的超声斑点。将位移估计值作为一帧中的窗口与和其最相似的另一帧中的那个窗口之间的滞后给出。在前两级计算中使用二维归一化交叉相关性计算斑点相似度,在最后一级计算中用零相位计算斑点相似度。进行等距点跟踪,但在每一级后续计算中细化间距和缩小窗口尺寸。混合方法的质量引导方面使用跟踪点的位移来初始化尚未被更精确跟踪的点初始化的相邻点的位移。
然后可以使用中值滤波法(M.M.Doyley,J.C.Bamber,F.Fuechsel,N.L.Bush.Afreehand elastographic imaging approach for clinical breast imaging:systemdevelopment and performance evaluation.Ultrasound in Medicine & Biology,vol.27,no.10,pp.1347-1357,2001,其通过引用整体并入本文以供所有目的使用)和去噪算法(A.Chambolle,An algorithm for total variation minimization and applications.JMath Imaging Vis,vol.20,pp.89-97,2004;X.Bresson,T.Chan.Fast dualminimization of the vectorial total variation norm and applications to color imageprocessing.CAM Report 07-25,2007,两者通过引用整体并入本文以供所有目的使用)来修平估计的位移。使用线性最小二乘法由位移估计应变(F.Kallel and J.Ophir,.Aleast-squares strain estimator for elastography,Ultrasonic Imaging,vol.19,no.3,pp.195-208,July 1997,其通过引用整体并入本文以供所有目的使用)。
结构方程
准确的结构方程,以及合适的几何因素和边界条件,使由应变评估值准确估算应力成为可能。有许多用于软组织的结构模型,但对于这些结构模型中的大多数而言,常见的有下列四个条件:软组织(几乎)不可压缩、有粘性、有超弹性,以及粘性和弹性应力项是附加项。Y.C.Fung的生物力学:活组织的力学性能(Y.C.Fung,Biomechanics:MechanicalProperties of Living Tissues,2nd Edition.New York:Springer-Verlag,1993,242-314)可为这些条件提供更详细的解释。尤其是肌肉,其行为与大多数软组织相似,但还可产生收缩力。有收缩性的组成部分可与被动组成部分平行建模(CITE)。依据这些条件得到的应力通用方程为四项之和:被动和主动组成部分的粘性和弹性项。
在此示例性的实施方式中,Pioletti等及Veronda和Westmann的现象逻辑模型分别用于被动组成部分的粘度和弹性项(D.P.Pioletti,L.R.Rakotomanana,J.-F.Benvenuti,P.-F.Leyvraz,.Viscoelastic constitutive law in large deformations:application to human knee ligaments and tendons,.Journal of Biomechanics,vol.31,no.8,pp.753.757,August 1998;D.R.Veronda,R.A.Westmann,Mechanicalcharacterization of skin-finite deformations,.Journal of Biomechanics,vol.3,no.1,pp.111.124,January 1970,两者通过引用整体并入本文以供所有目的使用)。在Veronda和Westmann的弹性项中有一个未知的罚函数,与(几乎)不可压缩性条件相关。在此示例性的实施方式中,使用由Sainte-Marie等在他们的心脏和另一个软组织模型中提出的罚函数(J.Sainte-Marie,D.Chapellea,R.Cimrmanc,M.Sorinea,Modelingand estimation of the cardiac electromechanical activity,Computers Structures,vol.84,no.28,pp.1743.1759,November 2006,其通过引用整体并入本文以供所有目的使用)。主动组成部分使用广义的超弹性模型。其对应力的贡献程度通过平衡方程获得。
几何模型
子宫几何模型使用这些教导有三个条件,这些条件是常见的用于数学壳理论条件,见SS Antman的弹性的非线性问题(S.S.Antman,Nonlinear Problems of Elasticiy.NewYork:Springer-Verlag,1995,353-383)中所描述的。根据这篇关于壳理论的论文,使用数学框架来导出下面将使用的对估算应力很重要的变形张量的表达式。此外,由于子宫壁与子宫壁到子宫腔中心的距离相比而言较小,因而在本示例性实施方式的技术描述中所解释的壳样条件非常准确。
捕捉超声波数据
将线性或矩阵阵列超声传感器放在孕妇的腹部。子宫的几何形状可以近似成关于对称轴对称,对称轴通过子宫腔中心几乎从远端延伸到近端。调整传感器,使得超声波A线阵列能线性通过这个对称轴。这种调整大致与将压电元件阵列放在腹部表面从而排列在矢状面内相对应。一旦传感器被固定在适当的位置上,需要用超声波标记对称轴。通过使用子宫两侧壁的超声影像比较其几何形状来定位此轴。每次子宫收缩时,捕捉超声射频数据序列,进行弹性成像分析结果。将获得的第一个超声射频数据序列的第一个矩阵阵列作为参考帧,与每个序列的首个矩阵阵列进行对比。进行这种区分比较是因为在每次子宫收缩前子宫内都可能有残留应变。
应变估计
在此示例性的实施方式中,使用前述混合方法评估子宫收缩期间的子宫壁位移。使用分割技术勾画子宫壁的外形。如上所述,用线性最小二乘法由位移确定应变。
子宫收缩建模
对此示例性实施方式中使用的子宫的几何形状和变形介绍如下。三个主要条件如下:子宫是轴对称的,变形是轴对称的,并且子宫的活动方式类似于壳体。将这些条件定义如下:
轴对称的子宫有对称轴,被称做
Figure BDA00002562010800071
从而在子宫绕此轴旋转任何角度后,子宫的几何形状与初始的几何形状相比无法区分。
当子宫发生轴对称变形时,在绕
Figure BDA00002562010800081
轴固定角度的平面内的物质仍然保持在该平面内,该平面内可见的变形与对应于任何其他固定角度的平面的内变形相比无法区分。
当最初沿子宫外法线的分布的材料在变形后仍保持在一直线上时,子宫的活动类似于壳体。
前两个条件允许在子宫收缩过程中使用二维传感器测量子宫的整体变形情况。第三个条件如此命名是因为其在数学的壳理论中被使用。当壳体的厚度比壳体的其他长度尺寸小很多时,能适当描述几何形状和变形情况。此条件是准确的,因为此误差与弹性成像技术变形评估的误差数量级相同。因而,可以假设“壳样”条件。
子宫模型的结构方程如下所述。下列条件用于结构方程:
子宫主动组成部分导致的应力可与被动组成部分导致的应力平行建模。
在子宫肌肉的各个点,平滑肌细胞(SMC)的长轴方向与子宫切平面平行,并与对称轴
Figure BDA00002562010800082
在同一平面内。
SMC内的肌节沿SMC的长轴分布。
子宫几乎不可压缩。
第一个条件用在许多肌肉的结构模型中。最大优势是使用弹性成像技术计算的变形能用于独立计算两个并行分支中的应力。对于被动矩阵,使用由Pioletti提出的粘弹性模型(D.P.Pioletti,L.R.Rakotomanana,J.-F.Benvenuti,P.-F.Leyvraz,.Viscoelasticconstitutive law in large deformations:application to human knee ligaments andtendons,.Journal of Biomechanics,vol.31,no.8,pp.753.757,August 1998,其通过引用整体并入本文以供所有目的使用)和Veronda-Westmann弹性模型(D.R.Veronda,R.A.Westmann,.Mechanical characterization of skin-finite deformations,Journalof Biomechanics,vol.3,no.1,pp.111-124,January 1970,其通过引用整体并入本文以供所有目的使用)。子宫被认为是几乎不可压缩的,而非不可压缩(J.Sainte-Marie,D.Chapellea,R.Cimrmanc,M.Sorinea,Modeling and estimation of the cardiacelectromechanical activity,Computers Structures,vol.84,no.28,pp.1743-1759,November 2006,其通过引用整体并入本文以供所有目的使用),这允许不使用边界条件而直接计算变形的被动组成部分中的应力,以发现不可压缩材料所用方程的拉格朗日系数。
主动组成部分导致的应力包含方程中的未知因素,因为计算应力时需要已知动作电位、钙的浓度和任何对抗剂。如下所示,可使用平衡方程发现此未知因素。一旦获知未知数值进而是已经获得的实际应力值,就可以计算子宫壁上的压力。并且还可以依据应力确定材料的硬度。
几何形状和变形
符号:小写粗体字指示
Figure BDA00002562010800091
的向量,大写粗体字指示
Figure BDA00002562010800092
的张量。带上标的向量为单位法向,qp
Figure BDA00002562010800093
的简写。
将子宫视为是对称的。换言之,存在一个轴,用
Figure BDA00002562010800094
指代,从而子宫绕轴旋转后的几何形状与原几何形状难以区分。定义和
Figure BDA00002562010800095
Figure BDA00002562010800096
,使
Figure BDA00002562010800097
符合右手规则。使用
Figure BDA00002562010800098
的三曲线坐标系来定义子宫上的每个点。与
Figure BDA00002562010800099
轴的夹角使用轴测定,显示为y是沿某个基准表面外法线的距离,s是沿着同一基准表面的距离。点
Figure BDA000025620108000912
关于于
Figure BDA000025620108000913
坐标的参考配置X由以下公式给出:
Figure BDA000025620108000914
其中b0是子宫表面的单元外法线,是从原点指到基准表面的向量。
轴对称变形意味着子宫内的任何物质点在子宫变形过程中不会围绕
Figure BDA000025620108000916
旋转。用
Figure BDA000025620108000917
作为单位向量,其指向在径向上远离
Figure BDA000025620108000918
并定位成与
Figure BDA000025620108000919
的夹角是
Figure BDA000025620108000920
因而,所有
Figure BDA000025620108000921
平面内的点变形后仍在该平面内。此外,此平面内的变形对于每个
Figure BDA000025620108000922
都相等。
“壳样”条件表明变形前沿外法线的子宫组织在变形后仍分布在一直线上。因而,此条件及变形是对称性变形的条件允许将变形分解成三个部分:基准表面的变形
Figure BDA000025620108000923
向外法线在
Figure BDA000025620108000924
平面θ(s)内旋转的角度,以及沿外法线的变形δ(s,r,y)。因此变形的一般形式可写成如下公式:
Figure BDA000025620108000925
其中
Figure BDA000025620108000926
是变形后初始外法线的单元方向,
Figure BDA000025620108000927
由下式给出:
Figure BDA00002562010800101
Figure BDA00002562010800102
此外,
Figure BDA00002562010800103
Figure BDA00002562010800104
Figure BDA00002562010800105
Figure BDA00002562010800106
Figure BDA00002562010800107
Figure BDA00002562010800108
注意有三个附加的坐标系统:
Figure BDA00002562010800109
Figure BDA000025620108001010
此处
Figure BDA000025620108001011
Figure BDA000025620108001012
分别是最初垂直于子宫表面的单位方向和此向量在变形之后的单位方向。
变形梯度F可通过使用链式法则给出:
Figure BDA000025620108001013
其中
Figure BDA000025620108001015
的逆变基础。
还有
X , s = ( 1 - θ 0 , s y ) a ^ 0
X , y = b ^ 0
Figure BDA000025620108001018
x , s = ( u , s cos θ + w , s sin θ - δ , s ) a ^ + ( - u , s sin θ + w , s cos θ - δ , s ) b ^
x , y = δ , y b ^
注意有三个附加的坐标系统:
Figure BDA00002562010800113
Figure BDA00002562010800114
此处
Figure BDA00002562010800115
Figure BDA00002562010800116
分别是最初垂直于子宫表面的单位方向和此向量在变形之后的单位方向。
变形梯度F如下:
Figure BDA00002562010800117
其中
Figure BDA00002562010800118
的逆变基础。使用上述定义,
X , s = ( 1 - θ 0 , s y ) a ^ 0
X , y = b ^ 0
Figure BDA000025620108001112
x , s = ( u , s cos θ + w , s sin θ - δ , s ) a ^ + ( - u , s sin θ + w , s cos θ - δ , s ) b ^
Figure BDA000025620108001114
x , y = δ , y b ^
由于,
Figure BDA000025620108001116
Figure BDA000025620108001117
Figure BDA000025620108001118
然后
g s = 1 1 - θ 0 , s y a ^ 0 ;
Figure BDA000025620108001120
g y = b ^ 0 .
因此变形梯度F变为:
Figure BDA00002562010800121
为了简化表达式,我们已经设定:
α=u,scosθ+w,ssinθ-δθ,s
α0=1-θ0,s
β=-u,ssinθ+w,scosθ-δ,s
γ=δ,y
χ=u-δsinθ
χ0=u0-ysinθ.
方程式的左和右柯西-格林张量(Cauchy-Green Tensor)B和C按照上述表达为:
Figure BDA00002562010800123
SMC的几何形状和位置条件暗示其长轴方向在变形前后分别是
Figure BDA00002562010800124
Figure BDA00002562010800125
因此,
Figure BDA00002562010800126
是沿SMC长轴的拉伸,γ是沿外法线方向的拉伸,而
Figure BDA00002562010800127
是环绕
Figure BDA00002562010800128
轴测得的环状拉伸。
结构模型
下面使用的结构模型考虑了因子宫中的收缩因素造成的子宫肌肉承受的总应力的应力总和,以及由有收缩因素嵌入其中的被动矩阵所产成的应力。此条件相当于说明主动组成部分和被动组成部分平行分布。该条件的直接结果是各分支的变形是整个子宫肌肉的变形,依据该变形,各平行分支中的应力被解偶。
被动矩阵
在时刻t时的被动矩阵被分解成两个附加组成部分:
σp(t)=σe(t)+σv(t)
方程式右侧的两项分别是与变形对应的弹性(e)和粘性(v)。
本示例的实施方式中,弹性响应利用Veronda和Westmann的弹性应变能(D.R.Veronda,R.A.Westmann,.Mechanical characterization of skin-finite deformations,Journalof Biomechanics,vol.3,no.1,pp.111-124,January 1970,其通过引用整体并入本文以供所有目的使用)。然而,由于组织被认为是(几乎)不可压缩的,因而应变能被分为体积部分(volumetric part)和等容部分:
W e = c 1 exp ( c 2 ( I ‾ 1 - 3 ) ) + c 3 ( I ‾ 2 - 3 ) + g ( J )
其中
J ≡ det ( F ) = αγχ α 0 χ 0
I ‾ 1 ≡ J - 2 3 trace ( C ) = ( α 2 + β 2 α 0 2 + γ 2 + χ 2 χ 0 2 ) ( α 0 χ 0 αγχ ) 2 / 3
I ‾ 2 ≡ 1 2 J - 4 / 3 ( I 1 2 - trace ( C · C ) ) = ( α 2 γ 2 α 0 2 + χ 2 χ 0 2 ( α 2 + β 2 α 0 2 + γ 2 ) ) ( α 0 χ 0 αγχ ) 4 / 3
此处,g是强加几乎不可压缩行为的罚函数。大多数软组织的分析模型均假定不可压缩性,即I3=1。相应的弹性应力σe涉及未知的拉格朗日乘数p,为了使用边界条件而求解该拉格朗日乘数p。数字上看,假定组织接近不可压缩常常比较容易,从而g(I3)是|I3-1|>0的罚函数。如同Sainte-Marie等的著作中(Sainte-Marie)所述,g(x)=c4(x-1-ln(x))
其中c4>>1。通过使用罚函数,无需知道任何边界条件即用B可明确求解σe
σ e = 2 I 3 - 1 ∂ W e ∂ B B = 2 I 3 - 1 ( ∂ W e ∂ I 1 + ∂ W e ∂ I 2 ( I 1 1 - B ) + ∂ W e ∂ I 3 B - 1 ) B
= 2 I 3 - 1 ( ( c 1 c 2 exp ( c 2 ( I 1 - 3 ) ) - c 3 I 1 ) B - c 3 B 2 + c 4 ( 1 - I 3 - 1 ) 1 )
请注意
J , C = 1 2 JC - 1
I ‾ 1 , C = ( J - 2 3 trace ( C ) ) , C = J - 2 / 3 1 - 2 3 J - 5 / 3 trace ( C ) J , C = J - 2 / 3 ( 1 - trace ( C ) 3 C - 1 )
I ‾ 2 , C = 1 2 ( J - 4 / 3 trace ( C ) 2 - J - 4 / 3 trace ( C · C ) ) , C
= J - 4 3 ( trace ( C ) 1 - C ) - 4 3 J - 7 3 ( trace ( C ) 2 - trace ( C · C ) ) J , C
= J - 4 3 ( ( trace ( C ) 1 - C ) - 2 3 ( trace ( C ) 2 - trace ( C · C ) ) C - 1 ) ,
和,
W e , J = c 4 ( 1 - 1 J ) ; W e , I ‾ 1 = c 1 c 2 exp ( c 2 ( I ‾ 1 - 3 ) ) ; W e , I ‾ 2 = c 3 .
求解弹性应力得出:
σ e = J - 1 F · ( 2 c 1 c 2 exp ( c 2 ( I ‾ 1 - 3 ) ) J - 2 3 ( 1 - trace ( C ) 3 C - 1 ) . . .
+ c 3 J - 4 3 ( ( trace ( C ) 1 - C ) - 2 3 ( trace ( C ) 2 - trace ( C · C ) ) C - 1 ) . . .
+ c 4 ( 1 - 1 x ) JC - 1 ) · F T
= 2 c 1 c 2 exp ( c 2 ( I ‾ 1 - 3 ) ) J 5 / 3 ( B - trace ( C ) 3 1 )
+ c 3 J 7 / 3 ( ( trace ( C ) B - B · B ) - 2 3 ( trace ( C ) 2 - trace ( C · C ) ) 1 )
+ c 4 ( 1 - 1 J ) 1
Pioletti等人研究的粘性表达式等包括较少的项目。为简化表达,仅考虑了上述表达式的第一项:
σ v = 2 I 3 - 1 ∂ W v ∂ B B = 2 I 3 - 1 ( ∂ W v ∂ J 1 ∂ J 1 ∂ B ) B
其中
W v = η 4 J 1 2 ( I 1 - 3 )
J 1 = trace ( dB dt )
而η是严格的正参数。此伪应变能函数满足Clausius-Duhem不等式,因为
Figure BDA00002562010800154
计算σv得出:
σ v = η ( I ‾ 1 - 3 ) J 1 2 J B
活性组成部分
在子宫肌肉中的各点,平滑肌细胞(SMC)的长轴,进而是肌节的长轴被认为是位于由子宫表面的切线和对称轴所确定的一个平面内,且垂直于该切线。变形后,“壳样”条件暗示SMC的长轴和肌节位于
Figure BDA00002562010800156
的方向上。
平滑肌细胞内的各肌节的动态模型,如Hai和Murphy建立的四态模型,依赖细丝和粗丝之间的相对滑动(速率v,和位移d,)以及其他状态变量,如钙浓度、拮抗剂的浓度和动作电位。这些状态变量定义为向量q(t)。相对滑动被认定是肌节和q方向上的拉伸和变形速度的函数,然后σa=αa(q,I4,J2),其中
Figure BDA00002562010800157
Figure BDA00002562010800158
认为主动组成部分是超弹性的:
σ a = 2 J - 1 F · ( Ψ a , C + Ψ a , C ) · F T = 2 J - 1 F · ( Ψ a , I 4 I 4 , C + Ψ a , J 2 J 2 , C ) · F T
= 2 J - 1 F · ( Ψ e , I 4 a ^ 0 α ^ 0 + Ψ e , I 4 a ^ 0 α ^ 0 ) · F T
= 2 J - 1 ( Ψ e , I 4 + Ψ e , J 2 ) ( α α 0 a ^ + β α 0 b ^ ) ( α α 0 a ^ + β α 0 b ^ )
= ξ ( α 2 a ^ a ^ + αβ b ^ b ^ + β 2 b ^ b ^ )
其中
ξ = 2 ( Ψ e , I 4 + Ψ e , J 2 ) J α 0 2
在上述方程中,最后的参量ξ未知,可使用平衡方程求解。
整个结构方程
结合主动和被动组成部分的方程,忽略长期历史效应,子宫内应力的最终方程是:
σ = σ e + σ v + σ a
= 2 c 1 c 2 exp ( c 2 ( I ‾ 1 - 3 ) ) J 5 / 3 ( B - trace ( C ) 3 1 )
+ c 3 J 7 / 3 ( ( trace ( C ) B - B · B ) - 2 3 ( trace ( C ) 2 - trace ( C · C ) ) 1 )
+ c 4 ( 1 - 1 J ) 1 + η ( I ‾ 1 - 3 ) J 1 2 J B + ξ ( α 2 + a ^ a ^ + αβ ( a ^ b ^ + b ^ a ^ ) + β 2 b ^ b ^ )
= d 1 a ^ a ^ + d 2 ( a ^ b ^ + b ^ a ^ ) + d 3 b ^ b ^ + d 4 e ^ 2 e ^ 2 + ξ ( α 2 a ^ a ^ + αβ ( a ^ b ^ + b ^ a ^ ) + β 2 b ^ b ^ )
其中
a 1 ≡ ( α 0 χ 0 αγχ ) 5 3 2 c 1 c 2 exp ( c 2 ( ( α 2 + β 2 α 0 2 + γ 2 + χ 2 χ 0 2 ) ( α 0 χ 0 αγχ ) 2 3 1 - 3 ) )
a 2 ≡ c 4 ( 1 - α 0 χ 0 αγχ )
a 3 ≡ α 0 χ 0 η 2 αγχ ( ( α 2 + β 2 α 0 2 + γ 2 + χ 2 χ 0 2 ) ( α 0 χ 0 αγχ ) 2 3 1 - 3 )
d 1 ≡ a 1 3 ( 2 α 2 - β 2 α 0 2 - γ 2 - χ 2 χ 0 2 ) + a 2 3 ( α 2 γ 2 α 0 2 + ( α 2 - 2 β 2 ) χ 2 α 0 2 χ 0 2 - 2 γ 2 χ 2 χ 0 2 )
+ c 3 ( α 0 χ 0 αγχ ) 7 / 3 + a 3 α 2 α 0 2
d 2 ≡ a 1 ( αβ α 0 2 - 1 3 ( α 2 + β 2 α 0 2 + γ 2 + χ 2 χ 0 2 ) )
+ a 2 ( αβ χ 2 α 0 2 χ 0 2 - 2 3 ( α 2 γ 2 α 0 2 + χ 2 χ 0 2 ( α 2 + β 2 α 0 2 + γ 2 ) ) )
d 3 ≡ a 1 3 ( 2 β 3 - α 2 α 0 2 + 2 γ 2 - χ 2 χ 0 2 ) + a 2 3 ( ( β 2 - 2 α 2 ) χ 2 α 0 2 χ 0 2 + γ 2 χ 2 χ 0 2 - 2 α 2 γ 2 α 0 2 )
+ c 3 ( α 0 χ 0 αγχ ) 7 / 3 + a 4 ( β 2 α 0 2 + γ 2 )
d 4 ≡ a 1 3 ( 2 α 2 - β 2 α 0 2 - γ 2 - χ 2 χ 0 2 ) + a 2 3 ( χ 2 χ 0 2 ( α 2 + β 2 α 0 2 ) + γ 2 χ 2 χ 0 2 - 2 α 2 γ 2 α 0 2 )
+ c 3 ( α 0 χ 0 αγχ ) 7 / 3 + a 3 χ 2 χ 0 2
平衡方程
若给定
Figure BDA00002562010800179
此时认定体力为零并且惯性条件可以忽略不计,则满足线性动量守恒。由于
Figure BDA000025620108001710
Figure BDA000025620108001711
的标准正交基,则平衡方程是
( ▿ · σ ) · a ^ = 0
( ▿ · σ ) · b ^ = 0
( ▿ · σ ) · e ^ 2 = 0
下面的标示用于任何向量c。
注意
∂ x ∂ s = α a ^ + α b ^
Figure BDA00002562010800181
∂ x ∂ y = y b ^
并且
a ^ , s = θ , s b ^ ; b ^ , s = - θ , s a ^
Figure BDA00002562010800185
Figure BDA00002562010800186
e ^ 2 , s = a ^ , y = b ^ , y = e ^ 2 , y = 0
因此对于任何标量函数l
▿ · ( l a ^ ) = l ( θ , s β - sin θ χ ) + l , s α
▿ · ( l b ^ ) = l ( θ , s α - cos θ χ ) + l , s β + l , y γ
Figure BDA000025620108001810
因此,线性动量的平衡方程是:
( d 1 θ , s + d 2 , s β + d 2 θ , s + d 1 , s α - d 1 sin θ + d 2 cos θ χ + d 2 , y γ )
= ξ ( - α 2 θ , s + ( αβ ) , s β - β θ , s - 2 α , s + α 2 sin θ + αβ cos θ χ - ( αβ ) , y γ ) - 2 ξ , s α ξ , y αβ γ
( d 2 θ , s + d 3 , s β + d 3 θ , s + d 2 , s α - d 2 sin θ + d 3 cos θ χ + d 3 , y γ )
= ξ ( - α θ , s - 2 β , s - β 2 θ , s + ( αβ ) , s α + αβ sin θ + β 2 cos θ χ - 2 ββ , y γ ) - 2 ξ , s β - ξ , y β 2 γ
Figure BDA000025620108001815
满足第三个方程,因为在轴对称条件下
Figure BDA00002562010800191
前两个方程可以结合起来,消除对ξ,s和ξ,y的依赖:
( ( d 1 - d 3 ) θ , s + ( α 2 - β 2 ) d 2 θ , s + d 1 , s β 2 - d 3 , s α 2 αβ )
+ - β ( d 1 sin θ + d 2 cos θ ) + α ( d 2 sin θ + d 3 cos θ ) χ + β d 2 , y - α d 3 , y γ )
= ξ ( 2 ( α β , s - α , s β ) + β ( α β , y - α , y β ) γ )
因此,可以由平衡条件解决未知变量ξ。
参数
为获得准确的应力,提供准确的η,c1,c2,c3,和c4值非常重要。前四个参数为材料参数,它们取决于子宫肌肉的性质。可以从软组织的相关文献获得这些材料参数的初始值。然而,通过机械测试被动子宫肌肉可以获得精确值。此外,如果参数因患者而异,对每个患者进行超声弹性成像检测可以获得该患者的估计值。可通过准静态地或借助机械振动而向被动材料施加变形,并通过迭代反演法估算参数。
为简化计算并保持数值稳定,引入最后一个参数c4。该参数的引入允许将受限的问题变成不受限的问题。认为c4>>1,从而获得近似解。
体内应用
在任何缩之前捕获的超声波RF数据允许获得初始动力学参量和向量,它们就是上文所标注的下标为零的参数。然后,在收缩过程中确定位移以后,确定动力学参量α,β,χ,和γ以及向量
Figure BDA00002562010800195
Figure BDA00002562010800196
使用上文描述的方法获得ξ,σ,和子宫壁压力的大小:
| p | = | n ^ σ n ^ |
其中在子宫壁处测得应力,
Figure BDA00002562010800198
是变形的子宫壁的外法线。
从测得数值获得参数与子宫准备分娩的情况的相关性。其他诊断技术依赖于电活动、宫腔内压力和材料硬度来预测分娩,本教导获得的变量得出一些相同或相关的参量。变量ξ与电活动相关,压力ν是宫腔内压力,而材料硬度涉及应力σ关于动力学参量的导函数。还计算了上述参量以外的参量,从而能够更好地进行相关性计算。
为了由参量获得最准确的预测,还要对患者进行检测,并记录分娩的结果。可以使用机器学习技术来确定能预测子宫准备分娩情况的参量值。可能的机器学习方法可以是训练Random ForestsTM(Breiman L.Random forests,Machine Learning,vol.45,pp.5-32,2001,其通过引用整体并入本文)来跟据测量值预测准备情况,除了Random ForestsTM,还有其他机器学习技术可以应用,例如,但不限于,支持向量机(Corinna Cortes and V.Vapnik,″Support-Vector Networks″,Machine Learning,20,1995,其通过引用整体并入本文)。各种其他机器学习技术也在这些教导的范围内。利用患者的检测值和由应力组织位移信息所获得的参量之间的关系,可确定子宫准备分娩或早产的情况(或为医师做决定提供信息)。
由于上述示例性实施方式是为澄清和解释这些教导,因而可能有多种其他实施方式。上述示例性实施方式中描述的某些可能变化包括以下内容,但这些教导不仅限于这些变化。
1.无简化假设的变形替代公式
2.被动组成部分的替代性结构方程。示例包括但不限于:
a.不同的粘度方程
b.包含纤维方向的横观各向同性结构方程
c.粘度应力和无附加的弹性行为
d.不同的弹性方程
3.主动组成部分的替代性结构方程。示例包括但不限于:
a.肌节方向的不同分布,无论其是离散的还是连续的
b.主动组成部分的应力是除
Figure BDA00002562010800201
Figure BDA00002562010800202
以外的其他动力学变量的函数
4.主动和被动组成部分不平行的替代性结构方程
5.使用平衡方程获得与收缩强度相关的信息。可能需要为了跟上述结构方程不同的结构方程而算出数字PDE。
6.用其他技术,如肌电图、分娩力计或外界振动,来补充超声弹性成像技术,以便估算子宫收缩过程中的应力或材料参数。
7.替代弹性成像算法。示例包括但不限于:
a.基于超声多普勒的方法
b.包括惩罚因子以增强平滑性的方法
c.使用小波的方法
d.包括全局或局部拉伸斑点的方法
8.将传感器放在多个位点,以便获得子宫或宫颈的RF数据
9.利用子宫收缩过程中通过超声弹性成像技术获得的参量来减少妊娠期出现并发症的可能性。
应当指出的是虽然上述揭示的示例性实施方式使用超声影像,但是这些教导不限于示例性实施方式。
为了描述和定义本教导的目的,应当指出此处使用的术语“基本上”来表现固有程度的不确定性,这种不确定性可归因于任何参量对比、数值、检测值或其他代表性数据。本文中还用术语“基本上”表现在不导致所讨论主题的基本功能变化的情况下,参量表现可能不同于规定的参照的程度。
为了实施相同的功能,可能进一步将本文所述的元件和组成部分划分为附加的组成部分或合并在一起而形成更少的组成部分。
每个计算机程序都可以使用任何编程语言来编写,如汇编语言、机器语言、高层次的过程编程语言或面向对象的编程语言。编程语言可以是编译的或解释的编程语言。
每个计算机程序可以在具体表现为计算机可读存储设备的计算机程序产品中通过计算机处理器执行来实施。可以计算机处理器执行在计算机可读介质上具体表现的程序来执行本发明的方法步骤,即通过进行输入和产生结果输出。
计算机可读介质的通常形式包括,例如,软盘,软磁盘,硬盘,磁带,或任何其他磁介质,CD-ROM,任何其他光学介质,任何具有孔图案的物理介质,RAM,PROM和EPROM,FLASH-EPROM,任何其他存储器芯片或卡带。如美国专利商标局2005年对专利申请的专利主题审查的临时指导原则中所述,1300 Off.Gaz.Pat.Office 142(Nov.22,2005)。“另一方面,从技术角度来看,由功能性描述素材编码的信号类似于由功能性描述素材编码的计算机可读存储器,它们都与计算机创建了功能性交互关系。换言之,计算机能执行编码功能,无论其格式是磁盘或是信号。”
虽然已参照多个实施方式对教导进行了描述,但应该认识到在所附权利要求限定的本发明的主旨和范围的情况下,这些教导还能提供多种的进一步的和其他的实施方式。

Claims (17)

1.一种确定子宫准备分娩的情况的方法,该方法包括:
由至少部分子宫壁在收缩前后的影像获得组织位移信息;以及
利用由所述组织位移信息获得的参量来确定子宫准备分娩的情况。
2.如权利要求1所述的方法,其特征在于,进一步包括:
由所述组织位移信息获得导致组织位移的应力;以及
由所述组织位移信息获得的参量包括由应力获得的参量。
3.如权利要求2所述的方法,其特征在于,由所述应力获得的参量包括子宫壁处的压力。
4.如权利要求2所述的方法,其特征在于,由所述应力获得的参量包括材料的硬度。
5.如权利要求2所述的方法,其特征在于,获得应力的步骤包括使用平衡关系的步骤。
6.如权利要求1所述的方法,其特征在于,所述至少部分子宫壁在收缩前后的影像是超声影像。
7.一种用于确定子宫准备分娩的情况的系统,该系统包括:
至少一个处理器;和
至少一种计算机可用介质,所述计算机可用介质中具有计算机可读代码,所述计算机可读代码能使所述至少一个处理器
由至少部分子宫壁在收缩前后的影像获得组织位移信息;以及
利用由组织位移信息获得的参量来确定子宫准备分娩的情况。
8.如权利要求7所述的系统,其特征在于,所述计算机可读代码还能使所述至少一个处理器
由所述组织位移信息获得导致组织位移的应力;以及
其中由所述组织位移信息获得的参量包括由应力获得的参量。
9.如权利要求8所述的系统,其特征在于,所述由应力获得的参量包括子宫壁处的压力。
10.如权利要求8所述的系统,其特征在于,所述由应力获得的参量包括材料的硬度。
11.如权利要求8所述的系统,其特征在于,在使所述至少一个处理器获得导致组织位移的应力时,所述计算机可读代码使所述至少一个处理器利用平衡关系。
12.一种计算机程序产品,包括具有计算机可读代码的计算机可用介质,所述计算机可读代码用于使计算机系统:
由至少部分子宫壁在收缩前后的影像获得组织位移信息;以及
利用由组织位移信息获得的参量来确定子宫准备分娩的情况。
13.如权利要求12所述的计算机程序产品,其特征在于,所述计算机可读代码可进一步使计算机系统:
由所述组织位移信息获得导致组织位移的应力;以及
所述由组织位移信息获得的参量包括由应力获得的参量。
14.如权利要求13所述的计算机程序产品,其特征在于,在使所述至少一个处理器获得导致组织位移的应力时,所述计算机可读代码使所述至少一个处理器使用平衡关系。
15.如权利要求13所述的计算机程序产品,其特征在于,所述由应力获得的参量包括子宫壁处的压力。
16.如权利要求13所述的计算机程序产品,其特征在于,所述由应力获得的参量包括材料的硬度。
17.如权利要求12所述的计算机程序产品,其特征在于,所述影像是超声影像。
CN201180028764.6A 2010-04-13 2011-04-11 确定子宫准备分娩的情况的方法和系统 Expired - Fee Related CN102933156B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US32341310P 2010-04-13 2010-04-13
US61/323,413 2010-04-13
PCT/US2011/031889 WO2011130143A1 (en) 2010-04-13 2011-04-11 Method and systems for determining preparedness of the uterus for delivery

Publications (2)

Publication Number Publication Date
CN102933156A true CN102933156A (zh) 2013-02-13
CN102933156B CN102933156B (zh) 2015-07-15

Family

ID=44798976

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201180028764.6A Expired - Fee Related CN102933156B (zh) 2010-04-13 2011-04-11 确定子宫准备分娩的情况的方法和系统

Country Status (3)

Country Link
US (1) US20130218289A1 (zh)
CN (1) CN102933156B (zh)
WO (1) WO2011130143A1 (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140007885A1 (en) * 2012-07-06 2014-01-09 Edward M. Gillis Systems and methods for treatment of sleep apnea
JP2013526392A (ja) 2010-05-21 2013-06-24 レベント メディカル インコーポレイテッド 睡眠時無呼吸を治療するためのシステム及び方法
US8628554B2 (en) 2010-06-13 2014-01-14 Virender K. Sharma Intragastric device for treating obesity
US10010439B2 (en) 2010-06-13 2018-07-03 Synerz Medical, Inc. Intragastric device for treating obesity
US10420665B2 (en) 2010-06-13 2019-09-24 W. L. Gore & Associates, Inc. Intragastric device for treating obesity
US9526648B2 (en) 2010-06-13 2016-12-27 Synerz Medical, Inc. Intragastric device for treating obesity
WO2012018545A2 (en) 2010-07-26 2012-02-09 Revent Medical, Inc. Systems and methods for treatment of sleep apnea
US11026826B2 (en) * 2016-01-29 2021-06-08 Gi Dynamics, Inc. Flanged gastrointestinal devices and methods of use thereof
US10779980B2 (en) 2016-04-27 2020-09-22 Synerz Medical, Inc. Intragastric device for treating obesity

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6066095A (en) * 1998-05-13 2000-05-23 Duke University Ultrasound methods, systems, and computer program products for determining movement of biological tissues
US6270459B1 (en) * 1998-05-26 2001-08-07 The Board Of Regents Of The University Of Texas System Method for estimating and imaging of transverse displacements, transverse strains and strain ratios
EP1161921A2 (en) * 2000-06-07 2001-12-12 GE Marquette Medical Systems, Inc. Uterine contraction detector and frequency trender
US20030114779A1 (en) * 1997-05-05 2003-06-19 Yoav Paltieli Method and apparatus for monitoring the progress of labor
US6749571B2 (en) * 2002-09-19 2004-06-15 Wisconsin Alumni Research Foundation Method and apparatus for cardiac elastography
CN1846233A (zh) * 2003-08-06 2006-10-11 特里格医疗有限公司 监测分娩参数的方法和装置
US20070015994A1 (en) * 2005-07-14 2007-01-18 Hyundae Hong In-vivo measurement of biomechanical properties of internal tissues
US20070167819A1 (en) * 2005-12-06 2007-07-19 Osborn Thomas W Iii Method for in-vivo measurement of biomechanical properties of internal tissues
CN101049245A (zh) * 2007-04-30 2007-10-10 陆尧胜 一种产程监护装置及其产程监护方法
US7297116B2 (en) * 2003-04-21 2007-11-20 Wisconsin Alumni Research Foundation Method and apparatus for imaging the cervix and uterine wall
US20080167581A1 (en) * 2007-01-10 2008-07-10 Yoav Paltieli Determining parameters associated with a female pelvis and cervix
CN101454792A (zh) * 2006-05-24 2009-06-10 梁伦锡 分娩控制系统及控制方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6149595A (en) * 1998-07-02 2000-11-21 Seitz; Walter S. Noninvasive apparatus and method for the determination of cardiac valve function
JP4269623B2 (ja) * 2002-10-07 2009-05-27 株式会社 東北テクノアーチ 血流可視化診断装置
US20070016391A1 (en) * 2005-07-14 2007-01-18 Ryo Minoguchi Reverse finite element analysis and modeling of biomechanical properties of internal tissues
WO2008073491A1 (en) * 2006-12-11 2008-06-19 University Of Florida Research Foundation, Inc. System and method for analyzing progress of labor and preterm labor

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030114779A1 (en) * 1997-05-05 2003-06-19 Yoav Paltieli Method and apparatus for monitoring the progress of labor
US6669653B2 (en) * 1997-05-05 2003-12-30 Trig Medical Ltd. Method and apparatus for monitoring the progress of labor
US6066095A (en) * 1998-05-13 2000-05-23 Duke University Ultrasound methods, systems, and computer program products for determining movement of biological tissues
US6270459B1 (en) * 1998-05-26 2001-08-07 The Board Of Regents Of The University Of Texas System Method for estimating and imaging of transverse displacements, transverse strains and strain ratios
EP1161921A2 (en) * 2000-06-07 2001-12-12 GE Marquette Medical Systems, Inc. Uterine contraction detector and frequency trender
US6749571B2 (en) * 2002-09-19 2004-06-15 Wisconsin Alumni Research Foundation Method and apparatus for cardiac elastography
US7297116B2 (en) * 2003-04-21 2007-11-20 Wisconsin Alumni Research Foundation Method and apparatus for imaging the cervix and uterine wall
CN1846233A (zh) * 2003-08-06 2006-10-11 特里格医疗有限公司 监测分娩参数的方法和装置
US20070015994A1 (en) * 2005-07-14 2007-01-18 Hyundae Hong In-vivo measurement of biomechanical properties of internal tissues
US20070167819A1 (en) * 2005-12-06 2007-07-19 Osborn Thomas W Iii Method for in-vivo measurement of biomechanical properties of internal tissues
CN101454792A (zh) * 2006-05-24 2009-06-10 梁伦锡 分娩控制系统及控制方法
US20080167581A1 (en) * 2007-01-10 2008-07-10 Yoav Paltieli Determining parameters associated with a female pelvis and cervix
CN101049245A (zh) * 2007-04-30 2007-10-10 陆尧胜 一种产程监护装置及其产程监护方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
M. HOUSE ET AL.: "The cervix as a biomechanical structure", 《ULTRASOUND IN OBSTETRICS & GYNECOLOGY》, vol. 28, no. 6, 30 November 2006 (2006-11-30), pages 745 - 749 *

Also Published As

Publication number Publication date
WO2011130143A1 (en) 2011-10-20
US20130218289A1 (en) 2013-08-22
CN102933156B (zh) 2015-07-15

Similar Documents

Publication Publication Date Title
CN102933156B (zh) 确定子宫准备分娩的情况的方法和系统
Hashemi et al. Global time-delay estimation in ultrasound elastography
Konofagou et al. Myocardial elastography—A feasibility study in vivo
Shiina et al. WFUMB guidelines and recommendations for clinical use of ultrasound elastography: Part 1: basic principles and terminology
CN101901335B (zh) 用于自动识别3d数据集中图像视图的方法和设备
Nenadic et al. Ultrasound elastography for biomedical applications and medicine
McCann et al. Multidimensional ultrasonic imaging for cardiology
US8118744B2 (en) Methods, systems and computer program products for ultrasound shear wave velocity estimation and shear modulus reconstruction
Barbosa et al. Fast and fully automatic 3-d echocardiographic segmentation using b-spline explicit active surfaces: Feasibility study and validation in a clinical setting
EP3231369A1 (en) Ultrasound diagnostic device and elasticity evaluation method
Kuo et al. Elastic properties of tendon measured by two different approaches
ES2964619T3 (es) Procedimiento para obtener un parámetro funcional de un músculo
EP3700429A1 (en) Method and apparatus for ultrasound measurement and imaging of biological tissue elasticity in real time
Zervantonakis et al. A novel, view-independent method for strain mapping in myocardial elastography: eliminating angle and centroid dependence
US9033883B2 (en) Flow quantification in ultrasound using conditional random fields with global consistency
Caenen et al. Investigating shear wave physics in a generic pediatric left ventricular model via in vitro experiments and finite element simulations
Lu et al. A real time displacement estimation algorithm for ultrasound elastography
Jiang et al. 6f-3 a regularized real-time motion tracking algorithm using dynamic programming for ultrasonic strain imaging
Dorosz et al. Three-dimensional visual guidance improves the accuracy of calculating right ventricular volume with two-dimensional echocardiography
Leung et al. A real-time intrasubject elastic registration algorithm for dynamic 2-D ultrasound images
Hjertaas et al. A left ventricular phantom for 3D echocardiographic twist measurements
Hu et al. A dynamic ultrasound simulation of a pulsating three‐layered CCA for validation of two‐dimensional wall motion and blood velocity estimation algorithms
Palmeri RSNA/QIBA efforts to standardize shear wave speed as a biomarker for liver fibrosis staging
Jia et al. 10B-5 4D Elasticity Imaging of PVA LV Phantom Integrated with Pulsatile Circulation System Using 2D Phased Array
Duann et al. Assessment of left ventricular cardiac shape by the use of volumetric curvature analysis from 3D echocardiography

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: 20150715

Termination date: 20160411

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