本发明涉及视频处理领域，特别涉及一种基于超声视频的对子宫微蠕动进行视频处理的方法。
背景技术：
子宫内膜随子宫内膜下肌层及子宫平滑肌层的非同步性收缩引起的不均匀性子宫腔内压力，而出现似肠道蠕动波一样的机械运动，称为子宫内膜波状运动，即子宫内膜蠕动波。子宫内膜蠕动和微蠕动能够反映子宫内膜容受性。子宫内膜蠕动波的方向是多种多样的，有宫颈到宫底的、有宫底到宫颈的、还有随机蠕动的等。子宫内膜蠕动波与精子运输、经血排出、胚胎植入及维持妊娠有关，蠕动波的方向和频率还会随着时间的变化而变化。除此之外，子宫内膜下肌层及子宫平滑肌层的非同步性收缩与血管的微小搏动都可能导致出现子宫微蠕动。
目前了解子宫蠕动波和微蠕动的方法是通过超声视频进行肉眼观察。由于蠕动的变化微弱，随机性强，传统超声设备采集的子宫视频中往往需要等待5～20分钟，才能通过仔细对比觉察到蠕动波，而对于子宫的微蠕动，往往更加不容易在超声视频中通过肉眼观察到的。
技术实现要素：
为了让子宫微蠕动变得可视化并能对子宫微蠕动的模式进行特征定量参考,本发明将超声视频中的子宫区域通过空间多尺度分解和相位放大进行微蠕动凸显，得到子宫微蠕动凸显视频；进而，基于预设的81种时空动态变化的标准频域特征字典，统计子宫微蠕动凸显视频中符合预设标准频域特征的频数分布，并将频数分布作为子宫的微蠕动模式；最终将子宫微蠕动凸显视频和微蠕动模式由显示模块输出展现。本发明的技术方案如下：
一种子宫超声微蠕动视频处理方法，如图1与图2所示，由以下模块组成：
1)超声采集模块，采集原始超声视频数据v1；
2)预处理模块，对原始超声视频数据进行视频去噪滤波和视觉标准化预处理，并将预处理后的超声视频数据v2发送给子宫区域检测模块；
3)子宫区域检测模块，对v2中的每一帧图像进行滑窗区域识别，将判断为子宫的区域按顺序整合为子宫视频数据v3，并将v3发送给子宫微蠕动凸显模块；
4)子宫微蠕动凸显模块，针对v3进行子宫微蠕动凸显得到v4；
5)微蠕动模式自动提取模块，负责提取v4中的微蠕动模式y；
6)显示模块，负责将v4和微蠕动模式y展现给用户。
优选地，为了能够采集超声视频数据，所述的超声采集模块的采集模式为b型超声或谐波成像超声，采集部位为子宫，采集得到包含子宫的原始超声视频数据v1。
优选地，为了满足对子宫微蠕动频率进行分析的奈奎斯特采样定律的要求，所述的超声采集模块的视频数据的采样频率不小于5hz。
优选地，为了去除噪声对视频质量的影响，所述的预处理模块将原始超声视频数据v1进行视频滤波去噪，视频滤波去噪的方法可以为高斯滤波，通过离散化窗口滑窗卷积来实现高斯滤波。
优选地，为了保证不同视频的图像一致性，所述的预处理模块将原始超声视频数据v1进行视觉标准化处理，视觉标准化处理可以为直方图均衡化。
优选地，为了降低非子宫区域对子宫区域的相位噪声，所述的子宫区域检测模块针对v2的每一帧进行子宫区域分割，分割方法可以为手工图像分割、或基于机器学习方法训练得到的识别模型，将每一帧的分割结果作为该帧的子宫区域，将仅有子宫区域的视频数据提取出来作为检测出的子宫视频数据v3。
优选地，为了实现对子宫微蠕动的凸显，即让肉眼不可见的子宫微蠕动变为肉眼可见，所述的子宫微蠕动凸显模块由尺度变换组、滤波器、相位放大器、逆尺度变换组构成；尺度变换组将子宫视频数据v3转换为不同尺度的视频数据，尺度变换方法可以为小波变换、高斯金字塔分解、可操纵方向金字塔变换；信号滤波器接收不同尺度的视频数据，并基于预设的带通范围对不同尺度的视频数据进行信号滤波后发送给相位放大器；相位放大器接收所述的信号滤波结果，并对其相位进行乘性扩增后发送给逆尺度变换组；逆尺度变换组接收所述的不同尺度的相位放大结果，并对其进行尺度逆向重建后得到子宫微蠕动凸显视频数据v4，最终将子宫微蠕动凸显视频数据v4发送给微蠕动模式自动提取模块。
优选地，为了提取子宫微蠕动的时空域特征模式，所述的微蠕动模式自动提取模块由频域特征编码器b、标准频域特征字典d、频数统计器构成；子宫微蠕动凸显视频数据v4的每一帧图像都通过傅立叶变换转换为频域数据fi(i＝1,2,…,k),k为总帧数；所有帧的频域数据f通过频域特征编码器b和特征简并方法j得到q种标准频域特征(0<q<＝81,且q为整数)；频数统计器按照标准频域特征字典的排列顺序统计上述q种标准频域特征的频数分布，将频数分布作为微蠕动模式y发送给显示模块。
优选地，为了对频域的时空特征进行描述，所述的频域特征编码器b将接收的频域数据f划分为一系列大小相同的子空间；子空间大小为a×b×c，a的范围是0.05～0.5m,b的范围是0.05～0.5n,c的范围是0.05～0.5k,a、b、c均为四舍五入后的正整数；视频v4的画面大小为m×n，总帧数为k，m、n、k均为正整数；每个子空间的最大实数和最大虚数为此子空间内所有点中的最大实数与最大虚数；每个子空间的最小实数和最小虚数为此子空间内所有点中的最小实数与最小虚数；每个子空间的中值实数和中值虚数为此子空间内所有点的中值实数与中值虚数；任意子空间与其邻域子空间构成3×3×3的三维矩阵，矩阵中每一个位置都包含3个复数，分别是最大实数与最大虚数、最小实数与最小虚数、中值实数与中值虚数；频域特征编码器将任意子空间与其邻域子空间的所有频域复数按顺序重整为一个1×162(162＝3×3×3×3×2)的一维向量，实数与对应的虚数紧挨排列；对这1×162的一维向量的每一个数值，大于0的被标记为1，小于等于0的被标记为0，从而得到1×162的只含1和0的一维向量，此一维向量即为二进制频域特征编码。
以三维3×3×3矩阵举例，给27个子空间分别按顺序编号为1、2、3、……、26、27，最中间的子空间就是第14个子空间，邻域子空间就是包围最中间的子空间的其余26个子空间。
优选地，为了对2162个二进制频域特征编码进行精简，所述的特征简并方法j将162位的二进制频域特征编码看作首尾相连的圆环，从逆时针方向遍历圆环上的每个二进制数，将出现从0到1或1到0的跳变总次数记为z；将具有相同z值的二进制频域特征编码简并为同一种标准频域特征，其中z的取值范围是0～80；所有81种频域特征组成标准频域特征字典。
以4位二进制频域特征编码为例来说明一下。比如1000、0110、1010、0101都是二进制频域特征编码；1000、0110都只有2次从0到1或从1到0的跳变，它们就被简并归为同一类标准频域特征，设定为具有相同的标准频域特征，编号为2；1010、0101都只有4次从0到1或从1到0的跳变，它们就归为另一个标准频域特征，编号为4；依次类推。之所以这样设定，是因为相似的频域特征具有相同的从0到1或从1到0的跳变。为了尽量去除频域特征的冗余性，同时也是为了减少计算量与存储空间，需要对频域特征进行简并。
本发明具有如下有益效果：
1)本发明将肉眼不可见的子宫超声微蠕动变为肉眼可见；
2)本发明能够提取子宫微蠕动的模式，为子宫微蠕动的模式分析提供了特征定量参考。
附图说明
图1是本发明提取子宫微蠕动模式的概括处理流程图；
图2是本发明提取子宫微蠕动模式的详细处理流程图；
图3是本发明实施例中子宫微蠕动凸显前后的时空切片对比图；
图4是本发明实施例中提取典型女性的子宫微蠕动模式图，横轴为频域特征序号，纵轴为相应频域特征的频数。
具体实施方式
以下通过实施例对本发明做进一步说明，以便更好地理解本发明的技术方案，但本发明并不局限于此。下述实施例中所用的超声设备、耦合剂等，如无特殊说明，均可从商业途径得到。
1.超声视频采集
筛选1名典型女性，在签署知情同意书后，纳入本实施例。采用普通超声诊断仪与阴道探头采集典型女性的子宫超声视频。采集模式为谐波成像超声，采集部位为子宫，采集帧率为30hz，探头中心频率为10mhz。
2.预处理
对原始超声视频数据进行视频去噪滤波和视觉标准化预处理。预处理模块对原始超声视频进行视频滤波去噪，视频滤波去噪的方法选择高斯滤波，通过离散化窗口滑窗卷积来实现高斯滤波，滑窗大小选择为5×5的矩阵。预处理模块对原始超声视频进行视频滤波去噪后再进行视觉标准化处理，视觉标准化处理方法选择直方图均衡化。
3.子宫区域检测
在预处理后的超声视频检测子宫区域并得到仅包含子宫区域的视频。针对预处理后的超声视频的每一帧进行子宫区域分割，分割方法基于机器学习方法训练得到的检测模型，将每一帧的分割结果作为该帧的子宫区域，将仅有子宫区域的视频数据提取出来作为检测出的子宫视频数据。此基于机器学习方法训练得到的子宫检测模型由滑窗核组、子宫分类器、聚类筛选器组成；滑窗核组由8个不同尺寸的滑窗核组成，从上述视频中的每帧图像中各自裁取不同尺寸的一小块图像作为感兴趣区域图像，直至遍历完所有图像数据；子宫分类器为一个人工神经网络，由15-30个卷积层构成，每个卷积层由100-250000个神经元组成，底层神经元读取滑窗核组传来的多组感兴趣区域图像，并将顶层神经元输出的多组识别数据发送给聚类筛选器；聚类筛选器为mean-shift聚类，mean-shift搜索半径为6个像素，对所述的多组分类数据进行聚类，将每一帧的聚类中心作为该帧的子宫中心，最终将聚类中心周围的子宫视频提取出来作为检测出的子宫视频。
4.子宫微蠕动凸显
子宫微蠕动凸显模块由尺度变换组、滤波器、相位放大器、逆尺度变换组构成；尺度变换组采用小波变换将子宫视频转换为不同尺度的视频；信号滤波器接收不同尺度的视频数据，并基于预设的带通范围0.1～14.9hz对不同尺度的视频数据进行信号滤波后发送给相位放大器；相位放大器接收所述的信号滤波结果，并对其相位进行10倍乘性扩增后发送给逆尺度变换组；逆尺度变换组接收所述的不同尺度的相位放大结果，并对其进行尺度逆向重建后得到子宫微蠕动凸显视频，并发送给显示模块。
图3是本发明实施例中子宫微蠕动凸显前后的时空切片对比图。图3.a是一名典型女性的原始子宫超声视频，虚线为时空切片采集的纵向单列位置；图3.b是相应的子宫微蠕动凸显视频，虚线为与图3.a中相同的时空切片采集的纵向单列位置；在视频中取一条纵切的像素在时间轴上展开并叠加为一幅图，就能得到时空切片图。时空切片图的横轴为视频的第1、2、……、第k帧，纵轴为对应的1、2、……、第k帧的相同位置的单列像素，其中k为视频总帧数。图3.c为原始子宫超声视频的子宫超声视频的时空切片图；图3.d为相应子宫微蠕动凸显视频的时空切片图。
5.微蠕动模式自动提取
微蠕动模式自动提取模块由频域特征编码器、标准频域特征字典、频数统计器构成；子宫微蠕动凸显视频数据v4的每一帧图像都通过傅立叶变换转换为频域数据fi(i＝1,2,…,300),300为视频总帧数；所有帧的频域数据f通过频域特征编码器b和特征简并方法j得到40种标准频域特征；频数统计器按照标准频域特征字典的排列顺序统计上述40种标准频域特征的频数分布，将频数分布作为微蠕动模式y发送给显示模块。
频域特征编码器接收所述的频域数据f，并划分为一系列的相同大小的子空间；视频画面大小为200×200，总帧数为300；子空间大小为20×20×30。每个子空间的最大实数和最大虚数为此子空间内所有点中的最大实数与最大虚数；每个子空间的最小实数和最小虚数为此子空间内所有点中的最小实数与最小虚数；每个子空间的中值实数和中值虚数为此子空间内所有点的中值实数与中值虚数；任意子空间与其邻域子空间构成3×3×3的三维矩阵，矩阵中每一个位置都包含3个复数，分别是最大实数与最大虚数、最小实数与最小虚数、中值实数与中值虚数。
频域特征编码器将任意子空间与其邻域子空间的所有频域复数按顺序重整为一个1×162(162＝3×3×3×3×2)的一维向量，实数与对应的虚数紧挨排列；对这1×162的一维向量的每一个数值，大于0的被标记为1，小于等于0的被标记为0，从而得到1×162的只含1和0的一维向量，此一维向量即为二进制频域特征编码。
本实施例对1名典型女性进行子宫超声视频采集，并根据本发明的处理方法对超声视频数据进行处理后，得到如图3和图4所示的结果。图3是本发明实施例中子宫微蠕动凸显前后的时空切片对比图，表明本发明能够对子宫微蠕动进行有效凸显。图4是本发明实施例中提取典型女性的子宫微蠕动模式图，表明本发明能够有效提取子宫微蠕动模式，为进行微蠕动模式的比较提供定量特征参考。