半全息曲面阵列快速成像方法与流程

1.本发明涉及光学成像、微波成像、雷达探测、声呐、超声成像以及基于声、光、电、磁等媒介的目标探测、成像识别、无线通信技术领域，具体涉及一种半全息曲面阵列快速成像方法及其在上述各领域中的应用。


背景技术：

2.当采用阵列单元同时收发信号的一维曲线阵，围绕某一直线进行旋转扫描成像或步进扫描成像时，此时在空间形成曲面探测阵列，此探测阵列在不同的方向具有不同的成像扫描属性，因而构成了一种不同于常规阵列的半全息曲面阵列系统。此外，对于机载、舰载和某些特殊应用场景，需要采用曲面阵列或共形阵列，当此曲面阵列沿不同方向采用不同的成像扫描体制时，也会构成半全息曲面阵列系统。此半全息曲面阵列系统，既无法使用传统的数字波束合成技术，又无法使用现有的全息成像技术，成像难度大，工程实现困难。
3.本发明人曾开发出了相关的快速成像技术，如适用于被动成像与主动成像的快速成像方法(中国专利申请号：202111123446x)、曲面阵列快速成像方法(中国专利申请号：2021112640085)等。
4.需要说明的是，在上述“曲面阵列快速成像方法”中，我们给出了一种曲面阵列快速成像方法，其中：
5.φf为聚焦相位加权系数：
[0006][0007]
φc为曲面相位补偿系数：
[0008]
φc＝ηkδz；
[0009]
其中，r为目标斜距，η为阵列属性参数，对于被动成像系统、半主动成像系统、常规相控阵系统，取η＝1；对于主动全息成像系统，取η＝2。
[0010]
然而，在上述发明申请中，仅涉及了1个阵列属性参数η，给出了具有单一阵列属性的曲面阵列快速成像实现方案。但是，当此曲面阵列在不同方向具有不同的阵列属性时，上述快速成像方法则无法完成有效成像。
[0011]
因此，针对半全息曲面阵列成像系统，需要开发一种兼容性好、成像效果出色的新型成像方法。


技术实现要素：

