CN110659583B - 信号处理方法、装置及相关产品 - Google Patents
信号处理方法、装置及相关产品 Download PDFInfo
- Publication number
- CN110659583B CN110659583B CN201910813881.1A CN201910813881A CN110659583B CN 110659583 B CN110659583 B CN 110659583B CN 201910813881 A CN201910813881 A CN 201910813881A CN 110659583 B CN110659583 B CN 110659583B
- Authority
- CN
- China
- Prior art keywords
- state change
- signals
- signal
- cross
- power spectrum
- 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.)
- Active
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 30
- 230000008859 change Effects 0.000 claims abstract description 201
- 238000001514 detection method Methods 0.000 claims abstract description 83
- 238000001228 spectrum Methods 0.000 claims abstract description 70
- 230000000737 periodic effect Effects 0.000 claims abstract description 43
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 238000001914 filtration Methods 0.000 claims description 88
- 238000000034 method Methods 0.000 claims description 21
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000000875 corresponding effect Effects 0.000 description 129
- 230000000694 effects Effects 0.000 description 12
- 238000000354 decomposition reaction Methods 0.000 description 9
- 238000006243 chemical reaction Methods 0.000 description 8
- 238000004422 calculation algorithm Methods 0.000 description 7
- 230000009466 transformation Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000000605 extraction Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000003321 amplification Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000003199 nucleic acid amplification method Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000004590 computer program Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/40—Scenes; Scene-specific elements in video content
- G06V20/41—Higher-level, semantic clustering, classification or understanding of video scenes, e.g. detection, labelling or Markovian modelling of sport events or news items
- G06V20/42—Higher-level, semantic clustering, classification or understanding of video scenes, e.g. detection, labelling or Markovian modelling of sport events or news items of sport video content
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种信号处理方法、装置及相关产品。信号处理方法包括获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱对应的多个状态变化信号，所述状态变化信号为时域信号；分析所述多个状态变化信号，以获取所述多个状态变化信号中的非周期信号；去除所述多个状态变化信号中的非周期信号。这样从检测视频中获取到的状态变化信号中，有用信息更多，从而可以使得根据检测视频得到的待检测装置的运行状况更加准确。
Description
技术领域
本发明涉及检测技术领域，尤其涉及了一种信号处理方法、装置及相关产品。
背景技术
在相关技术中，振动可以反应某些机械结构的运行状况，可以拍摄振动装置运行时的视频，然后从视频中获取待检测装置的振动情况，提取振动信号时，可通过先提取视频中各个像素点的状态变化信号，再根据状态变化信号得到待检测装置的振动情况，但是由于拍摄视频时，会受环境影响而引入各种各样的噪声信息，例如光照也会导致视频中各个像素点的状态变化，得到的待检测装置的振动情况也包含噪声信息，这样导致无法准确地获得装置的运行状况。
发明内容
本发明实施例提供一种信号处理方法、装置及相关产品。
第一方面，本发明实施例提供了一种信号处理方法，包括：
获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱对应的多个状态变化信号，所述状态变化信号为时域信号；
分析所述多个状态变化信号，以获取所述多个状态变化信号中的非周期信号；
去除所述多个状态变化信号中的非周期信号。
在某些实施例中，所述分析所述多个状态变化信号，以获取所述多个状态变化信号中的非周期信号包括：
从第一状态变化信号段，并获取各第一状态变化信号段的频率；
将所述多个状态变化信号中，对应的多个第一状态变化信号段的频率不一致的状态变化信号作为所述非周期信号。
在某些实施例中，所述分析所述多个状态变化信号，以获取所述多个状态变化信号中的非周期信号包括：
从各状态变化信号中提取预设长度的目标状态变化信号段,并获取所述目标状态变化信号段的目标频率；
根据各状态变化信号对应的目标频率设定对应的滑窗,并将各状态变化信号发送至对应的滑窗；
将不能通过对应的滑窗的状态变化信号作为所述非周期信号。
在某些实施例中，所述获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱对应的多个状态变化信号包括：
获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱，所述交叉互功率谱为频域信号；
对所述交叉互功率谱进行反傅里叶变换得到多个状态变化信号。
在某些实施例中，所述对所述交叉互功率谱进行反傅里叶变换得到多个状态变化信号之前，所述获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱对应的多个状态变化信号还包括：
根据所述交叉互功率谱中的各相关峰在所述检测视频中对应的位置及相关峰的频段确定各相关峰对应的滤波策略；
按照各相关峰对应的滤波策略对所述交叉互功率谱中的各相关峰进行滤波处理；
所述对所述交叉互功率谱进行反傅里叶变换得到多个状态变化信号包括：
对进行滤波处理之后的交叉互功率谱进行反傅里叶变换得到多个状态变化信号。
在某些实施例中，所述滤波策略包括滤波带宽，所述根据所述交叉互功率谱中的各相关峰在所述检测视频中对应的位置及相关峰的频段确定各相关峰对应的滤波策略包括：
根据所述交叉互功率谱中的各相关峰对应的频段确定各相关峰对应的初始滤波带宽；
根据各相关峰在所述检测视频中对应的位置确定所述相关峰对应的带宽调整系数；
将各相关峰对应的初始滤波带宽和带宽调整系数的乘积作为各相关峰对应的滤波带宽。
第二方面，本发明实施例还提供一种信号处理装置，包括：
信号获取模块，用于获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱对应的多个状态变化信号，所述状态变化信号为时域信号；
分析模块，用于分析所述多个状态变化信号，以获取所述多个状态变化信号中的非周期信号；
信号处理模块，用于去除所述多个状态变化信号中的非周期信号。
第三方面，本发明实施例还提供一种电子装置，包括处理器、存储器、以及存储在所述存储器上并可被所述处理器执行的信息推荐程序，其中所述信息推荐程序被所述处理器执行时，实现上述任一项实施例所述的信号处理方法中的步骤的指令。
第四方面，本发明还提供一种计算机可读存储介质，计算机可读存储介质上存储有信号处理程序，其中所述信号处理程序被处理器执行时，实现上述任一项所述的信号处理方法。
第五方面，本发明还提供一种检测设备，用于检测待检测装置的运行状况，所述检测设备包括上述实施例的信号处理装置或上述实施例的电子装置。
本发明实施例的技术方案，获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱对应的多个状态变化信号，状态变化信号为时域信号；分析多个状态变化信号，以获取多个状态变化信号中的非周期信号；去除多个状态变化信号中的非周期信号。这样从检测视频中获取到的状态变化信号中，有用信息更多，从而可以使得根据检测视频得到的待检测装置的运行状况更加准确。
附图说明
下面将对本发明实施例涉及的一些附图进行说明。
图1为本发明实施例的电子装置的硬件结构示意图；
图2为本发明实施例的信号处理方法的流程示意图；
图3为本发明实施例的信号处理方法用于检测待检测装置的运行状况时涉及的放大检测视频的过程示意图；
图4为本发明实施例的信号处理方法的另一流程示意图；
图5为本发明实施例的信号处理方法的又一流程示意图；
图6为本发明实施例的信号处理方法的再一流程示意图；
图7为本发明实施例的信号处理方法的再另一流程示意图；
图8为本发明实施例的信号处理方法的再又一流程示意图。
具体实施方式
下面结合本发明实施例中的附图对本发明实施例进行描述。
请参阅图1，图1是本发明实施例提供的电子装置100的硬件结构示意图。电子装置100 包括处理器101、存储器102、输入输出接口103，以及一个或多个程序，一个或多个程序被存储在存储器102中，并且被配置由处理器101执行，程序包括用以下任一实施例的信号处理方法的步骤的指令。电子装置100可以是服务器装置，也可以是终端装置。存储器102可以是高速RAM存储器，也可以是稳定的存储器(non-volatile memory)，例如磁盘存储器，存储器102可选的还可以是独立于前述处理器101的存储装置。输入输出接口103可选的可以包括USB接口、标准的有线接口、无线接口(如WI-FI接口)。
请参阅图2，图2是本发明实施例提供的一种信号处理方法的流程示意图，这种方法可包括但不限于如下步骤：
01、获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱对应的多个状态变化信号，状态变化信号为时域信号；
本发明实施例的滤波方法可用于但不限于检测待检测装置的振动情况。待检测装置的振动情况可以反应待检测装置的运行状况。检测视频可为由成像模组拍摄待检测装置得到的视频。可通过分析检测视频，从检测视频中提取出待检测装置的状态变化信息，根据状态变化信息得到振动信息，然后分析该振动信息得到待检测装置的运行状况。
在一个优选实施例中，检测视频为将成像模组拍摄待检测装置得到的视频采用欧拉算法对视频中的振动进行放大之后的振动放大视频。待检测装置包括但不限于各类机械设备、建筑结构等。得到检测视频之后，对检测视频对应的帧序列进行相位相关计算得到交叉互功率谱，交叉互功率谱中包括一个或多个相关峰。这样得到的交叉互功率谱为频域信号，将交叉互功率谱进行反傅里叶变换之后即可得到各相关峰对应的状态变化信号。每个状态变化信号可以反应检测视频中的某个位置的状态变化情况，对检测视频对应的帧序列进行相位相关计算得到交叉互功率谱可以理解为，从检测视频中提取出视频画面中的状态变化信息。状态变化信息包括振动信息和其他噪声信息，例如光照的变化也会导致视频画面中的状态变化，振动信息可以反应待检测装置的运行状况。本发明实施例中，通可通过本发明实施例的滤波方法，对交叉互功率谱中的各相关峰对应的状态变化信号进行处理，以使得根据检测视频得到的待检测装置的运行状况更加准确。
02、分析多个状态变化信号，以获取多个状态变化信号中的非周期信号；
可以理解，待检测装置的振动是一种周期往复运动，振动引起的状态变化也是周期性的。很多噪声信息虽然会导致检测视频中各像素点的状态变化，但是，噪声引起的状态变化常常不是周期性的，例如在室外场景下，由于室外光纤的明暗变化导致的状态变化。而且根据待检测装置的振动分析待检测装置的运行状况时，呈周期性的振动才能用于反应待检测装置的运行状况，因为非周期性的振动常常是由外界环境导致的，而不是由待检测装置自身引起的，那这部分非周期性的信号也不能用于分析待检测装置的运行状况。那么可通过获取状态变化信号中的非周期信号获取不是由自身振动引起的噪声信号。
03、去除多个状态变化信号中的非周期信号。
获得多个状态变化信号中的非周期信号之后，由于非周期信号常常是对分析待检测装置的运行状况作用不大或者没有作用甚至会有干扰的信号，那么可去除这部分非周期信号，以使得从检测视频中获取到的状态变化信号中，有用信息更多。
本发明实施例的信号处理方法，获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱对应的多个状态变化信号，状态变化信号为时域信号；分析多个状态变化信号，以获取多个状态变化信号中的非周期信号；去除多个状态变化信号中的非周期信号。这样从检测视频中获取到的状态变化信号中，有用信息更多，从而可以使得根据检测视频得到的待检测装置的运行状况更加准确。
检测待检测装置的运行状况时，需要从包含有检测装置运行时的画面的检测视频中提取振动信息，在检测视频中，振动信息可由状态变化信息体现。如图3所示，可先将检测视频中的帧序列由RGB颜色空间转换到YIQ颜色空间，分离视频帧的亮度信息和色度信息。RGB 和YIQ的转换关系为：
Y＝0.299*R+0.587*G+0.114*B；
I＝0.596*R–0.275*G–0.321*B；
Q＝0.212*R-0.523*G+0.311*B。
然后保持I、Q通道不变，对Y通道进行FFT操作，再利用欧拉运动放大算法对视频数据行放大处理。具体包括：先将FFT变换后的Y通道图像进行复数可操纵金子塔空域分解，得到由多个不同空间分辨率的子图像组成的金字塔形结构；对金字塔型结构中的多个子图像中每个子图像进行时域带通滤波处理，得到目标频带对应的变换信号。可以理解，在视频画面中，振动可由视频帧序列的亮度反映，那么可通过分析检测视频中的Y通道获取检测视频中的振动信息。然后放大时域带通滤波后得到的目标频带对应的变换信号并提取感兴趣的运动信息，对感兴趣的运动信息进行复数可操纵金字塔重建，得到放大后的Y通道图像；最后将重建的Y通道图像与原来的I、Q通道图像相加，再转化为RGB色彩空间，得到输出视频。
再对视频运动放大处理后的帧序列采用相位相关算法计算帧序列间的交叉互功率谱。相位相关算法采用如下的公式计算交叉互功率谱。
上式中，Fa为a帧图像的傅立叶变换，为b帧图像的傅里叶变换的共轭信号，为两个傅里叶变换的信号的相关积的模。R为本步骤的计算结果交叉互功率谱(包含频域噪音)。这样得到的交叉互功率谱为频域信号，将交叉互功率谱进行反傅立叶变换，逐相位比较，便可得到多个状态变化信号，可以反应检测视频中各处的状态变化信息。
最后按照本发明实施例的上述滤波方法对多个状态变化信号进行处理，去除非周期信号，提高了有用信号的信噪比，使得从包含有检测装置运行时的画面的检测视频中提取出的信息，有用信号的信噪比更高。
具体地，对目标振动视频的图像帧序列进行金字塔分解后，获得金字塔型结构。金字塔型结构包括多层子图像，且从上到下的图像分辨率依次降低，空间频率依次降低。对每层图像进行时域带通滤波处理，是为了得到目标频带，使得目标频带的子图像分辨率能够清晰表达图像的运动特征，同时不会因为分辨率过高造成计算量过大。因此，按照金字塔结构从下到上获取每个子图像在不同空间频率的频带，并且按照获取顺序与标准频带进行对比，当确定某一层子图像的频带与标准频带匹配成功时，就不需要再对该层子图像之上的一层或多层子图像进行空间频率的频带获取和匹配，提升了时域带通滤波的效率。另外，金字塔分解的最大层数确定：log2(min(xres,yres))，其中xres为图像的宽度像素值，yres为图像的高度像素值。
进一步地，对金字塔型结构中的多个子图像中每个子图像进行时域带通滤波处理，得到目标频带对应的变换信号，具体包括：
获取金字塔结构的总层数N，其中金字塔结构由下到上层数依次增大，对应的子图像分辨率依次升高；
对金字塔结构中的多个子图像进行分组，其中第一组包括2n+1层对应的子图像，第二组包括2(n+1)层对应的子图像，其中n为大于或等于0的整数，N≥max(2(n+1),2n+1)；
将第一组子图像采用第一处理器按照层数从小到大进行时域带通滤波处理，将第二组子图像采用第二处理器按照层数从小到大进行时域带通滤波处理，时域带通滤波处理为将子图像对应的不同空间频率的频带与标准频带进行匹配，第一处理器和第二处理器为独立运行的处理器；
当第一处理器或第二处理器确定频带与标准频带匹配成功时，确定频带为目标频带；
获取目标频带对应的子图像为第一目标子图像，获取第一目标子图像上面一层子图像为第二目标子图像；
确定第一目标子图像和第二目标子图像为目标频带对应的变换信号。
具体地，在对目标振动视频的多帧图像组成的帧序列进行空域金字塔分解后，可以获得金字塔结构的总层数N，且金字塔结构从下到上1～N层对应子图像的分辨率依次升高。然后对N层子图像进行分组，包括基数组和偶数组，基数组采用第一处理器按照1,3,5…2n+1的顺序进行时域带通滤波处理，偶数组采用第二处理器按照2,4,6,…2(n+1)的顺序进行时域带通滤波处理，两个处理器可以同时开始运行，也可以具有一定的处理时间间隔，另外，时域带通滤波处理是为了获取目标频带，获取方法可以是将子图像对应的不同空间频率的频带与标准频带进行匹配。当第一处理器或第二处理器确定频带与标准频带匹配成功时，确定所述频带为目标频带。
因为目标频带表示能反应图像振动信息的最低分辨率子图像的空间频率，那么确定目标频带对应的子图像为第一目标子图像，然后获取目标频带对应子图像上一层的子图像作为第二目标子图像，将第一目标子图像和第二目标子图像为目标频带对应的变换信号，进行后续放大变换。这样，可以更准确地对目标振动视频对应的帧图像进行放大，获得更准确的运动信息，同时也避免了对更高分辨率图像进行放大的需要增加的运算量。
在本申请实施例中，对金字塔结构对应的子图像进行分组，然后通过两个独立运行的处理器对不同组的子图像进行时域带通滤波处理，可以提升滤波处理效率，同时在获取到目标频带对应层的子图像作为第一目标子图像后，获取其上一层的子图像作为第二目标子图像，然后将第一目标子图像和第二目标子图像为目标频带对应的变换信号，这样可以获得更准确的放大运动信息，同时也避免了对更高分辨率图像进行放大的需要增加的运算量。
请参阅图4，基于上述实施例，在某些实施例中，步骤02包括但不限于以下步骤：
021、从每个状态变化信号中提取多个第一状态变化信号段，并获取各第一状态变化信号段的频率；
状态变化信号为时域信号，可从各状态变化信号中提取不同时段的第一状态变化信号段，然后通过数学方法获得各第一状态变化信号段的频率。从每个状态变化信号中提取多个第一状态变化信号段可理解为，从每个状态变化信号中提取出至少2个第一状态变化信号段，然后获取各个第一状态变化信号段的频率。多个可以为2个、3个、4个或其他任意大于4的个数。
022、将多个状态变化信号中，对应的多个第一状态变化信号段的频率不一致的状态变化信号作为非周期信号。
可将一个状态变化信号对应的多个第一状态变化信号段的频率进行比较，判断多个第一状态变化信号段的频率是否一致。如果一个状态变化信号对应的多个第一状态变化信号段的频率不一致，例如，从一个状态变化信号中提取出了4个第一状态变化信号段，这4个第一状态变化信号段中，有两个第一状态变化信号段的频率为a，有两个第一状态变化信号段的频率为b，则该状态变化信号段对应的第一状态变化信号段的频率不一致，说明该状态变化信号至少存在两个不同频率值的频段，则可将该状态变化信号作为非周期信号。
请参阅图5，在某些实施例中，步骤02包括但不限于以下步骤：
023、从各状态变化信号中提取预设长度的目标状态变化信号段,并获取目标状态变化信号段的目标频率；
判断各状态变化信号是否为周期信号时，可先提取预设长度的目标状态变化信号段，获取目标变化信号段的目标频率，然后将状态变化信号中的其他部分的频率与目标频率比对，如果状态变化信号中的其他部分的频率与目标频率不一致，则可认为该状态变化信号为非周期信号。预设长度可以由用户设定一个确定的值，也可以在信号处理过程中根据信号的长度自行适配，例如可将预设长度设定为状态变化信号长度的1/10。
024、根据各状态变化信号对应的目标频率设定对应的滑窗,并将各状态变化信号发送至对应的滑窗；
获得目标状态变换信号段的目标频率之后，根据目标频率设定滑窗的窗口大小，例如可将滑窗的窗口大小设定为与目标频率一致，这样只有频率与目标频率一致的信号才能通过滑窗，而与目标频率不一致的信号则不能通过滑窗。
025、将不能通过对应的滑窗的状态变化信号作为非周期信号。
如果状态变化信号不能通过对应的滑窗，则说明该状态变化信号中存在频率与目标频率不一致的信号段，也即该状态变化信号为非周期信号。
本实施例中，通过滑窗的方式判断状态变化信号中的其他部分的频率与目标频率是否一致，这样可以方便快捷地得到结论，计算量更小。
请参阅图6，在某些实施例中，步骤01包括：
011、获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱，交叉互功率谱为频域信号；
012、对交叉互功率谱进行反傅里叶变换得到多个状态变化信号。
这样实现将频域信号转换为时域信号，便于对识别出非周期信号，并去除非周期信号，提高了有用信号的信噪比，使得从包含有检测装置运行时的画面的检测视频中提取出的信息，有用信号的信噪比更高。
进一步地，请参阅图7，步骤012之前，步骤01还包括：
013、根据交叉互功率谱中的各相关峰在检测视频中对应的位置及相关峰的频段确定各相关峰对应的滤波策略；
步骤013在步骤011之后、步骤012之前执行。可分析各相关峰得到各相关峰在检测视频中对应的位置，相关峰在检测视频中对应的位置可以理解为该相关峰所反映的状态变化对应的检测视频中的位置。例如，检测视频中有一块位置的亮度随着时间呈周期性变化，对检测视频的对应的帧序列进行相位相关计算得到的一个或多个相关峰中，则包含至少一个相关峰可以体现该位置的亮度变化。
由于检测视频中，不同位置的振动对检测待检测装置的运行状况的作用大小不一样。例如，处于检测视频的画面中的边缘位置的振动相比较于中间位置的振动，对检测待检测装置的运行状况的作用更小，那么可将各相关峰在检测视频中对应的位置作为依据，确定各相关峰的滤波策略。每个相关峰对应的频率范围(即频段)不一样，根据相关峰对应的频段确定各相关峰的滤波策略，这样使得滤波策略与每个相关峰更加适配，从而有助于从各相关峰中提取出有用的状态变化信号。滤波策略例如可包括但不限于滤波器的带宽，滤波器的类型等。
014、按照各相关峰对应的滤波策略对交叉互功率谱中的各相关峰进行滤波处理；
获得各相关峰的滤波策略之后，按照对应的滤波策略对各相关峰进行滤波处理，以使得从检测视频中获得的信号的信噪比更高。较佳地，可以按照各相关峰对应的滤波策略对各相关峰通过插值滤波的方式进行滤波处理，这样滤波效果更好。
步骤012包括：
0121、对进行滤波处理之后的交叉互功率谱进行反傅里叶变换得到多个状态变化信号。
这样在进行反傅里叶变换得到状态变化信号之前，先对交叉互功率谱进行滤波处理，这样可以使得后续识别非周期信号的过程运算量更小。而且，通过根据各相关峰的位置及对应的频段选择合适的滤波策略进行滤波处理，在去除非周期信号之前，先进行一次去燥，减少各相关峰中的噪声信息，进一步地使得从检测视频中获得的状态变化信号的信噪比更高，从而使得根据检测视频得到的待检测装置的运行状况更加准确。
更进一步地，在根据交叉互功率谱中的各相关峰在检测视频中对应的位置及相关峰的频段确定各相关峰对应的滤波策略之前，还可以先利用主成分分解方法(PCA)对各相关峰进行降维处理，这样在步骤013中，对降维处理后的相关峰进行计算，得到各相关峰对应的滤波策略，这样可降低滤波过程的计算量，提高信号处理效率。主成分分解(PCA)用来将相关峰数据在振动检测的两个主方向上对数据维数约减。
再进一步地，在利用主成分分解方法对各相关峰进行降维处理之前对相关峰进行筛选，去除峰值不在预设范围内的相关峰，对峰值在预设峰值范围内的相关峰进行降维处理，过滤不在预设峰值范围内的相关峰，这样一方面可以提高从检测视频中获得的状态变化信号的信噪比，有用信号更多，另一方面由于去除了部分的噪声信号，从而可以减少后续的过滤的计算量，提升过滤效率。
请参阅图8，基于上述实施例，在某些实施例中，滤波策略包括滤波带宽，步骤013包括：
0131、根据交叉互功率谱中的各相关峰对应的频段确定各相关峰对应的初始滤波带宽；
可以理解，对检测视频对应的帧序列进行相位相关计算得到的相关峰中，包括检测视频对应的画面中位置相近的多个像素点的相似的多个状态变化信息，从而可以反应多个像素点所在的位置的状态变化。那么检测视频对应的画面中，不同位置的状态变化情况不同，对应的相关峰的频率分布也不同。如果用相同的滤波带宽进行滤波，那么为了避免有用信号被过滤，需要将滤波带宽设为较大值，但是这样有会导致噪声信号无法被过滤。而本申请的滤波方法中，根据各相关峰对应的频段确定个相关峰对应的初始滤波带宽，相关峰对应的频段与初始滤波带宽正相关，相关峰对应的频段越大，对应的初始滤波带宽越宽，可使得滤波带宽更好地与各相关峰适配，既能够有效地过滤噪声信号，也有助于避免有用信号被过滤，从而可以使得过滤效果更好。
例如，可以将B*(1-a)作为初始滤波带宽，B为各相关峰对应的频段，a可以根据实际需求确定，当需要得到更好的去燥效果时，可将a设为较大值，当需要获取到更多的有用信号时，可将a设定为较小值。当然，在其他实施例中，初始滤波带宽的设定方法不限于上述举例，在此不做限定。
0132、根据各相关峰在检测视频中对应的位置确定相关峰对应的带宽调整系数；
检测视频中不同位置的振动对检测待检测装置的运行状况的作用大小不一样。例如画面边缘的振动与画面中心的振动相比，画面边缘的作用对检测待检测装置的运行状况的作用更小。那么可根据各相关峰在检测视频中对应的位置确定相关峰对应的带宽调整系数。利用带宽调整系数来调整初始滤波带宽，使得滤波带宽更合理。
例如，可根据相关峰在检测视频中对应的位置的振动对检测待检测装置的运行状况的作用大小，将检测视频对应的画面分为多个区域，并为每个区域预先设定预设带宽调整系数，这样在步骤022中确定带宽调整系数时，可直接获取各相关峰在检测视频中对应的位置所在的区域对应的预设带宽调整系数，并将预设带宽调整系数作为相关峰对应的带宽调整系数。这样可以快速准确的获得带宽调整系数。
进一步地，用户可根据待检测装置的结构，对检测视频对应的画面中的各个区域设定带宽加权系数，然后将预设带宽调整系数与带宽加权系数的乘积作为对应的区域的相关的带宽调整系数，这样使得得到的各个相关峰的滤波带宽更加合理。例如，待检测装置中某些关键部位的振动相对其他部位更加能够反应待检测装置的运行状况，那么可将该关键部位在检测视频的画面中对应的区域的带宽加权系数设定为较大值，使得滤波带宽更大，从而可以提取到更多的有用信息。
0133、将各相关峰对应的初始滤波带宽和带宽调整系数的乘积作为各相关峰对应的滤波带宽。
这样得到的滤波带宽，综合考虑了相关峰在检测视频中对应的位置及相关峰对应的频段，使得滤波带宽与每个相关峰更加适配，从而有助于从交叉互功率谱中提取出有用的状态变化信号。
需要说明的是，滤波带宽的计算方法不限于上述方法，在其他实施例中也可以根据各相关峰对应的频段和各相关峰在检测视频中对应的位置选择合适的计算方法来确定滤波带宽。
本发明实施例还提供一种信号处理装置，包括：
信号获取模块，用于获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱对应的多个状态变化信号，状态变化信号为时域信号；
分析模块，用于分析多个状态变化信号，以获取多个状态变化信号中的非周期信号；
信号处理模块，用于去除多个状态变化信号中的非周期信号。
本发明实施例的信号处理装置，获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱对应的多个状态变化信号，状态变化信号为时域信号；分析多个状态变化信号，以获取多个状态变化信号中的非周期信号；去除多个状态变化信号中的非周期信号。这样从检测视频中获取到的状态变化信号中，有用信息更多，从而可以使得根据检测视频得到的待检测装置的运行状况更加准确。需要说明的是，上述信号处理方法各实施例的解释说明及技术效果也适用于本实施例的信号处理装置，为避免冗余，在此不再赘述。
在某些实施例中，分析模块包括：
第一信号提取单元，用于从每个状态变化信号中提取多个第一状态变化信号段，并获取各第一状态变化信号段的频率；
第一执行单元，用于将多个状态变化信号中，对应的多个第一状态变化信号段的频率不一致的状态变化信号作为非周期信号。
在某些实施例中，分析模块包括：
第二信号提取单元，用于从各状态变化信号中提取预设长度的目标状态变化信号段,并获取目标状态变化信号段的目标频率；
信号发送单元，用于根据各状态变化信号对应的目标频率设定对应的滑窗,并将各状态变化信号发送至对应的滑窗；
第二执行单元，用于将不能通过对应的滑窗的状态变化信号作为非周期信号。
在某些实施例中，信号获取模块包括：
信号获取单元，用于获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱，交叉互功率谱为频域信号；
信号转换单元，用于对交叉互功率谱进行反傅里叶变换得到多个状态变化信号。
在某些实施例中，信号获取模块还包括：
滤波策略确定单元，用于根据交叉互功率谱中的各相关峰在检测视频中对应的位置及相关峰的频段确定各相关峰对应的滤波策略；
滤波单元，用于按照各相关峰对应的滤波策略对交叉互功率谱中的各相关峰进行滤波处理；
信号转换单元还用于：
对进行滤波处理之后的交叉互功率谱进行反傅里叶变换得到多个状态变化信号。
在某些实施例中，滤波策略包括滤波带宽，滤波策略确定单元包括：
初始滤波带宽确定子单元，用于根据交叉互功率谱中的各相关峰对应的频段确定各相关峰对应的初始滤波带宽；
带宽调整系数确定子单元，用于根据各相关峰在检测视频中对应的位置确定相关峰对应的带宽调整系数；
滤波带宽确定子单元，用于将各相关峰对应的初始滤波带宽和带宽调整系数的乘积作为各相关峰对应的滤波带宽。
其中，上述信号处理装置中各个模块的功能实现及技术效果与上述信号处理方法实施例中各步骤相对应，其功能和实现过程在此处不再一一赘述。
本发明还提供一种计算机可读存储介质，计算机可读存储介质上存储有信号处理程序，其中信号处理程序被处理器执行时，实现上述任一实施例的信号处理方法的步骤。
其中，信号处理程序被执行时所实现的方法及对应的技术效果可参照本发明信号处理方法的各个实施例，此处不再赘述。
本发明实施例还提供一种检测设备，用于检测待检测装置的运行状况，检测设备包括上述实施例的信号处理装置或上述实施例的电子装置。
本发明实施例的检测设备，获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱对应的多个状态变化信号，状态变化信号为时域信号；分析多个状态变化信号，以获取多个状态变化信号中的非周期信号；去除多个状态变化信号中的非周期信号。这样从检测视频中获取到的状态变化信号中，有用信息更多，从而可以使得根据检测视频得到的待检测装置的运行状况更加准确。
检测待检测装置的运行状况时，需要从包含有检测装置运行时的画面的检测视频中提取振动信息，在检测视频中，振动信息可由状态变化信息体现。如图3所示，可先将检测视频中的帧序列由RGB颜色空间转换到YIQ颜色空间，分离视频帧的亮度信息和色度信息。RGB 和YIQ的转换关系为：
Y＝0.299*R+0.587*G+0.114*B；
I＝0.596*R–0.275*G–0.321*B；
Q＝0.212*R-0.523*G+0.311*B。
然后保持I、Q通道不变，对Y通道进行FFT操作，再利用欧拉运动放大算法对视频数据行放大处理。具体包括：先将FFT变换后的Y通道图像进行复数可操纵金子塔空域分解，将Y通道空域分解后的不同尺度的图像进行时域带通滤波，可以理解，在视频画面中，振动可由视频帧序列的亮度反映，那么可通过分析检测视频中的Y通道获取检测视频中的振动信息。然后放大时域带通滤波后感兴趣的运动信息，对感兴趣的运动信息进行复数可操纵金字塔重建，得到放大后的Y通道图像；最后将重建的Y通道图像与原来的I、Q通道图像相加，再转化为RGB色彩空间，得到输出视频。
再对视频运动放大处理后的帧序列采用相位相关算法计算帧序列间的交叉互功率谱。相位相关算法采用如下的公式计算交叉互功率谱。
上式中，Fa为a帧图像的傅立叶变换，为b帧图像的傅里叶变换的共轭信号，为两个傅里叶变换的信号的相关积的模。R为本步骤的计算结果交叉互功率谱(包含频域噪音)。这样得到的交叉互功率谱为频域信号，将交叉互功率谱进行反傅立叶变换，逐相位比较，便可得到多个状态变化信号，可以反应检测视频中各处的状态变化信息。
最后按照本发明实施例的上述滤波方法对多个状态变化信号进行处理，去除非周期信号，提高了有用信号的信噪比，使得从包含有检测装置运行时的画面的检测视频中提取出的信息，有用信号的信噪比更高。
在上述实施例中，可全部或部分地通过软件、硬件、固件、或其任意组合来实现。当使用软件实现时，可以全部或部分地以计算机程序产品的形式实现。所述计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行所述计算机程序指令时，全部或部分地产生按照本发明实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中，或者从一个计算机可读存储介质向另一个计算机可读存储介质传输，例如，所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线)或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输。所述计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质(例如软盘、硬盘、磁带)、光介质(例如光盘)、或者半导体介质(例如固态硬盘)等。在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。
在本发明所提供的几个实施例中，应该理解到，所揭露的装置，也可以通过其它的方式实现。例如以上所描述的装置实施例仅仅是示意性的，例如所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可结合或者可以集成到另一个系统，或一些特征可以忽略或不执行。另一点，所显示或讨论的相互之间的间接耦合或者直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者，也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例的方案的目的。
另外，在本发明各实施例中的各功能单元可集成在一个处理单元中，也可以是各单元单独物理存在，也可两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，或者也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质例如可包括：U盘、移动硬盘、只读存储器(ROM，Read-OnlyMemory)、随机存取存储器(RAM，Random Access Memory)、磁碟或光盘等各种可存储程序代码的介质。
Claims (9)
1.一种信号处理方法，其特征在于，包括：
获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱对应的多个状态变化信号，所述状态变化信号为时域信号，包括：获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱，所述交叉互功率谱为频域信号，所述交叉互功率谱包括多个相关峰，对所述交叉互功率谱包进行反傅里叶变换，得到所述多个相关峰对应的多个状态变化信号；
分析所述多个状态变化信号，以获取所述多个状态变化信号中的非周期信号；
去除所述多个状态变化信号中的非周期信号。
2.根据权利要求1所述的信号处理方法，其特征在于，所述分析所述多个状态变化信号，以获取所述多个状态变化信号中的非周期信号包括：
从每个状态变化信号中提取多个第一状态变化信号段，并获取各第一状态变化信号段的频率；
将所述多个状态变化信号中，对应的多个第一状态变化信号段的频率不一致的状态变化信号作为所述非周期信号。
3.根据权利要求1所述的信号处理方法，其特征在于，所述分析所述多个状态变化信号，以获取所述多个状态变化信号中的非周期信号包括：
从各状态变化信号中提取预设长度的目标状态变化信号段,并获取所述目标状态变化信号段的目标频率；
根据各状态变化信号对应的目标频率设定对应的滑窗,并将各状态变化信号发送至对应的滑窗；
将不能通过对应的滑窗的状态变化信号作为所述非周期信号。
4.根据权利要求3所述的信号处理方法，其特征在于，所述对所述交叉互功率谱进行反傅里叶变换得到多个状态变化信号之前，所述获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱对应的多个状态变化信号还包括：
根据所述交叉互功率谱中的各相关峰在所述检测视频中对应的位置及相关峰的频段确定各相关峰对应的滤波策略；
按照各相关峰对应的滤波策略对所述交叉互功率谱中的各相关峰进行滤波处理；
所述对所述交叉互功率谱进行反傅里叶变换得到多个状态变化信号包括：
对进行滤波处理之后的交叉互功率谱进行反傅里叶变换得到多个状态变化信号。
5.根据权利要求4所述的信号处理方法，其特征在于，所述滤波策略包括滤波带宽，所述根据所述交叉互功率谱中的各相关峰在所述检测视频中对应的位置及相关峰的频段确定各相关峰对应的滤波策略包括：
根据所述交叉互功率谱中的各相关峰对应的频段确定各相关峰对应的初始滤波带宽；
根据各相关峰在所述检测视频中对应的位置确定所述相关峰对应的带宽调整系数；
将各相关峰对应的初始滤波带宽和带宽调整系数的乘积作为各相关峰对应的滤波带宽。
6.一种信号处理装置，其特征在于，包括：
信号获取模块，用于获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱对应的多个状态变化信号，所述状态变化信号为时域信号，
所述信号获取模块还用于，获取对检测视频对应的帧序列进行相位相关计算得到的交叉互功率谱，所述交叉互功率谱为频域信号，所述交叉互功率谱包括多个相关峰，对所述交叉互功率谱包进行反傅里叶变换，得到所述多个相关峰对应的多个状态变化信号；
分析模块，用于分析所述多个状态变化信号，以获取所述多个状态变化信号中的非周期信号；
信号处理模块，用于去除所述多个状态变化信号中的非周期信号。
7.一种电子装置，其特征在于，包括处理器、存储器、以及存储在所述存储器上并可被所述处理器执行的信号处理程序，其中所述信号处理程序被所述处理器执行时，实现权利要求1至5中任一项所述的信号处理方法的步骤的指令。
8.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有信号处理程序，其中所述信号处理程序被处理器执行时，实现权利要求1至5中任一项所述的信号处理方法的步骤。
9.一种检测设备，用于检测待检测装置的运行状况，其特征在于，所述检测设备包括权利要求7所述的信号处理装置或权利要求7所述的电子装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2020/105591 WO2021036662A1 (zh) | 2019-04-26 | 2020-07-29 | 信号处理方法、装置及相关产品 |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2019103460300 | 2019-04-26 | ||
| CN201910346030 | 2019-04-26 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110659583A CN110659583A (zh) | 2020-01-07 |
| CN110659583B true CN110659583B (zh) | 2023-11-28 |
Family
ID=69036713
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910813881.1A Active CN110659583B (zh) | 2019-04-26 | 2019-08-30 | 信号处理方法、装置及相关产品 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN110659583B (zh) |
| WO (1) | WO2021036662A1 (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110659583B (zh) * | 2019-04-26 | 2023-11-28 | 深圳市豪视智能科技有限公司 | 信号处理方法、装置及相关产品 |
| CN110674697B (zh) * | 2019-04-26 | 2023-04-25 | 深圳市豪视智能科技有限公司 | 滤波方法、装置及相关产品 |
| CN113569696B (zh) * | 2021-07-22 | 2023-06-06 | 福建师范大学 | 基于视频的人体微小震颤信号提取的方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101858939A (zh) * | 2009-04-10 | 2010-10-13 | 华为技术有限公司 | 一种谐波信号检测方法及装置 |
| CN103714522A (zh) * | 2014-01-06 | 2014-04-09 | 海信集团有限公司 | 图像降噪方法和图像降噪系统 |
| WO2016146265A1 (en) * | 2015-03-17 | 2016-09-22 | Zynaptiq Gmbh | Methods for extending frequency transforms to resolve features in the spatio-temporal domain |
| CN106491117A (zh) * | 2016-12-06 | 2017-03-15 | 上海斐讯数据通信技术有限公司 | 一种基于ppg心率测量技术的信号处理方法及装置 |
| CN106846371A (zh) * | 2016-12-23 | 2017-06-13 | 深圳大学 | 一种基于视频图像的鸟类翅膀振动频率测量方法及装置 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3404302B2 (ja) * | 1998-11-30 | 2003-05-06 | 株式会社栗本鐵工所 | 土木建築構造物の振動測定方法 |
| CN102034022B (zh) * | 2010-12-03 | 2012-10-03 | 深圳市理邦精密仪器股份有限公司 | 一种基于倍频分析的信号处理方法及系统 |
| ES2525517B2 (es) * | 2013-05-23 | 2015-05-22 | Universidad De Alicante | Método y sistema para medir vibraciones |
| CN104614064B (zh) * | 2015-02-13 | 2018-01-12 | 福州大学 | 一种基于条纹靶的高速多维度振动测量装置及方法 |
| US10567655B2 (en) * | 2016-09-19 | 2020-02-18 | Triad National Security, Llc | System and method for automated extraction of high resolution structural dynamics from video |
| CN106989812B (zh) * | 2017-05-03 | 2020-11-03 | 湖南科技大学 | 基于摄影测量技术的大型风机叶片模态测试方法 |
| CN110068388A (zh) * | 2019-03-29 | 2019-07-30 | 南京航空航天大学 | 一种基于视觉和盲源分离的振动检测方法 |
| CN110659583B (zh) * | 2019-04-26 | 2023-11-28 | 深圳市豪视智能科技有限公司 | 信号处理方法、装置及相关产品 |
-
2019
- 2019-08-30 CN CN201910813881.1A patent/CN110659583B/zh active Active
-
2020
- 2020-07-29 WO PCT/CN2020/105591 patent/WO2021036662A1/zh active Application Filing
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101858939A (zh) * | 2009-04-10 | 2010-10-13 | 华为技术有限公司 | 一种谐波信号检测方法及装置 |
| CN103714522A (zh) * | 2014-01-06 | 2014-04-09 | 海信集团有限公司 | 图像降噪方法和图像降噪系统 |
| WO2016146265A1 (en) * | 2015-03-17 | 2016-09-22 | Zynaptiq Gmbh | Methods for extending frequency transforms to resolve features in the spatio-temporal domain |
| CN106491117A (zh) * | 2016-12-06 | 2017-03-15 | 上海斐讯数据通信技术有限公司 | 一种基于ppg心率测量技术的信号处理方法及装置 |
| CN106846371A (zh) * | 2016-12-23 | 2017-06-13 | 深圳大学 | 一种基于视频图像的鸟类翅膀振动频率测量方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| 数字信号处理技术在汽轮机组振动测量中的应用;柴庆芬等;《仪器仪表用户》;20120608(第03期);全文 * |
| 非平稳强噪条件下的视频相位放大技术研究;杨金钊等;《机械与电子》;20180624(第06期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110659583A (zh) | 2020-01-07 |
| WO2021036662A1 (zh) | 2021-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110659583B (zh) | 信号处理方法、装置及相关产品 | |
| CN110674697B (zh) | 滤波方法、装置及相关产品 | |
| CN110595745B (zh) | 针对设备的固定螺丝异常的检测方法及相关产品 | |
| JP4453648B2 (ja) | 画像処理装置および撮像装置 | |
| US9324153B2 (en) | Depth measurement apparatus, image pickup apparatus, depth measurement method, and depth measurement program | |
| US7672528B2 (en) | Method of processing an image to form an image pyramid | |
| JP5983373B2 (ja) | 画像処理装置、情報処理方法及びプログラム | |
| JP2012256202A (ja) | 画像処理装置、および画像処理方法、並びにプログラム | |
| US10043244B2 (en) | Image processing device, image processing method, image processing program, and recording medium | |
| KR20160102524A (ko) | 이미지의 역 톤 매핑을 위한 방법 | |
| US9319602B2 (en) | Method and apparatus for processing image | |
| KR101705605B1 (ko) | 코드화된 조명을 이용하는 촬상 장치 및 이미지 처리 장치와 그 방법 | |
| US20130329004A1 (en) | Method of and Apparatus for Image Enhancement | |
| US20100092086A1 (en) | Method and system for image deblurring | |
| US11272146B1 (en) | Content adaptive lens shading correction method and apparatus | |
| JP4936686B2 (ja) | 画像処理 | |
| Ponomarenko et al. | Adaptive visually lossless JPEG-based color image compression | |
| JP2021043874A (ja) | 画像処理装置、画像処理方法、およびプログラム | |
| JP5765893B2 (ja) | 画像処理装置、撮像装置および画像処理プログラム | |
| US7269295B2 (en) | Digital image processing methods, digital image devices, and articles of manufacture | |
| US20100227428A1 (en) | Method and apparatus for manufacturing color image pickup device including a random array pattern and a regular array pattern | |
| JP3743389B2 (ja) | 画像圧縮装置および画像圧縮プログラム | |
| US20130016905A1 (en) | Method and apparatus for correcting color distortion | |
| JP2014146872A5 (zh) | ||
| JP5088361B2 (ja) | 画像処理装置および撮像装置 |
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 |