WO2009110096A1 - 合成開口処理システム及び合成開口処理方法 - Google Patents
合成開口処理システム及び合成開口処理方法 Download PDFInfo
- Publication number
- WO2009110096A1 WO2009110096A1 PCT/JP2008/054165 JP2008054165W WO2009110096A1 WO 2009110096 A1 WO2009110096 A1 WO 2009110096A1 JP 2008054165 W JP2008054165 W JP 2008054165W WO 2009110096 A1 WO2009110096 A1 WO 2009110096A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- data
- model data
- correlation value
- synthetic aperture
- correlation
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8902—Side-looking sonar
- G01S15/8904—Side-looking sonar using synthetic aperture techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9004—SAR image acquisition techniques
- G01S13/9017—SAR image acquisition techniques with time domain processing of the SAR signals in azimuth
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/06—Systems determining the position data of a target
- G01S15/08—Systems for measuring distance only
- G01S15/10—Systems for measuring distance only using transmission of interrupted, pulse-modulated waves
- G01S15/102—Systems for measuring distance only using transmission of interrupted, pulse-modulated waves using transmission of pulses having some particular characteristics
- G01S15/104—Systems for measuring distance only using transmission of interrupted, pulse-modulated waves using transmission of pulses having some particular characteristics wherein the transmitted pulses use a frequency- or phase-modulated carrier wave
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/10—Systems for measuring distance only using transmission of interrupted, pulse modulated waves
- G01S13/26—Systems for measuring distance only using transmission of interrupted, pulse modulated waves wherein the transmitted pulses use a frequency- or phase-modulated carrier wave
- G01S13/28—Systems for measuring distance only using transmission of interrupted, pulse modulated waves wherein the transmitted pulses use a frequency- or phase-modulated carrier wave with time compression of received pulses
- G01S13/282—Systems for measuring distance only using transmission of interrupted, pulse modulated waves wherein the transmitted pulses use a frequency- or phase-modulated carrier wave with time compression of received pulses using a frequency modulated carrier wave
Definitions
- the present invention relates to a synthetic aperture processing system and a synthetic aperture processing method such as radar and sonar.
- Sonar and radar can survey a wide range of remote targets, but have the disadvantage that the resolution in the direction of travel (azimuth direction) of the platform on which sonar and radar are mounted is reduced in long-distance surveys.
- a synthetic aperture process is known that improves the azimuth resolution by virtually extending the receiving array using a short receiving array that is easy to handle.
- a one-dimensional algorithm and a two-dimensional algorithm are known (for example, see Non-Patent Document 1).
- the one-dimensional algorithm first, the reflected wave is pulse-compressed in the time direction based on a conventional pulse compression technique (this is called range compression).
- range compression the deviation of the recording time of the reflected wave due to the influence of the platform movement called range curvature is corrected.
- the waveform amplitudes of the plurality of reflected waves are added and output as a sonar image or radar image with improved resolution in the range direction and azimuth direction.
- a two-dimensional correlation calculation is performed between a plurality of received reflected waves and a reference signal created based on an ideal received wave when reflected at a plurality of predetermined positions within the receivable range of the reflected wave.
- a correlation value is calculated for each of a plurality of predetermined positions within the receivable range. Based on these correlation value data, a sonar image or radar image with improved resolution in the range direction and azimuth direction is created.
- the processing in the azimuth direction is a simple summation operation, so the calculation load is small.
- the data interpolation of the reflected wave is necessary for the range curvature correction processing, and the accuracy error depends on the conditions. There is a problem of growing.
- the reflected wave is used as it is, and the reference signal and the two-dimensional correlation calculation are performed, so that the synthetic aperture processing can be performed with high accuracy.
- the correlation calculation with a large number of reference signals is necessary, There is a problem that becomes larger.
- the present invention provides a synthetic aperture processing system and a synthetic aperture processing method that can solve the above-described problems, perform synthetic aperture processing with high accuracy, and reduce the calculation load. Objective.
- the synthetic aperture processing system of the present invention generates a frequency-modulated wave whose frequency changes with time, and transmits the irradiation region to the irradiation region from a plurality of measurement positions, and the irradiation unit irradiates the irradiation region.
- Receiving means for receiving a plurality of reflected waves by a plurality of frequency-modulated waves, and range compression processing means for performing range compression on each of the plurality of reflected waves received by the receiving means and generating received data comprising a plurality of sink functions And a plurality of model data composed of a plurality of sink functions obtained by range-compressing an ideal received wave when the frequency-modulated wave is reflected at a plurality of predetermined positions in the irradiation region, and a reception generated by the range compression processing means
- Correlation value calculating means for calculating correlation value data indicating the degree of correlation between each of the model data and the received data based on the data, and a correlation value calculating means It comprises output means for outputting the correlation values calculated, the by.
- the synthetic aperture processing method of the present invention generates a frequency-modulated wave whose frequency changes with time, and transmits the irradiation region from a plurality of measurement positions.
- a compression processing step a plurality of model data composed of a plurality of sink functions obtained by range-compressing an ideal received wave when a frequency-modulated wave is reflected at a plurality of predetermined positions in the irradiation region, and a range compression processing step.
- correlation value data indicating the degree of correlation between each of the model data and the received data is calculated. It includes a correlation value calculation step, and outputting the correlation values calculated in the correlation value calculation step, the that.
- the range compression is performed, and the correlation value data is calculated using the reception data composed of the sync function having a smaller number of data than the reflected wave, so that the calculation load can be reduced. Further, since the correlation value data between the received data and the model data is calculated as in the conventional two-dimensional algorithm, it is not necessary to modify the range curvature of the conventional one-dimensional algorithm, and an increase in accuracy error can be avoided.
- the correlation value data is an index of reliability indicating whether or not the frequency modulation wave is actually reflected at a plurality of predetermined positions in the irradiation region of the frequency modulation wave. For example, certain model data and sample data A large correlation value indicates that an object is actually present at a predetermined position corresponding to the model data, and the possibility that the frequency-modulated wave is reflected at this position is high.
- the synthetic aperture processing system of the present invention is a model data calculation unit that calculates model data used by the correlation value calculation unit for each of a plurality of predetermined positions in the irradiation region in accordance with the frequency modulation wave generated by the transmission unit. It is preferable to further include
- model data is generated according to the frequency modulation wave that is actually used, so that synthetic aperture processing can be performed reliably.
- the model data calculation means extracts two points of data at a predetermined interval across the peak for each sync function constituting the model data, and these extracted data sets Is calculated as model data, and the correlation value calculation means, for each of the plurality of model data calculated by the model data calculation means, according to the position information of the two points of data for each sink function constituting the model data, A selection means for selecting two points of data as sample data from among the sink functions constituting the reception data, a sample data of the reception data selected by the selection means, and a plurality of model data sets, respectively And a cross-correlation calculating means for performing cross-correlation calculation and calculating correlation value data. It is.
- sample data is extracted from the received data and cross-correlation calculation is performed, so that the number of data for calculating correlation value data is further reduced, and the calculation load can be further reduced.
- the model data calculation means extracts two points of data at a predetermined interval across the peak for each sync function constituting the model data, and these extracted data sets Is calculated as model data, and the correlation value calculation means, for each of the plurality of model data calculated by the model data calculation means, according to the position information of the two points of data for each sink function constituting the model data, A selecting means for selecting two points of data as sample data from the sync functions constituting the received data, and a sum operation for calculating the correlation value data by calculating the sum of the sample data of the received data selected by the selecting means And means.
- the synthetic aperture processing can be performed with high accuracy and the calculation load can be reduced.
- the conventional sonar 1 is mounted on a platform 2 such as a ship or a submarine, and irradiates a pulsed sound wave to an irradiation region 3 such as the sea bottom while the platform 2 travels in one direction. Then, by measuring how many seconds later the reflected wave returns, a sonar image is generated, and an object in the direction of pulsed sound wave emission (range direction) is detected.
- a chirp pulse (chirp wave), which is a kind of frequency-modulated wave, is mainly used as the pulse sound wave to be irradiated. Thereby, the resolution of the sonar image in the range direction can be improved by applying a known pulse compression technique to the measured reflected wave.
- the synthetic aperture sonar 4 synthesizes a plurality of reflected waves received by the conventional sonar 1 at a plurality of positions in the traveling direction (azimuth direction) of the platform 2 to virtually move the receiving array of the sonar 1 in the azimuth direction. It was extended. That is, as shown in FIG. 1, the synthetic aperture sonar 4 is equivalent to having a larger antenna than the conventional sonar 1, and the irradiation areas 3a to 3a at a plurality of positions in the azimuth direction by the conventional sonar 1 are provided. The resolution in the azimuth direction can be improved for the common part 5 of 3c.
- the synthetic aperture system according to the present invention refers to the synthetic aperture sonar 4 shown in FIG.
- the “azimuth direction” used in the following description refers to the traveling direction of the platform 2 as shown in FIG.
- the “range direction” is orthogonal to the azimuth direction and refers to a direction in which pulsed sound waves are emitted from the platform 2.
- FIG. 2 shows the configuration of the synthetic aperture processing system according to the first embodiment of the present invention.
- the synthetic aperture processing system 10 of this embodiment includes a measurement condition input unit 11, a signal transmission unit (transmission unit) 12, an antenna (transmission unit, reception unit) 13, and a signal reception unit (reception unit). 14, range compression processing unit (range compression processing unit) 15, model data calculation unit (model data calculation unit) 16, sample data selection unit (correlation value calculation unit, selection unit) 17, cross-correlation calculation unit (correlation value calculation unit) , A cross-correlation calculating means) 18 and an image output unit (output means) 19.
- Components other than the antenna 13 of the synthetic aperture processing system 10 are physically configured as a computer system connected to the antenna 13.
- the computer system includes a CPU (Central Processing Unit), a RAM (Random Access Memory) that is a main storage device, a ROM (Read Only Memory), an auxiliary storage device such as a hard disk, an input device such as a keyboard and a mouse that are input devices, and a display. Output devices and the like.
- the image output unit 19 reads predetermined computer software on hardware such as a CPU and a RAM to operate the communication module, the input device, and the output device under the control of the CPU, and the RAM and the auxiliary storage device. This is realized by reading and writing data in.
- each functional block will be described based on the functional blocks shown in FIG.
- the measurement condition input unit 11 is a part for inputting various measurement conditions.
- the measurement conditions may be in a format input by an operator, for example, or in a format in which data stored in advance in a storage device or the like is read.
- the measurement conditions include the platform traveling speed, the chirp wave output interval, the chirp wave speed, the number of data in the azimuth direction for performing synthetic aperture processing, the reception sampling frequency, the transmission chirp rate, and the like. Data of these measurement conditions is transmitted to the signal transmission unit 12, the signal reception unit 14, and the model data calculation unit 16.
- the signal transmission unit 12 is a part that generates a frequency-modulated wave (chirp wave) whose frequency changes with time and irradiates the irradiation region from a plurality of measurement positions.
- the signal transmission unit 12 creates a chirp wave according to conditions such as the speed of the chirp wave among various measurement conditions received from the measurement condition input unit 11. Then, irradiation is performed from the antenna 13 at a plurality of positions in the azimuth direction according to conditions such as the platform traveling speed and the interval at which chirp waves are output.
- the synthetic aperture processing is performed with five positions n1 to n5 at equal intervals along the azimuth direction as one set.
- FIG. 3 shows that chirp waves c1 to c5 are irradiated from five positions n1 to n5 in the azimuth direction in the range direction (right direction in FIG. 3), respectively.
- the common portions 5 of the irradiation regions 3a to 3c at a plurality of positions in the azimuth direction are to be detected.
- the signal receiving unit 14 is a part that receives a plurality of reflected waves by a plurality of chirp waves irradiated by the signal transmitting unit 12.
- the signal receiving unit 14 performs filtering processing, amplification processing, digital processing, and the like on the reflected wave received by the antenna 13 according to conditions such as a reception sampling frequency among various measurement conditions received from the measurement condition input unit 11, These are transmitted to the range compression processing unit 15.
- FIG. 4 shows an example of the reflected waves r1 to r5 for the chirp waves c1 to c5 illustrated in FIG.
- the chirp waves c1 to c5 are reflected by the object and received by the signal receiving unit 14.
- Reflected waves r1 to r5 are shown.
- the signal receiving unit 14 acquires data indicating the reflected wave as data on a predetermined coordinate system.
- the reflected waves r1 to r5 are a data string composed of data including position coordinates in the range direction and azimuth direction, and signal intensity.
- the interval in the range direction of the data string is determined by conditions such as the reception sampling frequency set in the measurement condition input unit 11.
- the position coordinates in the azimuth direction of the reflected waves r1 to r5 correspond to the irradiation positions n1 to n5 of the chirp waves c1 to c5.
- the position coordinates in the range direction of the reflected waves r1 to r5 are determined by the signal receiving unit 14 according to the time difference from the transmission of the chirp waves c1 to c5 to the reception of the reflected waves r1 to r5. The smaller the time difference and the shorter the distance from the measurement position to the object, the closer to the irradiation position. In the example of FIG.
- the reflected wave r3 at the position n3 where the distance from the measurement position to the object is the shortest is arranged closest to the irradiation position, and the reflected waves r2 and r4 at the positions n2 and n4 are closer to the irradiation position than the reflected wave r3.
- the reflection waves r1 and r5 at the positions n1 and n5 are arranged farther from the irradiation position.
- the range compression processing unit 15 is a part that performs range compression on each of the plurality of reflected waves r1 to r5 received by the signal receiving unit 14 and generates reception data including a plurality of sync functions s1 to s5. Specifically, the range compression processing unit 15 uses a conventional pulse compression method such as autocorrelation calculation to compress the reflected waves r1 to r5 along the range direction to reduce the number of data to be received as a sync function. Convert to data.
- FIG. 5 shows received data waveforms s1 to s5 after the range of the reflected wave illustrated in FIG. 4 is range-compressed. As shown in FIG.
- the received data is waveforms s1 to s5 having a sinc function reflecting the time difference between the reflected waves r1 to r5 in FIG.
- These sinc function-like waveforms s1 to s5 are also data strings composed of data including the position coordinates in the range direction and the azimuth direction and the signal intensity, like the reflected waves r1 to r5.
- the generated reception data is transmitted to the sample data selection unit 17.
- the model data calculation unit 16 converts model data including a plurality of sink functions obtained by range-compressing an ideal received wave when the chirp signal is reflected at a plurality of predetermined positions in the irradiation region into a plurality of predetermined positions in the irradiation region. It is a part to calculate for each.
- “a plurality of predetermined positions in the irradiation region” are positions to be detected by the synthetic aperture processing, and in this embodiment, as shown in FIG. 6, a set of azimuth directions for performing the synthetic aperture processing. Of the positions n1 to n5, the areas a to e are divided into five detection target areas A in the range direction from the center position n3.
- the model data is calculated according to each of the areas a to e.
- the “ideal received wave” refers to a chirp wave that is assumed not to be scattered or refracted even after being reflected by an object. These measurement positions and received waves are set according to various measurement conditions received from the measurement condition input unit 11.
- these reflected waves are processed in the same manner as the reflected waves r1 to r5 acquired by the signal receiving unit 14 according to conditions such as the reception sampling frequency among various measurement conditions received from the measurement condition input unit 11.
- these reflected waves are converted into five sync function waveforms by the range compression processing. Similar to the waveforms s 1 to s 5 created by the range compression processing unit 15, these sync function-like waveforms are data strings including data including position coordinates in the range direction and azimuth direction, and signal intensity.
- the sync function has a single pulse peak as shown in FIGS. 5 and 6, and two points having the same interval as the pulse width are extracted from each sync function data string. Specifically, for example, two adjacent points with the maximum sum of signal strengths are extracted from the data string.
- the reception sampling frequency and the transmission chirp rate are adjusted in advance in the measurement condition input unit 11 so that two points having the same interval as the pulse width are obtained across the pulse peak of the sync function.
- model data (a) to (e) corresponding to each of the five areas a to e in the range direction at the position n3 in the azimuth direction are calculated.
- FIG. 6 shows an example of model data (a) to (e). As shown in FIG. 6, the model data (a) to (e) have different arrangements in the range direction according to the positions of the corresponding areas a to e in the range direction.
- the calculated model data (a) to (e) are transmitted to the sample data selection unit 17 and the cross correlation calculation unit 18.
- the sample data selection unit 17 responds to the position information of two points of data for each of the sink functions constituting the model data for each of the plurality of model data (a) to (e) calculated by the model data calculation unit 16. Thus, two points of data are selected as sample data from the sync functions s1 to s5 constituting the received data.
- the sample data selection unit 17 selects 10 points of data from the received data calculated by the range compression processing unit 15 according to the data arrangement in the range direction of each model data (a) to (e). .
- FIG. 7 shows an example of sample data selected according to the model data (a) to (e).
- the areas indicated by dotted lines in FIG. 7 and marked with a to e indicate the positions in the range direction of the model data (a) to (e) shown in FIG. 6 at the positions n1 to n5 in the azimuth direction. Yes.
- s1 to s5 of the received data two points of data included in the areas a to e are selected corresponding to the model data (a) to (e), respectively.
- 10 points of data selected at positions n1 to n5 in the azimuth direction are collected for each of the regions a to e, and 10 points of data for one model data are selected as sample data.
- the selected sample data is transmitted to the cross correlation calculation unit 18.
- the cross-correlation calculator 18 performs a cross-correlation calculation between the sample data of the received data selected by the sample data selector 17 and a plurality of model data datasets calculated by the model data calculator 16. This is a part for calculating correlation value data.
- one of the plurality of model data (a) to (e) is selected, and for each data of the selected model data, for example, the data near the origin at the position n1 is the first and the other data is A data string is created with the second, the data near the origin at position n5 being the ninth and the other data being the tenth.
- the sample data is also created by setting the data near the origin at position n1 as the first, the other data as the second, the data near the origin at position n5 as the ninth, and the other data as the tenth. Then, a cross-correlation calculation is performed using the data string of model data and sample data.
- the correlation values are calculated by multiplying the signal intensities of the data with the data in the same order of the data strings, and calculating the sum of the multiplied values. Such processing is performed for all model data (a) to (e), and in this embodiment, five correlation value data corresponding to the areas a to e in FIG. 6 are calculated as shown in FIG.
- FIG. 8 shows an example of correlation value data regarding the areas a to e divided in the range direction of the azimuth position n3 shown in FIG.
- Symbols a to e in FIG. 8 represent correlation values in the regions a to e in FIG.
- These correlation values a to e indicate the degree of correlation between the model data and the sample data, and become larger as the signal intensity patterns of the respective data are closer.
- a large correlation value between a certain model data and sample data indicates that there is a high possibility that an object actually exists in a corresponding region of the model data.
- the correlation value of the area c is the largest, which indicates that the possibility that an object exists in the area c is the highest.
- the image output unit 19 is a part that outputs a sonar image based on the correlation value data calculated by the cross correlation calculation unit 18. Specifically, as described above, at any position in the azimuth direction, the processing for calculating the correlation value data in the range direction of the position using the reception data up to two positions before and after the position is performed. The process is repeated while moving one point at a time, and correlation value data at each position in the azimuth direction is acquired. These correlation value data are arranged along the azimuth direction, and each correlation value is two-dimensionally arranged with respect to the azimuth direction and the range direction. Colors are assigned according to the magnitudes of these arranged correlation values. For example, in the case of black and white gradation, the color is specified so that the larger the correlation value, the closer to white. The image output unit 19 creates a sonar image based on the color information designated in this way.
- FIG. 9 shows an example of the sonar image (a) obtained by the conventional sonar 1 and the sonar image (b) output by the method of the present embodiment.
- the horizontal axes of the sonar images (a) and (b) indicate the position coordinates in the azimuth direction
- the vertical axis indicates the position coordinates in the range direction.
- the correlation value data calculated by the cross-correlation calculation unit 19 as shown in FIG. 8 is arranged in the vertical direction, and a plurality of correlation value data is arranged along the azimuth direction on the horizontal axis. They are lined up and arranged two-dimensionally.
- the sonar image (b) output by the method of the present embodiment has improved sharpness and contrast compared to the sonar image (a) obtained by the conventional sonar 1, and the resolution of the sonar image is improved.
- a chirp wave is created by the signal transmission unit 12 and irradiated from the antenna 13 at a predetermined timing (S102: transmission step).
- model data when the chirp signal is reflected at a plurality of predetermined positions in the irradiation area is calculated for each of the plurality of predetermined positions in the irradiation area according to the measurement conditions (S103).
- a plurality of reflected waves by the irradiated chirp wave are received by the signal receiving unit 14 (S104: reception step), and the range compression processing unit 15 performs range compression on each of the plurality of reflected waves.
- Received data composed of the sync function is generated (S105: range compression processing step).
- one sample data is selected from the model data by the sample data selection unit 17 (S106), and two points of data are respectively selected from the sink functions constituting the reception data according to the selected model data. And a total of 10 points are selected from the received data as sample data (S107).
- the cross-correlation calculation unit 18 performs cross-correlation calculation between the model data selected in step S106 and the sample data selected in step S107, and calculates correlation value data (S108: correlation). Value calculation step).
- the sample data selection unit 17 confirms whether or not all model data calculated in step S103 has already been selected (S109). If there is still unselected model data, the process returns to step S106, and the processes of steps S106 to S108 are repeated. If all the model data has been selected, the process proceeds to step S110.
- step S109 If it is determined in step S109 that all model data has been selected, the image output unit 19 outputs a sonar image based on the correlation value data (S110: output step).
- the model data calculation step (S103) may be executed in any order as long as it is after the measurement condition input step (S101) and before the model data selection step (S106).
- the range compression processing unit 15 performs range compression, and using reception data including a sink function having a smaller number of data than reflected waves, Since the cross correlation calculation unit 18 calculates correlation value data, the calculation load can be reduced. Further, as with the conventional two-dimensional algorithm, the cross-correlation calculation unit 18 calculates correlation value data between the received data and the model data, so that the range curvature correction of the conventional one-dimensional algorithm is not required and the accuracy error is large. Can be avoided.
- the model data calculation unit 16 creates a waveform similar to the reception data composed of the waveform of the sync function generated by the range compression processing unit 15, and these syncs are generated.
- Model data is created by extracting two points having the same interval as the pulse width across the pulse peak of the functional waveform. As shown in FIG. 11, these two points P1 and P2 are referred to as “two points P1 and P2 sandwich the pulse peak” and “the distance L1 between the two points is the same as the pulse width (L / 2)”.
- the sum of the signal intensity L2 at the point P1 and the signal intensity L3 at the point P2 is always substantially the same as the signal intensity L4 at the pulse peak (L2 + L3 ⁇ L4).
- the received data may have a situation where the data string does not include the data of the pulse peak position of the waveform of the sync function due to the sampling frequency.
- the correlation calculation is performed after the data of the pulse peak position is generated by performing oversampling and interpolation processing in the frequency domain.
- the synthetic aperture processing system 10 uses the data of the position of the pulse peak if the correlation calculation with the sample data (received data) is performed using two points having the same interval as the pulse width across the pulse peak. It is possible to obtain a correlation tendency that is substantially equivalent to that of the above case.
- the correlation between the two can be expressed without strictly matching the position of the pulse peak of the sync function included in the received data or the model data, and the processing for obtaining the data of the conventional pulse peak position Since it becomes unnecessary, the calculation load can be further reduced. Furthermore, since the model data is generated according to the information of the chirp wave actually used, the synthetic aperture processing can be reliably performed.
- FIG. 12 shows the configuration of a synthetic aperture processing system according to the second embodiment of the present invention.
- the synthetic aperture processing system 20 of the present embodiment is different from the synthetic aperture processing system 10 of FIG. 2 without extracting sample data from the reception data generated by the range compression processing unit 15. A cross-correlation calculation with model data is performed using all the received data.
- the functional blocks of the measurement condition input unit 11, the signal transmission unit 12, the antenna 13, the signal reception unit 14, the range compression processing unit 15, and the image output unit 19 are the synthetic aperture processing illustrated in FIG. It is the same as that of the system 10.
- the model data calculation unit 26 has five sync function-like shapes in each of the areas a to e obtained by dividing the detection target area A in the range direction from the position n3 in the azimuth direction into five. Until the waveform is created, the same processing as that of the model data calculation unit 16 in FIG. 2 is performed. After that, the model data calculation unit 26 collects the data strings that form the five sync function waveforms without selecting the data of the two points sandwiching the pulse peak from each of the sync function waveforms. Calculate as data. Therefore, the number of model data calculated by the model data calculation unit 26 is larger than that by the model data calculation unit 16 (10 points) in FIG. The calculated model data is transmitted to the cross-correlation calculator 28.
- the cross-correlation calculation unit 28 performs cross-correlation calculation between the reception data generated by the range compression processing unit 15 and the plurality of model data calculated by the model data calculation unit 26, thereby calculating correlation value data. .
- the signal transmission unit 12 when the signal transmission unit 12 generates a chirp wave and irradiates the irradiation region from a plurality of measurement positions, the signal reception unit 14 performs the plurality of irradiated chirp waves. Receive multiple reflected waves.
- the range compression processing unit 15 performs range compression on each of the plurality of received reflected waves, and generates reception data including a plurality of sync functions.
- the model data calculation unit 26 calculates a plurality of model data including a plurality of sink functions obtained by performing range compression on an ideal received wave when the chirp wave is reflected at a plurality of predetermined positions in the irradiation region.
- the cross-correlation calculating unit 28 performs cross-correlation calculation between the reception data generated by the range compression processing unit 15 and a plurality of model data sets, and calculates correlation value data.
- the image output unit 19 outputs the correlation value data calculated by the cross correlation calculation unit 18.
- the range compression is performed, and the correlation value data is calculated using the reception data composed of the sync function having a smaller number of data than the reflected wave, so that the calculation load can be reduced. Further, since the correlation value data between the received data and the model data is calculated as in the conventional two-dimensional algorithm, it is not necessary to modify the range curvature of the conventional one-dimensional algorithm, and an increase in accuracy error can be avoided. Further, since the model data is generated according to the information of the chirp wave actually used, the synthetic aperture processing can be surely performed.
- FIG. 13 shows a configuration of a synthetic aperture processing system according to the third embodiment of the present invention.
- the synthetic aperture processing system 30 of this embodiment is different from the synthetic aperture processing system 10 of FIG. 2 in that the sample data selected by the sample data selection unit 17 is subjected to cross-correlation calculation with model data.
- the correlation value data is calculated by simply performing a summation operation.
- the functional blocks of the measurement condition input unit 11, signal transmission unit 12, antenna 13, signal reception unit 14, range compression processing unit 15, model data calculation unit 16, and image output unit 19 in the synthetic aperture processing system 30 are shown in FIG. The same as that of the synthetic aperture processing system 10 shown in FIG.
- the summation calculation unit 38 performs a summation operation on the sample data of the received data selected by the sample data selection unit 17 and calculates correlation value data. Specifically, all the signal intensities of the sample data are added together and calculated as correlation value data.
- the signal transmission unit 12 when the signal transmission unit 12 generates a chirp wave and irradiates the irradiation region from a plurality of measurement positions, the signal reception unit 14 performs the plurality of irradiated chirp waves. Receive multiple reflected waves.
- the range compression processing unit 15 performs range compression on each of the plurality of received reflected waves, and generates reception data including a plurality of sync functions.
- the model data calculation unit 16 calculates a plurality of model data composed of a plurality of sink functions obtained by range-compressing an ideal received wave when the chirp wave is reflected at a plurality of predetermined positions in the irradiation region.
- the model data calculation unit 16 extracts two points of data at a predetermined interval across the peak for each sync function constituting the model data, and calculates these extracted data sets as model data. Then, the sample data selection unit 17 converts the received data for each of the plurality of model data calculated by the model data calculation unit 16 according to the positional information of the two points of data for each sync function constituting the model data. Two points of data are selected as sample data from the sink functions to be constructed. The summation calculation unit 38 calculates the sum of the sampled data of the received data selected by the sample data selection unit 17 to calculate correlation value data. The image output unit 19 outputs the correlation value data calculated by the cross correlation calculation unit 18.
- the range compression is performed, and the correlation value data is calculated using the reception data composed of the sync function having a smaller number of data than the reflected wave, so that the calculation load can be reduced.
- the correlation value data between the received data and the model data is calculated as in the conventional two-dimensional algorithm, it is not necessary to modify the range curvature of the conventional one-dimensional algorithm, and an increase in accuracy error can be avoided.
- the sample data is extracted from the received data and the cross-correlation calculation is performed, the number of data for calculating the correlation value data is further reduced, and the calculation load can be further reduced.
- the correlation value data is calculated simply by performing the sum calculation of these data after extracting the sample data, the calculation load can be further reduced.
- the model data is generated according to the information of the chirp wave actually used, the synthetic aperture processing can be reliably performed.
- the present invention is not limited to the above-described embodiment, and various modifications can be made.
- the number of positions in the azimuth direction where the synthetic aperture processing is performed may be a number other than the five positions shown in the above embodiment.
- the model data calculation units 16 and 26 create model data according to various measurement conditions input by the measurement condition input unit 11, but model data according to various measurement conditions in advance. May be prepared in a storage device or the like and appropriately selected according to the measurement conditions.
- the synthetic aperture sonar has been described, but the present invention is also applicable to a synthetic aperture radar. Furthermore, the present invention can also be applied to conventional pulse compression processing that does not perform synthetic aperture processing.
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
Description
大内和夫、"リモートセンシングのための合成開口レーダの基礎"、東京電機大学出版局、2004
まず、図1を参照して、本発明に係る合成開口システムが適用される合成開口ソーナー(Synthetic Aperture Soner:SAS)の概要について説明する。
次に、本発明に係る合成開口システムの第1の実施形態について説明する。図2に本発明の第1の実施形態に係る合成開口処理システムの構成を示す。図2に示すように、本実施形態の合成開口処理システム10は、計測条件入力部11、信号送信部(送信手段)12、アンテナ(送信手段、受信手段)13、信号受信部(受信手段)14、レンジ圧縮処理部(レンジ圧縮処理手段)15、モデルデータ算出部(モデルデータ算出手段)16、サンプルデータ選択部(相関値算出手段、選択手段)17、相互相関演算部(相関値算出手段、相互相関演算手段)18及び画像出力部(出力手段)19を含んで構成されている。
次に、本発明に係る合成開口システムの第2の実施形態について説明する。図12に本発明の第2の実施形態に係る合成開口処理システムの構成を示す。図12に示すように、本実施形態の合成開口処理システム20は、図2の合成開口処理システム10とは異なり、レンジ圧縮処理部15により生成された受信データからサンプルデータを抽出せずに、受信データの全てのデータを用いてモデルデータとの相互相関演算を行う。
次に、本発明に係る合成開口システムの第3の実施形態について説明する。図13に本発明の第3の実施形態に係る合成開口処理システムの構成を示す。図13に示すように、本実施形態の合成開口処理システム30は、図2の合成開口処理システム10とは異なり、サンプルデータ選択部17により選択されたサンプルデータに、モデルデータとの相互相関演算を行わず、単なる総和演算を行って相関値データを算出する。
Claims (5)
- 周波数が時間変化する周波数変調波を生成して、複数の計測位置から照射領域に照射する送信手段と、
前記送信手段により照射された複数の周波数変調波による複数の反射波を受信する受信手段と、
前記受信手段により受信された複数の反射波のそれぞれにレンジ圧縮を行い、複数のシンク関数からなる受信データを生成するレンジ圧縮処理手段と、
前記周波数変調波が前記照射領域内の複数の所定位置で反射した場合の理想的な受信波をレンジ圧縮した複数のシンク関数からなる複数のモデルデータと、前記レンジ圧縮処理手段により生成された前記受信データとに基づいて、前記モデルデータのそれぞれと前記受信データとの相関度合いを示す相関値データを算出する相関値算出手段と、
前記相関値算出手段により算出された相関値データを出力する出力手段と、
を備える合成開口処理システム。 - 前記送信手段により生成された前記周波数変調波に応じて、前記相関値算出手段で用いられる前記モデルデータを前記照射領域内の複数の所定位置ごとに算出するモデルデータ算出手段をさらに備える、請求項1に記載の合成開口処理システム。
- 前記モデルデータ算出手段が、前記モデルデータを構成するシンク関数ごとに、ピークを挟んだ所定の間隔の2点のデータを抽出し、これらの抽出されたデータセットを前記モデルデータとして算出し、
前記相関値算出手段が、
前記モデルデータ算出手段により算出された複数のモデルデータのそれぞれについて、該モデルデータを構成するシンク関数ごとの2点のデータの位置情報に応じて、前記受信データを構成するシンク関数の中からそれぞれ2点のデータをサンプルデータとして選択する選択手段と、
前記選択手段により選択された前記受信データのサンプルデータと、前記複数のモデルデータのデータセットとの間でそれぞれ相互相関演算を行い、相関値データを算出する相互相関演算手段と、
をさらに備える、請求項1または2に記載の合成開口処理システム。 - 前記モデルデータ算出手段が、前記モデルデータを構成するシンク関数ごとに、ピークを挟んだ所定の間隔の2点のデータを抽出し、これらの抽出されたデータセットを前記モデルデータとして算出し、
前記相関値算出手段が、
前記モデルデータ算出手段により算出された複数のモデルデータのそれぞれについて、該モデルデータを構成するシンク関数ごとの2点のデータの位置情報に応じて、前記受信データを構成するシンク関数の中からそれぞれ2点のデータをサンプルデータとして選択する選択手段と、
前記選択手段により選択された前記受信データのサンプルデータの総和演算を行い、相関値データを算出する総和演算手段と、
をさらに備える、請求項1または2記載の合成開口処理システム。 - 周波数が時間変化する周波数変調波を生成して、複数の計測位置から照射領域に照射する送信ステップと、
前記送信ステップにおいて照射された複数の周波数変調波による複数の反射波を受信する受信ステップと、
前記受信ステップにおいて受信された複数の反射波のそれぞれにレンジ圧縮を行い、複数のシンク関数からなる受信データを生成するレンジ圧縮処理ステップと、
前記周波数変調波が前記照射領域内の複数の所定位置で反射した場合の理想的な受信波をレンジ圧縮した複数のシンク関数からなる複数のモデルデータと、前記レンジ圧縮処理ステップにおいて生成された前記受信データとに基づいて、前記モデルデータのそれぞれと前記受信データとの相関度合いを示す相関値データを算出する相関値算出ステップと、
前記相関値算出ステップにおいて算出された相関値データを出力する出力ステップと、
を備える合成開口処理方法。
Priority Applications (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
NO08721583A NO2256515T3 (ja) | 2008-03-07 | 2008-03-07 | |
US12/920,195 US8421669B2 (en) | 2008-03-07 | 2008-03-07 | Synthetic aperture processing system and synthetc aperture processing method |
EP08721583.6A EP2256515B1 (en) | 2008-03-07 | 2008-03-07 | Synthetic aperture processing system and synthetic aperture processing method |
CA2716958A CA2716958C (en) | 2008-03-07 | 2008-03-07 | Synthetic aperture processing system and synthetic aperture processing method |
JP2010501748A JP5182676B2 (ja) | 2008-03-07 | 2008-03-07 | 合成開口処理システム及び合成開口処理方法 |
PCT/JP2008/054165 WO2009110096A1 (ja) | 2008-03-07 | 2008-03-07 | 合成開口処理システム及び合成開口処理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2008/054165 WO2009110096A1 (ja) | 2008-03-07 | 2008-03-07 | 合成開口処理システム及び合成開口処理方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2009110096A1 true WO2009110096A1 (ja) | 2009-09-11 |
Family
ID=41055675
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2008/054165 WO2009110096A1 (ja) | 2008-03-07 | 2008-03-07 | 合成開口処理システム及び合成開口処理方法 |
Country Status (6)
Country | Link |
---|---|
US (1) | US8421669B2 (ja) |
EP (1) | EP2256515B1 (ja) |
JP (1) | JP5182676B2 (ja) |
CA (1) | CA2716958C (ja) |
NO (1) | NO2256515T3 (ja) |
WO (1) | WO2009110096A1 (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013128878A1 (ja) * | 2012-02-28 | 2013-09-06 | 日本電気株式会社 | 信号処理装置、物体検知装置、物体検知機能付き装置および物体検知方法 |
JP6011746B1 (ja) * | 2015-12-03 | 2016-10-19 | 三菱電機株式会社 | 合成開口レーダ装置および信号処理装置 |
WO2017195297A1 (ja) * | 2016-05-11 | 2017-11-16 | 三菱電機株式会社 | レーダ処理装置 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2949868B1 (fr) * | 2009-09-04 | 2012-05-04 | Thales Sa | Procede d'amelioration de la localisation des points composant une image de radar de type sar. |
WO2011036721A1 (ja) * | 2009-09-25 | 2011-03-31 | トヨタ自動車株式会社 | 車載用レーダ装置 |
GB201322803D0 (en) * | 2013-12-20 | 2014-02-05 | Hydrason Solutions Ltd | Improvements in or relating to sonar apparatus |
CN113960602B (zh) * | 2021-12-22 | 2022-05-20 | 中科星睿科技(北京)有限公司 | 轨道误差信息生成方法、装置、电子设备和可读介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002311136A (ja) * | 2001-04-10 | 2002-10-23 | Nec Corp | 埋没物体探知方法及び埋没物体探知ソーナー |
JP2003028955A (ja) * | 2001-07-11 | 2003-01-29 | Mitsubishi Electric Corp | 合成開口レーダ装置 |
JP2003161779A (ja) * | 2001-11-26 | 2003-06-06 | Nec Corp | 合成開口レーダシステムおよびその高精細画像化処理方法 |
JP2003222678A (ja) * | 2002-01-30 | 2003-08-08 | Mitsubishi Heavy Ind Ltd | 水中音響探査方法及び水中音響探査システム |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4045795A (en) * | 1975-06-23 | 1977-08-30 | Nasa | Charge-coupled device data processor for an airborne imaging radar system |
EP0083107B1 (en) * | 1981-12-28 | 1989-03-29 | Nec Corporation | Synthetic aperture radar image processing system |
US5530447A (en) * | 1995-01-13 | 1996-06-25 | Delco Electronics Corp. | Blind-zone target discrimination method and system for road vehicle radar |
JPH09281299A (ja) | 1996-04-12 | 1997-10-31 | Ishikawajima Harima Heavy Ind Co Ltd | ガラス固化体閉じ込め検査装置 |
JP3256655B2 (ja) * | 1996-04-12 | 2002-02-12 | 大阪瓦斯株式会社 | 埋設物の探査方法および装置 |
-
2008
- 2008-03-07 WO PCT/JP2008/054165 patent/WO2009110096A1/ja active Application Filing
- 2008-03-07 CA CA2716958A patent/CA2716958C/en not_active Expired - Fee Related
- 2008-03-07 US US12/920,195 patent/US8421669B2/en not_active Expired - Fee Related
- 2008-03-07 JP JP2010501748A patent/JP5182676B2/ja not_active Expired - Fee Related
- 2008-03-07 EP EP08721583.6A patent/EP2256515B1/en not_active Not-in-force
- 2008-03-07 NO NO08721583A patent/NO2256515T3/no unknown
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002311136A (ja) * | 2001-04-10 | 2002-10-23 | Nec Corp | 埋没物体探知方法及び埋没物体探知ソーナー |
JP2003028955A (ja) * | 2001-07-11 | 2003-01-29 | Mitsubishi Electric Corp | 合成開口レーダ装置 |
JP2003161779A (ja) * | 2001-11-26 | 2003-06-06 | Nec Corp | 合成開口レーダシステムおよびその高精細画像化処理方法 |
JP2003222678A (ja) * | 2002-01-30 | 2003-08-08 | Mitsubishi Heavy Ind Ltd | 水中音響探査方法及び水中音響探査システム |
Non-Patent Citations (1)
Title |
---|
KAZUO OUCHI: "Rimoto Sensingu no tameno Gousei Kaikoh Rehdah no Kiso", 2004, TOKYO DENKI UNIVERSITY PRESS |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013128878A1 (ja) * | 2012-02-28 | 2013-09-06 | 日本電気株式会社 | 信号処理装置、物体検知装置、物体検知機能付き装置および物体検知方法 |
JPWO2013128878A1 (ja) * | 2012-02-28 | 2015-07-30 | 日本電気株式会社 | 信号処理装置、物体検知装置、物体検知機能付き装置および物体検知方法 |
US9945945B2 (en) | 2012-02-28 | 2018-04-17 | Nec Corporation | Signal processing apparatus, object detecting apparatus, apparatus provided with object detecting function, and object detecting method |
JP6011746B1 (ja) * | 2015-12-03 | 2016-10-19 | 三菱電機株式会社 | 合成開口レーダ装置および信号処理装置 |
WO2017094157A1 (ja) * | 2015-12-03 | 2017-06-08 | 三菱電機株式会社 | 合成開口レーダ装置および信号処理装置 |
WO2017195297A1 (ja) * | 2016-05-11 | 2017-11-16 | 三菱電機株式会社 | レーダ処理装置 |
JPWO2017195297A1 (ja) * | 2016-05-11 | 2018-06-14 | 三菱電機株式会社 | レーダ処理装置 |
Also Published As
Publication number | Publication date |
---|---|
JPWO2009110096A1 (ja) | 2011-07-14 |
EP2256515A4 (en) | 2015-05-06 |
CA2716958C (en) | 2014-07-15 |
NO2256515T3 (ja) | 2018-09-29 |
CA2716958A1 (en) | 2009-09-11 |
EP2256515B1 (en) | 2018-05-02 |
US20110032142A1 (en) | 2011-02-10 |
JP5182676B2 (ja) | 2013-04-17 |
EP2256515A1 (en) | 2010-12-01 |
US8421669B2 (en) | 2013-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5182676B2 (ja) | 合成開口処理システム及び合成開口処理方法 | |
EP1743190B1 (en) | Generating three-dimensional images using impulsive radio frequency signals | |
US7894303B2 (en) | Detection device, detection program and detection method | |
NL9201060A (nl) | Volumetrische en terreinbeeldsonar. | |
Klaucke | Sidescan sonar | |
JP2012168122A (ja) | ソーナーシステム、送信装置、受信装置、ソーナー用目標物特定方法、及びそのプログラム | |
EP3418770B1 (en) | Synthetic-aperture-radar signal processing device | |
Khodabandeloo et al. | Nonlinear crosstalk in broadband multi-channel echosounders | |
Hiroji | Extracting sonar relative beam patterns for multi-sector multibeam sonar | |
JP5664869B2 (ja) | 測定装置、測定システム、測定方法、及びプログラム | |
US20190265343A1 (en) | Method of compressing sonar data | |
KR101331333B1 (ko) | 바닥 지형을 측량하는 방법 및 장치 | |
Behar et al. | Optimization of sparse synthetic transmit aperture imaging with coded excitation and frequency division | |
CN106918809B (zh) | 快速干涉合成孔径声纳原始回波时域仿真方法 | |
Durofchalk et al. | Analysis of the ray-based blind deconvolution algorithm for shipping sources | |
KR101702581B1 (ko) | 레이다 스캔 영상 처리 장치 및 그 처리 방법 | |
RU2510045C2 (ru) | Фазовый гидролокатор бокового обзора | |
WO2011058527A1 (en) | Method and apparatus for processing sonar signals | |
KR20220112547A (ko) | 다중경로 환경에서 레이더의 표적 고각 측정 방법 및 장치 | |
JP2020106392A (ja) | ソーナー装置、方法、及びプログラム | |
RU2581898C1 (ru) | Способ измерения угловых координат цели | |
JP3223897B2 (ja) | 地中レーダ信号処理装置 | |
KR101928799B1 (ko) | 재질정보 생성을 위한 측심 소나의 후방산란 음압 정보 처리 장치 | |
CN113820715B (zh) | 一种采用阵元级多基地数据融合的目标定位方法 | |
Aleixo Miguel Cândido | Implementing a reference backscatter calibration technique on a multi-sector multibeam echosounder |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 08721583 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2010501748 Country of ref document: JP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2716958 Country of ref document: CA |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2008721583 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 12920195 Country of ref document: US |