CN109298417B - Building internal structure detection method and device based on radar signal processing - Google Patents

Building internal structure detection method and device based on radar signal processing Download PDF

Info

Publication number
CN109298417B
CN109298417B CN201811156010.9A CN201811156010A CN109298417B CN 109298417 B CN109298417 B CN 109298417B CN 201811156010 A CN201811156010 A CN 201811156010A CN 109298417 B CN109298417 B CN 109298417B
Authority
CN
China
Prior art keywords
building
domain data
internal structure
pixel point
image domain
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811156010.9A
Other languages
Chinese (zh)
Other versions
CN109298417A (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.)
Huanuo Xingkong Technology Co ltd
Hunan Huanuo Xingkong Electronic Technology Co ltd
Original Assignee
HUNAN NOVASKY ELECTRONIC 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 HUNAN NOVASKY ELECTRONIC TECHNOLOGY CO LTD filed Critical HUNAN NOVASKY ELECTRONIC TECHNOLOGY CO LTD
Priority to CN201811156010.9A priority Critical patent/CN109298417B/en
Publication of CN109298417A publication Critical patent/CN109298417A/en
Application granted granted Critical
Publication of CN109298417B publication Critical patent/CN109298417B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/887Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons
    • G01S13/888Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons through wall detection
    • 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

Abstract

The invention discloses a method and a device for detecting an internal structure of a building based on radar signal processing, wherein the method comprises the following steps: s1, fixedly arranging a radar system on the outer side of a building to be detected, constructing an equivalent array by the radar system through a multi-transmitting multi-receiving antenna array to form an equivalent aperture, and detecting a target area through the arranged radar system to obtain a multi-channel echo signal of the target area; s2, converting all the channel echo domain signals acquired in the step S1 into image domain data; s3, converting the image domain data obtained in the step S2 into a characteristic binary image, and performing inversion to obtain the internal structure of the building; the device comprises a multi-channel echo signal acquisition module, an image domain data conversion module and an internal structure inversion module. The invention has the advantages of simple realization operation, easy realization, high detection precision and efficiency, strong real-time property and the like.

Description

