JP7102907B2 - Synthetic aperture sonar, its signal processing method, and signal processing program - Google Patents

Synthetic aperture sonar, its signal processing method, and signal processing program Download PDF

Info

Publication number
JP7102907B2
JP7102907B2 JP2018085371A JP2018085371A JP7102907B2 JP 7102907 B2 JP7102907 B2 JP 7102907B2 JP 2018085371 A JP2018085371 A JP 2018085371A JP 2018085371 A JP2018085371 A JP 2018085371A JP 7102907 B2 JP7102907 B2 JP 7102907B2
Authority
JP
Japan
Prior art keywords
synthetic aperture
reference table
processing
signal processor
sound velocity
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
JP2018085371A
Other languages
Japanese (ja)
Other versions
JP2019191045A (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.)
NEC Corp
Original Assignee
NEC Corp
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 NEC Corp filed Critical NEC Corp
Priority to JP2018085371A priority Critical patent/JP7102907B2/en
Priority to FR1903989A priority patent/FR3080686B1/en
Publication of JP2019191045A publication Critical patent/JP2019191045A/en
Application granted granted Critical
Publication of JP7102907B2 publication Critical patent/JP7102907B2/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
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8902Side-looking sonar
    • G01S15/8904Side-looking sonar using synthetic aperture techniques
    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/523Details of pulse systems
    • G01S7/526Receivers
    • G01S7/53Means for transforming coordinates or for evaluating data, e.g. using computers

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Description

本開示は、水中航走体に搭載される合成開口ソーナー(Synthetic Aperture Sonar:SAS)及びその信号処理方法、並びに信号処理プログラムに関する。 The present disclosure relates to a Synthetic Aperture Sonar (SAS) mounted on an underwater vehicle, a signal processing method thereof, and a signal processing program.

合成開口ソーナーは、低周波音響に固有の問題である空間分解能の低さを、合成開口処理により改善したソーナーであり、海面等の精密な画像を取得するために使用されている。合成開口ソーナーの合成開口処理は、例えば、特許文献1に開示されている。 The synthetic aperture sonar is a sonar that improves the low spatial resolution, which is a problem peculiar to low-frequency acoustics, by synthetic aperture processing, and is used to acquire a precise image of the sea surface or the like. The synthetic opening treatment of the synthetic opening sonar is disclosed in, for example, Patent Document 1.

特開2007-333473号公報Japanese Unexamined Patent Publication No. 2007-333473

しかし、合成開口ソーナーの合成開口処理は、計算量が比較的多い。そのため、合成開口処理を行う信号処理器の処理負荷が大きく、処理時間を要し、リアルタイム化の阻害要因となっていた。 However, the synthetic aperture processing of the synthetic aperture sonar requires a relatively large amount of calculation. Therefore, the processing load of the signal processor that performs the synthetic aperture processing is large, the processing time is required, and it is an obstacle to real-time operation.

また、合成開口処理の処理時間を短縮するために、信号処理器に高速・高性能なCPU(Central Processing Unit)を採用すると、消費電力増大を招いていた。密閉空間である水中航走体の内部で合成開口処理を実行する場合、消費電力の増大は、内部温度の上昇につながるため、合成開口ソーナーの安定運用の観点からも看過し得ない要因となっていた。 Further, if a high-speed and high-performance CPU (Central Processing Unit) is adopted for the signal processor in order to shorten the processing time of the combined aperture processing, the power consumption is increased. When the synthetic aperture treatment is performed inside the underwater vehicle, which is a closed space, the increase in power consumption leads to an increase in the internal temperature, which is a factor that cannot be overlooked from the viewpoint of stable operation of the synthetic aperture sonar. Was there.

本開示の目的は、上述した課題を解決し、合成開口処理を行う信号処理器の処理負荷を低減し、処理速度を高速化することができる合成開口ソーナー及びその信号処理方法、並びに信号処理プログラムを提供することにある。 An object of the present disclosure is a synthetic aperture sonar, a signal processing method thereof, and a signal processing program capable of solving the above-mentioned problems, reducing the processing load of a signal processor that performs synthetic aperture processing, and increasing the processing speed. Is to provide.

一態様による合成開口ソーナーは、
音波の送受波を行う送受波器と、
前記送受波器で受波された受波音波のデータに対して合成開口処理を行う信号処理器と、を備え、
前記信号処理器は、
参照テーブルを選択し、
前記データに対して、2次元FFTを実行し、
前記2次元FFTの出力から、前記選択した参照テーブルに含まれるインデックスに対応する要素を抽出し、
前記抽出した要素に対して、前記選択した参照テーブルに含まれるテーブル係数を乗算し、
前記乗算の結果に対して、2次元IFFTを実行し、
前記2次元IFFTの出力を画像化する。
Synthetic aperture sonar according to one aspect
A transmitter / receiver that transmits / receives sound waves,
A signal processor that performs synthetic aperture processing on the data of the received sound wave received by the transmitter / receiver is provided.
The signal processor
Select a lookup table and
A two-dimensional FFT is executed on the data,
From the output of the two-dimensional FFT, the element corresponding to the index included in the selected reference table is extracted.
Multiply the extracted elements by the table coefficients contained in the selected reference table.
A two-dimensional Fourier is executed on the result of the multiplication,
The output of the two-dimensional Fourier is imaged.

一態様による合成開口ソーナーの信号処理方法は、
送受波器で受波された受波音波のデータに対して合成開口処理を行う合成開口ソーナーの信号処理方法であって、
参照テーブルを選択し、
前記データに対して、2次元FFTを実行し、
前記2次元FFTの出力から、前記選択した参照テーブルに含まれるインデックスに対応する要素を抽出し、
前記抽出した要素に対して、前記選択した参照テーブルに含まれるテーブル係数を乗算し、
前記乗算の結果に対して、2次元IFFTを実行し、
前記2次元IFFTの出力を画像化する。
The signal processing method of the synthetic aperture sonar according to one aspect is
It is a signal processing method of a synthetic aperture sonar that performs synthetic aperture processing on the data of the received sound waves received by the transmitter / receiver.
Select a lookup table and
A two-dimensional FFT is executed on the data,
From the output of the two-dimensional FFT, the element corresponding to the index included in the selected reference table is extracted.
Multiply the extracted elements by the table coefficients contained in the selected reference table.
A two-dimensional Fourier is executed on the result of the multiplication,
The output of the two-dimensional Fourier is imaged.

一態様による信号処理プログラムは、
コンピュータを、
送受波器で受波された受波音波のデータに対して合成開口処理を行う信号処理手段として機能させる信号処理プログラムであって、
前記信号処理手段は、
参照テーブルを選択し、
前記データに対して、2次元FFTを実行し、
前記2次元FFTの出力から、前記選択した参照テーブルに含まれるインデックスに対応する要素を抽出し、
前記抽出した要素に対して、前記選択した参照テーブルに含まれるテーブル係数を乗算し、
前記乗算の結果に対して、2次元IFFTを実行し、
前記2次元IFFTの出力を画像化する。
The signal processing program according to one aspect is
Computer,
A signal processing program that functions as a signal processing means that performs synthetic aperture processing on the received sound wave data received by the transmitter / receiver.
The signal processing means
Select a lookup table and
A two-dimensional FFT is executed on the data,
From the output of the two-dimensional FFT, the element corresponding to the index included in the selected reference table is extracted.
Multiply the extracted elements by the table coefficients contained in the selected reference table.
A two-dimensional Fourier is executed on the result of the multiplication,
The output of the two-dimensional Fourier is imaged.

上述の態様によれば、合成開口処理を行う信号処理器の処理負荷を低減し、処理速度を高速化することができる合成開口ソーナー及びその信号処理方法、並びに信号処理プログラムを提供できるという効果が得られる。 According to the above aspect, there is an effect that a synthetic aperture sonar capable of reducing the processing load of a signal processor performing synthetic aperture processing and increasing the processing speed, a signal processing method thereof, and a signal processing program can be provided. can get.

実施の形態に係る合成開口ソーナーの構成例を示す図である。It is a figure which shows the structural example of the synthetic opening sonar which concerns on embodiment. 実施の形態に係る参照テーブルの例を示す図である。It is a figure which shows the example of the reference table which concerns on embodiment. 関連技術に係る合成開口処理のフローの例を示すフロー図である。It is a flow figure which shows the example of the flow of the synthetic opening processing which concerns on a related technique. 図3のステップS104のレプリカ及び参照関数の乗算処理の例を示す図である。It is a figure which shows the example of the multiplication process of the replica and the reference function of step S104 of FIG. 図3のステップS105~S107のStolt補間処理の例を示す図である。It is a figure which shows the example of the total interpolation processing of steps S105 to S107 of FIG. 図3のステップS105のインデックスの計算処理の例を示す図である。It is a figure which shows the example of the calculation process of the index of step S105 of FIG. 図3のステップS105の補間係数の計算処理の例を示す図である。It is a figure which shows the example of the calculation process of the interpolation coefficient of step S105 of FIG. 図3のステップS106の要素の抽出処理の例を示す図である。It is a figure which shows the example of the extraction process of the element of step S106 of FIG. 図3のステップS107の補間係数の乗算処理の例を示す図である。It is a figure which shows the example of the multiplication process of the interpolation coefficient of step S107 of FIG. 図3のステップS105~S107のStolt補間処理に適用可能な最近傍法の例を示す図である。It is a figure which shows the example of the nearest neighbor method applicable to the Salt interpolation processing of steps S105 to S107 of FIG. 実施の形態に係る合成開口ソーナーによる合成開口処理のフローの例を示すフロー図である。It is a flow figure which shows the example of the flow of the synthetic opening processing by the synthetic opening sonar which concerns on embodiment. 図11のステップS203の参照テーブルの読込処理の例を示す図である。It is a figure which shows the example of the reading process of the reference table of step S203 of FIG. 図11のステップS205の要素の抽出処理の例を示す図である。It is a figure which shows the example of the extraction process of the element of step S205 of FIG. 図11のステップS206のテーブル係数の乗算処理の例を示す図である。It is a figure which shows the example of the multiplication process of the table coefficient of step S206 of FIG. 実施の形態に係る合成開口ソーナーによる参照テーブルの読込方法の例を示す図である。It is a figure which shows the example of the reading method of the reference table by the synthetic opening sonar which concerns on embodiment. 実施の形態を概念的に示した合成開口ソーナーの構成例を示す図である。It is a figure which shows the structural example of the synthetic opening sonar which conceptually showed embodiment.

