CN107942328A - Terahertz aperture coding three-dimensional target scanning imaging method - Google Patents

Terahertz aperture coding three-dimensional target scanning imaging method Download PDF

Info

Publication number
CN107942328A
CN107942328A CN201711127463.4A CN201711127463A CN107942328A CN 107942328 A CN107942328 A CN 107942328A CN 201711127463 A CN201711127463 A CN 201711127463A CN 107942328 A CN107942328 A CN 107942328A
Authority
CN
China
Prior art keywords
mrow
dimensional imaging
directions
lens
terahertz
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.)
Granted
Application number
CN201711127463.4A
Other languages
Chinese (zh)
Other versions
CN107942328B (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201711127463.4A priority Critical patent/CN107942328B/en
Publication of CN107942328A publication Critical patent/CN107942328A/en
Application granted granted Critical
Publication of CN107942328B publication Critical patent/CN107942328B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The two-dimensional imaging technology under the current aperture coding imaging system is not suitable for three-dimensional targets any more, and aiming at the defects in the prior art, the invention provides a terahertz aperture coding three-dimensional target scanning imaging method. The invention divides the whole three-dimensional imaging area where the three-dimensional target is located into three-dimensional imaging units with the same size, calculates and updates the phase modulation factor of the lens according to the position of a single three-dimensional imaging unit aiming at each three-dimensional imaging unit, can realize the complete coverage of the three-dimensional imaging units and the complete scanning of the three-dimensional imaging area. The invention can realize forward-looking high-resolution imaging on a close-range three-dimensional target, effectively improve the imaging range and speed of terahertz aperture coding imaging, improve the imaging efficiency, reduce the requirement on the configuration of a computer memory, and can be applied to the close-range imaging fields of security inspection, anti-terrorism, target detection, identification and the like.

Description

A kind of Terahertz aperture encodes objective scan imaging method
Technical field
It is particularly a kind of based on Terahertz aperture coded imaging the invention belongs to objective blocked scan technical field Objective blocked scan method.
Background technology
With the development of society, radar high-resolution imaging is ensuring national strategy safety and the side of promoting national economic development Face plays an increasingly important role.Optical radar can forword-looking imaging, wavelength is short, and high resolution, image taking speed is fast, but according to Rely in target emanation, it is poor to penetration capacitys such as cigarette, dirt, mist and barriers, easily by such environmental effects.And microwave radar can be actively Detection, penetration capacity is strong, but since microwave frequency is low, wavelength length, angular resolution is low, and due to the limitation of image-forming principle, needs Integration time is imaged, can not realize the high frame frequency of forward sight, high-resolution imaging.Wherein synthetic aperture radar (SAR) and inverse synthesis hole For footpath radar (ISAR) although imaging can obtain the high-resolution in transverse direction by synthetic aperture, the two all relies on radar With the relative motion of target, can not forword-looking imaging, and real array of apertures radar with phased-array radar due to needing array element to be used Quantity is more, complicated, builds high with maintenance cost.Microwave relevance imaging technology can realize forward sight, stare under the conditions of High-resolution imaging.It is used as transmitting signal by the array signal that build time is uncorrelated, space is mutually orthogonal, is pushed away by calculating Drill to obtain the detectable signal of target region, and mesh is obtained by the association process between detectable signal and target echo signal Mark information.But this method needs to construct fairly large aerial array in transmitting terminal, and it is difficult to realize effective real-time wave beam It is directed toward regulation and control.
Terahertz aperture coded imaging uses for reference microwave relevance imaging thought, by array code aperture to terahertz wave beam Modulated in real time substitutes the radar array in microwave relevance imaging, so as to form the radiation field of Space Time two-dimensional random distribution, finally Realize high-resolution, forward sight using echo and radiation field reference signal matrix is detected by way of matrix equation solution and stare into Picture, compensate for the deficiency that synthetic aperture high-resolution imaging relies on target movement, realize more complicated various space wave modulation.Relatively In conventional radar, THz wave has higher frequency and shorter wavelength so that Terahertz radar is capable of providing the absolute band of bigger Width, aperture coding techniques is combined under same apertures antenna conditions, is more also easy to produce multifarious irradiation mode and faster pattern Switch speed, the more various then free degree of irradiation mode is higher, and the target information carried in echo is abundanter, and mesh is carried out using echo The potentiality of absolute altitude resolution imaging are also bigger.Meanwhile system receiving-transmitting chain single-shot list is received, miniaturization is easily realized.The system at present Under objective imaging method rarely have people's research, be imaged primarily directed to two dimension target, two-dimensional target imaging method master Be divided into two kinds, for low resolution small area two dimension target can single fraction irradiation direct imaging, for high-resolution a wide range of two Dimension target can carry out blocked scan joining image-forming.It is based on Terahertz aperture coded imaging system y directions two dimension as shown in Figure 1 Target high resolution scanning is imaged scheme schematic diagram, and capitalization A-D represents Terahertz dual-mode antenna, array code respectively in Fig. 1 Aperture, two-dimensional imaging plane and focal plane.
Array code aperture in the vertical direction is uniform-distribution with P array element, and array code aperture is equal in the horizontal direction Even to be dispersed with Q array element, whole array code aperture includes P × Q array element.In view of the imaging system in horizontally and vertically side Upward operating mode has symmetry, is illustrated by taking vertical direction as an example.In the vertical direction, array code aperture pth The lens phase modulation factor of a array element loading (1) formula controls terahertz wave beam to irradiate n-th of target scan areas:
Wherein, fyFor the focal length of lens, k=2 π fc/ c, fcFor the centre frequency of THz wave, c is the light velocity, ypFor array The ordinate of p-th of array element central point on the vertical direction of code aperture, p=1,2 ... P,For lens phase on array code aperture Ordinate at the phase centre location of position modulation factor.The focal plane of two-dimensional imaging is parallel with imaging plane, and the focal length of lens can With by Terahertz dual-mode antenna and array code aperture, the relative position between array code aperture and focal plane determines.Sweep During retouching, the focal length of lens in lens phase modulation factor immobilizes, can be with by varying phase centre location therein Change the scanning area of terahertz wave beam.As shown in the figure, loading N number of different phase modulation factor respectively, terahertz wave beam can To realize that n times scan on imaging plane, and each scanning area size is identical.It is seamless spliced between adjacent scan areas, Such as the scanning area 1 and 2 in Fig. 1.
But the Overlap-scanning method in two-dimensional imaging is not properly suited for objective imaging.By in Fig. 1 based on too The two dimension target high resolution scanning imaging scheme of hertz aperture coded imaging system is transplanted in objective imaging, in three-dimensional point During block scan, regard whole three-dimensional imaging region right flank as two-dimensional imaging plane, pass through loading and two dimension target high score Distinguish lens phase modulation factor identical in blocked scan imaging scheme, it is possible to achieve to whole three-dimensional imaging region right flank Seamless traverse scanning, but there are drawback for this method.As shown in Fig. 2, the Terahertz aperture that Fig. 2 is Y-direction encodes objective High-resolution blocked scan imaging problem illustrates schematic diagram.In the coded imaging system of Terahertz aperture, radiation field reference signal square Battle array scale is influenced by grid cell quantity, and grid cell quantity is determined by grid cell size and three-dimension object scope. Therefore it is true according to Terahertz aperture coded imaging systems parameter first to reduce the imaging system to calculator memory configuration requirement Its fixed imaging resolution size, is then uniformly split into M by imaging resolution and objective size by three-dimensional imaging region again × N number of three-dimensional imaging unit, wherein N number of three-dimensional imaging unit is split on the y directions in three-dimensional imaging region, in its x direction On be split into M three-dimensional imaging unit, the three-dimensional imaging unit number in whole three-dimensional imaging region is exactly M × N.Terahertz wave beam When irradiating three-dimensional imaging unit n, whole unit can be all irradiated to;But when terahertz wave beam irradiation three-dimensional imaging unit During 1,2 and N, the white space in targeted scans unit is not irradiated to, thus in echo the scanning area target information not It is complete.Similarly, the plane using the left surface in three-dimensional imaging region as two-dimensional imaging, two-dimentional blocked scan scheme equally can not Some three-dimensional imaging units are completely covered.
Itd is proposed similar to two-dimensional imaging planar mesh, in microwave relevance imaging to three-dimensional imaging region first along it is horizontal, Vertical and range direction mesh generation, then constructs reference signal matrix, then the method for being associated imaging, the same party Method is also applied for Terahertz aperture coded imaging.But this method be directed to microwave section low resolution imaging, mesh generation compared with To be sparse, reference signal matrix size is relatively small.And relatively large volume of objective high-resolution imaging is directed to, such as pacify Human body in inspection imaging, the reference signal matrix size of required construction expand at double, on the one hand the Terahertz in encoded aperture Wave beam is difficult to the whole human body of single fraction irradiation, on the other hand in large scale for the radiation field R-matrix of whole human body target, into As process is too high to calculator memory requirement.Therefore, there is an urgent need for a kind of new objective based on aperture coded imaging system into Image space method.
The content of the invention
Two-dimensional imaging technique under current aperture coded imaging system is no longer desirable for objective, in addition can be used at present 3 Dimension Image Technique imaging resolution it is low, areas imaging is narrow, high to system for computer request memory.Deposited for the prior art Drawbacks described above, the present invention provides a kind of Terahertz aperture coding objective scan imaging method.The present invention is by three-dimensional mesh Whole three-dimensional imaging region where mark is uniformly split into the identical three-dimensional imaging unit of size, respectively for each three-dimensional imaging Unit is scanned imaging, and the imaging results of all three-dimensional imaging units then are spliced to form whole three-dimensional imaging result.This Invention can realize forward sight high-resolution imaging to closely objective, can be applied to safety check and anti-terrorism, target detection and identification etc. Low coverage imaging field.The present invention can improve three-dimensional imaging resolution ratio, expand three-dimensional imaging scope, reduce the memory to computer Configuration requirement, improves computer operational efficiency.
To achieve the above object, the technical scheme is that:
A kind of Terahertz aperture encodes objective scan imaging method, comprises the following steps:
The first step, determine Terahertz aperture coded imaging systems parameter and carry out subdivision to three-dimensional imaging region;
Known Terahertz aperture coded imaging systems parameter is as follows:The vertical height in the y directions in array code aperture is h, The lateral length v in the x directions in array code aperture;The distance between Terahertz dual-mode antenna and array code aperture are a, array The distance of code aperture and three-dimensional imaging region left surface is b, and the thickness in three-dimensional imaging region is c.
According to the imaging resolution size of Terahertz aperture coded imaging systems and objective size by three-dimensional imaging Region is uniformly split into M × N number of three-dimensional imaging unit, wherein being split into N number of three-dimensional imaging on the y directions in three-dimensional imaging region Unit, is split into M three-dimensional imaging unit, the three-dimensional imaging unit number in whole three-dimensional imaging region is exactly M on its x direction ×N.The step is identical with subdivision method in background technology.
Second step, using the right flank in three-dimensional imaging region as two-dimensional imaging plane carry out the two-dimentional mesh of Terahertz aperture coding Blocked scan imaging is marked, its y directions (vertical side is tentatively solved respectively for each three-dimensional imaging unit on three-dimensional imaging region To) lens phase modulation factor and x directions (horizontal direction) lens phase modulation factor, including the y directions numeral focal length of lens fy,mnWith y directions phase centre locationAnd x directions numeral focal length of lens fx,mnWith x directions phase centre locationIts Middle m and n represent that m rows n-th arrange a three-dimensional imaging unit, m=1,2 ... M, n=1,2 ... N on three-dimensional imaging region.
3rd step, judge whether the corresponding lens phase modulation factor of each three-dimensional imaging unit needs in three-dimensional imaging region Correct;
4th step, to judging to need in the 3rd step to carry out the modified three-dimensional imaging unit of lens phase modulation factor to carry out The amendment of y directions (vertical direction) lens phase modulation factor solves and x directions (horizontal direction) lens phase modulation factor is repaiied It is positive to solve, determine the digital focal length of lens f in its revised y directiony,mnWith the phase centre location in revised y directionsAnd the digital focal length of lens f ' in revised x directionsx,mnWith the phase centre location in revised x directions
5th step, lens phase modulation factor comprehensive modification
To make three-dimensional imaging unit all completely be covered by terahertz wave beam in x directions and y directions, more revised x side To digital focal length of lens f 'x,mnWith the digital focal length of lens f ' in revised y directionsy,mn, select wherein larger conduct numeral The focal length of lens:
F '=max (f 'x,mn,f′y,mn) (14)
When Terahertz aperture coded imaging systems carry out blocked scan to three-dimensional imaging region, Terahertz aperture coded imaging The pth row in array code aperture in system, the lens phase modulation factor F of q row array element loading (15) formulaspq,mnTo control too M rows n-th arrange a three-dimensional imaging unit on hertz beam three-dimensional imaging region:
Wherein:P=1,2 ... P, q=1,2 ... Q, k=2 π fc/ c, fcFor the centre frequency of THz wave, c is the light velocity.Array P array element is uniform-distribution with code aperture in the vertical direction, that is, y directions, array code aperture is x side in the horizontal direction Q array element is uniform-distribution with upwards, whole array code aperture includes P × Q array element.(xpq,ypq) for array code aperture the P rows, the coordinate position of q row array element central points, p=1,2 ... P.
In second step of the present invention, a three-dimensional is arranged for each three-dimensional imaging unit such as m rows n-th on three-dimensional imaging region Imaging unit, tentatively solves its y directions (vertical direction) lens phase modulation factor, including y directions numeral focal length of lens fy,mn With y directions phase centre locationMethod is as follows:
The 2.1 preliminary digital focal length of lens f for solving y directionsy,mn
The two-dimentional mesh of Terahertz aperture coding in y directions is carried out using the right flank in three-dimensional imaging region as two-dimensional imaging plane Blocked scan imaging is marked, determines that the terahertz wave beam that the Terahertz dual-mode antenna in the coded imaging systems of Terahertz aperture is sent exists The coverage S of each three-dimensional imaging unit right flankh
According to the focal plane under each three-dimensional imaging unit right flank of triangle similarity relation calculating and two-dimensional imaging system (i.e. The right flank in three-dimensional imaging region carries out corresponding focal plane when Terahertz aperture encodes two-dimensional imaging as two-dimensional imaging plane) Between distance d:
Then, calculated according to lens imaging formula and carry out y using the right flank in three-dimensional imaging region as two-dimensional imaging plane During the Terahertz aperture coding two dimension target blocked scan imaging in direction, m rows n-th arrange its y direction of three-dimensional imaging unit and need The digital focal length of lens f to be loadedy,mn
The 2.2 preliminary phase centre locations for solving y directions
The two-dimentional mesh of Terahertz aperture coding in y directions is carried out using the right flank in three-dimensional imaging region as two-dimensional imaging plane When marking blocked scan imaging, to realize seamless spliced, the phase in two dimension target blocked scan imaging process between each scanning area The lens phase modulation factor phase centre location stepping length of loading is needed between adjacent three-dimensional imaging unitFor:
According to stepping lengthThree-dimensional imaging region m rows n-th can be calculated arrange its y direction of three-dimensional imaging unit to add The phase centre location of load
The imaging system of Terahertz aperture coded imaging systems is in x directions (i.e. horizontal direction) and y directions (i.e. vertical side To) on operating mode be it is identical, for each three-dimensional imaging unit on three-dimensional imaging region can use with 2.1 and 2.2 identical methods are completed tentatively to solve its x directions (horizontal direction) lens phase modulation factor, including the numeral in x directions is saturating Mirror focal length fx,mnWith the phase centre location in x directions
In 3rd step of the invention, the corresponding lens phase modulation factor of each three-dimensional imaging unit in three-dimensional imaging region is judged Whether need modified method as follows:
Before judging whether the corresponding lens phase modulation factor of each three-dimensional imaging unit needs amendment, one is first determined Gray area:If sending parallel beam along the z-axis direction from whole array code aperture, gray area is exactly covering for the parallel beam Lid scope, the vertical height in gray area y directions are the vertical height h, the x of gray area in the y directions in array code aperture The lateral length in direction is the lateral length in the x directions in array code aperture, and the length in gray area z directions is array code Aperture is to the distance between three-dimensional imaging region right flank (the i.e. distance b of array code aperture and three-dimensional imaging region left surface The sum of with the thickness c in three-dimensional imaging region);
When three-dimensional imaging unit is completely in gray area, array code aperture needs the lens phase modulation loaded The factor need not be corrected, and the lens phase modulation factor that array code aperture needs load is identical with two-dimensional target imaging, Jiao Ping Face remains as the perpendicular parallel with three-dimensional imaging region.
When there is the three-dimensional imaging unit for being not entirely in gray area, the lens under loading two-dimensional scan imaging system Phase modulation factor cannot completely cover three-dimensional imaging region, it is necessary to correct lens phase modulation factor so that terahertz wave beam is shone It is mapped to whole three-dimensional imaging region.
In 4th step of the invention, for judge to need in the 3rd step to carry out lens phase modulation factor it is modified it is three-dimensional into As unit (there is the three-dimensional imaging unit for being not entirely in gray area), as m rows n-th arrange a three on three-dimensional imaging region Dimension imaging unit is not entirely in gray area, then the lens phase modulation factor of the three-dimensional imaging unit is modified, its In the amendment including the digital focal length of lens to y directions, the phase center to y directions are modified to lens phase modulation factor The amendment of position, the amendment to the digital focal length of lens in x directions and the amendment of the phase centre location to x directions, determine it The digital focal length of lens f ' in revised y directionsy,mnWith the phase centre location in revised y directionsIt is and revised The digital focal length of lens f ' in x directionsx,mnWith the phase centre location in revised x directionsDiscussed in the 4th step of the invention Three-dimensional imaging unit all referring to the three-dimensional imaging unit for being not entirely in gray area, be not entirely in grey area for any The three-dimensional imaging unit in domain such as m rows n-th arrange three-dimensional imaging unit, using following steps to its lens phase modulation factor into Row is corrected.
The three-dimensional imaging unit that gray area is not entirely on 4.1 pairs of three-dimensional imaging regions carries out the numeral in y directions thoroughly The amendment of mirror focal length, solves the digital focal length of lens f ' in revised y directionsy,mn
The left surface in whole three-dimensional imaging region is regarded as to the Terahertz aperture coding two in two-dimensional imaging plane progress y directions Dimension target segment scanning imagery, the three-dimensional imaging unit for being not entirely in gray area on any three-dimensional imaging region, such as The m rows n-th for being not entirely in gray area arrange three-dimensional imaging unit, the Terahertz in the coded imaging systems of Terahertz aperture The terahertz wave beam that dual-mode antenna is sent (is not entirely in grey in the three-dimensional imaging unit being not entirely in gray area The m rows n-th in region arrange three-dimensional imaging unit) coverage of right flank is Sh, terahertz wave beam is not entirely at this It is s that the m rows n-th of gray area, which arrange a coverage for three-dimensional imaging unit left surface,h
If the coverage of the three-dimensional imaging unit left surface obtained after lens phase modulation factor amendment is s "h, There are relation s "h=s 'h+ Δ s, Δ s can be solved according to the following formula:
Wherein hmnIt is the focus before amendment (before not carrying out lens phase modulation factor amendment, with three-dimensional imaging region Right flank is carried out when Terahertz aperture coding two dimension target blocked scan is imaged on three-dimensional imaging region as two-dimensional imaging plane M rows n-th arrange the focus that a three-dimensional imaging unit obtains) y-axis coordinate, can be tried to achieve according to triangle similarity relation:
Composite type (6), (7) and (8), can be in the hope of s "h
By triangle similarity relation obtain carry out lens phase modulation factor amendment after with the right flank in three-dimensional imaging region The focus that is obtained during the coding two dimension target blocked scan imaging of Terahertz aperture is carried out to three-dimensional imaging as two-dimensional imaging plane The distance of region right flank is d '
The digital focal length of lens f ' in revised y directions is calculated according to lens imaging formulay,mn
The amendment of the phase centre location in 4.2 pairs of y directions, solves the phase centre location in revised y directions
Correct forward and backward terahertz wave beam and arrange a three-dimensional imaging unit left side with being not entirely in the m rows n-th of gray area The intersection point that side upper end is intersected is respectively respectively h with the most short vertical distance of gray area1With h '1, similar to s 'hWith s "h, two There are relation h ' by person1=h1+ Δ s, wherein Δ s can be determined by formula (7).In addition, h1It can be tried to achieve by following formula:
Carried out after carrying out lens phase modulation factor amendment using the right flank in three-dimensional imaging region as two-dimensional imaging plane M rows n-th arrange a three-dimensional imaging unit and obtain on three-dimensional imaging region during the coding two dimension target blocked scan imaging of Terahertz aperture The y-axis coordinate h ' of the focus arrivedmnFor
, can be in the hope of the phase centre location in revised y directions according to triangle similarity relation
The imaging system of 4.3 Terahertz aperture coded imaging systems is (i.e. vertical in x directions (i.e. horizontal direction) and y directions Direction) on operating mode be identical, need to carry out lens phase modulation factor for each on three-dimensional imaging region to repair Positive three-dimensional imaging unit, the digital focal length of lens in x directions is modified using the method identical with 4.1 and 4.2 and The phase centre location in x directions is modified, determines the digital focal length of lens f ' in revised x directionsx,mnAnd phase center Position
Further, during three-dimensional blocked scan of the invention, it is completely covered in three-dimension object unit same When, adjacent non-targeted three-dimensional imaging unit can be also irradiated to by terahertz wave beam, and the echo information so obtained is redundancy. The radiation field reference signal matrix that redundancy can be included by construction carries out echo information reconstruct, then again by redundancy The method of removal reconstructs three-dimension object unit, and the reconstruction result of so each unit is free of redundancy, finally again will The reconstruction result combination of all units obtains complete target information.
The method have the benefit that:
The present invention is for the coding objective imaging of Terahertz aperture, the system of proposition blocked scan splicing, blocked scan Imaging can be reduced to calculate to calculator memory configuration requirement, in addition each three-dimensional imaging cell imaging process can be with parallel processing, most Spliced again afterwards, the areas imaging and speed of Terahertz aperture coded imaging can be effectively improved under the premise of high-resolution, is carried High imaging efficiency, reduces to calculator memory configuration requirement, effectively increases computer operational efficiency.
The present invention is during blocked scan, different from general two-dimentional blocked scan method, according to single 3 D into As the position of unit, renewal lens phase modulation factor is calculated, it is possible to achieve the complete covering to three-dimensional imaging unit, realization pair The complete scanning in three-dimensional imaging region.
During the present invention is for the problem that three-dimensional imaging Cell Reconstruction, information redundancy, the appropriate radiation field that expands is with reference to letter Number matrix size, redundancy is first reconstructed and rejected again, obtain pure three-dimensional imaging unit information then in conjunction with echo vector, Reduce influence of the redundancy to three-dimensionalreconstruction process.
Brief description of the drawings
Fig. 1 is the two dimension target high resolution scanning imaging scheme schematic diagram based on Terahertz aperture coded imaging system.
Fig. 2 is that Terahertz aperture coding objective high-resolution blocked scan imaging problem illustrates schematic diagram.
Fig. 3 is the y direction schematic diagrams of Terahertz aperture coding objective high-resolution blocked scan imaging scheme of the present invention.
Fig. 4 is the x direction schematic diagrams of Terahertz aperture coding objective high-resolution blocked scan imaging scheme of the present invention.
Fig. 5 is the flow chart of specific embodiment provided by the invention.
Embodiment
In order to which technical scheme and advantage is more clearly understood, with reference to the accompanying drawings and embodiments, to this hair It is bright to be further elaborated.It should be appreciated that specific embodiment described herein is only used for explaining the present invention, it is not used to Limit the present invention.
With reference to Fig. 3, Fig. 3 is the y side of Terahertz aperture coding objective high-resolution blocked scan imaging scheme of the present invention To schematic diagram.Gray area y directions in Fig. 3 are highly that the vertical height in array code aperture y directions is h, gray area z side To length be array code aperture to the distance between three-dimensional imaging region right flank (i.e. array code aperture and three-dimensional imaging The sum of the distance b of region left surface and the thickness c in three-dimensional imaging region).Between Terahertz dual-mode antenna and array code aperture Distance be a, the distance of array code aperture and three-dimensional imaging region left surface is b, and three-dimensional imaging area thickness is c, three-dimensional The distance of imaging region right flank and the focal plane under two-dimensional imaging system is d.When three-dimensional imaging cell distribution is in gray area When interior, the y directions lens phase modulation factor that array code aperture needs load is identical with two-dimensional target imaging, and focal plane is still For the perpendicular parallel with three-dimensional imaging region.When the upper end of a certain three-dimensional imaging unit is higher than or lower end is less than gray area When (when being not entirely in there are three-dimensional imaging unit in gray area), the y directions under loading two-dimensional scan imaging system Lens phase modulation factor the three-dimensional imaging unit cannot be completely covered in y directions, it is necessary to which correcting lens phase modulation factor makes The three-dimensional imaging unit can be completely irradiated in y directions by obtaining terahertz wave beam.As shown in figure 3, the lens phase tune in y directions Before the factor processed is not corrected, terahertz wave beam is not entirely in the covering model of the three-dimensional imaging unit right flank in gray area at this Enclose for Sh, and the coverage of its left surface is s 'h.The y directions lens phase modulation factor of the three-dimensional imaging unit is repaiied After just so that the beam coverage of the three-dimensional imaging unit left surface is s "h, such THz wave can be irradiated in y directions The whole three-dimensional imaging unit, the distance of revised focus to the three-dimensional imaging unit right flank is d '.hmn、h′mnRespectively Correct the y-axis coordinate of forward and backward focus.h1、h′1Respectively correct forward and backward terahertz wave beam and be not entirely in gray area with this The interior intersection point of three-dimensional imaging unit left surface and the minimum vertical distance of gray area.What is provided in Fig. 3 is three-dimensional imaging list First upper end is higher than the situation of gray area.
With reference to Fig. 4, Fig. 4 is the x side of Terahertz aperture coding objective high-resolution blocked scan imaging scheme of the present invention To schematic diagram.Gray area x directions length in Fig. 4 is the lateral length v in the x directions in array code aperture, gray area z side To length be array code aperture to the distance between three-dimensional imaging region right flank (i.e. array code aperture and three-dimensional imaging The sum of the distance b of region left surface and the thickness c in three-dimensional imaging region).Between Terahertz dual-mode antenna and array code aperture Distance be a, the distance of array code aperture and three-dimensional imaging region left surface is b, and three-dimensional imaging area thickness is c, three-dimensional The distance of imaging region right flank and the focal plane under two-dimensional imaging system is d.When three-dimensional imaging cell distribution is in gray area When interior, the lens phase modulation factor in the x directions that array code aperture needs load is identical with two-dimensional target imaging, and focal plane is still It is so the perpendicular parallel with three-dimensional imaging region.When a certain three-dimensional imaging unit is when being not entirely in gray area, The lens phase modulation factor in the x directions under loading two-dimensional scan imaging system cannot completely cover the three-dimensional imaging in x directions Unit enables terahertz wave beam to be completely irradiated to the three-dimensional imaging unit, it is necessary to correct lens phase modulation factor.Such as Shown in Fig. 4, before the lens phase modulation factor in x directions is not corrected, terahertz wave beam is not entirely in gray area at this The coverage of three-dimensional imaging unit right flank is Sv, and the coverage of its left surface is s 'v.To the x of the three-dimensional imaging unit After direction lens phase modulation factor is modified so that the beam coverage of the three-dimensional imaging unit left surface is s "v, this Sample THz wave is irradiated to the whole three-dimensional imaging unit x directions, and revised focus is right to the three-dimensional imaging unit The distance of side is d '.vmn、v′mnRespectively correct the z-axis coordinate of forward and backward focus.v1、v′1Respectively correct forward and backward terahertz Hereby wave beam is not entirely in the intersection point of three-dimensional imaging unit left surface and the minimum level of gray area in gray area with this Distance.What is provided in Fig. 4 is the situation that three-dimensional imaging unit upper end is higher than gray area.
Lens phase modulation factor is determined by the digital focal length of lens and phase centre location loaded, therefore core is to ask Go out the revised digital focal length of lens and phase centre location.
With reference to Fig. 5, for the flow chart of the present invention, comprise the following steps:
The first step, determine Terahertz aperture coded imaging systems parameter and carry out subdivision to three-dimensional imaging region;
Known Terahertz aperture coded imaging systems parameter is as follows:The vertical height in the y directions in array code aperture is h, The lateral length v in the x directions in array code aperture;The distance between Terahertz dual-mode antenna and array code aperture are a, array The distance of code aperture and three-dimensional imaging region left surface is b, and the thickness in three-dimensional imaging region is c.
According to the imaging resolution size of Terahertz aperture coded imaging systems and objective size by three-dimensional imaging Region is uniformly split into M × N number of three-dimensional imaging unit, wherein being split into N number of three-dimensional imaging on the y directions in three-dimensional imaging region Unit, is split into M three-dimensional imaging unit, the three-dimensional imaging unit number in whole three-dimensional imaging region is exactly M on its x direction ×N.The step is identical with subdivision method in background technology.
Second step, using the right flank in three-dimensional imaging region as two-dimensional imaging plane carry out the two-dimentional mesh of Terahertz aperture coding Blocked scan imaging is marked, its y directions (vertical side is tentatively solved respectively for each three-dimensional imaging unit on three-dimensional imaging region To) lens phase modulation factor and x directions (horizontal direction) lens phase modulation factor, including the y directions numeral focal length of lens fy,mnWith y directions phase centre locationAnd x directions numeral focal length of lens fx,mnWith x directions phase centre locationIts Middle m and n represent that m rows n-th arrange a three-dimensional imaging unit, m=1,2 ... M, n=1,2 ... N on three-dimensional imaging region.
The imaging system of Terahertz aperture coded imaging systems is in x directions (i.e. horizontal direction) and y directions (i.e. vertical side To) on operating mode be identical.Therefore preliminary y direction (vertical direction) the lens phase modulation factors that solve are solved with preliminary The method of x directions (horizontal direction) lens phase modulation factor is identical.Below tentatively to solve y directions (vertical direction) lens phase Exemplified by the modulation factor of position, illustrate its specific solution procedure:
A three-dimensional imaging unit is arranged for any one three-dimensional imaging unit such as m rows n-th on three-dimensional imaging region, is tentatively asked Solve its y directions (vertical direction) lens phase modulation factor, including y directions numeral focal length of lens fy,mnWith y directions phase center PositionMethod is as follows:
The 2.1 preliminary digital focal length of lens f for solving y directionsy,mn
The two-dimentional mesh of Terahertz aperture coding in y directions is carried out using the right flank in three-dimensional imaging region as two-dimensional imaging plane Blocked scan imaging is marked, determines that the terahertz wave beam that the Terahertz dual-mode antenna in the coded imaging systems of Terahertz aperture is sent exists The coverage S of each three-dimensional imaging unit right flankh
According to the focal plane under each three-dimensional imaging unit right flank of triangle similarity relation calculating and two-dimensional imaging system (i.e. The right flank in three-dimensional imaging region carries out corresponding focal plane when Terahertz aperture encodes two-dimensional imaging as two-dimensional imaging plane) Between distance d:
Then, calculated according to lens imaging formula and carry out y using the right flank in three-dimensional imaging region as two-dimensional imaging plane During the Terahertz aperture coding two dimension target blocked scan imaging in direction, m rows n-th arrange its y direction of three-dimensional imaging unit and need The digital focal length of lens f to be loadedy,mn
The 2.2 preliminary phase centre locations for solving y directions
The two-dimentional mesh of Terahertz aperture coding in y directions is carried out using the right flank in three-dimensional imaging region as two-dimensional imaging plane When marking blocked scan imaging, to realize seamless spliced, the phase in two dimension target blocked scan imaging process between each scanning area The lens phase modulation factor phase centre location stepping length of loading is needed between adjacent three-dimensional imaging unitFor:
According to stepping lengthThree-dimensional imaging region m rows n-th can be calculated arrange its y direction of three-dimensional imaging unit to add The phase centre location of load
For each three-dimensional imaging unit on three-dimensional imaging region the method identical with 2.1 and 2.2 can be used complete Into its x directions (horizontal direction) lens phase modulation factor is tentatively solved, include the digital focal length of lens f in x directionsx,mnWith x side To phase centre location
3rd step, judge whether the corresponding lens phase modulation factor of each three-dimensional imaging unit needs in three-dimensional imaging region Correct.
Before judging whether the corresponding lens phase modulation factor of each three-dimensional imaging unit needs amendment, one is first determined Gray area:If sending parallel beam along the z-axis direction from whole array code aperture, gray area is exactly covering for the parallel beam Lid scope, the vertical height in gray area y directions are the vertical height h, the x of gray area in the y directions in array code aperture The lateral length in direction is the lateral length in the x directions in array code aperture, and the length in gray area z directions is array code Aperture is to the distance between three-dimensional imaging region right flank (the i.e. distance b of array code aperture and three-dimensional imaging region left surface The sum of with the thickness c in three-dimensional imaging region);
When three-dimensional imaging unit is completely in gray area, array code aperture needs the lens phase modulation loaded The factor need not be corrected, and the lens phase modulation factor that array code aperture needs load is identical with two-dimensional target imaging, Jiao Ping Face remains as the perpendicular parallel with three-dimensional imaging region.
When there is the three-dimensional imaging unit for being not entirely in gray area, the lens under loading two-dimensional scan imaging system Phase modulation factor cannot completely cover three-dimensional imaging region, it is necessary to correct lens phase modulation factor so that terahertz wave beam is shone It is mapped to whole three-dimensional imaging region.
4th step, to judging to need in the 3rd step to carry out the modified three-dimensional imaging unit of lens phase modulation factor to carry out The amendment of y directions (vertical direction) lens phase modulation factor solves and x directions (horizontal direction) lens phase modulation factor is repaiied It is positive to solve, determine the digital focal length of lens f ' in its revised y directiony,mnWith the phase centre location in revised y directionsAnd the digital focal length of lens f ' in revised x directionsx,mnWith the phase centre location in revised x directions
For judging to need to carry out the modified three-dimensional imaging unit of lens phase modulation factor in the 3rd step (i.e. in the presence of not It is completely in the three-dimensional imaging unit of gray area), as to arrange three-dimensional imaging unit endless for m rows n-th on three-dimensional imaging region Gray area is in entirely, then the lens phase modulation factor of the three-dimensional imaging unit is modified, wherein to lens phase tune The factor processed is modified the amendment including the digital focal length of lens to y directions, the amendment to the phase centre location in y directions, to x The amendment of the digital focal length of lens in direction and the amendment of the phase centre location to x directions, determine its revised y direction Digital focal length of lens f 'y,mnWith the phase centre location in revised y directionsAnd the digital lens in revised x directions Focal length f 'x,mnWith the phase centre location in revised x directions
The three-dimensional imaging unit that gray area is not entirely on 4.1 pairs of three-dimensional imaging regions carries out the numeral in y directions thoroughly The amendment of mirror focal length, solves the digital focal length of lens f ' in revised y directionsy,mn
The left surface in whole three-dimensional imaging region is regarded as to the Terahertz aperture coding two in two-dimensional imaging plane progress y directions Dimension target segment scanning imagery, the three-dimensional imaging unit for being not entirely in gray area on any three-dimensional imaging region, such as The m rows n-th for being not entirely in gray area arrange three-dimensional imaging unit, the Terahertz in the coded imaging systems of Terahertz aperture The terahertz wave beam that dual-mode antenna is sent is not entirely in the covering model of the three-dimensional imaging unit right flank in gray area at this Enclose for Sh, terahertz wave beam arranges a covering for three-dimensional imaging unit left surface in the m rows n-th for being not entirely in gray area Scope is s 'h
If the coverage of the three-dimensional imaging region left surface obtained after lens phase modulation factor amendment is s "h, There are relation s "h=s 'h+ Δ s, Δ s can be solved according to the following formula:
Wherein hmnIt is the focus before amendment (before not carrying out lens phase modulation factor amendment, with three-dimensional imaging region Right flank is carried out when Terahertz aperture coding two dimension target blocked scan is imaged on three-dimensional imaging region as two-dimensional imaging plane M rows n-th arrange the initial focus that a three-dimensional imaging unit obtains) y-axis coordinate, can be tried to achieve according to triangle similarity relation:
Composite type (6), (7) and (8), can be in the hope of s "h
By triangle similarity relation obtain carry out lens phase modulation factor amendment after with the right flank in three-dimensional imaging region The focus that is obtained during the coding two dimension target blocked scan imaging of Terahertz aperture is carried out to three-dimensional imaging as two-dimensional imaging plane The distance of region right flank is d '
The digital focal length of lens f ' in revised y directions is calculated according to lens imaging formulay,mn
The amendment of the phase centre location in 4.2 pairs of y directions, solves the phase centre location in revised y directions
Correct forward and backward terahertz wave beam and arrange a three-dimensional imaging unit left side with being not entirely in the m rows n-th of gray area The intersection point that side upper end is intersected is respectively respectively h with the most short vertical distance of gray area1With h '1, similar to s 'hWith s "h, two There are relation h ' by person1=h1+ Δ s, wherein Δ s can be determined by formula (7).In addition, h1It can be tried to achieve by following formula:
Carried out after carrying out lens phase modulation factor amendment using the right flank in three-dimensional imaging region as two-dimensional imaging plane M rows n-th arrange a three-dimensional imaging unit and obtain on three-dimensional imaging region during the coding two dimension target blocked scan imaging of Terahertz aperture The y-axis coordinate h ' of the focus arrivedmnFor
, can be in the hope of the phase centre location in revised y directions according to triangle similarity relation
The imaging system of Terahertz aperture coded imaging systems is in x directions (i.e. horizontal direction) and y directions (i.e. vertical side To) on operating mode be identical, need to carry out lens phase modulation factor amendment for each on three-dimensional imaging region Three-dimensional imaging unit, the digital focal length of lens in x directions is modified and to x using the method identical with 4.1 and 4.2 The phase centre location in direction is modified, and determines the digital focal length of lens f ' in revised x directionsx,mnAnd phase centre location
5th step, lens phase modulation factor comprehensive modification
To make three-dimensional imaging unit all completely be covered by terahertz wave beam in x directions and y directions, more revised x side To digital focal length of lens f 'x,mnWith the digital focal length of lens f ' in revised y directionsy,mn, select wherein larger conduct numeral The focal length of lens:
F '=max (f 'x,mn,f′y,mn) (14)
When Terahertz aperture coded imaging systems carry out blocked scan to three-dimensional imaging region, Terahertz aperture coded imaging The pth row in array code aperture in system, the lens phase modulation factor F of q row array element loading (15) formulaspq,mnTo control too M rows n-th arrange a three-dimensional imaging unit on hertz beam three-dimensional imaging region:
Wherein:P=1,2 ... P, q=1,2 ... Q, k=2 π fc/ c, fcFor the centre frequency of THz wave, c is the light velocity.Array P array element is uniform-distribution with code aperture in the vertical direction, that is, y directions, array code aperture is x side in the horizontal direction Q array element is uniform-distribution with upwards, whole array code aperture includes P × Q array element.(xpq,ypq) for array code aperture the P rows, the coordinate position of q row array element central points, p=1,2 ... P.
Further, during three-dimensional blocked scan of the invention, it is completely covered in three-dimension object unit same When, adjacent non-targeted three-dimensional imaging unit can be also irradiated to by terahertz wave beam, and the echo information so obtained is redundancy. The radiation field reference signal matrix that redundancy can be included by construction carries out echo information reconstruct, then again by redundancy The method of removal reconstructs three-dimension object unit, and the reconstruction result of so each unit is free of redundancy, finally again will The reconstruction result combination of all units obtains complete target information.
Below for Terahertz aperture coded imaging systems, given system parameter, obtains this by way of computer sim- ulation Lens phase modulation factor core parameter in invention, verifies the practicality of the method for the invention.
Determine the height h=0.50m on the vertical direction of array code aperture, the width in horizontal direction is also 0.5m, is erected Nogata includes 25 row array elements upwards, and 25 arrays member is included in horizontal direction, shares 625 array elements.Terahertz dual-mode antenna and battle array The distance b=0.75m of the level interval a=0.25m of code aperture, array code aperture and three-dimensional imaging region left surface are arranged, The thickness in three-dimensional imaging region is c=0.25m, and drop shadow spread of the three-dimensional imaging region in xoy planes is 1.8m × 1.8m, if Determine three-dimensional imaging unit right flank size s=0.1m.
Prime minister by three-dimensional imaging region transversely with longitudinally split into 18 × 18=324 three-dimensional imaging unit, for the 9th The three-dimensional imaging unit that row the 10th arranges, the method according to the invention determine its lens phase modulation factor core parameter numeral lens Focal length and phase centre location.Judge to understand the three-dimensional imaging unit in gray area, numeral can be obtained thoroughly by directly calculating Mirror focal length and phase centre location are respectively 0.2065m and (- 0.0109m, 0.0109m).
For the three-dimensional imaging unit of the 4th row the 5th row, judge to understand the three-dimensional imaging unit not in gray area, point Not Ji Suan x directions and y directions the digital focal length of lens and phase centre location, obtain f 'x=0.2127m, f 'y=0.2153m and (0.0978m,0.1196m).Again according to step 2 select the larger focal length of lens can determine the revised digital focal length of lens and Phase centre location is 0.2153m and (0.0978m, 0.1196m).
The above method provided according to the present invention can calculate the corresponding lens phase modulation of each three-dimensional imaging unit respectively The factor, removes redundancy after blocked scan by method proposed by the present invention again, and splicing obtains complete three-dimensional imaging mesh Mark.
In conclusion although the present invention is disclosed above with preferred embodiment, so it is not limited to the present invention, any Those of ordinary skill in the art, without departing from the spirit and scope of the present invention, when can make it is various change and retouch, therefore this hair Bright protection domain is subject to the scope defined depending on claims.

