CN117607971A - 一种成像方法、成像设备、装置及存储介质 - Google Patents

一种成像方法、成像设备、装置及存储介质 Download PDF

Info

Publication number
CN117607971A
CN117607971A CN202311390015.9A CN202311390015A CN117607971A CN 117607971 A CN117607971 A CN 117607971A CN 202311390015 A CN202311390015 A CN 202311390015A CN 117607971 A CN117607971 A CN 117607971A
Authority
CN
China
Prior art keywords
imaging
receiving array
receiving
array
transmitting
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.)
Pending
Application number
CN202311390015.9A
Other languages
English (en)
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.)
Suzhou Weimo Electronic Information Technology Co ltd
Original Assignee
Suzhou Weimo Electronic Information 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 Suzhou Weimo Electronic Information Technology Co ltd filed Critical Suzhou Weimo Electronic Information Technology Co ltd
Priority to CN202311390015.9A priority Critical patent/CN117607971A/zh
Publication of CN117607971A publication Critical patent/CN117607971A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/12Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V8/00Prospecting or detecting by optical means
    • G01V8/005Prospecting or detecting by optical means operating with millimetre waves, e.g. measuring the black losey radiation

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Electromagnetism (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及光学成像、微波成像、声成像以及磁成像等技术领域,具体涉及一种成像方法、成像设备、装置及存储介质。本发明方法基于近场成像技术,通过设置位置固定的多个发射单元,采用合成孔径接收阵列或实孔径接收阵列接收目标回波信号,对每个发射单元对应的回波信号进行近场成像和叠加处理,最终形成聚合成像结果。本发明方法具有发射单元、接收单元数量少、硬件成本低、成像速度快等优点,可广泛应用于光学成像、微波成像、声成像、磁成像等领域。

Description

一种成像方法、成像设备、装置及存储介质
技术领域
本发明涉及光学成像、微波成像、声成像、磁成像等技术领域,具体涉及一种成像方法、成像设备、装置及存储介质。
背景技术
现有的微波人体安检成像设备,大部分都采用主动式微波全息成像技术,该技术虽然具有较好的成像效果,但设备复杂,需要设置大量间距很小的密集的发射和接收单元,甚至需要采用极为复杂的开关矩阵对发射、接收单元进行控制,从而造成设备复杂、硬件成本较高。
在主动式微波成像领域,目前应用较为广泛的是微波全息成像技术,如“一种近场微波全息成像方法、装置、电子设备及存储介质(中国专利申请号:202210388580.0)”、“一种用于人体安全检查的毫米波全息成像设备(中国专利申请号:2013103568638)”,相关技术都是基于微波全息成像技术,该技术理论上要求发射/接收天线单元的间距不大于1/4波长,并且发射/接收单元按1:1配置,且单元数量庞大,成本极其高昂。
因此,针对微波人体安检成像设备及类似设备,开发一种效果相近、成本更低的成像方法具有重大的应用价值。
发明内容
本发明所要解决的问题是:在保证成像效果基本相同的情况下,大幅降低成像系统的硬件成本。
传统微波主动全息成像技术中,发射单元数量和接收单元数量是1:1设置的,并且为了防止成像混叠,需要设置发射/接收单元的间距不大于1/4波长,此外,由于不同的发射——接收对在时间上错开依次工作,因而还需要设置复杂的开关矩阵。若能够将原来按接收单元数量1:1设置的发射单元数量降低为数量有限的几个单元,则可以将发射单元的数量降低90%以上,若再将接收单元的间距增加到1/2波长,则接收单元数量可降低50%,则显然可以极大的压缩硬件成本,提高产品的市场竞争力。基于上述设想,经过技术可行性验证后,本发明提供了一种成像方法、成像设备、装置及存储介质:
第一方面,本发明提供了一种成像方法,包括:
步骤一,设置至少2个工作时空间位置相对固定的发射单元;
步骤二,设置具有多个接收单元的接收阵列;
步骤三,设置其中一个发射单元发射成像探测信号;
步骤四,接收阵列接收目标回波信号;
步骤五,采用近场反演方法,对接收阵列接收的目标回波信号进行近场成像;
步骤六,对成像结果进行叠加求和处理;
步骤七,变更发射单元,重复步骤3过程,直至完成全部发射单元的发射和接收成像处理。
第二方面,本发明提供了一种成像设备,包括:
所述成像设备包括壳体、发射单元、接收阵列、信号处理模块、显示装置,其中,发射单元和接收阵列沿待探测区域设置,发射单元用来对待探测区域发射探测信号,接收阵列用来接收来自待探测区域的目标反射信号,信号处理模块用来对接收阵列收到的目标反射信号采用发明1中方法进行近场成像,显示装置用来显示反演成像结果;所述成像设备设置有多个发射单元,发射单元的工作位置固定不变,不随着接收阵列的位置改变而改变。
第三方面,本发明提供了一种装置,包括:
所述装置包含处理器、可编程器件、存储器和通信接口;其中,存储器、可编程器件存储一个或多个程序,该一个或多个程序包含可执行指令;该装置运行时,处理器、可编程器件执行存储的可执行指令,以使该装置执行发明1中的方法。
第四方面,本发明提供了一种计算机可读存储介质,包括:
所述计算机可读存储介质存储有可执行指令,当所述可执行指令在计算上运行时,使得计算机执行发明1中的方法。
本发明方案中,将原来和接收单元1:1设置的发射单元减少为有限的几个,同时去掉了复杂的开关矩阵,并且将接收单元的数量减少了50%,显然能够极大的降低硬件成本。在成像处理时采用近场成像技术,保证了基本相同的成像效果。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对本发明实施例描述中所需要使用的附图作简要介绍,显而易见地,以下附图仅仅是本发明中记载的一些实施例,对于本领域技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明方法的流程图。
图2为本发明方法的坐标系示意图。
图3为本发明方法的成像效果验证结果,其中:左图为目标模型,右图为本发明方法成像结果。
图4为本发明成像设备主要组成框图。
图5为本发明平面或曲面形式实孔径阵列收发单元分布示意图,其中:左图为平面阵列,右图为曲面阵列。
图6为本发明一维直线阵直线扫描的平面合成孔径示意图,其中:左图为竖直一维直线阵水平扫描,右图为水平一维直线阵垂直扫描。
图7为本发明一维线阵的曲面合成孔径示意图其中:左图为一维曲线阵直线扫描,右图为一维直线阵曲线扫描。
图8为本发明装置的主要组成框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合具体实施例及相应的附图对本发明的技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例,本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。
同时,应理解,本发明的保护范围并不局限于下述特定的具体实施方案;还应当理解,本发明实施例中使用的术语是为了描述特定的具体实施方案,而不是为了限制本发明的保护范围。
实施例1:一种成像方法及成像验证结果
附图1给出成像方法的一般流程,其具体步骤如下:
步骤(1),设置至少2个工作时空间位置相对固定的发射单元:
根据所要探测目标的位置、尺寸等因素,合理设置若干数量的发射单元,通过不同位置的发射单元的照射,从而获得无明显盲区的目标全景像;
设置不同位置的多个发射单元,其目的是“照亮”目标,因而对发射单元的间距没有确定要求,只要能够“照亮”目标的全部即可,可以通过合理布局,尽量用最少的发射单元;
也可以采用1个发射单元,通过运动机构在几个相对固定的位置“照亮”目标即可。
步骤(2),设置具有多个接收单元的接收阵列:
接收阵列由大量接收单元组成,对于常规阵列,接收单元的间距不大于1/2波长,对于稀疏阵,则接收单元的间距可大于1/2波长;
根据成像要求,可以选择合成孔径阵列,也可以选择实孔径阵列。
步骤(3),设置其中一个发射单元发射成像探测信号:
发射单元根据需要,选择合适的发射探测信号,通常可选择线性调频信号、频率步进信号等。
步骤(4),接收阵列接收目标回波信号:
对应着每个不同的发射单元,接收阵接收所对应的目标回波信号,对回波信号进行放大、变频、滤波,进行A/D变换后将数字回波信号送信号处理系统;
步骤(5),采用近场反演方法,对接收阵列接收的目标回波信号进行近场成像:
近场成像方法,包括且不限于基于快速傅里叶变换(Fast Fourier Transform,FFT)的口径场快速反演方法、基于等效磁流法的快速反演方法、时间反演法(TimeReversal,TR)、后向投影法(Back Projection,BP)、虚拟透镜成像方法(Virtual LensImaging,VLI)、近场波束合成法;
附图2给出了本发明成像方法的坐标系示意图,其中σ为待成像的目标,位于xy平面上;接收阵列位于XY平面上,目标所在平面到接收阵列所在平面的距离为R,某个发射单元对目标进行照射,接收阵列接收目标的回波信号;
成像时可采用基于FFT的口径场快速反演方法进行成像,方法如下:
其中,σ(x,y,z)为待求解的目标的像,(x,y,z)为目标的坐标,IFT3表示三维快速傅里叶逆变换,FT2表示二维快速傅里叶变换,E(X,Y,ω)表示接收阵列接收的目标反射信号,(X,Y)为接收单元的坐标,ω为角频率,R为目标到接收阵列的距离,为波矢量/>沿坐标z方向的分量,kx为波矢量/>沿坐标x方向的分量,ky为波矢量/>沿坐标z方向的分量,/>为波数,e为欧拉常数,j为虚数单位。
步骤(6),对成像结果进行叠加求和处理:
针对每个不同发射单元的成像结果,其直接累加处理时会出现相互干涉现象,可能导致成像效果变差,为了防止该现象并改善成像效果,采用取绝对值后进行叠加求和:
σn(x,y,z)=σn-1(x,y,z)+|σ(x,y,z)|;
其中,σn(x,y,z)为当前成像结果,σn-1(x,y,z)为前期成像结果,符号“||”表示取绝对值运算。
步骤(7),变更发射单元,重复步骤3过程,直至完成全部发射单元的发射和接收成像处理。
成像验证:对本实施例方法进行验证:目标为“F”形理想点目标集合,点目标之间的间距为15mm,点目标集合位于阵列法线方向,探测信号频率为30GHz,接收阵列规模为256×256,单元间距为1/2波长,发射单元有5个,分别位于接收阵列的中心以及4个角处。点目标集合距离阵列中心0.4m。附图3左图为“F”形理想点目标集合的模型,右图为本发明方法的成像结果,由图中能够明显分辨出其中的点目标。
需要说明的是,上述步骤仅给出了本发明方法的一般实现流程,在实际应用中可能会调整顺序或增加/减少一些步骤。
实施例2:一种成像设备
附图4给出一种成像设备的主要组成框图,所述成像设备包括壳体、发射单元、接收阵列、信号处理模块、显示装置,其中,发射单元和接收阵列沿待探测区域设置,发射单元用来对待探测区域发射探测信号,接收阵列用来接收来自待探测区域的目标反射信号,信号处理模块用来对接收阵列收到的目标反射信号采用实施例1中方法进行近场成像,显示装置用来显示成像结果。
所述成像设备设置有多个发射单元,发射单元的工作位置固定不变,不随着接收阵列的位置改变而改变。
本实施例成像设备的接收阵列可以采用实孔径阵列,附图5给出了平面形式和曲面形式的实孔径收发单元分布示例,其中发射单元也可设置在接收阵列之外。
本实施例成像设备也可以采用一维线阵扫描形成合成孔径阵列。附图6给出了一维直线阵直线扫描的平面合成孔径示意图,图7给出了一维直线阵曲线扫描或一维曲线阵直线扫描的曲面合成孔径示意图。
需要注意的是,本实施例成像设备的接收阵列的单元间距不大于1/2波长。
实施例4:一种装置
所述装置组成结构见附图8,包含处理器、可编程器件、存储器和通信接口;其中,存储器、可编程器件存储一个或多个程序,该一个或多个程序包含可执行指令;该装置运行时,处理器、可编程器件执行存储的可执行指令,以使该装置执行实施例1中任一项方法或其混合方法。
实施例5:一种存储介质
所述存储介质存储有可执行指令,当所述可执行指令在计算上运行时,使得计算机执行实施例1中任一项方法或其混合方法。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到本发明可以借助软件加必要的硬件平台的方式来实现。基于这样的理解,上述实施方案可以以软件产品的形式体现出来,该软件产品可以存储在一个非易失性存储介质(包含光盘,U盘,EPROM,硬盘等)中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者嵌入式设备等)执行本发明实施方案所述的方法。
最后应当说明的是:以上叙述仅用以说明本发明的基本技术方案,而非对其限制;尽管参照上述内容对本发明进行了详细的说明,任何熟悉本领域的技术人员可以对前述的技术方案进行修改完善或适当删减,或者对其中部分技术特征进行等同替换,而这些修改或者替换,都应涵盖在本发明的保护范围之内。