以下、図面を参照して本開示の実施の形態について説明する。なお、以下の記載及び図面は、説明の明確化のため、適宜、省略及び簡略化がなされている。また、以下の各図面において、同一の要素には同一の符号が付されており、必要に応じて重複説明は省略されている。 Hereinafter, embodiments of the present disclosure will be described with reference to the drawings. The following descriptions and drawings have been omitted or simplified as appropriate for the purpose of clarifying the explanation. Further, in each of the following drawings, the same elements are designated by the same reference numerals, and duplicate explanations are omitted as necessary.

<実施の形態の構成>
まず、図1を参照して、本実施の形態に係る合成開口ソーナー1の構成例について説明する。本実施の形態に係る合成開口ソーナー1は、RAWデータを取得しながら航走する水中航走体の内部で、リアルタイムに合成開口処理を実行するものであるとする。また、本実施の形態に係る合成開口ソーナー1は、左舷と右舷とで送波波形を変えて音波を送受波することで、左右両舷のSAS画像を生成するものであるとする。また、本実施の形態に係る合成開口ソーナー1は、水中音速を計測するセンサ(後述の音速センサ13)を備え、RAWデータ取得時の水中音速を参照できるものであるとする。
<Structure of Embodiment>
First, a configuration example of the synthetic opening sonar 1 according to the present embodiment will be described with reference to FIG. It is assumed that the synthetic opening sonar 1 according to the present embodiment executes synthetic opening processing in real time inside an underwater vehicle that is navigating while acquiring RAW data. Further, it is assumed that the synthetic aperture sonar 1 according to the present embodiment generates SAS images on both the left and right sides by transmitting and receiving sound waves by changing the wave transmission waveform on the port side and the starboard side. Further, it is assumed that the synthetic opening sonar 1 according to the present embodiment includes a sensor for measuring the underwater sound velocity (sound velocity sensor 13 described later), and can refer to the underwater sound velocity at the time of acquiring RAW data.

図1に示されるように、本実施の形態に係る合成開口ソーナー1は、送受波器11L,11Rと、送受信器12と、音速センサ13と、信号処理器14と、外部記憶装置15と、を備えている。 As shown in FIG. 1, the synthetic aperture sonar 1 according to the present embodiment includes transmitters / receivers 11L and 11R, a transmitter / receiver 12, a sound velocity sensor 13, a signal processor 14, an external storage device 15, and the like. It has.

送受波器11L,11Rは、音波の送受波を行うものである。送受波器11Lは、水中航走体の左舷側に設けられ、送波波形1の音波を送波する。また、送受波器11Rは、水中航走体の右舷側に設けられ、送波波形2の音波を送波する。 The transmitter / receiver 11L and 11R transmit and receive sound waves. The transmitter / receiver 11L is provided on the left side of the underwater vehicle and transmits the sound wave of the transmission waveform 1. Further, the transmitter / receiver 11R is provided on the right side of the underwater vehicle and transmits the sound wave of the transmission waveform 2.

送受信器12は、左舷の送受波器11Lで受波された受波音波をRAWデータに変換する共に、右舷の送受波器11Rで受波された受波音波をRAWデータに変換し、左舷及び右舷のRAWデータを、それぞれ独立に信号処理器14に出力する。 The transmitter / receiver 12 converts the received sound wave received by the left-side transmitter / receiver 11L into RAW data, and also converts the received sound wave received by the right-side transmitter / receiver 11R into RAW data, and converts the received sound wave into RAW data. The right-side RAW data is independently output to the signal processor 14.

音速センサ13は、水中音速を計測するセンサであり、計測した水中音速を信号処理器14に出力する。 The sound velocity sensor 13 is a sensor that measures the underwater sound velocity, and outputs the measured underwater sound velocity to the signal processor 14.

信号処理器14は、CPU及びメモリからハードウェア構成され、そのメモリ上に、後述のアルゴリズムを有する信号処理プログラムがソフトウェア実装される。すなわち、信号処理器14の機能は、CPUが、メモリに記憶された信号処理プログラムを読み出し実行することにより、実現される。 The signal processor 14 is composed of a CPU and a memory in hardware, and a signal processing program having an algorithm described later is mounted on the memory by software. That is, the function of the signal processor 14 is realized by the CPU reading and executing the signal processing program stored in the memory.

信号処理器14は、HDD(Hard Disc Drive)等の外部記憶装置15に接続され、外部記憶装置15との間でファイルを入出力可能な構成となっている。この外部記憶装置15に予め生成した複数の参照テーブルをファイル形式で保存しておく。 The signal processor 14 is connected to an external storage device 15 such as an HDD (Hard Disc Drive), and has a configuration capable of inputting / outputting files to / from the external storage device 15. A plurality of reference tables generated in advance are stored in the external storage device 15 in a file format.

参照テーブルの内容は、FFT(Fast Fourier Transform)ポイント数、送波波形、水中音速、内挿方式などの処理パラメータに依存した内容となる。ただし、水中音速以外の処理パラメータは固定とすることで、参照テーブルの総数を節約できる。 The contents of the reference table depend on the processing parameters such as the number of FFT (Fast Fourier Transform) points, the wave waveform, the underwater sound velocity, and the interpolation method. However, by fixing the processing parameters other than the underwater sound velocity, the total number of reference tables can be saved.

水中音速は、最小値/最大値及びステップが定められ、信号処理器14は、音速センサ13にて計測された水中音速を、例えば、C(-M)、C(-M+1)、・・・、C(0)、C(1)、・・・、C(M)の合計2M+1個の離散値で近似する。参照テーブルは、水中音速以外の処理パラメータを固定する場合、合計2M+1個の水中音速毎に予め生成し、外部記憶装置15に保存しておく。本実施の形態では、左舷と右舷とで送波波形が異なるため、図2に示されるように、参照テーブルの総数は、2×(2M+1)個となる。 The minimum / maximum value and step of the underwater sound velocity are determined, and the signal processor 14 sets the underwater sound velocity measured by the sound velocity sensor 13, for example, C (-M), C (-M + 1), ... , C (0), C (1), ..., C (M), totaling 2M + 1 discrete values. When the processing parameters other than the underwater sound velocity are fixed, the reference table is generated in advance for each of a total of 2M + 1 underwater sound velocity and stored in the external storage device 15. In the present embodiment, since the wave transmission waveforms are different between the port side and the starboard side, the total number of reference tables is 2 × (2M + 1) as shown in FIG.

また、図2に示されるように、参照テーブルは、インデックス及びテーブル係数が含まれている。インデックス及びテーブル係数はいずれも、後述の2D-FFT(2次元FFT)のFFTポイント数と同じ要素数の行列形式となっている。 Also, as shown in FIG. 2, the reference table includes an index and table coefficients. Both the index and the table coefficient are in a matrix format having the same number of elements as the number of FFT points of the 2D-FFT (two-dimensional FFT) described later.

信号処理器14は、左舷及び右舷のRAWデータのそれぞれに対して、合成開口処理を行う。その際、信号処理器14は、合成開口処理毎に、外部記憶装置15に保存された複数の参照テーブルの中から、処理パラメータに応じた最適な参照テーブルを選択し、選択した参照テーブルをメモリに読込んで使用する。水中音速以外の処理パラメータを固定とする場合は、信号処理器14は、音速センサ13にて計測された水中音速を近似した水中音速Cに応じた参照テーブルを選択する。なお、本実施の形態に係る参照テーブルを使用した合成開口処理については後述する。 The signal processor 14 performs synthetic aperture processing on each of the port and starboard RAW data. At that time, the signal processor 14 selects the optimum reference table according to the processing parameters from the plurality of reference tables stored in the external storage device 15 for each synthetic aperture processing, and stores the selected reference table in the memory. Read and use in. When the processing parameters other than the underwater sound velocity are fixed, the signal processor 14 selects a reference table corresponding to the underwater sound velocity C that approximates the underwater sound velocity measured by the sound velocity sensor 13. The synthetic opening process using the reference table according to the present embodiment will be described later.

<実施の形態の動作>
以下、本実施の形態に係る合成開口ソーナー1による合成開口処理について説明する。ここでは、対比のために、関連技術に係る合成開口処理についても合わせて説明する。なお、以下では、Wavenumber Domainアルゴリズムによる合成開口処理を例に挙げて説明する。また、水中音速以外の処理パラメータは固定するものとする。
<Operation of the embodiment>
Hereinafter, the synthetic opening treatment by the synthetic opening sonar 1 according to the present embodiment will be described. Here, for comparison, the synthetic opening treatment according to the related technique will also be described. In the following, the synthetic aperture processing by the Wavenumber Domine algorithm will be described as an example. In addition, processing parameters other than the underwater sound velocity shall be fixed.

