JPS6217654A - Image pick-up apparatus - Google Patents

Image pick-up apparatus

Info

Publication number
JPS6217654A
JPS6217654A JP60155136A JP15513685A JPS6217654A JP S6217654 A JPS6217654 A JP S6217654A JP 60155136 A JP60155136 A JP 60155136A JP 15513685 A JP15513685 A JP 15513685A JP S6217654 A JPS6217654 A JP S6217654A
Authority
JP
Japan
Prior art keywords
spectrum
fourier transform
memory
imaging
measuring device
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
Application number
JP60155136A
Other languages
Japanese (ja)
Other versions
JPH0664019B2 (en
Inventor
Fuminobu Takahashi
高橋 文信
Masahiro Koike
正浩 小池
Satoshi Ogura
聰 小倉
Izumi Yamada
泉 山田
Kazunori Koga
古賀 和則
Tadashi Sonobe
園部 正
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.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP60155136A priority Critical patent/JPH0664019B2/en
Publication of JPS6217654A publication Critical patent/JPS6217654A/en
Publication of JPH0664019B2 publication Critical patent/JPH0664019B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Length-Measuring Devices Using Wave Or Particle Radiation (AREA)

Abstract

PURPOSE:To make it possible to display a sharp detailed examination image, by allowing a measuring apparatus to determine the resolving power inherent to said apparatus and converting a measured spectrum corresponding to resolving power. CONSTITUTION:The contents of a reference shape memory 1 and a spectrum memory 2 having a measured result stored therein receive functional conversion by a Fourier transform device 9 corresponding to outputs '1', '0' of a change-over controller 3 to be separated into a real number part and imaginary number part while both parts are stored in a reference spectrum memory 10 and a gamma-spectrum memory 11 and the content of the memory 11 is further stored in an R-spectrum memory 12 through a spectrum divider 14. The contents of the memories 11, 12 are stored in an I-spectrum memory 13 through the divider 14 and subsequently converted by an inverse Fourier transform device 15, operated by a power operator 17 through a memory 16, converted by an address converter 19 through a power spectrum memory 18, again converted by a power converter 21 through a shift spectrum memory 20 and processed by an address converter 22 to be applied to an image display device 23 as a vertical/ horizontal deflection signal with predetermined voltage and an image having no blur is displayed.

Description

【発明の詳細な説明】 〔発明の利用分野〕 本発明は、撮像装置に係り、特に分解能の低い測定装置
、例えば、超音波探傷装置、放射線CT表装置レーダ撮
像装置などの映像化を目的とした測定装置に接続するも
ので、高解像の映像を表示するのに好適な撮像装置に関
する。
[Detailed Description of the Invention] [Field of Application of the Invention] The present invention relates to an imaging device, and is particularly aimed at imaging a measuring device with a low resolution, such as an ultrasonic flaw detection device, a radiation CT table device, or a radar imaging device. The present invention relates to an imaging device that is connected to a measuring device and is suitable for displaying high-resolution images.

〔発明の背景〕[Background of the invention]

従来の撮像装置は、測定装置の空間分解能が正規分布で
近似できる測定装置についてのみ適用でき、測定された
スペクトルを逆ワイアストラス変換により分解能を向上
させたスペクトルを導出していた。このため、空間分解
能が未知、あるいは、正規分布と異なる分解能を有する
測定装置には使用することが出来ないという問題点があ
った。
Conventional imaging devices can only be applied to measuring devices whose spatial resolution can be approximated by a normal distribution, and a spectrum with improved resolution is derived from the measured spectrum by inverse Wire Strass transformation. For this reason, there was a problem in that it could not be used for measurement devices with unknown spatial resolution or resolutions different from normal distribution.

尚、これに関連する先行技術としては、特頴昭49−4
2228号明細書に開示されたものが存在する。
In addition, as prior art related to this, Tokkokuho 49-4
There is one disclosed in the No. 2228 specification.

〔発明の目的〕[Purpose of the invention]

本発明は上記した従来技術の問題点に鑑みなされたもの
で、測定装置に接続し、測定装置固有の分解能を装置自
身で把握させ、分解能によって解像度の劣化した測定ス
ペクトルを自動的に解像度の高いスペクトルに変換して
、鮮明で詳細な検査映像を表示できる撮像装置を提供す
ることを目的としている。
The present invention was made in view of the above-mentioned problems of the prior art.The present invention is connected to a measuring device, allows the device itself to grasp the resolution specific to the measuring device, and automatically converts a measured spectrum whose resolution has deteriorated depending on the resolution into a high-resolution one. The purpose of the present invention is to provide an imaging device that can display a clear and detailed inspection image by converting it into a spectrum.

〔発明の概要〕[Summary of the invention]

本発明の撮像装置は、既知形状の試験体を撮像し、測定
した試験体の形状のフーリエ変換関数を真の試験体の形
状のフーリエ変換関数で除することによって、11!l
定装置の分解能をフーリエ変換関数として保持し1次に
測定対象たる試験体を撮像して測定した形状のフーリエ
変換関数を上記フーリエ変換関数で除し、これによって
得られたスペクトルを逆フーリエ変換し、更に逆フーリ
エ変換によって得られた変換スペクトルのパワースペク
トルを該スペクトル幅を中心として座標を折り返したス
ペクトルに変換し、このスペクトルを画像として表示す
ることを特徴としている。
The imaging device of the present invention images a test specimen with a known shape and divides the Fourier transform function of the measured shape of the test specimen by the Fourier transform function of the true shape of the test specimen, thereby obtaining 11! l
The resolution of the fixed device is maintained as a Fourier transform function, and the Fourier transform function of the shape measured by first imaging the test object to be measured is divided by the above Fourier transform function, and the spectrum obtained thereby is inversely Fourier transformed. The present invention is further characterized in that the power spectrum of the transformed spectrum obtained by inverse Fourier transform is transformed into a spectrum whose coordinates are folded around the spectral width, and this spectrum is displayed as an image.

即ち、本発明の撮像装置は、測定した空間スペクトルの
フーリエ変換関数と真の空間スペクトルのフーリエ変換
関数との比が、測定装置固有の空間分解能のフーリエ変
換関数に一致することに着目し、既知形状の試験体を測
定した形状と真の形状から測定装置の空間分解能を表わ
すフーリエ変換関数を決定し、得られた分解能のフーリ
エ変換関数を用いて、未知の試験体を測定した得たぼけ
のあるまたは測定誤差によるバラツキの多いスペクトル
から、高分解能のスペクトルを表示するものである。
That is, the imaging device of the present invention focuses on the fact that the ratio of the Fourier transform function of the measured spatial spectrum to the Fourier transform function of the true spatial spectrum matches the Fourier transform function of the spatial resolution specific to the measuring device, and A Fourier transform function representing the spatial resolution of the measuring device is determined from the measured shape and the true shape of the test specimen, and the Fourier transform function of the obtained resolution is used to calculate the blur obtained by measuring the unknown test specimen. This is to display a high-resolution spectrum from a spectrum that has a lot of variation due to measurement errors.

次に、本発明の原理について説明する。内部欠陥、ある
いは内部構造の形状を画像で表示する非破壊検査装置、
たとえば、超音波探傷機器、放射線C,Tやマイクロ波
レーダ探査装置には、超音波、放射線あるいはマイクロ
波などの空間的なひろがり、それらのセンサ受信部の開
口の大きさにより、画像にぼけが生じる。このぼけは、
測定装置固有の空間分解能とみなせる。空間分解能は、
煮物体も画像で表示した時、画像の空間的なひろがりで
ある。装置固有の空間分解能を把握し、その空間分解能
によって生じた画像のぼけを修正して真の形状に近い画
像を表示するのが、本発明の骨子である。まず、空間分
解能が既知であれば、ぼけを有する画像からぼけのない
画像を導出する原理について述べる。
Next, the principle of the present invention will be explained. Non-destructive inspection equipment that displays images of internal defects or the shape of internal structures;
For example, in ultrasonic flaw detection equipment, radiation C, T, and microwave radar detection devices, images may be blurred due to the spatial spread of ultrasonic waves, radiation, or microwaves, and the size of the aperture of the sensor receiving section. arise. This blur is
It can be regarded as the spatial resolution specific to the measuring device. The spatial resolution is
When a boiled object is also displayed as an image, it is the spatial spread of the image. The gist of the present invention is to grasp the spatial resolution unique to the device, correct the blurring of the image caused by the spatial resolution, and display an image close to the true shape. First, we will describe the principle of deriving an unblurred image from a blurred image if the spatial resolution is known.

まず簡単のために、1次元のぼけを修正する数学的な処
理を考える。実際の画像は2次元であるが、以下に述べ
る1次元のぼけを修正する数学的な処理は、2次元の場
合にも全く同様に適用できるものである。なぜなら、超
音波探傷機器等の測定装置は一次元の画像信号を出力し
、主走査方向と副走査方向の走査により2次元の画像を
得る様に構成されているからである。
First, for the sake of simplicity, we will consider a mathematical process for correcting one-dimensional blur. Although the actual image is two-dimensional, the mathematical processing for correcting one-dimensional blur described below can be applied in exactly the same way to a two-dimensional image. This is because a measuring device such as an ultrasonic flaw detector is configured to output a one-dimensional image signal and obtain a two-dimensional image by scanning in the main scanning direction and the sub-scanning direction.

今、測定装置からの画像信号を便宜上スペクトルと呼ぶ
ことにする。スペクトルの横軸は空間位置、縦軸は強度
とする。測定装置で測定したスペクトル(ぼけを有する
もの)をスペクトル0と呼び第2図に示す。一方、装置
のぼけを表わす空間分解能は、スペクトルRと呼び、第
3図に示す。
For convenience, the image signal from the measuring device will be referred to as a spectrum. The horizontal axis of the spectrum is spatial position, and the vertical axis is intensity. The spectrum (with blur) measured by the measuring device is called spectrum 0 and is shown in FIG. On the other hand, the spatial resolution representing the blur of the device is called spectrum R and is shown in FIG.

第2図に示すスペクトルOおよび第3図のスペクトルR
をそれぞれ関数0(X)、R(X)で表わし、ぼけのな
い時のスペクトルをスペクトル■と呼び関数I (X)
で表わすと次の式(1)の関係がある。
Spectrum O shown in Figure 2 and spectrum R shown in Figure 3
are expressed by the functions 0(X) and R(X), respectively, and the spectrum when there is no blur is called the spectrum ■ and the function I (X)
When expressed as , there is a relationship expressed by the following equation (1).

ここで、未知の関数I (X)を導出するには、以下に
説明する関係を利用して導出する。
Here, the unknown function I (X) is derived using the relationship described below.

まず、式(1)の両辺をフーリエ変換する。すなわち、
左辺は、 と定義する。ここで、jは虚数、ωは周波数である。同
様に、I (X)のフーリエ変換関数をT(ω)。
First, both sides of equation (1) are Fourier transformed. That is,
The left side is defined as . Here, j is an imaginary number and ω is a frequency. Similarly, the Fourier transform function of I (X) is T(ω).

R(X)のフーリエ変換関数をR(ω)とし、下式(3
)、(4)で定義する。
Let the Fourier transform function of R(X) be R(ω), and use the following formula (3
), defined in (4).

式(1)の右辺をフーリエ変換すると次の様になる。When the right side of equation (1) is Fourier transformed, it becomes as follows.

・・・(5) ところで、関数f(X−xo)のフーリエ変換は次の様
に求まる。すなわち、X−x、をtで変数変換し、 ここで、f(ω)は!<  )のフーリエ変換関数であ
る0式(6)の関係から1式(5)中の〔〕内の部分は
次の様になる。
...(5) By the way, the Fourier transform of the function f(X-xo) can be found as follows. That is, convert X−x into variables using t, where f(ω) is! From the relationship of equation (6), which is the Fourier transform function of

