Traitement en cours

Veuillez attendre...

Paramétrages

Paramétrages

Aller à Demande

1. WO2021042908 - PROCÉDÉ ET APPAREIL DE DÉTECTION DE VIBRATIONS DE VOIE, ET DISPOSITIF DE DÉTECTION DE VIBRATIONS

Document

说明书

发明名称 0001   0002   0003   0004   0005   0006   0007   0008   0009   0010   0011   0012   0013   0014   0015   0016   0017   0018   0019   0020   0021   0022   0023   0024   0025   0026   0027   0028   0029   0030   0031   0032   0033   0034   0035   0036   0037   0038   0039   0040   0041   0042   0043   0044   0045   0046   0047   0048   0049   0050   0051   0052   0053   0054   0055   0056   0057   0058   0059   0060   0061   0062   0063   0064   0065   0066   0067   0068   0069   0070   0071   0072   0073   0074   0075   0076   0077   0078   0079   0080   0081   0082   0083   0084   0085   0086   0087   0088   0089   0090   0091   0092   0093   0094   0095   0096   0097   0098   0099   0100   0101   0102   0103   0104   0105   0106   0107   0108   0109   0110   0111   0112   0113   0114   0115   0116   0117   0118   0119   0120   0121   0122   0123   0124   0125   0126   0127   0128   0129   0130   0131   0132   0133   0134   0135   0136   0137   0138   0139   0140   0141   0142   0143   0144   0145   0146   0147   0148   0149   0150  

权利要求书

1   2   3   4   5   6   7   8   9   10  

附图

1   2A   2B   2C   3   4   5   6  

说明书

发明名称 : 轨道振动检测方法及装置、振动检测设备

[0001]
本申请要求于2019年09月04日提交中国专利局、申请号为2019108358154、发明名称为“轨道振动检测方法及装置、振动检测设备”的中国专利申请的优先权,其全部内容通过引用结合在本申请中。

技术领域

[0002]
本申请涉及计算机技术领域,具体涉及一种轨道振动检测方法及装置、振动检测设备。

背景技术

[0003]
轨道是指用条形的钢轨铺成的供火车、电车等行驶的路线。其中,维护人员需要定期对轨道进行维护与检查,以保证交通运输的安全。在轨道维护与检查工作中,轨道振动检测是一个重要的环节。
[0004]
在传统振动检测中,维护人员往往通过个人经验或者简单振动检测设备排查轨道接头、轨道扣件和道床等是否存在安全隐患,这也导致轨道振动检测的检测效率低以及检测结果不精准等问题。此外,在相关技术中,通过相关设备拍摄轨道振动时的视频,并对轨道振动视频进行处理和分析,以减小人工成本和检测时间,但是轨道振动视频也存在振动部位的振动幅度细微,肉眼难以观测,导致振动信息的提取难度大。
[0005]
发明内容
[0006]
本申请提供了一种轨道振动检测方法及装置、振动检测设备,以期望提升振动信息提取的准确性和应用的广泛性,进而提高振动检测的效率、准确和可靠性。
[0007]
第一方面,本申请实施例提供一种轨道振动检测方法,应用于振动检测设备,所述方法包括:
[0008]
获取待测轨道的轨道振动视频,并根据所述待测轨道的轨道振动视频选取第一RGB图像;
[0009]
将所述第一RGB图像由RGB颜色空间线性转换到YIQ颜色空间,获得第一YIQ图像;
[0010]
通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列;
[0011]
根据所述特征图像序列提取所述待测轨道的实际振动数据;
[0012]
根据所述实际振动数据确定所述待测轨道的稳定性。
[0013]
第二方面,本申请实施例提供一种轨道振动检测装置,应用于振动检测设备,所述装置包括处理单元和通信单元,其中:
[0014]
所述处理单元,用于获取待测轨道的轨道振动视频,并根据所述待测轨道的轨道振动视频选取第一RGB图像;以及用于将所述第一RGB图像由RGB颜色空间线性转换到YIQ颜色空间,获得第一YIQ图像;以及用于通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列;以及用于根据所述特征图像序列提取所述待测轨道的实际振动数据;以及用于根据所述实际振动数据确定所述待测轨道的稳定性。
[0015]
第三方面,本申请实施例提供一种振动检测设备,包括处理器、存储器、通信接口,以及一个或多个程序,所述一个或多个程序被存储在所述存储器中,并且被配置由所述处理器执行,所述程序包括用于执行本申请实施例第一方面任一方法中的步骤的指令。
[0016]
第四方面,本申请实施例提供一种计算机可读存储介质,其中,上述计算机可读存储介质存储用于电子数据交换的计算机程序,并且上述计算机程序使得计算机执行如本申请实施例第一方面任一方法中所描述的部分或全部步骤。
[0017]
第五方面,本申请实施例提供一种计算机程序产品,其中,上述计算机程序产品包括存储了计算机程序的非瞬时性计算机可读存储介质,并且上述计算机程序可操作来使计算机执行如本申请实施例第一方面任一方法中所描述的部分或全部步骤。该计算机程序产品可以为一个软件安装包。
[0018]
可以看出,本申请实施例公开了一种轨道振动检测的方法及装置,通过获取待测轨道的轨道振动视频,并根据所述待测轨道的轨道振动视频选取第一RGB图像;将所述第一RGB图像由RGB颜色空间线性转换到YIQ颜色空间,获得第一YIQ图像;通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列;根据所述特征图像序列提取所述待测轨道的实际振动数据;根据所述实际振动数据确定所述待测轨道的稳定性,提升振动数据提取的准确性以及工程应用的广泛性,进而提升轨道振动检测的效率、准确和可靠性。

附图说明