<関連技術に係る合成開口処理のフロー>
まず、図3を参照して、関連技術に係る合成開口処理のフローについて説明する。
図3に示されるように、まず、水中音速の平均を算出し、これを水中音速Cとする(ステップS101)。
続いて、レプリカR(i)及び参照関数φ(i,j)を計算する(ステップS102)。ここで、レプリカR(i)は、送波波形をFFTしたスペクトルである。また、参照関数φ(i,j)は、r-u領域を2D-FFTしたkr(i)-ku(j)領域での空間チャープ参照関数である。また、rはレンジ方向距離、kr(i)は、レンジ方向波数であり、i≦Nである(N=レンジ方向FFTポイント数)。また、uはアジマス方向距離、ku(j)は、アジマス方向波数であり、j≦Nである(N=アジマス方向FFTポイント数)。
<Flow of synthetic aperture processing related to related technology>
First, with reference to FIG. 3, the flow of synthetic aperture processing according to the related technique will be described.
As shown in FIG. 3, first, the average of the underwater sound velocity is calculated, and this is defined as the underwater sound velocity C (step S101).
Subsequently, the replica R (i) and the reference function φ (i, j) are calculated (step S102). Here, the replica R (i) is a spectrum obtained by FFTing the transmitted waveform. Further, the reference function φ (i, j) is a spatial chirp reference function in the kr (i) -ku (j) region in which the ru-u region is 2D-FFT. Further, r is the distance in the range direction, kr (i) is the wave number in the range direction, and i ≦ NR ( NR = number of FFT points in the range direction). Further, u is the distance in the azimuth direction, and ku (j) is the wave number in the azimuth direction, and j ≦ NA ( NA = the number of FFT points in the azimuth direction).

続いて、RAWデータを、2D-FFTして、r-u領域からkr(i)-ku(j)領域に変換する(ステップS103)。RAWデータの2D-FFTにより、図4に示されるように、kr(i)-ku(j)座標上のスペクトルA(i,j)が得られる。 Subsequently, the RAW data is 2D-FFT and converted from the ru-u region to the kr (i) -ku (j) region (step S103). The 2D-FFT of the RAW data gives the spectrum A0 (i, j) on the kr (i) -ku (j) coordinates, as shown in FIG.

続いて、スペクトルA(i,j)にレプリカR(i)を乗算し、さらに、参照関数φ(i,j)を乗算する(ステップS104)。レプリカR(i)及び参照関数φ(i,j)の乗算により、最終的に、図4に示されるように、スペクトルA(i,j)が得られる。 Subsequently, the spectrum A 0 (i, j) is multiplied by the replica R (i), and further multiplied by the reference function φ (i, j) (step S104). By multiplying the replica R (i) and the reference function φ (i, j), the spectrum A (i, j) is finally obtained, as shown in FIG.

続いて、不等間隔リサンプリング処理であるStolt補間を行う。Stolt補間は、図5に示されるように、座標変換[kr(i),ku(j)]→[kx(i),ky(j)]を適用した後のスペクトルB(i,j)を求める処理である。 Subsequently, Retro interpolation, which is an unequal-interval resampling process, is performed. As shown in FIG. 5, the Salt interpolation applies the spectrum B (i, j) after applying the coordinate transformation [kr (i), ku (j)] → [kx (i), ky (j)]. This is the desired process.

kr(i)はkx座標のグリッド点上に写像されるとは限らないため、kx座標のグリッド点上のスペクトルB(i,j)は、補間により求める必要がある。この補間方法はいくつかの例があるが、以下では、最も一般的な線形補間法を例に挙げて説明する。 Since kr (i) is not always mapped on the grid point of the kx coordinate, the spectrum B (i, j) on the grid point of the kx coordinate needs to be obtained by interpolation. There are several examples of this interpolation method, but the most common linear interpolation method will be described below as an example.

まず、インデックス及び補間係数αを計算する(ステップS105)。具体的には、図6に示されるように、水中音速Cに基づく写像[kr,ku]→[kx’,ky’]により、kx’(i)インデックス(i=1~N)を得る。続いて、図7に示されるように、kx(i)の前後のkx’を抽出し(この例では、kx’(i),kx’(i+1))、補間係数αを計算する。 First, the index and the interpolation coefficient α are calculated (step S105). Specifically, as shown in FIG. 6, the kx'(i) index (i = 1 to NR ) is obtained by the mapping [kr, ku] → [kx', ky'] based on the underwater sound velocity C. .. Subsequently, as shown in FIG. 7, kx'before and after kx (i) is extracted (in this example, kx'(i), kx'(i + 1)), and the interpolation coefficient α is calculated.

続いて、図8に示されるように、スペクトルA(i,j)から、kx’インデックスに対応する要素(この例では、A(i,j),A(i+1,j))を抽出する(ステップS106)。
続いて、図9に示されるように、補間係数αと要素A(i,j),A(i+1,j)との乗算を行い、以下のスペクトルB(i,j)を得る(ステップS107)。
B(i,j)=α・A(i,j)+(1-α)・A(i+1,j)
以上でStolt補間処理が終了する。
Subsequently, as shown in FIG. 8, the elements corresponding to the kx'index (A (i, j), A (i + 1, j) in this example) are extracted from the spectrum A (i, j) (in this example, A (i, j), A (i + 1, j)) ( Step S106).
Subsequently, as shown in FIG. 9, the interpolation coefficient α is multiplied by the elements A (i, j) and A (i + 1, j) to obtain the following spectrum B (i, j) (step S107). ..
B (i, j) = α · A (i, j) + (1-α) · A (i + 1, j)
This completes the Salt interpolation process.

続いて、スペクトルB(i,j)を、2D-IFFT(Inverse Fast Fourier Transform)して、kx(i)-ky(j)領域からx-y領域に変換する(ステップS108)。
その後、2D-IFFTの出力を、SAS画像に画像化する(ステップS109)。
Subsequently, the spectrum B (i, j) is subjected to 2D-IFFT (Inverse Fast Fourier Transform) to transform the kx (i) -ky (j) region into the xy region (step S108).
After that, the output of the 2D-IFFT is imaged into a SAS image (step S109).

以上の同じ処理を、左舷のRAWデータと右舷のRAWデータについて独立に行う。ただし、本実施の形態のように、左舷と右舷とで送波波形が異なる場合は、左舷と右舷とで処理パラメータは異なる。具体的には、上記のレプリカが左舷と右舷とで異なる。 The same process as described above is performed independently for the port RAW data and the starboard RAW data. However, when the wave transmission waveform is different between the port side and the starboard side as in the present embodiment, the processing parameters are different between the port side and the starboard side. Specifically, the above replica differs between port and starboard.

なお、Stolt補間処理に適用される、線形補間法以外の補間方法としては、最近傍法及びN次補間法等が挙げられる。
最近傍法は、図10に示されるように、kxに最も近いkx’に対応する要素で置き換えるものである。この例では、kx(i)に最も近いのはkx’(i+1)なので、B(i,j)=A(i+1,j)となる。
Examples of the interpolation method other than the linear interpolation method applied to the halt interpolation processing include the nearest neighbor method and the Nth-order interpolation method.
The nearest neighbor method replaces with the element corresponding to kx'closest to kx, as shown in FIG. In this example, kx'(i + 1) is the closest to kx (i), so B (i, j) = A (i + 1, j).

一方、N次補間法は、kx(i)を区間内に含むN+1点のkx’及びAを使ってB(i,j)を計算するものである。なお、N次補間法の図示及び具体的な計算方法の説明は省略する。 On the other hand, in the Nth-order interpolation method, B (i, j) is calculated using kx'and A of N + 1 points including kx (i) in the interval. The illustration of the Nth-order interpolation method and the description of a specific calculation method will be omitted.

このように、関連技術に係る合成開口処理では、Stolt補間などの不等間隔リサンプリング処理を計算で行っているため、計算量が比較的多かった。そのため、合成開口処理を行う信号処理器の処理負荷が大きく、処理時間を要し、リアルタイム化の阻害要因となっていた。 As described above, in the synthetic aperture processing according to the related technique, the amount of calculation is relatively large because the non-equidistant resampling processing such as the Salt interpolation is performed by calculation. Therefore, the processing load of the signal processor that performs the synthetic aperture processing is large, the processing time is required, and it is an obstacle to real-time operation.

<実施の形態に係る合成開口処理のフロー>
続いて、図11を参照して、本実施の形態に係る合成開口ソーナー1による合成開口処理のフローについて説明する。
図11に示されるように、まず、信号処理器14は、音速センサ13にて計測された水中音速の平均を算出し、算出した水中音速の平均を、上述の合計2M+1個の離散値のいずれかで近似し、これを水中音速Cとする(ステップS201)。
<Flow of synthetic aperture processing according to the embodiment>
Subsequently, with reference to FIG. 11, the flow of the synthetic opening processing by the synthetic opening sonar 1 according to the present embodiment will be described.
As shown in FIG. 11, first, the signal processor 14 calculates the average of the underwater sound velocity measured by the sound velocity sensor 13, and the calculated average of the underwater sound velocity is one of the above-mentioned total 2M + 1 discrete values. Approximate with, and let this be the underwater sound velocity C (step S201).