・・・(7) 従って、式(5)は、 =R(ω)・I(ω)          ・・・(8
)に変形できる。
...(7) Therefore, equation (5) is: =R(ω)・I(ω) ...(8
) can be transformed into

式(1)の関係が成り立ているとき、式(8)で明らか
な如く、フーリエ変換関数の積で表わす次の式(9)の
関係を得る。
When the relationship in Equation (1) holds true, as is clear from Equation (8), the following relationship in Equation (9) expressed by the product of Fourier transform functions is obtained.

従って、0(X)、R(X)が既知のスペクトルであれ
ば、未知のスペクトルエ(X)は、次式のフーリエ変換
関数の形で導出できる。
Therefore, if 0(X) and R(X) are known spectra, the unknown spectrum E(X) can be derived in the form of the Fourier transform function of the following equation.

さて、フーリエ変換関数I(ω)を実際の空間スペクト
ルI()に戻すには、次式で示す逆フーリエ変換の手順
で可能である。
Now, in order to return the Fourier transform function I(ω) to the actual spatial spectrum I(), it is possible to use the inverse Fourier transform procedure shown in the following equation.

= I (X)             ・・・(1
1)式(10)、 (11)から、I (X)を次式で
求メル。
= I (X) ... (1
1) From formulas (10) and (11), find I (X) using the following formula.

・・・(12) さて、式(12)において、0(ω)は、測定したスペ
クトル0()をフーリエ変換して得られるが、R(ω)
の原関数であるスペクトルR(X)を取得する方法が課
題となる。スペクトルR(X)を測定袋の空間分解能と
すれば、点状の物体を測定装置で画像化し、その画像の
空間的なひろがりをスペクトルR(X)とみなすことが
できる、しかし、現実には、点状の物体は存在し得ない
し、位置によつて空間分解能が変化する可能性もある。
...(12) Now, in equation (12), 0(ω) is obtained by Fourier transforming the measured spectrum 0(), but R(ω)
The challenge is how to obtain the spectrum R(X), which is the original function of . If the spectrum R(X) is the spatial resolution of the measurement bag, it is possible to image a point-like object using a measurement device and consider the spatial spread of the image as the spectrum R(X). However, in reality, , point-like objects cannot exist, and the spatial resolution may change depending on the position.

そこで、形状が既知の物体を測定装置で画像化し、真の
形状の画像との比較から、測定装置の空間分解能を導出
する。これによって導出した分解能を用いて、式(12
)により、物体の真の画像を導出すればよい。
Therefore, an object whose shape is known is imaged using a measuring device, and the spatial resolution of the measuring device is derived from a comparison with an image of the true shape. Using the resolution derived from this, equation (12
) to derive the true image of the object.

まず、分解能、すなわち測定装置のスペクトルRを導出
するには、真の形状を表わすスペクトルエ6、既知の物
体の画像を表わすスペクトル0゜とすれば1式(9)の
関係より、 R(ω)=O,(ω)/Ia(ω)     ・・・(
13)を得る。ここでR(ω)はスペクトルRの、0.
(ω)はスペクトル0.の、およびIa(ω)はスペク
トル1.の、それぞれフーリエ変換スペクトルである。
First, to derive the resolution, that is, the spectrum R of the measuring device, if the spectrum E representing the true shape is 6, and the spectrum representing the image of a known object is 0°, then from the relationship of Equation 1 (9), R(ω )=O, (ω)/Ia(ω) ...(
13) is obtained. Here, R(ω) is 0.
(ω) is the spectrum 0. , and Ia(ω) are spectra 1. are respectively Fourier transform spectra.

