JP3290845B2  Radar image processing device  Google Patents
Radar image processing deviceInfo
 Publication number
 JP3290845B2 JP3290845B2 JP07487395A JP7487395A JP3290845B2 JP 3290845 B2 JP3290845 B2 JP 3290845B2 JP 07487395 A JP07487395 A JP 07487395A JP 7487395 A JP7487395 A JP 7487395A JP 3290845 B2 JP3290845 B2 JP 3290845B2
 Authority
 JP
 Japan
 Prior art keywords
 image
 radar
 image cell
 radar image
 cell
 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.)
 Expired  Fee Related
Links
 230000001133 acceleration Effects 0.000 claims description 51
 281000190407 TrackingPoint companies 0.000 claims description 8
 238000002592 echocardiography Methods 0.000 claims description 4
 238000010586 diagrams Methods 0.000 description 6
 238000000034 methods Methods 0.000 description 6
 238000003672 processing method Methods 0.000 description 2
 230000001131 transforming Effects 0.000 description 2
 280000125387 Around The Point companies 0.000 description 1
 280000468441 Strongly, Ltd. companies 0.000 description 1
 230000005540 biological transmission Effects 0.000 description 1
 239000000460 chlorine Substances 0.000 description 1
 239000003086 colorants Substances 0.000 description 1
 230000018109 developmental process Effects 0.000 description 1
 238000003384 imaging method Methods 0.000 description 1