続いて、信号処理器14は、メモリに現在読込んでいる参照テーブルを変更するか否かを判断する(ステップS202)。ここでは、水中音速Cに応じた参照テーブルをメモリに現在読込んでいなければ、参照テーブルを変更すると判断する。 Subsequently, the signal processor 14 determines whether or not to change the reference table currently being read into the memory (step S202). Here, if the reference table corresponding to the underwater sound velocity C is not currently read into the memory, it is determined that the reference table is changed.

信号処理器14は、参照テーブルを変更する場合(ステップS202のYes)、外部記憶装置15から新たに参照テーブルを読込む(ステップS203)。なお、本実施の形態に係る参照テーブルの読込方法の詳細については後述する。 When the signal processor 14 changes the reference table (Yes in step S202), the signal processor 14 newly reads the reference table from the external storage device 15 (step S203). The details of the reference table reading method according to the present embodiment will be described later.

ここで、水中音速Cに応じた参照テーブルには、水中音速Cに応じたkx’インデックスと、水中音速Cに応じた、補間係数に替わるテーブル係数T(i,j)と、が含まれている。そのため、以降に、kx’インデックスに対応する要素を抽出すれば、補間計算が可能となる。 Here, the reference table corresponding to the underwater sound velocity C includes a kx'index corresponding to the underwater sound velocity C and a table coefficient T (i, j) corresponding to the underwater sound velocity C instead of the interpolation coefficient. There is. Therefore, if the element corresponding to the kx'index is extracted thereafter, the interpolation calculation can be performed.

直線補間の場合、1点を補間計算に必要なサンプルは2個である。よって、図12に示されるように、テーブル係数T(i,j)も2個必要である(この例では、T(i,j)、T(i,j)のペアでB(i,j)を計算する)。 In the case of linear interpolation, two samples are required for one point interpolation calculation. Therefore, as shown in FIG. 12, two table coefficients T (i, j) are also required (in this example, a pair of T 1 (i, j) and T 2 (i, j) is B (i). , J) is calculated).

テーブル係数T(i,j)は、2D-FFTの出力であるA(i,j)からB(i,j)を算出するための係数であり、関連技術に係るフローに基づき、以下の定義で予め計算し保存しておく。
(i,j)=α・A(i,j)/A(i,j)=α・R(i)・φ(i,j)
(i,j)=(1-α)・A(i+1,j)/A(i+1,j)
=(1-α)・R(i+1)・φ(i+1,j)
このように、テーブル係数T(i,j)は、補間係数αに対し、使用するレプリカR(i)及び参照関数φ(i,j)を予め乗算した形式の係数となる。
The table coefficient T (i, j) is a coefficient for calculating B (i, j) from A 0 (i, j), which is the output of the 2D-FFT, and is as follows based on the flow related to the related technology. Calculate in advance by definition and save.
T 1 (i, j) = α · A (i, j) / A 0 (i, j) = α · R (i) · φ (i, j)
T 2 (i, j) = (1-α) · A (i + 1, j) / A 0 (i + 1, j)
= (1-α) ・ R (i + 1) ・ φ (i + 1, j)
As described above, the table coefficient T (i, j) is a coefficient in the form of the interpolation coefficient α multiplied in advance by the replica R (i) to be used and the reference function φ (i, j).

続いて、信号処理器14は、送受信器12から得られたRAWデータを、2D-FFTして、r-u領域からkr(i)-ku(j)領域に変換する(ステップS204)。 Subsequently, the signal processor 14 performs 2D-FFT of the RAW data obtained from the transmitter / receiver 12 to convert the ru region into the kr (i) -ku (j) region (step S204).

続いて、信号処理器14は、関連技術に係るレプリカ及び参照関数の乗算処理(ステップS104)及びStolt補間処理(ステップS105~S107)に相当する処理を行い、kx(i),ky(j)のスペクトルB(i,j)を求める。以下では、補間方法として、関連技術に係るフローと同様に、線形補間法を例に挙げて説明する。 Subsequently, the signal processor 14 performs a process corresponding to the replica and reference function multiplication process (step S104) and the Salt interpolation process (steps S105 to S107) related to the related technique, and kx (i), ky (j). The spectrum B (i, j) of In the following, as the interpolation method, a linear interpolation method will be described as an example, as in the flow related to the related technique.

まず、信号処理器14は、図13に示されるように、2D-FFTの出力であるA(i,j)から、kx’インデックスに対応する要素(この例では、A(i,j),A(i+1,j))を抽出する(ステップS205)。
続いて、信号処理器14は、図14に示されるように、テーブル係数T(i,j),T(i,j)と要素A(i,j),A(i+1,j)との乗算を行い、以下のスペクトルB(i,j)を得る(ステップS206)。
B(i,j)=T(i,j)・A(i,j)+T(i,j)・A(i+1,j)
以上で、関連技術に係るレプリカ及び参照関数の乗算処理及びStolt補間処理に相当する処理が終了する。
First, as shown in FIG. 13, the signal processor 14 has elements corresponding to the kx'index from A 0 (i, j), which is the output of the 2D-FFT (in this example, A 0 (i, j). ), A 0 (i + 1, j)) (step S205).
Subsequently, as shown in FIG. 14, the signal processor 14 has table coefficients T 1 (i, j), T 2 (i, j) and elements A 0 (i, j), A 0 (i + 1, j). ) To obtain the following spectrum B (i, j) (step S206).
B (i, j) = T 1 (i, j) · A 0 (i, j) + T 2 (i, j) · A 0 (i + 1, j)
This completes the processing corresponding to the replica and reference function multiplication processing and the Salt interpolation processing related to the related techniques.

続いて、信号処理器14は、スペクトルB(i,j)を、2D-IFFTして、kx(i)-ky(j)領域からx-y領域に変換する(ステップS207)。
その後、信号処理器14は、2D-IFFTの出力を、SAS画像に画像化し(ステップS208)、SAS画像を外部記憶装置15に保存する。
Subsequently, the signal processor 14 performs 2D-IFFT of the spectrum B (i, j) to convert the kx (i) -ky (j) region into the xy region (step S207).
After that, the signal processor 14 images the output of the 2D-IFFT into a SAS image (step S208), and stores the SAS image in the external storage device 15.

以上の同じ処理を、左舷のRAWデータと右舷のRAWデータについて独立に行う。ただし、本実施の形態では、左舷と右舷とで送波波形が異なるため、左舷と右舷とで処理パラメータは異なる。具体的には、上記のレプリカが左舷と右舷とで異なる。 The same process as described above is performed independently for the port RAW data and the starboard RAW data. However, in the present embodiment, since the wave transmission waveform is different between the port side and the starboard side, the processing parameters are different between the port side and the starboard side. Specifically, the above replica differs between port and starboard.

<実施の形態に係る参照テーブルの読込方法>
続いて、図15を参照して、本実施の形態に係る合成開口ソーナー1による参照テーブルの読込方法について説明する。
図15に示されるように、水中音速Cは、時間経過に伴い変化する場合がある(図中の折れ線を参照)。そこで、信号処理器14は、水中音速Cの変化に対応するため、現在の水中音速Cに応じて選択した参照テーブルだけでなく、その参照テーブルの前及び/又は後の参照テーブルも、外部記憶装置15から予めメモリに読込み、保持しておく。
<Method of reading the reference table according to the embodiment>
Subsequently, with reference to FIG. 15, a method of reading the reference table by the synthetic aperture sonar 1 according to the present embodiment will be described.
As shown in FIG. 15, the underwater sound velocity C may change over time (see the polygonal line in the figure). Therefore, in order to respond to the change in the underwater sound velocity C, the signal processor 14 externally stores not only the reference table selected according to the current underwater sound velocity C but also the reference tables before and / or after the reference table. It is read from the device 15 into the memory in advance and stored.

この例は、3つの参照テーブルをメモリに常に保持する例である。初回は、現在の水中音速Cに応じて選択した参照テーブル及びその前後の参照テーブルをメモリに読込んで保持し、選択した参照テーブルを処理に適用する。次回以降は、現在の水中音速が、参照テーブルを保持している水中音速の範囲内であれば、参照テーブルの新規の読込みは行わず(処理に適用する参照テーブルは適宜変更)、範囲外となったら、現在の水中音速に応じて選択した参照テーブルを新規にメモリに読込み、処理に適用する。 This example is an example of always holding three reference tables in memory. At the first time, the reference table selected according to the current underwater sound velocity C and the reference tables before and after it are read and held in the memory, and the selected reference table is applied to the processing. From the next time onward, if the current underwater sound velocity is within the range of the underwater sound velocity that holds the reference table, the reference table will not be newly read (the reference table applied to the processing will be changed as appropriate), and it will be out of range. When it becomes, the reference table selected according to the current underwater sound velocity is newly read into the memory and applied to the processing.

具体的には、この例では、初回(1回目)の水中音速はC(0)である。そのため、信号処理器14は、C(0)及びその前後のC(-1),C(1)の参照テーブルをメモリに読込み、C(0)の参照テーブルを処理に適用する。
続く2回目の水中音速もC(0)である。そのため、信号処理器14は、参照テーブルの新規の読込みは行わず、C(0)の参照テーブルを継続して処理に適用する。
Specifically, in this example, the first (first) underwater sound velocity is C (0). Therefore, the signal processor 14 reads the reference tables of C (0) and C (-1) and C (1) before and after it into the memory, and applies the reference table of C (0) to the processing.
The second underwater sound velocity that follows is also C (0). Therefore, the signal processor 14 does not read a new reference table, and continuously applies the reference table of C (0) to the processing.