この様にすれば、測定装置の空間分解能をR(ω)の形
で導出でき、以後の測定で得たぼけのある画像を、R(
ω)を用いて式(12)の手順で修正し真の画像を表わ
すスペクトルエ(X)を算出できる。
In this way, the spatial resolution of the measuring device can be derived in the form of R(ω), and the blurred images obtained in subsequent measurements can be calculated using R(ω).
ω) can be used to correct the procedure of equation (12) and calculate the spectrum E(X) representing the true image.

第1図は、本発明の第1の実施例を示すブロック図であ
る。
FIG. 1 is a block diagram showing a first embodiment of the present invention.

先ず、測定装置のスペクトルi(ω)を求める場合につ
いて、第1図に示す実施例の動作について説明する。測
定装置のR(ω)スペクトルを求めるときは、切り替え
制御器3の出力を1411+とする。
First, the operation of the embodiment shown in FIG. 1 will be described in the case of determining the spectrum i(ω) of the measuring device. When obtaining the R(ω) spectrum of the measuring device, the output of the switching controller 3 is set to 1411+.

これによって、スイッチ5,6.7の端子aとCが導通
状態になる0次に5選択指示回路24の出力を“1″と
すると、アンド回路25の出方がパ1”となり、スイッ
チ4の端子aとCが導通状態になる。これで、測定装置
のスペクトルR(ω)を求めることが可能となる。尚、
スイッチ5,6゜7は、切り替え制御器3から“1”が
出方されると、端子すとCが導通状態になる。また、ス
イッチ4は、アンド回路25から“1”が出方される場
合に限って、その端子aとCが導通状態になり。
As a result, the terminals a and C of the switches 5, 6.7 become conductive. When the output of the 0th and 5 selection instruction circuit 24 is set to "1", the output of the AND circuit 25 becomes "P1", and the switch 4 Terminals a and C become conductive.This makes it possible to obtain the spectrum R(ω) of the measuring device.
When "1" is output from the switching controller 3, the terminals C of the switches 5, 6 and 7 become conductive. Further, the terminals a and C of the switch 4 are brought into conduction only when the AND circuit 25 outputs "1".

その他の場合には端子すとCが導通状態になる。In other cases, the terminal C becomes conductive.

今、スイッチ4の端子aとCが導通状態になっているた
め、基準形状メモリ1に記憶されている内容(既知形状
の空間位置を番地に対応させ1強度を記憶内容とする。
Now, since the terminals a and C of the switch 4 are in a conductive state, the contents stored in the reference shape memory 1 (the spatial position of the known shape corresponds to the address and the 1 intensity is set as the stored content).

)が、スイッチ4を介してフーリエ変換器9に入力され
、前記した式(3)に従って、フーリエ変換関数工(ω
)が算出される。
) is input to the Fourier transformer 9 via the switch 4, and the Fourier transform function (ω
) is calculated.

フーリエ変換関数工(ω)は、スイッチ5を介して、基
準スペクトルメミリ10の実数部と虚数部に格納される
The Fourier transform function (ω) is stored in the real part and imaginary part of the reference spectrum memory 10 via the switch 5.

次に、選択指示器24の出力を“0”にして、基準形状
の試料を測定して得られるスペクトルを測定装置より入
力し、Ws定スペクトルメモリ2に格納する。8I定ス
ペクトルメモリ2の内容は、スイッチ4を介してフーリ
エ変換器9に入力され、ここで式(4)に従ってスペク
トル0を算出し、実数部、虚数部に分けてOスペクトル
メモリ11に格納する。0スペクトルメモリ11の内容
は、ベクトル除算器14のボートAに入力され、また、
基準スペクトルメモリ10の内容はスイッチ6を介して
ベクトル除算器14のボートBに入力される。ベクトル
除算器14では、ボートAに入力された値A、(実数部
)、A1(虚数部)とボートBに入力された値B、(実
数部)、B、(虚数部)を用いて次式に従って実数値c
2.虚数値c1.を算出する。
Next, the output of the selection indicator 24 is set to "0", and the spectrum obtained by measuring the sample of the reference shape is input from the measuring device and stored in the Ws constant spectrum memory 2. The contents of the 8I constant spectrum memory 2 are inputted to the Fourier transformer 9 via the switch 4, where the spectrum 0 is calculated according to equation (4), and the result is divided into a real part and an imaginary part and stored in the O spectrum memory 11. . The contents of the zero spectrum memory 11 are input to port A of the vector divider 14, and
The contents of the reference spectrum memory 10 are input to port B of the vector divider 14 via switch 6. The vector divider 14 uses the values A, (real part), A1 (imaginary part) input to boat A and the values B, (real part), B, (imaginary part) input to boat B to calculate the following: The real value c according to the formula
2. Imaginary value c1. Calculate.

上式で算出したC、、Ctは、スイッチ7を介してそれ
ぞれRスペクトルメモリ12の実数部。
C, , Ct calculated by the above formula are sent to the real part of the R spectrum memory 12 via the switch 7.

虚数部に格納される。これによって、測定装置の空間が
分解能に関するスペクトルRが得られる。
Stored in the imaginary part. This results in a spectrum R related to the spatial resolution of the measuring device.

次に、切り替え制御器3の出力を“0”にする。Next, the output of the switching controller 3 is set to "0".

そして、未知の試料を測定した結果を測定装置からスイ
ッチ4を介してフーリエ変換器9に入力し。
Then, the results of measuring the unknown sample are input from the measuring device to the Fourier transformer 9 via the switch 4.

式(2)に従って、スペクトルざを算出し、実数部。According to equation (2), calculate the spectral ratio and calculate the real part.

虚数部に分け、スイッチ5を介してOスペクトルメモリ
11に記録する。0スペクトルメモリ11の内容は、ス
イッチ6を介してベクトル除算器14のボートAに、R
スペクトルメモリ12の内容はスイッチ6を介して除算
器14のボートBに入力され、除算器14で式(14)
、 (IS)に従った演算を実行し、その結果C,,C
,をスペクトルエとして、スイッチ7を介してニスベク
トルメモリ13の実数部および虚数部にそれぞれ記録す
る。
It is divided into imaginary parts and recorded in the O spectrum memory 11 via the switch 5. The contents of the zero spectrum memory 11 are transferred to the port A of the vector divider 14 via the switch 6.
The contents of the spectrum memory 12 are input to the port B of the divider 14 via the switch 6, and the divider 14 calculates the equation (14).
, executes the operation according to (IS), and the result is C,,C
, are recorded in the real part and imaginary part of the varnish vector memory 13 via the switch 7, respectively, as a spectrum.

■スペクトルメモリ13の内容は、逆フーリエ変換器1
5に入力され、ここで、式(11)に従って。
■The contents of the spectrum memory 13 are stored in the inverse Fourier transformer 1.
5, where according to equation (11).

スペクトルIが算出され、その結果は、Iスペクトルメ
モリ16内に実数部、虚数部にわけて格納される。パワ
ー演算器17は、■スペクトルメモリ16の実数部の値
I1、虚数部の値I、を用いて、I、を下式に従って算
出する。
The spectrum I is calculated, and the result is stored in the I spectrum memory 16 separately into a real part and an imaginary part. The power calculator 17 uses the real part value I1 and the imaginary part value I of the spectrum memory 16 to calculate I according to the following formula.

I、=(I、”+Iz”)”        ・・・(
16)この値は、パワースペクトルメモリ18に記録さ
れる。
I,=(I,"+Iz")"...(
16) This value is recorded in the power spectrum memory 18.

番地変換器19は、メモリ18の内容を番地をシフトし
、シフトスペクトルメモリ20に配列しなおして記録す
る1番地のシフトは次の様に行なわれる0番地が2Nま
であり、パワースペクトルメモリ18の番地にの内容を
エア(k)で表わすと、次式に従って番地mのメモリ2
0の内容1.、(m)が決定する。
The address converter 19 shifts the contents of the memory 18 by address and rearranges and records them in the shift spectrum memory 20.The shift of the 1st address is performed as follows. If the content at the address is expressed as air (k), then the memory 2 at address m according to the following formula:
Contents of 01. , (m) are determined.