Claims (5)

1. a kind of Terahertz aperture encodes objective scan imaging method, it is characterised in that comprises the following steps:
The first step, determine Terahertz aperture coded imaging systems parameter and carry out subdivision to three-dimensional imaging region;
Second step, using the right flank in three-dimensional imaging region as two-dimensional imaging plane carry out Terahertz aperture coding two dimension target point Block scan is imaged, its y directions lens phase tune is tentatively solved respectively for each three-dimensional imaging unit on three-dimensional imaging region The factor processed and x directions lens phase modulation factor, including y directions numeral focal length of lens fy,mnWith y directions phase centre locationAnd x directions numeral focal length of lens fx,mnWith x directions phase centre locationWherein m and n represents three-dimensional imaging area M rows n-th arrange a three-dimensional imaging unit, m=1,2 ... M, n=1,2 ... N on domain;
3rd step, judge whether the corresponding lens phase modulation factor of each three-dimensional imaging unit needs to repair in three-dimensional imaging region Just;
4th step, to judging to need in the 3rd step to carry out the modified three-dimensional imaging unit of lens phase modulation factor to carry out y side Solved to the amendment of lens phase modulation factor and the lens phase modulation factor amendment of x directions solves, determine its revised y side To digital focal length of lens f 'y,mnWith the phase centre location in revised y directionsAnd the numeral in revised x directions Focal length of lens f 'x,mnWith the phase centre location in revised x directions
5th step, lens phase modulation factor comprehensive modification;
To make three-dimensional imaging unit all completely be covered by terahertz wave beam in x directions and y directions, more revised x directions Digital focal length of lens f 'x,mnWith the digital focal length of lens f ' in revised y directionsy,mn, select wherein larger conduct numeral lens Focal length:
F '=max (f 'x,mn,f′y,mn) (14)
When Terahertz aperture coded imaging systems carry out blocked scan to three-dimensional imaging region, Terahertz aperture coded imaging systems The pth row in middle array code aperture, the lens phase modulation factor F of q row array element loading (15) formulaspq,mnTo control Terahertz M rows n-th arrange a three-dimensional imaging unit on beam three-dimensional imaging region:
<mrow> <msub> <mi>F</mi> <mrow> <mi>p</mi> <mi>q</mi> <mo>,</mo> <mi>m</mi> <mi>n</mi> </mrow> </msub> <mo>=</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>i</mi> <mfrac> <mi>k</mi> <mrow> <mn>2</mn> <msup> <mi>f</mi> <mo>&amp;prime;</mo> </msup> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>p</mi> <mi>q</mi> </mrow> </msub> <mo>-</mo> <msubsup> <mover> <mi>x</mi> <mo>^</mo> </mover> <mrow> <mi>m</mi> <mi>n</mi> </mrow> <mo>&amp;prime;</mo> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>p</mi> <mi>q</mi> </mrow> </msub> <mo>-</mo> <msubsup> <mover> <mi>y</mi> <mo>^</mo> </mover> <mrow> <mi>m</mi> <mi>n</mi> </mrow> <mo>&amp;prime;</mo> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
Wherein:P=1,2 ... P, q=1,2 ... Q, k=2 π fc/ c, fcFor the centre frequency of THz wave, c is the light velocity;Array code P array element is uniform-distribution with aperture in the vertical direction, that is, y directions, array code aperture is on x directions in the horizontal direction Q array element is uniform-distribution with, whole array code aperture includes P × Q array element;(xpq,ypq) it is array code aperture pth row, The coordinate position of q row array element central points, p=1,2 ... P.
2. Terahertz aperture according to claim 1 encodes objective scan imaging method, it is characterised in that:The first step In:Known Terahertz aperture coded imaging systems parameter is as follows:The vertical height in the y directions in array code aperture is h, and array is compiled The lateral length in the x directions in code holes footpath;The distance between Terahertz dual-mode antenna and array code aperture are a, array code hole The distance of footpath and three-dimensional imaging region left surface is b, and the thickness in three-dimensional imaging region is c;
According to the imaging resolution size of Terahertz aperture coded imaging systems and objective size by three-dimensional imaging region M × N number of three-dimensional imaging unit uniformly is split into, wherein being split into N number of three-dimensional imaging list on the y directions in three-dimensional imaging region Member, is split into M three-dimensional imaging unit on its x direction, the three-dimensional imaging unit number in whole three-dimensional imaging region be exactly M × N。
3. Terahertz aperture according to claim 2 encodes objective scan imaging method, it is characterised in that:Second step In, a three-dimensional imaging unit is arranged for each three-dimensional imaging unit such as m rows n-th on three-dimensional imaging region, tentatively solves its y Direction lens phase modulation factor, including y directions numeral focal length of lens fY, mnWith y directions phase centre locationMethod is such as Under:
The 2.1 preliminary digital focal length of lens f for solving y directionsy,mn
The Terahertz aperture coding two dimension target point in y directions is carried out using the right flank in three-dimensional imaging region as two-dimensional imaging plane Block scan is imaged, and determines terahertz wave beam that the Terahertz dual-mode antenna in the coded imaging systems of Terahertz aperture is sent each three Tie up the coverage S of imaging unit right flankh
Calculated according to triangle similarity relation between the focal plane under each three-dimensional imaging unit right flank and two-dimensional imaging system away from From d:
<mrow> <mi>d</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>s</mi> <mi>h</mi> </msub> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mi>b</mi> <mo>+</mo> <mi>c</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>h</mi> <mo>-</mo> <msub> <mi>s</mi> <mi>h</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
Then, calculated according to lens imaging formula and carry out y directions using the right flank in three-dimensional imaging region as two-dimensional imaging plane The coding two dimension target blocked scan imaging of Terahertz aperture when, m rows n-th arrange its y direction of three-dimensional imaging unit need plus The digital focal length of lens f of loady,mn
<mrow> <msub> <mi>f</mi> <mrow> <mi>y</mi> <mo>,</mo> <mi>m</mi> <mi>n</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mi>a</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mi>b</mi> <mo>+</mo> <mi>c</mi> <mo>+</mo> <mi>d</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>a</mi> <mo>+</mo> <mi>b</mi> <mo>+</mo> <mi>c</mi> <mo>+</mo> <mi>d</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
The 2.2 preliminary phase centre locations for solving y directions
The Terahertz aperture coding two dimension target point in y directions is carried out using the right flank in three-dimensional imaging region as two-dimensional imaging plane It is seamless spliced between each scanning area in two dimension target blocked scan imaging process to realize when block scan is imaged, adjacent three The lens phase modulation factor phase centre location stepping length of loading is needed between dimension imaging unitFor:
<mrow> <mi>&amp;Delta;</mi> <mover> <mi>y</mi> <mo>^</mo> </mover> <mo>=</mo> <mfrac> <mrow> <mi>a</mi> <mo>&amp;CenterDot;</mo> <mi>s</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mi>b</mi> <mo>+</mo> <mi>c</mi> <mo>+</mo> <mi>d</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>(</mo> <mi>a</mi> <mo>+</mo> <mi>b</mi> <mo>+</mo> <mi>c</mi> <mo>+</mo> <mi>d</mi> <mo>)</mo> <mo>&amp;CenterDot;</mo> <mo>(</mo> <mi>b</mi> <mo>+</mo> <mi>c</mi> <mo>)</mo> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
According to stepping lengthThree-dimensional imaging region m rows n-th can be calculated and arrange the phase that its y direction of three-dimensional imaging unit need to load Position center
<mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mrow> <mi>m</mi> <mi>n</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mi>N</mi> <mo>+</mo> <mn>1</mn> <mo>-</mo> <mn>2</mn> <mi>n</mi> </mrow> <mn>2</mn> </mfrac> <mo>&amp;CenterDot;</mo> <mi>&amp;Delta;</mi> <mover> <mi>y</mi> <mo>^</mo> </mover> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
Operating mode of the imaging system of Terahertz aperture coded imaging systems on x directions and y directions is identical, for three It is saturating that each three-dimensional imaging unit on dimension imaging region uses the method completion identical with 2.1 and 2.2 tentatively to solve its x direction Mirror phase modulation factor, includes the digital focal length of lens f in x directionsx,mnWith the phase centre location in x directions
4. Terahertz aperture according to claim 3 encodes objective scan imaging method, it is characterised in that:3rd step In, judge whether the corresponding lens phase modulation factor of each three-dimensional imaging unit needs modified method such as in three-dimensional imaging region Under:
Before judging whether the corresponding lens phase modulation factor of each three-dimensional imaging unit needs amendment, a grey is first determined Region:If sending parallel beam along the z-axis direction from whole array code aperture, gray area is exactly the covering model of the parallel beam Enclose, the vertical height in gray area y directions is the vertical height h in the y directions in array code aperture, the x directions of gray area Lateral length be array code aperture x directions lateral length, the length in gray area z directions is array code aperture The distance between to three-dimensional imaging region right flank;
When three-dimensional imaging unit is completely in gray area, array code aperture needs the lens phase modulation factor loaded It need not correct, the lens phase modulation factor that array code aperture needs load is identical with two-dimensional target imaging;
When there is the three-dimensional imaging unit for being not entirely in gray area, the lens phase under loading two-dimensional scan imaging system Modulation factor cannot completely cover three-dimensional imaging region, it is necessary to correct lens phase modulation factor so that terahertz wave beam is irradiated to Whole three-dimensional imaging region.
5. Terahertz aperture according to claim 4 encodes objective scan imaging method, it is characterised in that:3rd step In, for judging to need in the 3rd step to carry out, lens phase modulation factor is modified any is not entirely in the three of gray area Dimension imaging unit such as m rows n-th arrange a three-dimensional imaging unit, its lens phase modulation factor is repaiied using following steps Just:
The three-dimensional imaging unit that gray area is not entirely on 4.1 pairs of three-dimensional imaging regions carries out digital lens Jiao in y directions Away from amendment, solve the digital focal length of lens f ' in revised y directionsy,mn
The left surface in whole three-dimensional imaging region is regarded as to the two-dimentional mesh of Terahertz aperture coding in two-dimensional imaging plane progress y directions Blocked scan imaging is marked, for any three-dimensional imaging unit such as m rows for being not entirely in gray area on three-dimensional imaging region A n-th row three-dimensional imaging unit, the terahertz wave beam that the Terahertz dual-mode antenna in the coded imaging systems of Terahertz aperture is sent It is S in the coverage for the three-dimensional imaging unit right flank that this is not entirely in gray areah, terahertz wave beam is endless at this It is s that the full m rows n-th in gray area, which arrange a coverage for three-dimensional imaging unit left surface,h′;
<mrow> <msubsup> <mi>s</mi> <mi>h</mi> <mo>&amp;prime;</mo> </msubsup> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <mi>c</mi> <mo>+</mo> <mi>d</mi> <mo>)</mo> <mo>&amp;CenterDot;</mo> <msub> <mi>s</mi> <mi>h</mi> </msub> </mrow> <mi>d</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
If the coverage of the three-dimensional imaging unit left surface obtained after lens phase modulation factor amendment is sh", there is pass It is sh"=sh'+Δ s, Δ s are solved according to the following formula:
<mrow> <mi>&amp;Delta;</mi> <mi>s</mi> <mo>=</mo> <mfrac> <mrow> <mi>c</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mo>|</mo> <msub> <mi>h</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> </msub> <mo>|</mo> <mo>-</mo> <mi>h</mi> <mo>/</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>b</mi> <mo>+</mo> <mi>c</mi> <mo>+</mo> <mi>d</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
Wherein hmnFor the y-axis coordinate of the focus before amendment, can be tried to achieve according to triangle similarity relation:
<mrow> <msub> <mi>h</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mrow> <mi>m</mi> <mi>n</mi> </mrow> </msub> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mi>a</mi> <mo>+</mo> <mi>b</mi> <mo>+</mo> <mi>c</mi> <mo>+</mo> <mi>d</mi> <mo>)</mo> </mrow> </mrow> <mi>a</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
Composite type (6), (7) and (8), tries to achieve s "h
By triangle similarity relation obtain carry out lens phase modulation factor amendment after using the right flank in three-dimensional imaging region as Two-dimensional imaging plane carries out the focus that is obtained during the coding two dimension target blocked scan imaging of Terahertz aperture to three-dimensional imaging region The distance of right flank is d '
<mrow> <msup> <mi>d</mi> <mo>&amp;prime;</mo> </msup> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>s</mi> <mi>h</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msubsup> <mo>&amp;CenterDot;</mo> <mi>b</mi> <mo>+</mo> <msubsup> <mi>s</mi> <mi>h</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msubsup> <mo>&amp;CenterDot;</mo> <mi>c</mi> <mo>-</mo> <mi>h</mi> <mo>&amp;CenterDot;</mo> <mi>c</mi> </mrow> <mrow> <mi>h</mi> <mo>-</mo> <msubsup> <mi>s</mi> <mi>h</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msubsup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
The digital focal length of lens f in revised y directions is calculated according to lens imaging formulay,mn
<mrow> <msubsup> <mi>f</mi> <mrow> <mi>y</mi> <mo>,</mo> <mi>m</mi> <mi>n</mi> </mrow> <mo>&amp;prime;</mo> </msubsup> <mo>=</mo> <mfrac> <mrow> <mi>a</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mi>b</mi> <mo>+</mo> <mi>c</mi> <mo>+</mo> <msup> <mi>d</mi> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> </mrow> <mrow> <mi>a</mi> <mo>+</mo> <mi>b</mi> <mo>+</mo> <mi>c</mi> <mo>+</mo> <msup> <mi>d</mi> <mo>&amp;prime;</mo> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
The amendment of the phase centre location in 4.2 pairs of y directions, solves the phase centre location in revised y directions
Correct forward and backward terahertz wave beam and arrange a three-dimensional imaging unit left surface with being not entirely in the m rows n-th of gray area The intersection point that upper end is intersected is respectively respectively h with the most short vertical distance of gray area1And h1', similar to sh' and sh", Liang Zhecun In relation h1'=h1+ Δ s, wherein Δ s can be determined by formula (7);In addition, h1It can be tried to achieve by following formula:
<mrow> <msub> <mi>h</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mi>b</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mo>|</mo> <msub> <mi>h</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> </msub> <mo>|</mo> <mo>-</mo> <mi>h</mi> <mo>/</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>b</mi> <mo>+</mo> <mi>c</mi> <mo>+</mo> <mi>d</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
Terahertz is carried out using the right flank in three-dimensional imaging region as two-dimensional imaging plane after carrying out lens phase modulation factor amendment M rows n-th arrange what a three-dimensional imaging unit obtained on three-dimensional imaging region during hereby aperture coding two dimension target blocked scan imaging The y-axis coordinate h ' of focusmnFor
<mrow> <msubsup> <mi>h</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> <mo>&amp;prime;</mo> </msubsup> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mfrac> <mrow> <msubsup> <mi>h</mi> <mn>1</mn> <mo>&amp;prime;</mo> </msubsup> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mrow> <mi>b</mi> <mo>+</mo> <mi>c</mi> <mo>+</mo> <msup> <mi>d</mi> <mo>&amp;prime;</mo> </msup> </mrow> <mo>)</mo> </mrow> </mrow> <mi>b</mi> </mfrac> <mo>+</mo> <mi>h</mi> <mo>/</mo> <mn>2</mn> </mrow> </mtd> <mtd> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>n</mi> </msub> <mo>&gt;</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mfrac> <mrow> <msubsup> <mi>h</mi> <mn>1</mn> <mo>&amp;prime;</mo> </msubsup> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mrow> <mi>b</mi> <mo>+</mo> <mi>c</mi> <mo>+</mo> <msup> <mi>d</mi> <mo>&amp;prime;</mo> </msup> </mrow> <mo>)</mo> </mrow> </mrow> <mi>b</mi> </mfrac> <mo>-</mo> <mi>h</mi> <mo>/</mo> <mn>2</mn> </mrow> </mtd> <mtd> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>n</mi> </msub> <mo>&gt;</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
, can be in the hope of the phase centre location in revised y directions according to triangle similarity relation
<mrow> <msubsup> <mover> <mi>y</mi> <mo>^</mo> </mover> <mrow> <mi>m</mi> <mi>n</mi> </mrow> <mo>&amp;prime;</mo> </msubsup> <mo>=</mo> <mfrac> <mrow> <mi>a</mi> <mo>&amp;CenterDot;</mo> <msubsup> <mi>h</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> <mo>&amp;prime;</mo> </msubsup> </mrow> <mrow> <mi>a</mi> <mo>+</mo> <mi>b</mi> <mo>+</mo> <mi>c</mi> <mo>+</mo> <msup> <mi>d</mi> <mo>&amp;prime;</mo> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
Operating mode of the imaging system of 4.3 Terahertz aperture coded imaging systems on x directions and y directions be it is identical, it is right Need to carry out the modified three-dimensional imaging unit of lens phase modulation factor in each on three-dimensional imaging region, use with 4.1 with And 4.2 identical method the digital focal length of lens in x directions is modified and the phase centre location in x directions is repaiied Just, the digital focal length of lens f in revised x directions is determinedx,mnAnd phase centre location
CN201711127463.4A 2017-11-15 2017-11-15 Terahertz aperture coding three-dimensional target scanning imaging method Active CN107942328B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711127463.4A CN107942328B (en) 2017-11-15 2017-11-15 Terahertz aperture coding three-dimensional target scanning imaging method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711127463.4A CN107942328B (en) 2017-11-15 2017-11-15 Terahertz aperture coding three-dimensional target scanning imaging method