続く3回目の水中音速はC(-1)に変化しているが、C(-1)は、参照テーブルを保持しているC(-1)~C(1)の範囲内である。そのため、信号処理器14は、参照テーブルの新規の読込みは行わず、処理に適用する参照テーブルをC(-1)の参照テーブルに変更する。
続く4回目の水中音速もC(-1)である。そのため、信号処理器14は、参照テーブルの新規の読込みは行わず、C(-1)の参照テーブルを継続して処理に適用する。
The subsequent third underwater sound velocity changes to C (-1), which is in the range of C (-1) to C (1) holding the reference table. Therefore, the signal processor 14 does not read a new reference table, but changes the reference table applied to the processing to the reference table of C (-1).
The fourth underwater sound velocity that follows is also C (-1). Therefore, the signal processor 14 does not read a new reference table, and continuously applies the reference table of C (-1) to the processing.

続く5回目の水中音速はC(-2)に変化しており、C(-2)は、参照テーブルを保持しているC(-1)~C(1)の範囲外である。そのため、信号処理器14は、C(-2)の参照テーブルを新規にメモリに読込み、C(1)の参照テーブルをメモリから削除する。そして、信号処理器14は、処理に適用する参照テーブルをC(-2)の参照テーブルに変更する。
6回目以降の処理についての説明は省略する。
The subsequent fifth underwater sound velocity changes to C (-2), which is outside the range of C (-1) to C (1) holding the reference table. Therefore, the signal processor 14 newly reads the reference table of C (-2) into the memory, and deletes the reference table of C (1) from the memory. Then, the signal processor 14 changes the reference table applied to the processing to the reference table of C (-2).
The description of the sixth and subsequent processes will be omitted.

上述したように、本実施の形態によれば、信号処理器14は、参照テーブルを選択し、RAWデータに対して、2D-FFTを実行し、2D-FFTの出力から、参照テーブルに含まれるインデックスに対応する要素を抽出し、抽出した要素に対して、参照テーブルに含まれるテーブル係数を乗算し、乗算の結果に対して、2D-IFFTを実行し、2D-IFFTの出力を画像化する。 As described above, according to the present embodiment, the signal processor 14 selects the reference table, executes the 2D-FFT on the RAW data, and is included in the reference table from the output of the 2D-FFT. The element corresponding to the index is extracted, the extracted element is multiplied by the table coefficient included in the reference table, the 2D-IFFT is executed on the result of the multiplication, and the output of the 2D-IFFT is imaged. ..

このように、本実施の形態によれば、関連技術では処理負荷の大きかった、Stolt補間などの不等間隔リサンプリング処理を、計算ではなく、参照テーブルの参照に置き換えている。これにより、合成開口処理の計算量が減るため、信号処理器14の処理負荷を低減し、処理速度を高速化することができる。その結果、リアルタイム化に寄与し得る。また、信号処理器14に高速・高性能なCPUを採用する必要がないため、信号処理器14の低消費電力化が可能となる。 As described above, according to the present embodiment, the non-equidistant resampling process such as the Salt interpolation, which has a large processing load in the related technique, is replaced with the reference of the reference table instead of the calculation. As a result, the amount of calculation for the synthetic aperture processing is reduced, so that the processing load of the signal processor 14 can be reduced and the processing speed can be increased. As a result, it can contribute to real-time operation. Further, since it is not necessary to adopt a high-speed and high-performance CPU for the signal processor 14, the power consumption of the signal processor 14 can be reduced.

また、本実施の形態によれば、参照テーブルとして、処理パラメータに基づき生成した複数の参照テーブルを予め生成し、外部記憶装置15に保存しておき、信号処理器14は、合成開口処理毎に、処理パラメータに応じた最適な参照テーブルを選択し、メモリに読み込んで使用する。これにより、信号処理器14の処理速度と共に画像性能も合わせて担保することが可能となる。 Further, according to the present embodiment, as the reference table, a plurality of reference tables generated based on the processing parameters are generated in advance and stored in the external storage device 15, and the signal processor 14 performs each synthetic aperture processing. , Select the most suitable reference table according to the processing parameters, read it into the memory, and use it. As a result, it is possible to guarantee the image performance as well as the processing speed of the signal processor 14.

また、本実施の形態によれば、水中音速以外の処理パラメータは固定とすることで、参照テーブルの総数を節約できる。
また、本実施の形態によれば、信号処理器14は、水中音速に応じて選択した参照テーブルと共に、その参照テーブルの前及び/又は後の参照テーブルも、外部記憶装置15からメモリに読込んで保持しておく。これにより、水中音速が時間経過に伴い変化する場合においても、参照テーブルを外部記憶装置15から読込む頻度や負荷を低減することが可能となる。
Further, according to the present embodiment, the total number of reference tables can be saved by fixing the processing parameters other than the underwater sound velocity.
Further, according to the present embodiment, the signal processor 14 reads the reference table before and / or after the reference table selected according to the underwater sound velocity into the memory from the external storage device 15. Keep it. This makes it possible to reduce the frequency and load of reading the reference table from the external storage device 15 even when the underwater sound velocity changes with the passage of time.

<実施の形態の概念>
続いて、図16を参照して、上述の実施の形態を概念的に示した合成開口ソーナー2の構成例について説明する。図16に示されるように、合成開口ソーナー2は、送受波器101及び信号処理器102を備えている。
<Concept of embodiment>
Subsequently, a configuration example of the synthetic aperture sonar 2 conceptually showing the above-described embodiment will be described with reference to FIG. As shown in FIG. 16, the synthetic aperture sonar 2 includes a transmitter / receiver 101 and a signal processor 102.

送受波器101は、音波の送受波を行う。送受波器101は、送受波器11L,11Rに対応する。
信号処理器102は、送受波器101で受波された受波音波のデータに対して合成開口処理を行う。信号処理器102は、信号処理器14に対応する。そのため、信号処理器102は、信号処理器14による合成開口処理と同様の合成開口処理を行う。
The transmitter / receiver 101 transmits / receives sound waves. The transmitter / receiver 101 corresponds to the transmitter / receiver 11L, 11R.
The signal processor 102 performs a synthetic aperture process on the data of the received sound wave received by the transmitter / receiver 101. The signal processor 102 corresponds to the signal processor 14. Therefore, the signal processor 102 performs the same synthetic aperture processing as the synthetic aperture processing by the signal processor 14.

すなわち、信号処理器102は、参照テーブルを選択し、送受波器101で受波された受波音波のデータに対して、2次元FFTを実行し、2次元FFTの出力から、選択した参照テーブルに含まれるインデックスに対応する要素を抽出し、抽出した要素に対して、選択した参照テーブルに含まれるテーブル係数を乗算し、乗算の結果に対して、2次元IFFTを実行し、2次元IFFTの出力を画像化する。 That is, the signal processor 102 selects the reference table, executes the two-dimensional FFT on the data of the received sound wave received by the transmitter / receiver 101, and selects the reference table from the output of the two-dimensional FFT. The element corresponding to the index included in is extracted, the extracted element is multiplied by the table coefficient contained in the selected reference table, the result of the multiplication is subjected to the two-dimensional Fourier, and the two-dimensional Fourier is executed. Image the output.

以上、実施の形態を参照して本開示を説明したが、本開示は上記の実施の形態に限定されるものではない。本開示の構成や詳細には、本開示のスコープ内で当業者が理解し得る様々な変更をすることができる。 Although the present disclosure has been described above with reference to the embodiments, the present disclosure is not limited to the above embodiments. Various changes that can be understood by those skilled in the art can be made to the structure and details of the present disclosure within the scope of the present disclosure.

例えば、上記の信号処理プログラムは、様々なタイプの非一時的なコンピュータ可読媒体(non-transitory computer readable medium)を用いて格納され、コンピュータに供給することができる。非一時的なコンピュータ可読媒体は、様々なタイプの実体のある記録媒体(tangible storage medium)を含む。非一時的なコンピュータ可読媒体の例は、磁気記録媒体(例えばフレキシブルディスク、磁気テープ、ハードディスクドライブ)、光磁気記録媒体(例えば光磁気ディスク)、CD-ROM(Compact Disc-Read Only Memory)、CD-R(CD-Recordable)、CD-R/W(CD-ReWritable)、半導体メモリ(例えば、マスクROM、PROM(Programmable ROM)、EPROM(Erasable PROM)、フラッシュROM、RAM(Random Access Memory))を含む。また、上記の信号処理プログラムは、様々なタイプの一時的なコンピュータ可読媒体(transitory computer readable medium)によってコンピュータに供給されても良い。一時的なコンピュータ可読媒体の例は、電気信号、光信号、及び電磁波を含む。一時的なコンピュータ可読媒体は、電線及び光ファイバ等の有線通信路、又は無線通信路を介して、プログラムをコンピュータに供給できる。 For example, the signal processing program described above can be stored and supplied to a computer using various types of non-transitory computer readable media. Non-transitory computer-readable media include various types of tangible storage media. Examples of non-temporary computer-readable media include magnetic recording media (eg flexible disks, magnetic tapes, hard disk drives), photomagnetic recording media (eg optomagnetic discs), CD-ROMs (Compact Disc-Read Only Memory), CDs. -R (CD-Recordable), CD-R / W (CD-ReWritable), semiconductor memory (for example, mask ROM, PROM (Programmable ROM), EPROM (Erasable PROM), flash ROM, RAM (Random Access Memory)) include. The signal processing program described above may also be supplied to the computer by various types of transient computer readable media. Examples of temporary computer-readable media include electrical, optical, and electromagnetic waves. The temporary computer-readable medium can supply the program to the computer via a wired communication path such as an electric wire and an optical fiber, or a wireless communication path.