・・・(17) パワー変換器21は、番地mのメモリ20の内容I、、
(m)に比例した電圧値vvを1番地換算器22はメモ
リ20の番地mに比例した電圧値vxをそれぞれ画像表
示器23の垂直偏向および水平偏向信号として出力し、
画像表示器23において。
...(17) The power converter 21 stores the contents I of the memory 20 at address m, .
(m), the 1 address converter 22 outputs the voltage value vx proportional to the address m of the memory 20 as a vertical deflection signal and a horizontal deflection signal of the image display 23, respectively.
In the image display 23.

ぼけのない画像が表示される。An unblurred image is displayed.

次に、物体の厚み放射線透過法で測定し、これによって
得られたスペクトルを第1図に示す装置を用いて処現し
た例について述べる。第1!Iは、放射線透過法による
試験体厚さ測定装置の一例を示す図である。放射線源1
00から放射される放射線101はコリメータ103で
コリメートされ。
Next, an example will be described in which the thickness of an object is measured by a radiographic transmission method and the spectrum obtained thereby is processed using the apparatus shown in FIG. 1st! I is a diagram showing an example of a test specimen thickness measuring device using a radiographic method. radiation source 1
Radiation 101 emitted from 00 is collimated by a collimator 103.

検出器104で検出される。線源100とコリメータ1
03を備えた検出器104は、X軸方向に並進走査され
る。そして、線源100と検出器104の間に試験体1
02が配置される。厚み測定装!!105は、走査軸X
における試験体102のY軸方向厚みを記録する。走査
位WtXでの放射線101検出器104による検出強度
をφ(X)とする、試験体102がない場所の検出強度
をφ。
It is detected by the detector 104. Source 100 and collimator 1
03 is translated and scanned in the X-axis direction. The test object 1 is placed between the radiation source 100 and the detector 104.
02 is placed. Thickness measuring device! ! 105 is the scanning axis
The thickness of the test specimen 102 in the Y-axis direction is recorded. Let φ(X) be the detection intensity by the radiation 101 detector 104 at the scanning position WtX, and φ is the detection intensity at a place where there is no test object 102.

とし、試験体102のY軸方向厚みをd()とすると、
検出強度φ(X)は次の様に表わせる。
If the thickness of the test specimen 102 in the Y-axis direction is d(),
The detection intensity φ(X) can be expressed as follows.

φ(X)=φ*exp(−p d (X ))    
 −(18)ここで、μは試験体102の放射線101
に対する線吸収係数である。
φ(X)=φ*exp(-p d (X ))
-(18) Here, μ is the radiation 101 of the specimen 102
is the linear absorption coefficient for

式(18)の関係を変形すると次式の様になる。When the relationship in equation (18) is transformed, it becomes as shown in the following equation.

* (X)/ $a=exp(−p d (X))  
  −(19)式(19)の両辺の対数をとると d (X)= −Q n (+6/ φ(X))   
 −(20)μ を得る。厚み測定装置は、式(20)の関係を用い、透
過放射線の検出強度φ(X)に基づいて、厚み分布d(
’)を測定する。
* (X)/ $a=exp(-p d (X))
-(19) Taking the logarithm of both sides of equation (19), d (X) = -Q n (+6/ φ(X))
−(20) μ is obtained. The thickness measuring device uses the relationship of equation (20) to calculate the thickness distribution d(
') to measure.

!@5図(a)は、試験体102を測定して得られた厚
み分布を示す図であり、縦軸dが厚みを示し。
! @5 Figure (a) is a diagram showing the thickness distribution obtained by measuring the test specimen 102, and the vertical axis d indicates the thickness.

横軸Xが位置を示している。第5図(b)は、試験体1
02の真の厚み分布を示す図である。第5図(b)は、
試験体102の真の厚み分布を示す図である。第5図(
a)に示す測定により得られた厚み分布は、第5図(b
)に示す真の厚み分布と比較して、なめらかな分布、即
ちぼけのある分布となっている。
The horizontal axis X indicates the position. Figure 5(b) shows test specimen 1.
02 is a diagram showing the true thickness distribution. Figure 5(b) shows
3 is a diagram showing the true thickness distribution of the test specimen 102. FIG. Figure 5 (
The thickness distribution obtained by the measurement shown in a) is shown in Figure 5 (b).
) is a smooth distribution, that is, a blurred distribution compared to the true thickness distribution shown in ().

第4図に示す測定装置の分解能は、第6図(a)に示す
分布となり、この分布をスペクトルRとして、第1図に
示す装置を用いてぼけのない厚み分布を画像表示すると
、第6図(b)に示すスペクトル分布が得られる。第5
図(b)に示す真の厚み分布と第6図(b)に示すスペ
クトルをくらべると。
The resolution of the measuring device shown in FIG. 4 becomes the distribution shown in FIG. A spectral distribution shown in Figure (b) is obtained. Fifth
Comparing the true thickness distribution shown in Figure (b) with the spectrum shown in Figure 6 (b).

はとんど同じ厚み分布が得られていることがわかる。It can be seen that almost the same thickness distribution is obtained.

さて、第1図で示した装置では、測定装置が測定したス
ペクトルに測定値のバラツキが大きいと、そのバラツキ
を強調してしまう傾向がある。従って、測定スペクトル
中の測定誤差によるバラツキを低下させて、鮮明なスペ
クトルエを導出することが必要となる。測定スペクトル
中に測定誤差によるバラツキが多く含まれる場合、第1
図で示す装置で得られるスペクトルエの像は、第7図で
示すスペクトルの様になる。第7図から明らかな様に、
特定周波数成分がリップル状に発生している。
Now, in the apparatus shown in FIG. 1, if there is a large variation in the measured values in the spectrum measured by the measurement apparatus, there is a tendency to emphasize the variation. Therefore, it is necessary to reduce the variation due to measurement errors in the measured spectrum and derive a clear spectrum. If the measured spectrum contains many variations due to measurement errors, the first
The spectral image obtained by the apparatus shown in the figure looks like the spectrum shown in FIG. As is clear from Figure 7,
A specific frequency component appears in a ripple pattern.

このノイズを低減する第2の手段は、スペクトルIを平
滑化して、高い空間周波数成分を除去することである。
A second means of reducing this noise is to smooth the spectrum I to remove high spatial frequency components.

第7図に示すスペクトルエに対し次式に示す様な平滑化
処理を行う、つまり、各点XIに関して、その周囲の点
X1.にのスペクトル強度Iap(Xt*h)の平均値
を、新しいスペクトル強度1.、’ (x、)とする、
この処理は1次式による。
Smoothing processing as shown in the following equation is performed on the spectrum shown in FIG. The average value of the spectral intensity Iap(Xt*h) of 1. is the new spectral intensity 1. ,' (x,),
This process is based on a linear equation.

1、、’  (X、)=(Σ I□(X、やk))/(
2Q+ 1 )・・・(21) 式(21)により第7図のスペクトルを平滑化した新し
いスペクトルを第8図に示す、第7図に示すスペクトル
Iから第8図に示すスペクトルを減算し、第9図に示す
スペクトルを得る。第9図のスペクトルは第7図のスペ
クトルに含まれているリップル状ノイズになっている。
1,,' (X,)=(Σ I□(X, or k))/(
2Q+ 1 )...(21) A new spectrum obtained by smoothing the spectrum in FIG. 7 using equation (21) is shown in FIG. 8. Subtracting the spectrum shown in FIG. 8 from spectrum I shown in FIG. 7, The spectrum shown in FIG. 9 is obtained. The spectrum in FIG. 9 is ripple-like noise included in the spectrum in FIG. 7.

