CN110687535B - 一种快速微波成像方法 - Google Patents

一种快速微波成像方法 Download PDF

Info

Publication number
CN110687535B
CN110687535B CN201910911179.9A CN201910911179A CN110687535B CN 110687535 B CN110687535 B CN 110687535B CN 201910911179 A CN201910911179 A CN 201910911179A CN 110687535 B CN110687535 B CN 110687535B
Authority
CN
China
Prior art keywords
matrix
function
current
green
inversion
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
Application number
CN201910911179.9A
Other languages
English (en)
Other versions
CN110687535A (zh
Inventor
徐魁文
刘林春
赵鹏
王高峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Fan Li Technology Co ltd
Original Assignee
Hangzhou Fan Li Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hangzhou Fan Li Technology Co ltd filed Critical Hangzhou Fan Li Technology Co ltd
Priority to CN201910911179.9A priority Critical patent/CN110687535B/zh
Publication of CN110687535A publication Critical patent/CN110687535A/zh
Application granted granted Critical
Publication of CN110687535B publication Critical patent/CN110687535B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开一种快速微波成像方法。在DBIM算法的每一次迭代过程中，通过调用这个高效的正问题求解器算正问题来提高反演速度。同时为了使整个反演算法稳定地进行迭代反演并且，采用了截断奇异值分解(TSVD)技术，这避免了大量的数值测试来确定正则化项。为了进一步的减少计算量，在反演算法的成本函数中，用入射场来代替总场。因此该发明能快速，稳定地进行反演，保证了图像的反演速度和质量。

Description

一种快速微波成像方法
技术领域
本发明属于微波成像技术领域，提出了一种快速微波成像方法。
背景技术
微波成像技术以微波在各种复杂媒质中的传播和散射的研究为基础，通过测量被测媒质外部的散射场数据，重建被测媒质内部的复介电常数图像。被测的散射场携带大量有关散射体的信息，利用关于散射目标的先验知识，经过适当的数学处理之后可以提取出散射体本身所具有的某些特性，如散射体的形状，介电常数的分布等。由于微波成像能实现对目标的无损检测，能同时对目标进行几何成像和物理成像，且具有较高分辨率，因此微波成像技术在遥感，医学成像，雷达目标识别等领域有广泛的用途。
微波成像技术一般包含两部分：一部分是正问题计算，即计算给定模型的电场分布；另一部分是逆问题的求解，即根据给定的测量场重构电场的分布。当未知目标的物理尺寸比拟，甚至大于入射波的波长，逆问题就变得非线性，病态，并且计算量巨大。这些困难在处理三维微波成像问题中变得尤其棘手。根据等效电流是否包含在内，反演方法可大致分为两类：场型反演方法和源型反演方法。场型反演方法主要包括变型波恩迭代法(DBIM)和牛顿迭代法等。在DBIM方法中，每次需要更新非均匀背景下的格林函数，并且采用了波恩近似用入射场来代替总场。源型反演方法主要包括对比源方法(CSI)和子空间优化方法(SOM)等。在SOM方法中，等效感应电流被划分为确定性电流和非确定性电流。在SOM方法中，通过对确定性电流进行一定的处理，提出了一种基于子空间的DBIM方法(S-DBIM)来加速迭代收敛。在S-DBIM 的迭代过程中，需要更新背景非均匀格林函数和总场，此时总场不再由入射场代替，而是入射场和散射场的总和。对总场赋予更精确的物理含义，这使得S-DBIM收敛更快。但是在S-DBIM中，非均匀格林函数和总场都需要更新，S-DBIM在处理三维反演问题时计算量变得很大。
自从周永祖教授首次提出DBIM反演算法，人们提出了各种改进的DBIM算法来解决不同的成像问题，比如乳腺成像、密度成像。由于DBIM算法反演速度不够快，我们将其运用到实际问题上还是存在诸多的困难，特别是在三维的微波成像问题上。在DBIM算法中，每次迭代我们都需要计算一次正问题。反演的速度在很大程度上取决于正问题的计算速度。因此人们将不同的高效正问题求解器集成到DBIM 反演算法中来加快正问题的求解速度。稳定双共轭梯度快速傅立叶变换(BCGS-FFT) 曾作为正问题求解器来反演多层介质问题。但是在这种方法中，均匀离散化的必要性使得这种方法处理复杂几何物体变得很吃力。在2009年，国外学者提出将多层快速多级子(MLFMA)和DBIM相结合来。然而MLFMA这种高效正问题求解器在处理多层介质问题的反演时，变得很棘手。另一方面，之后报导的一种叫作多层格林函数插值方法(MLGFIM)被用来解决各种电磁正问题。这种方法采用了BCGS-FFT的插值思想和MLFMA的多层结构。因此这个算法作为一个可行性高的正问题求解器 MLGFIM-Slover是一个很好的选择。
发明内容
本发明所要解决的技术问题是传统的三维微波成像方法反演速度慢的问题，提出了一种快速微波成像方法，这种方法将修正的MS-DBIM,基于变形波恩近似(DBIM) 反演算法的基础上，集成了一个高效的正问题求解器，该求解器采用了一种与积分核无关的插值算法，即多层格林函数插值算法(MLGFIM)，可以将这个求解器命名为 MLGFIM-Solver.其优点是在DBIM算法的每一次迭代过程中，通过调用这个高效的正问题求解器，我们能相对快速地计算正问题来提高反演速度。同时为了使整个反演算法稳定地进行迭代反演并且，采用了截断奇异值分解(TSVD)技术，这避免了大量的数值测试来确定正则化项。为了进一步的减少计算量，在反演算法的成本函数中，用入射场来代替总场。因此该发明能快速，稳定地进行反演，保证了图像的反演速度和质量。
本发明的技术方案：
本发明设计方法利用发射天线照射电磁波到探测区域，用接收天线获取到一定量的散射场数据后，在电磁场积分方程的基础上通过离散化探测区域构建相对应的矩阵方程和成本函数方程，用发明的MS-DBIM-MLGFIM方法来快速迭代求解成本函数的最小值，同时为了使整个算法在迭代过程能稳定地进行反演，我们采用了截断奇异值分解技术来代替繁琐的正则化参数选取，具体过程如下:
本发明快速微波成像方法，包括以下步骤：
步骤(1)、探测区域外设置Ninc个发射天线和Nr个接收天线，发射天线发射频率f的电磁波到有未知散射体的探测区域，接收天线接收对应的测量散射场
Figure BDA0002214747460000021
步骤(2)、首先将探测区域离散为M个足够小的离散方块，并根据已知的电磁波频率f、探测区域位置、发射天线和接收天线的位置和极化方式、测量散射场
Figure BDA0002214747460000031
计算格林函数矩阵
Figure BDA0002214747460000032
和离散后入射场
Figure BDA0002214747460000033
并构建MLGFIM-solve求解器。
格林函数矩阵
Figure BDA0002214747460000034
为离散后的格林函数
Figure BDA0002214747460000035
的积分算子；
Figure BDA0002214747460000036
为离散后的入射场矩阵其中并矢格林函数
Figure BDA0002214747460000037
为并矢格林函数,表示探测区域内点r′处的点源对其接收天线rs处的作用，是一个3×3的矩阵；标量格林函数 g(rs,r′)＝exp(ik0|rs-r′|)/(4π|rs-r′|)，k0为背景介质的波数。
Figure BDA0002214747460000038
是一个3×3的单位对角矩阵，
Figure BDA0002214747460000039
表示旋度算子，i表示复数的虚部。
MLGFIM-solve求解器是对
Figure BDA00022147474600000310
(矩阵向量相乘操作进行加速；其中
Figure BDA00022147474600000311
代表任意一个向量)
Figure BDA00022147474600000312
为离散化的格林函数
Figure BDA00022147474600000313
的积分算子；并矢格林函数
Figure BDA00022147474600000314
表示的是探测区域内点r′对探测区域内另一点r的作用；
Figure BDA00022147474600000315
是一个3M×3M的矩阵，可由9个M×M的矩阵
Figure BDA00022147474600000316
组成,即
Figure BDA00022147474600000317
当m≠n,(n,m＝1,2,…,M)
Figure BDA00022147474600000318
当m＝n，
Figure BDA00022147474600000319
其中w分别是电磁波的角频率；μ0是背景介质的磁导率；ρ(u-v)表示冲激函数，当u＝v时ρ(u-v)＝1，当u≠v时ρ(u-v)＝0； Rn,m＝|rn-rm|，rn＝(r1:n,r2；n,r3；n)，rm＝(r1:m,r2；m,r3；m)分别是第n,m个离散方块的中心坐标,其中1，2，3分别代表三维坐标下的x,y,z方向上的分量。
当未知数M很大时,
Figure BDA00022147474600000320
会变成非常大的密集矩阵，如果通过传统的矩量法去填充
Figure BDA00022147474600000321
的每个元素来完成
Figure BDA00022147474600000322
矩阵向量乘法操作计算量很大并且耗费的时间也比较长，因此我们通过构建MLGFIM-solver求解器来完成操作。
步骤(3)、定义新的成本函数
在原始的DBIM算法中，成本函数可以定义为：
Figure BDA0002214747460000041
其中
Figure BDA0002214747460000042
是二次入射电场；
Figure BDA0002214747460000043
是上一迭代后的最新离散后的格林函数
Figure BDA0002214747460000044
的积分算子；δ是散射张量系数,
Figure BDA0002214747460000045
为散射强度张量，为一个3M×3M的对角矩阵,可以表示为
Figure BDA0002214747460000046
当m≠M,2M,3M，则q＝mod(m,M)，mod表示的是取余符号；其他情况下 q＝M；Vq,∈r；q分别是第q个离散方块的体积和介电常数。γ为Tikhonov正则化参数，主要的作用是为了保证反演算法在迭代过程中的稳定性，但是这个参数需要大量的数值测试来确定。
在本发明的MS-DBIM-MLGFIM算法中，采用了截断奇异值分解(TSVD)技术,可以代替正则化参数的作用。对
Figure BDA0002214747460000047
进行奇异值分解：
Figure BDA0002214747460000048
其中
Figure BDA0002214747460000049
表示对
Figure BDA00022147474600000410
奇异值分解得到的左边矩阵的第j列，
Figure BDA00022147474600000411
表示第j个电流基；
Figure BDA00022147474600000412
表示
Figure BDA00022147474600000413
的共轭；将奇异值σj的绝对值按照从大到小的顺序排，取前L个较大奇异值σj对应的电流基构成一个电流子空间，这个电流子空间构成了电流的主要组成部分。即L 是一个临界值，它满足条件|σj|≥阈值，j≤L。上述电流子空间即为感应电流的确定性部分，可以表示为：
Figure BDA00022147474600000414
其中
Figure BDA00022147474600000415
表示
Figure BDA00022147474600000416
的共轭。利用TSVD技术，将归一化处理后的探测区域内的测量散射场
Figure BDA00022147474600000417
和感应电流确定性部分的叠加后得到新的成本函数：
Figure BDA00022147474600000418
其中
Figure BDA00022147474600000419
为每次更新得到的第p个入射场对应的感应电流确定性部分。
步骤(4)、初始化：设定迭代次数
Figure BDA00022147474600000420
步骤(5)、利用散射强度张量通过公式(8)-(9)分别更新二次入射电场
Figure BDA00022147474600000421
和非均匀格林函数
Figure BDA0002214747460000051
将步骤(2)MLGFIM-solver嵌入到共轭梯度迭代法(CG)中来求解
Figure BDA0002214747460000052
Figure BDA00022147474600000510
其中
Figure BDA0002214747460000053
是3M×3M的单位对角矩阵；
步骤(6)、对更新后的格林函数矩阵
Figure BDA0002214747460000054
进行奇异值分解
Figure BDA0002214747460000055
则感应电流相对应的确定性部分可以通过下式得到
Figure BDA0002214747460000056
步骤(7)、然后将
Figure BDA0002214747460000057
代入到公式(7)得到当前迭代的成本函数f,若f小于阈值，则停止迭代,跳转至步骤(9)，否则跳转至步骤8。
步骤(8)、重置n＝n+1通过求公式(7)的最小值来更新
Figure BDA0002214747460000058
根据公式(12) 得到新的散射强度张量；然后跳转至步骤5继续迭代：
Figure BDA0002214747460000059
步骤(9)、根据当前更新得到的散射强度张量数据，在计算机上编程可视化进行成像。
本发明的有益效果：在三维微波成像中，提出了一种快速成像方法，在DBIM 反演迭代过程中，每次计算正问题通过调用MLGFIM-Solver这个高效的正问题求解器。并且在反演算法MS-DBIM中，采用了截断奇异值分解技术，避免了繁琐的正则化参数选取过程,为了进一步的减少计算量，在反演成本函数中，用入射场来代替总场,避免了散射场那部分的计算量。本发明所提出的成像方法有很高的反演速度，并且构造出的图像质量高。
附图说明
图1是本发明方法流程图。
图2是三维微波成像的实验测量装置图。
图3是MLGFIM-solver中二维插值示意图。
图4是实施例1的实物图。
图5是MS-DBIM-MLGFIM对实施例1的反演结果，其中(a1)、(a2)分别为z＝0 截面的实部、虚部成像图；(b1)、(b 2)分别为x＝0截面的实部、虚部成像图；(c 1)、(c2)分别为y＝-0.5λ截面的实部、虚部成像图。
图6是实施例2的实物图。
图7是MS-DBIM-MLGFIM对仿真例2的反演结果；其中(a1)、(a 2)分别为x＝0 截面的实部、虚部成像图；(b1)、(b2)分别为z＝-0.4λ截面的实部、虚部成像图；(c1)、 (c2)分别为y＝-0.5λ截面的实部、虚部成像图。
具体实施方式
下面结合附图和仿真测试例子对本技术做进一步详细说明。
首先我们对三维电磁场正问题进行描述，如图2所示。在三维坐标系下，一个在长方体探测区域内D内的非磁性散射体被Ninc个平面波照射。
假设三维坐标系下任何一点可以表示为r＝(r1,r2,r3)。对于每一个入射天线而言，有Nr个接收天线位于r′n＝(r′1:l,r′2；l,r′3；l),l＝1,2,...,Nr来测量三个方向上的矢量场。在r处的散射场满足下列积分方程：
Figure BDA0002214747460000061
其中
Figure BDA0002214747460000062
为对比度函数，∈(r′)表示离散方块r′的介电常数，∈0表示空气中的介电常数，
Figure BDA0002214747460000063
为三维并矢格林函数，表示一个探测区域r′处的点源对另一点r处的作用，是一个3×3的矩阵；标量格林函数g(r,r′)＝ exp(ik0|r-r′|)/(4π|rs-r′|)，k0为电磁波的波数；Etot(r′)表示了位于r′的总场。 Etot＝Einc+Esca代表了总场是入射场和散射场之和，所以可得下列方程：
Figure BDA0002214747460000064
考虑到计算的简便性和MLGFIM的应用，可将探测区域离散为M个足够小的小立方体，这些小立方体的中心坐标分别为rm＝(r1；m,r2；m,r3；m)，m＝1,2,3,...,M。因此(14)可以转换成一个矩阵方程：
Figure BDA0002214747460000065
其中
Figure BDA0002214747460000066
由
Figure BDA0002214747460000067
组成， T表示对一个矩阵进行转置，这表示入射场在M个小立方体上k方向上的场，1,2,3分别代表三维坐标下的x,y,z方向上的分量。
Figure BDA0002214747460000068
和
Figure BDA0002214747460000069
有着相同的结构，式(15)中其他变量在发明步骤中已经详细阐明，这里就不做说明。
在探测区域外的Nr个接收天线上的散射场可以通过下式得到：
Figure BDA0002214747460000071
其中
Figure BDA0002214747460000072
由
Figure BDA0002214747460000073
组成，这表示散射场在Nr个接收天线上k方向上的场。
当未知数M很大时，
Figure BDA0002214747460000074
将会变成一个很大的密集矩阵，所以我们构建了MLGFIM-solver求解器来降低这个复杂度。MLGFIM-solve采用了分层的思想，首先我们将正方体探测区域D作为MLGFIM-solver第一层的最大立方体，然后把它划分为8个更小的立方体，每个子立方体又继续往下划分，直到划分后的最低层立方体满足条件。
MLGFIM-solver用插值方法来快速地估计
Figure BDA0002214747460000075
矩阵中的元素，如图3所示，立方体m有一源点r′，在立方体n中有一场点r，如果立方体m和立方体 n的中心点距离大于其2倍的边长，白色的空心点为插值点。我们可以通过插值的方法来计算这两点的作用GD；uv(r,r′)，即：
Figure BDA0002214747460000076
其中wm,p，wn,q分别为场立方体m的第p个插值基函数和源立方体n的第q个插值基函数，rm,p,rn,q分别为场立方体m的第p个插值点和源立方体n的第q个插值点，K为插值点数。
通过(17)，场立方体m中所有点和源立方体n中所有点之间的作用可由矩阵
Figure BDA0002214747460000077
表示，其具体表达式为：
Figure BDA0002214747460000078
其中
Figure BDA0002214747460000079
是插值函数矩阵，由wm,p(r)组成；
Figure BDA00022147474600000710
是作用函数矩阵，由 GD；uv(rm,p,rn,q)组成。
假设立方体m中有M个点，立方体n中有N个点，那么用矩量法(MOM)直接计算
Figure BDA00022147474600000711
的时间复杂度为O(MN)，用插值公式计算的时间复杂度为O(MK+KK+KN)，其中K为插值点的个数，插值点的数目是有限的，而且很少的插值点就可以达到很高的精度，即K<<min(M,N),所以时间复杂度又可以写成O(MK+KN)。可见 MLGFIM-solver比MOM效率要好得多。
具体的迭代流程如图1所示，本发明设计的一种快速微波成像方法具体实施方法包括以下步骤：
步骤1：探测区域外设置Ninc个发射天线和Nr个接收天线，发射天线发射频率f 的电磁波到有未知散射体的探测区域，接收天线接收对应的测量散射场
Figure BDA0002214747460000081
步骤2：根据公式(1)-(2)计算格林函数矩阵
Figure BDA0002214747460000082
入射场矩阵
Figure BDA0002214747460000083
由平面波的数学表达式得到，根据公式(17)-(18)来构建MLGFIM-solver正问题求解器。
步骤3：根据截断奇异值分解技术定义新的成本函数(7)，并且进行初始化：设定迭代次数n＝0，
Figure BDA0002214747460000084
并且利用波恩近似得到
Figure BDA0002214747460000085
步骤4：更新二次入射电场
Figure BDA0002214747460000086
和非均匀格林函数
Figure BDA0002214747460000087
在共轭梯度迭代法中 CG中嵌入构建好的MLGFIM-solver求解器来求解方程(8)-(9)。
步骤5：对更新后的格林函数矩阵
Figure BDA0002214747460000088
进行奇异值分解，然后通过公式(11)得到确定性电流的主导部分
Figure BDA0002214747460000089
步骤6：计算成本函数(7)的值，判断是否小于阈值，如果满足要求则停止迭代，如果不满足要求，通过求(7)的最小值更新
Figure BDA00022147474600000810
通过(12)得到新的散射强度张量然后转到步骤4继续执行。
本发明可用于快速安检成像，将以往的x光成像换成微波成像，降低了用途成本，并且成像速度快，实验装置结构图如图2所示，实施例的探测区域D为一个边长为2λ， (λ为在自由空间介质中的入射波波长)中心为原点的立方体。用48个发射天线去照射探测区域，这48个发射天线分别位于xoy,yoz,xoz的三个平面圆上，每个圆半径为 3λ并且均匀地分布16个发射天线。在xoy,yoz,xoz平面上的天线极化方式分别为z极化,x极化,y极化。接收天线的位置和个数和入射天线是一样的。为了得到实验仿真数据，对探测区域进行剖分为60×60×60个网格，我们采用偶极耦合法(CDM)正向求解器来得到散射场数据。
实施例1
如图4，实施例1是两个半径为0.5λ的球，两个球的相对介电常数为1.5，它们的圆心坐标分别为(0,0.5λ,0)，(0,0.5λ,0)。在MS-DBIM-MLGFIM的迭代过程中，将成像区域离散为32×32×32个网格，通过18次的迭代，图5为所发明方法的成像结果，图5中的(a),(b),(c)分别为反演结果的z＝0,x＝0,y＝-0.5λ的截面图。可以看出，所发明的反演算法MS-DBIM-MLGFIM能精确的反演出两个球的位置和尺寸，并且能估计它们的相对介电常数大约为1.5。
仿真例2
为了进一步验证所提出方法的有效性，用一个更复杂的设计目标来测试该方法。这个未知由两个球体和U形结构组成，如图6所示。两个球体的球心分别位于(0, -0.5λ,0.4λ),(0,0.5λ,0.4λ)。U型结构的最低点位于(0,0,-0.9λ)，整个结构关于x,y 两个坐标轴对称。将成像区域均匀划分为48×48×48均匀小立方体，每个小立方体的边长为0.03125m，图6物体中的各参数r＝0.35λ,H＝0.5λ,W＝0.4λ,L＝0.1912m,S＝0.2λ。经过16次的迭代，图7为所发明方法的成像结果，图7中的(a),(b),(c)分别为反演结果的x＝0,z＝0.4λ,y＝-0.5λ的截面图。我们可以观察到，U形结构部分的实部介电常数略低于1.5，但散射体的形状和位置仍然得到了很好的重构。
上述两实例仅仅只是例证本发明方法，并非是对于本发明的限制，本发明也并非仅限于上述实例，只要符合本发明方法的要求，均属于本发明方法的保护范围。

Claims (1)

1.一种快速微波成像方法，其特征在于包括以下步骤：
步骤(1)、探测区域外设置Ninc个发射天线和Nr个接收天线，发射天线发射频率f的电磁波到有未知散射体的探测区域，接收天线接收对应的测量散射场
Figure FDA0003200954290000011
步骤(2)、首先将探测区域离散为M个离散方块，并根据已知的电磁波频率f、探测区域位置、发射天线和接收天线的位置和极化方式、测量散射场
Figure FDA0003200954290000012
计算格林函数矩阵
Figure FDA0003200954290000013
和离散后入射场
Figure FDA0003200954290000014
并构建MLGFIM-solve求解器；
格林函数矩阵
Figure FDA0003200954290000015
为离散后的格林函数
Figure FDA0003200954290000016
的积分算子；
Figure FDA0003200954290000017
为离散后的入射场矩阵，其中并矢格林函数
Figure FDA0003200954290000018
表示探测区域内点r′处的点源对其接收天线rs处的作用，是一个3×3的矩阵；标量格林函数g(rs，r′)＝exp(ik0|rs-r′|)/(4π|rs-r′|)，k0为背景介质的波数；I是一个3×3的单位对角矩阵，
Figure FDA0003200954290000019
表示旋度算子，i表示复数的虚部；
MLGFIM-solve求解器是对
Figure FDA00032009542900000110
矩阵向量相乘操作进行加速；其中
Figure FDA00032009542900000111
代表任意一个向量，
Figure FDA00032009542900000112
为离散化的格林函数
Figure FDA00032009542900000113
的积分算子；并矢格林函数
Figure FDA00032009542900000114
表示的是探测区域内点r′对另一点r的作用；
Figure FDA00032009542900000115
是一个3M×3M的矩阵，可由9个M×M的矩阵
Figure FDA00032009542900000116
组成，即
Figure FDA00032009542900000117
当m≠n，(n，m＝1，2，...，M)
Figure FDA00032009542900000118
当m＝n，
Figure FDA00032009542900000119
其中w分别是电磁波的角频率；μ0是背景介质的磁导率；ρ(u-v)表示冲激函数，当u＝v时ρ(u-v)＝1，当u≠v时ρ(u-v)＝0；Rn，m＝|rn-rm|，rn＝(r1：n，r2；n，r3；n)，rm＝(r1：m，r2；m，r3；m)分别是第n，m个离散方块的中心坐标，其中1，2，3分别代表三维坐标下的x，y，z方向上的分量；
步骤(3)、定义新的成本函数
采用了截断奇异值分解(TSVD)技术，对
Figure FDA0003200954290000021
进行奇异值分解：
Figure FDA0003200954290000022
其中
Figure FDA0003200954290000023
表示对
Figure FDA0003200954290000024
奇异值分解得到的左边矩阵的第j列，
Figure FDA0003200954290000025
表示第j个电流基；
Figure FDA0003200954290000026
表示
Figure FDA0003200954290000027
的共轭，将奇异值σj的绝对值按照从大到小排序，取前L个较大奇异值σj对应的电流基构成一个电流子空间，这个电流子空间构成了电流的主要组成部分；上述电流子空间即为感应电流的确定性部分，可以表示为：
Figure FDA0003200954290000028
其中
Figure FDA0003200954290000029
表示
Figure FDA00032009542900000210
的共轭；
Figure FDA00032009542900000211
是二次入射电场；
Figure FDA00032009542900000212
为散射强度张量，为一个3M×3M的对角矩阵，可以表示为
Figure FDA00032009542900000213
当m≠M，2M，3M，则q＝mod(m，M)，mod表示的是取余符号；其他情况下q＝M；Vq，∈r；q分别是第q个离散方块的体积和介电常数；
利用TSVD技术，将归一化处理后的探测区域内的测量散射场
Figure FDA00032009542900000214
和感应电流确定性部分的叠加后得到新的成本函数：
Figure FDA00032009542900000215
其中
Figure FDA00032009542900000216
为每次更新得到的第p个入射场对应的感应电流确定性部分；δ是散射张量系数；
步骤(4)、初始化：设定迭代次数n＝0，
Figure FDA00032009542900000217
步骤(5)、通过公式(8)-(9)分别更新二次入射电场
Figure FDA00032009542900000218
和非均匀格林函数将步骤(2)MLGFIM-solver嵌入到共轭梯度迭代法(CG)中来求解；
Figure FDA00032009542900000219
Figure FDA00032009542900000220
步骤(6)、对更新后的格林函数矩阵
Figure FDA0003200954290000031
根据公式(10)进行奇异值分解，则感应电流相对应的确定性部分
Figure FDA0003200954290000032
通过公式(11)得到；
Figure FDA0003200954290000033
Figure FDA0003200954290000034
步骤(7)、然后将
Figure FDA0003200954290000035
代入到公式(7)得到当前迭代的成本函数f，若f小于阈值，则停止迭代输出当前散射强度张量
Figure FDA0003200954290000036
跳转至步骤(9)，否则跳转至步骤8；
步骤(8)、重置n＝n+1通过求公式(7)的最小值来更新
Figure FDA0003200954290000037
根据公式(12)得到新的散射强度张量
Figure FDA0003200954290000038
然后跳转至步骤5继续迭代：
Figure FDA0003200954290000039
步骤(9)、根据当前更新得到的散射强度张量进行成像。
CN201910911179.9A 2019-09-25 2019-09-25 一种快速微波成像方法 Active CN110687535B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910911179.9A CN110687535B (zh) 2019-09-25 2019-09-25 一种快速微波成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910911179.9A CN110687535B (zh) 2019-09-25 2019-09-25 一种快速微波成像方法

Publications (2)

Publication Number Publication Date
CN110687535A CN110687535A (zh) 2020-01-14
CN110687535B true CN110687535B (zh) 2021-10-08

Family

ID=69110101

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910911179.9A Active CN110687535B (zh) 2019-09-25 2019-09-25 一种快速微波成像方法

Country Status (1)

Country Link
CN (1) CN110687535B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111428407B (zh) * 2020-03-23 2023-07-18 杭州电子科技大学 一种基于深度学习的电磁散射计算方法
CN115184927B (zh) * 2022-07-25 2023-04-07 北京众智信安信息技术研究院 一种微波无损成像目标检测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101581779A (zh) * 2008-05-14 2009-11-18 中国科学院电子学研究所 一种层析合成孔径雷达三维成像原始回波信号生成方法
CN101581780A (zh) * 2008-05-14 2009-11-18 中国科学院电子学研究所 一种用于侧视层析合成孔径雷达的三维聚焦成像方法
CN108182316A (zh) * 2017-12-27 2018-06-19 杭州泛利科技有限公司 一种基于人工智能的电磁仿真方法及其电磁大脑

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101581779A (zh) * 2008-05-14 2009-11-18 中国科学院电子学研究所 一种层析合成孔径雷达三维成像原始回波信号生成方法
CN101581780A (zh) * 2008-05-14 2009-11-18 中国科学院电子学研究所 一种用于侧视层析合成孔径雷达的三维聚焦成像方法
CN108182316A (zh) * 2017-12-27 2018-06-19 杭州泛利科技有限公司 一种基于人工智能的电磁仿真方法及其电磁大脑

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Fast inverse scattering solutions using the distorted Born iterative method and the multilevel fast multipole algorithm;Andrew J. Hesford等;《Acoustical Society of America》;20100831;全文 *
Mert Hidayetoglu等.Parallel Solutions of Inverse Multiple Scattering Problems with Born-type Fast Solvers.《2016 Progress In Electromagnetic Research Symposium (PIERS)》.2016, *
Parallel Solutions of Inverse Multiple Scattering Problems with Born-type Fast Solvers;Mert Hidayetoglu等;《2016 Progress In Electromagnetic Research Symposium (PIERS)》;20160831;全文 *
Subspace-Based Distorted-Born Iterative Method for Solving Inverse Scattering Problems;Xiuzhu Ye等;《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》;20171231;第65卷(第12期);全文 *
基于频率综合技术的新型微波介质检测方法;陈世昌等;《2018年全国微波毫米波会议文集（上集）》;20181231;全文 *

Also Published As

Publication number Publication date
CN110687535A (zh) 2020-01-14

Similar Documents

Publication Publication Date Title
Abubakar et al. Application of the multiplicative regularized Gauss–Newton algorithm for three-dimensional microwave imaging
Joachimowicz et al. Inverse scattering: An iterative numerical method for electromagnetic imaging
Guo et al. Physics embedded deep neural network for solving full-wave inverse scattering problems
Zhou et al. An improved deep learning scheme for solving 2-D and 3-D inverse scattering problems
Chen et al. Microwave imaging of buried inhomogeneous objects using parallel genetic algorithm combined with FDTD method
CN111610374B (zh) 一种基于卷积神经网络的散射场相位恢复方法
Rekanos et al. Inverse scattering using the finite-element method and a nonlinear optimization technique
CN110687535B (zh) 一种快速微波成像方法
Carpio et al. Bayesian approach to inverse scattering with topological priors
Li et al. A three-dimensional model-based inversion algorithm using radial basis functions for microwave data
He et al. Microwave imaging of 3D dielectric-magnetic penetrable objects based on integral equation method
Chai et al. Compressive sensing for monostatic scattering from 3-D NURBS geometries
Brown et al. On the use of electromagnetic inversion for near-field antenna measurements: A review
Zhang et al. Image reconstruction of electrical capacitance tomography based on an efficient sparse Bayesian learning algorithm
Fedeli Microwave tomography with LSTM-based processing of the scattered field
Zdun et al. Fast MPI reconstruction with non-smooth priors by stochastic optimization and data-driven splitting
Guo et al. A nonlinear model compression scheme based on variational autoencoder for microwave data inversion
Narendra et al. Phaseless Gauss-Newton inversion for microwave imaging
Lobel et al. Microwave imaging: Reconstructions from experimental data using conjugate gradient and enhancement by edge‐preserving regularization
Kılıç et al. Solution of 3D inverse scattering problems by combined inverse equivalent current and finite element methods
Franza et al. SICS: A sensor interaction compensation scheme for microwave imaging
Martin et al. Experimental microwave imaging system calibration via cycle-GAN
Caorsi et al. Equivalent current density reconstruction for microwave imaging purposes
Ilvonen et al. Numerical assessment of induced ELF currents in the human head due to the battery current of a digital mobile phone
Wang et al. Underwater active electrosense: A scattering formulation and its application

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