上記の実施の形態の一部又は全部は、以下の付記のようにも記載されうるが、以下には限られない。
(付記1)
音波の送受波を行う送受波器と、
前記送受波器で受波された受波音波のデータに対して合成開口処理を行う信号処理器と、を備え、
前記信号処理器は、
参照テーブルを選択し、
前記データに対して、2次元FFTを実行し、
前記2次元FFTの出力から、前記選択した参照テーブルに含まれるインデックスに対応する要素を抽出し、
前記抽出した要素に対して、前記選択した参照テーブルに含まれるテーブル係数を乗算し、
前記乗算の結果に対して、2次元IFFTを実行し、
前記2次元IFFTの出力を画像化する、
合成開口ソーナー。
(付記2)
前記信号処理器は、
メモリを備え、
複数の前記参照テーブルが保存された外部記憶装置に接続され、
前記外部記憶装置に保存された複数の前記参照テーブルの中から参照テーブルを選択し、前記選択した参照テーブルを前記メモリに読込んで使用する、
付記1に記載の合成開口ソーナー。
(付記3)
前記外部記憶装置には、処理パラメータ毎に生成された複数の前記参照テーブルが保存されており、
前記信号処理器は、前記処理パラメータに応じて、前記外部記憶装置に保存された複数の前記参照テーブルの中から参照テーブルを選択する、
付記2に記載の合成開口ソーナー。
(付記4)
水中音速を計測する音速センサを備え、
前記処理パラメータは、水中音速以外は固定とされ、
前記外部記憶装置には、前記水中音速毎に生成された複数の前記参照テーブルが保存されており、
前記信号処理器は、前記音速センサにて計測された水中音速に応じて、前記外部記憶装置に保存された複数の前記参照テーブルの中から参照テーブルを選択する、
付記3に記載の合成開口ソーナー。
(付記5)
前記信号処理器は、前記音速センサにて計測された水中音速に応じて選択した参照テーブルと共に、当該参照テーブルの前及び/又は後の参照テーブルを、前記メモリに読込んで保持する、
付記4に記載の合成開口ソーナー。
(付記6)
前記テーブル係数は、補間係数に対し、前記合成開口処理に使用するレプリカ及び参照関数を予め乗算した形式の係数である、
付記1から5のいずれか1項に記載の合成開口ソーナー。
(付記7)
前記参照テーブルに含まれるインデックス及びテーブル係数は、前記2次元FFTのFFTポイント数と同じ要素数の行列形式である、
付記1から6のいずれか1項に記載の合成開口ソーナー。
(付記8)
前記送受波器は、水中航走体の左右両舷にそれぞれ設けられ、
前記信号処理器は、左右両舷の前記送受波器で受波された受波音波のデータのそれぞれに対して、前記合成開口処理を行う、
付記1から7のいずれか1項に記載の合成開口ソーナー。
(付記9)
送受波器で受波された受波音波のデータに対して合成開口処理を行う合成開口ソーナーの信号処理方法であって、
参照テーブルを選択し、
前記データに対して、2次元FFTを実行し、
前記2次元FFTの出力から、前記選択した参照テーブルに含まれるインデックスに対応する要素を抽出し、
前記抽出した要素に対して、前記選択した参照テーブルに含まれるテーブル係数を乗算し、
前記乗算の結果に対して、2次元IFFTを実行し、
前記2次元IFFTの出力を画像化する、
合成開口ソーナーの信号処理方法。
(付記10)
外部記憶装置に保存された複数の前記参照テーブルの中から参照テーブルを選択し、前記選択した参照テーブルをメモリに読込んで使用する、
付記9に記載の合成開口ソーナーの信号処理方法。
(付記11)
前記外部記憶装置には、処理パラメータ毎に生成された複数の前記参照テーブルが保存されており、
前記処理パラメータに応じて、前記外部記憶装置に保存された複数の前記参照テーブルの中から参照テーブルを選択する、
付記10に記載の合成開口ソーナーの信号処理方法。
(付記12)
前記処理パラメータは、水中音速以外は固定とされ、
前記外部記憶装置には、前記水中音速毎に生成された複数の前記参照テーブルが保存されており、
音速センサにて計測された水中音速に応じて、前記外部記憶装置に保存された複数の前記参照テーブルの中から参照テーブルを選択する、
付記11に記載の合成開口ソーナーの信号処理方法。
(付記13)
前記音速センサにて計測された水中音速に応じて選択した参照テーブルと共に、当該参照テーブルの前及び/又は後の参照テーブルを、前記メモリに読込んで保持する、
付記12に記載の合成開口ソーナーの信号処理方法。
(付記14)
前記テーブル係数は、補間係数に対し、前記合成開口処理に使用するレプリカ及び参照関数を予め乗算した形式の係数である、
付記9から13のいずれか1項に記載の合成開口ソーナーの信号処理方法。
(付記15)
前記参照テーブルに含まれるインデックス及びテーブル係数は、前記2次元FFTのFFTポイント数と同じ要素数の行列形式である、
付記9から14のいずれか1項に記載の合成開口ソーナーの信号処理方法。
(付記16)
前記送受波器は、水中航走体の左右両舷にそれぞれ設けられており、
左右両舷の前記送受波器で受波された受波音波のデータのそれぞれに対して、前記合成開口処理を行う、
付記9から15のいずれか1項に記載の合成開口ソーナーの信号処理方法。
(付記17)
コンピュータを、
送受波器で受波された受波音波のデータに対して合成開口処理を行う信号処理手段として機能させる信号処理プログラムであって、
前記信号処理手段は、
参照テーブルを選択し、
前記データに対して、2次元FFTを実行し、
前記2次元FFTの出力から、前記選択した参照テーブルに含まれるインデックスに対応する要素を抽出し、
前記抽出した要素に対して、前記選択した参照テーブルに含まれるテーブル係数を乗算し、
前記乗算の結果に対して、2次元IFFTを実行し、
前記2次元IFFTの出力を画像化する、
信号処理プログラム。
(付記18)
前記信号処理手段は、外部記憶装置に保存された複数の前記参照テーブルの中から参照テーブルを選択し、前記選択した参照テーブルをメモリに読込んで使用する、
付記17に記載の信号処理プログラム。
(付記19)
前記外部記憶装置には、処理パラメータ毎に生成された複数の前記参照テーブルが保存されており、
前記信号処理手段は、前記処理パラメータに応じて、前記外部記憶装置に保存された複数の前記参照テーブルの中から参照テーブルを選択する、
付記18に記載の信号処理プログラム。
(付記20)
前記処理パラメータは、水中音速以外は固定とされ、
前記外部記憶装置には、前記水中音速毎に生成された複数の前記参照テーブルが保存されており、
前記信号処理手段は、音速センサにて計測された水中音速に応じて、前記外部記憶装置に保存された複数の前記参照テーブルの中から参照テーブルを選択する、
付記19に記載の信号処理プログラム。
(付記21)
前記信号処理手段は、前記音速センサにて計測された水中音速に応じて選択した参照テーブルと共に、当該参照テーブルの前及び/又は後の参照テーブルを、前記メモリに読込んで保持する、
付記20に記載の信号処理プログラム。
(付記22)
前記テーブル係数は、補間係数に対し、前記合成開口処理に使用するレプリカ及び参照関数を予め乗算した形式の係数である、
付記17から21のいずれか1項に記載の信号処理プログラム。
(付記23)
前記参照テーブルに含まれるインデックス及びテーブル係数は、前記2次元FFTのFFTポイント数と同じ要素数の行列形式である、
付記17から22のいずれか1項に記載の信号処理プログラム。
(付記24)
前記送受波器は、水中航走体の左右両舷にそれぞれ設けられており、
前記信号処理手段は、左右両舷の前記送受波器で受波された受波音波のデータのそれぞれに対して、前記合成開口処理を行う、
付記17から23のいずれか1項に記載の信号処理プログラム。
Some or all of the above embodiments may also be described, but not limited to:
(Appendix 1)
A transmitter / receiver that transmits / receives sound waves,
A signal processor that performs synthetic aperture processing on the data of the received sound wave received by the transmitter / receiver is provided.
The signal processor
Select a lookup table and
A two-dimensional FFT is executed on the data,
From the output of the two-dimensional FFT, the element corresponding to the index included in the selected reference table is extracted.
Multiply the extracted elements by the table coefficients contained in the selected reference table.
A two-dimensional Fourier is executed on the result of the multiplication,
Imaging the output of the two-dimensional Fourier,
Synthetic opening sonar.
(Appendix 2)
The signal processor
Equipped with memory
A plurality of the reference tables are connected to the stored external storage device, and the reference table is connected to the storage device.
A reference table is selected from the plurality of reference tables stored in the external storage device, and the selected reference table is read into the memory and used.
The synthetic opening sonar according to Appendix 1.
(Appendix 3)
A plurality of the reference tables generated for each processing parameter are stored in the external storage device.
The signal processor selects a reference table from a plurality of the reference tables stored in the external storage device according to the processing parameters.
The synthetic opening sonar according to Appendix 2.
(Appendix 4)
Equipped with a sound velocity sensor that measures the speed of sound underwater
The processing parameters are fixed except for the underwater sound velocity.
The external storage device stores a plurality of the reference tables generated for each underwater sound velocity.
The signal processor selects a reference table from a plurality of the reference tables stored in the external storage device according to the underwater sound velocity measured by the sound velocity sensor.
The synthetic opening sonar according to Appendix 3.
(Appendix 5)
The signal processor reads and holds a reference table before and / or after the reference table in the memory together with a reference table selected according to the underwater sound velocity measured by the sound velocity sensor.
The synthetic opening sonar according to Appendix 4.
(Appendix 6)
The table coefficient is a coefficient in the form of pre-multiplying the interpolation coefficient by the replica and the reference function used for the synthetic aperture processing.
The synthetic opening sonar according to any one of Appendix 1 to 5.
(Appendix 7)
The indexes and table coefficients included in the reference table are in the form of a matrix having the same number of elements as the number of FFT points of the two-dimensional FFT.
The synthetic opening sonar according to any one of Appendix 1 to 6.
(Appendix 8)
The transmitter / receiver is provided on both the starboard and starboard sides of the underwater vehicle.
The signal processor performs the synthetic aperture processing on each of the received sound wave data received by the transmitter / receiver on both the starboard and starboard sides.
The synthetic opening sonar according to any one of Appendix 1 to 7.
(Appendix 9)
It is a signal processing method of a synthetic aperture sonar that performs synthetic aperture processing on the data of the received sound waves received by the transmitter / receiver.
Select a lookup table and
A two-dimensional FFT is executed on the data,
From the output of the two-dimensional FFT, the element corresponding to the index included in the selected reference table is extracted.
Multiply the extracted elements by the table coefficients contained in the selected reference table.
A two-dimensional Fourier is executed on the result of the multiplication,
Imaging the output of the two-dimensional Fourier,
Signal processing method for synthetic aperture sonar.
(Appendix 10)
A reference table is selected from a plurality of the reference tables stored in the external storage device, and the selected reference table is read into the memory and used.
The signal processing method for a synthetic aperture sonar according to Appendix 9.
(Appendix 11)
A plurality of the reference tables generated for each processing parameter are stored in the external storage device.
A reference table is selected from the plurality of reference tables stored in the external storage device according to the processing parameters.
The signal processing method for a synthetic aperture sonar according to Appendix 10.
(Appendix 12)
The processing parameters are fixed except for the underwater sound velocity.
The external storage device stores a plurality of the reference tables generated for each underwater sound velocity.
A reference table is selected from the plurality of reference tables stored in the external storage device according to the underwater sound velocity measured by the sound velocity sensor.
The signal processing method for a synthetic aperture sonar according to Appendix 11.
(Appendix 13)
Along with the reference table selected according to the underwater sound velocity measured by the sound velocity sensor, the reference table before and / or after the reference table is read and held in the memory.
The signal processing method for a synthetic aperture sonar according to Appendix 12.
(Appendix 14)
The table coefficient is a coefficient in the form of pre-multiplying the interpolation coefficient by the replica and the reference function used for the synthetic aperture processing.
The signal processing method for a synthetic aperture sonar according to any one of Appendix 9 to 13.
(Appendix 15)
The indexes and table coefficients included in the reference table are in the form of a matrix having the same number of elements as the number of FFT points of the two-dimensional FFT.
The signal processing method for a synthetic aperture sonar according to any one of Appendix 9 to 14.
(Appendix 16)
The transmitter / receiver is provided on both the starboard and starboard sides of the underwater vehicle.
The synthetic aperture processing is performed on each of the received sound wave data received by the transmitters and receivers on both the starboard and starboard sides.
The signal processing method for a synthetic aperture sonar according to any one of Appendix 9 to 15.
(Appendix 17)
Computer,
A signal processing program that functions as a signal processing means that performs synthetic aperture processing on the received sound wave data received by the transmitter / receiver.
The signal processing means
Select a lookup table and
A two-dimensional FFT is executed on the data,
From the output of the two-dimensional FFT, the element corresponding to the index included in the selected reference table is extracted.
Multiply the extracted elements by the table coefficients contained in the selected reference table.
A two-dimensional Fourier is executed on the result of the multiplication,
Imaging the output of the two-dimensional Fourier,
Signal processing program.
(Appendix 18)
The signal processing means selects a reference table from a plurality of the reference tables stored in the external storage device, and reads the selected reference table into the memory for use.
The signal processing program according to Appendix 17.
(Appendix 19)
A plurality of the reference tables generated for each processing parameter are stored in the external storage device.
The signal processing means selects a reference table from a plurality of the reference tables stored in the external storage device according to the processing parameters.
The signal processing program according to Appendix 18.
(Appendix 20)
The processing parameters are fixed except for the underwater sound velocity.
The external storage device stores a plurality of the reference tables generated for each underwater sound velocity.
The signal processing means selects a reference table from a plurality of the reference tables stored in the external storage device according to the underwater sound velocity measured by the sound velocity sensor.
The signal processing program according to Appendix 19.
(Appendix 21)
The signal processing means reads and holds a reference table before and / or after the reference table in the memory together with a reference table selected according to the underwater sound velocity measured by the sound velocity sensor.
The signal processing program according to Appendix 20.
(Appendix 22)
The table coefficient is a coefficient in the form of pre-multiplying the interpolation coefficient by the replica and the reference function used for the synthetic aperture processing.
The signal processing program according to any one of Appendix 17 to 21.
(Appendix 23)
The indexes and table coefficients included in the reference table are in the form of a matrix having the same number of elements as the number of FFT points of the two-dimensional FFT.
The signal processing program according to any one of Appendix 17 to 22.
(Appendix 24)
The transmitter / receiver is provided on both the starboard and starboard sides of the underwater vehicle.
The signal processing means performs the combined aperture processing on each of the received sound wave data received by the transmitters and receivers on both the starboard and starboard sides.
The signal processing program according to any one of Appendix 17 to 23.