従って、第9図のスペクトルをフーリエ変換することに
より、リップル状ノイズの特徴的な周波数を抽出できる
。第10図は、第9図に示すスペクトルをフーリエ変換
したスペクトルの例を示す図である。リップル状ノイズ
の空間周波数は高いので、第10図に示すスペクトルに
おいて、ω軸の中央部に形成されるピークが、リップル
状ノイズに起因するノイズである。そこで、第10図に
示すスペクトルのうちω軸で中央部幅Wの範囲内のスペ
クトル(第11図に示す、)だけを抽出し、第1図に示
すニスベクトルメモリ13に記憶されているスペクトル
Iから第11図のスペクトルを減算すれば、スペクトル
I中で、リップル状ノイズを形成する周波数成分を除去
できる。
Therefore, by Fourier transforming the spectrum shown in FIG. 9, characteristic frequencies of ripple-like noise can be extracted. FIG. 10 is a diagram showing an example of a spectrum obtained by Fourier transforming the spectrum shown in FIG. 9. Since the spatial frequency of the ripple-like noise is high, in the spectrum shown in FIG. 10, the peak formed at the center of the ω-axis is the noise caused by the ripple-like noise. Therefore, from the spectrum shown in FIG. 10, only the spectrum within the range of the center width W on the ω axis (shown in FIG. 11) is extracted, and the spectrum stored in the varnish vector memory 13 shown in FIG. By subtracting the spectrum shown in FIG. 11 from I, frequency components that form ripple-like noise can be removed from spectrum I.

第12図は、上記した一連の処理(第7図から第11図
に関する処理)を実行する本発明の第2の実施例を示す
ブロック図である。第12図において、上記した一連の
処理を実行するのは、点線枠内の機器であり、点線枠外
の機器は、第1図に示した機器と同一であり、点線枠内
の機器を第1図に示す装置に増設した構成を有している
0図示する様に、ニスベクトルメモリ13と逆フーリエ
変換器15との間にスイッチ201を挿入し、かつパワ
ースペクトルメモリ18と番地変換器19との間にスイ
ッチ202を挿入する。リップル除去指示器200の出
力を“1″にすると、スイッチ201、スイッチ202
の端子Cと端子aとが導通状態となり、リップル除去指
示器200の出力を“0″にするとスイッチ201,2
02の端子Cと端子すが導通になる。まず、リップル除
去指示器200の出力レベルを“1”にする、Iスペク
トルメモリ13の内容はスイッチ201を介して逆フー
リエ変換器15に入力され、以下、第1図の装置で説明
した如く、パワースペクトルメモリ18に第7図で示し
た様なリップルのあるスペクトルが格納される。このス
ペクトルは、スイッチ202を介して平滑演算器203
に入力され、式(21)に従った処理が施され、第8図
に示す平滑化されたスペクトルに変換され平滑スペクト
ルメモリ204に格納される0次に減算器205は。
FIG. 12 is a block diagram showing a second embodiment of the present invention that executes the series of processes described above (processes related to FIGS. 7 to 11). In FIG. 12, the devices inside the dotted line frame execute the series of processes described above, and the devices outside the dotted line frame are the same as those shown in FIG. As shown in the figure, a switch 201 is inserted between the varnish vector memory 13 and the inverse Fourier transformer 15, and a switch 201 is inserted between the power spectrum memory 18 and the address converter 19. A switch 202 is inserted between them. When the output of the ripple removal indicator 200 is set to "1", the switches 201 and 202
Terminal C and terminal a become conductive, and when the output of the ripple removal indicator 200 is set to "0", the switches 201 and 2
Terminal C of 02 and terminal S become conductive. First, the contents of the I-spectrum memory 13, which sets the output level of the ripple removal indicator 200 to "1", are input to the inverse Fourier transformer 15 via the switch 201, and as explained below with respect to the apparatus of FIG. A spectrum with ripples as shown in FIG. 7 is stored in the power spectrum memory 18. This spectrum is transmitted to the smoothing calculator 203 via the switch 202.
The zero-order subtracter 205 processes the spectrum according to equation (21), converts it into a smoothed spectrum shown in FIG. 8, and stores it in the smoothed spectrum memory 204.

パワースペクトルメモリ18の内容、すなわちリップル
を有するスペクトルから平滑スペクトルメモリ204の
内容、すなわち平滑スペクトルを減算し、第9図に示す
様なリップルのスペクトルを減算スペクトルメモリ20
6に格納される。フーリエ変換器207は、減算スペク
トルメモリ206の内容をフーリエ変換し、第10図に
示すリップルの周波数スペクトルを実数部と虚数部にわ
けてノイズスペクトルメモリ208に格納する。ノイズ
スペクトルメモリ208の内容は、第11図に示す様に
、第10図のスペクトルの中央部分の幅Wのスペクトル
だけを保存し、幅W以外は0にしたスペクトルに変換さ
れ、実数部と虚数部とに分けてリップルスペクトルメモ
リ210に格納される。
The contents of the smooth spectrum memory 204, that is, the smooth spectrum, are subtracted from the contents of the power spectrum memory 18, that is, the spectrum with ripples, and the spectrum with ripples as shown in FIG. 9 is subtracted from the spectrum memory 20.
6. The Fourier transformer 207 performs a Fourier transform on the contents of the subtraction spectrum memory 206, divides the frequency spectrum of the ripple shown in FIG. As shown in FIG. 11, the contents of the noise spectrum memory 208 are converted into a spectrum in which only the spectrum with a width W at the center of the spectrum in FIG. The ripple spectrum is stored in the ripple spectrum memory 210 separately.

減算器211は、図示する様に、ニスベクトルメモリ1
3の内容からりップルスペクトルメモリ210の内容を
減じて、リップル除去スペクトルメモリ212に実数部
と虚数部とをわけて格納する0次に、リップル除去指示
器200の出力を110”にすると、リップル除去スペ
クトルメモリ212の内容がスイッチ201を介して逆
フーリエ変換器15に入力され、以下、第1図の装置で
説明した如く、パワースペクトルメモリ18までスペク
トルが格納され、さらにスイッチ202を介して番地変
換器19に格納内容が流れ、リップルの除去された画像
が表示される。
The subtracter 211 is connected to the varnish vector memory 1 as shown in the figure.
Subtract the contents of the ripple spectrum memory 210 from the contents of 3 and store the real part and the imaginary part separately in the ripple removal spectrum memory 212. Next, if the output of the ripple removal indicator 200 is set to 110'', The contents of the ripple-removed spectrum memory 212 are inputted to the inverse Fourier transformer 15 via the switch 201, and the spectra are stored up to the power spectrum memory 18 as explained in connection with the apparatus shown in FIG. The stored contents flow to the address converter 19, and an image with ripples removed is displayed.

次に、第12図に示す装置を、超音波探傷装置による画
像のぼけの修正に用いた場合の具体例について説明する
Next, a specific example will be described in which the apparatus shown in FIG. 12 is used to correct image blur caused by an ultrasonic flaw detector.

第13図は、超音波探傷装置の一例を示す図であり、第
12図に示す装置を用いて、試験体表面の凹凸の映像化
を図るものである。第13図において、探触子300を
走査装置301によりX軸方向に走査しながら、超音波
ビーム307を発信し、水槽302の水中303に設置
した試験体304の表面から反射した超音波を受信する
。探傷器305では、位置における超音波発信がら受信
までの時間間隔t(X)を測定し、次式により試料の厚
さd (X)を測定する。
FIG. 13 is a diagram showing an example of an ultrasonic flaw detection device, and the device shown in FIG. 12 is used to visualize the irregularities on the surface of a test specimen. In FIG. 13, a probe 300 is scanned in the X-axis direction by a scanning device 301 while transmitting an ultrasonic beam 307 and receiving ultrasonic waves reflected from the surface of a test specimen 304 placed underwater 303 in a water tank 302. do. The flaw detector 305 measures the time interval t(X) between ultrasonic transmission and reception at a position, and measures the thickness d(X) of the sample using the following equation.

