CN110674697B - Filtering method, device and related product - Google Patents
Filtering method, device and related product Download PDFInfo
- Publication number
- CN110674697B CN110674697B CN201910815821.3A CN201910815821A CN110674697B CN 110674697 B CN110674697 B CN 110674697B CN 201910815821 A CN201910815821 A CN 201910815821A CN 110674697 B CN110674697 B CN 110674697B
- Authority
- CN
- China
- Prior art keywords
- correlation peak
- filtering
- correlation
- peak
- video
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/40—Scenes; Scene-specific elements in video content
- G06V20/46—Extracting features or characteristics from the video content, e.g. video fingerprints, representative shots or key frames
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
- G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Multimedia (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
Description
技术领域technical field
本发明涉及检测技术领域,尤其涉及了一种滤波方法、装置及相关产品。The invention relates to the technical field of detection, in particular to a filtering method, device and related products.
背景技术Background technique
在相关技术中,振动可以反应某些机械结构的运行状况,可以拍摄振动装置运行时的视频,然后从视频中提取出振动信号,再根据振动信号获得装置的运行状况,但是由于拍摄视频时,会受环境影响而引入各种各样的噪声信息,从视频中提取出的振动信号也包含有大量的噪声信号,这样导致无法准确地获得装置的运行状况。In related technologies, vibration can reflect the operating conditions of certain mechanical structures. It is possible to take a video of the vibration device during operation, then extract the vibration signal from the video, and then obtain the operating status of the device according to the vibration signal. However, when shooting the video, A variety of noise information will be introduced due to the influence of the environment, and the vibration signal extracted from the video also contains a large number of noise signals, which makes it impossible to accurately obtain the operating status of the device.
发明内容Contents of the invention
本发明实施例提供一种滤波方法、装置及相关产品。Embodiments of the present invention provide a filtering method, device and related products.
第一方面,本发明实施例提供了一种滤波方法,包括:In a first aspect, an embodiment of the present invention provides a filtering method, including:
获取对检测视频对应的帧序列进行相位相关计算得到的一个或多个相关峰,所述相关峰为频域信号;Obtaining one or more correlation peaks obtained by performing phase correlation calculation on the frame sequence corresponding to the detected video, where the correlation peaks are frequency domain signals;
根据各相关峰在所述检测视频中对应的位置及相关峰的频段确定各相关峰对应的滤波策略;Determine the filtering strategy corresponding to each correlation peak according to the corresponding position of each correlation peak in the detection video and the frequency band of the correlation peak;
按照各相关峰对应的滤波策略对各相关峰进行滤波处理。Each correlation peak is filtered according to the filtering strategy corresponding to each correlation peak.
在某些实施例中,所述滤波策略包括滤波带宽,所述根据各相关峰在所述检测视频中对应的位置及所述相关峰对应的频段确定各相关峰对应的滤波策略包括:In some embodiments, the filtering strategy includes a filtering bandwidth, and determining the filtering strategy corresponding to each correlation peak according to the corresponding position of each correlation peak in the detected video and the frequency band corresponding to the correlation peak includes:
根据各相关峰对应的频段确定各相关峰对应的初始滤波带宽;Determine the initial filtering bandwidth corresponding to each correlation peak according to the frequency band corresponding to each correlation peak;
根据各相关峰在所述检测视频中对应的位置确定所述相关峰对应的带宽调整系数;Determine the bandwidth adjustment coefficient corresponding to the correlation peak according to the corresponding position of each correlation peak in the detection video;
将各相关峰对应的初始滤波带宽和带宽调整系数的乘积作为各相关峰对应的滤波带宽。The product of the initial filter bandwidth corresponding to each correlation peak and the bandwidth adjustment coefficient is used as the filter bandwidth corresponding to each correlation peak.
在某些实施例中,所述根据各相关峰在所述检测视频中对应的位置确定所述相关峰对应的带宽调整系数包括:In some embodiments, the determining the bandwidth adjustment coefficient corresponding to the correlation peak according to the corresponding position of each correlation peak in the detection video includes:
获取各相关峰在所述检测视频中对应的位置所在的区域对应的预设带宽调整系数,并将所述预设带宽调整系数作为所述相关峰对应的带宽调整系数。Obtain the preset bandwidth adjustment coefficient corresponding to the region where the corresponding position of each correlation peak in the detected video is located, and use the preset bandwidth adjustment coefficient as the bandwidth adjustment coefficient corresponding to the correlation peak.
在某些实施例中,所述根据各相关峰在所述检测视频中对应的位置及相关峰对应的频段确定各相关峰对应的滤波策略之前,所述滤波方法还包括:In some embodiments, before determining the filtering strategy corresponding to each correlation peak according to the corresponding position of each correlation peak in the detection video and the frequency band corresponding to the correlation peak, the filtering method further includes:
利用主成分分解方法对各相关峰进行降维处理;Use the principal component decomposition method to reduce the dimensionality of each correlation peak;
所述根据各相关峰在所述检测视频中对应的位置及相关峰对应的频段确定各相关峰对应的滤波策略包括:According to the corresponding position of each correlation peak in the detection video and the frequency band corresponding to the correlation peak, the filtering strategy corresponding to each correlation peak is determined to include:
根据各降维处理后的相关峰在所述检测视频中对应的位置及相关峰对应的频段确定各相关峰对应的滤波策略。The filtering strategy corresponding to each correlation peak is determined according to the corresponding position of each correlation peak after the dimensionality reduction processing in the detection video and the frequency band corresponding to the correlation peak.
在某些实施例中,所述利用主成分分解方法对各相关峰进行降维处理之前,所述滤波方法还包括:In some embodiments, before performing dimensionality reduction processing on each correlation peak by using the principal component decomposition method, the filtering method further includes:
将各相关峰的峰值与预设峰值范围比对,分别判断各相关峰的峰值是否在预设峰值范围内;Comparing the peak value of each correlation peak with the preset peak range, respectively judging whether the peak value of each correlation peak is within the preset peak range;
所述利用主成分分解方法对各相关峰进行降维处理包括:The method of utilizing the principal component decomposition method to carry out dimensionality reduction processing on each correlation peak includes:
对峰值在预设峰值范围内的相关峰进行降维处理。Dimensionality reduction is performed on related peaks whose peaks are within the preset peak range.
在某些实施例中,所述按照各相关峰对应的滤波策略对各相关峰进行滤波处理包括:In some embodiments, the filtering processing of each correlation peak according to the filtering strategy corresponding to each correlation peak includes:
按照各相关峰对应的滤波策略对各相关峰通过插值滤波的方式进行滤波处理。According to the filtering strategy corresponding to each correlation peak, filter processing is performed on each correlation peak by way of interpolation filtering.
第二方面,本发明还提供一种滤波装置,包括:In a second aspect, the present invention also provides a filtering device, including:
第一获取模块,用于获取对检测视频对应的帧序列进行相位相关计算得到的一个或多个相关峰,所述相关峰为频域信号;The first acquisition module is used to acquire one or more correlation peaks obtained by performing phase correlation calculation on the frame sequence corresponding to the detected video, and the correlation peaks are frequency domain signals;
滤波策略确定模块,用于根据各相关峰在所述检测视频中对应的位置及所述相关峰对应的频段确定各相关峰对应的滤波策略;A filter strategy determination module, configured to determine the filter strategy corresponding to each correlation peak according to the corresponding position of each correlation peak in the detection video and the frequency band corresponding to the correlation peak;
滤波模块,用于按照各相关峰对应的滤波策略对各相关峰进行滤波处理。The filtering module is configured to perform filtering processing on each correlation peak according to a filtering strategy corresponding to each correlation peak.
第三方面,本发明实施例还提供一种电子装置,包括处理器、存储器、以及存储在所述存储器上并可被所述处理器执行的信息推荐程序,其中所述信息推荐程序被所述处理器执行时,实现上述任一项实施例所述的滤波方法中的步骤的指令。In the third aspect, the embodiment of the present invention also provides an electronic device, including a processor, a memory, and an information recommendation program stored on the memory and executable by the processor, wherein the information recommendation program is executed by the When executed by a processor, the instructions to implement the steps in the filtering method described in any one of the above embodiments.
第四方面,本发明还提供一种计算机可读存储介质,计算机可读存储介质上存储有滤波程序,其中所述滤波程序被处理器执行时,实现上述任一项所述的滤波方法。In a fourth aspect, the present invention further provides a computer-readable storage medium, on which a filtering program is stored, wherein when the filtering program is executed by a processor, the filtering method described in any one of the above is implemented.
第五方面,本发明还提供一种检测设备,用于检测待检测装置的运行状况,所述检测设备包括上述实施例的滤波装置或上述实施例的电子装置。In a fifth aspect, the present invention further provides a detection device for detecting the operation status of the device to be detected, the detection device includes the filter device of the above embodiment or the electronic device of the above embodiment.
本发明实施例的滤波方法,获取对检测视频对应的帧序列进行相位相关计算得到的一个或多个相关峰;根据各相关峰在检测视频中对应的位置及相关峰对应的频段确定各相关峰对应的滤波策略;按照各相关峰对应的滤波策略对各相关峰进行滤波处理。通过根据各相关峰的位置及对应的频段选择合适的滤波策略进行滤波处理,可减少各相关峰中的噪声信息,使得从检测视频中获得的状态变化信号的信噪比更高,从而使得根据检测视频得到的待检测装置的运行状况更加准确。The filtering method of the embodiment of the present invention obtains one or more correlation peaks obtained by performing phase correlation calculation on the frame sequence corresponding to the detection video; and determines each correlation peak according to the corresponding position of each correlation peak in the detection video and the frequency band corresponding to the correlation peak Corresponding filtering strategy: filter each correlation peak according to the filtering strategy corresponding to each correlation peak. By selecting an appropriate filtering strategy for filtering according to the position of each correlation peak and the corresponding frequency band, the noise information in each correlation peak can be reduced, so that the signal-to-noise ratio of the state change signal obtained from the detection video is higher, so that according to The operation status of the device to be detected obtained by detecting the video is more accurate.
附图说明Description of drawings
下面将对本发明实施例涉及的一些附图进行说明。Some drawings related to the embodiments of the present invention will be described below.
图1为本发明实施例的电子装置的结构示意图;FIG. 1 is a schematic structural diagram of an electronic device according to an embodiment of the present invention;
图2为本发明实施例的滤波方法的流程示意图;FIG. 2 is a schematic flow chart of a filtering method according to an embodiment of the present invention;
图3为本发明实施例的滤波方法用于检测待检测装置的运行状况时涉及的放大检测视频的过程示意图;FIG. 3 is a schematic diagram of the process of zooming in and detecting video involved when the filtering method according to an embodiment of the present invention is used to detect the operating status of the device to be detected;
图4为本发明实施例的滤波方法的另一流程示意图;FIG. 4 is another schematic flowchart of a filtering method according to an embodiment of the present invention;
图5为本发明实施例的滤波方法的又一流程示意图;FIG. 5 is another schematic flowchart of a filtering method according to an embodiment of the present invention;
图6为本发明实施例的滤波方法的再一流程示意图。FIG. 6 is another schematic flowchart of a filtering method according to an embodiment of the present invention.
具体实施方式Detailed ways
下面结合本发明实施例中的附图对本发明实施例进行描述。Embodiments of the present invention will be described below with reference to the drawings in the embodiments of the present invention.
请参阅图1,图1是本发明实施例提供的电子装置100的硬件结构示意图。电子装置100包括处理器101、存储器102、输入输出接口103,以及一个或多个程序,一个或多个程序被存储在存储器102中,并且被配置由处理器101执行,程序包括用以下任一实施例的滤波方法的步骤的指令。电子装置100可以是服务器设备,也可以是终端设备。存储器102可以是高速RAM存储器,也可以是稳定的存储器(non-volatile memory),例如磁盘存储器,存储器102可选的还可以是独立于前述处理器101的存储装置。输入输出接口103可选的可以包括USB接口、标准的有线接口、无线接口(如WI-FI接口)等。Please refer to FIG. 1 . FIG. 1 is a schematic diagram of a hardware structure of an
请参阅图2,图2是本发明实施例提供的一种滤波方法的流程示意图,这种方法可包括但不限于如下步骤:Please refer to FIG. 2. FIG. 2 is a schematic flowchart of a filtering method provided by an embodiment of the present invention. This method may include but not limited to the following steps:
01、获取对检测视频对应的帧序列进行相位相关计算得到的一个或多个相关峰,相关峰为频域信号;01. Obtain one or more correlation peaks obtained by performing phase correlation calculation on the frame sequence corresponding to the detected video, and the correlation peaks are frequency domain signals;
本发明实施例的滤波方法可用于但不限于检测待检测装置的振动情况。待检测装置的振动情况可以反应待检测装置的运行状况。检测视频可为由成像模组拍摄待检测装置得到的视频。可通过分析检测视频,从检测视频中提取出待检测装置的振动信息,然后分析该振动信息得到待检测装置的运行状况。The filtering method of the embodiment of the present invention can be used for but not limited to detecting the vibration of the device to be detected. The vibration condition of the device to be detected can reflect the operation status of the device to be detected. The inspection video can be a video obtained by shooting the device to be inspected by the imaging module. By analyzing the detection video, the vibration information of the device to be detected can be extracted from the detection video, and then the vibration information can be analyzed to obtain the operation status of the device to be detected.
在一个优选实施例中,检测视频为将成像模组拍摄待检测装置得到的视频采用欧拉算法对视频中的振动进行放大之后的振动放大视频。待检测装置包括但不限于各类机械设备、建筑结构等。得到检测视频之后,对检测视频对应的帧序列进行相位相关计算得到交叉互功率谱,交叉互功率谱中包括一个或多个相关峰。可以理解,相关峰为频域信号,每个相关峰可以反应检测视频中的某个位置的某个频段的状态变化剧烈程度,对检测视频对应的帧序列进行相位相关计算得到交叉互功率谱可以理解为,从检测视频中提取出视频画面中的状态变化信息。状态变化信息包括振动信息和其他噪声信息,振动信息可以反应待检测装置的运行状况。可通过本发明实施例的滤波方法,对交叉互功率谱中的各相关峰进行滤波处理,以减少各相关峰中的噪声信息,以使得根据检测视频得到的待检测装置的运行状况更加准确。In a preferred embodiment, the detection video is a vibration amplified video after the vibration in the video is amplified by using the Euler algorithm to capture the video of the device to be detected by the imaging module. The devices to be tested include but are not limited to various types of mechanical equipment, building structures, etc. After the detected video is obtained, phase correlation calculation is performed on the frame sequence corresponding to the detected video to obtain a cross-cross power spectrum, and the cross-cross power spectrum includes one or more correlation peaks. It can be understood that the correlation peak is a frequency domain signal, and each correlation peak can reflect the severity of the state change of a certain frequency band at a certain position in the detection video, and the phase correlation calculation is performed on the frame sequence corresponding to the detection video to obtain the cross cross power spectrum. It is understood that the state change information in the video picture is extracted from the detected video. The status change information includes vibration information and other noise information, and the vibration information can reflect the operating status of the device to be detected. The filtering method in the embodiment of the present invention can be used to filter each correlation peak in the cross-cross power spectrum to reduce noise information in each correlation peak, so that the operating status of the device to be detected obtained from the detection video is more accurate.
02、根据各相关峰在检测视频中对应的位置及相关峰对应的频段确定各相关峰对应的滤波策略;02. Determine the filtering strategy corresponding to each correlation peak according to the corresponding position of each correlation peak in the detection video and the frequency band corresponding to the correlation peak;
可分析各相关峰得到各相关峰在检测视频中对应的位置,相关峰在检测视频中对应的位置可以理解为该相关峰所反映的状态变化对应的检测视频中的位置。例如,检测视频中有一块位置的亮度随着时间呈周期性变化,对检测视频的对应的帧序列进行相位相关计算得到的一个或多个相关峰中,则包含至少一个相关峰可以体现该位置的亮度变化。Each correlation peak can be analyzed to obtain the corresponding position of each correlation peak in the detection video, and the corresponding position of the correlation peak in the detection video can be understood as the position in the detection video corresponding to the state change reflected by the correlation peak. For example, the brightness of a position in the detection video changes periodically with time, and one or more correlation peaks obtained by performing phase correlation calculation on the corresponding frame sequence of the detection video include at least one correlation peak to reflect the position brightness changes.
由于检测视频中,不同位置的振动对检测待检测装置的运行状况的作用大小不一样。例如,处于检测视频的画面中的边缘位置的振动相比较于中间位置的振动,对检测待检测装置的运行状况的作用更小,那么可将各相关峰在检测视频中对应的位置作为依据,确定各相关峰的滤波策略。每个相关峰对应的频率范围(即频段)不一样,根据相关峰对应的频段确定各相关峰的滤波策略,这样使得滤波策略与每个相关峰更加适配,从而有助于从各相关峰中提取出有用的状态变化信号。滤波策略例如可包括但不限于滤波器的带宽,滤波器的类型等。In the detection video, vibrations at different locations have different effects on detecting the operating conditions of the device to be detected. For example, the vibration at the edge position in the picture of the detection video is less effective in detecting the operating condition of the device to be detected than the vibration at the middle position, so the corresponding positions of each correlation peak in the detection video can be used as the basis, Determine the filtering strategy for each correlation peak. The frequency range (ie, frequency band) corresponding to each correlation peak is different, and the filtering strategy of each correlation peak is determined according to the frequency band corresponding to the correlation peak, so that the filtering strategy is more suitable for each correlation peak, thus helping to obtain the corresponding frequency range from each correlation peak. Extract useful state change signals from . The filtering strategy may include, but not limited to, the bandwidth of the filter, the type of the filter, and the like, for example.
03、按照各相关峰对应的滤波策略对各相关峰进行滤波处理。03. Filter each correlation peak according to the filtering strategy corresponding to each correlation peak.
获得各相关峰的滤波策略之后,按照对应的滤波策略对各相关峰进行滤波处理,以使得从检测视频中获得的信号的信噪比更高。较佳地,可以按照各相关峰对应的滤波策略对各相关峰通过插值滤波的方式进行滤波处理,这样滤波效果更好。After the filtering strategy of each correlation peak is obtained, each correlation peak is filtered according to the corresponding filtering strategy, so that the signal-to-noise ratio of the signal obtained from the detection video is higher. Preferably, each correlation peak can be filtered by means of interpolation filtering according to the filtering strategy corresponding to each correlation peak, so that the filtering effect is better.
本发明实施例的滤波方法,获取对检测视频对应的帧序列进行相位相关计算得到的一个或多个相关峰;根据各相关峰在检测视频中对应的位置及相关峰对应的频段确定各相关峰对应的滤波策略;按照各相关峰对应的滤波策略对各相关峰进行滤波处理。通过根据各相关峰的位置及对应的频段选择合适的滤波策略进行滤波处理,可减少各相关峰中的噪声信息,使得从检测视频中获得的状态变化信号的信噪比更高,从而使得根据检测视频得到的待检测装置的运行状况更加准确。The filtering method of the embodiment of the present invention obtains one or more correlation peaks obtained by performing phase correlation calculation on the frame sequence corresponding to the detection video; and determines each correlation peak according to the corresponding position of each correlation peak in the detection video and the frequency band corresponding to the correlation peak Corresponding filtering strategy: filter each correlation peak according to the filtering strategy corresponding to each correlation peak. By selecting an appropriate filtering strategy for filtering according to the position of each correlation peak and the corresponding frequency band, the noise information in each correlation peak can be reduced, so that the signal-to-noise ratio of the state change signal obtained from the detection video is higher, so that according to The operation status of the device to be detected obtained by detecting the video is more accurate.
检测待检测装置的运行状况时,需要从包含有检测装置运行时的画面的检测视频中提取振动信息。如图3所示,可先将检测视频中的帧序列由RGB颜色空间转换到YIQ颜色空间,分离视频帧的亮度信息和色度信息。RGB和YIQ的转换关系为:When detecting the running status of the device to be detected, vibration information needs to be extracted from the detection video including the picture of the detection device running. As shown in Figure 3, the frame sequence in the detected video can be converted from the RGB color space to the YIQ color space, and the brightness information and chrominance information of the video frame can be separated. The conversion relationship between RGB and YIQ is:
Y=0.299*R+0.587*G+0.114*B;Y=0.299*R+0.587*G+0.114*B;
I=0.596*R–0.275*G–0.321*B;I=0.596*R-0.275*G-0.321*B;
Q=0.212*R-0.523*G+0.311*B。Q=0.212*R−0.523*G+0.311*B.
然后保持I、Q通道不变,对Y通道进行FFT操作,再利用欧拉运动放大算法对视频数据行放大处理。具体包括:将FFT变换后的Y通道图像进行复数可操纵金子塔空域分解,得到由多个不同空间分辨率的子图像组成的金字塔形结构;对金字塔型结构中的多个子图像中每个子图像进行时域带通滤波处理,得到目标频带对应的变换信号。可以理解,在视频画面中,振动可由视频帧序列的亮度反应,那么可通过分析检测视频中的Y通道获取检测视频中的振动信息。然后放大时域带通滤波后得到的目标频带对应的变换信号并提取感兴趣的运动信息,对感兴趣的运动信息进行复数可操纵金字塔重建,得到放大后的Y通道图像;最后将重建的Y通道图像与原来的I、Q通道图像相加,再转化为RGB色彩空间,得到输出视频。Then keep the I and Q channels unchanged, perform FFT operation on the Y channel, and then use the Euler motion amplification algorithm to amplify the video data. It specifically includes: decomposing the Y-channel image after FFT transformation into a complex number controllable pyramid space domain to obtain a pyramid structure composed of multiple sub-images with different spatial resolutions; for each sub-image in the multiple sub-images in the pyramid structure The time-domain band-pass filtering process is performed to obtain the transformed signal corresponding to the target frequency band. It can be understood that in a video picture, vibration can be reflected by the brightness of the video frame sequence, then the vibration information in the detected video can be obtained by analyzing the Y channel in the detected video. Then amplify the converted signal corresponding to the target frequency band obtained after time-domain bandpass filtering and extract the motion information of interest, perform complex manipulation pyramid reconstruction on the motion information of interest, and obtain the enlarged Y channel image; finally, the reconstructed Y The channel image is added to the original I and Q channel images, and then converted into RGB color space to obtain the output video.
再对视频运动放大处理后的帧序列采用相位相关算法计算帧序列间的交叉互功率谱。相位相关算法采用如下的公式计算交叉互功率谱。Then, the phase correlation algorithm is used to calculate the cross power spectrum between the frame sequences after video motion amplification processing. The phase correlation algorithm uses the following formula to calculate the cross-cross power spectrum.
上式中,Fa为a帧图像的傅立叶变换,为b帧图像的傅里叶变换的共轭信号,除式的下边为两个傅里叶变换的信号的相关积的模。R为本步骤的计算结果交叉互功率谱(包含频域噪音)。In the above formula, Fa is the Fourier transform of the frame a image, is the conjugate signal of the Fourier transform of the b-frame image, and the lower side of the division formula is the modulus of the correlation product of the two Fourier transformed signals. R is the cross-cross power spectrum (including frequency domain noise) of the calculation result of this step.
最后按照本发明实施例的上述滤波方法对交叉互功率谱中的各相关峰进行滤波处理,提高了有用信号的信噪比,使得从包含有检测装置运行时的画面的检测视频中提取出的信息,有用信号的信噪比更高。Finally, according to the above-mentioned filtering method of the embodiment of the present invention, the correlation peaks in the cross-cross power spectrum are filtered to improve the signal-to-noise ratio of the useful signal, so that Information, the signal-to-noise ratio of the useful signal is higher.
具体地,对目标振动视频的图像帧序列进行金字塔分解后,获得金字塔型结构。金字塔型结构包括多层子图像,且从上到下的图像分辨率依次降低,空间频率依次降低。对每层图像进行时域带通滤波处理,是为了得到目标频带,使得目标频带的子图像分辨率能够清晰表达图像的运动特征,同时不会因为分辨率过高造成计算量过大。因此,按照金字塔结构从下到上获取每个子图像在不同空间频率的频带,并且按照获取顺序与标准频带进行对比,当确定某一层子图像的频带与标准频带匹配成功时,就不需要再对该层子图像之上的一层或多层子图像进行空间频率的频带获取和匹配,提升了时域带通滤波的效率。另外,金字塔分解的最大层数确定:log2(min(xres,yres)),其中xres为图像的宽度像素值,yres为图像的高度像素值。Specifically, after performing pyramid decomposition on the image frame sequence of the target vibration video, a pyramid structure is obtained. The pyramid structure includes multiple sub-images, and the image resolution and spatial frequency decrease successively from top to bottom. The time-domain band-pass filtering process is performed on each layer of images to obtain the target frequency band, so that the sub-image resolution of the target frequency band can clearly express the motion characteristics of the image, and at the same time, the calculation amount will not be too large due to the high resolution. Therefore, according to the pyramid structure, the frequency bands of each sub-image at different spatial frequencies are obtained from bottom to top, and compared with the standard frequency bands in the order of acquisition. The frequency band acquisition and matching of the spatial frequency of one or more layers of sub-images above the layer of sub-images improves the efficiency of time-domain band-pass filtering. In addition, the maximum number of layers for pyramid decomposition is determined: log2(min(xres,yres)), where xres is the width pixel value of the image, and yres is the height pixel value of the image.
进一步地,对金字塔型结构中的多个子图像中每个子图像进行时域带通滤波处理,得到目标频带对应的变换信号,具体包括:Further, performing time-domain band-pass filtering processing on each sub-image in the multiple sub-images in the pyramid structure to obtain a transformed signal corresponding to the target frequency band, specifically including:
获取金字塔结构的总层数N,其中金字塔结构由下到上层数依次增大,对应的子图像分辨率依次升高;Obtain the total number of layers N of the pyramid structure, wherein the number of layers of the pyramid structure increases sequentially from bottom to top, and the resolution of the corresponding sub-images increases sequentially;
对金字塔结构中的多个子图像进行分组,其中第一组包括2n+1层对应的子图像,第二组包括2(n+1)层对应的子图像,其中n为大于或等于0的整数,N≥max(2(n+1),2n+1);Group multiple sub-images in the pyramid structure, where the first group includes sub-images corresponding to 2n+1 layers, and the second group includes sub-images corresponding to 2(n+1) layers, where n is an integer greater than or equal to 0 , N≥max(2(n+1),2n+1);
将第一组子图像采用第一处理器按照层数从小到大进行时域带通滤波处理,将第二组子图像采用第二处理器按照层数从小到大进行时域带通滤波处理,时域带通滤波处理为将子图像对应的不同空间频率的频带与标准频带进行匹配,第一处理器和第二处理器为独立运行的处理器;The first group of sub-images is processed by the first processor according to the number of layers from small to large, and the second group of sub-images is processed by the second processor according to the number of layers from small to large. The time-domain band-pass filtering process is to match the frequency bands of different spatial frequencies corresponding to the sub-images with the standard frequency bands, and the first processor and the second processor are independent processors;
当第一处理器或第二处理器确定频带与标准频带匹配成功时,确定频带为目标频带;When the first processor or the second processor determines that the frequency band matches the standard frequency band successfully, determine that the frequency band is the target frequency band;
获取目标频带对应的子图像为第一目标子图像,获取第一目标子图像上面一层子图像为第二目标子图像;Obtaining the sub-image corresponding to the target frequency band as the first target sub-image, and obtaining the sub-image above the first target sub-image as the second target sub-image;
确定第一目标子图像和第二目标子图像为目标频带对应的变换信号。Determine the first target sub-image and the second target sub-image as transformed signals corresponding to the target frequency band.
具体地,在对目标振动视频的多帧图像组成的帧序列进行空域金字塔分解后,可以获得金字塔结构的总层数N,且金字塔结构从下到上1~N层对应子图像的分辨率依次升高。然后对N层子图像进行分组,包括基数组和偶数组,基数组采用第一处理器按照1,3,5…2n+1的顺序进行时域带通滤波处理,偶数组采用第二处理器按照2,4,6,…2(n+1)的顺序进行时域带通滤波处理,两个处理器可以同时开始运行,也可以具有一定的处理时间间隔,另外,时域带通滤波处理是为了获取目标频带,获取方法可以是将子图像对应的不同空间频率的频带与标准频带进行匹配。当第一处理器或第二处理器确定频带与标准频带匹配成功时,确定所述频带为目标频带。Specifically, after performing spatial pyramid decomposition on the frame sequence composed of multiple frames of images of the target vibration video, the total number of layers N of the pyramid structure can be obtained, and the resolutions of the corresponding sub-images of the pyramid structure from bottom to top 1-N layers are sequentially raised. Then group the N-layer sub-images, including the base group and the even group. The base group uses the first processor to perform time-domain band-pass filtering in the order of 1, 3, 5...2n+1, and the even group uses the second processor. Perform time-domain band-pass filter processing in the order of 2, 4, 6, ... 2(n+1), the two processors can start running at the same time, or have a certain processing time interval. In addition, time-domain band-pass filter processing In order to acquire the target frequency band, the acquisition method may be to match the frequency bands of different spatial frequencies corresponding to the sub-images with the standard frequency band. When the first processor or the second processor determines that the frequency band matches the standard frequency band successfully, determine that the frequency band is the target frequency band.
因为目标频带表示能反应图像振动信息的最低分辨率子图像的空间频率,那么确定目标频带对应的子图像为第一目标子图像,然后获取目标频带对应子图像上一层的子图像作为第二目标子图像,将第一目标子图像和第二目标子图像为目标频带对应的变换信号,进行后续放大变换。这样,可以更准确地对目标振动视频对应的帧图像进行放大,获得更准确的运动信息,同时也避免了对更高分辨率图像进行放大的需要增加的运算量。Because the target frequency band represents the spatial frequency of the lowest resolution sub-image that can reflect image vibration information, then determine the sub-image corresponding to the target frequency band as the first target sub-image, and then obtain the sub-image above the sub-image corresponding to the target frequency band as the second For the target sub-image, the first target sub-image and the second target sub-image are used as transformed signals corresponding to the target frequency band, and subsequently enlarged and transformed. In this way, the frame image corresponding to the target vibration video can be amplified more accurately, and more accurate motion information can be obtained, and at the same time, the increased computational load required for amplifying a higher-resolution image can be avoided.
在本申请实施例中,对金字塔结构对应的子图像进行分组,然后通过两个独立运行的处理器对不同组的子图像进行时域带通滤波处理,可以提升滤波处理效率,同时在获取到目标频带对应层的子图像作为第一目标子图像后,获取其上一层的子图像作为第二目标子图像,然后将第一目标子图像和第二目标子图像为目标频带对应的变换信号,这样可以获得更准确的放大运动信息,同时也避免了对更高分辨率图像进行放大的需要增加的运算量。In the embodiment of the present application, the sub-images corresponding to the pyramid structure are grouped, and then the sub-images of different groups are processed by time-domain band-pass filtering through two independently operating processors, which can improve the filtering processing efficiency, and at the same time, the acquired After the sub-image of the layer corresponding to the target frequency band is used as the first target sub-image, the sub-image of the upper layer is obtained as the second target sub-image, and then the first target sub-image and the second target sub-image are transformed signals corresponding to the target frequency band , so that more accurate zoom-in motion information can be obtained, and at the same time, the increased computation required for zooming in on higher-resolution images can be avoided.
请参阅图4,基于上述实施例,在某些实施例中,滤波策略包括滤波带宽,根据各相关峰在检测视频中对应的位置及相关峰对应的频段确定各相关峰对应的滤波策略包括但不限于以下步骤:Please refer to FIG. 4 , based on the above embodiments, in some embodiments, the filtering strategy includes filtering bandwidth, and the filtering strategy corresponding to each correlation peak is determined according to the corresponding position of each correlation peak in the detection video and the frequency band corresponding to the correlation peak, including but Not limited to the following steps:
021、根据各相关峰对应的频段确定各相关峰对应的初始滤波带宽;021. Determine the initial filter bandwidth corresponding to each correlation peak according to the frequency band corresponding to each correlation peak;
可以理解,对检测视频对应的帧序列进行相位相关计算得到的相关峰中,包括检测视频对应的画面中位置相近的多个像素点的相似的多个状态变化信息,从而可以反应多个像素点所在的位置的状态变化。那么检测视频对应的画面中,不同位置的状态变化情况不同,对应的相关峰的频率分布也不同。如果用相同的滤波带宽进行滤波,那么为了避免有用信号被过滤,需要将滤波带宽设为较大值,但是这样有会导致噪声信号无法被过滤。而本申请的滤波方法中,根据各相关峰对应的频段确定个相关峰对应的初始滤波带宽,相关峰对应的频段与初始滤波带宽正相关,相关峰对应的频段越大,对应的初始滤波带宽越宽,可使得滤波带宽更好地与各相关峰适配,既能够有效地过滤噪声信号,也有助于避免有用信号被过滤,从而可以使得过滤效果更好。It can be understood that the correlation peak obtained by performing phase correlation calculation on the frame sequence corresponding to the detected video includes similar multiple state change information of multiple pixels with similar positions in the picture corresponding to the detected video, so that multiple pixel points can be reflected The status of the location changes. Then, in the picture corresponding to the detection video, the state changes of different positions are different, and the frequency distributions of the corresponding correlation peaks are also different. If the same filter bandwidth is used for filtering, in order to avoid useful signals from being filtered, the filter bandwidth needs to be set to a larger value, but this may cause noise signals to be unable to be filtered. And in the filtering method of the present application, determine the corresponding initial filter bandwidth of each correlation peak according to the frequency band corresponding to each correlation peak, the frequency band corresponding to the correlation peak is positively correlated with the initial filter bandwidth, the larger the frequency band corresponding to the correlation peak, the corresponding initial filter bandwidth The wider it is, the better the filtering bandwidth can be adapted to each correlation peak, which can not only effectively filter the noise signal, but also help prevent the useful signal from being filtered, thus making the filtering effect better.
例如,可以将B*(1-a)作为初始滤波带宽,B为各相关峰对应的频段,a可以根据实际需求确定,当需要得到更好的去燥效果时,可将a设为较大值,当需要获取到更多的有用信号时,可将a设定为较小值。当然,在其他实施例中,初始滤波带宽的设定方法不限于上述举例,在此不做限定。For example, B*(1-a) can be used as the initial filter bandwidth, B is the frequency band corresponding to each correlation peak, and a can be determined according to actual needs. When a better denoising effect is required, a can be set to a larger value value, when more useful signals need to be acquired, a can be set to a smaller value. Of course, in other embodiments, the method for setting the initial filtering bandwidth is not limited to the above example, and no limitation is made here.
022、根据各相关峰在检测视频中对应的位置确定相关峰对应的带宽调整系数;022. Determine the bandwidth adjustment coefficient corresponding to the correlation peak according to the corresponding position of each correlation peak in the detection video;
检测视频中不同位置的振动对检测待检测装置的运行状况的作用大小不一样。那么可根据各相关峰在检测视频中对应的位置确定相关峰对应的带宽调整系数。利用带宽调整系数来调整初始滤波带宽,使得滤波带宽更合理。Detecting vibrations at different positions in the video has different effects on detecting the operating status of the device to be detected. Then, the bandwidth adjustment coefficient corresponding to the correlation peak can be determined according to the corresponding position of each correlation peak in the detection video. The bandwidth adjustment coefficient is used to adjust the initial filtering bandwidth to make the filtering bandwidth more reasonable.
例如,可根据相关峰在检测视频中对应的位置的振动对检测待检测装置的运行状况的作用大小,将检测视频对应的画面分为多个区域,并为每个区域设定预设带宽调整系数,这样在步骤022中确定带宽调整系数时,可直接获取各相关峰在检测视频中对应的位置所在的区域对应的预设带宽调整系数,并将预设带宽调整系数作为相关峰对应的带宽调整系数。这样可以快速准确的获得带宽调整系数。For example, according to the effect of the vibration of the position corresponding to the correlation peak in the detection video on the detection of the operating status of the device to be detected, the screen corresponding to the detection video can be divided into multiple regions, and the preset bandwidth adjustment can be set for each region coefficient, when determining the bandwidth adjustment coefficient in
进一步地,用户可根据待检测装置的结构,对检测视频对应的画面中的各个区域设定带宽加权系数,然后将预设带宽调整系数与带宽加权系数的乘积作为对应的区域的相关的带宽调整系数,这样使得得到的各个相关峰的滤波带宽更加合理。例如,待检测装置中某些关键部位的振动相对其他部位更加能够反应待检测装置的运行状况,那么可将该关键部位在检测视频的画面中对应的区域的带宽加权系数设定为较大值,使得滤波带宽更大,从而可以提取到更多的有用信息。Further, according to the structure of the device to be detected, the user can set the bandwidth weighting coefficient for each area in the screen corresponding to the detection video, and then use the product of the preset bandwidth adjustment coefficient and the bandwidth weighting coefficient as the relevant bandwidth adjustment of the corresponding area coefficient, which makes the filtering bandwidth of each correlation peak more reasonable. For example, the vibration of some key parts in the device to be detected is more able to reflect the operating conditions of the device to be detected than other parts, then the bandwidth weighting coefficient of the corresponding area of the key part in the detection video screen can be set to a larger value , so that the filtering bandwidth is larger, so that more useful information can be extracted.
023、将各相关峰对应的初始滤波带宽和带宽调整系数的乘积作为各相关峰对应的滤波带宽。023. Use the product of the initial filtering bandwidth corresponding to each correlation peak and the bandwidth adjustment coefficient as the filtering bandwidth corresponding to each correlation peak.
这样得到的滤波带宽,综合考虑了相关峰在检测视频中对应的位置及相关峰对应的频段,使得滤波带宽与每个相关峰更加适配,从而有助于从各相关峰中提取出有用的状态变化信号。The filtering bandwidth obtained in this way comprehensively considers the corresponding position of the correlation peak in the detection video and the frequency band corresponding to the correlation peak, so that the filtering bandwidth is more suitable for each correlation peak, thereby helping to extract useful information from each correlation peak. State change signal.
需要说明的是,滤波带宽的计算方法不限于上述方法,在其他实施例中也可以根据各相关峰对应的频段和各相关峰在检测视频中对应的位置选择合适的计算方法来确定滤波带宽。It should be noted that the calculation method of the filter bandwidth is not limited to the above method, and in other embodiments, an appropriate calculation method may be selected according to the frequency band corresponding to each correlation peak and the corresponding position of each correlation peak in the detection video to determine the filter bandwidth.
请参阅图5,基于上述实施例,在某些实施例中,滤波方法还包括步骤:Please refer to FIG. 5, based on the above embodiments, in some embodiments, the filtering method further includes the steps of:
04、利用主成分分解方法对各相关峰进行降维处理;04. Use the principal component decomposition method to reduce the dimensionality of each correlation peak;
步骤04在步骤01之后,步骤02之前执行。主成分分解(PCA)用来将相关峰数据在振动检测的两个主方向上对数据维数约减。可以理解,从相关峰中获得振动信息时,对相关峰进行滤波的过程中,需要对各个相关峰进行计算,而从检测视频中提取出的相关峰包括检测视频对应的画面的多个维度的状态变化状况信息,这会导致滤波处理过程中计算繁琐,效率低。利用主成分分解方法将各相关峰降维处理之后,再执行其他关于滤波处理的计算,可大大降低滤波过程的计算量。
步骤02包括:
024、根据各降维处理后的相关峰在检测视频中对应的位置及相关峰对应的频段确定各相关峰对应的滤波策略。024. Determine the filtering strategy corresponding to each correlation peak according to the corresponding position of each correlation peak after dimensionality reduction processing in the detection video and the frequency band corresponding to the correlation peak.
将相关峰数据降维之后,再对降维处理后的相关峰进行计算,得到各相关峰对应的滤波策略,这样可降低滤波过程的计算量,提高滤波效率。After reducing the dimensionality of the correlation peak data, the correlation peaks after dimensionality reduction processing are calculated to obtain the filtering strategy corresponding to each correlation peak, which can reduce the calculation amount of the filtering process and improve the filtering efficiency.
进一步地,请参阅图6,在步骤04之前,滤波方法还包括:Further, referring to Fig. 6, before
05、将各相关峰的峰值与预设峰值范围比对,分别判断各相关峰的峰值是否在预设峰值范围内;05. Compare the peak value of each correlation peak with the preset peak range, and judge whether the peak value of each correlation peak is within the preset peak range;
步骤05在步骤01之后、步骤04之前执行。在成像模组拍摄待检测装置的检测视频时,一些像素点的微小变化也可以形成相关峰,但是频率的变化程度不够大,相关峰的峰值较小,这些信号对分析待检测装置的运行状况的作用很小甚至是没有作用的,可以理解为噪声信号。那么可通过预先设定一个预设峰值范围,在滤波时将各相关峰的峰值与预设峰值范围比对,过滤掉一些常见的噪声信号。
步骤04包括:
041、对峰值在预设峰值范围内的相关峰进行降维处理。041. Perform dimensionality reduction processing on related peaks whose peak values are within the preset peak range.
对峰值在预设峰值范围内的相关峰进行降维处理,过滤不在预设峰值范围内的相关峰,这样一方面可以提高从检测视频中获得的状态变化信号的信噪比,有用信号更多,另一方面由于去除了部分的噪声信号,从而可以减少后续的过滤的计算量,提升过滤效率。Dimensionality reduction processing is performed on the correlation peaks whose peak value is within the preset peak range, and the correlation peaks not within the preset peak range are filtered, so that on the one hand, the signal-to-noise ratio of the state change signal obtained from the detection video can be improved, and there are more useful signals , on the other hand, since part of the noise signal is removed, the calculation amount of subsequent filtering can be reduced, and the filtering efficiency can be improved.
进一步地,待检测装置的振动是一种周期往复运动,振动引起的状态变化也是周期性的。很多噪声信息虽然会导致检测视频中各像素点的状态变化,但是,噪声引起的状态变化常常不是周期性的,例如在室外场景下,由于室外光纤的明暗变化导致的状态变化。而且根据待检测装置的振动分析待检测装置的运行状况时,呈周期性的振动才能用于反应待检测装置的运行状况,因为非周期性的振动常常是由外界环境导致的,而不是由待检测装置自身引起的,那这部分非周期性的信号也不能用于分析待检测装置的运行状况。那么,在按照各相关峰对应的滤波策略对各相关峰进行滤波处理之后,还可以对相关峰进行反傅里叶变换,将频域信号转换为时域信号,得到多个状态变化信号,然后去除状态变化信号中的非周期信号,以进一步地去除噪声,使得根据检测视频得到的待检测装置的运行状况更加准确。Further, the vibration of the device to be detected is a periodic reciprocating motion, and the state change caused by the vibration is also periodic. Although a lot of noise information can lead to the detection of the state changes of each pixel in the video, the state changes caused by noise are often not periodic, for example, in outdoor scenes, the state changes caused by the light and dark changes of outdoor optical fibers. Moreover, when analyzing the operating conditions of the device to be detected based on the vibration of the device to be detected, only periodic vibrations can be used to reflect the operating status of the device to be detected, because non-periodic vibrations are often caused by the external environment, not by the device to be detected. If it is caused by the detection device itself, this part of the aperiodic signal cannot be used to analyze the operating status of the device to be detected. Then, after filtering each correlation peak according to the filtering strategy corresponding to each correlation peak, the correlation peak can also be inversely Fourier transformed to convert the frequency domain signal into a time domain signal to obtain multiple state change signals, and then The aperiodic signal in the state change signal is removed to further remove the noise, so that the operation status of the device to be detected obtained according to the detection video is more accurate.
本发明实施例还提供一种滤波装置,包括:An embodiment of the present invention also provides a filtering device, including:
第一获取模块,用于获取对检测视频对应的帧序列进行相位相关计算得到的一个或多个相关峰,相关峰为频域信号;The first acquisition module is used to acquire one or more correlation peaks obtained by performing phase correlation calculation on the frame sequence corresponding to the detection video, and the correlation peaks are frequency domain signals;
滤波策略确定模块,用于根据各相关峰在检测视频中对应的位置及相关峰对应的频段确定各相关峰对应的滤波策略;The filter strategy determination module is used to determine the filter strategy corresponding to each correlation peak according to the corresponding position of each correlation peak in the detection video and the corresponding frequency band of the correlation peak;
滤波模块,用于按照各相关峰对应的滤波策略对各相关峰进行滤波处理。The filtering module is configured to perform filtering processing on each correlation peak according to a filtering strategy corresponding to each correlation peak.
本发明实施例的滤波装置,获取对检测视频对应的帧序列进行相位相关计算得到的一个或多个相关峰,相关峰为频域信号;根据各相关峰在检测视频中对应的位置及相关峰对应的频段确定各相关峰对应的滤波策略;按照各相关峰对应的滤波策略对各相关峰进行滤波处理。通过根据各相关峰的位置及对应的频段选择合适的滤波策略进行滤波处理,可减少各相关峰中的噪声信息,使得从检测视频中获得的状态变化信号的信噪比更高,从而使得根据检测视频得到的待检测装置的运行状况更加准确。需要说明的是,上述滤波方法各实施例的解释说明及技术效果也适用于本实施例的滤波装置,为避免冗余,在此不再赘述。The filter device of the embodiment of the present invention obtains one or more correlation peaks obtained by performing phase correlation calculation on the frame sequence corresponding to the detection video, and the correlation peaks are frequency domain signals; according to the corresponding positions and correlation peaks of each correlation peak in the detection video The corresponding frequency band determines the filtering strategy corresponding to each correlation peak; performs filtering processing on each correlation peak according to the filtering strategy corresponding to each correlation peak. By selecting an appropriate filtering strategy for filtering according to the position of each correlation peak and the corresponding frequency band, the noise information in each correlation peak can be reduced, so that the signal-to-noise ratio of the state change signal obtained from the detection video is higher, so that according to The operation status of the device to be detected obtained by detecting the video is more accurate. It should be noted that the explanations and technical effects of the above filtering method embodiments are also applicable to the filtering device of this embodiment, and will not be repeated here to avoid redundancy.
在某些实施例中,滤波策略包括滤波带宽滤波策略确定模块包括:In some embodiments, the filtering strategy includes a filtering bandwidth. The filtering strategy determination module includes:
初始滤波带宽确定单元,用于根据各相关峰对应的频段确定各相关峰对应的初始滤波带宽;An initial filter bandwidth determination unit, configured to determine the initial filter bandwidth corresponding to each correlation peak according to the frequency band corresponding to each correlation peak;
带宽调整系数确定单元,用于根据各相关峰在检测视频中对应的位置确定相关峰对应的带宽调整系数;A bandwidth adjustment coefficient determination unit is used to determine the bandwidth adjustment coefficient corresponding to the correlation peak according to the corresponding position of each correlation peak in the detection video;
滤波带宽确定单元,用于将各相关峰对应的初始滤波带宽和带宽调整系数的乘积作为各相关峰对应的滤波带宽。The filter bandwidth determination unit is configured to use the product of the initial filter bandwidth corresponding to each correlation peak and the bandwidth adjustment coefficient as the filter bandwidth corresponding to each correlation peak.
在某些实施例中,带宽调整系数确定单元用于获取各相关峰在检测视频中对应的位置所在的区域对应的预设带宽调整系数,并将预设带宽调整系数作为相关峰对应的带宽调整系数。In some embodiments, the bandwidth adjustment coefficient determination unit is configured to obtain the preset bandwidth adjustment coefficient corresponding to the area where the corresponding position of each correlation peak in the detection video is located, and use the preset bandwidth adjustment coefficient as the bandwidth adjustment corresponding to the correlation peak coefficient.
在某些实施例中,滤波装置还包括:In some embodiments, the filtering device also includes:
降维模块,用于利用主成分分解方法对各相关峰进行降维处理;A dimensionality reduction module is used to perform dimensionality reduction processing on each correlation peak by using a principal component decomposition method;
滤波策略确定模块用于根据各降维处理后的相关峰在检测视频中对应的位置及相关峰对应的频段确定各相关峰对应的滤波策略。The filtering strategy determination module is used to determine the corresponding filtering strategy of each correlation peak according to the corresponding position in the detection video and the corresponding frequency band of each correlation peak after dimensionality reduction processing.
在某些实施例中,滤波装置还包括:In some embodiments, the filtering device also includes:
判断模块,用于将各相关峰的峰值与预设峰值范围比对,分别判断各相关峰的峰值是否在预设峰值范围内;A judging module, configured to compare the peak value of each correlation peak with a preset peak range, and determine whether the peak value of each correlation peak is within the preset peak range;
降维模块用于对峰值在预设峰值范围内的相关峰进行降维处理。The dimensionality reduction module is used to perform dimensionality reduction processing on related peaks whose peak value is within a preset peak range.
在某些实施例中,滤波模块用于按照各相关峰对应的滤波策略对各相关峰通过插值滤波的方式进行滤波处理。In some embodiments, the filtering module is configured to filter each correlation peak through interpolation filtering according to a filtering strategy corresponding to each correlation peak.
其中,上述滤波装置中各个模块的功能实现及技术效果与上述滤波方法实施例中各步骤相对应,其功能和实现过程在此处不再一一赘述。Wherein, the function implementation and technical effect of each module in the above filtering device correspond to each step in the above filtering method embodiment, and the functions and implementation processes thereof will not be repeated here.
本发明还提供一种计算机可读存储介质,计算机可读存储介质上存储有滤波程序,其中滤波程序被处理器执行时,实现上述任一实施例的滤波方法的步骤。The present invention also provides a computer-readable storage medium, on which a filtering program is stored, wherein when the filtering program is executed by a processor, the steps of the filtering method in any of the foregoing embodiments are implemented.
其中,滤波程序被执行时所实现的方法及对应的技术效果可参照本发明滤波方法的各个实施例,此处不再赘述。Wherein, the method and the corresponding technical effect realized when the filtering program is executed can refer to the various embodiments of the filtering method of the present invention, which will not be repeated here.
本发明实施例还提供一种检测设备,用于检测待检测装置的运行状况,检测设备包括上述实施例的滤波装置或上述实施例的电子装置。An embodiment of the present invention also provides a detection device for detecting the operation status of the device to be detected, and the detection device includes the filter device of the above embodiment or the electronic device of the above embodiment.
本发明实施例的检测设备,获取对检测视频对应的帧序列进行相位相关计算得到的一个或多个相关峰;根据各相关峰在检测视频中对应的位置及相关峰对应的频段确定各相关峰对应的滤波策略;按照各相关峰对应的滤波策略对各相关峰进行滤波处理。通过根据各相关峰的位置及对应的频段选择合适的滤波策略进行滤波处理,可减少各相关峰中的噪声信息,使得从检测视频中获得的状态变化信号的信噪比更高,从而使得根据检测视频得到的待检测装置的运行状况更加准确。The detection device in the embodiment of the present invention obtains one or more correlation peaks obtained by performing phase correlation calculation on the frame sequence corresponding to the detection video; and determines each correlation peak according to the corresponding position of each correlation peak in the detection video and the frequency band corresponding to the correlation peak Corresponding filtering strategy: filter each correlation peak according to the filtering strategy corresponding to each correlation peak. By selecting an appropriate filtering strategy for filtering according to the position of each correlation peak and the corresponding frequency band, the noise information in each correlation peak can be reduced, so that the signal-to-noise ratio of the state change signal obtained from the detection video is higher, so that according to The operation status of the device to be detected obtained by detecting the video is more accurate.
检测设备检测待检测装置的运行状况时,可利用成像模组拍摄待检测装置运行时的检测视频,然后从包含有检测装置运行时的画面的检测视频中提取状态变化信息,状态变化信息可以反应待检测装置的振动信息,可以根据待检测装置的振动信息获得待检测装置的运行状况。如图3所示,可先将检测视频中的帧序列由RGB颜色空间转换到YIQ颜色空间,分离视频帧的亮度信息和色度信息。RGB和YIQ的转换关系为:When the detection equipment detects the operating status of the device to be detected, the imaging module can be used to shoot the detection video of the device to be detected when it is running, and then the state change information can be extracted from the detection video containing the picture of the detection device running. The state change information can reflect The vibration information of the device to be detected can obtain the operation status of the device to be detected according to the vibration information of the device to be detected. As shown in Figure 3, the frame sequence in the detected video can be converted from the RGB color space to the YIQ color space, and the brightness information and chrominance information of the video frame can be separated. The conversion relationship between RGB and YIQ is:
Y=0.299*R+0.587*G+0.114*B;Y=0.299*R+0.587*G+0.114*B;
I=0.596*R–0.275*G–0.321*B;I=0.596*R-0.275*G-0.321*B;
Q=0.212*R-0.523*G+0.311*B。Q=0.212*R−0.523*G+0.311*B.
然后保持I、Q通道不变,对Y通道进行FFT操作,再利用欧拉运动放大算法对视频数据行放大处理。具体包括:将FFT变换后的Y通道图像进行复数可操纵金子塔空域分解。将Y通道空域分解后的不同尺度的图像进行时域带通滤波,可以理解,在视频画面中,振动可由视频帧序列的亮度反应,那么可通过分析检测视频中的Y通道获取检测视频中的振动信息。;然后放大时域带通滤波后感兴趣的运动信息,对感兴趣的运动信息进行复数可操纵金字塔重建,得到放大后的Y通道图像;最后将重建的Y通道图像与原来的I、Q通道图像相加,再转化为RGB色彩空间,得到输出视频。Then keep the I and Q channels unchanged, perform FFT operation on the Y channel, and then use the Euler motion amplification algorithm to amplify the video data. It specifically includes: decomposing the Y-channel image after the FFT transformation into a complex number steerable pyramid space domain. The images of different scales decomposed in the Y channel space domain are subjected to time domain bandpass filtering. It can be understood that in the video picture, the vibration can be reflected by the brightness of the video frame sequence, so the Y channel in the detection video can be obtained by analyzing and detecting the Y channel in the video. vibration information. ; Then amplify the motion information of interest after time-domain bandpass filtering, perform complex manipulation pyramid reconstruction on the motion information of interest, and obtain an enlarged Y channel image; finally combine the reconstructed Y channel image with the original I, Q channel The images are added and converted to RGB color space to obtain the output video.
再对视频运动放大处理后的帧序列采用相位相关算法计算帧序列间的交叉互功率谱。相位相关算法采用如下的公式计算交叉互功率谱。Then, the phase correlation algorithm is used to calculate the cross power spectrum between the frame sequences after video motion amplification processing. The phase correlation algorithm uses the following formula to calculate the cross-cross power spectrum.
上式中,Fa为a帧图像的傅立叶变换,为b帧图像的傅里叶变换的共轭信号,除式的下边为两个傅里叶变换的信号的相关积的模。R为本步骤的计算结果交叉互功率谱(包含频域噪音)。交叉互功率谱中包括一个或多个相关峰,相关峰为频域信号。可以理解,每个相关峰可以反应检测视频中的某个位置的状态变化,对检测视频对应的真序列进行相位相关计算得到交叉互功率谱可以理解为,从检测视频中提取出视频画面中的状态变化信息。状态变化信息包括振动信息和其他噪声信息。可通过本发明实施例的滤波方法,对交叉互功率谱中的各相关峰进行滤波处理,以减少各相关峰中的噪声信息,使得从包含有检测装置运行时的画面的检测视频中提取出的信息中,有用信号的信噪比更高。In the above formula, Fa is the Fourier transform of the frame a image, is the conjugate signal of the Fourier transform of the b-frame image, and the lower side of the division formula is the modulus of the correlation product of the two Fourier transformed signals. R is the cross-cross power spectrum (including frequency domain noise) of the calculation result of this step. The cross-cross power spectrum includes one or more correlation peaks, and the correlation peaks are frequency-domain signals. It can be understood that each correlation peak can reflect the state change of a certain position in the detection video, and the phase correlation calculation of the true sequence corresponding to the detection video to obtain the cross-cross power spectrum can be understood as extracting the Status change information. State change information includes vibration information and other noise information. Through the filtering method of the embodiment of the present invention, each correlation peak in the cross-cross power spectrum can be filtered to reduce the noise information in each correlation peak, so that it can be extracted from the detection video containing the picture when the detection device is running. In the information, the signal-to-noise ratio of the useful signal is higher.
在上述实施例中,可全部或部分地通过软件、硬件、固件、或其任意组合来实现。当使用软件实现时,可以全部或部分地以计算机程序产品的形式实现。所述计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行所述计算机程序指令时,全部或部分地产生按照本发明实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中,或者从一个计算机可读存储介质向另一个计算机可读存储介质传输,例如,所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线)或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输。所述计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质(例如软盘、硬盘、磁带)、光介质(例如光盘)、或者半导体介质(例如固态硬盘)等。在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。In the above-mentioned embodiments, all or part may be implemented by software, hardware, firmware, or any combination thereof. When implemented using software, it may be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on the computer, the processes or functions according to the embodiments of the present invention will be generated in whole or in part. The computer can be a general purpose computer, a special purpose computer, a computer network, or other programmable devices. The computer instructions may be stored in or transmitted from one computer-readable storage medium to another computer-readable storage medium, for example, the computer instructions may be transmitted from a website, computer, server or data center Transmission to another website site, computer, server, or data center by wired (eg, coaxial cable, optical fiber, DSL) or wireless (eg, infrared, wireless, microwave, etc.) means. The computer-readable storage medium may be any available medium that can be accessed by a computer, or a data storage device such as a server or a data center integrated with one or more available media. The available medium may be a magnetic medium (such as a floppy disk, a hard disk, or a magnetic tape), an optical medium (such as an optical disk), or a semiconductor medium (such as a solid-state disk), and the like. In the foregoing embodiments, the descriptions of each embodiment have their own emphases, and for parts not described in detail in a certain embodiment, reference may be made to relevant descriptions of other embodiments.
在本发明所提供的几个实施例中,应该理解到,所揭露的装置,也可以通过其它的方式实现。例如以上所描述的装置实施例仅仅是示意性的,例如所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可结合或者可以集成到另一个系统,或一些特征可以忽略或不执行。另一点,所显示或讨论的相互之间的间接耦合或者直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性或其它的形式。In the several embodiments provided by the present invention, it should be understood that the disclosed device can also be implemented in other ways. For example, the device embodiments described above are only illustrative, for example, the division of the units is only a logical function division, and there may be other division methods in actual implementation, for example, multiple units or components can be combined or integrated to another system, or some features may be ignored or not implemented. In another point, the shown or discussed mutual indirect coupling or direct coupling or communication connection may be through some interfaces, and the indirect coupling or communication connection of devices or units may be in electrical or other forms.
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者,也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例的方案的目的。The units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or may be distributed to multiple network units . Part or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.
另外,在本发明各实施例中的各功能单元可集成在一个处理单元中,也可以是各单元单独物理存在,也可两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,或者也可以采用软件功能单元的形式实现。In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, each unit may exist separately physically, or two or more units may be integrated into one unit. The above-mentioned integrated units can be implemented in the form of hardware, or can also be implemented in the form of software functional units.
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质例如可包括:U盘、移动硬盘、只读存储器(ROM,Read-OnlyMemory)、随机存取存储器(RAM,Random Access Memory)、磁碟或光盘等各种可存储程序代码的介质。If the integrated unit is realized in the form of a software function unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the essence of the technical solution of the present invention or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium. Several instructions are included to make a computer device (which may be a personal computer, server or network device, etc.) execute all or part of the steps of the methods described in the various embodiments of the present invention. And the aforementioned storage medium can include, for example: U disk, mobile hard disk, read-only memory (ROM, Read-OnlyMemory), random access memory (RAM, Random Access Memory), magnetic disk or optical disk, etc., which can store program codes. medium.
Claims (9)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2020/105583 WO2021036661A1 (en) | 2019-04-26 | 2020-07-29 | Filtering method and apparatus, and related product |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2019103460141 | 2019-04-26 | ||
| CN201910346014 | 2019-04-26 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110674697A CN110674697A (en) | 2020-01-10 |
| CN110674697B true CN110674697B (en) | 2023-04-25 |
Family
ID=69075853
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910815821.3A Expired - Fee Related CN110674697B (en) | 2019-04-26 | 2019-08-30 | Filtering method, device and related product |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN110674697B (en) |
| WO (1) | WO2021036661A1 (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110674697B (en) * | 2019-04-26 | 2023-04-25 | 深圳市豪视智能科技有限公司 | Filtering method, device and related product |
| CN114504327B (en) * | 2021-12-28 | 2024-05-17 | 深圳大学 | Electroencephalogram noise processing method and device and computer equipment |
| CN117576091B (en) * | 2024-01-15 | 2024-04-09 | 深圳昱拓智能有限公司 | A cooling tower fan vibration detection method and system based on video detection |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2187059A (en) * | 1986-02-20 | 1987-08-26 | British Broadcasting Corp | Motion detector for video signal prediction systems |
| CN103714522A (en) * | 2014-01-06 | 2014-04-09 | 海信集团有限公司 | Image denoising method and image denoising system |
| CN107844779A (en) * | 2017-11-21 | 2018-03-27 | 重庆邮电大学 | A kind of video key frame extracting method |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103024248B (en) * | 2013-01-05 | 2016-01-06 | 上海富瀚微电子股份有限公司 | The video image noise reducing method of Motion Adaptive and device thereof |
| US10667723B2 (en) * | 2016-02-19 | 2020-06-02 | Covidien Lp | Systems and methods for video-based monitoring of vital signs |
| CN110659583B (en) * | 2019-04-26 | 2023-11-28 | 深圳市豪视智能科技有限公司 | Signal processing methods, devices and related products |
| CN110674697B (en) * | 2019-04-26 | 2023-04-25 | 深圳市豪视智能科技有限公司 | Filtering method, device and related product |
-
2019
- 2019-08-30 CN CN201910815821.3A patent/CN110674697B/en not_active Expired - Fee Related
-
2020
- 2020-07-29 WO PCT/CN2020/105583 patent/WO2021036661A1/en not_active Ceased
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2187059A (en) * | 1986-02-20 | 1987-08-26 | British Broadcasting Corp | Motion detector for video signal prediction systems |
| CN103714522A (en) * | 2014-01-06 | 2014-04-09 | 海信集团有限公司 | Image denoising method and image denoising system |
| CN107844779A (en) * | 2017-11-21 | 2018-03-27 | 重庆邮电大学 | A kind of video key frame extracting method |
Non-Patent Citations (1)
| Title |
|---|
| 微小运动放大及其在超声波传播中的应用;轩建平等;《华中科技大学学报(自然科学版)》;20180309(第02期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2021036661A1 (en) | 2021-03-04 |
| CN110674697A (en) | 2020-01-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110595745B (en) | Detection method for abnormality of fixing screw of equipment and related product | |
| US9324153B2 (en) | Depth measurement apparatus, image pickup apparatus, depth measurement method, and depth measurement program | |
| CN110659583B (en) | Signal processing methods, devices and related products | |
| CN110674697B (en) | Filtering method, device and related product | |
| CN101345891B (en) | No-reference image quality assessment method based on information entropy and contrast | |
| US11272146B1 (en) | Content adaptive lens shading correction method and apparatus | |
| JP5983373B2 (en) | Image processing apparatus, information processing method, and program | |
| US9153013B2 (en) | Image processing apparatus, image processing method and computer readable medium | |
| US10217193B2 (en) | Image processing apparatus, image capturing apparatus, and storage medium that stores image processing program | |
| US20070097236A1 (en) | Method and apparatus for reducing noise of image sensor | |
| Ponomarenko et al. | Adaptive visually lossless JPEG-based color image compression | |
| JP2015166979A5 (en) | Image processing method, image processing apparatus, imaging apparatus, image processing program, and recording medium | |
| JP2022130888A (en) | Image correction system, image correction method, and program | |
| CN109672829B (en) | Image brightness adjusting method and device, storage medium and terminal | |
| JP2014146872A5 (en) | ||
| WO2025025690A1 (en) | Method and apparatus for determining accuracy of vision system, and computer readable storage medium | |
| Temel et al. | BLeSS: Bio-inspired low-level spatiochromatic similarity assisted image quality assessment | |
| US20130108163A1 (en) | Image evaluation apparatus, image evaluation method, and program | |
| US10249020B2 (en) | Image processing unit, imaging device, computer-readable medium, and image processing method | |
| WO2016139871A1 (en) | Pulse measuring device and control method therefor | |
| Burns | Estimation error in image quality measurements | |
| US12106451B2 (en) | Block-based noise reduction intensity determination method and device | |
| CN103325091B (en) | Low-frequency spectra data padding method image acquiring method and system | |
| JP5836878B2 (en) | Image processing apparatus, method, and program | |
| Sadou et al. | No reference image quality assessment: Feature fusion using relevance vector machine |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | 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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240909 Address after: Room 21A, Junhao Ge, No. 2 Xiangxuan Road, Donghai Community, Xiangmihu Street, Futian District, Shenzhen City, Guangdong Province 518040 Patentee after: SHENZHEN GUANGNING INDUSTRIAL Co.,Ltd. Country or region after: China Address before: 518000 7th floor, building 9, phase II software park, Yuehai street, Nanshan District, Shenzhen City, Guangdong Province Patentee before: Shenzhen Haoxi Intelligent Technology Co.,Ltd. Country or region before: China Patentee before: SHENZHEN GUANGNING INDUSTRIAL Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20241224 Address after: 518000 floor 5, building C3, Nanshan Zhiyuan, No. 1001, Xueyuan Avenue, Nanshan District, Shenzhen, Guangdong Patentee after: SHENZHEN SECOM TECHNOLOGY Co.,Ltd. Country or region after: China Address before: Room 21A, Junhao Ge, No. 2 Xiangxuan Road, Donghai Community, Xiangmihu Street, Futian District, Shenzhen City, Guangdong Province 518040 Patentee before: SHENZHEN GUANGNING INDUSTRIAL Co.,Ltd. Country or region before: China |
|
| TR01 | Transfer of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20230425 |