[0019]
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根 据这些附图获得其他的附图。
[0020]
图1是本申请实施例提供的一种轨道振动检测的系统示意图;
[0021]
图2a为本申请实施例提供的一种轨道振动检测方法的流程示意图;
[0022]
图2b是本申请实施例提供的一种可能的处理轨道振动视频的示意图;
[0023]
图2c是本申请实施例提供的一种可能的轨道扣件的平面结构示意图;
[0024]
图3是本申请实施例提供的另一种轨道振动检测方法的流程示意图;
[0025]
图4是本申请实施例提供的另一种轨道振动检测方法的流程示意图
[0026]
图5是本申请实施例提供的一种振动检测设备的结构示意图;
[0027]
图6是本申请实施例提供的一种轨道振动检测装置的功能单元组成框图。

具体实施方式

[0028]
为了使本技术领域的人员更好地理解本申请方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
[0029]
本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别不同对象,而不是用于描述特定顺序。此外,术语“包括”和“具有”以及它们任何变形,意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元,而是可选地还包括没有列出的步骤或单元,或可选地还包括对于这些过程、方法、产品或设备固有的其他步骤或单元。
[0030]
在本文中提及“实施例”意味着,结合实施例描述的特定特征、结构或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例,也不是与其它实施例互斥的独立的或备选的实施例。本领域技术人员显式地和隐式地理解的是,本文所描述的实施例可以与其它实施例相结合。
[0031]
下面对本申请实施例进行详细介绍。
[0032]
请参阅图1,图1是本申请实施例提供的一种轨道振动检测的系统示意图。图1包括轨道100和振动检测设备110。
[0033]
具体的,本申请的轨道100主要包括钢轨101、轨枕102、轨道接头103、轨道扣件104、道床105等等。其中,轨枕102包括木枕1021和混凝土轨枕1022,其承受来自钢轨101 传来的垂直力和水平力,并将这些力分布于床道105上,同时有效地保持轨道100的轨距、方向和位置。因此,轨枕102需要有一定的坚固性、耐久性和弹性,同时具有足够的阻力,以避免在列车作用下发生横向移动。
[0034]
具体的,本申请的钢轨连接结构分为中间连接结构和接头连接结构两类。其中,中间连接为钢轨101与轨枕102之间的连接,中间连接结构通称轨道扣件104,其主要作用是阻止钢轨101作相对与轨枕102的纵横向移动,并保持钢轨101的稳固位置。轨道扣件104包括木枕扣件1041和混泥土轨枕扣件1042。木枕扣件1041包括挡板、道钉垫板和普通道钉,混泥土轨枕扣件1042包括挡板、道钉垫板和普通道钉。其次,接头连接为钢轨101与钢轨101的连接。接头连接由夹板、螺栓、弹簧垫圈等组成,通称轨道接头103,其作用是在接头处把钢轨101连接起来,使钢轨接头部分具有与钢轨101一样的整体性,以抵抗弯曲和位移。
[0035]
具体的,本申请的道床105的作用是把轨枕102传来的力均匀地传布到路基上,固定轨枕102位置,保持轨道100的稳定性,排除路基水分,保持轨道100的弹性,并调整轨道100的平面及纵断面。
[0036]
具体的,本申请的振动检测设备110包括具备视频处理能力的电子设备111、摄像头113和传感器112。其中,电子设备111可以包括各种具有视频处理功能的手持设备、车载设备、可穿戴设备、计算设备或连接到无线调制解调器的其他视频处理设备,以及各种形式的智能终端设备(intelligent terminal device)等等。摄像头112可以包括红外摄像头和可见光摄像头,可见光摄像头也可以包括普通摄像头或广角摄像头,在此不作限定。摄像头113主要用于录制轨道100的一段视频,并上传给电子设备111,电子设备111中的轨道振动检测装置再对上传的一段进行振动检测。传感器112包括超声波传感器,用于检测轨道100。
[0037]
下面具体描述一种轨道振动检测的执行步骤,请参阅图2a。图2a为本申请实施例提供的一种轨道振动检测方法的流程示意图,应用于振动检测设备,所述方法包括:
[0038]
S201、获取待测轨道的轨道振动视频,并根据所述待测轨道的轨道振动视频选取第一RGB图像;
[0039]
具体的,所述待测轨道的轨道振动视频可以通过以下操作获得:根据所述振动检测设备拍摄待测轨道,并得到所述待测轨道的第一视频和第二视频,其中第一视频和第二视频为所述待测轨道在相同时间内拍摄的不同视频;获取述第一视频对应的第一图像和第二视 频对应的第二图像;对第一图像和第二图像进行重叠,再清除第一图像和第二图像不能重叠的像素点,并获得所述待测轨道的轨道振动视频。
[0040]
其中,在根据所述振动检测设备拍摄待测轨道时,存在一些外界干扰因素,例如摄像头晃动、摄像头故障等,导致所述待测轨道的轨道振动视频与真实情况存在一定的偏差。因此,通过采用不同摄像头在同一时间内拍摄所述待测轨道的不同视频,获得所述待测轨道的第一视频和第二视频,并且对第一视频对应的第一图像和第二视频对应的第二图像进行重叠。在摄像头不存在上述干扰因素的情况下,所述第一图像和所述第二图像对应的像素点可以完全重叠,通过清除所述第一图像和所述第二图像不能重叠的像素点,即不同摄像头中的拍摄偏差像素点,可以获得噪声更少的图像。同样,可以对所述待测轨道拍摄更多的同一时间段内的不同视频,以进一步减少图像噪声。
[0041]
S202、将所述第一RGB图像由RGB颜色空间线性转换到YIQ颜色空间,获得第一YIQ图像;
[0042]
具体的,RGB是工业界的一种由红基色(Red)、绿基色(Green)和蓝基色(Blue)三个颜色分量表示的颜色标准。RGB颜色空间利用三个颜色分量的线性组合来表示颜色,而YIQ是美国全国电视标准委员会(National Television Standards Committee,NTSC)的电视系统标准。Y表示提供黑白电视和彩色电视的亮度分量,I(In-phase)表示色彩从橙色到青色的色度分量,Q(Quadrature-phase)表示色彩从紫色到黄绿色的色度分量。YIQ颜色空间能够将视频图像中的亮度分量分离提取出来,并且YIQ颜色空间与RGB颜色空间之间是线性变换的关系,具有计算量小和聚类特性好等优点,可以适应光照强度不断变化的场合。
[0043]
具体的,RGB颜色空间和YIQ颜色空间的线性转换关系为:
[0044]
Y=a1*R+a2*G+a3*B;
[0045]
I=b1*R+b2*G+b3*B;
[0046]
Q=c1*R+c2*G+c3*B;
[0047]
其中,a1、a2、a3、b1、b2、b3、c1、c2和c3取值范围为[-1,1]。
[0048]
优选的,RGB颜色空间和YIQ颜色空间的线性转换关系为:
[0049]
Y=0.299*R+0.587*G+0.114*B;
[0050]
I=0.596*R-0.275*G-0.321*B;
[0051]
Q=0.212*R-0.523*G+0.311*B。
[0052]
举例中,请参阅图2b,图2b是本申请实施例提供的一种可能的处理轨道振动视频的示意图。其中,振动检测设备110拍摄得到一段轨道100的轨道振动视频,并采集到轨道振动视频中轨道扣件104振动的图像。通过RGB颜色空间分别提取所述第一图像的红基色分量、绿基色分量和蓝基色分量。根据RGB颜色空间和YIQ颜色空间的线性转换关系,得到第一YIQ图像的Y图像分量、Q图像分量和I图像分量。
[0053]
S203、通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列;
[0054]
具体的,所述视频放大算法可以包括以下至少一种:拉普拉斯运动放大算法、欧拉运动放大算法、复数相位运动放大算法、里斯金字塔运动放大算法。
[0055]
S204、根据所述特征图像序列提取所述待测轨道的实际振动数据;
[0056]
S205、根据所述实际振动数据确定所述待测轨道的稳定性。
[0057]
具体的,通过将实际振动数据与待测轨道中的轨道接头103、轨道扣件104和道床105的振动数据进行比较,确定所述待测轨道的稳定性。
[0058]
可以看出,本申请实施例中所描绘的轨道振动检测方法,应用于轨道振动设备。所述轨道振动设备通过视频放大算法从轨道振动视频中提取振动数据,以提升了振动数据提取的准确性,进而提高轨道振动检测的效率、准确和可靠性。
[0059]
在一个可能的示例中,所述通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列,包括:对所述第一YIQ图像中的Y通道图像进行傅里叶变换;以及对傅里叶变换后的Y通道图像进行图像金字塔分解;以及对图像金字塔分解后的Y通道图像进行归一化、时空滤波以及线性放大处理;将处理后的Y通道图像与所述YIQ图像中的I通道图像和Q通道图像进行合成,并形成特征图像序列。
[0060]
其中,Y通道图像表示亮度信息,通过对Y通道图像进行傅里叶变换,以将时域的亮度信息转化为频域的相位变化。
[0061]
具体的,图像金字塔是一种空域对图像进行多分辨率处理的方法。物体图像的尺寸和对比度的不同,通过对物体图像进行多分辨率的分析,能有利于分析物体图像和提取物体图像的特征参数。此外,通过时频小波分解,将物体图像分解为高频和低频分量,再经过二抽样取得小波分解系数。根据小波分解系数将物体图像分解为多个尺度和方向,以获得多种物体图像的子带图像。
[0062]
具体的,图像归一化是一种对图像进行一系列标准变换处理,并使其变换为标准图像 的过程,所述标准图像称作归一化图像。对图像进行处理,并得到多种子带图像,所述多种子带图像在经过相同参数的图像归一化处理后能够得到相同形式的标准图像。进一步的,图像归一化就是通过一系列变换将图像转换成相应的唯一标准形式,所述该标准形式图像对平移、旋转、缩放等仿射变换具有不变特性。
[0063]
可选的,所述通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列,包括:根据所述振动检测设备拍摄待测轨道,并得到所述待测轨道的轨道振动视频;获取所述轨道振动视频的多帧的待处理图像,并将所述待处理图像进行分区;将每个分区中的像素点作为初始特征点,并基于最小差值平方和SSD匹配,以计算所述初始特征点的流向量;根据所述初始特征点对应的流向量计算所述初始特征点的偏移距离;对多个初始特征点对应的多个偏移距离采用k-means聚类算法进行聚类,获得多个聚类类簇;确定所述多个聚类类簇中每个聚类类簇中的偏移距离平均值是否处于预设范围;若是,则确定所述聚类类簇中的初始特征点对应的所述分区为运动分区;保留所述目标视频的多帧待处理图像中的运动分区,形成特征图像序列。
[0064]
可选的,所述通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列,包括:对所述第一YIQ图像中的Y通道图像进行傅里叶变换;对傅里叶变换后的Y通道图像进行下采样处理,得到第一子图像,所述下采样处理为降低图像的分辨率;对所述第一子图像进行上采样处理,得到第二子图像,所述上采样处理为提升图像的分辨率;将傅里叶变换后的Y通道图像与所述第二子图进行像素处理,得到第三子图;将所述第三子图进行时域滤波处理,得到目标频带;根据所述目标频段与所述第一YIQ图像中的Y通道图像确定多个Y通道图像信号;对所述多个Y通道图像信号进行放大,得到放大后的多个Y通道图像信号;根据拉普拉斯金字塔重建算法对放大后的多个Y通道图像信号进行合成;获取所述第一YIQ图像中的I通道图像与Q通道图像,并将合成后的多个Y通道图像与所述I通道图像和所述Q通道图像进行相加,形成特征图像序列。
[0065]
可选的,所述通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列,包括:对所述第一YIQ图像中的Y通道图像进行傅里叶变换;对傅里叶变换后的Y通道图像进行空间滤波,得到不同空间分辨率的Y通道图像;对不同空间分辨率的Y通道图像进行时域滤波处理,获取目标频段;根据所述目标频段与所述不同空间分辨率的Y通道图像确定多个Y通道图像信号;对所述多个Y通道图像信号进行放大,得到放大后的多个Y通道图像信号;根据复数可操纵金字塔重建算法对放大后的多个Y通道图像 信号进行合成;获取所述第一YIQ图像中的I通道图像与Q通道图像,并将合成后的多个Y通道图像与所述I通道图像和所述Q通道图像进行相加,形成特征图像序列。
[0066]
可选的,所述通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列,包括:对所述第一YIQ图像中的Y通道图像进行傅里叶变换;根据拉普拉斯金字塔对傅里叶变换后的Y通道图像进行分解;将分解后的Y通道图像进行里斯变换;对里斯变换后的Y通道图像进行正交和相位处理;对处理后的Y通道图像进行空间时域滤波;对空间时域滤波后的Y通道图像进行方法和相移,形成特征图像序列。
[0067]
可以看出,所述轨道振动设备通过图像金字塔对第一YIQ图像中的Y通道图像进行分解,并对图像金字塔分解后的Y通道图像进行归一化、时空滤波以及线性放大处理,以放大轨道振动视频,提升振动数据提取的准确性,进而提高轨道振动检测的效率、准确和可靠性。
[0068]
在一个可能的示例中,所述根据所述特征图像序列提取所述待测轨道的实际振动数据,包括:根据相位相关公式计算所述特征图像序列的交叉互功率谱;以及根据所述交叉互动功率谱提取所述特征图像序列中像素的振动信息;以及根据所述振动信息确定所述待测轨道的实际振动数据。
[0069]
具体的,所述相位相关算法采用如下的公式计算所述特征图像序列的交叉互功率谱:
[0070]
[0071]
其中,F a为a帧图像的傅立叶变换, 为b帧图像的傅里叶变换的共轭信号, 表示F a的相关积的模,R表示交叉互功率谱(包含频域噪音)。
[0072]
可以看出,振动检测设备通过对视频放大后的特征图像序列采用相位相关算法计算交叉互功率谱,以提高轨道振动检测的准确性。
[0073]
在一个可能示例中,所述根据所述交叉互动功率谱提取所述特征图像序列中像素的振动信息,包括:根据所述交叉互功率谱的峰值位置选择自适应滤波器;以及根据所述自适应滤波器对所述交叉互功率谱进行滤波;以及对滤波后的交叉互功率谱进行反傅里叶变换;以及根据反傅里叶变换后的交叉互功率谱提取所述特征图像序列中像素的振动信息。
[0074]
具体的,由于待测轨道的轨道振动是有方向性的,对处理后的Y通道图像进行归一化、时空滤波以及线性放大处理,得到的特征图像序列包含有冗余信号。对冗余信号都进行差值、加权和重构,将带来很高的运算量,而且其对振动信息的分析没有作用。通过采用主 成分分解对特征图像序列进行降维处理,再用自适应滤波器对降维后的特征图像序列进行滤波,能有效减小运算处理时间和噪声信号的影响。
[0075]
可以看出,振动检测设备通过自适应滤波器对特征图像序列进行滤波,以减小轨道振动检测的运算量,以及提高轨道振动检测的效率,有效保证轨道的稳定性。
[0076]
在一个可能的示例中,所述根据所述实际振动数据确定所述待测轨道的稳定性,包括:根据所述实际振动数据确定所述待测轨道的振动部位,所述振动部位包括轨道扣件或轨道接头或道床;以及根据所述待测轨道的振动部位和所述实际振动数据确定所述待测轨道的稳定性。
[0077]
可以看出,振动检测设备通过待测轨道的振动部位和所述实际振动数据确定所述待测轨道的稳定性,以保证轨道振动检测的准确和可靠性。
[0078]
在一个可能的示例中,所述振动部位包括轨道扣件或轨道接头时,所述根据所述待测轨道的振动部位和所述实际振动数据确定所述待测轨道的稳定性,包括:采集所述待测轨道的轨道扣件或轨道接头的N段振动视频;以及统计所述N段振动视频的平均振动频率、平均振幅和平均振动周期;以及比较所述平均振动频率、所述平均振幅和所述平均振动周期与所述实际振动数据,并根据比较结果确定所述待测轨道的稳定性。
[0079]
具体的,图2c是本申请实施例提供的一种可能的轨道扣件的平面结构示意图。其中,混泥土轨枕扣件1042包括挡板301、普通道钉302和道钉垫板303。振动检测设备采集混泥土轨枕扣件1042的N段振动视频,并通过N端振动视频统计挡板301和普通道钉302在上下振动和左右振动方向的平均振动频率、平均振幅和平均振动周期。
[0080]
可以看出,振动检测设备通过统计待测轨道的轨道扣件的平均振动数据,并将所述平均振动数据比较实际振动数据,以提高轨道振动检测的准确和可靠性。
[0081]
在一个可能的示例中,所述振动部位包括道床时,所述根据所述待测轨道的振动部位和所述实际振动数据确定所述待测轨道的稳定性,包括:获取道床的道床系数;以及根据所述道床系数确定所述道床的振动刚度和振动弹性系数;以及比较所述振动刚度和所述振动弹性系数与所述实际振动数据,并根据比较结果确定所述待测轨道的稳定性。
[0082]
可以看出,振动检测设备通过获取待测轨道的道床的振动刚度和振动弹性系数,并将所述振动刚度和所述振动弹性系数比较实际振动数据,以提高轨道振动检测的准确和可靠性。
[0083]
与上述图2a所述的实施例一致,请参阅图3。图3是本申请实施例提供的另一种轨道 振动检测方法的流程示意图,应用于振动检测设备,所述方法包括:
[0084]
S301、获取待测轨道的轨道振动视频,并根据所述待测轨道的轨道振动视频选取第一RGB图像;
[0085]
具体的,所述待测轨道的轨道振动视频可以通过以下操作获得:根据所述振动检测设备拍摄待测轨道,并得到所述待测轨道的第一视频和第二视频,其中第一视频和第二视频为所述待测轨道在相同时间内拍摄的不同视频;获取述第一视频对应的第一图像和第二视频对应的第二图像;对第一图像和第二图像进行重叠,再清除第一图像和第二图像不能重叠的像素点,并获得所述待测轨道的轨道振动视频。
[0086]
S302、将所述第一RGB图像由RGB颜色空间线性转换到YIQ颜色空间,获得第一YIQ图像;
[0087]
具体的,RGB颜色空间和YIQ颜色空间的线性转换关系为:
[0088]
Y=0.299*R+0.587*G+0.114*B;
[0089]
I=0.596*R-0.275*G-0.321*B;
[0090]
Q=0.212*R-0.523*G+0.311*B。
[0091]
S303、对所述第一YIQ图像中的Y通道图像进行傅里叶变换,并对傅里叶变换后的Y通道图像进行图像金字塔分解;
[0092]
其中,Y通道图像表示亮度信息,通过对Y通道图像进行傅里叶变换,以将时域的亮度信息转化为频域的相位变化。
[0093]
其中,图像金字塔包括高斯金字塔、拉普拉斯金字塔、复数可操纵金字塔、里斯金字塔等等。
[0094]
S304、对图像金字塔分解后的Y通道图像进行归一化、时空滤波以及线性放大处理;
[0095]
S305、将处理后的Y通道图像与所述YIQ图像中的I通道图像和Q通道图像进行合成,并形成特征图像序列;
[0096]
S306、根据相位相关公式计算所述特征图像序列的交叉互功率谱,并根据所述交叉互功率谱的峰值位置选择自适应滤波器;
[0097]
S307、根据所述自适应滤波器对所述交叉互功率谱进行滤波,并对滤波后的交叉互功率谱进行反傅里叶变换;
[0098]
S308、根据反傅里叶变换后的交叉互功率谱提取所述特征图像序列中像素的振动信息,并根据所述振动信息确定所述待测轨道的实际振动数据;
[0099]
S309、根据所述实际振动数据确定所述待测轨道的振动部位,所述振动部位包括轨道扣件或轨道接头或道床;
[0100]
S310、采集所述待测轨道的轨道扣件或轨道接头的N段振动视频,并统计所述N段振动视频的平均振动频率、平均振幅和平均振动周期;
[0101]
S311、比较所述平均振动频率、所述平均振幅和所述平均振动周期与所述实际振动数据,并根据比较结果确定所述轨道扣件或轨道接头的稳定性。
[0102]
可以看出,本申请实施例中所描绘的轨道振动检测方法,应用于轨道振动设备。所述轨道振动设备通过视频放大算法对轨道振动视频进行放大,并根据相位相关公式和反傅里叶变换提取轨道振动视频的实际振动数据,再比较轨道扣件的平均振动频率、平均振幅和平均振动周期与所述实际振动数据,以提升了振动数据提取的准确性,进而提高轨道振动检测的效率、准确和可靠性。
[0103]
请参阅图4,图4是本申请实施例提供的另一种轨道振动检测方法的流程示意图。
[0104]
S401、获取待测轨道的轨道振动视频,并根据所述待测轨道的轨道振动视频选取第一RGB图像;
[0105]
S402、处理所述第一RGB图像,并得到多个分辨率对应的多个子带图像序列;
[0106]
S403、按照预设的分区灰度值筛选策略从所述多个子带图像序列筛选出用于放大处理的至少一个子带图像序列;
[0107]
具体的,所述按照预设的分区灰度值筛选策略从所述多个子带图像序列筛选出用于放大处理的至少一个子带图像序列包括:确定所述多个子带图像序列的前景图像和后景图像,所述前景图像包含所述被检测产品的发生往复运动的区域影像,所述后景图像为除所述被检测产品的影像之外的影像;以及确定所述前景图像在子带图像中的面积占比;以及根据所述面积占比和预设的子分区计算公式确定所述前景图像的子分区数量,并按照所述子分区数量划分所述前景图像为多个前景子分区;以及针对每个子带图像序列,执行以下操作(1)-(6)以得到所述每个子带图像序列的灰度值变化频率:(1)确定当前处理的子带图像序列中每个子带图像的每个前景子分区的被测像素点;(2)根据每个被测像素点在当前子带图像序列包含的多个子带图像的灰度值生成每个被测像素点的灰度值时域变化波形图;(3)针对所述每个前景子分区执行如下(a)(b)(c)操作:(a)根据当前处理的前景子分区包含的多个被测像素点的灰度值时域变化波形图,确定所述当前处理的前景子分区是否包含灰度值周期性变化的被测像素点;(b)若是,则标记所述当前处理的前景子分区为 被选择的前景子分区;(c)若否,则标记所述当前处理的前景子分区为未被选择的前景子分区;(4)针对标记后的被选择的多个前景子分区,按照区域关联性将具有相邻关系的前景子分区拼接为振动参考区域;(5)确定所述振动参考区域中多个像素点中灰度值发生周期性变化的多个参考像素点,以及确定每个参考像素点的灰度值变化频率;(6)加权计算所述振动参考区域中所述多个参考像素点的灰度值变化频率,得到所述当前处理的子带图像序列的灰度值变化频率;以及根据所述每个子带图像序列的灰度值变化频率筛选出符合预设参考振动频率的至少一个子带图像序列。
[0108]
其中,预设的子分区计算公式为:
[0109]
[0110]
其中,x为面积占比,y为子分区数量,x大于0且小于或等于1。
[0111]
其中,所述确定当前处理的子带图像序列中每个子带图像的每个前景子分区的被测像素点包括以下至少一种:所述每个前景子分区的边缘像素点;所述每个前景子分区的中间区域的像素点;所述每个前景子分区的随机筛选出来的多个像素点。
[0112]
S404、放大处理所述至少一个子带图像序列,并得到放大处理后的至少一个子带图像序列;
[0113]
S405、将所述放大处理后的至少一个子带图像序列和所述多个子带图像序列中除所述至少一个子带图像序列之外的子带图像序列进行融合处理,并形成特征图像序列;
[0114]
S406、根据相位相关公式计算所述特征图像序列的交叉互功率谱,并根据所述交叉互功率谱的峰值位置选择自适应滤波器;
[0115]
S407、根据所述自适应滤波器对所述交叉互功率谱进行滤波,并对滤波后的交叉互功率谱进行反傅里叶变换;
[0116]
S408、根据反傅里叶变换后的交叉互功率谱提取所述特征图像序列中像素的振动信息,并根据所述振动信息确定所述待测轨道的实际振动数据;
[0117]
S409、根据所述实际振动数据确定所述待测轨道的振动部位,所述振动部位包括轨道扣件或轨道接头或道床;
[0118]
S410、采集所述待测轨道的轨道扣件或轨道接头的N段振动视频,并统计所述N段振动视频的平均振动频率、平均振幅和平均振动周期;
[0119]
S411、比较所述平均振动频率、所述平均振幅和所述平均振动周期与所述实际振动数 据,并根据比较结果确定所述轨道扣件或轨道接头的稳定性。
[0120]
可以看出,本申请实施例中所描绘的轨道振动检测方法,应用于轨道振动设备。所述轨道振动设备通过分区灰度值筛选策略对轨道振动视频进行处理,并从处理后的轨道振动视频中提取振动数据,以提升了振动数据提取的准确性,进而提高轨道振动检测的效率、准确和可靠性。
[0121]
与上述图2a和图3所述的实施例一致,请参阅图5。图5是本申请实施例提供的一种振动检测设备500的结构示意图,所述振动检测设备500包括应用处理器510、存储器520、通信接口530以及一个或多个程序521,其中,所述一个或多个程序521被存储在上述存储器520中,并且被配置由上述应用处理器510执行,所述一个或多个程序521包括用于执行以下步骤的指令:获取待测轨道的轨道振动视频,并根据所述待测轨道的轨道振动视频选取第一RGB图像;以及将所述第一RGB图像由RGB颜色空间线性转换到YIQ颜色空间,获得第一YIQ图像;以及通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列;以及根据所述特征图像序列提取所述待测轨道的实际振动数据;以及根据所述实际振动数据确定所述待测轨道的稳定性。
[0122]
可以看出,振动检测设备通过获取待测轨道的轨道振动视频,并根据所述待测轨道的轨道振动视频选取第一RGB图像;将所述第一RGB图像由RGB颜色空间线性转换到YIQ颜色空间,获得第一YIQ图像;通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列;根据所述特征图像序列提取所述待测轨道的实际振动数据;根据所述实际振动数据确定所述待测轨道的稳定性,提升振动数据提取的准确性以及工程应用的广泛性,进而提升轨道振动检测的效率、准确和可靠性。
[0123]
在一个可能示例中,在所述通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列方面,所述程序中的指令具体用于执行以下操作:对所述第一YIQ图像中的Y通道图像进行傅里叶变换;以及对傅里叶变换后的Y通道图像进行图像金字塔分解;以及对图像金字塔分解后的Y通道图像进行归一化、时空滤波以及线性放大处理;以及将处理后的Y通道图像与所述YIQ图像中的I通道图像和Q通道图像进行合成,并形成特征图像序列。
[0124]
在一个可能示例中,在所述根据所述特征图像序列提取所述待测轨道的实际振动数据方面,所述程序中的指令具体用于执行以下操作:根据相位相关公式计算所述特征图像序列的交叉互功率谱;以及根据所述交叉互动功率谱提取所述特征图像序列中像素的振动信 息;以及根据所述振动信息确定所述待测轨道的实际振动数据。
[0125]
在一个可能示例中,在所述根据所述交叉互动功率谱提取所述特征图像序列中像素的振动信息方面,所述程序中的指令具体用于执行以下操作:根据所述交叉互功率谱的峰值位置选择自适应滤波器;以及根据所述自适应滤波器对所述交叉互功率谱进行滤波;以及对滤波后的交叉互功率谱进行反傅里叶变换;以及根据反傅里叶变换后的交叉互功率谱提取所述特征图像序列中像素的振动信息。
[0126]
在一个可能示例中,在所述根据所述实际振动数据确定所述待测轨道的稳定性方面,所述程序中的指令具体用于执行以下操作:根据所述实际振动数据确定所述待测轨道的振动部位,所述振动部位包括轨道扣件或轨道接头或道床;以及根据所述待测轨道的振动部位和所述实际振动数据确定所述待测轨道的稳定性。
[0127]
在一个可能示例中,所述振动部位包括轨道扣件或轨道接头时,在所述根据所述待测轨道的振动部位和所述实际振动数据确定所述待测轨道的稳定性方面,所述程序中的指令具体用于执行以下操作:采集所述待测轨道的轨道扣件或轨道接头的N段振动视频;以及统计所述N段振动视频的平均振动频率、平均振幅和平均振动周期;以及比较所述平均振动频率、所述平均振幅和所述平均振动周期与所述实际振动数据,并根据比较结果确定所述待测轨道的稳定性。
[0128]
在一个可能示例中,所述振动部位包括道床时,在所述根据所述待测轨道的振动部位和所述实际振动数据确定所述待测轨道的稳定性方面,所述程序中的指令具体用于执行以下操作:获取道床的道床系数;以及根据所述道床系数确定所述道床的振动刚度和振动弹性系数;以及比较所述振动刚度和所述振动弹性系数与所述实际振动数据,并根据比较结果确定所述待测轨道的稳定性。
[0129]
上述主要从方法侧执行过程的角度对本申请实施例的方案进行了介绍。可以理解的是,振动检测设备为了实现上述功能,其包含了执行各个功能相应的硬件结构和/或软件模块。本领域技术人员应该很容易意识到,结合本文中所提供的实施例描述的各示例的单元及算法步骤,本申请能够以硬件或硬件和计算机软件的结合形式来实现。某个功能究竟以硬件还是计算机软件驱动硬件的方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用使用不同方法来实现所描述的功能,但是这种实现不应认为超出本申请的范围。
[0130]
本申请实施例可以根据上述方法示例对振动检测设备进行功能单元的划分,例如,可 以对应各个功能划分各个功能单元,也可以将两个或两个以上的功能集成在一个处理单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。需要说明的是,本申请实施例中对单元的划分是示意性的,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式。
[0131]
图6是本申请实施例提供的一种轨道振动检测装置600的功能单元组成框图。所述轨道振动检测装置600应用于振动检测设备,所述装置包括处理单元601和通信单元602。
[0132]
其中,所述处理单元601,用于执行如上述方法实施例中的任一步骤,且在执行诸如发送等数据传输时,可选择的调用所述通信单元602来完成相应操作,下面进行详细说明。
[0133]
所述处理单元601具体用于:获取待测轨道的轨道振动视频,并根据所述待测轨道的轨道振动视频选取第一RGB图像;以及将所述第一RGB图像由RGB颜色空间线性转换到YIQ颜色空间,获得第一YIQ图像;以及通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列;以及根据所述特征图像序列提取所述待测轨道的实际振动数据;以及根据所述实际振动数据确定所述待测轨道的稳定性。
[0134]
可以看出,本申请实施例公开了一种轨道振动检测装置,通过获取待测轨道的轨道振动视频,并根据所述待测轨道的轨道振动视频选取第一RGB图像;将所述第一RGB图像由RGB颜色空间线性转换到YIQ颜色空间,获得第一YIQ图像;通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列;根据所述特征图像序列提取所述待测轨道的实际振动数据;根据所述实际振动数据确定所述待测轨道的稳定性,提升振动数据提取的准确性以及工程应用的广泛性,进而提升轨道振动检测的效率、准确和可靠性。
[0135]
在一个可能的示例中,在所述通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列方面,所述处理单元601具体用于:对所述第一YIQ图像中的Y通道图像进行傅里叶变换;以及对傅里叶变换后的Y通道图像进行图像金字塔分解;以及对图像金字塔分解后的Y通道图像进行归一化、时空滤波以及线性放大处理;以及将处理后的Y通道图像与所述YIQ图像中的I通道图像和Q通道图像进行合成,并形成特征图像序列。
[0136]
在一个可能的示例中,在所述根据所述特征图像序列提取所述待测轨道的实际振动数据方面,所述处理单元601具体用于:根据相位相关公式计算所述特征图像序列的交叉互功率谱;以及根据所述交叉互动功率谱提取所述特征图像序列中像素的振动信息;以及根 据所述振动信息确定所述待测轨道的实际振动数据。
[0137]
在一个可能的示例中,在所述根据所述交叉互动功率谱提取所述特征图像序列中像素的振动信息方面,所述处理单元601具体用于:根据所述交叉互功率谱的峰值位置选择自适应滤波器;以及根据所述自适应滤波器对所述交叉互功率谱进行滤波;以及对滤波后的交叉互功率谱进行反傅里叶变换;以及根据反傅里叶变换后的交叉互功率谱提取所述特征图像序列中像素的振动信息。
[0138]
在一个可能的示例中,在所述根据所述实际振动数据确定所述待测轨道的稳定性方面,所述处理单元601具体用于:根据所述实际振动数据确定所述待测轨道的振动部位,所述振动部位包括轨道扣件或轨道接头或道床;以及根据所述待测轨道的振动部位和所述实际振动数据确定所述待测轨道的稳定性。
[0139]
在一个可能的示例中,所述振动部位包括轨道扣件或轨道接头时,在所述根据所述待测轨道的振动部位和所述实际振动数据确定所述待测轨道的稳定性方面,所述处理单元601具体用于:采集所述待测轨道的轨道扣件或轨道接头的N段振动视频;以及统计所述N段振动视频的平均振动频率、平均振幅和平均振动周期;以及比较所述平均振动频率、所述平均振幅和所述平均振动周期与所述实际振动数据,并根据比较结果确定所述待测轨道的稳定性。
[0140]
在一个可能的示例中,所述振动部位包括道床时,在所述根据所述待测轨道的振动部位和所述实际振动数据确定所述待测轨道的稳定性方面,所述处理单元601具体用于:获取道床的道床系数;以及根据所述道床系数确定所述道床的振动刚度和振动弹性系数;以及比较所述振动刚度和所述振动弹性系数与所述实际振动数据,并根据比较结果确定所述待测轨道的稳定性。
[0141]
本申请实施例还提供一种计算机存储介质,其中,该计算机存储介质存储用于电子数据交换的计算机程序,该计算机程序使得计算机执行如上述方法实施例中记载的任一方法的部分或全部步骤,上述计算机包括振动检测设备。
[0142]
本申请实施例还提供一种计算机程序产品,上述计算机程序产品包括存储了计算机程序的非瞬时性计算机可读存储介质,上述计算机程序可操作来使计算机执行如上述方法实施例中记载的任一方法的部分或全部步骤。该计算机程序产品可以为一个软件安装包,上述计算机包括振动检测设备。
[0143]
需要说明的是,对于前述的各方法实施例,为了简单描述,故将其都表述为一系列的 动作组合,但是本领域技术人员应该知悉,本申请并不受所描述的动作顺序的限制,因为依据本申请,某些步骤可以采用其他顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作和模块并不一定是本申请所必须的。
[0144]
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。
[0145]
在本申请所提供的几个实施例中,应该理解到,所揭露的装置,可通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如上述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性或其它的形式。
[0146]
上述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
[0147]
另外,在本申请各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
[0148]
上述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储器中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储器中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本申请各个实施例上述方法的全部或部分步骤。而前述的存储器包括:U盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
[0149]
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读存储器中,存储器可以包括:闪存盘、只读存储器(英文:Read-Only Memory,简称:ROM)、随机存取器(英文: Random Access Memory,简称:RAM)、磁盘或光盘等。
[0150]
以上对本申请实施例进行了详细介绍,本文中应用了具体个例对本申请的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本申请的方法及其核心思想;同时,对于本领域的一般技术人员,依据本申请的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本申请的限制。