Building internal structure detection method and device based on radar signal processing
Technical Field
The invention relates to the technical field of radar detection, in particular to a method and a device for detecting an internal structure of a building based on radar signal processing.
Background
By utilizing the penetrability of electromagnetic waves, the through-wall imaging radar can detect human eyes such as wall bodies or optically invisible shielding areas, and can display an imaging interested target area after processing echo information, so that the through-wall imaging radar is widely applied to the military and civil fields of urban street fighting, anti-terrorism fighting, post-disaster rescue and the like. By combining the structure inversion in the building with the target detection result, the panoramic information of the detection range can be obtained, and the method has important practical values for fighting reference in the building, accurate terrorist offenders, target relative position determination and the like, so that the method for determining the building layout by inverting the structure in the building also becomes an important direction for current research and development.
The method is characterized in that a building internal structure is detected based on a radar detection mode, at present, a synthetic aperture array is generally adopted, multiple aperture echo data are obtained by moving an antenna, then the multiple aperture echo data are processed, certain building internal structure characteristics are obtained, namely, the imaging of a building wall is realized by moving the radar position and utilizing the synthetic aperture principle, typically, a through-wall radar adopting a transmit-receive or time division multiplexing transmit-receive antenna array mode needs to change the position by moving the antenna for many times in each visual angle so as to obtain the multiple aperture data. However, the method of using the synthetic aperture radar to detect the structure in the building is complex in operation and difficult to implement, the antenna needs to be moved for multiple times to change positions, the real-time performance is poor, real-time signal processing is not facilitated, radar signal processing can be performed after detection of each position is completed, the internal structure and the target state of the building cannot be reflected in real time, and the method cannot be applied to occasions with high real-time requirements.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: aiming at the technical problems in the prior art, the invention provides the method and the device for detecting the internal structure of the building based on radar signal processing, which are simple to operate, easy to realize, high in detection precision and efficiency and strong in real-time performance.
In order to solve the technical problems, the technical scheme provided by the invention is as follows:
a building internal structure detection method based on radar signal processing comprises the following steps:
s1, fixedly arranging a radar system on the outer side of a building to be detected, constructing an equivalent array by adopting a multi-transmitting multi-receiving antenna array to form an equivalent aperture, and detecting a target area through the arranged radar system to obtain a multi-channel echo signal of the target area;
s2, converting all the channel echo domain signals acquired in the step S1 into image domain data;
and S3, transforming the image domain data obtained in the step S2 into a characteristic binary image, and performing inversion to obtain the internal structure of the building.
As a further improvement of the method of the present invention, in step S2, coherent superposition processing is performed on all channel echo domain signals by using a back projection BP algorithm, and image domain data is obtained through conversion.
As a further improvement of the method of the present invention, the step of performing coherent superposition processing by using a back projection BP algorithm includes:
s21, dividing an imaging area into M multiplied by N pixel points;
s22, calculating the two-way time delay between each pixel point and all the receiving and transmitting antenna combinations;
s23, calculating the phase compensation of each pixel point according to the calculated double-pass time delay;
and S24, calculating data of each aperture domain in the equivalent aperture according to the phase compensation, accumulating the data of each aperture domain to obtain a calculation coherence stack result, and using the calculation coherence stack result as the image domain data obtained by conversion.
As a further improvement of the method of the present invention, the step S22 specifically obtains each pixel point (x) by following the following formulan,ym) With two-way time delay tau between all transceiver antenna combinationsk
Figure BDA0001818959040000021
Figure BDA0001818959040000022
Figure BDA0001818959040000023
Wherein, ymAnd xnThe coordinate values of the distance-direction and the azimuth-direction pixel points are respectively represented, M is 1,2T(k),yT(k) Is the transmit antenna coordinate, (x)R(k),yR(k) Is a coordinate of the receiving antenna, c represents the speed of light, dTDistance of pixel point to transmitting antenna, dRThe distance between the pixel point and the receiving antenna;
in the step S23, each pixel point (x) is calculatedn,ym) Phase compensation phase ofkThe expression of (a) is:
phasek(xn,ym)=exp(j2πfc·(dT+dR)/c)
wherein f iscCarrier frequency, c is speed of light;
in the step S24, a coherent superposition result I is obtained by calculation according to the following formulas(x, y) and as said image domain data;
Figure BDA0001818959040000024
wherein z isk(x,y)=echos(t-τk(x,y))×phasek(x, y) denotes the k-th aperture domain data, echos(t-τk(x, y)) is the time delay corrected channel echo signal.
As a further improvement of the method, in step S3, two-dimensional CFAR detection is performed on the image domain data to obtain a binary image of the internal structure characteristics of the building to be detected, and the internal structure of the building is obtained through inversion; the step of performing two-dimensional CFAR detection comprises: determining a detection threshold value in advance according to the difference between the human body echo and the building echo in the image domain, taking the unit to be detected as the center, comparing the image domain data obtained in the step S2 with the detection threshold value, and obtaining a binary image of the internal structure characteristics of the building to be detected according to the comparison result.
As a further improvement of the method of the present invention, the step of determining the detection threshold value includes: respectively setting a protection window for avoiding the self-shielding effect of the object, a background reference window for estimating the unit to be detected and a corresponding detection probability PAAccording to the set protection window, reference window and detection probability PAThe detection threshold value is determined to be obtained.
As a further improvement of the method of the present invention, in step S2, all the channel echo domain signals are converted into image domain data by using a range-doppler-based imaging algorithm.
As a further improvement of the method, in step S3, the image domain data is processed by using a maximum threshold method to obtain a binary image of the internal structure characteristics of the building to be detected, and the internal structure of the building is obtained by inversion, which specifically includes: and searching the maximum value max (x, y) in the array of the image domain data, taking the designated multiple a.max (x, y) of the maximum value as a comparison threshold value, wherein 0< a <1, comparing each pixel point value with the comparison threshold value, taking the pixel point value greater than the comparison threshold value as 1, and otherwise, taking 0 to obtain the internal structural feature binary image of the building to be detected.
A building internal structure detection device based on radar signal processing comprises:
the system comprises a multi-channel echo signal acquisition module, a multi-channel echo signal acquisition module and a multi-channel echo signal acquisition module, wherein the multi-channel echo signal acquisition module is used for fixedly arranging a radar system on the outer side of a building to be detected, the radar system adopts a multi-transmitting multi-receiving antenna array to construct an equivalent array so as to form an equivalent aperture, and a target area is detected through the arranged radar system so as to acquire a multi-channel echo signal of the target area;
an image domain data conversion module, configured to convert all the channel echo domain signals obtained in step S1 into image domain data;
and an internal structure inversion module, configured to transform the image domain data obtained in step S2 into a characteristic binary image, and perform inversion to obtain an internal structure of the building.
As a further improvement of the device, the image domain data conversion module adopts an image domain data conversion module based on a BP algorithm, and is used for performing coherent superposition processing on all channel echo domain signals by adopting a back projection BP algorithm to obtain image domain data through conversion;
the internal structure inversion module is based on CFAR detection and is used for carrying out two-dimensional CFAR detection on the image domain data to obtain a binary image of the internal structure characteristics of the building to be detected, and the internal structure of the building is obtained through inversion.
Compared with the prior art, the invention has the advantages that:
1. the method and the device for detecting the internal structure of the building based on radar signal processing acquire multi-channel echo signals by adopting an equivalent aperture mode at a fixed position, convert echo domain data into image domain data, and then invert the image domain data to obtain the internal structure of the building.
2. According to the method and the device for detecting the internal structure of the building based on radar signal processing, disclosed by the invention, the multi-channel echo signal can be obtained through one-time radar detection, the real-time detection of the internal structure of the building can be realized, and further, the method and the device can be directly fused with a target detection signal processing mode in the building, so that the real-time state positions of the internal structure of the building and a target can be accurately judged.
3. The building internal structure detection method and device based on radar signal processing further utilize equivalent aperture multi-channel echo signal data at a fixed position to carry out coherent BP imaging, multi-channel echo data can be obtained by one-time detection, real-time inversion can be carried out on the building internal structure, two-dimensional CFAR detection is carried out on image domain data according to human body echoes and building structure characteristics after BP imaging, the building internal structure can be accurately inverted, and the method and device can be conveniently and directly fused with target detection in a building.
Drawings
Fig. 1 is a schematic flow chart of an implementation of the method for detecting the internal structure of a building based on radar signal processing according to the embodiment.
Fig. 2 is a schematic diagram of an actual detection scenario in an embodiment of the present invention.
Fig. 3 is a schematic diagram of a specific implementation flow for implementing the detection of the internal structure of the building according to the embodiment.
Fig. 4 is a schematic diagram of a simulation result of measured data in a specific application embodiment.
Fig. 5 is a schematic diagram of the result obtained by the method of the present invention applied to the measured data in the specific application embodiment.
Fig. 6 is a schematic diagram of the detection results obtained by combining the movement of a single object in a building with the method of the present invention.
Detailed Description
The invention is further described below with reference to the drawings and specific preferred embodiments of the description, without thereby limiting the scope of protection of the invention.
As shown in fig. 1, the method for detecting the internal structure of a building based on radar signal processing in the present embodiment includes the steps of:
s1, fixedly arranging a radar system on the outer side of a building to be detected, constructing an equivalent array by the radar system through a multi-transmitting multi-receiving antenna array to form an equivalent aperture, and detecting a target area through the arranged radar system to obtain a multi-channel echo signal of the target area.
S2, converting all the channel echo domain signals acquired in the step S1 into image domain data;
and S3, transforming the image domain data obtained in the step S2 into a characteristic binary image, and performing inversion to obtain the internal structure of the building.
The radar system is specifically ultra-wideband through-wall radar equipment, and the interior of a building to be detected is detected by using a fixedly arranged through-wall radar, and a radar echo signal is returned; because the equivalent aperture is formed by adopting the array antenna design and structure on the premise of determining the working frequency band, the array element density can be increased to avoid grating lobes, and the imaging resolution can be improved, the embodiment utilizes the characteristics, the radar adopts the multi-transmitting multi-receiving antenna array structure to form the equivalent aperture, the multi-channel echo signal is obtained by adopting the equivalent aperture mode of a fixed position, the echo domain data is converted into image domain data, the image domain data is converted into a characteristic binary image, and then the internal structure of the building is obtained by inversion, thereby realizing the internal structure detection of the building by utilizing the equivalent aperture mode based on the fixed radar, compared with the traditional method of adopting a mobile radar to realize the detection by utilizing a synthetic aperture mode, a radar device does not need to be moved, the internal structure detection of the building can be realized by adopting one-time radar detection, and the operation complexity can be greatly reduced, the method improves the detection efficiency, improves the detection real-time property, and is particularly suitable for application occasions with higher real-time requirement.
The method can acquire multi-channel echo signals through one-time radar detection, can realize real-time detection of the internal structure of the building, can be executed independently and parallelly with the detection process of the targets in the building, and cannot influence the detection of the targets in the building, so that the detection of the internal structure of the building can be directly fused with the processing mode of the target detection signals in the building, and is favorable for accurately judging the real-time state positions of the internal structure and the targets of the building.
In the specific application example, when the multi-channel echo signal is extracted in step S1, the ultra-wideband radar apparatus is first fixed to the wall by using a tripod or other device, the ultra-wideband radar device adopts a multi-transmitting multi-receiving antenna array to construct an equivalent array, thereby forming an equivalent aperture, all channels sequentially transmit and receive electromagnetic wave signals, the multi-channel echo signal of the target area is obtained by directly carrying out one-time detection on the target area, the actual detection scene in the specific application embodiment is shown in figure 2, wherein the room is about 27.3m long and more than 8m wide, the radar is arranged to be close to the wall with a certain height from the ground, the center of the room is 4.4m away from the left window, 3 bearing columns of 0.5m by 0.5m are arranged in the room at positions of 7m, 14m and 21m, the left windows are 4m, 11m and 15m, 3 windows are arranged on the rear side of the room, and no wall exists between 19m and 22.65 m.
In step S2 of this embodiment, coherent superposition processing is performed on all the channel echo domain signals acquired in step S1 by using a back projection BP algorithm, and image domain data is obtained through conversion. The backward projection BP algorithm is adopted to convert the echo domain data into the image domain data, the operation is simple, the real-time performance is strong, and the detection real-time performance can be further improved.
In this embodiment, the step of performing coherent superposition processing by using a back projection BP algorithm includes:
s21, dividing an imaging area into M multiplied by N pixel points, wherein y ism(M-1, 2, …, M) and xn(m ═ 1,2, …, N) represents coordinate values of pixel points in the distance direction and the azimuth direction, respectively;
s22, calculating each pixel point (x)n,ym) A two-way time delay with all the transmitting and receiving antenna combinations;
s23, calculating the phase compensation of each pixel point according to the calculated double-pass time delay;
and S24, calculating data of each aperture domain in the equivalent aperture according to phase compensation, accumulating the data of each aperture domain to obtain a calculation coherence stack result, and using the calculation coherence stack result as the image domain data obtained by conversion.
Through the steps, the multi-channel echo data acquired by the equivalent aperture can be converted into image domain data, through coherent superposition processing, incoherent clutter, noise and side lobes among channels can be suppressed, and the suppression effect is good.
In step S22, each pixel point (x) is calculated according to the following formulas (1), (2) and (3)n,ym) Two-way time delay with all transmit-receive antenna combinations:
Figure BDA0001818959040000061
Figure BDA0001818959040000062
Figure BDA0001818959040000063
wherein, ym(M ═ 1,2,. multidot., M) and xnN is a coordinate value representing a pixel point in a distance direction and an orientation direction, k is 1,2,.., NUM is the channel label, NUM is the equivalent aperture number, (x)T(k),yT(k) Is the transmit antenna coordinate, (x)R(k),yR(k) Is a coordinate of the receiving antenna, c represents the speed of light, dTDistance of pixel point to transmitting antenna, dRThe distance between the pixel point and the receiving antenna;
in step S23, each pixel point (x) is calculatedn,ym) Phase compensation phase ofkThe expression of (a) is:
phasek(xn,ym)=exp(j2πfc·(dT+dR)/c) (4)
wherein f iscCarrier frequency, c is speed of light;
in step S24, a coherent addition result I is obtained by calculation according to the following formula (5)s(x, y) and as image domain data;
Figure BDA0001818959040000064
wherein z isk(x,y)=echos(t-τk(x,y))×phasek(x, y) denotes the k-th aperture domain data, echos(t-τk(x, y)) is the time delay corrected channel echo signal.
It can be understood that, the above-mentioned conversion of the echo domain signals of each channel into the image domain data may also adopt other methods according to actual requirements to reduce implementation complexity, etc.
In this embodiment, in step S3, performing two-dimensional CFAR detection on the image domain data to obtain a two-value graph of internal structure characteristics of the building to be detected, and performing inversion to obtain the internal structure of the building; the step of performing two-dimensional CFAR detection comprises: and determining a detection threshold value in advance according to the difference between the human body echo and the building echo in the image domain, taking the unit to be detected as a center, comparing the image domain data obtained in the step S2 with the detection threshold value, and obtaining a two-value image of the internal structure characteristic of the building to be detected according to the comparison result, wherein if the pixel value of the image domain data is greater than the detection threshold value, the two-value image is 1, otherwise, the two-value image is 0, and the two-value image is obtained. The difference between the human body echo and the building echo in the image domain comprises the difference between the distance width and the area, namely, the two-dimensional CFAR detection is carried out by utilizing the different characteristics of the distance width and the area of the human body echo and the building echo in the image domain.
The image domain data are converted into the characteristic binary image through the two-dimensional CFAR detection, the target can be accurately detected by utilizing the difference between the human body echo and the building echo in the image domain, so that the internal structure characteristics of the building are accurately inverted, and the threshold can be subjected to adaptive change according to the echo intensity by only selecting the size and the width of an appropriate window in the CFAR detection, so that the optimal detection in each local area is ensured, and the target detection probability can be improved.
In this embodiment, the step of determining the detection threshold includes: respectively setting a protection window for avoiding the self-shielding effect of the object, a background reference window for estimating the unit to be detected and a corresponding detection probability PAAccording to the set protection window, reference window and detection probability PAThe threshold value is determined according to the statistical value, and the specific expression is as follows:
Figure BDA0001818959040000071
in which alpha is a threshold factor, PAFor detecting the probability, N is the total number of the protection window and the reference window participating in the average, and T is the finally obtained threshold value.
After the threshold value is determined by the method, the image domain data obtained in the step S2 are respectively compared with the threshold value to be converted into a characteristic binary image, so that the target self-shielding effect can be fully considered, and the detection accuracy can be further improved.
It can be understood that, the obtaining of the binary image of the internal structure feature of the building to be detected from the image domain data may also adopt other methods according to actual requirements, for example, a maximum threshold method is used, that is, the image domain data obtained in step S2 is processed by using the maximum threshold method to obtain the binary image of the internal structure feature of the building to be detected, and the internal structure of the building is obtained by inversion, and the specific steps include: and searching the maximum value max (x, y) in the array of the image domain data, taking the designated multiple a.max (x, y) of the maximum value as a comparison threshold value, wherein 0< a <1, comparing each pixel point value with the comparison threshold value, taking the pixel point value greater than the comparison threshold value as 1, and otherwise, taking 0 to obtain the internal structural feature binary image of the building to be detected.
As shown in fig. 3, in the embodiment, when the detection of the structure in the building is realized based on the radar signal processing, specifically, after the position of the ultra-wideband radar is fixed, a target area is directly detected once, so as to obtain a multi-channel echo signal of the target area; then, carrying out coherent superposition processing on all the acquired channel echo domain signals by adopting a back projection BP algorithm to obtain image domain data; and then, carrying out two-dimensional CFAR detection on the obtained image domain data according to the structural characteristics of the building, and further performing inversion on the internal structure of the building. The method utilizes the equivalent aperture multi-channel echo signal data of the fixed position to carry out coherent BP imaging, can obtain the multi-channel echo data by one-time detection, thereby carrying out real-time inversion on the internal structure of the building, carrying out two-dimensional CFAR detection on image domain data according to the human body echo and the structural characteristics of the building after BP imaging, being capable of accurately inverting the internal structure of the building, having strong real-time performance, simple actual operation and high efficiency, and being capable of being directly fused with target detection in the building.
In order to verify the effectiveness of the invention, the actual measurement data shown in fig. 2 are processed by simulation processing and the method of the invention, fig. 4 is the result of the simulation processing of the actual measurement data shown in fig. 2, fig. 5 is the result of the detection of the internal structure of the building by the method of the invention, and as can be seen from fig. 4 and 5, the internal structure of the building detected by the method of the invention is basically matched with the actual detection scene, the protruding wall parts of the back wall, the column and the left window edge can be detected, and the internal structure of the building can be obtained by accurate inversion; the detection result obtained by further fusing the movement of the single target in the building with the method of the invention is shown in fig. 6, that is, the movement result of the single target in the building is fused with the detection structure of the internal structure of the building, wherein the curve represents the target state track, and the oblique line frame filling area represents the detection result of the structure of the building.
This embodiment is based on device of radar detection inversion building inner structure includes:
the system comprises a multi-channel echo signal acquisition module, a multi-channel echo signal acquisition module and a multi-channel echo signal acquisition module, wherein the multi-channel echo signal acquisition module is used for fixedly arranging a radar system on the outer side of a building to be detected, the radar system adopts a multi-transmitting multi-receiving antenna array to construct an equivalent array so as to form an equivalent aperture, and a target area is detected through the arranged radar system so as to acquire a multi-channel echo signal of the target area;
the image domain data conversion module is used for converting all the channel echo domain signals acquired by the multi-channel echo signal acquisition module into image domain data;
and the internal structure inversion module is used for transforming the image domain data obtained by the image domain data conversion module into a characteristic binary image and performing inversion to obtain the internal structure of the building.
In this embodiment, the image domain data conversion module adopts an image domain data conversion module based on a BP algorithm, and is configured to perform coherent superposition processing on all channel echo domain signals by using a back projection BP algorithm, so as to obtain image domain data through conversion;
the internal structure inversion module is based on CFAR detection and is used for carrying out two-dimensional CFAR detection on the image domain data to obtain a two-value image of the internal structure characteristics of the building to be detected, and the internal structure of the building is obtained through inversion.
The above-mentioned devices of this embodiment correspond to the above-mentioned building internal structure detection method based on radar signal processing one by one, and are not described herein again one by one.
The foregoing is considered as illustrative of the preferred embodiments of the invention and is not to be construed as limiting the invention in any way. Although the present invention has been described with reference to the preferred embodiments, it is not intended to be limited thereto. Therefore, any simple modification, equivalent change and modification made to the above embodiments according to the technical spirit of the present invention should fall within the protection scope of the technical scheme of the present invention, unless the technical spirit of the present invention departs from the content of the technical scheme of the present invention.