Claims (10)

1.一种成像方法,其特征在于,包括以下步骤:
步骤一,设置至少2个工作时坐标位置相对固定的发射单元;
步骤二,设置具有多个接收单元的接收阵列;
步骤三,设置其中一个发射单元发射成像探测信号;
步骤四,接收阵列接收目标回波信号;
步骤五,采用近场成像方法,对接收阵列接收的目标回波信号进行近场成像;
步骤六,对成像结果进行叠加求和处理;
步骤七,变更发射单元,重复步骤三过程,直至完成全部发射单元的发射和接收成像处理。
2.根据权利要求1所述的方法,其特征在于,步骤五中所述采用近场成像方法,对接收阵列接收的目标回波信号进行近场成像,包括:
所述近场成像方法,包括且不限于基于快速傅里叶变换(Fast Fourier Transform,FFT)的口径场快速反演方法、基于等效磁流法的快速反演方法、时间反演法(TimeReversal,TR)、后向投影法(Back Projection,BP)、虚拟透镜成像方法(Virtual LensImaging,VLI)、近场波束合成法。
3.根据权利要求2所述的方法,其特征在于,包括:
所述基于FFT的口径场快速反演方法,包括:
其中,σ(x,y,z)为待求解的目标的像,(x,y,z)为目标的坐标,IFT3表示三维快速傅里叶逆变换,FT2表示二维快速傅里叶变换,E(X,Y,ω)表示接收阵列接收的目标反射信号,(X,Y)为接收单元的坐标,ω为角频率,R为目标到接收阵列的距离,为波矢量沿坐标z方向的分量,kx为波矢量/>沿坐标x方向的分量,ky为波矢量/>沿坐标z方向的分量,/>为波数,e为欧拉常数,j为虚数单位。
4.根据权利要求1所述的方法,其特征在于,步骤六中所述对成像结果进行叠加求和处理,包括:
将本次成像结果与前期成像结果进行叠加求和处理:
σn(x,y,z)=σn-1(x,y,z)+|σ(x,y,z)|;
其中,σn(x,y,z)为当前成像结果,σn-1(x,y,z)为前期成像结果,符号“||”表示取绝对值运算。
5.一种成像设备,包括:
所述成像设备包括壳体、发射单元、接收阵列、信号处理模块、显示装置,发射单元和接收阵列沿待探测区域设置,发射单元用来对待探测区域发射探测信号,接收阵列用来接收来自待探测区域的目标反射信号,显示装置用来显示成像结果,其特征在于:设置有多个发射单元,发射单元的工作位置固定不变,不随着接收阵列的位置改变而改变;信号处理模块采用权利要求1~4任一项所述方法对接收阵列收到的目标反射信号进行近场成像。
6.根据权利要求5所述的成像设备,其特征在于,所述接收阵列的单元间距不大于1/2波长。
7.根据权利要求5所述的成像设备,其特征在于,所述接收阵列为一维线阵,通过一维线阵的直线或曲线扫描,形成合成孔径平面或曲面二维接收阵列。
8.根据权利要求5所述的成像设备,其特征在于,所述接收阵列为平面或曲面实孔径二维阵列。
9.一种装置,其特征在于,包括:
所述装置包含处理器、可编程器件、存储器和通信接口;其中,存储器、可编程器件存储一个或多个程序,该一个或多个程序包含可执行指令;该装置运行时,处理器、可编程器件执行存储的可执行指令,以使该装置执行权利要求1~4中任一项方法或其混合方法。
10.一种存储介质,其特征在于,包括:
所述存储介质存储有可执行指令,当所述可执行指令在计算上运行时,使得计算机执行权利要求1~4中任一项方法或其混合方法。
CN202311390015.9A 2023-10-25 2023-10-25 一种成像方法、成像设备、装置及存储介质 Pending CN117607971A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311390015.9A CN117607971A (zh) 2023-10-25 2023-10-25 一种成像方法、成像设备、装置及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311390015.9A CN117607971A (zh) 2023-10-25 2023-10-25 一种成像方法、成像设备、装置及存储介质

