JPS62247391A - Image forming system - Google Patents
Image forming systemInfo
- Publication number
- JPS62247391A JPS62247391A JP9012986A JP9012986A JPS62247391A JP S62247391 A JPS62247391 A JP S62247391A JP 9012986 A JP9012986 A JP 9012986A JP 9012986 A JP9012986 A JP 9012986A JP S62247391 A JPS62247391 A JP S62247391A
- Authority
- JP
- Japan
- Prior art keywords
- data
- amplitude
- correction function
- phase correction
- hologram
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000012937 correction Methods 0.000 claims abstract description 96
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 238000000034 method Methods 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012552 review Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 241001123248 Arma Species 0.000 description 1
- 241000287107 Passer Species 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Landscapes
- Holo Graphy (AREA)
Abstract
Description
(産業上の利用分野)
本発明は、超音波映像システムやソーナシステム等に使
用される高解像度の映像形成方式に関するものである。
(従来の技術)
音響映像の分野における目標物の映像は、音波を目標物
に送波してその反射波を受波器アレイでとらえ受波デー
タから目標物強度を求めて得られる。この場合、目標物
と受波器あるいは受渡器アレイとの間での音波放射によ
る伝搬は空間情報の変換をなしている。横方向の分解能
をもつどのような映像装置も、像を再生するためにはこ
の変換の逆変換を必要とする。このため、従来は、遠距
離音場やフレネル領域ではフーリエ変換を用い、近距離
音場ではコンボリューション変換を用いて映像再生がな
されていた。以下従来の映像再生方法について詳細に述
べる。
まず、受波器アレイ上のホログラムデータと目標物像と
の関係を述べる。第2図に受波器アレイの存在するホロ
グラム平面と目標物の存在する目標物平面の概略図を示
す。同図に示す様に受波器かホログラム平面座標(Xi
、O)、(i=1.・・・。
N、但しNは受波器数)に存在し、それぞれの受波器位
置で受波される音波信号(以後ホログラムデータと呼ぶ
)をS□(X 、、O) 、(i = 1 、 ・・・
。
N)とする。一方、目標物は同図に示す様に目標物平面
座標(x n、 z n)に存在し、目標物の反射係数
(以後目標物像と呼ぶ)をρ(x (1,z n)と表
わす。このとき、ホログラムデータと目標物像との関係
は、近距離音場、フレネル領域、遠距離音場のすべての
領域の音場を対象として、音場の波動理論から、
で示される。ここに2.、は目標物平面−ホログラム平
面量比w1(以後、単に距離と呼ぶ)である。
またKは波数であり、目標物に送波する音波の送2πf
信用波数f、媒質中の音速■を用いてに=−■
と示される。なお、jは虚数単位である。又、。−・K
〜η;;:フTア1−−フー/、へx、−)で−ゴー−
ニジー1−1− は目標物位置(x no z n)か
ら球面拡散波として放射される音波13号の伝搬関数で
ある。
ここで、遠距離音場やフレネル領域でフーリエ変換を用
いて映像再生する方法について詳細に述べる。この方法
は、たとえば、ウルトラソニックス シンポジウム プ
ロシーディング(Ultra−sonics Symp
、 Proc、)、「超音波映像再生のためのA RM
A処理(ARMA Processing for U
ltrasonicleconstructive I
maging)」ヨー ハン バオ(Yoh、 Han
Pao) 、1982年、(A%9.27−29、I
EEE。
ニューヨーク(1982年))に記載されている。
l Xi−X+l<zn (i=1.−、N)が成立す
る場合、すなわち、目標物平面上における目標物の存在
する範囲や、アレイ開口が距gznに比べて充分小さい
という場合には、フレネル近似が成立する。そこで、フ
レネル近似を用いて前記伝搬関数の分母(x 、−X
、)”+z 、””z nと置き、伝搬関数の分子の位
相項−jに (x、−X工” +z、”−−jにz 、
(1+ (x n−X 、)2/ (2z 、”))と
置く。そして、目標物像とホログラムの関係式(1)を
次の様に近似する。
・−(2)
次に、ホログラムS、(Xi、O)に振幅と位相の補正
を行い、補正されたホログラム、
5H(X+ 、0) :SH(L 、0) Z、 e
JK”11”◆Xi /+2z、” ))、、、 (3
)を算出する。そこで、弐(2)と式(3)から補正さ
れたホログラムSH(Xi、O)と目標物像ρ(xn。
zn)の関係が、
のように求まる。ここで補正されたホログラムS、4(
Xi、O)と、目標物像ρ(xn、zn)とe−JKX
i” /+2z・1の積は、フーリエ変換の関係にある
。従って補正されたホログラムSH(Xi、O)をフー
リエ逆変換し、その値をe−JKXi’ / +2z・
)で割ると、目標物像゛ρ(xn、zn)が次の様に得
られここに、Lは開口長さであり、目標物の高解像度化
をはかるためには大きくとらなければならない値である
。従来はフレネル近似の成立する範囲で開口りを拡張し
解像度を上げていた。
次に、近距離音場でコンボリューション変換を用いて映
像再生する方法について詳細に述べる。
この方法は、たとえば、「波長解像度制限のないホログ
ラフ映像(tlolographic Imaging
wiLhoutthe Wavelength Re
5olution Lim1t)」 アールジー ウィ
リアムズ(Earl G、 Williams)、フイ
ジカ)L/ レビュー レターズ(Physical
ReviewLetters)、旦、1980年に記
載されている。
前記式(1)はコンボリューションの形であり、フーリ
エ変換すると、
ン(S、4(ξ))−Eマひ(ρ(ζ))び(P(ξ)
)−(6)と示される。ここに、7(S、(E))はホ
ログラムs、O)のフーリエ変換9.y(ρ(ξ))は
目標物像ρ(x n、z n)のフーリエ変換、V7(
P(ξ))は。−J KF糖マ/ t 24 z 、
2のフーリエ変換であり、あり、
と示される。ここにHQ”は第2種バンケル関数であり
、に。は第2種変形ベッセル関数である。
目標物像ρ(x n、z n)は、式(7) 、 (8
)を式(6)に代入し逆フーリエ変換することによって
求められ、
ρ(xn、zo) = 3−’ (ρ(E ))と得
られる。ここにJ−1はフーリエ逆変換を表わす。
通常は、式(9)の一番目の式)−’[,7(S□(ξ
))/ (−j 、Ho(2J znE”τ]))]の
みを用いることにより像の再生が行なわれる。この場合
の再生像は高解像度な像ではない。そこで更に解像度を
上げるために5.?−’[,7(S□(ζ))/(−j
πHO(2)(z n5下]コ))]の情情報性にン一
’[JG(SH(ξ))/(2にo<z nπ1コQ>
>]の情報も用いられる。ここでに。はznfの増大に
ともなって急激に小さくなる値である。
そこで従来はこうしたに。の性質により、3−” [
J (S)+(ξ))/ (2K O(Z 、n>
)コが発散するのを押さえるため、距離2゜を極めて小
さく、1波長程度以内(超近距離音場)に設定し、1ξ
1〉になる範囲でξを大きくとれる様にして、フーリエ
逆変換を行ない解像度の高い像を得ていた。
(発明が解決しようとする問題点)
しかしながら、上記方法のうち、近距離音場においてフ
ーリエ変換を用いる場合は、フレネル近似が成立しない
。そのため訂記伝搬関数のフレネル近似による誤差が増
大し、実測ホログラムデータの訂記振幅と位相の補正に
も誤差が生じ再生像が劣化するという問題点があった。
一方、コンボリューション変換を用いる場合は、解像度
を上げようとすると、眞記式(9)の二番目の情報、l
i’−’[,7(SH(ξ))/(2にo (z 、E
「1)]も用いなければならない。しかしこの情報が距
離2、の増大にともなって急激に増大するため、近距離
場においてもznの増大に対して、ホログラム信号中の
ノイズを急激に増幅させることになる。従って、再生像
の劣化が生じるという問題点、があった。
本発明は、以上述べた近距離音場において目標l物の存
在する範囲のみが距離に比べて充分小さいという条件で
高解像度な目標物像を得る場合、大きなアレイ開口でフ
レネル近似を使ったときの近似上の誤差により再生像が
劣化する問題点や、コンボリューション変換を用いた場
合のノイズ増幅により再生像が劣化する問題点を解決し
、安定で高解像度な再生像が得られる映像形成方式を提
供することを目的とする。
(問題点を解決するための手段)
本発明は、目標物に音波を送波してその反射波をN個の
受波器から成る受渡器アレイで受波し、受波した実測ホ
ログラムデータに基づいて目標物像の再生を行う映像形
成方式に係るもので、面記従来技術の問題点を解消する
ため、ホログラム平面と目標物平面間の距離データzn
と受渡器の位置データXiから振幅補正開数 Xi”+
z、”を算出する振幅補正関数算出手段と:受波したホ
ログラムデータ5H(Xi、0)を、振幅補正関数算出
手段の出力により振幅補正し、振幅補正されたホログラ
ムデータS□O) 、ら(,24−乙]を得る振幅補
正手段と;目標物に送波する音波の周波数データfと音
速データVから波数K(=2πf/V)を算出する波数
算出手段と:振幅補正関数算出手段の出力と波数算出手
段の出力から位相補正関数e J K tを算出する位
相補正関数算出手段と:振幅補正手段の出力する振幅補
正されたホログラムデータを、位相補正関数算出手段の
出力により位相補正し、振幅補正及び位相補正されたホ
ログラムデータS1.I(Xi、O)mゴe・に5耳T
を得る位相補正手段と:受波器の位置データと振幅補正
関数算出手段の出力から、新たな変数でl (=xr
/ X % + Z n”)を求める変数変換手段と
二波数算出手段、位相補正手段及び変数変換手段の各出
力並びに目標物平面座標上のサンプリング位置データX
。を受取り、位相補正手段の出力にフーリエ逆変換を施
し、目標物像データ(ただしSH(τ1.0)は位相補
正手段の出力と同値で該出力をτ1 で表わしたデータ
)を得るフーリエ逆変換手段とを設けたものである。
る。
(作 用)
本発明では、以上のように映像形成方式を構成したので
各技術手段は以下の通り作用する。
振幅補正関数算出手段は、距離に比べて大きさの制限の
ない受波器アレイ開口全体で受波して得られるホログラ
ムデータに対する振幅補正開数Xi’+z、’を算出し
、振幅補正手段、位相補正関数算出手段及び変数変換手
段に出力する。これは、従来フレネル近似を使って伝搬
関数分母のπ7τX万”+z、’を〜2.と近似したか
わりに、距離znに比べて目標物の存在する範囲のみが
充分小さいという条件から成立するl 2 XrI X
if(<X;’ + z 、”の関係にもとづいて、
Rマ巧こゴ75(Industrial Application Field) The present invention relates to a high-resolution image forming method used in ultrasonic imaging systems, sonar systems, and the like. (Prior Art) An image of a target in the field of audiovisuals is obtained by transmitting sound waves to the target, capturing the reflected waves with a receiver array, and determining the intensity of the target from the received wave data. In this case, the propagation of the acoustic radiation between the target and the receiver or passer array constitutes a transformation of spatial information. Any imaging device with lateral resolution requires the inverse of this transformation to reproduce the image. For this reason, conventionally, video has been played using Fourier transform in the far sound field or Fresnel region, and convolution transform in the near sound field. The conventional video playback method will be described in detail below. First, the relationship between the hologram data on the receiver array and the target object image will be described. FIG. 2 shows a schematic diagram of a hologram plane where a receiver array exists and a target plane where a target exists. As shown in the figure, the receiver or hologram plane coordinates (Xi
, O), (i=1...N, where N is the number of receivers), and the sound wave signals (hereinafter referred to as hologram data) received at each receiver position are S□ (X,,O),(i=1,...
. N). On the other hand, as shown in the figure, the target exists at the target plane coordinates (x n, z n), and the reflection coefficient of the target (hereinafter referred to as the target image) is expressed as ρ(x (1, z n)). At this time, the relationship between the hologram data and the target object image is expressed by the following from the wave theory of the sound field, targeting the sound field in all regions of the near sound field, Fresnel region, and far sound field. Here, 2. is the target object plane - hologram plane amount ratio w1 (hereinafter simply referred to as distance). Also, K is the wave number, 2πf of the sound wave transmitted to the target object, the reliable wave number f, and the wave number in the medium. Using the sound speed ■, it is shown as =-■. Note that j is an imaginary unit. Also, -・K
~η;;: FuTa1--fu/, hex,-)de-go-
1-1- is the propagation function of the sound wave No. 13 emitted from the target position (x no z n) as a spherical diffused wave. Here, a method for reproducing images using Fourier transform in a far sound field or Fresnel domain will be described in detail. This method is described, for example, in the Ultrasonics Symposium Proceedings.
, Proc, ), “A RM for Ultrasound Image Reproduction
A processing (ARMA Processing for U
ltrasonicleconstructive I
Yoh, Han Bao
Pao), 1982, (A%9.27-29, I
EEE. New York (1982)). If l Xi−X+l<zn (i=1.−, N) holds, that is, if the range where the target exists on the target plane or the array aperture is sufficiently small compared to the distance gzn, Fresnel approximation holds. Therefore, using the Fresnel approximation, the denominator of the propagation function (x, −X
,)"+z,""z n, and the phase term of the numerator of the propagation function -j is (x, -X"+z,"--j is z,
(1+ (x n - Amplitude and phase correction is performed on S, (Xi, O), and the corrected hologram, 5H (X+, 0): SH (L, 0) Z, e
JK”11”◆Xi /+2z,” )),,, (3
) is calculated. Therefore, from (2) and equation (3), the relationship between the corrected hologram SH (Xi, O) and the target object image ρ (xn. zn) is determined as follows. Here the corrected hologram S,4(
Xi, O), target object image ρ(xn, zn) and e-JKX
The product of i''/+2z・1 is in a Fourier transform relationship. Therefore, the corrected hologram SH(Xi, O) is inversely Fourier transformed and its value is expressed as e-JKXi'/+2z・
), the target image ρ(xn, zn) is obtained as follows, where L is the aperture length, which must be large in order to achieve high resolution of the target. It is. Conventionally, the resolution was increased by expanding the aperture within the range where Fresnel approximation holds. Next, a method for reproducing video using convolution transformation in a near-field sound field will be described in detail. This method can be used, for example, for holographic imaging without wavelength resolution limitations.
wiLoutthe Wavelength Re
5olution Lim1t)" Earl G Williams (Earl G, Williams), Fujika) L/Review Letters (Physical
Review Letters), Dan, 1980. The above equation (1) is in the form of convolution, and when Fourier transform is performed, it becomes
)−(6). Here, 7(S, (E)) is the Fourier transform of the hologram s, O)9. y(ρ(ξ)) is the Fourier transform of the target object image ρ(x n, z n), V7(
P(ξ)) is. -J KF glycoma/t24z,
is the Fourier transform of 2, and is denoted by . Here, HQ'' is the Wankel function of the second kind, and . is the modified Bessel function of the second kind.
) is substituted into equation (6) and subjected to inverse Fourier transform, and ρ(xn, zo) = 3−' (ρ(E)) is obtained. Here, J-1 represents the inverse Fourier transform. Usually, the first equation in equation (9))−'[,7(S□(ξ
))/ (-j , Ho(2J znE”τ]))]. The reconstructed image in this case is not a high-resolution image. Therefore, in order to further increase the resolution, .?-'[,7(S□(ζ))/(-j
πHO(2)(z n5 lower] ko))
>] information is also used. Here. is a value that rapidly decreases as znf increases. So this is how it used to be. Due to the property of 3-” [
J (S) + (ξ))/ (2K O(Z, n>
) In order to suppress the divergence of ko, the distance 2° is set extremely small, within about 1 wavelength (ultra short distance sound field), and the distance 2°
A high-resolution image was obtained by performing inverse Fourier transform while allowing ξ to be large within the range of 1>. (Problems to be Solved by the Invention) However, among the above methods, when Fourier transform is used in the near field, Fresnel approximation does not hold. Therefore, errors due to Fresnel approximation of the correction propagation function increase, errors also occur in the correction of correction amplitude and phase of the measured hologram data, and the reproduced image deteriorates. On the other hand, when using convolution transformation, if you try to increase the resolution, the second information in Makki equation (9), l
i'-'[,7(SH(ξ))/(2 to o (z, E
"1)] must also be used. However, since this information increases rapidly as the distance 2 increases, the noise in the hologram signal will be rapidly amplified even in the near field as zn increases. Therefore, there was a problem in that the reproduced image deteriorated.The present invention is based on the condition that only the range in which the target object exists in the near-field sound field described above is sufficiently small compared to the distance. When obtaining a high-resolution target image, there is a problem that the reproduced image deteriorates due to approximation errors when using Fresnel approximation with a large array aperture, and the reproduced image deteriorates due to noise amplification when using convolution transformation. It is an object of the present invention to provide an image forming method that can obtain stable and high-resolution reproduced images. This relates to an image forming method in which the reflected waves are received by a receiver array consisting of N receivers, and a target object image is reproduced based on the received measured hologram data. In order to solve the problem, distance data zn between the hologram plane and the target plane
From the position data Xi of the transfer device and the amplitude correction numerical value Xi”+
amplitude correction function calculation means for calculating the amplitude correction function calculation means: the received hologram data 5H (Xi, 0) is amplitude-corrected by the output of the amplitude correction function calculation means, and the amplitude-corrected hologram data S□O), et al. (,24-B); wave number calculation means for calculating the wave number K (=2πf/V) from the frequency data f of the sound wave transmitted to the target object and the sound velocity data V; and the amplitude correction function calculation means and a phase correction function calculation means for calculating a phase correction function e J K t from the output of the amplitude correction means and the output of the wave number calculation means. Then, the amplitude-corrected and phase-corrected hologram data S1.
phase correction means for obtaining: l (=xr
/
. and performs inverse Fourier transform on the output of the phase correction means to obtain target object image data (where SH (τ1.0) is the same value as the output of the phase correction means and the output is expressed as τ1). The means are provided. Ru. (Function) In the present invention, since the image forming system is configured as described above, each technical means functions as follows. The amplitude correction function calculation means calculates an amplitude correction numerical aperture Xi'+z,' for the hologram data obtained by receiving the wave with the entire receiver array aperture, which has no size limit compared to the distance, and Output to phase correction function calculation means and variable conversion means. This holds true from the condition that only the range where the target exists is sufficiently small compared to the distance zn, instead of approximating the propagation function denominator π7τ 2 XrI
Based on the relationship “if(<X;’ + z,”
R makogo 75
【3−% 、 24 z 、2と近似す
るものである。
振幅補正手段は、受波したホログラムデータS□(Xi
、O)に対し、振幅補正関数算出手段で求めた振幅補正
関数〜X 、24 z 、”を乗じることにより、ホロ
グラムデータの振幅の補正を行ない、振幅補正されたホ
ログラムデータS 、(X 、、O)nを位相補正手段
に出力する。従来
のように (x 、−x 1)” + z 、”yz
nと近似した場合は、アレイ開口が大きいとき、すなわ
ちX。
が大のとき、誤差が大きくなる。しかし本発明では、
5τフ7〒z 、2〜 Xi”+z、”と近似すること
により誤差が小さくなる。したがって、受波器アレイで
受波したホログラムデータを振幅補正手段に通すことに
より、球面拡散波として放射される音波43号の振幅に
精度の良い補正が行なわれる。
波数算出手段は、位相補正関数の算出とフーリエ逆変換
に使用する波数を算出し、波数データを位相補正関数算
出手段とフーリエ逆変換手段とに供給する。
位相補正関数算出手段は、振幅補正手段の出力する振幅
補正されたホログラムデータに位相補正をするため位相
補正関数B J K J口1マを算出し、位相補正手段
に出力する。これは、従来フレネル近似を使って、伝搬
関数分子の位相項一 J K (X n−X +)2
を” −j K zn(1+ (Xn−Xl)
2/(2Z、”))と近似したかわりに−jK(t
x n x i/ (x 、2+z n’))と近似す
るものである。
位相補正手段は、振幅補正されたホログラムデータSH
(Xi、O) 、1”>こ7;7コに対し、位相補正
関数算出手段で求めた位相補正関数eJ K 2を乗じ
ることによりホログラムデータの位相の補正を行ない、
振幅補正及び位相補正されたホログラムデータ5o(X
t、0) Xi”+z、’6バ5ア7σをフーリエ逆
変換手段に出力する。振幅補正されたホログラムデータ
を位相補正手段に通すことにより、球面拡散波として放
射される音波信号の位相に粒度の良い位相補正が行なわ
れる。すなわち、従来の様に−jK#=x、−X)24
z 、2QjL−j Kzn(1+ (xn−Xi)
2)/ (2zn”))と近似した場合は、アレイの開
口が大きいとき、すなわち、Xiが大のとき、この誤差
が大きくなるが、本発明では−jKzへ−(、−Xコ2
+ Z 、2−L−jに X%+ z、2 (I
Xixi / (Xi”+Z、2))と近似することに
より、誤差は小となる。
したがって、振幅の補正されたホログラムデータ5H(
Xt、O)5τT7に位相の補正をしたテ’) S H
(X r 、 0 ) fi e ”’ア”下は積度
の良い補正データとなるのである。
変数変換手段はフーリエ逆変換における演算のためにτ
、=Xi/ Xi’+z、’なる式にしたがい変数変
換を行ない、そのデータをフーリエ逆変換手段に出力す
る。
フーリエ逆変換手段は、振幅及び位相の補正されたホロ
グラムデータにフーリエ逆変換を施し、目標物像データ
ρ(x yl、z n)を得る。位相補正手段より出力
される振幅及び位相の補正されたホログラムデータは目
標物像データとの関係においてフーリエ変換
SH(で1.0)は、補正されたホログラムデータ5n
(Xl、O) 5丁フeaKfiア7フと同値で、τ、
=Xi/シロ5T「;覆−】の変数変換により、T1の
関数として表わした補正ホログラムデータである。そこ
で、振幅・位相の補正されたホログラムデータをフーリ
エ逆変換手段に通すことにより、目標物像ρ(xn、z
n)が再生される。
したがって、本発明によれば、アレイ開口がフレネル近
似の成立しない様な大きな場合でも全開口ホログラムデ
ータに積度の良い振幅及び位相の補正がなされ、高解像
度な目標物像が得られるようになり、前記従来技術の問
題点が解決される。
(実施例)
以下、本発明の実施例を詳細に説明する。
第1図は本実施例の映像形成システムの構成を示すブロ
ック図である。この映像形成システムは、目標物に音波
を送波してその反射波を受渡器アレイで受波して得た実
測ホログラムデータから目標物像の再生を行うもので、
振幅補正関数算出手段1、振幅補正手段2、波数算出手
段3、位相補正関数算出手段4、位相補正手段5、変数
変換手段6及びフーリエ逆変換手段7より構成される。
ここでは、−例として受波器アレイの受渡器の個数をN
とする。また第1図中、T1はホログラム平面と目標物
平面間の距離データzn(第2図参照)の入力端子、T
2は受波器の位置データXi(i=1. ・・・、N)
(7)入力端子、T3は受渡器で受波したホログラムデ
ータSH(Xi、O) 、(i=1.−、N)の入力端
子、T4は目標物に送波する音波の送信周波数データf
の入力端子、T5は音速データVの入力端子、T6は目
標物平面座標上のサンプリングポイントデータxnの入
力端子、T7は目標物像データρ(Xn、zn)の出力
端子である。各要素の接続関係について述べると、入力
端子T、は振幅補正関数算出手段1に接続される。入力
端子T2は振幅補正関数算出手段1及び変数変換手段6
に接続される。入力端子T3は振幅補正手段2に接続さ
れる。入力端子T4.T5は波数算出手段3に接続され
る。入力端子T6はフーリエ逆変換手段7に接続される
。振幅補正関数算出手段1の出力は振幅補正手段2、位
相補正関数算出手段4及び変数変換手段6の各人力に接
続される。振幅補正手段2の出力は位相補正手段5の人
力に接続され、位相補正手段5の出力はフーリエ逆変換
手段7の入力に接続される。波数算出手段3の出力は位
相補正関数算出手段4及びフーリエ逆変換手段7の各入
力に接続される。位相補正関数算出手段4の出力は位相
補正手段5の人力に接続される。変数変換手段6の出力
はフーリエ逆変換手段7の人力に接続される。更にフー
リエ逆変換手段7は目標物像データを出力する目標物像
出力端子T、に接続されている。
次に動作について説明する。先ず入力端子T1から距離
データznが振幅補正関数算出手段1に送られ記憶され
る。入力端子T6から目標物平面座標サンプリングポイ
ントデータX。が入力され、フーリエ逆変換手段7に送
られ記憶される。
入力端子T2から受波器位置データX1(i=1゜−、
N)が人力され、振幅補正関数算出手段1と変数変換手
段6に送られそれぞれ記憶される。入力端子T4から送
信周波数データfが人力され波数算出手段3に送られ記
憶される。入力端子T。
から音速データVが入力され波数算出手段3に送られ記
憶される。振幅補正関数算出下段1では、入力端子TI
から送られてきた距離データz、と入力端子T2から送
られてきた受波器位置データXi(i=1.・−、N)
とを用いて振幅補正関数データ 、5こフロE覆−フ、
(i=1.−、N)を求める(+3i算が行なわれる。
そして求められ振幅補正関数データパーT;7)、
(i=1.・−、N)は、振幅補正手段2と、位相補正
関数算出手段4と変数変換手段6に送られ、それぞれ記
憶される。波数算出手段3では、入力端子T4から入力
された送信周波数データfと、入力端子T5から人力さ
れた音速データVを用いて、波数データに=2πf/V
を求める演算が行なわれる。算出された波数データには
、位相補正関数算出手段4とフーリエ逆変換手段7に送
られそれぞれ記憶される。位相補正関数算出手段4では
、振幅補正関数算出下段1から送られてきた振幅補正関
数データ5亜7フ、(i=1.−、N)と波数算出手段
3から送られてきた波数データKを用いて、位相補正関
数チーJ e ””T”’ 、 (i = 1、−、
N )を求める演算が行なわれる。そして、このデータ
は、位相補正手段5に送られ記憶される。変数変換手段
6では、入力端子T2から送られてきた受波器位置デー
タXI、(i =1.−ms、 N)と振幅補正関数算
出手段1から送られてきた振幅補正関数データ、bε、
”+=フ−7,(i=1.−、N)を用いて、受波器位
置の変数変換されたデータτ1=x+/(>ε、24
z、”、(i = 1、−、 N )を求める演算が行
なわれる。そして、このデータはフーリエ逆変換手段7
に送られ記憶される。
次に、入力端子T3から実測ホログラムデータ5H(X
t、0 ) 、(i = 1 、 m、 N )が入力
され、振幅補正手段2に送られる。ここでは、振幅補正
関数算出手段1から送られてきた振幅補正関数データに
T+ Z n”、 (1=1 、 ”・、 N )を実
測ホログラムデータに乗゛じる演算S□(Xi、O)〜
x 、24 z 、”。
(i=1. ・・・、N)が行なわれ、実測ホログラム
データの振幅補正がなされる。振幅補正されたホログラ
ムデータは、位相補正手段5に送られる。
ここでは、位相補正関数算出手段4から送られてきた位
相補正関数データeバム77茫、(1=1゜−、N)を
振幅の補正されたホログラムデータに乗じる演算5s(
Xt、0)’4こ、 i 4 z 、26 J K −
〃y7寥。
(i= 1 、−−−、 N )が行なわれ、振幅補正
されたホログラムデータに位相の補正がなされる。この
様にして、振幅及び位相の補正されたホログラムデータ
は、フーリエ逆変換手段7に送られる。ここでは、振幅
と位相の補正されたホログラムデータ5HO)5戸+Z
n’eバ1が、変数変換手段6より送られてきた受波器
位置データ0)に置き換えられる。更に波数算出手段3
から送られてきた波数にと入力端子T6から送られてき
た目標物平面座標サンプリングポイントデータx、とS
H(τ1.O)を用いて、目標物像データρ(x ll
+ z n)をフーリエ逆変換、
ρ(xn、zn) =ΣS、(τi 、o) e −J
K”aτ・によって算出する。算出された目標物像デー
タρ(x n、 z n)は出力端子T?から再生像と
して出力される。
(発明の効果)
以上詳細に説明したように、本発明によれば次のような
利点がある。
■目標物像の高解像度化をはかるため、フレネル近似が
成立しないような大きさにアレイ開口を拡張した場合で
も開口全体の実測ホログラムデータに振幅と位相の精度
の良い補正がなされるので、安定で高解像度な目標物像
が得られる。
■フレネル近似の成立しない近距離音場でも目標物像の
再生が可能なので、音波信号の伝搬時間が短縮され、処
理効率が上がる。[3-%, 24 z, approximates 2. The amplitude correction means receives the received hologram data S□(Xi
, O) is multiplied by the amplitude correction function ~ O) Output n to the phase correction means. As in the conventional case, (x, -x 1)" + z, "yz
When approximated to n, when the array aperture is large, that is, X. When is large, the error becomes large. However, in the present invention,
The error is reduced by approximating 5τf7〒z,2~Xi''+z,''. Therefore, by passing the hologram data received by the receiver array through the amplitude correction means, the amplitude of the sound wave No. 43 emitted as a spherical diffused wave can be accurately corrected. The wave number calculation means calculates a wave number used for calculation of the phase correction function and inverse Fourier transform, and supplies wave number data to the phase correction function calculation means and the inverse Fourier transform means. The phase correction function calculation means calculates a phase correction function BJKJ and outputs it to the phase correction means in order to perform phase correction on the amplitude-corrected hologram data output from the amplitude correction means. This can be calculated by using the conventional Fresnel approximation to calculate the phase term of the numerator of the propagator 1 J K (X n - X +) 2
” −j K zn(1+ (Xn-Xl)
2/(2Z,”)) instead of approximating −jK(t
x n x i/(x, 2+z n')). The phase correction means uses amplitude-corrected hologram data SH.
(Xi, O), 1''>7; The phase of the hologram data is corrected by multiplying 7 by the phase correction function eJ K 2 obtained by the phase correction function calculation means,
Hologram data 5o (X
t, 0) Xi"+z, '6ba5a7σ is output to the Fourier inverse transform means. By passing the amplitude-corrected hologram data through the phase correction means, the phase of the sound wave signal emitted as a spherical diffuse wave is adjusted. Fine-grained phase correction is performed.In other words, as in the conventional case, −jK#=x, −X)24
z , 2QjL-j Kzn(1+ (xn-Xi)
2)/ (2zn”)), this error increases when the array aperture is large, that is, when Xi is large, but in the present invention, −jKz to −(, −Xko2
+ Z, 2-L-j to X% + z, 2 (I
By approximating Xixi / (Xi"+Z, 2)), the error becomes small. Therefore, the amplitude-corrected hologram data 5H (
Xt, O)5τT7 with phase correction Te') S H
(X r , 0) fi e The lower part of "'a" becomes correction data with good integration. The variable transformation means is τ for the operation in the Fourier inverse transform.
,=Xi/Xi'+z,' Variable transformation is performed according to the formula, and the data is output to the inverse Fourier transform means. The inverse Fourier transform means performs inverse Fourier transform on the hologram data whose amplitude and phase have been corrected to obtain target object image data ρ(x yl, z n). The amplitude and phase corrected hologram data output from the phase correction means has a Fourier transform SH (at 1.0) in relation to the target object image data, which is the corrected hologram data 5n.
(Xl, O) Same value as 5th feaKfi a7f, τ,
This is corrected hologram data expressed as a function of T1 by the variable transformation of =Xi/shiro 5T ";reverse-".Then, by passing the hologram data whose amplitude and phase have been corrected through an inverse Fourier transform means, the target object is Image ρ(xn, z
n) is played. Therefore, according to the present invention, even when the array aperture is so large that the Fresnel approximation does not hold, the full-aperture hologram data is subjected to well-integrated amplitude and phase correction, and a high-resolution target image can be obtained. , the problems of the prior art are solved. (Example) Examples of the present invention will be described in detail below. FIG. 1 is a block diagram showing the configuration of the image forming system of this embodiment. This image forming system reproduces the image of the target object from measured hologram data obtained by transmitting sound waves to the target object and receiving the reflected waves with a receiver array.
It is composed of amplitude correction function calculation means 1, amplitude correction means 2, wave number calculation means 3, phase correction function calculation means 4, phase correction means 5, variable conversion means 6, and Fourier inverse transformation means 7. Here, as an example, the number of transfer devices in the receiver array is N
shall be. In addition, in FIG. 1, T1 is an input terminal for distance data zn (see FIG. 2) between the hologram plane and the target plane;
2 is the receiver position data Xi (i=1....,N)
(7) Input terminal, T3 is the input terminal for the hologram data SH (Xi, O), (i=1.-, N) received by the delivery device, T4 is the transmission frequency data f of the sound wave sent to the target object
T5 is an input terminal for sound velocity data V, T6 is an input terminal for sampling point data xn on the target object plane coordinates, and T7 is an output terminal for target object image data ρ(Xn, zn). Regarding the connection relationship of each element, the input terminal T is connected to the amplitude correction function calculation means 1. Input terminal T2 is connected to amplitude correction function calculation means 1 and variable conversion means 6.
connected to. Input terminal T3 is connected to amplitude correction means 2. Input terminal T4. T5 is connected to the wave number calculation means 3. Input terminal T6 is connected to inverse Fourier transform means 7. The output of the amplitude correction function calculation means 1 is connected to each of the amplitude correction means 2, the phase correction function calculation means 4, and the variable conversion means 6. The output of the amplitude correction means 2 is connected to the input of the phase correction means 5, and the output of the phase correction means 5 is connected to the input of the inverse Fourier transform means 7. The output of the wave number calculation means 3 is connected to each input of the phase correction function calculation means 4 and the inverse Fourier transform means 7. The output of the phase correction function calculation means 4 is connected to the human power of the phase correction means 5. The output of the variable transformation means 6 is connected to the input of the inverse Fourier transformation means 7. Furthermore, the inverse Fourier transform means 7 is connected to a target image output terminal T, which outputs target image data. Next, the operation will be explained. First, distance data zn is sent from the input terminal T1 to the amplitude correction function calculation means 1 and stored therein. Target plane coordinate sampling point data X is input from input terminal T6. is input and sent to the inverse Fourier transform means 7 and stored. Receiver position data X1 (i=1°-,
N) is input manually and sent to the amplitude correction function calculation means 1 and the variable conversion means 6 and stored therein. The transmission frequency data f is manually inputted from the input terminal T4 and sent to the wave number calculation means 3, where it is stored. Input terminal T. The sound velocity data V is inputted from the wave number calculation means 3 and is stored therein. In the lower stage 1 of amplitude correction function calculation, the input terminal TI
distance data z sent from the input terminal T2, and receiver position data Xi (i=1..-, N) sent from the input terminal T2.
Using the amplitude correction function data,
(i = 1.-, N) is calculated (+3i calculation is performed. Then, the obtained amplitude correction function data par T; 7),
(i=1..-, N) is sent to the amplitude correction means 2, the phase correction function calculation means 4, and the variable conversion means 6, and is stored respectively. The wave number calculation means 3 uses the transmission frequency data f input from the input terminal T4 and the sound speed data V manually input from the input terminal T5, and converts the wave number data to =2πf/V.
A computation is performed to find . The calculated wave number data is sent to the phase correction function calculation means 4 and the inverse Fourier transform means 7 and stored therein. The phase correction function calculation means 4 uses the amplitude correction function data 5, (i=1.-, N) sent from the lower amplitude correction function calculation stage 1 and the wave number data K sent from the wave number calculation means 3. Using the phase correction function Qi J ``”T”', (i = 1, -
An operation is performed to obtain N ). This data is then sent to the phase correction means 5 and stored. The variable conversion means 6 converts the receiver position data XI, (i = 1.-ms, N) sent from the input terminal T2 and the amplitude correction function data, bε, sent from the amplitude correction function calculation means 1, to
”+=Fu-7, (i=1.-, N), the variable-converted data of the receiver position τ1=x+/(>ε, 24
An operation is performed to obtain z,'', (i = 1, -, N). Then, this data is transferred to the Fourier inverse transform means 7.
is sent to and stored. Next, the measured hologram data 5H (X
t, 0) and (i = 1, m, N) are input and sent to the amplitude correction means 2. Here, an operation S□(Xi,O ) ~
x, 24 z,". (i=1....,N) is performed, and the amplitude of the measured hologram data is corrected. The amplitude-corrected hologram data is sent to the phase correction means 5. Here, , an operation 5s (5s) of multiplying the amplitude-corrected hologram data by the phase correction function data eBum77, (1=1°-,N) sent from the phase correction function calculation means 4.
Xt, 0) '4 ko, i 4 z, 26 J K -
〃y7寥. (i=1, ---, N) is performed, and the amplitude-corrected hologram data is subjected to phase correction. The hologram data whose amplitude and phase have been corrected in this manner is sent to the inverse Fourier transform means 7. Here, the amplitude and phase corrected hologram data 5HO) 5+Z
n'e bar 1 is replaced with the receiver position data 0) sent from the variable conversion means 6. Furthermore, wave number calculation means 3
, the target plane coordinate sampling point data x sent from the input terminal T6, and S
Using H(τ1.O), target object image data ρ(x ll
+ z n) by inverse Fourier transform, ρ(xn, zn) = ΣS, (τi, o) e −J
The calculated target object image data ρ(x n, z n) is output as a reproduced image from the output terminal T? (Effects of the Invention) As explained in detail above, the present invention According to the invention, there are the following advantages: ■In order to increase the resolution of the target object image, even if the array aperture is expanded to a size that does not satisfy the Fresnel approximation, the amplitude and amplitude of the measured hologram data of the entire aperture are Since the phase is corrected with high precision, a stable and high-resolution target object image can be obtained. ■The target object image can be reproduced even in a near-field sound field where Fresnel approximation does not hold, so the propagation time of the sound wave signal is shortened. processing efficiency.
第1図は本発明の一実施例の映像形成システムの構成を
示すブロック図、第2図は目標物平面とホログラム平面
との関係を示す図である。
1・・・振幅補正関数算出手段 2・・・振幅補正手段
3・・・波数算出手段 4・・・位相補正関数算出手段
5・−・位相補正手段 6・−変数変換手段7・・−フ
ーリエ逆変換手段
ffメffJPII 9(b′t $ UヴζムQlk
J ta%l’jlf、54fl!1第2図FIG. 1 is a block diagram showing the configuration of an image forming system according to an embodiment of the present invention, and FIG. 2 is a diagram showing the relationship between a target object plane and a hologram plane. 1...Amplitude correction function calculation means 2...Amplitude correction means 3...Wave number calculation means 4...Phase correction function calculation means 5...Phase correction means 6.-Variable conversion means 7...-Fourier Inverse conversion meansffmeffJPII 9(b't $UvζmQlk
J ta%l'jlf, 54fl! 1Figure 2
Claims (1)
成る受波器アレイで受波し、受波した実測ホログラムデ
ータに基づいて目標物像の再生を行う映像形成方式にお
いて、 ホログラム平面と目標物平面間の距離データZ_nと受
波器の位置データX_iから振幅補正関数√(X_i^
2+Z_n^2)を算出する振幅補正関数算出手段と、
受波したホログラムデータS_H(X_i、O)を、振
幅補正関数算出手段の出力により振幅補正し、振幅補正
されたホログラムデータS_H(X_i、O)√(X_
i^2+Z_n^2)を得る振幅補正手段と、目標物に
送波する音波の周波数データfと音速データVから波数
K(=2πf/V)を算出する波数算出手段と、 振幅補正関数算出手段の出力と波数算出手段の出力から
位相補正関数e^j^K^√^{^(X_j)^2^+
^(Z_n)^2^}を算出する位相補正関数算出手段
と、 振幅補正手段の出力する振幅補正されたホログラムデー
タを、位相補正関数算出手段の出力により位相補正し、
振幅補正及び位相補正されたホログラムデータS_H(
X_i、O)√X_i^2+Z_n^2e^j^K^√
^{^(X_i)^2^+^(Z_n)^2^}を得る
位相補正手段と、受波器の位置データと振幅補正関数算
出手段の出力から、新たな変数τ_i(=Xi/√(X
_i^2+Z_n^2)を求める変数変換手段と、 波数算出手段、位相補正手段及び変数変換手段の各出力
並びに目標物平面座標上のサンプリング位置データX_
nを受取り、位相補正手段の出力にフーリエ逆変換を施
し、目標物像データρ(X_n、Zn)=Σ^N_(_
i_=_1_)■_H(τi、O)e^−^j^K^(
X_n)^(τ_i)(ただし、■_H(τ_i、O)
は位相補正手段の出力と同値で該出力をτ_iで表わし
たデータ)を得るフーリエ逆変換手段とを具備すること
を特徴とする映像形成方式。[Claims] A method of transmitting sound waves to a target object, receiving the reflected waves by a receiver array consisting of N receivers, and reproducing the image of the target object based on the received actually measured hologram data. In an image forming method that performs this, an amplitude correction function √(X_i^
2+Z_n^2);
The amplitude of the received hologram data S_H(X_i, O) is corrected by the output of the amplitude correction function calculating means, and the amplitude-corrected hologram data S_H(X_i, O)√(X_
i^2+Z_n^2), a wave number calculating means for calculating a wave number K (=2πf/V) from frequency data f and sound velocity data V of a sound wave transmitted to a target object, and an amplitude correction function calculating means. From the output of
a phase correction function calculation means for calculating ^(Z_n)^2^^; and a phase correction of the amplitude-corrected hologram data output from the amplitude correction means by the output of the phase correction function calculation means;
Amplitude-corrected and phase-corrected hologram data S_H (
X_i, O)√X_i^2+Z_n^2e^j^K^√
A new variable τ_i (=Xi/√ (X
_i^2+Z_n^2), each output of the wave number calculation means, phase correction means, and variable conversion means, and sampling position data X_ on the target plane coordinates.
n, performs inverse Fourier transform on the output of the phase correction means, and obtains target object image data ρ(X_n, Zn)=Σ^N_(_
i_=_1_)■_H(τi, O)e^-^j^K^(
X_n)^(τ_i) (However, ■_H(τ_i, O)
and an inverse Fourier transform means for obtaining data whose output is equal to the output of the phase correction means and represented by τ_i.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP9012986A JPS62247391A (en) | 1986-04-21 | 1986-04-21 | Image forming system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP9012986A JPS62247391A (en) | 1986-04-21 | 1986-04-21 | Image forming system |
Publications (2)
Publication Number | Publication Date |
---|---|
JPS62247391A true JPS62247391A (en) | 1987-10-28 |
JPH0217115B2 JPH0217115B2 (en) | 1990-04-19 |
Family
ID=13989898
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP9012986A Granted JPS62247391A (en) | 1986-04-21 | 1986-04-21 | Image forming system |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPS62247391A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8701491B2 (en) | 2008-04-25 | 2014-04-22 | Stichting Voor De Technische Wetenschapen | Acoustic holography |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59208572A (en) * | 1983-04-29 | 1984-11-26 | コミサリヤ・ア・レネルジ・アトミク | Method and apparatus for sonic holography using converged ultrasonic beam |
-
1986
- 1986-04-21 JP JP9012986A patent/JPS62247391A/en active Granted
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59208572A (en) * | 1983-04-29 | 1984-11-26 | コミサリヤ・ア・レネルジ・アトミク | Method and apparatus for sonic holography using converged ultrasonic beam |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8701491B2 (en) | 2008-04-25 | 2014-04-22 | Stichting Voor De Technische Wetenschapen | Acoustic holography |
Also Published As
Publication number | Publication date |
---|---|
JPH0217115B2 (en) | 1990-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPH01207042A (en) | Method and apparatus for reducing repeated adaptation of phase differential effect | |
CN109581275B (en) | Two-dimensional underwater DOA estimation method and device based on non-circular signal and three-dimensional orthogonal array | |
US4484477A (en) | Variable delay memory system | |
CN113702909A (en) | Sound source positioning analytic solution calculation method and device based on sound signal arrival time difference | |
JP3367462B2 (en) | Active sonar and target detection method thereof | |
JP6696575B2 (en) | Moving target detecting system and moving target detecting method | |
Cho et al. | Adaptive near-field beamforming techniques for sound source imaging | |
Thong-un et al. | A linearization-based method of simultaneous position and velocity measurement using ultrasonic waves | |
JPS62247391A (en) | Image forming system | |
US5029144A (en) | Synthetic aperture active underwater imaging system | |
Ma et al. | Low-frequency sound source localization in enclosed space based on time reversal method | |
Sampson et al. | FPGA implementation of low-power 3D ultrasound beamformer | |
Johansson et al. | Acoustic direction of arrival estimation, a comparison between root-music and SRP-PHAT | |
CN112263277B (en) | Ultrasonic Doppler blood flow imaging method, device, equipment and computer medium | |
JPS6391684A (en) | Video forming system | |
CN110613476B (en) | Ultrasonic signal processing method and device | |
WO1996024053A1 (en) | Acoustic probing apparatus | |
KR101956456B1 (en) | Method of transmitting and receiving ultrasound signal using convex probe capable of securing wide imaging area, apparatus using the same, and system using the same | |
JPH0824258A (en) | Ultrasonic transmitting and receiving method and apparatus using dispersion compression method | |
Suzuki et al. | 2.5 dimensional sound field reproduction based on mode matching and equivalent sources considering primary reflections | |
Hasegawa et al. | Self-Localization of Single Microphone from Ultrasound Field Generated by Sources with Subwavelength Displacement | |
CN115825968B (en) | Side-scan image generation method and system and side-scan sonar device | |
CN103126720B (en) | The ultrasonic imaging method of many spread speeds and device | |
JP3136635B2 (en) | Exploration method and device | |
JPH0560B2 (en) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
EXPY | Cancellation because of completion of term |