Description
[0001]
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a radar image processing apparatus for creating a radar image reflecting a target shape from a radar echo from a radar target.
[0002]
2. Description of the Related Art A conventional radar image processing apparatus mainly uses a processing method which is easy to implement in hardware, that is, obtains twodimensional shape information on a rangeDoppler plane from a radar echo from a radar target. Met.
On the other hand, as a processing method for obtaining target threedimensional shape information, “DA” is completely different from the processing for obtaining a radar image on the rangeDoppler plane.
LE A. AUSHERMAN, ADAM KOZMA, et al. Developments in
Radar Imaging, IEEE Transaction on Aerospace and
Electronic Systems, AES20 (Jul 1984), 363400 ".
[0004] However, in the method according to the abovementioned literature,
It is necessary to collect various data on the aspect direction of the radar as viewed from the radar target, which is strongly limited to the movement of the target. For this reason, it is necessary to collect data for a long time, and it is difficult to apply the method to a target other than a target that performs a motion capable of collecting various aspect directions.
For this reason, it is extremely difficult to obtain a threedimensional radar image in actual operation, and a technique for obtaining information on a threedimensional shape from a radar image on a rangeDoppler plane has been conventionally required. .
[0006]
As described above, there has been a demand for a technique for acquiring information on a threedimensional shape from a radar image on a rangeDoppler plane. SUMMARY OF THE INVENTION The present invention has been made to solve the abovementioned problem, and it is not necessary to collect data in various aspect directions of a radar when viewed from a radar target. It is an object of the present invention to provide a radar image processing apparatus capable of obtaining target threedimensional shape information using a radar image obtained.
[0007]
In order to achieve the above object, a radar image processing apparatus according to the present invention provides a twodimensional radar image consisting of a rangeDoppler plane created from radar echoes in an area of a designated size. Image cell selecting means for selecting, for each area, an image cell having the maximum amplitude in each area including the target portion, and for each image cell selected by this means, a radar image one unit time before An image cell tracking unit that outputs a position and a complex value of the point of the image cell having the maximum amplitude as a new selected image cell within an area around the point of the image cell that has been selected, Image cell acceleration calculating means for calculating the phase of the complex value of the image cell point sequence, calculating the second derivative of the phase with respect to time, converting the calculation result into acceleration, and outputting the acceleration; A threedimensional quantity (range, Doppler frequency, acceleration) obtained by inputting an image and acceleration, and adding an acceleration to a twodimensional position (range, Doppler frequency) on the rangeDoppler plane of the image cell tracking point and the image cell A threedimensional radar image creating unit that outputs the amount of (range position, Doppler frequency, acceleration, intensity) as the threedimensional radar image of the image cell tracking point by adjusting the intensity of the tracking point; 3D radar image display means which regards a range position, Doppler frequency, acceleration) as a coordinate point and displays a point indicating its intensity on an orthogonal threedimensional coordinate system represented by a range axis, a Doppler axis, and an acceleration axis. It is configured so that
[0008]
In the radar image processing apparatus having the above configuration,
The twodimensional radar image is divided into regions of a designated size, and an image cell having a maximum amplitude in each region including a target portion is selected for each region. For each of the selected image cells, one unit time In the area centered on the image cell point selected for the previous radar image, the position of the image cell having the largest amplitude is determined as a new selected image cell, and its position and complex value are obtained. Calculate the phase of the complex value of the cell point sequence, calculate the second derivative of the phase with respect to time,
The calculation result is converted into acceleration, and a threedimensional amount (range, Doppler frequency) obtained by applying acceleration to a twodimensional position (range, Doppler frequency) on the rangeDoppler plane of the image cell tracking point from the twodimensional radar image and acceleration. (Frequency, acceleration) and the intensity of the image cell tracking point
Doppler frequency, acceleration, and intensity) are created as a threedimensional radar image of the image cell tracking point, and the (range position, Doppler frequency, acceleration) of the threedimensional radar image is regarded as a coordinate point, and a range axis, a Doppler axis, A point indicating the intensity is displayed on a threedimensional coordinate system represented by an acceleration axis.
[0009]
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS An embodiment of the present invention will be described below in detail with reference to the drawings. FIG. 1 shows the configuration of a radar image processing apparatus according to the present invention. A twodimensional radar image generating apparatus 11 sequentially receives radar echoes from a radar receiver (not shown) and outputs the echoes on a rangeDoppler plane. A dimensional radar image is created continuously at regular intervals. The radar image created here has a complex value consisting of phase and amplitude for each image cell.
As the twodimensional radar image creating device 11, a conventional device can be used as it is. The radar image created by the device 11 is sent to an image cell selection device 12, an image cell tracking device 13, and an image cell acceleration calculation device 14.
The image cell selecting device 12 sequentially selects image cells to be used for the subsequent processing from the input radar image, and outputs a set of the position on the rangeDoppler plane and a complex value to the image cell tracking device 13. I do.
The image cell tracking device 13 selects each of the input radar images selected one unit time ago based on the set of the position on the rangeDoppler plane and the complex value from the image cell selection device 12. The image cell of each target part is extracted from the image cell, the movement over the image cell of the target part is tracked, and the image cell including the target part to be targeted is selected. The image cell acceleration calculation device 14 and the threedimensional radar image creation Output to the device 15.
The image cell acceleration calculating device 14 calculates the phase of the complex value of each selected image cell from the image cell tracking device 13, and calculates the phase of the phase for each selected image cell.
The differential is calculated, the calculation result is converted into acceleration, and the acceleration is output to the threedimensional radar image creation device 15.
The threedimensional radar image creating device 15 has two
A threedimensional radar image is created from one radar image from the threedimensional radar image creation device 11, the selected image cell from the image tracking device 13, and the acceleration from the image cell acceleration calculation device 14, and the threedimensional radar image display device 1
Output to 6.
The threedimensional radar image display device 16 displays the input threedimensional radar image on a rectangular coordinate system or a skew coordinate system, and changes or rotates the scale factor of the coordinate system according to an instruction from an operator as needed. Do.
The processing procedure of the above configuration will be described below with reference to FIGS. First, the twodimensional radar image creation device 11 continuously creates a radar image on a rangeDoppler plane at regular time intervals from radar echoes from radar targets. The radar image created here has a complex value consisting of phase and amplitude for each image cell.
In the image cell selecting device 12, the input 2
The dimensional radar image is divided into regions of a designated size as shown in FIG. 2, and for example, an image that exceeds a preset threshold for the amplitude of the image cell is detected as an image cell including a target part. The image cell having the largest amplitude in each region including the target portion is selected for each region. For example, a twodimensional radar image is displayed on a screen of a display device,
An image cell may be selected according to an instruction of the operator.
The image cell selected in this manner is represented by a set of a position on the rangeDoppler plane, a complex value, and a position of an area to which the set belongs. Image cell selection device 12
Is used only once at the start for one processing.
Subsequently, a point P with respect to the area A1 shown in FIG.
An example in which 1 is selected will be described with reference to FIG. In the image cell tracking device 13, the point P of the image cell selected for the radar image one unit time ago is selected.
In the area A2 centered at 1, the position and the complex value of the point P2 of the image cell having the maximum amplitude are output as a new selected image cell. Alternatively, a new selected image cell may be output by an operator inputting a target rotational motion or using target rotational motion information from a device not included in the present configuration.
This processing is performed for each of the radar images at fixed time intervals, and the selected image cell rows (P1, P2, P
3, ...). The image cell acceleration calculation device 14 calculates the phase of the complex value of the selected image cell point sequence, calculates the second derivative of the phase with respect to time, converts it into acceleration, and outputs it. For example, the phase is calculated for each of the previously selected image cell rows, and the phase row (PH1, PH2, PH3,
…).
Next, subtraction of adjacent phases, for example, (phase PH2phase PH1) is performed, and a phase difference sequence (PH2
PH1, PH3PH2, PH4PH3, ...). Assuming that each term of this phase difference sequence is Y and each time is X,
For example, the gradient of the phase difference with respect to time is obtained by the least square method (see FIG. 4). This inclination is converted into acceleration by the following equation.
Acceleration = inclination / (2π) · λ / 2 · 1 / Δt ^{2} [m / s ** 2] Here, Δt represents a time interval at which a radar image is created, and λ represents a radar transmission wavelength.
The threedimensional radar image creating device 15 receives one radar image, for example, a radar image at the start of processing, and acceleration as inputs, and moves the point P1 to a twodimensional position (range, Doppler frequency) on the rangeDoppler plane. The threedimensional amount (range, Doppler frequency, acceleration) to which the acceleration is applied and the intensity of the point P1 are combined to obtain (the range position of the point P1, the point P1).
Doppler frequency, acceleration, intensity of point P1)
As a threedimensional radar image.
The above processing is performed for all the regions in FIG. 2 to obtain a threedimensional radar image of the selected image cell in each region. This situation will be described with reference to FIGS.
The threedimensional radar image display device 16 regards the (range position, Doppler frequency, acceleration) of the threedimensional radar image as a coordinate point, and displays it on an orthogonal threedimensional coordinate system represented by a range axis, a Doppler axis, and an acceleration axis. Points indicating the intensity, for example, points of different colors depending on the intensity, are displayed. In addition, a display using a threedimensional skewed coordinate system according to an instruction of an operator,
The rotation of the coordinate system and the change of the scale factor of each axis may be performed.
Further, the Doppler frequency [Hz] → distance [m] or acceleration [m / s ** 2
] → It may be displayed on an orthogonal threedimensional coordinate system using the distance unit [m] as a unit.
Here, a threedimensional coordinate system Σ (X, Y, Z) fixed to a target rotating at a constant angular velocity vector ω_ (an underline represents a vector) together with the target is considered. However, for the sake of simplicity, it is assumed that the line of sight of the radar coincides with the Y axis and the angular velocity vector ω_ is included in the YZ plane.
This assumption does not lose generality.
Under these conditions, coordinates (R ', dR / dt, d ^{2} R / dt ^{2} ), which are observable by radar, and the coordinates shown in FIG. In the system, for example, a velocity vector v_ of one wing tip r_ is expressed as follows.
[0029]
(Equation 1) Further, an acceleration vector a_ due to rotation is expressed as follows.
[0030]
(Equation 2) The first derivative dR / dt with respect to the time of the distance R from the radar to the observable target is expressed as a distance R >>  r_ 
Since it matches the projection component of the velocity vector v_ in the line of sight of the radar, it is expressed as follows.
[0031]
(Equation 3) This is modified using the above velocity vector v_ to obtain the following equation.
[0032]
(Equation 4) This is the Doppler axis component of a conventional twodimensional radar image. Similarly, the second derivative d ^{2} R with respect to the time of the distance R
/ Dt ^{2} is expressed as follows.
[0033]
(Equation 5) Further, the distance R is expressed as follows.
[0034]
(Equation 6) This is the range axis component of a conventional twodimensional radar image. Since r_ reflects the target shape, the second term of the above equation will be considered below.
[0035]
(Equation 7) Here, the following two vectors are considered as definitions.
[0036]
(Equation 8) From equations (1) to (3), if the three vectors ey_, γ_, and ζ_ are linearly independent, r_ is a point (R ′,
dR / dt, d ^{2} R / dt ^{2} ).
Here, as shown in FIG. 5, γ_ is parallel to the X axis of the threedimensional coordinate system Σ (X, Y, Z), and as shown in the figure, the rotational angular velocity vector ω_ is not parallel to the Z axis. As far as possible, ζ_ includes a component in the Zaxis direction. This is the traditional range
When the radar image on the Doppler plane is obtained, the same conditions as in the case where information on the Doppler direction is lost when the rotational angular velocity vector ω_ is parallel to the range axis, that is, the Y axis.
That is, if ω is not parallel to the Y axis and the Z axis, the threedimensional coordinate system Σ (X, Y, Z) becomes the skewed coordinate system Σ ′ (γ, Y, ζ).
Skewed coordinate system Σ '
A point above the point (R ', dR / dt, d 2 R / dt 2) is observed by the radar as.
Therefore, by using the points (R ′, dR / dt, d ^{2} R / dt ^{2} ) consisting of the amounts observable by the radar, each axis of the threedimensional coordinate system Σ (X, Y, Z) is obtained. Is obtained, for example, by calculating the relative ratio of each point for each axis (FIG. 6).
(B)), or the rotational angular velocity vector ω_
, The inverse transformation of the above coordinate transformation is obtained, and a point on the oblique coordinate system Σ ′ is transformed into a point on the threedimensional coordinate system Σ. Since the technique for obtaining the inverse transform is a wellknown technique, its description is omitted here.
On the other hand, the radar image obtained on a normal rangeDoppler plane is the above (R, dR / d
Since the image is composed of the components of t), the target shape X
A projection component onto the Y plane, and information in the Zaxis direction has been lost (see FIG. 6A).
Since an arbitrary portion of a target on a radar image obtained on a normal rangeDoppler plane can be made threedimensional, there is no need to collect data in various aspects of the radar viewed from the radar target. .
Hereinafter, the processing operation will be described with reference to FIG. Each target portion included in the image cell selected by the image cell selecting device 12 is corrected by the image cell acceleration calculating device 14 for d ^{2} while the image cell tracking device 13 corrects the movement across the image cells between radar images.
R / dt ^{2} is calculated using the phase of the image cell, and this is combined with the position (R ′, dR / dt) of each image cell selected by the image cell selection device 12 to obtain a point ( R ′, dR / dt, d ^{2} R / dt ^{2} ).
The threedimensional radar image display device 16 uses a set of obtained points (R ', dR / dt, d ^{2} R / dt ^{2} ) in a threedimensional coordinate system using (ey_, γ_, ζ_) as base vectors. ο (ey _, γ_, ζ_). At the time of display, (ey_, γ_, ζ_) may be replaced with three orthogonal or linearly independent unit vectors. Alternatively, a set of points (R ′, dR / dt, d ^{2} R / dt ^{2} ) is converted into a point on the threedimensional coordinate system Σ (X, Y, Z) by giving information about the rotational angular velocity vector ω_ from outside. Later, it may be displayed.
As described above, according to the configuration of this embodiment, the data on the rangeDoppler plane for which the technology has been established can be obtained without collecting data in various aspects of the radar when viewed from the radar target. The threedimensional shape information of the target can be obtained from the radar image obtained in step (1).
[0044]
As described above, according to the present invention, a range in which a technique is established without collecting data in various aspects of the radar when viewed from a radar target.
It is possible to provide a radar image processing apparatus capable of obtaining target threedimensional shape information using a radar image obtained on a Doppler plane.
FIG. 1 is a block diagram showing a configuration of an embodiment of a radar image processing apparatus according to the present invention.
FIG. 2 is a diagram showing a twodimensional radar image for explaining the image cell selection processing of the embodiment.
FIG. 3 is a diagram showing a twodimensional radar image for explaining the image cell tracking processing of the embodiment.
FIG. 4 is a diagram showing a state in which a gradient of a phase difference with respect to time is obtained in the image cell acceleration calculation processing of the embodiment.
FIG. 5 is a vector diagram for explaining the threedimensional radar image processing of the embodiment.
FIG. 6 is a diagram for describing threedimensional radar image image processing of the embodiment.
DESCRIPTION OF SYMBOLS 11 ... 2D radar image creation apparatus 12 ... Image cell selection apparatus 13 ... Image cell tracking apparatus 14 ... Image cell acceleration calculation apparatus 15 ... 3D radar image creation apparatus 16 ... 3D radar image display apparatus
Continuation of the front page (58) Field surveyed (Int.Cl. ^{7} , DB name) G01S ^{7} /007/42 G01S 13/0013/95
Claims (8)
3D radar image creating means for outputting the amount of intensity as a 3D radar image of the image cell tracking point; (Range position, Doppler frequency, acceleration) of the 3D radar image are regarded as coordinate points, Doppler axis,
A radar image processing apparatus comprising: threedimensional radar image display means for displaying a point indicating the intensity on a threedimensional coordinate system represented by an acceleration axis.
A radar image processing apparatus as described in the above.
2. The radar image processing apparatus according to claim 1, further comprising a function of using a onedimensional rectangular coordinate system or a threedimensional skew coordinate system to further change and display a scale factor of each axis.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

JP07487395A JP3290845B2 (en)  19950331  19950331  Radar image processing device 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

JP07487395A JP3290845B2 (en)  19950331  19950331  Radar image processing device 
Publications (2)
Publication Number  Publication Date 

JPH08271612A JPH08271612A (en)  19961018 
JP3290845B2 true JP3290845B2 (en)  20020610 
Family
ID=13559900
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

JP07487395A Expired  Fee Related JP3290845B2 (en)  19950331  19950331  Radar image processing device 
Country Status (1)
Country  Link 

JP (1)  JP3290845B2 (en) 
Families Citing this family (3)
Publication number  Priority date  Publication date  Assignee  Title 

US7212149B2 (en) *  20040617  20070501  The Boeing Company  System, method and computer program product for detecting and tracking a moving ground target having a single phase center antenna 
JP5424667B2 (en) *  20080602  20140226  三菱電機株式会社  Image radar device 
JP6235557B2 (en) *  20130312  20171122  古野電気株式会社  Radar device 

1995
 19950331 JP JP07487395A patent/JP3290845B2/en not_active Expired  Fee Related
Also Published As
Publication number  Publication date 

JPH08271612A (en)  19961018 
Similar Documents
Publication  Publication Date  Title 

Wang et al.  Threedimensional ISAR imaging of maneuvering targets using three receivers  
US8107722B2 (en)  System and method for automatic stereo measurement of a point of interest in a scene  
JP5393318B2 (en)  Position and orientation measurement method and apparatus  
Livingston et al.  Magnetic tracker calibration for improved augmented reality registration  
US8320709B2 (en)  Information processing method and apparatus for calculating information regarding measurement target on the basis of captured images  
Kang et al.  3D scene data recovery using omnidirectional multibaseline stereo  
DE19882939C2 (en)  Method for determining the relative camera orientation position to produce visible 3D images  
JP4079690B2 (en)  Object tracking apparatus and method  
US7193636B2 (en)  Image processing device and method therefor and program codes, storing medium  
CN101901481B (en)  Image mosaic method  
CA1212166A (en)  Range/azimuth/elevation angle ship imaging for ordnance control  
US6570566B1 (en)  Image processing apparatus, image processing method, and program providing medium  
JP2016057108A (en)  Arithmetic device, arithmetic system, arithmetic method and program  
US6456731B1 (en)  Optical flow estimation method and image synthesis method  
EP1724592B1 (en)  System for estimating the speed of an aircraft and its application to the detection of obstacles  
Johnson et al.  Overview of terrain relative navigation approaches for precise lunar landing  
EP1505408B1 (en)  A method for SAR processing without INS data  
Zhang et al.  Focusing parallel bistatic SAR data using the analytic transfer function in the wavenumber domain  
You et al.  Hybrid inertial and vision tracking for augmented reality registration  
US4084158A (en)  Method of operating synthetic aperture radar  
JP3338618B2 (en)  Display method and display device for real space image and virtual space image  
US7911503B2 (en)  Information processing apparatus and information processing method  
US6445815B1 (en)  Measurement of depth image considering time delay  
Chen et al.  Simulation of ISAR imaging of moving targets  
US7800645B2 (en)  Image display method and image display apparatus 
Legal Events
Date  Code  Title  Description 

FPAY  Renewal fee payment (prs date is renewal date of database) 
Free format text: PAYMENT UNTIL: 20080322 Year of fee payment: 6 

FPAY  Renewal fee payment (prs date is renewal date of database) 
Free format text: PAYMENT UNTIL: 20090322 Year of fee payment: 7 

FPAY  Renewal fee payment (prs date is renewal date of database) 
Free format text: PAYMENT UNTIL: 20100322 Year of fee payment: 8 

FPAY  Renewal fee payment (prs date is renewal date of database) 
Free format text: PAYMENT UNTIL: 20100322 Year of fee payment: 8 

FPAY  Renewal fee payment (prs date is renewal date of database) 
Free format text: PAYMENT UNTIL: 20110322 Year of fee payment: 9 

FPAY  Renewal fee payment (prs date is renewal date of database) 
Free format text: PAYMENT UNTIL: 20120322 Year of fee payment: 10 

FPAY  Renewal fee payment (prs date is renewal date of database) 
Free format text: PAYMENT UNTIL: 20130322 Year of fee payment: 11 

FPAY  Renewal fee payment (prs date is renewal date of database) 
Free format text: PAYMENT UNTIL: 20130322 Year of fee payment: 11 

FPAY  Renewal fee payment (prs date is renewal date of database) 
Free format text: PAYMENT UNTIL: 20140322 Year of fee payment: 12 

LAPS  Cancellation because of no payment of annual fees 