Publications (2)

Publication Number Publication Date
CN107942328A true CN107942328A (en) 2018-04-20
CN107942328B CN107942328B (en) 2018-11-13

Family

ID=61932256

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711127463.4A Active CN107942328B (en) 2017-11-15 2017-11-15 Terahertz aperture coding three-dimensional target scanning imaging method

Country Status (1)

Country Link
CN (1) CN107942328B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109739268A (en) * 2018-12-18 2019-05-10 深圳市太赫兹科技创新研究院 Regulating device is imaged
CN110609297A (en) * 2019-09-19 2019-12-24 中国人民解放军军事科学院国防科技创新研究院 Three-dimensional target imaging method and device
CN116087954A (en) * 2023-04-06 2023-05-09 中国科学院空天信息创新研究院 Real aperture system close-range out-of-focus synthetic aperture imaging method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012000075A1 (en) * 2010-06-28 2012-01-05 Institut National D'optique Synthetic aperture imaging interferometer
CN106680796A (en) * 2017-01-11 2017-05-17 中国人民解放军国防科学技术大学 Method for reconstructing three-dimensional surface of planar holographic array target on basis of frequency interference
CN106772367A (en) * 2016-11-30 2017-05-31 中国人民解放军国防科学技术大学 Long distance staring imaging device and method in the coding high-resolution of Terahertz frequency range aperture
CN106950555A (en) * 2017-05-03 2017-07-14 中国人民解放军国防科学技术大学 A kind of Area Objects imaging method based on Terahertz aperture coded imaging system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012000075A1 (en) * 2010-06-28 2012-01-05 Institut National D'optique Synthetic aperture imaging interferometer
CN106772367A (en) * 2016-11-30 2017-05-31 中国人民解放军国防科学技术大学 Long distance staring imaging device and method in the coding high-resolution of Terahertz frequency range aperture
CN106680796A (en) * 2017-01-11 2017-05-17 中国人民解放军国防科学技术大学 Method for reconstructing three-dimensional surface of planar holographic array target on basis of frequency interference
CN106950555A (en) * 2017-05-03 2017-07-14 中国人民解放军国防科学技术大学 A kind of Area Objects imaging method based on Terahertz aperture coded imaging system

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109739268A (en) * 2018-12-18 2019-05-10 深圳市太赫兹科技创新研究院 Regulating device is imaged
CN109739268B (en) * 2018-12-18 2021-11-16 深圳市太赫兹科技创新研究院 Imaging adjustment device
CN110609297A (en) * 2019-09-19 2019-12-24 中国人民解放军军事科学院国防科技创新研究院 Three-dimensional target imaging method and device
CN116087954A (en) * 2023-04-06 2023-05-09 中国科学院空天信息创新研究院 Real aperture system close-range out-of-focus synthetic aperture imaging method