Publications (1)

Publication Number Publication Date
CN117607971A true CN117607971A (zh) 2024-02-27

Family

ID=89952320

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311390015.9A Pending CN117607971A (zh) 2023-10-25 2023-10-25 一种成像方法、成像设备、装置及存储介质

Country Status (1)

Country Link
CN (1) CN117607971A (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105877747A (zh) * 2016-03-30 2016-08-24 厦门大学 基于快速体积分方程和磁共振的人体电磁特性反演方法
CN108828592A (zh) * 2018-04-25 2018-11-16 中国科学院电子学研究所 基于mimo矩形平面阵列的方位向成像方法
US20200225337A1 (en) * 2018-03-22 2020-07-16 Panasonic Intellectual Property Management Co., Ltd. Radar apparatus, position estimation apparatus, and position estimation method
CN111505637A (zh) * 2020-04-29 2020-08-07 中国科学院国家空间科学中心 一种基于二单元扫描干涉仪的自定标近场成像方法及其系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105877747A (zh) * 2016-03-30 2016-08-24 厦门大学 基于快速体积分方程和磁共振的人体电磁特性反演方法
US20200225337A1 (en) * 2018-03-22 2020-07-16 Panasonic Intellectual Property Management Co., Ltd. Radar apparatus, position estimation apparatus, and position estimation method
CN108828592A (zh) * 2018-04-25 2018-11-16 中国科学院电子学研究所 基于mimo矩形平面阵列的方位向成像方法
CN111505637A (zh) * 2020-04-29 2020-08-07 中国科学院国家空间科学中心 一种基于二单元扫描干涉仪的自定标近场成像方法及其系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
周勇,郁佳蜻,陈挺: "基于CGFFT和等效磁流法的天线近场测量方法研究", 中国测试, vol. 46, no. 1, 31 January 2020 (2020-01-31), pages 18 - 23 *
彭虎编著: "超声成像算法导论", 31 October 2008, 中国科学技术大学出版社, pages: 114 *
李志;金添;陈波;周智敏;: "基于协同阵的超宽带MIMO成像阵列优化设计方法", 信号处理, no. 07, 25 July 2012 (2012-07-25) *
王祥;郭良辉;: "球坐标系密度界面反演方法及在华南大陆的应用", 物探与化探, no. 05, 29 September 2020 (2020-09-29) *

Similar Documents

Publication Publication Date Title
US11215730B2 (en) Close range microwave imaging method and system
US11561300B2 (en) Close-range microwave imaging method and system
US20150061926A1 (en) Target detection utilizing image array comparison
CN111650585B (zh) 近场毫米波稀疏mimo扫描阵列全聚焦成像方法和装置
WO2018196254A1 (zh) 一种微波成像系统的幅相校正方法及系统
CN110609282B (zh) 基于后向投影的太赫兹孔径编码三维成像方法及装置
CN106680796B (zh) 基于频率干涉的平面全息阵列目标三维表面重构方法
CN109884627B (zh) 任意线阵构型的近程毫米波快速三维成像方法
CN107390216A (zh) 基于波数域相干因子的高速超分辨率驻点扫描成像方法
CN113158485B (zh) 近场条件下电大尺寸目标的电磁散射仿真方法
CN107024693B (zh) 一种单发射体制的雷达关联成像方法
CN108008389A (zh) 一种基于gpu的快速频域后向投影三维成像方法
CN113156429B (zh) 基于毫米波的成像方法、装置、设备和存储介质
US11199612B2 (en) Direct wave suppression method and system for microwave imaging system
CN117607971A (zh) 一种成像方法、成像设备、装置及存储介质
CN112649806A (zh) 一种mimo雷达近场三维成像方法
CN111580112A (zh) 基于平面波的水下声呐传感器阵列成像方法
CN112099102B (zh) 一种天线阵列机械扫描式全息图像重建方法
CN115657017A (zh) 电大尺寸目标快速成像的增量长度绕射理论射线成像方法
CN114200448A (zh) 综合孔径辐射计波数域近场成像方法及设备
Berizzi et al. A new fast method for the reconstruction of 2-D microwave images of rotating objects
US20020138004A1 (en) Ultrasound imaging method using redundant synthetic aperture concepts
CN112666625A (zh) 用于毫米波安检的快速成像设备及其方法
Prechtel et al. Short-range MIMO radar system considerations
Xu et al. Ultra‐wideband imaging with an improved backward projection algorithm for far‐field applications

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