d (X)= d o  v t (X)/ 2   
   ・・・(22)式(22)においてdoは探触子
300から水槽302の底面までの距離、■は水303
での音速である。超音波ビーム307は、水303中で
ひろがるため、この時の測定系の空間分解能は、超音波
ビーム307の空間的な強度分布に近く、第14図(b
)に示すスペクトルRになる。
d (X) = d ov t (X)/2
(22) In equation (22), do is the distance from the probe 300 to the bottom of the water tank 302, and ■ is the water 303.
The speed of sound is . Since the ultrasonic beam 307 spreads in the water 303, the spatial resolution of the measurement system at this time is close to the spatial intensity distribution of the ultrasonic beam 307, as shown in Fig. 14 (b).
) is the spectrum R shown in FIG.

この空間分解能のため、第14図(a)に示す様な厚さ
分布を有する試験体304を測定すると、探傷器305
から第14図(c)に示すスペクトルが得られる。第1
4図(c)に示すスペクトルを用いて、第1に示す装置
で映像化を図ると、第15図(a)に示す様にリップル
の多い厚み分布になる。
Due to this spatial resolution, when measuring the specimen 304 having the thickness distribution as shown in FIG. 14(a), the flaw detector 305
From this, the spectrum shown in FIG. 14(c) is obtained. 1st
When the spectrum shown in FIG. 4(c) is visualized with the apparatus shown in FIG. 1, a thickness distribution with many ripples is obtained as shown in FIG. 15(a).

しかし、第1図に示す装置に第12図に示す装置を増設
すると、第15図(b)に示す様に、リップルのない真
の厚み分布に近いものが得られる。
However, if the device shown in FIG. 12 is added to the device shown in FIG. 1, a thickness distribution close to the true thickness without ripples can be obtained, as shown in FIG. 15(b).

上記した様に、第1図に示す装置に第12図に示す装置
を増設すると、測定装置から出力されるスペクトルにば
らつきが多く存在する場合にも、第15図(b)に示す
様なリップルのない鮮明なスペクトルを導出することが
できる。
As mentioned above, if the device shown in FIG. 12 is added to the device shown in FIG. It is possible to derive a clear spectrum without any

測定したスペクトルがバラツキの大きいスペクトルの場
合、リップルを低減して鮮明なスペクトルを導出する第
2の手段は、つぎのようにすることである。
When the measured spectrum is a spectrum with large variations, a second means of reducing ripples and deriving a clear spectrum is as follows.

X Xを乗じ、第16図(b)に示すスペクトルとし、ΔX を乗じ第16図(c)に示すスペクトルにする。第16
図(b)を実数部のスペクトル、第16図(c)を虚数
部のスペクトルとしてフーリエ変換すれば。
Multiply by X to obtain the spectrum shown in FIG. 16(b), and multiply by ΔX to obtain the spectrum shown in FIG. 16(c). 16th
If Figure 16(b) is the spectrum of the real part and Figure 16(c) is the spectrum of the imaginary part, Fourier transform is performed.

ノイズによる高調波成分の強度を低減することができる
The intensity of harmonic components due to noise can be reduced.

第17図は、上記の処理を実行する装置構成の一例を示
す図である。第17図の点線枠内の装置は、第1図に示
した装置に増設する部分である。
FIG. 17 is a diagram showing an example of the configuration of a device that executes the above processing. The device within the dotted line frame in FIG. 17 is a part that is added to the device shown in FIG.

第1図中の測定スペクトルメモリ2とスイッチ4の間に
挿入する1次に、増設する装置の動作内容を述べる。第
17図において、cos発生器400は、測定スペクト
ルメモリ2の番地に比例した産生器400の値と測定ス
ペクトルメモリ2の内容を乗算し、その値をキャリアス
ペクトルメモリ405の実数部に記録する。−万乗算器
403は。
The operation of the primary and additional equipment inserted between the measurement spectrum memory 2 and the switch 4 in FIG. 1 will be described. In FIG. 17, a cos generator 400 multiplies the value of the generator 400 proportional to the address of the measured spectrum memory 2 by the contents of the measured spectrum memory 2, and records the value in the real part of the carrier spectrum memory 405. - million multiplier 403.

発生器401の値と、測定スペクトルメモリ2の内容を
乗算した結果をキャリアスペクトルメモリ405の内容
は第1図で実数部だけの測定スペクトルの内容にかわっ
て、スイッチ4を介して以下、第1図の装置に入力され
処理される。この時に得られる画像は、第12図の装置
を増設した場合にくらべ、リップルが多少残る場合もあ
るが、第15図(a)に示す様なリップルは確実に低減
できる。
The content of the carrier spectrum memory 405 is the result of multiplying the value of the generator 401 by the content of the measured spectrum memory 2. Instead of the content of the measured spectrum of only the real part in FIG. It is input to the device shown in the figure and processed. The image obtained at this time may have some ripples compared to when the apparatus shown in FIG. 12 is added, but the ripples as shown in FIG. 15(a) can be definitely reduced.

第12図、第17図の装置をともに第1図の装置に増設
すれば、さらに効果が上がる。
If both the devices shown in FIGS. 12 and 17 are added to the device shown in FIG. 1, the effect will be further improved.

以上説明した本発明の装置は、測定装置の固有の空間分
解能を求め、その分解能によってぼけの生じた測定スペ
クトルを、ぼけのないスペクトル、すなわち、さらに高
い分解能を有する測定装置で測定したスペクトルと同等
なスペクトルを得るものである。
The device of the present invention as described above determines the inherent spatial resolution of the measuring device, and converts the measured spectrum that is blurred due to that resolution into a spectrum that is not blurred, that is, equivalent to a spectrum measured with a measuring device that has higher resolution. This will give you a spectrum.

第1図の本発明装置の動作を流れ図で簡単に表わすと、
第18図に示す様になる。第17図の装置を第1図に増
設した場合には、第18図の流れ図において、ステップ
S1の処理のかわりに第19図の流れ図、ステップS2
の処理の代わりに第20図の流れ図を挿入した図になる
The operation of the device of the present invention shown in FIG. 1 can be simply expressed in a flowchart as follows:
The result will be as shown in FIG. When the device shown in FIG. 17 is added to the device shown in FIG. 1, the process in step S1 in the flow chart shown in FIG.
The flowchart shown in FIG. 20 is inserted instead of the process shown in FIG.

第12図の装置を第1図の装置に増設した場合の動作の
流れ図は、第18図に示す流れ図で、ステップS3と8
4の処理の間に第21図の流れ図を挿入した形になる。
The flow chart of the operation when the device shown in FIG. 12 is added to the device shown in FIG. 1 is the flow chart shown in FIG.
The flowchart shown in FIG. 21 is inserted between the processes in step 4.

〔発明の効果〕〔Effect of the invention〕

本発明の撮像装置は、測定装置固有の空間分解能によっ
て劣化した測定スペクトルを測定装置の改造を何んら伴
なわずに、分解能を著じるしく向上させたスペクトルに
変換することが出来る・測定スペクトルに含まれる測定
誤差に起因するバラツキの影響も低減できる機能も実現
できた。
The imaging device of the present invention can convert a measurement spectrum degraded by the spatial resolution inherent to the measurement device into a spectrum with significantly improved resolution without any modification of the measurement device. We were also able to realize a function that can reduce the effects of variations caused by measurement errors contained in spectra.

また、本発明の撮像装置は測定装置固有の空間分解能を
測定するため、分解能が未知のものなど、どんな分解能
を有する測定装置に対しても、撮像装置として広く活用
できる。したがって、これまで分解能が悪いと一般にい
われていた超音波探傷法、レーダ撮像法や放射線撮像法
などで、高解像の映像を表示できる様になる。
Further, since the imaging device of the present invention measures the spatial resolution specific to the measurement device, it can be widely used as an imaging device for measurement devices having any resolution, including those with unknown resolution. Therefore, it is now possible to display high-resolution images using ultrasonic flaw detection methods, radar imaging methods, radiation imaging methods, etc., which were generally said to have poor resolution.