Also Published As

Publication number Publication date
CN107942328B (en) 2018-11-13

Similar Documents

Publication Publication Date Title
CN107942328B (en) Terahertz aperture coding three-dimensional target scanning imaging method
Alvarez et al. Application of deep learning image-to-image transformation networks to GPR radargrams for sub-surface imaging in infrastructure monitoring
Shen et al. Parallel sensing in metaverses: Virtual-real interactive smart systems for “6S” sensing
CN106772367B (en) Long distance staring imaging device and method in the coding high-resolution of Terahertz frequency range aperture
CN106324853B (en) A kind of double object distance lamination imaging methods of visible domain
CN105403889B (en) A kind of Terahertz frequency range aperture encodes high-resolution low coverage staring imaging device
CN107966710B (en) Three-dimensional target high-resolution imaging method based on terahertz aperture coding imaging
CN107918125B (en) Single-emitting multi-receiving terahertz aperture coding imaging device and method
CN106872978A (en) A kind of Electromagnetic Modeling emulation mode of complex scene
CN116229224A (en) Fusion perception method and device, electronic equipment and storage medium
Mergy et al. Vision-based neural scene representations for spacecraft
CN102096913A (en) Multi-strategy image fusion method under compressed sensing framework
Almanza-Medina et al. Deep learning architectures for navigation using forward looking sonar images
Hurt et al. Benchmark meta-dataset of high-resolution remote sensing imagery for training robust deep learning models in machine-assisted visual analytics
US20240054698A1 (en) Three-dimensional imaging method and apparatus and three-dimensional imaging device
CN107833180A (en) A kind of method using complex field neutral net rapid solving nonlinear electromagnetic inverse Problem
Sengupta et al. A review of recent advancements including machine learning on synthetic aperture radar using millimeter-wave radar
CN109917385A (en) Long-distance frequency domain rapid imaging method and device for ground SAR (synthetic aperture radar) scanned by rotary arm
Bittner et al. DSM building shape refinement from combined remote sensing images based on WNET-CGANS
Nazarczuk et al. AIM 2024 sparse neural rendering challenge: Methods and results
Saha et al. Turbulence strength C n2 estimation from video using physics-based deep learning
CN104931963A (en) Moving object microwave stare correlated imaging method
Perez-Eijo et al. A physical optics simulator for multireflector THz imaging systems
Hafner et al. Exploring the fusion of sentinel-1 sar and sentinel-2 msi data for built-up area mapping using deep learning
CN114973181B (en) Multi-view BEV (beam steering angle) visual angle environment sensing method, device, equipment and storage medium

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