Claims (8)

1. A building internal structure detection method based on radar signal processing is characterized by comprising the following steps:
s1, fixedly arranging a radar system on the outer side of a building to be detected, constructing an equivalent array by adopting a multi-transmitting multi-receiving antenna array to form an equivalent aperture, and detecting a target area through the arranged radar system to obtain a multi-channel echo signal of the target area;
s2, performing coherent superposition processing on all the channel echo domain signals acquired in the step S1 by adopting a back projection BP algorithm to convert the signals into image domain data, wherein the step of performing coherent superposition processing by adopting the back projection BP algorithm comprises the following steps:
s21, dividing an imaging area into M multiplied by N pixel points;
s22, calculating the two-way time delay between each pixel point and all the receiving and transmitting antenna combinations;
s23, calculating the phase compensation of each pixel point according to the distance from the pixel point to the transmitting antenna and the distance from the pixel point to the receiving antenna;
s24, calculating each aperture domain data in the equivalent aperture according to the phase compensation and the two-way time delay, accumulating each aperture domain data to obtain a calculation coherence stack result, and using the calculation coherence stack result as the image domain data obtained by conversion;
and S3, transforming the image domain data obtained in the step S2 into a characteristic binary image, and performing inversion to obtain the internal structure of the building.
2. The method for detecting the internal structure of a building based on radar signal processing as claimed in claim 1, wherein said step S22 is implemented by calculating the two-way time delay τ according to the following formulak
Figure FDA0003507374590000011
Figure FDA0003507374590000012
Wherein, ymAnd xnRespectively represent distancesCoordinate values of pixel points in the direction and the azimuth direction, M is 1,2, …, M, N is 1,2, …, N, k is 1,2, …, NUM is a channel label, NUM is an equivalent aperture number, (x) is a value of a coordinate value of a pixel point in the direction and the azimuth direction, and (M) is a value of a coordinate value of a pixel point in the direction and the azimuth direction, NUM is a channel label, NUM is a value of an equivalent aperture number, (x) is a value of a pixel point in the direction and the azimuth direction, and (M) is a value of a pixel point in the direction and the azimuth direction, NUM is a value of a pixel point in the direction and the azimuth direction, M is a pixel point in the direction and the azimuth direction are 1,2, …, N is a pixel point in the same as a pixel point in the azimuth direction, and the azimuth direction of the pixel point in the azimuth direction and the azimuth directionT(k),yT(k) Is the transmit antenna coordinate, (x)R(k),yR(k) Is a coordinate of the receiving antenna, c represents the speed of light, dTDistance of pixel point to transmitting antenna, dRThe distance between the pixel point and the receiving antenna;
in the step S23, each pixel point (x) is calculatedn,ym) Phase compensation phase ofkThe expression of (a) is:
phasek(xn,ym)=exp(j2πfc·(dT+dR)/c)
wherein f iscCarrier frequency, c is speed of light;
in the step S24, a coherent superposition result I is obtained by calculation according to the following formulas(x, y) and as said image domain data;
Figure FDA0003507374590000021
wherein z isk(x,y)=echos(t-τk(x,y))×phasek(x, y) denotes the k-th aperture domain data, echos(t-τk(x, y)) is the time delay corrected channel echo signal.
3. The building internal structure detection method based on radar signal processing according to any one of claims 1 to 2, wherein in step S3, two-dimensional CFAR detection is performed on image domain data to obtain a binary image of the internal structure characteristics of the building to be detected, and the internal structure of the building is obtained through inversion; the step of performing two-dimensional CFAR detection comprises: determining a detection threshold value in advance according to the difference between the human body echo and the building echo in the image domain, taking the unit to be detected as the center, comparing the image domain data obtained in the step S2 with the detection threshold value, and obtaining a binary image of the internal structure characteristics of the building to be detected according to the comparison result.
4. The radar signal processing-based building interior structure detection method according to claim 3, wherein the step of determining the detection threshold value comprises: respectively setting a protection window for avoiding the self-shielding effect of the object, a background reference window for estimating the unit to be detected and a corresponding detection probability PAAccording to the set protection window, reference window and detection probability PAThe detection threshold value is determined to be obtained.
5. The building internal structure detection method based on radar signal processing according to any one of claims 1 to 2, wherein in step S2, all channel echo domain signals are converted into image domain data by using a range-doppler-based imaging algorithm.
6. The building internal structure detection method based on radar signal processing according to any one of claims 1 to 2, wherein in step S3, image domain data is processed by using a maximum value threshold method to obtain a binary image of the internal structure characteristics of the building to be detected, and the building internal structure is obtained by inversion, and the specific steps include: and searching the maximum value max (x, y) in the array of the image domain data, taking the designated multiple a.max (x, y) of the maximum value as a comparison threshold value, wherein 0< a <1, comparing each pixel point value with the comparison threshold value, taking the pixel point value greater than the comparison threshold value as 1, and otherwise, taking 0 to obtain the internal structural feature binary image of the building to be detected.
7. A building inner structure detection device based on radar signal processing is characterized by comprising:
the system comprises a multi-channel echo signal acquisition module, a multi-channel echo signal acquisition module and a multi-channel echo signal acquisition module, wherein the multi-channel echo signal acquisition module is used for fixedly arranging a radar system on the outer side of a building to be detected, the radar system adopts a multi-transmitting multi-receiving antenna array to construct an equivalent array so as to form an equivalent aperture, and a target area is detected through the arranged radar system so as to acquire a multi-channel echo signal of the target area;
the image domain data conversion module is used for performing coherent superposition processing on all the channel echo domain signals acquired by the multi-channel echo signal acquisition module by adopting a back projection BP algorithm to convert the signals into image domain data;
the internal structure inversion module is used for transforming the image domain data obtained by the image domain data conversion module into a characteristic binary image and performing inversion to obtain the internal structure of the building;
the image domain data conversion module comprises:
dividing an imaging area into M multiplied by N pixel points;
calculating the two-way time delay between each pixel point and all the transceiving antenna combinations;
calculating the phase compensation of each pixel point according to the distance from the pixel point to the transmitting antenna and the distance from the pixel point to the receiving antenna;
and calculating each aperture domain data in the equivalent aperture according to the phase compensation and the two-way time delay, accumulating each aperture domain data to obtain a calculation coherence stack result, and using the calculation coherence stack result as the image domain data obtained by conversion.
8. The radar signal processing-based building interior structure detecting apparatus according to claim 7,
the internal structure inversion module is based on CFAR detection and is used for carrying out two-dimensional CFAR detection on the image domain data to obtain a binary image of the internal structure characteristics of the building to be detected, and the internal structure of the building is obtained through inversion.
CN201811156010.9A 2018-09-30 2018-09-30 Building internal structure detection method and device based on radar signal processing Active CN109298417B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811156010.9A CN109298417B (en) 2018-09-30 2018-09-30 Building internal structure detection method and device based on radar signal processing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811156010.9A CN109298417B (en) 2018-09-30 2018-09-30 Building internal structure detection method and device based on radar signal processing

