CN116342438A - 一种空间频谱自适应滤波方法和装置 - Google Patents
一种空间频谱自适应滤波方法和装置 Download PDFInfo
- Publication number
- CN116342438A CN116342438A CN202310627794.3A CN202310627794A CN116342438A CN 116342438 A CN116342438 A CN 116342438A CN 202310627794 A CN202310627794 A CN 202310627794A CN 116342438 A CN116342438 A CN 116342438A
- Authority
- CN
- China
- Prior art keywords
- target
- frequency
- spectrogram
- phase range
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 93
- 238000001914 filtration Methods 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 47
- 239000011159 matrix material Substances 0.000 claims abstract description 71
- 238000012545 processing Methods 0.000 claims abstract description 49
- 230000003044 adaptive effect Effects 0.000 claims abstract description 27
- 238000013519 translation Methods 0.000 claims abstract description 13
- 238000010276 construction Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 230000003595 spectral effect Effects 0.000 claims 1
- 230000008685 targeting Effects 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 9
- 238000004458 analytical method Methods 0.000 abstract description 7
- 238000005516 engineering process Methods 0.000 abstract description 6
- 238000004891 communication Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000010339 dilation Effects 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004556 laser interferometry Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010587 phase diagram Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration using non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
- G06T5/30—Erosion or dilatation, e.g. thinning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20004—Adaptive image processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20036—Morphological image processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
-
- 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
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明提供了一种空间频谱自适应滤波方法和装置，涉及激光干涉的技术领域，包括：获取待处理激光干涉图像，并对激光干涉图像进行傅里叶变换，得到频谱图；基于频谱图，构建目标矩阵，其中，目标矩阵用于表征频谱图的一倍频的频谱信号中目标像素，目标像素的幅值与频谱图的中心频谱的幅值之间的比值大于预设自适应阈值；对目标矩阵进行图像膨胀处理，确定出目标相位范围；对目标相位范围进行目标处理，得到自适应滤波后的相位，其中，目标处理包括：高斯低通滤波处理，平移处理和傅里叶逆变换处理，解决了现有的激光干涉条纹傅里叶分析技术在进行空间频谱空间滤波时的误差较大的技术问题。
Description
技术领域
本发明涉及激光干涉的技术领域，尤其是涉及一种空间频谱自适应滤波方法和装置。
背景技术
现有激光干涉多使用的条纹傅里叶分析技术，在频谱空间滤波时一般使用高斯滤波，而选择的频谱信号范围一般为一圆形，使用圆形直观上也相对较容易理解，这是因为在空间频谱图中，距离代表所相差的频率，圆形则代表所选择范围使用低通滤波，在物理上含义较为清楚。但是，使用圆形包裹空间频谱下的目标相位时，一方面，需避开零频噪声信息与二倍频干扰信息，另一方面，目标相位表现在频谱上的形状一般不规则，难以使用。由此，使用圆形包括目标相位进行滤波后多会出现信号的损失，引入较大噪声。
针对上述问题，还未提出有效的解决方案。
发明内容
有鉴于此，本发明的目的在于提供一种空间频谱自适应滤波方法和装置，以缓解了现有的激光干涉条纹傅里叶分析技术在进行频谱空间滤波时的误差较大的技术问题。
第一方面，本发明实施例提供了一种空间频谱自适应滤波方法，包括：获取待处理激光干涉图像，并对所述激光干涉图像进行傅里叶变换，得到频谱图；基于所述频谱图，构建目标矩阵，其中，所述目标矩阵用于表征所述频谱图的一倍频的频谱信号中目标像素，所述目标像素的幅值与所述频谱图的中心频谱的幅值之间的比值大于预设自适应阈值；对所述目标矩阵进行图像膨胀处理，确定出目标相位范围；对所述目标相位范围进行目标处理，得到自适应滤波后的相位，其中，所述目标处理包括：高斯低通滤波处理，平移处理和傅里叶逆变换处理。
进一步地，基于所述频谱图，构建目标矩阵，包括：确定出所述频谱图中的目标区域，其中，所述目标区域为零频信号和二倍频信号对应的区域；基于所述目标区域，滤除所述频谱图中的所述零频信号和所述二倍频信号，得到所述一倍频的频谱信号对应的区域；基于所述一倍频的频谱信号对应的区域，构建所述目标矩阵。
进一步地，基于所述一倍频的频谱信号对应的区域，构建所述目标矩阵，包括：基于所述一倍频的频谱信号对应的区域中各个像素的位置坐标，构建位置矩阵，其中，所述位置矩阵中的各个元素的元素值为0；确定出所述一倍频的频谱信号对应的区域中各个像素的幅值；基于所述一倍频的频谱信号对应的区域中各个像素的幅值和所述频谱图的中心频谱的幅值，确定出所述目标像素；将所述位置矩阵中与所述目标像素对应的元素的元素值修改为1，得到所述目标矩阵。
进一步地，对所述目标相位范围进行目标处理，得到自适应滤波后的相位，包括：确定出所述目标相位范围中幅值的最大值，利用所述最大值对所述目标相位范围进行高斯低通滤波，得到滤波后的目标相位范围；将所述滤波后的目标相位范围平移至所述频谱图的中心，得到平移后的目标相位范围；对所述平移后的目标相位范围进行傅里叶逆变换，得到所述自适应滤波后的相位。
第二方面，本发明实施例还提供了一种空间频谱自适应滤波装置，包括： 获取单元，用于获取待处理激光干涉图像，并对所述激光干涉图像进行傅里叶变换，得到频谱图；构建单元，用于基于所述频谱图，构建目标矩阵，其中，所述目标矩阵用于表征所述频谱图的一倍频的频谱信号中目标像素，所述目标像素的幅值与所述频谱图的中心频谱的幅值之间的比值大于预设自适应阈值；膨胀单元，用于对所述目标矩阵进行图像膨胀处理，确定出目标相位范围；处理单元，用于对所述目标相位范围进行目标处理，得到自适应滤波后的相位，其中，所述目标处理包括：高斯低通滤波处理，平移处理和傅里叶逆变换处理。
进一步地，所述构建单元，用于：确定出所述频谱图中的目标区域，其中，所述目标区域为零频信号和二倍频信号对应的区域；基于所述目标区域，滤除所述频谱图中的所述零频信号和所述二倍频信号，得到所述一倍频的频谱信号对应的区域；基于所述一倍频的频谱信号对应的区域，构建所述目标矩阵。
进一步地，所述构建单元，用于：基于所述一倍频的频谱信号对应的区域中各个像素的位置坐标，构建位置矩阵，其中，所述位置矩阵中的各个元素的元素值为0；确定出所述一倍频的频谱信号对应的区域中各个像素的幅值；基于所述一倍频的频谱信号对应的区域中各个像素的幅值和所述频谱图的中心频谱的幅值，确定出所述目标像素；将所述位置矩阵中与所述目标像素对应的元素的元素值修改为1，得到所述目标矩阵。
进一步地，所述处理单元，用于：确定出所述目标相位范围中幅值的最大值，利用所述最大值对所述目标相位范围进行高斯低通滤波，得到滤波后的目标相位范围；将所述滤波后的目标相位范围平移至所述频谱图的中心，得到平移后的目标相位范围；对所述平移后的目标相位范围进行傅里叶逆变换，得到所述自适应滤波后的相位。
第三方面，本发明实施例还提供了一种电子设备，包括存储器以及处理器，所述存储器用于存储支持处理器执行上述第一方面中所述方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。
第四方面，本发明实施例还提供了一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序。
在本发明实施例中，通过获取待处理激光干涉图像，并对所述激光干涉图像进行傅里叶变换，得到频谱图；基于所述频谱图，构建目标矩阵，其中，所述目标矩阵用于表征所述频谱图的一倍频的频谱信号中目标像素，所述目标像素的幅值与所述频谱图的中心频谱的幅值之间的比值大于预设自适应阈值；对所述目标矩阵进行图像膨胀处理，确定出目标相位范围；对所述目标相位范围进行目标处理，得到自适应滤波后的相位，其中，所述目标处理包括：高斯低通滤波处理，平移处理和傅里叶逆变换处理，达到了空间频谱滤波的降噪和信号保真的目的，进而解决了现有的激光干涉条纹傅里叶分析技术在进行频谱空间滤波时的误差较大的技术问题，从而实现了降低间频谱空间滤波时的误差的技术效果。
本发明的其他特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种空间频谱自适应滤波方法的流程图；
图2为本发明实施例提供的结构元素B的示意图；
图3为本发明实施例提供的一种空间频谱自适应滤波装置的示意图；
图4为本发明实施例提供的一种电子设备的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例一：
根据本发明实施例，提供了一种空间频谱自适应滤波方法的实施例，需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
图1是根据本发明实施例的一种空间频谱自适应滤波方法的流程图，如图1所示，该方法包括如下步骤：
步骤S102，获取待处理激光干涉图像，并对所述激光干涉图像进行傅里叶变换，得到频谱图；
步骤S104，基于所述频谱图，构建目标矩阵，其中，所述目标矩阵用于表征所述频谱图的一倍频的频谱信号中目标像素，所述目标像素的幅值与所述频谱图的中心频谱的幅值之间的比值大于预设自适应阈值；
步骤S106，对所述目标矩阵进行图像膨胀处理，确定出目标相位范围；
步骤S108，对所述目标相位范围进行目标处理，得到自适应滤波后的相位，其中，所述目标处理包括：高斯低通滤波处理，平移处理和傅里叶逆变换处理。
在本发明实施例中，通过获取待处理激光干涉图像，并对所述激光干涉图像进行傅里叶变换，得到频谱图；基于所述频谱图，构建目标矩阵，其中，所述目标矩阵用于表征所述频谱图的一倍频的频谱信号中目标像素，所述目标像素的幅值与所述频谱图的中心频谱的幅值之间的比值大于预设自适应阈值；对所述目标矩阵进行图像膨胀处理，确定出目标相位范围；对所述目标相位范围进行目标处理，得到自适应滤波后的相位，其中，所述目标处理包括：高斯低通滤波处理，平移处理和傅里叶逆变换处理，达到了空间频谱滤波的降噪和信号保真的目的，进而解决了现有的激光干涉条纹傅里叶分析技术在进行频谱空间滤波时的误差较大的技术问题，从而实现了降低间频谱空间滤波时的误差的技术效果。
下面对步骤S102进行详细说明。
待处理激光干涉图像一般可以写成如下形式：
其中
表示干涉图像的灰度强度，/>
表示参考波前与测试波前之间的相位差，/>
表示平均强度，/>
是在该位置处的调制强度。
条纹傅里叶分析技术利用快速傅里叶在频谱分离干涉图的背景项和相位项。利用复数形式，干涉图可以写成如下形式：
其中
表示线性载波条纹的函数，但一般并未规定载波条纹的方向。要获取
的值，一般的处理方法是先对式(2)进行傅里叶变化，在空间频域种提取所要求得的相位值。式(2)的傅里叶变换结果写作：
其中
是频域下的坐标，/>
是德尔塔函数，/>
和/>
是一对对应式(2)中等式右边第二项和第三项的复数共轭函数，线性载波函数
在频谱中将/>
和/>
两相偏移到了相反的两边，这一过程可以通过下述过程解释：
并且有
式(5)的结果实际上是对目标频谱区域叠加条纹载波后的结果。类似的，还有-1级条纹：
和
现有的空间频谱滤波方法在获取干涉图样后，对图像进行傅里叶变换，之后一般的处理方法是：以频谱图中1级一瓣的最大点为中心进行高斯低频滤波，选取合适的半径选取圆形信息进行平移，平移至图像中心，再进行傅里叶逆变换，得到目标相位。滤波过程如下式所示：
但在实际应用中，圆形区域的选取受限条件较多，难以应对复杂的相位环境，同时，该方法存在理论的最大解析梯度：
其中
是载波频率，x表示图像中的载波条纹方向。
在本发明实施例中，步骤S104包括如下步骤：
确定出所述频谱图中的目标区域，其中，所述目标区域为零频信号和二倍频信号对应的区域；
基于所述目标区域，滤除所述频谱图中的所述零频信号和所述二倍频信号，得到所述一倍频的频谱信号对应的区域；
基于所述一倍频的频谱信号对应的区域，构建所述目标矩阵。
具体的，基于所述一倍频的频谱信号对应的区域，构建所述目标矩阵，包括如下步骤：
基于所述一倍频的频谱信号对应的区域中各个像素的位置坐标，构建位置矩阵，其中，所述位置矩阵中的各个元素的元素值为0；
确定出所述一倍频的频谱信号对应的区域中各个像素的幅值；
基于所述一倍频的频谱信号对应的区域中各个像素的幅值和所述频谱图的中心频谱的幅值，确定出所述目标像素；
将所述位置矩阵中与所述目标像素对应的元素的元素值修改为1，得到所述目标矩阵。
在本发明实施例中，在得到频谱图之后，首先需要确定干涉图频谱中零频信号及二倍频信号的范围，并滤去相关信号，防止对所需相位信号造成噪声重叠，引起噪声。
接下来，确定目标相位范围，具体的，对于一倍频，即载波条纹的频率范围，以其坐标位置，构建一相同大小的零矩阵，称为位置矩阵，然后，分析频谱信号中的相位信息，对于信号大于中心频谱的幅值m%的信号，认为是目标相位，在位置矩阵中将此位置的数值置1。m的选取根据实际情况灵活试验调整。最终，所得到的矩阵中，数值为1的区域即为目标相位的范围。
下面对步骤S106进行详细说明。
对位置矩阵，使用半径为n的圆为模板进行图像膨胀处理，所得范围即为要进行平移的目标相位范围。膨胀操作的具体过程跟卷积操作类似，对于位置矩阵和圆模板结构元素B，结构元素B在位置矩阵上面移动，其中结构元素B定义其中心为锚点，计算B覆盖下位置矩阵的最大像素值用来替换锚点的像素，此处的B选取类似图2所示，对于方模板，圆内部分占大多数的正方形标识为黑色，即其矩阵数值为1，其他的为空白区域，矩阵数值为0。膨胀就是求局部最大值的操作，核B与位置矩阵卷积，即计算核B覆盖的区域的像素点的最大值，并把这个最大值赋值给参考点指定的像素。这样就会使图像中的高亮区域逐渐增长，在本发明实施例中，使感兴趣的目标相位范围增大。
在本发明实施例中，步骤S108包括如下步骤：
确定出所述目标相位范围中幅值的最大值，利用所述最大值对所述目标相位范围进行高斯低通滤波，得到滤波后的目标相位范围；
将所述滤波后的目标相位范围平移至所述频谱图的中心，得到平移后的目标相位范围；
对所述平移后的目标相位范围进行傅里叶逆变换，得到所述自适应滤波后的相位。
在本发明实施例中，将目标相位范围中幅值的最大值为中心对该范围进行高斯低通滤波后，再将该区域平移至图像频谱中心，进行傅里叶逆变换，即可得到经过自适应滤波后的相位。
本发明实施例中对空间频谱进行低通滤波及信号平移处理，能够避免目标信号被过度截断的问题，有效将测量相位保真传递，得到良好的相位复原结果。这是因为原有圆形信号选取方式没有考虑到实际信号的特点，容易将部分有用信号数据截断，或处理不当，选取信号与零频或二倍频信号产生重叠。
对干涉图样进行傅里叶变换获取的空间频谱，选择圆形滤波区域进行叶频谱是一种常见的方法。使用本发明的自适应滤波方法对包裹主要信息的椭圆中使用本方法确定滤波及平移数据的范围。在经过选定的滤波器并在频域中移动后，进行反傅立叶变换，得到干扰场的相位值。使用原本圆形滤波器重建出的相位图中出现了多处相位间断的现象。而使用本发明实施例中提出的自适应滤波后的相位，可以消除不连续的相位了。此外，用上一帧的相位减去本帧的相位，就可以得到相位差。比较圆形滤波器和自适应滤波器的相位差，能够明显发现自适应滤波器后的相位差结果背景噪声较小。此外，自适应滤波的相位差结果呈现出较大的变化充分说明了使用自适应滤波对干涉条纹处理产生的降噪及信号保真效果。
实施例二：
本发明实施例还提供了一种空间频谱自适应滤波装置，该空间频谱自适应滤波装置用于执行本发明实施例上述内容所提供的空间频谱自适应滤波方法，以下是本发明实施例提供的 装置的具体介绍。
如图3所示，图3为上述空间频谱自适应滤波装置的示意图，该空间频谱自适应滤波装置包括：
获取单元10，用于获取待处理激光干涉图像，并对所述激光干涉图像进行傅里叶变换，得到频谱图；
构建单元20，用于基于所述频谱图，构建目标矩阵，其中，所述目标矩阵用于表征所述频谱图的一倍频的频谱信号中目标像素，所述目标像素的幅值与所述频谱图的中心频谱的幅值之间的比值大于预设自适应阈值；
膨胀单元30，用于对所述目标矩阵进行图像膨胀处理，确定出目标相位范围；
处理单元40，用于对所述目标相位范围进行目标处理，得到自适应滤波后的相位，其中，所述目标处理包括：高斯低通滤波处理，平移处理和傅里叶逆变换处理。
在本发明实施例中，通过获取待处理激光干涉图像，并对所述激光干涉图像进行傅里叶变换，得到频谱图；基于所述频谱图，构建目标矩阵，其中，所述目标矩阵用于表征所述频谱图的一倍频的频谱信号中目标像素，所述目标像素的幅值与所述频谱图的中心频谱的幅值之间的比值大于预设自适应阈值；对所述目标矩阵进行图像膨胀处理，确定出目标相位范围；对所述目标相位范围进行目标处理，得到自适应滤波后的相位，其中，所述目标处理包括：高斯低通滤波处理，平移处理和傅里叶逆变换处理，达到了空间频谱滤波的降噪和信号保真的目的，进而解决了现有的激光干涉条纹傅里叶分析技术在进行频谱空间滤波时的误差较大的技术问题，从而实现了降低间频谱空间滤波时的误差的技术效果。
实施例三：
本发明实施例还提供了一种电子设备，包括存储器以及处理器，所述存储器用于存储支持处理器执行上述实施例一中所述方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。
参见图4，本发明实施例还提供一种电子设备100，包括：处理器50，存储器51，总线52和通信接口53，所述处理器50、通信接口53和存储器51通过总线52连接；处理器50用于执行存储器51中存储的可执行模块，例如计算机程序。
其中，存储器51可能包含高速随机存取存储器（RAM，Random Access Memory），也可能还包括非不稳定的存储器（non-volatilememory），例如至少一个磁盘存储器。通过至少一个通信接口53（可以是有线或者无线）实现该系统网元与至少一个其他网元之间的通信连接，可以使用互联网，广域网，本地网，城域网等。
总线52可以是ISA总线、PCI总线或EISA总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示，图4中仅用一个双向箭头表示，但并不表示仅有一根总线或一种类型的总线。
其中，存储器51用于存储程序，所述处理器50在接收到执行指令后，执行所述程序，前述本发明实施例任一实施例揭示的流过程定义的装置所执行的方法可以应用于处理器50中，或者由处理器50实现。
处理器50可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器50中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器50可以是通用处理器，包括中央处理器(CentralProcessingUnit，简称CPU)、网络处理器(NetworkProcessor，简称NP)等；还可以是数字信号处理器(Digital SignalProcessing，简称DSP)、专用集成电路(Application Specific Integrated Circuit，简称ASIC)、现成可编程门阵列(Field-Programmable Gate Array，简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器51，处理器50读取存储器51中的信息，结合其硬件完成上述方法的步骤。
实施例四：
本发明实施例还提供了一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，计算机程序被处理器运行时执行上述实施例一中所述方法的步骤。
另外，在本发明实施例的描述中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。
Claims (10)
1.一种空间频谱自适应滤波方法，其特征在于，包括：
获取待处理激光干涉图像，并对所述激光干涉图像进行傅里叶变换，得到频谱图；
基于所述频谱图，构建目标矩阵，其中，所述目标矩阵用于表征所述频谱图的一倍频的频谱信号中目标像素，所述目标像素的幅值与所述频谱图的中心频谱的幅值之间的比值大于预设自适应阈值；
对所述目标矩阵进行图像膨胀处理，确定出目标相位范围；
对所述目标相位范围进行目标处理，得到自适应滤波后的相位，其中，所述目标处理包括：高斯低通滤波处理，平移处理和傅里叶逆变换处理。
2.根据权利要求1所述的方法，其特征在于，基于所述频谱图，构建目标矩阵，包括：
确定出所述频谱图中的目标区域，其中，所述目标区域为零频信号和二倍频信号对应的区域；
基于所述目标区域，滤除所述频谱图中的所述零频信号和所述二倍频信号，得到所述一倍频的频谱信号对应的区域；
基于所述一倍频的频谱信号对应的区域，构建所述目标矩阵。
3.根据权利要求2所述的方法，其特征在于，基于所述一倍频的频谱信号对应的区域，构建所述目标矩阵，包括：
基于所述一倍频的频谱信号对应的区域中各个像素的位置坐标，构建位置矩阵，其中，所述位置矩阵中的各个元素的元素值为0；
确定出所述一倍频的频谱信号对应的区域中各个像素的幅值；
基于所述一倍频的频谱信号对应的区域中各个像素的幅值和所述频谱图的中心频谱的幅值，确定出所述目标像素；
将所述位置矩阵中与所述目标像素对应的元素的元素值修改为1，得到所述目标矩阵。
4.根据权利要求1所述的方法，其特征在于，对所述目标相位范围进行目标处理，得到自适应滤波后的相位，包括：
确定出所述目标相位范围中幅值的最大值，利用所述最大值对所述目标相位范围进行高斯低通滤波，得到滤波后的目标相位范围；
将所述滤波后的目标相位范围平移至所述频谱图的中心，得到平移后的目标相位范围；
对所述平移后的目标相位范围进行傅里叶逆变换，得到所述自适应滤波后的相位。
5.一种空间频谱自适应滤波装置，其特征在于，包括：
获取单元，用于获取待处理激光干涉图像，并对所述激光干涉图像进行傅里叶变换，得到频谱图；
构建单元，用于基于所述频谱图，构建目标矩阵，其中，所述目标矩阵用于表征所述频谱图的一倍频的频谱信号中目标像素，所述目标像素的幅值与所述频谱图的中心频谱的幅值之间的比值大于预设自适应阈值；
膨胀单元，用于对所述目标矩阵进行图像膨胀处理，确定出目标相位范围；
处理单元，用于对所述目标相位范围进行目标处理，得到自适应滤波后的相位，其中，所述目标处理包括：高斯低通滤波处理，平移处理和傅里叶逆变换处理。
6.根据权利要求5所述的装置，其特征在于，所述构建单元，用于：
确定出所述频谱图中的目标区域，其中，所述目标区域为零频信号和二倍频信号对应的区域；
基于所述目标区域，滤除所述频谱图中的所述零频信号和所述二倍频信号，得到所述一倍频的频谱信号对应的区域；
基于所述一倍频的频谱信号对应的区域，构建所述目标矩阵。
7.根据权利要求6所述的装置，其特征在于，所述构建单元，用于：
基于所述一倍频的频谱信号对应的区域中各个像素的位置坐标，构建位置矩阵，其中，所述位置矩阵中的各个元素的元素值为0；
确定出所述一倍频的频谱信号对应的区域中各个像素的幅值；
基于所述一倍频的频谱信号对应的区域中各个像素的幅值和所述频谱图的中心频谱的幅值，确定出所述目标像素；
将所述位置矩阵中与所述目标像素对应的元素的元素值修改为1，得到所述目标矩阵。
8.根据权利要求5所述的装置，其特征在于，所述处理单元，用于：
确定出所述目标相位范围中幅值的最大值，利用所述最大值对所述目标相位范围进行高斯低通滤波，得到滤波后的目标相位范围；
将所述滤波后的目标相位范围平移至所述频谱图的中心，得到平移后的目标相位范围；
对所述平移后的目标相位范围进行傅里叶逆变换，得到所述自适应滤波后的相位。
9.一种电子设备，其特征在于，包括存储器以及处理器，所述存储器用于存储支持处理器执行权利要求1至4任一项所述方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。
10.一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，其特征在于，计算机程序被处理器运行时执行上述权利要求1至4任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310627794.3A CN116342438B (zh) | 2023-05-31 | 2023-05-31 | 一种空间频谱自适应滤波方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310627794.3A CN116342438B (zh) | 2023-05-31 | 2023-05-31 | 一种空间频谱自适应滤波方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116342438A true CN116342438A (zh) | 2023-06-27 |
| CN116342438B CN116342438B (zh) | 2023-08-01 |
Family
ID=86879158
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310627794.3A Active CN116342438B (zh) | 2023-05-31 | 2023-05-31 | 一种空间频谱自适应滤波方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116342438B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101776489A (zh) * | 2010-02-05 | 2010-07-14 | 中国科学院上海光学精密机械研究所 | 飞秒激光等离子体通道干涉图相位和电子密度提取方法 |
| US20100189372A1 (en) * | 2009-01-23 | 2010-07-29 | National Taipei University Of Technology | Method for acquiring phase information and system for measuring three dimensional surface profiles |
| US20100288910A1 (en) * | 2009-05-14 | 2010-11-18 | Raytheon Company | Adaptive spatial-spectral processing (assp) |
| CN103812598A (zh) * | 2014-02-28 | 2014-05-21 | 北京邮电大学 | 一种基于空间频谱滤波的模式转换中频谱模场半径匹配方法 |
| CN104103048A (zh) * | 2014-08-01 | 2014-10-15 | 中国科学院遥感与数字地球研究所 | 一种合成孔径雷达干涉图频域滤波的方法 |
| US20190003959A1 (en) * | 2017-06-30 | 2019-01-03 | Guangdong University Of Technology | Blind separation based high accuracy perspective detection method for multilayer complex structure material |
| CN111383186A (zh) * | 2018-12-29 | 2020-07-07 | Tcl集团股份有限公司 | 一种图像处理方法、装置及终端设备 |
-
2023
- 2023-05-31 CN CN202310627794.3A patent/CN116342438B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100189372A1 (en) * | 2009-01-23 | 2010-07-29 | National Taipei University Of Technology | Method for acquiring phase information and system for measuring three dimensional surface profiles |
| US20100288910A1 (en) * | 2009-05-14 | 2010-11-18 | Raytheon Company | Adaptive spatial-spectral processing (assp) |
| CN101776489A (zh) * | 2010-02-05 | 2010-07-14 | 中国科学院上海光学精密机械研究所 | 飞秒激光等离子体通道干涉图相位和电子密度提取方法 |
| CN103812598A (zh) * | 2014-02-28 | 2014-05-21 | 北京邮电大学 | 一种基于空间频谱滤波的模式转换中频谱模场半径匹配方法 |
| CN104103048A (zh) * | 2014-08-01 | 2014-10-15 | 中国科学院遥感与数字地球研究所 | 一种合成孔径雷达干涉图频域滤波的方法 |
| US20190003959A1 (en) * | 2017-06-30 | 2019-01-03 | Guangdong University Of Technology | Blind separation based high accuracy perspective detection method for multilayer complex structure material |
| CN111383186A (zh) * | 2018-12-29 | 2020-07-07 | Tcl集团股份有限公司 | 一种图像处理方法、装置及终端设备 |
Non-Patent Citations (3)
| Title |
|---|
| 刘莹雪;王永红;高新亚;孙建飞;涂思琪;: "基于散斑干涉的相位测量技术研究", 工具技术, no. 09 * |
| 安海华;曹益平;李红梅;吴海涛;: "一种基于混合滤波窗口的复合光栅相位测量轮廓术", 中国激光, no. 06 * |
| 陈杰;阳雷;罗建东;: "面向小直径薄壁管表面质量检测的散斑干涉图像滤波方法研究", 科技视界, no. 15 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116342438B (zh) | 2023-08-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7483559B2 (en) | Method and apparatus for deblurring mask images | |
| EP2924972A1 (en) | Image processing apparatus and imaging system | |
| JP2016077904A (ja) | 撮像方法、画像処理装置、コンピュータ可読媒体、方法、装置およびシステム | |
| US20130202181A1 (en) | Integration of filters and phase unwrapping algorithms for removing noise in image reconstruction | |
| CN112034463A (zh) | 一种sar数据的高精度沉降监测方法和装置 | |
| CN114993452A (zh) | 基于宽带相位运动放大的结构微小振动测量方法与系统 | |
| Zhao et al. | An iterative Goldstein SAR interferogram filter | |
| CN112014349A (zh) | 太赫兹时域光谱成像复原方法、装置、存储介质及终端 | |
| US20160131767A1 (en) | Nonlinear processing for off-axis frequency reduction in demodulation of two dimensional fringe patterns | |
| CN116342438B (zh) | 一种空间频谱自适应滤波方法和装置 | |
| CN106651814B (zh) | 一种图像处理方法和装置 | |
| CN108445458B (zh) | 一种合成孔径雷达轨道误差消除方法和装置 | |
| CN114202476A (zh) | 一种红外图像增强方法、装置、设备及计算机可读介质 | |
| CN111723753B (zh) | 卫星遥感图像的去条带方法、装置和电子设备 | |
| CN110617775B (zh) | 一种基于视觉感知的桥梁挠度的监测方法、装置及系统 | |
| CN103607249A (zh) | 一种ds/fh混合扩频信号参数测试方法与系统 | |
| CN116336934B (zh) | 一种提高激光干涉测量精度的方法和装置 | |
| CN102967774B (zh) | 低雷达散射截面金属支架总散射场计算方法 | |
| CN115494303A (zh) | 一种emi接收机信号转化方法、装置及存储介质 | |
| CN114255179A (zh) | 地震图像噪声压制方法、装置及电子设备 | |
| Magnard et al. | GAMMA Technical Report: Offset estimation programs update | |
| CN110737024B (zh) | 地震波信号的处理方法、装置、设备及存储介质 | |
| CN117572464B (zh) | Gps测试方法、装置、电子设备及存储介质 | |
| EP3198554B1 (en) | Image processing | |
| CN116481416B (zh) | 基于北斗导航的桥梁挠度监测方法、电子设备及存储介质 |
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 |