【図面の簡単な説明】[Brief explanation of the drawing]

第1図は本発明の撮像装置の第1の実施例を示すブロッ
ク図、第2図は分解能の低い測定装置で測定して得られ
たスペクトバノ、〕−例を示す図、第3図は測定装置の
空間分解能の一例を示す図、第4図は放射線透過法によ
り試験体厚み分布を測定する装置の一例を示す図、第5
図(a)は、第4図に示す装置で測定した試験体厚み分
布の一例を示す図、第5図(b)は真の厚み分布を示す
図、第6図は第4図に示す装置の空間分解能を示す図、
第6図(b)は第1図に示す第1の実施例を用いて第5
図(、)に示す測定スペクトルより導出した試験体厚み
分布を示す図、第7図は測定して得られたスペクトルに
誤差によるバラツキがあるとき第1図に示す第1の実施
例により得られる出カスベクトルを示す図、第8図は第
7図に示すスペクトルを平滑化したスペクトルを示す図
、第9図は第7図のスペクトルから第8図のスペクトル
を差し引いたリップルノイズのスペクトルを示す図、第
10図は第9図のスペクトルをフーリエ変換したスペク
トルを示す図、第11図は第10図のスペクトル中央部
の幅Wの部分を抽出したスペクトルを示す図、第12図
はリップルノイズを低減するための本発明の撮像装置の
第2の実施例を示す図、第13図は超音波により試料厚
み測定を行なう装置の一例を示す図、第14図(a)は
第13図の条件で測定する試料の真の厚み分布を示す図
、第14図(b)は第13図の測定装置の空間分解能を
示す図、第14図(c)は測定した試料の厚み分布を示
す図、第15図(a)は第1図に示す第1の実施例を用
いて導出した試料の厚み分布を示す図、第15図(b)
は第12図に示す第2の実施例を用いて導出した試料の
厚み分布を示す図、第16図(a)はバラツキの多い測
定スペクトルの一例を示す図、第16図(b)は第16
図(a)に示すスペクトルにcos波を乗じたスペクト
ルを示す図、第16図(Q)は第16図(a)に示すス
ペクトルにsin波を乗じたスペクトルを示す図、第1
7図はリップルノイズを低減するための第3の実施例を
示す図、第18図は第1図に示す第1の実施例の動作を
説明した流れ図、第19図は第17図に示す第2の実施
例に係る装置を第1図に示す第1の実施例に増設したこ
とによって竿18図の流れ図に新たに加わる流れ図、第
20図も第19図と同様に新たに加わる動作の流れ図を
示した図、第21図は第12図に示す第2の実施例に係
る装置を第1図に示す第1の実施例に増設したときに第
18図に新たに加わる動作の流れ図を示した図である。 1・・・基準形状メモリ、2・・・測定スペクトルメモ
リ、3・・・切り替え制御器、4,5,6.7・・・ス
イッチ。 9・・・フーリエ変換器、10・・・基準スペクトルメ
モリ、11・・・Oスペクトロメモリ、12・・・Rス
ペクトルメモリ、13・・・ニスベクトルメモリ、14
・・・ベクトル除算器、15・・・逆フーリエ変換器、
16・・・ニスベクトルメモリ、17・・・パワー演算
器、18・・・パワースペクトルメモリ、19・・・番
地変換器、20・・・シフトスペクトルメモリ、21・
・・パワー換算器、22・・・番地換算器、23・・・
画像表示器、103・・・コリメータ、104・・・検
出器、105・・・厚み測定装置、200・・・リップ
ル除去指示器。 201.202・・・スイッチ、203・・・平滑演算
器、204・・・平滑スペクトルメモリ、205・・・
減算器、206・・・減算スペクトルメモリ、207・
・・フーリエ変換器、208・・・ノイズスペクトルメ
モリ、209・・・周波数ゲート、210・・・リップ
ルスペクトルメモリ、211・・・減算器、212・・
・リップル除去スペクトルメモリ、300・・・探触子
、301X 403・・・乗算器、405・・・キャリアスペクトル
メモリ。
Fig. 1 is a block diagram showing a first embodiment of the imaging device of the present invention, Fig. 2 is a diagram showing an example of spectrovano obtained by measuring with a measuring device with low resolution, and Fig. 3 is a diagram showing measurement Figure 4 shows an example of the spatial resolution of the device.
Figure (a) shows an example of the specimen thickness distribution measured with the apparatus shown in Figure 4, Figure 5 (b) shows the true thickness distribution, and Figure 6 shows the apparatus shown in Figure 4. A diagram showing the spatial resolution of
FIG. 6(b) shows the fifth embodiment using the first embodiment shown in FIG.
Figure 7 shows the specimen thickness distribution derived from the measured spectrum shown in Figure 7, which is obtained by the first embodiment shown in Figure 1 when there is variation due to error in the spectrum obtained by measurement. Figure 8 shows the spectrum shown in Figure 7, smoothed, and Figure 9 shows the ripple noise spectrum obtained by subtracting the spectrum in Figure 8 from the spectrum in Figure 7. Figure 10 shows the spectrum obtained by Fourier transforming the spectrum in Figure 9, Figure 11 shows the spectrum extracted from the width W at the center of the spectrum in Figure 10, and Figure 12 shows ripple noise. FIG. 13 is a diagram showing an example of an apparatus for measuring sample thickness using ultrasonic waves, and FIG. 14(a) is a diagram showing a second embodiment of the imaging device of the present invention for reducing Figure 14(b) is a diagram showing the true thickness distribution of the sample measured under the conditions, Figure 14(b) is a diagram showing the spatial resolution of the measuring device in Figure 13, and Figure 14(c) is a diagram showing the thickness distribution of the measured sample. , FIG. 15(a) is a diagram showing the thickness distribution of the sample derived using the first example shown in FIG. 1, FIG. 15(b)
is a diagram showing the thickness distribution of the sample derived using the second example shown in FIG. 12, FIG. 16
A diagram showing a spectrum obtained by multiplying the spectrum shown in FIG. 16(a) by a cosine wave, FIG. 16(Q) is a diagram showing a spectrum obtained by multiplying the spectrum shown in FIG.
7 is a diagram showing the third embodiment for reducing ripple noise, FIG. 18 is a flowchart explaining the operation of the first embodiment shown in FIG. 1, and FIG. 19 is a diagram showing the operation of the first embodiment shown in FIG. 17. By adding the device according to the second embodiment to the first embodiment shown in FIG. 1, a new flowchart is added to the flowchart in FIG. 18, and FIG. 20 is a flowchart of newly added operations similar to FIG. FIG. 21 shows a flowchart of operations newly added to FIG. 18 when the device according to the second embodiment shown in FIG. 12 is added to the first embodiment shown in FIG. This is a diagram. DESCRIPTION OF SYMBOLS 1... Reference shape memory, 2... Measured spectrum memory, 3... Switching controller, 4, 5, 6.7... Switch. 9... Fourier transformer, 10... Reference spectrum memory, 11... O spectro memory, 12... R spectrum memory, 13... Varnish vector memory, 14
...vector divider, 15...inverse Fourier transformer,
16... Varnish vector memory, 17... Power calculator, 18... Power spectrum memory, 19... Address converter, 20... Shift spectrum memory, 21...
...Power converter, 22...Address converter, 23...
Image display device, 103... collimator, 104... detector, 105... thickness measuring device, 200... ripple removal indicator. 201.202...Switch, 203...Smoothing calculator, 204...Smoothing spectrum memory, 205...
Subtractor, 206... Subtraction spectrum memory, 207.
...Fourier transformer, 208...Noise spectrum memory, 209...Frequency gate, 210...Ripple spectrum memory, 211...Subtractor, 212...
- Ripple removal spectrum memory, 300... Probe, 301X 403... Multiplier, 405... Carrier spectrum memory.