Publications (2)

Publication Number Publication Date
CN109298417A CN109298417A (en) 2019-02-01
CN109298417B true CN109298417B (en) 2022-04-22

Family

ID=65161371

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811156010.9A Active CN109298417B (en) 2018-09-30 2018-09-30 Building internal structure detection method and device based on radar signal processing

Country Status (1)

Country Link
CN (1) CN109298417B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110361722B (en) * 2019-07-17 2023-03-31 四川川大智胜软件股份有限公司 Multi-radar signal fusion method and system
CN110687534A (en) * 2019-10-28 2020-01-14 湖南华诺星空电子技术有限公司 Attitude identification method based on unmanned airborne life detector
CN110988872B (en) * 2019-12-25 2023-10-03 中南大学 Rapid identification method for detecting wall health state by unmanned aerial vehicle through-wall radar
CN114296074B (en) * 2021-12-29 2022-10-14 湖南华诺星空电子技术有限公司 Penetration type reconnaissance radar imaging method and radar system based on flight platform

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6177903B1 (en) * 1999-06-14 2001-01-23 Time Domain Corporation System and method for intrusion detection using a time domain radar array
US8040273B2 (en) * 2009-07-14 2011-10-18 Raytheon Company Radar for imaging of buildings
CN102176011B (en) * 2011-01-24 2014-01-01 陕西延长石油(集团)有限责任公司 Method for realizing three-dimensional coherent imaging by ground penetrating radar under near field condition
CN103197289B (en) * 2013-04-03 2015-02-18 电子科技大学 Multi-aperture phantom inhibition method applicable to through-the-wall radar imaging
CN104898115B (en) * 2015-05-11 2017-08-25 电子科技大学 Multi-object tracking method after a kind of through-wall radar imaging
CN105549011B (en) * 2015-12-16 2017-11-14 成都理工大学 A kind of unilateral 2 points of strabismus imaging method of building based on MIMO through-wall radars
CN106371093B (en) * 2016-08-16 2018-11-20 电子科技大学 Multi-target detection localization method based on building perspective radar imagery

