本公开涉及图像处理技术领域，具体地，涉及一种图像处理方法、装置、设备及计算机可读存储介质。
背景技术：
合成孔径雷达(syntheticapertureradar,sar)是一种可以对地球进行观测的高分辨率遥感探测雷达。由于星载sar能克服雨雾和暗夜条件的影响，实现全天时、全天候对地观测，因此，星载sar在农业、海洋、灾害监测以及3d绘图等领域具有广阔的应用前景。
在实际的星载sar体制中，目标的全部多普勒带宽极大，而受制于整体波位设计和数据率的限制，脉冲重复频率不可能覆盖全部带宽，因此在卫星接收到的回波中，超出脉冲重复频率的多普勒信号会出现反褶，最后在图像中形成方位模糊。目前还没有一种有效的方法可以完美地抑制方位模糊。
技术实现要素：
本公开的目的是提供一种图像处理方法、装置、设备及计算机可读存储介质，以抑制图像中的方位模糊。
为了实现上述目的，本公开第一方面提供一种图像处理方法，包括：
将待处理图像信号转换到频域，以得到初始频域信号；
根据原目标天线方向图函数和模糊目标天线方向图函数，获得非对称加权矩阵；
根据所述非对称加权矩阵，对所述初始频域信号进行非对称加权处理；
将经所述非对称加权处理后的频域信号转换到时域，以得到非对称加权的图像信号。
可选地，根据原目标天线方向图函数和模糊目标天线方向图函数，获得非对称加权矩阵，包括：
按照以下公式，建立m×1的加权数组f(m)：
其中，floor(·)为向下取整操作符，ceil(·)为向上取整操作符，m＝0,1,…,m-1，m为纵向索引值，m为所述待处理图像信号对应的待处理图像的纵向像素总数，ρ为频谱过采样系数，wmain(m)为所述原目标天线方向图函数，wamb(m)为所述模糊目标天线方向图函数，t为预设门限值；
对所述加权数组f(m)进行归一化，得到归一化数组f1(m)；
建立m×n的非对称加权矩阵wf(m,n)，所述非对称加权矩阵的每一列满足wf(m,n)＝f1(m),n＝0,1,…,n-1，n为横向索引值，n为所述待处理图像信号对应的待处理图像的横向像素总数。
可选地，对所述加权数组f(m)进行归一化，得到归一化数组f1(m)，包括：
按照以下公式确定归一化常数para：
其中，sum(·)为整体求和操作符；
根据所述加权数组f(m)和所述归一化常数para，按照以下公式，得到所述归一化数组f1(m)：
f1(m)＝f(m)·para。
可选地，所述方法还包括：
按照以下公式，生成所述原目标天线方向图函数wmain(m)：
按照以下公式，生成所述模糊目标天线方向图函数wamb(m)：
其中，vr为星地等效速度，la为天线长度，prf为脉冲重复频率，fdc为多普勒中心频率。
本公开第二方面提供一种图像处理装置，包括：
第一转换模块，被配置为将待处理图像信号转换到频域，以得到初始频域信号；
获得模块，被配置为根据原目标天线方向图函数和模糊目标天线方向图函数，获得非对称加权矩阵；
处理模块，被配置为根据所述非对称加权矩阵，对所述初始频域信号进行非对称加权处理；
第二转换模块，被配置为将经所述非对称加权处理后的频域信号转换到时域，以得到非对称加权的图像信号。
可选地，所述获得模块包括：
第一建立子模块，被配置为按照以下公式，建立m×1的加权数组f(m)：
其中，floor(·)为向下取整操作符，ceil(·)为向上取整操作符，m＝0,1,…,m-1，m为纵向索引值，m为所述待处理图像信号对应的待处理图像的纵向像素总数，ρ为频谱过采样系数，wmain(m)为所述原目标天线方向图函数，wamb(m)为所述模糊目标天线方向图函数，t为预设门限值；
归一化子模块，被配置为对所述加权数组f(m)进行归一化，得到归一化数组f1(m)；
第二建立子模块，被配置为建立m×n的非对称加权矩阵wf(m,n)，所述非对称加权矩阵的每一列满足wf(m,n)＝f1(m),n＝0,1,…,n-1，n为横向索引值，n为所述待处理图像信号对应的待处理图像的横向像素总数。
可选地，所述归一化子模块包括：
确定子模块，被配置为按照以下公式确定归一化常数para：
其中，sum(·)为整体求和操作符；
获得子模块，被配置为根据所述加权数组f(m)和所述归一化常数para，按照以下公式，得到所述归一化数组f1(m)：
f1(m)＝f(m)·para。
可选地，所述装置还包括：
第一生成模块，被配置为按照以下公式，生成所述原目标天线方向图函数wmain(m)：
第二生成模块，被配置为按照以下公式，生成所述模糊目标天线方向图函数wamb(m)：
其中，vr为星地等效速度，la为天线长度，prf为脉冲重复频率，fdc为多普勒中心频率。
本公开第三方面提供一种图像处理设备，包括：至少一个处理器、至少一个存储器以及存储在所述存储器中的计算机程序指令，当所述计算机程序指令被所述处理器执行时实现本公开第一方面所述的方法。
本公开第四方面提供一种计算机可读存储介质，其上存储有计算机程序指令，当所述计算机程序指令被处理器执行时实现如本公开第一方面所述的方法。
通过上述技术方案，通过非对称加权矩阵对频域信号进行处理，在主频谱占优的频点上提供较高的频域通过性，在模糊频谱占优的频点上降低通过性，既实现了模糊信号的抑制处理，又尽可能地保留原目标特性。采用本公开实施例提供的方法，只需要少量参数(即上文中的相关参数)即可完成处理，不需要主成像区域数据，可适用于各种模式的星载合成孔径雷达图像模糊抑制。并且，该方法计算速度较快，适合对处理实时性要求高的应用场景。
本公开的其他特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本公开的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本公开，但并不构成对本公开的限制。在附图中：
图1是本公开实施例提供的图像处理方法的流程图。
图2是本公开实施例中待处理图像信号对应的待处理图像的示意图。
图3是采用本公开实施例提供的图像处理方法得到的结果图像的示意图。
图4是本公开实施例提供的图像处理装置的示意图。
图5是本公开实施例提供的图像处理设备的示意图。
具体实施方式
以下结合附图对本公开的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本公开，并不用于限制本公开。
本公开实施例提供了一种图像处理方法，该方法利用原目标天线方向图函数和模糊目标天线方向图函数的不同，从频域信号入手，利用非对称加权矩阵最大化地抑制图像中的方位模糊。
图1是本公开实施例提供的图像处理方法的流程图，如图1所示，该方法包括以下步骤：
在步骤s11中，将待处理图像信号转换到频域，以得到初始频域信号；
在步骤s12中，根据原目标天线方向图函数和模糊目标天线方向图函数，获得非对称加权矩阵；
在步骤s13中，根据所述非对称加权矩阵，对所述初始频域信号进行非对称加权处理；
在步骤s14中，将经所述非对称加权处理后的频域信号转换到时域，以得到非对称加权的图像信号。
本公开实施例中，待处理图像信号对应的待处理图像即为具有方位模糊的图像，例如：仿真雷达卫星图像。图2是本公开实施例中待处理图像信号对应的待处理图像的示意图。由图2可见，整个区域完全被模糊占据，无法判别出目标。
本公开实施例提供的图像处理方法的实施过程如下：
第一步，读取待处理图像信号s(m,n)，其中，m＝0,1,…,m-1，m为纵向索引值，m为待处理图像信号对应的待处理图像的纵向像素总数，n＝0,1,…,n-1，n为横向索引值，n为待处理图像信号对应的待处理图像的横向像素总数。
第二步，按照公式(1)，将待处理图像信号转换到频域，以得到初始频域信号s(k,n)：
第三步，根据相关参数和天线情况，生成原目标天线方向图函数和模糊目标天线方向图函数。相关参数包括：星地等效速度vr，天线长度la，脉冲重复频率prf，多普勒中心频率fdc，频谱过采样系数ρ。示例地，相关参数对照terrasar-x真实波位，如表1所示。
表1相关参数
按照公式(2)，生成原目标天线方向图函数wmain(m)：
按照公式(3)，生成模糊目标天线方向图函数wamb(m)：
第四步，构建非对称加权矩阵，具体包括以下步骤：
步骤(1)，按照公式(4)，建立m×1的加权数组f(m)：
本公开实施例中，floor(·)为向下取整操作符，ceil(·)为向上取整操作符，t为预设门限值。一般情况下可设置t＝500。
具体地，首先，建立m×1的加权数组f(m)。
然后，针对以下索引位置：
令f(m)＝0；
针对以下索引位置：
若wamb(m)＝0，则令f(m)＝t；
若wamb(m)≠0且
若wamb(m)≠0且
步骤(2)，对加权数组f(m)进行归一化，得到归一化数组f1(m)。一种实施方式是：
首先，按照公式(5)确定归一化常数para：
本公开实施例中，sum(·)为整体求和操作符。
然后，按照公式(6)，得到归一化数组f1(m)：
f1(m)＝f(m)·para公式(6)
步骤(3)，建立m×n的非对称加权矩阵wf(m,n)，非对称加权矩阵wf(m,n)的每一列满足wf(m,n)＝f1(m)。
第五步，按照公式(7)，对初始频域信号s(k,n)进行非对称加权处理，得到经非对称加权处理后的频域信号ws(k,n)：
ws(k,n)＝s(k,n)·wf(m,n)公式(7)
其中，k＝0,1,…,n-1，k与n均为横向索引值。
第六步，按照公式(8)，将经非对称加权处理后的频域信号ws(k,n)转换到时域，以得到非对称加权的图像信号ws(m,n)：
非对称加权的图像信号对应的图像即为结果图像。图3是采用本公开实施例提供的图像处理方法得到的结果图像的示意图。如图3所示，被模糊区域覆盖的目标凸显出来，达到了模糊抑制效果。
采用上述技术方案，通过非对称加权矩阵对频域信号进行处理，在主频谱占优的频点上提供较高的频域通过性，在模糊频谱占优的频点上降低通过性，既实现了模糊信号的抑制处理，又尽可能地保留原目标特性。采用本公开实施例提供的方法，只需要少量参数(即上文中的相关参数)即可完成处理，不需要主成像区域数据，可适用于各种模式的星载合成孔径雷达图像模糊抑制。并且，该方法计算速度较快，适合对处理实时性要求高的应用场景。
基于同一发明构思，本公开实施例还提供一种图像处理装置。图4是本公开实施例提供的图像处理装置的示意图。如图4所示，该装置200包括：
第一转换模块201，被配置为将待处理图像信号转换到频域，以得到初始频域信号；
获得模块202，被配置为根据原目标天线方向图函数和模糊目标天线方向图函数，获得非对称加权矩阵；
处理模块203，被配置为根据所述非对称加权矩阵，对所述初始频域信号进行非对称加权处理；
第二转换模块204，被配置为将经所述非对称加权处理后的频域信号转换到时域，以得到非对称加权的图像信号。
可选地，所述获得模块包括：
第一建立子模块，被配置为按照以下公式，建立m×1的加权数组f(m)：
其中，floor(·)为向下取整操作符，ceil(·)为向上取整操作符，m＝0,1,…,m-1，m为纵向索引值，m为所述待处理图像信号对应的待处理图像的纵向像素总数，ρ为频谱过采样系数，wmain(m)为所述原目标天线方向图函数，wamb(m)为所述模糊目标天线方向图函数，t为预设门限值；
归一化子模块，被配置为对所述加权数组f(m)进行归一化，得到归一化数组f1(m)；
第二建立子模块，被配置为建立m×n的非对称加权矩阵wf(m,n)，所述非对称加权矩阵的每一列满足wf(m,n)＝f1(m),n＝0,1,…,n-1，n为横向索引值，n为所述待处理图像信号对应的待处理图像的横向像素总数。
可选地，所述归一化子模块包括：
确定子模块，被配置为按照以下公式确定归一化常数para：
其中，sum(·)为整体求和操作符；
获得子模块，被配置为根据所述加权数组f(m)和所述归一化常数para，按照以下公式，得到所述归一化数组f1(m)：
f1(m)＝f(m)·para。
可选地，所述装置还包括：
第一生成模块，被配置为按照以下公式，生成所述原目标天线方向图函数wmain(m)：
第二生成模块，被配置为按照以下公式，生成所述模糊目标天线方向图函数wamb(m)：
其中，vr为星地等效速度，la为天线长度，prf为脉冲重复频率，fdc为多普勒中心频率。
本公开实施例提供了一种图像处理设备，包括：至少一个处理器、至少一个存储器以及存储在存储器中的计算机程序指令，当计算机程序指令被处理器执行时实现如上述实施方式中图像处理方法。
本公开实施例提供了一种计算机可读存储介质，其上存储有计算机程序指令，当计算机程序指令被处理器执行时实现如上述实施方式中图像处理方法。
结合图1描述的本公开实施例的图像处理方法可以由图像处理设备来实现。图5示出了本公开实施例提供的图像处理设备的硬件结构示意图。
图像处理设备可以包括处理器301以及存储有计算机程序指令的存储器302。
具体地，上述处理器301可以包括中央处理器(cpu)，或者特定集成电路(applicationspecificintegratedcircuit，asic)，或者可以被配置成实施本公开实施例的一个或多个集成电路。
存储器302可以包括用于数据或指令的大容量存储器。举例来说而非限制，存储器302可包括硬盘驱动器(harddiskdrive，hdd)、软盘驱动器、闪存、光盘、磁光盘、磁带或通用串行总线(universalserialbus，usb)驱动器或者两个或更多个以上这些的组合。在合适的情况下，存储器302可包括可移除或不可移除(或固定)的介质。在合适的情况下，存储器302可在数据处理装置的内部或外部。在特定实施例中，存储器302是非易失性固态存储器。在特定实施例中，存储器302包括只读存储器(rom)。在合适的情况下，该rom可以是掩模编程的rom、可编程rom(prom)、可擦除prom(eprom)、电可擦除prom(eeprom)、电可改写rom(earom)或闪存或者两个或更多个以上这些的组合。
处理器301通过读取并执行存储器302中存储的计算机程序指令，以实现上述实施例中的任意一种图像处理方法。
在一个示例中，图像处理设备还可包括通信接口303和总线310。其中，如图5所示，处理器301、存储器302、通信接口303通过总线310连接并完成相互间的通信。
通信接口303，主要用于实现本公开实施例中各模块、装置、单元和/或设备之间的通信。
总线310包括硬件、软件或两者，将执行图像处理方法的设备的部件彼此耦接在一起。举例来说而非限制，总线可包括加速图形端口(agp)或其他图形总线、增强工业标准架构(eisa)总线、前端总线(fsb)、超传输(ht)互连、工业标准架构(isa)总线、无限带宽互连、低引脚数(lpc)总线、存储器总线、微信道架构(mca)总线、外围组件互连(pci)总线、pci-express(pci-x)总线、串行高级技术附件(sata)总线、视频电子标准协会局部(vlb)总线或其他合适的总线或者两个或更多个以上这些的组合。在合适的情况下，总线310可包括一个或多个总线。尽管本公开实施例描述和示出了特定的总线，但本公开考虑任何合适的总线或互连。
另外，结合上述实施例中的图像处理方法，本公开实施例可提供一种计算机可读存储介质来实现。该计算机可读存储介质上存储有计算机程序指令；该计算机程序指令被处理器执行时实现上述实施例中的任意一种图像处理方法。
需要明确的是，本公开并不局限于上文所描述并在图中示出的特定配置和处理。为了简明起见，这里省略了对已知方法的详细描述。在上述实施例中，描述和示出了若干具体的步骤作为示例。但是，本公开的方法过程并不限于所描述和示出的具体步骤，本领域的技术人员可以在领会本公开的精神后，作出各种改变、修改和添加，或者改变步骤之间的顺序。
以上所述的结构框图中所示的功能块可以实现为硬件、软件、固件或者它们的组合。当以硬件方式实现时，其可以例如是电子电路、专用集成电路(asic)、适当的固件、插件、功能卡等等。当以软件方式实现时，本公开的元素是被用于执行所需任务的程序或者代码段。程序或者代码段可以存储在机器可读介质中，或者通过载波中携带的数据信号在传输介质或者通信链路上传送。“机器可读介质”可以包括能够存储或传输信息的任何介质。机器可读介质的例子包括电子电路、半导体存储器设备、rom、闪存、可擦除rom(erom)、软盘、cd-rom、光盘、硬盘、光纤介质、射频(rf)链路，等等。代码段可以经由诸如因特网、内联网等的计算机网络被下载。
还需要说明的是，本公开中提及的示例性实施例，基于一系列的步骤或者装置描述一些方法或系统。但是，本公开不局限于上述步骤的顺序，也就是说，可以按照实施例中提及的顺序执行步骤，也可以不同于实施例中的顺序，或者若干步骤同时执行。
以上所述，仅为本公开的具体实施方式，所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，上述描述的系统、模块和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。应理解，本公开的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本公开揭露的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本公开的保护范围之内。
以上结合附图详细描述了本公开的优选实施方式，但是，本公开并不限于上述实施方式中的具体细节，在本公开的技术构思范围内，可以对本公开的技术方案进行多种简单变型，这些简单变型均属于本公开的保护范围。
另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合，为了避免不必要的重复，本公开对各种可能的组合方式不再另行说明。
此外，本公开的各种不同的实施方式之间也可以进行任意组合，只要其不违背本公开的思想，其同样应当视为本公开所公开的内容。