[0012]
为了解决半全息曲面阵列成像难题，我们进行了系统深入的研究，并获得了一种新的切实可行的解决方案。
[0013]
半全息成像系统具有两个互不相同的阵列属性参数η
x
和ηy，某一维度方向的阵列单元同时收发信号时对应的阵列属性参数值为1，某一维度方向的阵列单元依次收发信号时对应的阵列属性参数值为2。
[0014]
如附图1所示，建立成像系统的坐标系，p为目标，q为目标的像，假设曲面阵列的中
心位于x，y，z坐标系的原点，并且与xy平面相切，标记表示收发阵列单元。
[0015]
对于目标到阵列单元的传播相移φ1、等效透镜相移φ
l
以及阵列单元到成像平面的传播相移φ2，根据变量x、y和z进行分解，获得与阵列方向相关联的分量：
[0016][0017]
其中，为波数，λ为波长，u为物距，即目标到阵列平面的距离，v为像距，即成像平面到阵列平面的距离，f为透镜焦距；(ζ，ξ)为目标坐标，(x，y)为阵列单元坐标，(δ，σ)为像点坐标。
[0018]
为了便于在阵列两个维度方向进行分离，采用如下公式进行成像处理：
[0019][0020]
其中，为像场分布，为目标散射信号，j为虚数单位，e为欧拉常数。
[0021]
上式化简整理后得：
[0022][0023]
其中，其中，
[0024]
当满足成像条件时有此时有φ
x
＝0、φy＝0，代入上式化简得：
[0025][0026]
研究表明，阵面孔径越大，上式右边的二重积分结果越接近狄拉克函数(dirac)的形态，此时所获得的像与源场之间存在如下近似线性关系：
[0027][0028]
上式表明，此时半全息曲面阵列成像系统实现了目标的近似线性成像。
[0029]
对于实际的阵列成像系统，成像时则需要对阵列接收到的信号作如下处理：
[0030][0031]
其中，e为阵列单元接收到的目标信号，a为幅度加权系数，代入相关表达式后化简整理得：
[0032][0033]
其中，其中，
[0034]
当满足成像条件时，有ψ
x
＝0、ψy＝0，上式可简化为：
[0035][0036]
将上式离散化：
[0037][0038]
令xm＝x0 mδ
x
，yn＝y0 nδy，带入上式化简整理得：
[0039][0040]
其中，m、n分别为阵列单元x方向与y方向的序号，(x0，y0)为阵列起始单元的坐标，δ
x
、δy分别为x方向、y方向的阵列单元间距，e
mn
为阵列单元接收到的目标信号，a
mn
为阵列单元的幅度加权系数，m为x方向的阵列单元数量，n为y方向的阵列单元数量，单元的幅度加权系数，m为x方向的阵列单元数量，n为y方向的阵列单元数量，
[0041]
上式右边的系数满足反映了像场的空间波动特性，对成像基本无影响，可忽略。不考虑系数的影响，求和运算可用二维ifft进行快速求解，则像场计算公式为：
[0042][0043]
其中，为最终获得的谱域像，ifft表示二维快速傅里叶逆变换。
[0044]
在上述认识的基础上，本发明提出了一种半全息曲面阵列快速成像方法，本方法基于透镜成像原理，根据天线阵列接收到的目标信号，通过单元信号的幅度、相位加权，采用高效并行算法，获得目标对应的像场分布。
[0045]
具体而言，本发明提供了一种半全息曲面阵列快速成像方法，包括以下步骤：
[0046]
步骤一：确定半全息曲面阵列的两个互不相同的扫描属性参数值；
[0047]
步骤二：对阵列信号进行曲面相位补偿；
[0048]
步骤三：根据上述两个互不相同的扫描属性参数值对阵列信号进行聚焦相位加权；
[0049]
步骤四：根据上述两个互不相同的扫描属性参数值对阵列信号进行扫描相位加权；
[0050]
步骤五：对阵列信号进行幅度加权；
[0051]
步骤六：采用高效并行算法，对阵列信号进行快速成像处理；
[0052]
步骤七：根据上述两个互不相同的扫描属性参数值对像场进行坐标反演。
[0053]
进一步地，本发明半全息曲面阵列快速成像方法步骤一中所述确定半全息曲面阵列的两个互不相同的扫描属性参数值，包括：
[0054]
在半全息曲面阵列系统中，阵列在不同方向的扫描属性参数η
x
、ηy的取值不同；
[0055]
若阵列某一方向的阵列单元同时发射和接收目标信号，则选择该方向的扫描属性参数η
x
或ηy的值为1；
[0056]
若阵列某一方向的阵列单元依次发射和接收目标信号，则选择该方向的扫描属性参数η
x
或ηy的值为2。
[0057]
进一步地，本发明半全息曲面阵列快速成像方法步骤二中所述对阵列信号进行曲面相位补偿，其计算公式为：
[0058][0059]
其中，j为虚数单位，e为欧拉常数，s0为阵列回波信号，s1为曲面相位补偿后的阵列信号，φc为曲面相位补偿系数；
[0060]
其中，曲面相位补偿系数的计算公式为：
[0061]
φc＝2kz；
[0062]
其中，为波数，λ为波长，π为圆周率，z表示曲面阵列单元z方向的直角坐标。
[0063]
进一步地，本发明半全息曲面阵列快速成像方法步骤三中所述根据上述两个互不相同的扫描属性参数值对阵列信号进行聚焦相位加权，其计算公式为：
[0064][0065]
其中，s2为聚焦相位加权后的阵列信号，为阵列x方向的聚焦相位加权值，为阵列y方向的聚焦相位加权值；
[0066]
其中，阵列x方向的聚焦相位加权值的计算方法如下：
[0067][0068]
阵列y方向的聚焦相位加权值的计算方法如下：
[0069][0070]
其中，r为目标斜距，即目标到阵列中心的距离，(x，y)为曲面阵列单元的坐标。
[0071]
进一步地，本发明半全息曲面阵列快速成像方法步骤四中所述根据上述两个互不
相同的扫描属性参数值对阵列信号进行扫描相位加权，通过扫描相位加权调整成像系统中心视角方向，其计算公式为：
[0072][0073]
其中，s3为扫描相位加权后的阵列信号，为阵列x方向的扫描相位加权值，为阵列y方向的扫描相位加权值；
[0074]
其中，阵列x方向的扫描相位加权值的计算方法如下：
[0075][0076]
阵列y方向的扫描相位加权值的计算方法如下：
[0077][0078]
其中，(θ
ζ
，θ
ξ
)为成像系统的中心视角方向对应的角度坐标，符号sin表示正弦函数。
[0079]
进一步地，本发明半全息曲面阵列快速成像方法步骤五中所述对阵列信号进行幅度加权，其计算公式为：
[0080]
s4＝s3a；
[0081]
其中，s4为幅度加权后的阵列信号，a为幅度加权系数；
[0082]
幅度加权的方法包括但不限于均匀分布、余弦加权、汉明窗、taylor分布、切比雪夫分布及混合加权方法。
[0083]
进一步地，本发明半全息曲面阵列快速成像方法步骤六中所述采用高效并行算法，对阵列信号进行快速成像处理，其计算公式为：
[0084][0085]
其中，q
t
为快速成像处理获得的谱域像，ω
δ
、ω
σ
为谱坐标，符号表示高效并行算法函数；
[0086]
高效并行算法包括但不限于二维或三维fft、ifft、非均匀fft、稀疏fft方法；
[0087]
上述快速成像结果对应的谱坐标ω
δ
、ω
σ
取值范围为：ω
δ
∈[0，2π]、ω
σ
∈[0，2π]，进行fftshift运算后将ω
δ
、ω
σ
取值范围变换为：ω
δ
∈[-π，π]、ω
σ
∈[-π，π]，此时的像是符合实际分布的像：
[0088]
q(ω
δ
，ω
σ
)＝fftshift[q
t
(ω
δ
，ω
σ
)]。
[0089]
进一步地，本发明半全息曲面阵列快速成像方法步骤七中所述根据上述两个互不相同的扫描属性参数值对像场进行坐标反演，包括：
[0090]
对高效并行算法获得的像场进行坐标解算，并对像场进行坐标反演，获得真实目标的分布情况；其中：
[0091]
对于ifft类的高效并行算法，像场扫描角坐标计算公式为：
[0092][0093]
对于fft类的高效并行算法，像场扫描角坐标计算公式为：
[0094][0095]
其中，θ
δ
、θ
σ
分别为x方向和y方向的像场扫描角坐标，δ
x
、δy分别为阵列单元x方向和y方向的单元间距，符号sin-1
表示反正弦函数；
[0096]
像场的直角坐标计算公式为：
[0097][0098]
其中，δ、σ分别为x方向和y方向的像场直角坐标，符号tan表示正切函数，p为成像平面到阵列所在平面的距离；
[0099]
真实目标的直角坐标反演计算公式为：
[0100][0101]
进一步地，本发明半全息曲面阵列快速成像方法中设置阵列单元间距进一步地，本发明半全息曲面阵列快速成像方法中设置阵列单元间距以避免出现成像混叠现象；
[0102]
当阵列单元实际间距不满足上述条件时，采用插值方法对数据进行处理以保证满足成像无混叠条件。
[0103]
此外，本发明还提供了一种半全息曲面阵列快速成像方法，所述快速成像方法用于远距离成像，包括：
[0104]
适用于远距离成像的计算公式为：
[0105][0106]
采用前述的高效并行算法计算像场，通过一次运算获得宽视角范围内的目标分布情况。
[0107]
此外，本发明还涉及上述快速成像方法在光学成像、微波成像、雷达探测、声呐、超声成像以及基于声、光、电、磁的目标探测、成像识别、无线通信领域中的应用。
[0108]
综上，本发明方法解决了具有混合阵列属性的半全息曲面阵列的成像难题，具有良好的应用前景，可广泛应用于以声、光、电、磁等为媒介的目标探测及无线通信领域。当探测媒介为电磁波时，本技术适用于微波成像、雷达探测、无线通信、合成孔径雷达、逆合成孔径雷达；当探测媒介为声波、超声波时，本技术适用于声呐、超声成像、合成孔径声呐；当探测媒介为光时，本技术适用于光学成像、合成孔径光学成像等。
附图说明
[0109]
为了更清楚地说明本发明实施例的技术方案，下面将对本发明实施例描述中需要使用的附图作简要介绍，显而易见地，以下附图仅是本发明中记载的一些实施例，对于本领域技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0110]
图1为本发明成像方法的成像系统坐标系示意图。
[0111]
图2为本发明成像方法的方法流程图。
[0112]
图3为利用本发明成像方法进行的半全息曲面阵列快速成像结果图，其中：左图为目标模型，右图为本发明半全息曲面阵列快速成像结果。
具体实施方式
[0113]
为使本发明的目的、技术方案和优点更加清楚，下面将结合具体实施例及相应的附图对本发明的技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例，本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。
[0114]
同时，应理解，本发明的保护范围并不局限于下述特定的具体实施方案；还应当理解，本发明实施例中使用的术语是为了描述特定的具体实施方案，而不是为了限制本发明的保护范围。
[0115]
实施例1：一种半全息曲面阵列快速成像方法(参见附图1-2)，本方法包括下述步骤：
[0116]
步骤一：确定半全息曲面阵列的两个互不相同的扫描属性参数值，包括：
[0117]
在半全息曲面阵列系统中，阵列在不同方向的扫描属性参数η
x
、ηy的取值不同；
[0118]
若阵列某一方向的阵列单元同时发射和接收目标信号，则选择该方向的扫描属性参数η
x
或ηy的值为1；
[0119]
若阵列某一方向的阵列单元依次发射和接收目标信号，则选择该方向的扫描属性参数η
x
或ηy的值为2。
[0120]
步骤二：对阵列信号进行曲面相位补偿，其计算公式为：
[0121][0122]
其中，j为虚数单位，e为欧拉常数，s0为阵列回波信号，s1为曲面相位补偿后的阵列信号，φc为曲面相位补偿系数；
[0123]
其中，曲面相位补偿系数的计算公式为：
[0124]
φc＝2kz；
[0125]
其中，为波数，λ为波长，π为圆周率，z表示曲面阵列单元z方向的直角坐标。
[0126]
步骤三：根据上述两个互不相同的扫描属性参数值对阵列信号进行聚焦相位加权，其计算公式为：
[0127][0128]
其中，s2为聚焦相位加权后的阵列信号，为阵列x方向的聚焦相位加权值，为阵列y方向的聚焦相位加权值；
[0129]
其中，阵列x方向的聚焦相位加权值的计算方法如下：
[0130]
[0131]
阵列y方向的聚焦相位加权值的计算方法如下：
[0132][0133]
其中，r为目标斜距，即目标到阵列中心的距离，(x，y)为曲面阵列单元的坐标。
[0134]
步骤四：根据上述两个互不相同的扫描属性参数值对阵列信号进行扫描相位加权，通过扫描相位加权调整成像系统中心视角方向，其计算公式为：
[0135][0136]
其中，s3为扫描相位加权后的阵列信号，为阵列x方向的扫描相位加权值，为阵列y方向的扫描相位加权值；
[0137]
其中，阵列x方向的扫描相位加权值的计算方法如下：
[0138][0139]
阵列y方向的扫描相位加权值的计算方法如下：
[0140][0141]
其中，(θ
ζ
，θ
ξ
)为成像系统的中心视角方向对应的角度坐标，符号sin表示正弦函数。
[0142]
步骤五：对阵列信号进行幅度加权，其计算公式为：
[0143]
s4＝s3a；
[0144]
其中，s4为幅度加权后的阵列信号，a为幅度加权系数；
[0145]
幅度加权的方法包括但不限于均匀分布、余弦加权、汉明窗、taylor分布、切比雪夫分布及混合加权方法。
[0146]
步骤六：采用高效并行算法，对阵列信号进行快速成像处理，其计算公式为：
[0147][0148]
其中，q
t
为快速成像处理获得的谱域像，ω
δ
、ω
σ
为谱坐标，符号表示高效并行算法函数；
[0149]
高效并行算法包括但不限于二维或三维fft、ifft、非均匀fft、稀疏fft方法；
[0150]
上述快速成像结果对应的谱坐标ω
δ
、ω
σ
取值范围为：ω
δ
∈[0，2π]、ω
σ
∈[0，2π]，进行fftshift运算后将ω
δ
、ω
σ
取值范围变换为：ω
δ
∈[-π，π]、ω
σ
∈[-π，π]，此时的像是符合实际分布的像：
[0151]
q(ω
δ
，ω
σ
)＝fftshift[q
t
(ω
δ
，ω
σ
)]。
[0152]
步骤七：根据上述两个互不相同的扫描属性参数值对像场进行坐标反演，包括：对高效并行算法获得的像场进行坐标解算，并对像场进行坐标反演，获得真实目标的分布情况；其中：
[0153]
对于ifft类的高效并行算法，像场扫描角坐标计算公式为：
[0154]
[0155]
对于fft类的高效并行算法，像场扫描角坐标计算公式为：
[0156][0157]
其中，θ
δ
、θ
σ
分别为x方向和y方向的像场扫描角坐标，δ
x
、δy分别为阵列单元x方向和y方向的单元间距，符号sin-1
表示反正弦函数；
[0158]
像场的直角坐标计算公式为：
[0159][0160]
其中，δ、σ分别为x方向和y方向的像场直角坐标，符号tan表示正切函数，v为成像平面到阵列所在平面的距离；
[0161]
真实目标的直角坐标反演计算公式为：
[0162][0163]
此外，本发明半全息曲面阵列快速成像方法中设置阵列单元间距此外，本发明半全息曲面阵列快速成像方法中设置阵列单元间距以避免出现成像混叠现象；当阵列单元实际间距不满足上述条件时，采用插值方法对数据进行处理以保证满足成像无混叠条件。
[0164]
实施例2：本成像方法(实施例1方法)用于半全息曲面阵列快速成像效果验证试验
[0165]
试验条件：目标为“f”形理想点目标集合，位于阵列法线方向，距离阵列中心1m，点目标之间的间距为15mm。探测信号频率为10ghz，阵列单元位于半径为1m的球面上，其中，水平方向(x方向)天线单元间距为λ/4，竖直方向(y方向)天线单元间距为λ/2，阵列孔径为0.64m
×
0.64m，采用本发明方法进行半全息曲面阵列成像，仿真时选取η
x
＝1、ηy＝2，成像结果见附图3。
[0166]
实施例3：一种半全息曲面阵列快速成像方法，该方法用于远距离成像，选取r＝∞，适用于远距离成像的计算公式为：
[0167][0168]
采用实施例1中的高效并行算法计算像场，通过一次运算获得宽视角范围内的目标分布情况。
[0169]
本发明中的各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。
[0170]
以上所述仅为本发明的实施例而已，并不用于限制本发明。对于本领域技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原理之内所作的任何修改、替换等，均应包含在本发明的权利要求保护范围之内。