Also Published As

Publication number Publication date
CN109298417A (en) 2019-02-01

Similar Documents

Publication Publication Date Title
CN109298417B (en) Building internal structure detection method and device based on radar signal processing
CN106772365B (en) A kind of multipath based on Bayes&#39;s compressed sensing utilizes through-wall radar imaging method
CN109298418B (en) Radar detection false alarm suppression method and device based on building internal structure characteristics
CN109581352B (en) Super-resolution angle measurement system based on millimeter wave radar
CN109188387B (en) Target parameter estimation method for distributed coherent radar based on interpolation compensation
CN107561507B (en) Clutter cancellation method for external radiation source radar
CN104749570B (en) It is a kind of to move constant airborne biradical synthetic aperture radar target localization method
Wang et al. Through-wall imaging of moving targets using UWB random noise radar
CN113687356B (en) Airborne multichannel circular SAR moving target detection and estimation method
Jia et al. A novel approach to target localization through unknown walls for through-the-wall radar imaging
CN109471097B (en) Through-wall radar signal optimization processing method and device
CN110568434A (en) multi-channel uniform acceleration trajectory SAR moving target two-dimensional speed estimation method
CN103777178A (en) Synchronization error compensation method, device and system
CN108020834B (en) Improved EDPCA (enhanced data analysis component analysis) -based moving target detection method and device and electronic equipment
CN109188436B (en) Efficient bistatic SAR echo generation method suitable for any platform track
CN109143235A (en) A kind of biradical forward sight synthetic aperture radar Ground moving target detection method
Lv et al. Coherence improving algorithm for airborne multichannel SAR-GMTI
CN110082760A (en) A kind of random radiation radar three-dimensional high-resolution imaging method
CN113820678B (en) Coherent accumulation method for high repetition frequency hopping signals of external radiation source radar
CN108710117A (en) A kind of synthetic aperture radar clutter cancellation method
Gao et al. Static background removal in vehicular radar: Filtering in azimuth-elevation-doppler domain
CN114994676A (en) One-station fixed type double-station low-frequency ultra-wideband SAR moving target imaging method
CN111796277B (en) Through-wall radar rapid imaging method based on unmanned aerial vehicle platform
CN110045374B (en) Chirp Scaling-based multi-angle ultra-long bunching SAR super-resolution imaging algorithm
CN108983192B (en) Radar moving target parameter estimation method based on GPS radiation source

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
CP01 Change in the name or title of a patent holder

Address after: Building B7, Lugu Enterprise Plaza, No. 27 Wenxuan Road, High tech Development Zone, Changsha City, Hunan Province, 410221

Patentee after: Huanuo Xingkong Technology Co.,Ltd.

Address before: Building B7, Lugu Enterprise Plaza, No. 27 Wenxuan Road, High tech Development Zone, Changsha City, Hunan Province, 410221

Patentee before: Hunan Huanuo Xingkong Electronic Technology Co.,Ltd.

CP01 Change in the name or title of a patent holder
CP03 Change of name, title or address

Address after: Building B7, Lugu Enterprise Plaza, No. 27 Wenxuan Road, High tech Development Zone, Changsha City, Hunan Province, 410221

Patentee after: Hunan Huanuo Xingkong Electronic Technology Co.,Ltd.

Address before: 410205 B7 Lugu Enterprise Plaza, 27 Wenxuan Road, Changsha High-tech Development Zone, Hunan Province

Patentee before: HUNAN NOVASKY ELECTRONIC TECHNOLOGY Co.,Ltd.

CP03 Change of name, title or address