1 合成開口ソーナー
11L,11R 送受波器
12 送受信器
13 音速センサ
14 信号処理器
15 外部記憶装置
2 合成開口ソーナー
101 送受波器
102 信号処理器
1 Synthetic aperture sonar 11L, 11R Transmitter / receiver 12 Transmitter / receiver 13 Sound speed sensor 14 Signal processor 15 External storage device 2 Synthetic aperture sonar 101 Transmitter / receiver 102 Signal processor

Claims (9)

音波の送受波を行う送受波器と、
前記送受波器で受波された受波音波のデータに対して合成開口処理を行う信号処理器と、を備え、
前記信号処理器は、
参照テーブルを選択し、
前記データに対して、2次元FFTを実行し、
前記2次元FFTの出力から、前記選択した参照テーブルに含まれるインデックスに対応する要素を抽出し、
前記抽出した要素に対して、前記選択した参照テーブルに含まれるテーブル係数を乗算し、
前記乗算の結果に対して、2次元IFFTを実行し、
前記2次元IFFTの出力を画像化
前記テーブル係数は、補間係数に対し、前記合成開口処理に使用するレプリカ及び参照関数を予め乗算した形式の係数である、
合成開口ソーナー。
A transmitter / receiver that transmits / receives sound waves,
A signal processor that performs synthetic aperture processing on the data of the received sound wave received by the transmitter / receiver is provided.
The signal processor
Select a lookup table and
A two-dimensional FFT is executed on the data,
From the output of the two-dimensional FFT, the element corresponding to the index included in the selected reference table is extracted.
Multiply the extracted elements by the table coefficients contained in the selected reference table.
A two-dimensional Fourier is executed on the result of the multiplication,
The output of the two-dimensional Fourier is imaged and
The table coefficient is a coefficient in the form of pre-multiplying the interpolation coefficient by the replica and the reference function used for the synthetic aperture processing.
Synthetic opening sonar.
前記信号処理器は、
メモリを備え、
複数の前記参照テーブルが保存された外部記憶装置に接続され、
前記外部記憶装置に保存された複数の前記参照テーブルの中から参照テーブルを選択し、前記選択した参照テーブルを前記メモリに読込んで使用する、
請求項1に記載の合成開口ソーナー。
The signal processor
Equipped with memory
A plurality of the reference tables are connected to the stored external storage device, and the reference table is connected to the storage device.
A reference table is selected from the plurality of reference tables stored in the external storage device, and the selected reference table is read into the memory and used.
The synthetic aperture sonar according to claim 1.
前記外部記憶装置には、処理パラメータ毎に生成された複数の前記参照テーブルが保存されており、
前記信号処理器は、前記処理パラメータに応じて、前記外部記憶装置に保存された複数の前記参照テーブルの中から参照テーブルを選択する、
請求項2に記載の合成開口ソーナー。
A plurality of the reference tables generated for each processing parameter are stored in the external storage device.
The signal processor selects a reference table from a plurality of the reference tables stored in the external storage device according to the processing parameters.
The synthetic aperture sonar according to claim 2.
水中音速を計測する音速センサを備え、
前記処理パラメータは、水中音速以外は固定とされ、
前記外部記憶装置には、前記水中音速毎に生成された複数の前記参照テーブルが保存されており、
前記信号処理器は、前記音速センサにて計測された水中音速に応じて、前記外部記憶装置に保存された複数の前記参照テーブルの中から参照テーブルを選択する、
請求項3に記載の合成開口ソーナー。
Equipped with a sound velocity sensor that measures the speed of sound underwater
The processing parameters are fixed except for the underwater sound velocity.
The external storage device stores a plurality of the reference tables generated for each underwater sound velocity.
The signal processor selects a reference table from a plurality of the reference tables stored in the external storage device according to the underwater sound velocity measured by the sound velocity sensor.
The synthetic aperture sonar according to claim 3.
前記信号処理器は、前記音速センサにて計測された水中音速に応じて選択した参照テーブルと共に、当該参照テーブルの前及び/又は後の参照テーブルを、前記メモリに読込んで保持する、
請求項4に記載の合成開口ソーナー。
The signal processor reads and holds a reference table before and / or after the reference table in the memory together with a reference table selected according to the underwater sound velocity measured by the sound velocity sensor.
The synthetic aperture sonar according to claim 4.
前記参照テーブルに含まれるインデックス及びテーブル係数は、前記2次元FFTのFFTポイント数と同じ要素数の行列形式である、
請求項1からのいずれか1項に記載の合成開口ソーナー。
The indexes and table coefficients included in the reference table are in the form of a matrix having the same number of elements as the number of FFT points of the two-dimensional FFT.
The synthetic opening sonar according to any one of claims 1 to 5 .
前記送受波器は、水中航走体の左右両舷にそれぞれ設けられ、
前記信号処理器は、左右両舷の前記送受波器で受波された受波音波のデータのそれぞれに対して、前記合成開口処理を行う、
請求項1からのいずれか1項に記載の合成開口ソーナー。
The transmitter / receiver is provided on both the starboard and starboard sides of the underwater vehicle.
The signal processor performs the synthetic aperture processing on each of the received sound wave data received by the transmitter / receiver on both the starboard and starboard sides.
The synthetic opening sonar according to any one of claims 1 to 6 .
送受波器で受波された受波音波のデータに対して合成開口処理を行う合成開口ソーナーの信号処理方法であって、
参照テーブルを選択し、
前記データに対して、2次元FFTを実行し、
前記2次元FFTの出力から、前記選択した参照テーブルに含まれるインデックスに対応する要素を抽出し、
前記抽出した要素に対して、前記選択した参照テーブルに含まれるテーブル係数を乗算し、
前記乗算の結果に対して、2次元IFFTを実行し、
前記2次元IFFTの出力を画像化
前記テーブル係数は、補間係数に対し、前記合成開口処理に使用するレプリカ及び参照関数を予め乗算した形式の係数である、
合成開口ソーナーの信号処理方法。
It is a signal processing method of a synthetic aperture sonar that performs synthetic aperture processing on the data of the received sound waves received by the transmitter / receiver.
Select a lookup table and
A two-dimensional FFT is executed on the data,
From the output of the two-dimensional FFT, the element corresponding to the index included in the selected reference table is extracted.
Multiply the extracted elements by the table coefficients contained in the selected reference table.
A two-dimensional Fourier is executed on the result of the multiplication,
The output of the two-dimensional Fourier is imaged and
The table coefficient is a coefficient in the form of pre-multiplying the interpolation coefficient by the replica and the reference function used for the synthetic aperture processing.
Signal processing method for synthetic aperture sonar.
コンピュータを、
送受波器で受波された受波音波のデータに対して合成開口処理を行う信号処理手段として機能させる信号処理プログラムであって、
前記信号処理手段は、
参照テーブルを選択し、
前記データに対して、2次元FFTを実行し、
前記2次元FFTの出力から、前記選択した参照テーブルに含まれるインデックスに対応する要素を抽出し、
前記抽出した要素に対して、前記選択した参照テーブルに含まれるテーブル係数を乗算し、
前記乗算の結果に対して、2次元IFFTを実行し、
前記2次元IFFTの出力を画像化
前記テーブル係数は、補間係数に対し、前記合成開口処理に使用するレプリカ及び参照関数を予め乗算した形式の係数である、
信号処理プログラム。
Computer,
A signal processing program that functions as a signal processing means that performs synthetic aperture processing on the received sound wave data received by the transmitter / receiver.
The signal processing means
Select a lookup table and
A two-dimensional FFT is executed on the data,
From the output of the two-dimensional FFT, the element corresponding to the index included in the selected reference table is extracted.
Multiply the extracted elements by the table coefficients contained in the selected reference table.
A two-dimensional Fourier is executed on the result of the multiplication,
The output of the two-dimensional Fourier is imaged and
The table coefficient is a coefficient in the form of pre-multiplying the interpolation coefficient by the replica and the reference function used for the synthetic aperture processing.
Signal processing program.
JP2018085371A 2018-04-26 2018-04-26 Synthetic aperture sonar, its signal processing method, and signal processing program Active JP7102907B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2018085371A JP7102907B2 (en) 2018-04-26 2018-04-26 Synthetic aperture sonar, its signal processing method, and signal processing program
FR1903989A FR3080686B1 (en) 2018-04-26 2019-04-15 SONAR A OPENING SYNTHESIS, ASSOCIATED SIGNAL PROCESSING PROCESS, AND SIGNAL PROCESSING PROGRAM

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018085371A JP7102907B2 (en) 2018-04-26 2018-04-26 Synthetic aperture sonar, its signal processing method, and signal processing program