Claims (1)

【特許請求の範囲】 1、既知形状の試験体を測定装置で撮像して得られた形
状を示す信号をフーリエ変換し、かつ試験体の真の形状
を示す信号をフーリエ変換して、上記既知形状の試験体
の撮像により得られたフーリエ変換関数を上記真の形状
によるフーリエ変換関数で除し、上記測定装置の分解能
を示すフーリエ変換関数を求める第1の手段と、目的と
する試験体を上記測定装置で撮像して得られた形状を示
す信号をフーリエ変換し、上記目的とする試験体の撮像
により得られたフーリエ変換関数を上記測定装置の分解
能を示すフーリエ変換関数で除し、得られたフーリエ変
換関数を逆フーリエ変換したスペクトルを求める第2の
手段と、第2の手段により得られたスペクトルのパワー
スペクトルを、スペクトル幅を中心として座標を折り返
したスペクトルに変換して、画像として表示する第3の
手段とを備えて構成されていることを特徴とする撮像装
置。 2、前記第3の手段は、前記パワースペクトルを平滑化
し、平滑前のパワースペクトルから平滑後のパワースペ
クトルを減じてリップルスペクトルを得、上記リップル
スペクトルをフーリエ変換して得たスペクトルから、ス
ペクトル幅の中心から所定範囲内のスペクトルを抽出し
、前記第2の手段から出力されたスペクトルから上記抽
出されたスペクトルを減じる手段を備えていることを特
徴とする特許請求の範囲第1項記載の撮像装置。 3、前記第2の手段は、前記目的とする試験体の撮像に
より得られたフーリエ変換関数に、所定周期の余弦関数
と正弦関数を乗じ、これによつて得られたスペクトルを
前記測定装置の分解能を示すフーリエ変換関数で除する
手段を備えていることを特徴とする特許請求の範囲第1
項記載の撮像装置。
[Claims] 1. Fourier transform a signal representing the shape obtained by imaging a test specimen with a known shape using a measuring device, and Fourier transform a signal representing the true shape of the test specimen to obtain the known shape. A first means for obtaining a Fourier transform function indicative of the resolution of the measuring device by dividing a Fourier transform function obtained by imaging the shape of the test object by the Fourier transform function according to the true shape; A signal representing the shape obtained by imaging with the above measuring device is Fourier transformed, and the Fourier transform function obtained by imaging the above-mentioned target test object is divided by the Fourier transform function indicating the resolution of the above measuring device. A second means for obtaining a spectrum obtained by inverse Fourier transform of the Fourier transform function obtained, and converting the power spectrum of the spectrum obtained by the second means into a spectrum whose coordinates are folded around the spectral width as an image. An imaging device comprising: third means for displaying. 2. The third means smoothes the power spectrum, subtracts the power spectrum after smoothing from the power spectrum before smoothing to obtain a ripple spectrum, and calculates the spectral width from the spectrum obtained by Fourier transforming the ripple spectrum. The imaging system according to claim 1, further comprising means for extracting a spectrum within a predetermined range from the center of the image pickup apparatus, and subtracting the extracted spectrum from the spectrum output from the second means. Device. 3. The second means multiplies the Fourier transform function obtained by imaging the target test object by a cosine function and a sine function of a predetermined period, and transmits the resulting spectrum to the measuring device. Claim 1, further comprising means for dividing by a Fourier transform function indicating resolution.
The imaging device described in Section 1.
JP60155136A 1985-07-16 1985-07-16 Image display processor Expired - Lifetime JPH0664019B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP60155136A JPH0664019B2 (en) 1985-07-16 1985-07-16 Image display processor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP60155136A JPH0664019B2 (en) 1985-07-16 1985-07-16 Image display processor

Publications (2)

Publication Number Publication Date
JPS6217654A true JPS6217654A (en) 1987-01-26
JPH0664019B2 JPH0664019B2 (en) 1994-08-22

Family

ID=15599338

Family Applications (1)

Application Number Title Priority Date Filing Date
JP60155136A Expired - Lifetime JPH0664019B2 (en) 1985-07-16 1985-07-16 Image display processor

Country Status (1)

Country Link
JP (1) JPH0664019B2 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6486974A (en) * 1987-06-19 1989-03-31 Univ Harvard Drug compound and implant for curing meniscus tissue
WO1992003016A2 (en) * 1990-08-03 1992-02-20 Thomson-Csf Charge transfer type panning camera in a security system
JP2006275664A (en) * 2005-03-28 2006-10-12 Toshiba It & Control Systems Corp Radiographical inspection device and radiographical inspection method
JP2010513867A (en) * 2006-12-15 2010-04-30 オトクリトエ アクツィオネルノエ オブシェストボ ”ノルディンクラフト” Method for enhancing interference protection in ultrasonic testing and apparatus for carrying out this method
JP2013210268A (en) * 2012-03-30 2013-10-10 Nippon Telegr & Teleph Corp <Ntt> Two-dimensional image reconstruction method

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6486974A (en) * 1987-06-19 1989-03-31 Univ Harvard Drug compound and implant for curing meniscus tissue
JPH0311785B2 (en) * 1987-06-19 1991-02-18 Purejidento Ando Fueroozu Obu Haabaado Karetsuji
WO1992003016A2 (en) * 1990-08-03 1992-02-20 Thomson-Csf Charge transfer type panning camera in a security system
JP2006275664A (en) * 2005-03-28 2006-10-12 Toshiba It & Control Systems Corp Radiographical inspection device and radiographical inspection method
JP2010513867A (en) * 2006-12-15 2010-04-30 オトクリトエ アクツィオネルノエ オブシェストボ ”ノルディンクラフト” Method for enhancing interference protection in ultrasonic testing and apparatus for carrying out this method
JP2013210268A (en) * 2012-03-30 2013-10-10 Nippon Telegr & Teleph Corp <Ntt> Two-dimensional image reconstruction method

Also Published As

Publication number Publication date
JPH0664019B2 (en) 1994-08-22

Similar Documents

Publication Publication Date Title
JPH0233190B2 (en)
CN109087348B (en) Single-pixel imaging method based on adaptive area projection
CN101122580A (en) Inspection systems and methods of operation
Beyerer et al. Suppression of inhomogeneities in images of textured surfaces
JPS6217654A (en) Image pick-up apparatus
JP2528829B2 (en) Imaging method and imaging apparatus
US8520216B2 (en) Shape measuring apparatus
JPS6140678A (en) Image processor
JP2001174235A (en) Pattern analyzing method utilizing fourier transformation
CN116109962A (en) Vibration frequency measurement method and device based on video phase, electronic equipment and storage medium
JPH03262947A (en) Device and method for computer tomography and computer tomographic scanner
JPS63198862A (en) Eddy current flaw detecting device capable of inspecting thickness reduction
CN110764088B (en) Super-resolution standing point scanning real-time imaging algorithm
CN114488139A (en) Imaging method, imaging device, equipment and security check system for detection object
Klestov et al. Digital 3d x-ray microtomographic scanners for electronic equipment testing
JP7147894B2 (en) signal source estimator
JPH08212341A (en) Method and device for image processing of radiograph
EP0648359B1 (en) Correlation signal processing
JPS61751A (en) Ultrasonic medium characteristic measuring apparatus
JPS63278183A (en) Picture distortion correction device
Wu et al. EMI Detection on Time Domain–Practice through Cloud Computing
Popescu et al. Phantom-based point spread function estimation for terahertz imaging system
JP5085144B2 (en) Ultrasonic diagnostic equipment
Drobakhin et al. Multifrequency microwave measurements in the intermediate zone of radiation with inverse aperture synthesis for radio image forming
JP3492875B2 (en) Radar signal processing equipment