Method and apparatus for reconstructing dielectric image using electromagnetic waves

Patent metadata

Field	Value
Publication number	US-9519054-B2
Application number	US-201414254439-A
Country	US
Kind code	B2
Filing date	Apr 16, 2014
Priority date	Aug 9, 2013
Publication date	Dec 13, 2016
Grant date	Dec 13, 2016

Abstract

A method for reconstructing a dielectric image using electromagnetic waves, comprising: acquiring a measurement value of the electromagnetic waves; generating a matching system matrixes between the meshes; generating a smoothed dyadic Green's function matrix; generating an electromagnetic wave calculation value; calculating a misfit error between the acquired electromagnetic wave measurement value and the generated electromagnetic wave calculation value, and checking whether a change of the calculated misfit error satisfies a predetermined optimization determination condition; updating dielectric parameters at the meshes; and outputting a reconstructed dielectric image in the image reconstruction region.

First claim

What is claimed is:

1. A method for reconstructing an image using electromagnetic waves, comprising:
acquiring a measurement value of the electromagnetic waves that are radiated by a transmitting antenna and propagate through an image reconstruction region;
setting a plurality of predetermined meshes in the image reconstruction region;
generating a matching system matrixes using basis functions related to the plurality of predetermined meshes;
generating a smoothed Green's function matrix using the generated matching system matrix;
generating an electromagnetic wave calculation value through performing of a forward electromagnetic analysis with respect to the image reconstruction region on the basis of the generated smoothed dyadic Green's function matrix;
calculating a misfit between the acquired electromagnetic wave measurement value and the generated electromagnetic wave calculation value, and checking whether a the calculated misfit error satisfies a predetermined optimization determination condition;
updating the dielectric parameters at the meshes if the predetermined optimization determination condition is not satisfied; and
outputting a reconstructed image in the image reconstruction region when the predetermined optimization determination condition is satisfied.

2. The method of claim 1 , wherein the plurality of meshes include a fine mesh, a coarse mesh having a relatively larger spacing than the fine mesh, and an image reconstruction mesh having a relatively larger spacing than the coarse mesh.

3. The method of claim 2 , wherein the plurality of meshes are in any one of rectangular grids, or hexagonal grids.

4. The method of claim 2 , wherein the basis function corresponding to the fine mesh and the coarse mesh and the basis function corresponding to the coarse mesh and the image reconstruction mesh are in a sinc-type function, spline function type or wavelet type.

5. The method of claim 2 , wherein the matrix B c of matching system of the basis function, allocated at the fine mesh and centered at the coarse mesh nodes is generated as in an equation,
B c =[b c,1 ,b c,2 ,b c,3 , . . . ,b c,n , . . . ,b c,N c ]
where, b c,n denotes the vectors of basis function allocated to the fine mesh and centered at the n-th node of coarse mesh, and N c denotes the total number of coarse mesh nodes.

6. The method of claim 5 , wherein the matching system matrix satisfies an equation,
R c ≡B c ·{tilde over (B)} c T =[r c,1 ,r c,2 ,r c,3 , . . . ,r c,n , . . . ,r c,N f ]
where, R c denotes a spatial resolution matrix, {tilde over (B)} c T : denotes a transpose matrix, which is dual to B c , r c,n denotes a vector of a point spread function corresponding to the n-th fine mesh node, and N f denotes the total number of fine mesh nodes.

7. The method of claim 6 , wherein R c has property of a spatial filter for smoothing an image reconstruction value.

8. The method of claim 6 , wherein the matching system matrix simultaneously satisfies two conditions:
Condition 1: The effective width of the basis function is equal to or close to the grid spacing of the coarse mesh
Condition 2: All rows and columns of the spatial resolution matrix correspond to functions having the same shape and the same peak values but centered at different nodes of the fine mesh.

9. The method of claim 5 , wherein the matching system matrix B r of the basis functions allocated to the fine mesh and centered at the image reconstruction mesh nodes is generated as in an equation,
B r =[b r,1 ,b r,2 ,b r,3 , . . . ,b r,n , . . . ,b r,N r ]
where, b r,n denotes the basis function allocated to the fine mesh and centered at the n-th node of image reconstruction mesh, and N r denotes the total number of image reconstruction mesh nodes.

10. The method of claim 9 , wherein the generating the matching system matrix comprises generating a spatial resolution matrix using multiplication of the matching system matrix and a transpose matrix, which is dual to the matching system matrix.

11. The method of claim 2 , wherein the generating the smoothed Green's function matrix comprises generating the smoothed Green's function matrix by multiplying a Green's function matrix defined at the fine mesh by the matching system matrix and a transpose matrix, which is dual to the matching system matrix.

12. The method of claim 11 , wherein the generating the smoothed Green's function matrixes comprises:
calculating the dyadic Green's function matrix defined at all the fine mesh nodes in the image reconstruction; and
calculating the dyadic Green's function matrix for calculating scattered waves incident on receiving antennas.

13. The method of claim 2 , wherein the performing the forward electromagnetic analysis comprises analyzing electromagnetic wave scattering with respect to the coarse mesh and calculating the electromagnetic waves that can be acquired in positions of the receiving antennas.

14. The method of claim 1 , wherein the checking comprises determining that the predetermined optimization condition is satisfied if a difference value between a current misfit error and a previously calculated misfit error is smaller than a predetermined allowable difference value or exceeds a predetermined number of calculations.

15. The method of claim 2 , wherein the updating the dielectric parameters comprises:
generating a regularized Jacobian matrix using the generated matching system matrix with respect to the image reconstruction mesh, and the smoothed dyadic Green's function matrix; and
updating the dielectric parameters at the coarse mesh using the generated regularized Jacobian matrix.

16. The method of claim 15 , wherein the regularized Jacobian matrix J r is generated as in an equation,

J
r

⁡

(

Tx
/
Rx

)

=

1

V
s

⁡

(

Tx
/
Rx

)

·

∑

i
=
1

3

⁢

(

G

signal
,
c

⁡

(
Rx
)

·
A

Assignees

• Electronics & Telecommunications Res Inst

Inventors

• Kim Bo Ra
• Simonov Nikolai
• Son Seong-Ho
• Jeon Soon Ik
• Kim Hyuk Je
• Lee Jong Moon
• Choi Hyung Do

Classifications

• G01S13/888

  through wall detection · CPC title

• G01S13/88

  Radar or analogous systems specially adapted for specific applications (electromagnetic prospecting or detecting of objects, e.g. near-field detection, G01V3/00) · CPC title

• G01S13/887

  for detection of concealed objects, e.g. contraband or weapons · CPC title

• G01S13/89 Primary

  for mapping or imaging · CPC title

• G06T1/00

  General purpose image data processing · CPC title