权利要求书

[权利要求 1]
一种轨道振动检测方法,其特征在于,应用于振动检测设备,所述方法包括: 获取待测轨道的轨道振动视频,并根据所述待测轨道的轨道振动视频选取第一RGB图像; 将所述第一RGB图像由RGB颜色空间线性转换到YIQ颜色空间,获得第一YIQ图像; 通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列; 根据所述特征图像序列提取所述待测轨道的实际振动数据; 根据所述实际振动数据确定所述待测轨道的稳定性。
[权利要求 2]
根据权利要求1所述的方法,其特征在于,所述通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列,包括: 对所述第一YIQ图像中的Y通道图像进行傅里叶变换; 对傅里叶变换后的Y通道图像进行图像金字塔分解; 对图像金字塔分解后的Y通道图像进行归一化、时空滤波以及线性放大处理; 将处理后的Y通道图像与所述YIQ图像中的I通道图像和Q通道图像进行合成,并形成特征图像序列。
[权利要求 3]
根据权利要求1所述的方法,其特征在于,所述根据所述特征图像序列提取所述待测轨道的实际振动数据,包括: 根据相位相关公式计算所述特征图像序列的交叉互功率谱; 根据所述交叉互动功率谱提取所述特征图像序列中像素的振动信息; 根据所述振动信息确定所述待测轨道的实际振动数据。
[权利要求 4]
根据权利要求3所述的方法,其特征在于,所述根据所述交叉互动功率谱提取所述特征图像序列中像素的振动信息,包括: 根据所述交叉互功率谱的峰值位置选择自适应滤波器; 根据所述自适应滤波器对所述交叉互功率谱进行滤波; 对滤波后的交叉互功率谱进行反傅里叶变换; 根据反傅里叶变换后的交叉互功率谱提取所述特征图像序列中像素的振动信息。
[权利要求 5]
根据权利要求1所述的方法,其特征在于,所述根据所述实际振动数据确定所述待测轨道的稳定性,包括: 根据所述实际振动数据确定所述待测轨道的振动部位,所述振动部位包括轨道扣件或 轨道接头或道床; 根据所述待测轨道的振动部位和所述实际振动数据确定所述待测轨道的稳定性。
[权利要求 6]
根据权利要求5所述的方法,其特征在于,所述振动部位包括轨道扣件或轨道接头时,所述根据所述待测轨道的振动部位和所述实际振动数据确定所述待测轨道的稳定性,包括: 采集所述待测轨道的轨道扣件或轨道接头的N段振动视频; 统计所述N段振动视频的平均振动频率、平均振幅和平均振动周期; 比较所述平均振动频率、所述平均振幅和所述平均振动周期与所述实际振动数据,并根据比较结果确定所述待测轨道的稳定性。
[权利要求 7]
根据权利要求5所述的方法,其特征在于,所述振动部位包括道床时,所述根据所述待测轨道的振动部位和所述实际振动数据确定所述待测轨道的稳定性,包括: 获取道床的道床系数; 根据所述道床系数确定所述道床的振动刚度和振动弹性系数; 比较所述振动刚度和所述振动弹性系数与所述实际振动数据,并根据比较结果确定所述待测轨道的稳定性。
[权利要求 8]
一种轨道振动检测装置,其特征在于,应用于振动检测设备,所述装置包括处理单元和通信单元,其中: 所述处理单元,用于获取待测轨道的轨道振动视频,并根据所述待测轨道的轨道振动视频选取第一RGB图像;以及用于将所述第一RGB图像由RGB颜色空间线性转换到YIQ颜色空间,获得第一YIQ图像;以及用于通过视频放大算法对所述第一YIQ图像中的Y通道图像进行处理,获得特征图像序列;以及用于根据所述特征图像序列提取所述待测轨道的实际振动数据;以及用于根据所述实际振动数据确定所述待测轨道的稳定性。
[权利要求 9]
一种振动检测设备,其特征在于,包括处理器、存储器,所述存储器用于存储一个或多个程序,并且被配置由所述处理器执行,所述程序包括用于执行如权利要求1-7任一项所述的方法中的步骤的指令。
[权利要求 10]
一种计算机可读存储介质,其特征在于,所述计算机可读存储介质用于存储计算机程序,所述计算机程序被处理器执行,以实现如权利要求1-7任一项所述的方法。

附图

[ 图 1]  
[ 图 2A]  
[ 图 2B]  
[ 图 2C]  
[ 图 3]  
[ 图 4]  
[ 图 5]  
[ 图 6]