[发明专利]点扩散函数估计及高分辨处理方法、装置、设备及介质

权利要求书

1.一种点扩散函数估计及高分辨处理方法，其特征在于，包括：

对成像剖面进行构造属性分析与图像分块，获得局部成像块；

基于所述局部成像块，采用图像学处理方式对数据进行平滑滤波，获得对应的波数域点扩散函数；

基于所述波数域点扩散函数，通过高维反褶积算子进行成像剖面的高分辨率处理。

2.根据权利要求1所述的点扩散函数估计及高分辨处理方法，其中，分块后的局部成像块为：

其中，I_b(x_i,j)为局部成像块，I(x)表示叠后成像剖面，B(·)表示图像分块函数。

3.根据权利要求1所述的点扩散函数估计及高分辨处理方法，其中，所述波数域点扩散函数为：

PSF_b(x_i,j)＝AFilt(S(FK(I_b(x_i,j))))           (2)

其中，PSF_b(x_i,j)为波数域点扩散函数，I_b(x_i,j)为局部成像块，FK(.)为FK变换，S(.)为常规二维平滑算子，AFilt(.)为有照明角度范围约束的角度滤波算子。

4.根据权利要求3所述的点扩散函数估计及高分辨处理方法，其中，通过公式(3)计算角度滤波算子：

AFilt(x_i,j)＝AFilt(angle_l,angle_r,x_i,j,angle_k,f(x))    (3)

其中，angle_l、angle_r为局部成像块FK谱对应的点扩散函数的左右两方向的照明范围约束，angle_k为在照明范围内的角度数据，f(x)为滤波函数。

5.根据权利要求4所述的点扩散函数估计及高分辨处理方法，其中，通过f(x)对angle_k进行一维滤波。

6.根据权利要求5所述的点扩散函数估计及高分辨处理方法，其中，滤波函数f(x)为：

f(x)＝0.54-0.46*cos(2*π*(ntaper-x)/ntaper)        (4)

其中，ntaper表示滤波算子长度。

7.根据权利要求1所述的点扩散函数估计及高分辨处理方法，其中，通过公式(5)进行成像剖面的高分辨率处理：

其中，R(x_i,j)为高分辨率处理结果，FK^-1(.)为FK逆变换，α为阻尼因子，表示PSF_b(.)的共轭转置矩阵。

8.一种点扩散函数估计及高分辨处理装置，其特征在于，包括：

分块模块，对成像剖面进行构造属性分析与图像分块，获得局部成像块；

平滑滤波模块，基于所述局部成像块，采用图像学处理方式对数据进行平滑滤波，获得对应的波数域点扩散函数；

高分辨率处理模块，基于所述波数域点扩散函数，通过高维反褶积算子进行成像剖面的高分辨率处理。

9.一种电子设备，其特征在于，所述电子设备包括：

存储器，存储有可执行指令；

处理器，所述处理器运行所述存储器中的所述可执行指令，以实现权利要求1-7中任一项所述的点扩散函数估计及高分辨处理方法。

10.一种计算机可读存储介质，其特征在于，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现权利要求1-7中任一项所述的点扩散函数估计及高分辨处理方法。