Publications (2)

Publication Number Publication Date
JP2019191045A JP2019191045A (en) 2019-10-31
JP7102907B2 true JP7102907B2 (en) 2022-07-20

Family

ID=68342058

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018085371A Active JP7102907B2 (en) 2018-04-26 2018-04-26 Synthetic aperture sonar, its signal processing method, and signal processing program

Country Status (2)

Country Link
JP (1) JP7102907B2 (en)
FR (1) FR3080686B1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112099021A (en) * 2020-08-14 2020-12-18 杭州电子科技大学 Point target imaging method of synthetic aperture sonar

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007333473A (en) 2006-06-13 2007-12-27 Nec Corp Sonar system and its phase error correction method
JP2011089789A (en) 2009-10-20 2011-05-06 Mitsubishi Electric Corp Synthetic aperture radar system
JP2012108122A (en) 2010-10-28 2012-06-07 Univ Of Tokyo Submarine audio/video system
JP2012137305A (en) 2010-12-24 2012-07-19 Mitsubishi Heavy Ind Ltd Synthetic aperture exploratory device and underwater traveling body

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4953143A (en) * 1981-01-12 1990-08-28 The United States Of America As Represented By The Secretary Of The Navy Multiple frequency synthetic aperture sonar

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007333473A (en) 2006-06-13 2007-12-27 Nec Corp Sonar system and its phase error correction method
JP2011089789A (en) 2009-10-20 2011-05-06 Mitsubishi Electric Corp Synthetic aperture radar system
JP2012108122A (en) 2010-10-28 2012-06-07 Univ Of Tokyo Submarine audio/video system
JP2012137305A (en) 2010-12-24 2012-07-19 Mitsubishi Heavy Ind Ltd Synthetic aperture exploratory device and underwater traveling body

Also Published As

Publication number Publication date
FR3080686A1 (en) 2019-11-01
FR3080686B1 (en) 2021-05-21
JP2019191045A (en) 2019-10-31

Similar Documents

Publication Publication Date Title
KR102009274B1 (en) Fir filter coefficient calculation for beam forming filters
JP6533340B2 (en) Adaptive phase distortion free amplitude response equalization for beamforming applications
JP6254142B2 (en) Apparatus and method for calculating loudspeaker signals for multiple loudspeakers using delay in the frequency domain
US20120179043A1 (en) Subject diagnosis system, medical image system for providing diagnosis image of subject, and method of displaying diagnosis image of subject
JP7102907B2 (en) Synthetic aperture sonar, its signal processing method, and signal processing program
CN111768353B (en) Hole filling method and device for three-dimensional model
Choe et al. Deep point cloud reconstruction
JP6261839B1 (en) Synthetic aperture radar signal processor
CN114706076A (en) Millimeter wave near-field SAR (synthetic aperture radar) velocity imaging method based on improved range migration algorithm
JPH10503849A (en) Adaptive optimization method of ultrasonic measurement signal
CN115359771A (en) Underwater acoustic signal noise reduction method, system, equipment and storage medium
CN112669805B (en) Active noise control system of compressor in gas station based on equation error algorithm
CN116879952A (en) Calculation method, storage medium and equipment for point source elastic wave seabed reflection coefficient
JP6520276B2 (en) Noise suppression device, noise suppression method, and program
CN116359990A (en) Depth domain inversion method based on variable density acoustic wave equation
JP6997589B2 (en) Direction estimation device, direction estimation method and program
CN112729532B (en) Method, device, equipment and medium for acquiring near-field transfer function of shallow sea motion sound source
JP7375904B2 (en) Filter coefficient optimization device, latent variable optimization device, filter coefficient optimization method, latent variable optimization method, program
CN113917395A (en) Sound source positioning method and device based on reflector extension and electronic equipment
CN107329128A (en) A kind of high-precision low latitude angle-measuring method of the super-resolution based on Rapid storage technique
KR20220123655A (en) Segmentation methods, encoders, decoders and computer storage media
JP7173356B2 (en) PSD optimization device, PSD optimization method, program
JP7173355B2 (en) PSD optimization device, PSD optimization method, program
CN112666625A (en) Rapid imaging device and method for millimeter wave security inspection
CN113835067B (en) Method and device for resisting spectrum dispersion interference

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210303

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20211228

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220201

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220328

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20220607

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220620

R151 Written notification of patent or utility model registration

Ref document number: 7102907

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151