JP4521584B2 - Sonic velocity measuring method using ultrasonic microscope, sonic velocity measuring device, sonic velocity image acquiring method using ultrasonic microscope, and diagnostic imaging device thereof - Google Patents

Sonic velocity measuring method using ultrasonic microscope, sonic velocity measuring device, sonic velocity image acquiring method using ultrasonic microscope, and diagnostic imaging device thereof Download PDF

Info

Publication number
JP4521584B2
JP4521584B2 JP2004105213A JP2004105213A JP4521584B2 JP 4521584 B2 JP4521584 B2 JP 4521584B2 JP 2004105213 A JP2004105213 A JP 2004105213A JP 2004105213 A JP2004105213 A JP 2004105213A JP 4521584 B2 JP4521584 B2 JP 4521584B2
Authority
JP
Japan
Prior art keywords
measured
pulse
reflected wave
ultrasonic
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2004105213A
Other languages
Japanese (ja)
Other versions
JP2005291827A (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.)
Honda Electronics Co Ltd
Original Assignee
Honda Electronics Co 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 Honda Electronics Co Ltd filed Critical Honda Electronics Co Ltd
Priority to JP2004105213A priority Critical patent/JP4521584B2/en
Publication of JP2005291827A publication Critical patent/JP2005291827A/en
Application granted granted Critical
Publication of JP4521584B2 publication Critical patent/JP4521584B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04Analysing solids
    • G01N29/06Visualisation of the interior, e.g. acoustic microscopy
    • G01N29/0609Display arrangements, e.g. colour displays
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/028Material parameters
    • G01N2291/02854Length, thickness
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/028Material parameters
    • G01N2291/0289Internal structure, e.g. defects, grain size, texture

Landscapes

  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

この発明は、超音波顕微鏡を使用した音速測定方法、その音速測定装置、超音波顕微鏡を使用した音速像取得方法およびその画像診断装置に関する。詳しくは超音波顕微鏡として特に、広帯域のパルス励起型超音波顕微鏡を使用して生体組織などの被測定物の音速を測定するに当たり、超音波振動子が正対する試料載置板の機械的な傾き量を信号処理によって補正することで、試料載置板の傾き量に対する機械的な事前調整を不要にして、測定時間の大幅短縮を図ると共に、機械的な調整誤差による測定精度への影響を除去したものである。
The present invention relates to a sound velocity measuring method using an ultrasonic microscope, a sound velocity measuring device thereof, a sound velocity image acquiring method using an ultrasonic microscope, and an image diagnostic device thereof. Specifically, when measuring the speed of sound of an object to be measured such as a living tissue using a wide-band pulse excitation type ultrasonic microscope, the mechanical tilt of the sample mounting plate directly facing the ultrasonic transducer Correcting the amount by signal processing eliminates the need for mechanical pre-adjustment with respect to the amount of inclination of the sample mounting plate, greatly shortening the measurement time and eliminating the effect of measurement errors due to mechanical adjustment errors It is a thing.

さらに、生体組織などの性状を音速像として取得するに当たり、試料載置板の機械的な傾き量を電気的に補正することで、音速像の処理時間を短縮すると共に、高精度な音速像が得られるようにしたものである。   Furthermore, when acquiring the properties of biological tissue as a sonic image, the mechanical tilt amount of the sample mounting plate is electrically corrected, thereby reducing the processing time of the sonic image and producing a highly accurate sonic image. It is intended to be obtained.

超音波顕微鏡を使用して試料の表面状態を観察したり、その音速を測定する手法は以前から知られている(例えば、非特許文献1)。   A technique for observing the surface state of a sample using an ultrasonic microscope and measuring the speed of sound has been known (for example, Non-Patent Document 1).

また、一方では生体組織の観察及び診断に、この超音波顕微鏡などが応用されるようになってきた。例えば、悪性腫瘍を切除するような場合、切除部位の組織を摘出したり、切除部位近傍の生体組織を摘出して組織観察する生体組織診断が行われている。   On the other hand, this ultrasonic microscope has been applied to the observation and diagnosis of living tissue. For example, when a malignant tumor is excised, biological tissue diagnosis is performed in which the tissue at the excision site is removed or the biological tissue near the excision site is excised and observed.

この生体組織診断に使用される診断装置(観察装置)として従来から光学顕微鏡が使用されている。これは摘出した生体組織を染色した状態でその切除面を光学顕微鏡で観察しているが、完全な染色を行うには1日以上待機しなければならないため、光学顕微鏡を使用した組織観察は術後になる。   Conventionally, an optical microscope has been used as a diagnostic apparatus (observation apparatus) used for this biological tissue diagnosis. In this state, the excised surface is observed with an optical microscope in a state where the extracted living tissue is stained. However, in order to perform complete staining, it is necessary to wait for one day or more. Later.

また、染色する場合でも術中に観察が可能なように術中迅速診断法も開発されているが、この場合には、通常の染色による場合には組織輪郭の内部に「染色によって塗り分けられる」はずの組織情報が、術中迅速診断法では失われてしまうことが報告されている。   In addition, an intraoperative rapid diagnosis method has been developed so that observation can be performed during surgery even when staining is performed. In this case, in the case of normal staining, the tissue outline should be “colored by staining”. It has been reported that the tissue information is lost by intraoperative rapid diagnosis.

このようなことから、従来の染色方法を用いないでも、術中に摘出した生体組織の診断を行える術中診断法や、術中診断装置の開発が切望されている。この要望に応えるものとして、超音波顕微鏡を使用した診断装置の研究および開発がなされている。   For this reason, development of an intraoperative diagnostic method and an intraoperative diagnostic apparatus capable of diagnosing a living tissue extracted during surgery without using a conventional staining method is eagerly desired. In response to this demand, research and development of a diagnostic apparatus using an ultrasonic microscope has been conducted.

超音波顕微鏡とは、生体組織に超音波を照射し、その反射波の強度や位相などを演算処理して可視像化しながら観察することで、その生体組織を診断する手法である(非特許文献1参照)。   An ultrasonic microscope is a method of diagnosing a living tissue by irradiating the living tissue with ultrasonic waves and observing the reflected wave intensity and phase while performing a calculation process and visualizing it (non-patent document). Reference 1).

図7はこれを模式的に示した概念図であって、超音波顕微鏡1は送受波部10を有する。この送受波部10にはその先端部側に超音波振動子11が内蔵されている。送受波部10と対峙するようにガラス板などを使用したプレパラート(試料載置板)12が置かれ、このプレパラート12上に被測定物(例えば、生体組織切片)13が載置される。   FIG. 7 is a conceptual diagram schematically showing this, and the ultrasonic microscope 1 has a wave transmitting / receiving unit 10. An ultrasonic transducer 11 is built in the wave transmitting / receiving unit 10 on the tip side. A preparation (sample mounting plate) 12 using a glass plate or the like is placed so as to face the transmission / reception unit 10, and an object to be measured (for example, a biological tissue section) 13 is placed on the preparation 12.

そして、この被測定物13をターゲットとするように被測定物面に対して、超音波が水などの媒質17を介して照射される。超音波出射面10aは図示するように湾曲しており、超音波はビーム状に絞り込まれた状態で被測定物13の面に到達するようになされている。ビームの絞り込み方などで分解能が決まる。   Then, ultrasonic waves are applied to the surface of the object to be measured through the medium 17 such as water so as to target the object 13 to be measured. The ultrasonic wave emission surface 10a is curved as shown in the figure, and the ultrasonic wave reaches the surface of the object to be measured 13 in a state of being narrowed down into a beam shape. The resolution is determined by how the beam is narrowed down.

被測定物13の面上を二次元走査できるように、この例では二次元走査手段16が設けられ、この二次元走査手段16に上述したプレパラート12が載置される。図7の例では、二次元走査手段16としてX−Yステージを使用した場合を示す。X−Yステージ16はXステージ16XとYステージ16Yとで構成される。   In this example, a two-dimensional scanning unit 16 is provided so that the surface of the measurement target 13 can be two-dimensionally scanned, and the above-described preparation 12 is placed on the two-dimensional scanning unit 16. In the example of FIG. 7, a case where an XY stage is used as the two-dimensional scanning unit 16 is shown. The XY stage 16 includes an X stage 16X and a Y stage 16Y.

したがって、超音波出射側からの平面図は図8のようになり、また図9のようにx、x’方向(水平方向)への往復走査と、y方向への走査を行うことによって被測定物13に対して超音波を二次元的に走査することができる。   Therefore, the plan view from the ultrasonic wave emission side is as shown in FIG. 8, and as shown in FIG. 9, the measurement is performed by performing reciprocal scanning in the x and x ′ directions (horizontal direction) and scanning in the y direction. The object 13 can be scanned two-dimensionally with ultrasonic waves.

超音波顕微鏡1は音波の伝搬という定量化が容易な物理的パラメータを取り扱うのが特長である。具体的には反射波(反射波信号)の強度や位相をパラメータとして取り扱う。反射波を可視像として観察するため、図10のようにプレパラート12の面上に照射された超音波の反射波信号(参照反射波信号)Srと、被測定物13の面上に照射された超音波の反射波信号Soが利用される。   The acoustic microscope 1 is characterized by handling physical parameters that are easy to quantify such as propagation of sound waves. Specifically, the intensity and phase of the reflected wave (reflected wave signal) are handled as parameters. In order to observe the reflected wave as a visible image, the reflected ultrasonic wave signal (reference reflected wave signal) Sr irradiated on the surface of the preparation 12 and the surface of the object to be measured 13 are irradiated as shown in FIG. The reflected ultrasonic wave signal So is used.

ここに、反射波信号Soは、被測定物13の表面から反射した反射波信号Ssと、その裏面側つまりプレパラート12の面で反射した反射波信号Sdとの干渉波である。   Here, the reflected wave signal So is an interference wave between the reflected wave signal Ss reflected from the surface of the DUT 13 and the reflected wave signal Sd reflected from the back side thereof, that is, the surface of the preparation 12.

図11には超音波振動子11を従来のようにバースト信号で励起するのではなく、説明の都合上、広帯域のパルス信号で励起したときのそれぞれの波形が示されている。   FIG. 11 shows respective waveforms when the ultrasonic transducer 11 is not excited with a burst signal as in the prior art, but is excited with a broadband pulse signal for convenience of explanation.

図11Aは超音波振動子11を励起する励起パルスSiの波形である。ほぼ1サイクルのパルス波形で、その周波数特性は図12のようになる。この例では、中心周波数が80MHz程度で、40〜150MHz程度までの帯域を持つ励起パルスが使用される。図11Bはそのとき得られる参照反射波信号Srの波形である。また図11Cは反射波信号SsとSdとが干渉した干渉反射波信号Soの波形が示されている。   FIG. 11A shows a waveform of an excitation pulse Si that excites the ultrasonic transducer 11. FIG. 12 shows the frequency characteristics of a pulse waveform of almost one cycle. In this example, an excitation pulse having a center frequency of about 80 MHz and a band of about 40 to 150 MHz is used. FIG. 11B shows the waveform of the reference reflected wave signal Sr obtained at that time. FIG. 11C shows the waveform of the interference reflected wave signal So in which the reflected wave signals Ss and Sd interfere.

図示する干渉反射波信号Soは、被測定物13として生体組織切片を用いたときの波形で、励起パルスSi自身の波形(Si’)のほかに、反射波信号Suが得られる。超音波は超音波振動子11と被測定物13との間で繰り返し反射されるものであるから、最初に得られる反射波信号を一次反射波信号Su1とすると、数次に至った反射波信号Su2,Su3,・・・が得られる。   The interference reflected wave signal So shown in the figure is a waveform when a biological tissue section is used as the object 13 to be measured, and a reflected wave signal Su is obtained in addition to the waveform (Si ′) of the excitation pulse Si itself. Since the ultrasonic wave is repeatedly reflected between the ultrasonic transducer 11 and the object 13 to be measured, if the reflected wave signal obtained first is the primary reflected wave signal Su1, the reflected wave signal reaching several orders. Su2, Su3,... Are obtained.

ここで、このような構成をなす超音波顕微鏡にあっては、得られた反射波信号Sr,Soの強度(反射波強度)や位相を利用して試料である被測定物の音速が求められる。このときの測定精度を高めるためにはプレパラート12の傾きに影響されない測定方法や測定装置であることが望ましい。また、この超音波顕微鏡1を生体の画像診断として利用する場合には、被測定物13は摘出した生体組織切片となる。   Here, in the ultrasonic microscope having such a configuration, the sound velocity of the object to be measured, which is a sample, is obtained using the intensity (reflected wave intensity) and phase of the obtained reflected wave signals Sr and So. . In order to increase the measurement accuracy at this time, it is desirable that the measurement method or the measurement apparatus be not affected by the inclination of the preparation 12. When the ultrasonic microscope 1 is used for image diagnosis of a living body, the object to be measured 13 is an extracted living tissue section.

以下に画像診断について説明する。生体組織切片からの反射波の強度や位相に基づいて生体組織を評価することで、当該生体組織内における音響減衰評価が可能になる。   Image diagnosis will be described below. By evaluating the living tissue based on the intensity and phase of the reflected wave from the living tissue section, it is possible to evaluate the acoustic attenuation in the living tissue.

非特許文献2にも示されるように、バースト信号を利用した超音波顕微鏡による診断ではこのように音響減衰による評価が主で、したがって、従来では生体組織の音速評価は全くなされていない。   As shown in Non-Patent Document 2, the diagnosis based on the acoustic microscope using the burst signal is mainly based on the evaluation based on the acoustic attenuation as described above. Therefore, the sound velocity of the living tissue has not been evaluated at all.

生体組織内に超音波を照射した場合、超音波は生体組織の性状によって得られる音速が相違することが知られている。   It is known that when ultrasonic waves are irradiated into a living tissue, the ultrasonic waves have different sound speeds depending on the properties of the living tissue.

つまり、生体組織と音速とは相関性が高く、生体組織が正常な組織である場合と、病変部のように変質した組織の場合とでは、その物理的性質の違い、換言すれば生体組織の弾性の違いが音速の違いとなって検知できる。   In other words, there is a high correlation between living tissue and sound velocity, and the difference in physical properties between the case where the living tissue is a normal tissue and the case where the tissue is altered like a lesion, in other words, the living tissue. Differences in elasticity can be detected as differences in sound speed.

超音波を使用する場合では、音響特性の相違から生体組織の違いとして「塗り分ける」ことが可能であることが、最近の諸種の実験によって確認されている。したがって、超音波による生体組織の性状診断を行うには、生体組織の音響特性をミクロンレベルで計測することが必要になる。   In the case of using ultrasonic waves, it has been confirmed by various kinds of recent experiments that it can be “painted” as a difference in biological tissue from a difference in acoustic characteristics. Therefore, in order to perform a property diagnosis of a living tissue using ultrasonic waves, it is necessary to measure the acoustic characteristics of the living tissue at a micron level.

その一方で、術中診断を可能にするためには、
(1)摘出生体組織切片の迅速な可視像生成が可能であること
(2)できれば、超音波顕微鏡を使用した診断装置の手術室内への搬入が可能であること
(3)診断のための装置操作が簡易であること
(4)目視できること
などの諸種の課題をクリアしなければならない。
On the other hand, to enable intraoperative diagnosis,
(1) It is possible to quickly generate a visible image of an excised living tissue section (2) If possible, it is possible to carry a diagnostic device using an ultrasonic microscope into an operating room (3) For diagnosis It is necessary to clear various problems such as simple operation of the apparatus (4) being visible.

さらには、商業ベースを考えると、比較的安価に提供できることも必要である。あまり高価であると、装置設置者にとって負担が大きく、その負担の一部は患者への医療費負担に還元されてしまうからである。   Furthermore, considering a commercial base, it is also necessary to be able to provide it at a relatively low cost. This is because if the cost is too high, the burden on the apparatus installer is large, and a part of the burden is reduced to the medical cost burden on the patient.

従来の超音波顕微鏡を使用した診断装置は、上述したように単一周波数のバースト信号を超音波信号として利用し、反射した超音波信号の強度や位相を解析することで、生体組織の性状を観察している。   A diagnostic device using a conventional ultrasonic microscope uses a burst signal of a single frequency as an ultrasonic signal as described above, and analyzes the intensity and phase of the reflected ultrasonic signal, so that the properties of biological tissue can be determined. Observe.

プレパラート12から直接反射された参照反射波信号Srと生体組織切片13からの反射波信号Soを用いて生体組織切片13の厚みdや生体組織の音速を算出することができる。その詳細は後述する。   Using the reference reflected wave signal Sr directly reflected from the preparation 12 and the reflected wave signal So from the biological tissue section 13, the thickness d of the biological tissue section 13 and the sound speed of the biological tissue can be calculated. Details thereof will be described later.

ところで、上述した表面反射波信号Ssと裏面反射波信号Sdとの時間差よりも、超音波振動子11に加えるバースト信号の持続時間を短くすることが難しいため、通常は表面反射波信号Ssと裏面反射波信号Sdとが互いに干渉した干渉反射波信号Soを用いて解析せざるを得ない。   By the way, since it is difficult to shorten the duration of the burst signal applied to the ultrasonic transducer 11 rather than the time difference between the surface reflected wave signal Ss and the back surface reflected wave signal Sd described above, the surface reflected wave signal Ss and the back surface are usually used. It must be analyzed using the interference reflected wave signal So in which the reflected wave signal Sd interferes with each other.

干渉による影響をなくすため、干渉波を分離できるようにバースト信号の周波数を順次変化させて、そのときの反射波の強度と位相スペクトルを解析することで、生体組織の厚さと音速を算出している。   In order to eliminate the influence of interference, the burst signal frequency is sequentially changed so that the interference wave can be separated, and the intensity and phase spectrum of the reflected wave at that time are analyzed to calculate the thickness and sound speed of the living tissue. Yes.

研究室レベルで研究開発が進行しているこの超音波顕微鏡を使用した診断装置では、上述したようにバースト信号を使用した位相検波方式を採用している。上述したように得られた反射波の強度と位相スペクトルを解析して、生体組織の厚みや音速を算出し、得られた音響特性から生体組織の性状を診断するものであるから、可視像化するまでの処理時間が非常に長い。   In the diagnostic apparatus using the ultrasonic microscope, which is under research and development at the laboratory level, the phase detection method using the burst signal is adopted as described above. Analyzing the intensity and phase spectrum of the reflected wave obtained as described above, calculating the thickness and speed of sound of the living tissue, and diagnosing the properties of the living tissue from the obtained acoustic characteristics, so a visible image It takes a very long time to process.

例えば、2mm×2mmの生体組織切片を100×100ピクセルとして測定した場合でも10分以上費やしてしまう。したがってこれを300×300ピクセル程度まで測定精度を上げるためには90分程度の処理時間がかかることになり、病理判断を行う術中診断装置としては充分とは言えない。   For example, even when a 2 mm × 2 mm biological tissue section is measured as 100 × 100 pixels, it takes 10 minutes or more. Therefore, in order to increase the measurement accuracy up to about 300 × 300 pixels, it takes about 90 minutes of processing time, which is not sufficient as an intraoperative diagnostic apparatus for performing pathological judgment.

計測時間が長くかかるのは、さらに生体組織切片に対する走査速度が関与している。従来装置では生体組織切片をX−Yステージ16の上に載せ、これをx、y方向に間欠的に動かすことで、例えば100×100ピクセルの生体像を得ている。この間欠走査による生体像の取得が測定時間を遅延させる大きな要因となっている。さらに生体組織の切片から得られた反射波を処理して音速情報を得る信号処理速度の遅さも一因となっている。   The long measurement time further involves the scanning speed of the biological tissue section. In the conventional apparatus, a biological tissue section is placed on the XY stage 16 and is moved intermittently in the x and y directions to obtain a biological image of, for example, 100 × 100 pixels. Acquisition of a biological image by this intermittent scanning is a major factor that delays the measurement time. In addition, the slow signal processing speed for obtaining sound speed information by processing the reflected wave obtained from a slice of biological tissue is also a factor.

このような問題を解決するための超音波顕微鏡を使用した画像診断装置として、近年パルス励起による超音波顕微鏡を使用した画像診断装置の研究が進められている。   As an image diagnostic apparatus using an ultrasonic microscope for solving such a problem, research on an image diagnostic apparatus using an ultrasonic microscope based on pulse excitation has recently been advanced.

上述したようにほぼ1サイクル分の広帯域パルスを用いて超音波振動子を励起する方式である。このパルス励起型超音波顕微鏡を用いると共に、X−Yステージ16のx、y方向における走査速度を速め、さらに信号処理系を工夫することによって、超音波を照射してから最終的な診断画像を可視像として描画するまでの時間を、従来のバースト信号を使用した超音波顕微鏡よりも格段に短縮できることが確認された。   As described above, this is a method of exciting an ultrasonic transducer using a broadband pulse for approximately one cycle. By using this pulse excitation type ultrasonic microscope, the scanning speed in the x and y directions of the XY stage 16 is increased, and the signal processing system is further devised to irradiate the ultrasonic wave before final diagnostic images are obtained. It was confirmed that the time required for drawing as a visible image can be significantly reduced as compared with a conventional ultrasonic microscope using a burst signal.

本出願人らによる実験によれば、図13のように、参照面としてのプレパラート12を含めた300×300ピクセルの画像を描画する場合でも2分以下に短縮できることが判った。もちろん、このような画像処理時間の短縮は、同時に超音波を収束させて照射する超音波ビームを使用したこと、超音波の励起周波数(ほぼ中心周波数)を80MHz程度に選定することなどが大きく寄与していることは明らかである。   According to the experiments by the present applicants, it was found that even when an image of 300 × 300 pixels including the preparation 12 as a reference surface is drawn as shown in FIG. Of course, such a reduction in image processing time greatly contributes to the use of an ultrasonic beam that simultaneously converges and irradiates the ultrasonic wave, and the selection of the ultrasonic excitation frequency (approximately the center frequency) to about 80 MHz. Obviously.

図14に超音波顕微鏡を使用して作成した生体組織の画像(生体像)の一例を示す。プレパラート12の面が参照画像20となり、生体組織切片が生体像21となって得られる。生体像は減衰像や音速像が考えられる。   FIG. 14 shows an example of a biological tissue image (biological image) created using an ultrasonic microscope. The surface of the preparation 12 becomes the reference image 20 and the biological tissue section becomes the biological image 21. The biological image can be an attenuated image or a sound velocity image.

なお、パルス励振型の超音波顕微鏡や超音波治療装置としては特許文献1〜特許文献3などが知られている。   Patent Documents 1 to 3 are known as pulse excitation type ultrasonic microscopes and ultrasonic treatment apparatuses.

特許文献1は、パルス増幅器からのノイズを遮断するためのスイッチ回路を備えた超音波顕微鏡が開示されている。したがってこの発明とは直接関係のない技術である。   Patent Document 1 discloses an ultrasonic microscope provided with a switch circuit for blocking noise from a pulse amplifier. Therefore, the present invention is not directly related to the present invention.

特許文献2は、パルス励振型超音波治療装置に関する技術であって、特にパルスの焦点位置がずれないような治療装置に関するもので、これまたこの発明の内容とは直接関係がない。   Patent Document 2 is a technique related to a pulse excitation type ultrasonic therapeutic apparatus, and particularly relates to a therapeutic apparatus in which the focal position of a pulse is not shifted, and is not directly related to the contents of the present invention.

また、特許文献3は、三次元超音波顕微鏡に関する技術が開示されている。これは球面超音波によって三次元に対象物を観察できるようにしたもので、二次元観察を対象とするこの発明とはその技術が相違する。   Patent Document 3 discloses a technique related to a three-dimensional ultrasonic microscope. This is a technique in which an object can be observed three-dimensionally by spherical ultrasonic waves, and the technique is different from that of the present invention which is intended for two-dimensional observation.

「超音波技術とその応用:超音波顕微鏡の最新成果」 (「電子材料」1992年11月号(102〜107頁) 工業調査会発行)"Ultrasonic technology and its application: Latest results of ultrasonic microscope" ("Electronic materials" November 1992 issue (pp. 102-107), published by Industrial Research Council) 「医用超音波:パルス励振型超音波顕微鏡」 (「超音波TECHNO」 VOL.15 No.6(2003.11〜12)(101〜105頁)日本工業出版社発行)“Medical Ultrasound: Pulse-excitation Ultrasound Microscope” (“Ultrasonic TECHNO” VOL.15 No.6 (November 11-12, 2003) (101-105 pages), published by Nihon Kogyo Shuppansha) 特開2001−46370号公報JP 2001-46370 A 特開平7−303657号公報Japanese Patent Laid-Open No. 7-303657 特表2000−517414号公報Special Table 2000-517414

ところで、このようなパルス励起による超音波顕微鏡を使用した高速処理可能な音速測定にあっては、超音波を照射したときの反射波から、被測定物の音速を正確に測定しなければならないし、画像診断装置にあっては、照射点(測定点)での生体組織の音速を正確に算出しなければならない。そのためには、被測定物と超音波振動子(トランスデューサ)11が正対し、図10のように参照基板であるプレパラート12や被測定物13に対してそれぞれ直角に超音波が照射されることが必要である。   By the way, in such a sound speed measurement that can be processed at high speed using an ultrasonic microscope by pulse excitation, the sound speed of the object to be measured must be accurately measured from the reflected wave when the ultrasonic wave is irradiated. In the diagnostic imaging apparatus, the speed of sound of the living tissue at the irradiation point (measurement point) must be accurately calculated. For this purpose, the object to be measured and the ultrasonic transducer (transducer) 11 face each other, and ultrasonic waves are radiated at right angles to the preparation 12 and the object to be measured 13 which are reference substrates as shown in FIG. is necessary.

被測定物13はプレパラート12の上に載せられており、このプレパラート12はさらにX−Yステージ16に載っているものであるから、超音波振動子11とプレパラート12および被測定物13がそれぞれ正対の関係にあるとは限らず、微少ではあるが、僅かな傾きを以て正対する場合が殆どである。この傾きがゼロでないと、参照反射波信号Srおよび干渉反射波信号Soから厚みや音速を正確に算出することができず、僅かな計測誤差が発生してしまう。   Since the object to be measured 13 is placed on the preparation 12, and this preparation 12 is further placed on the XY stage 16, the ultrasonic transducer 11, the preparation 12, and the object 13 to be measured are respectively correct. The pair is not necessarily in a pair relationship, and although it is very small, it is almost the case where it faces directly with a slight inclination. If the inclination is not zero, the thickness and sound speed cannot be accurately calculated from the reference reflected wave signal Sr and the interference reflected wave signal So, and a slight measurement error occurs.

例えば、図15に示すように基準面refに対してプレパラート12の平面12aが傾きθをもっているときには、参照点Poと被測定物13が載置されているところの測定点Ppとの関係では、測定点Ppは、z軸方向(超音波出射方向)に対して参照点PoよりもΔdだけずれている。   For example, as shown in FIG. 15, when the plane 12 a of the preparation 12 has an inclination θ with respect to the reference plane ref, the relationship between the reference point Po and the measurement point Pp where the DUT 13 is placed is The measurement point Pp is shifted by Δd from the reference point Po with respect to the z-axis direction (ultrasound emission direction).

このz軸方向での測定点Ppのずれによって、測定された厚みが変わる。また厚みが変わると音速の算出誤差が生じ、この誤差による影響で、より精密で、きめ細かな画像診断の妨げとなるおそれがある。   The measured thickness changes due to the deviation of the measurement point Pp in the z-axis direction. In addition, if the thickness changes, a sound speed calculation error occurs, and the influence of this error may hinder more precise and detailed image diagnosis.

特に、被測定物として生体組織を観察するときは、この生体組織を極めて薄く切り落として使用する。通常4〜10μm程度の厚みとなるようにスライスする。これに対してプレパラート12の傾きは微調整後であっても1μm程度の調整誤差(傾き誤差)が発生しているのが常である。しかし、この調整誤差は生体組織の厚みからすると看過できない値であるから、極めて高精度に調整して、傾きを限りなくゼロに近づけなければならない。   In particular, when observing a living tissue as an object to be measured, the living tissue is cut out extremely thinly and used. Usually, it slices so that it may become a thickness of about 4-10 micrometers. On the other hand, the inclination of the preparation 12 usually has an adjustment error (tilt error) of about 1 μm even after fine adjustment. However, since this adjustment error is a value that cannot be overlooked from the thickness of the living tissue, it must be adjusted with extremely high accuracy and the inclination must be as close to zero as possible.

従来の測定では、このずれΔdを機械的に補正している。例えば図16に示すように生体組織切片13が載置されていないプレパラート12の平面12aの任意の点Poを基準点とし、これとx軸方向での任意の点P1およびy軸方向での任意の点P2をとり、基準点Poからみたプレパラート12のx軸方向の傾き(Δθx)がゼロとなるようにプレパラート12より、具体的にはXステージ16Xのx軸方向の傾きが微調整される。同様にしてy軸方向の傾き(Δθy)がゼロとなるようにYステージ16Yの傾きが微調整される。   In the conventional measurement, this deviation Δd is mechanically corrected. For example, as shown in FIG. 16, an arbitrary point Po on the plane 12a of the preparation 12 on which the biological tissue section 13 is not placed is used as a reference point, and an arbitrary point P1 in the x-axis direction and an arbitrary point in the y-axis direction. The tilt of the X-axis direction of the preparation 12 as viewed from the reference point Po (Δθx) is finely adjusted from the preparation 12, more specifically, the X-stage tilt of the X stage 16X. . Similarly, the tilt of the Y stage 16Y is finely adjusted so that the tilt (Δθy) in the y-axis direction becomes zero.

このような機械的な微調整を繰り返し行うことで、超音波の出射方向とプレパラート12および生体組織切片13のそれぞれが直交するようになるので、音速の算出精度が高まる。   By repeatedly performing such mechanical fine adjustment, the ultrasonic wave emission direction and each of the preparation 12 and the biological tissue section 13 are orthogonal to each other, so that the calculation accuracy of the sound speed is increased.

しかし、このX−Yステージ16の超音波出射方向に対する傾きθは上述したように極めて僅かな値であるから、これを正確に微調整するためには相当の熟練を要するため、被測定物に対する事前調整が面倒で、事前の調整時間に相当時間を費やしてしまう。結果として迅速で高精度な生体組織診断を阻害する要因となっている。しかも、この機械的な微調整は、熟練を要すると共に、生体組織を診断する都度実施しなければならないので、非常に面倒であるなどの問題も惹起している。   However, since the inclination θ of the XY stage 16 with respect to the ultrasonic wave emission direction is extremely small as described above, considerable skill is required to accurately fine-tune this, so Pre-adjustment is cumbersome and spends considerable time on pre-adjustment time. As a result, it is a factor that impedes rapid and highly accurate biological tissue diagnosis. In addition, this mechanical fine adjustment requires skill and must be carried out every time a living tissue is diagnosed, which causes problems such as being very troublesome.

このようなことから、超音波顕微鏡を用いて音速を測定する場合にあっては、超音波振動子に対する被測定物の傾き補正を短時間にかつ高精度に行うことが極めて重要になる。そのことが超音波顕微鏡を用いた画像診断装置に対する画像診断の精度に大きく影響を及ぼすことになる。   For this reason, when measuring the speed of sound using an ultrasonic microscope, it is extremely important to correct the inclination of the object to be measured with respect to the ultrasonic transducer in a short time and with high accuracy. This greatly affects the accuracy of image diagnosis for an image diagnosis apparatus using an ultrasonic microscope.

そこで、この発明はこのような従来の課題を解決したものであって、特に機械的調整によることなく、換言すれば電気的な処理のみで傾き補正が可能な超音波顕微鏡を使用した音速測定方法、音速測定装置、及びこの傾き補正を応用した画像診断方法およびその画像診断装置を提案するものである。   Therefore, the present invention solves such a conventional problem, and is not particularly based on mechanical adjustment, in other words, a sound speed measuring method using an ultrasonic microscope capable of correcting an inclination only by electrical processing. The present invention proposes a sound speed measuring device, an image diagnostic method applying the tilt correction, and an image diagnostic device thereof.

上述の課題を解決するため、請求項1に記載したこの発明に係る超音波顕微鏡を使用した音速測定方法は、超音波振動子をパルス励起することで被測定物の音速を測定するパルス励起型超音波顕微鏡を使用した音速測定方法であって、所定の間隔を有するパルス励起用のトリガー信号を生成するトリガー信号生成ステップと、トリガー信号が供給されて、トリガー信号よりも間隔の狭い励起パルスで超音波振動子を集中的に励起する集中パルス励起ステップと、この集中パルス励起期間中に得られる励起パルスからの反射波信号を平均化するステップと、超音波振動子と対峙する試料載置板の傾き量を検出するステップと、検出されたこの傾き量と、平均化した反射波信号とに基づいて被測定物の厚みを算出するステップと、傾き量によって電気的に補正された被測定物の厚みに基づいて被測定物の音速を算出するステップとを有することを特徴とする。
In order to solve the above-mentioned problem, a sound speed measuring method using an ultrasonic microscope according to the present invention described in claim 1 is a pulse excitation type that measures the sound speed of an object to be measured by pulse exciting an ultrasonic transducer. the speed-of-sound measurement method using the ultrasonic microscope, a trigger signal generating step of generating a trigger signal for pulse excitation with a predetermined interval, and a trigger signal is supplied, the interval than the trigger signal narrow I励 electromotive a centralized pulsed excitation step of intensively excited ultrasonic transducer in pulse, a step of averaging the reflected wave signal from the obtained that excitation pulses during the intensive pulsed excitation period, facing the ultrasonic vibrator The step of detecting the amount of inclination of the sample mounting plate, the step of calculating the thickness of the object to be measured based on the detected amount of inclination and the averaged reflected wave signal, and the amount of inclination Characterized by a step of calculating the sound speed of the object to be measured based on the thickness of the gas-to corrected DUT.

請求項6に記載したこの発明に係る超音波顕微鏡を使用した音速測定装置は、超音波振動子をパルス励起することで被測定物の音速を測定するパルス励起型超音波顕微鏡を使用した音速測定装置あって、所定の間隔を有するパルス励起用のトリガー信号を生成するトリガー信号生成部と、トリガー信号生成部で生成されたトリガー信号が供給されて、トリガー信号よりも間隔の狭い励起パルスで超音波振動子を集中的に励起する集中パルス励起機能を有したパルサー部と、集中パルス励起期間中に得られるパルサー部で励起された励起パルスからの反射波信号を平均化する信号演算手段と、超音波振動子と対峙する試料載置板の傾き量を検出する傾き量検出手段と、検出されたこの傾き量と、信号演算手段で平均化された反射波信号とに基づいて被測定物の厚みを算出する厚み算出手段と、傾き量によって電気的に補正された被測定物の厚みに基づいて、被測定物の音速を測定する音速測定手段とで構成されたことを特徴とする。
The sound speed measuring apparatus using the ultrasonic microscope according to the present invention described in claim 6 is a sound speed measurement using a pulse excitation type ultrasonic microscope that measures the sound speed of the object to be measured by pulse exciting the ultrasonic vibrator. there device, a trigger signal generator for generating a trigger signal for pulse excitation with a predetermined interval, is supplied with the trigger signal generated by trigger signal generator, the interval than the trigger signal narrow I励 electromotive pulse in signal calculation for averaging a pulser unit having a centralized pulsed excitation function to focus exciting the ultrasonic transducer, a reflected wave signal from the excited excitation pulse pulser portion obtained during intensive pulsed excitation period Means, an inclination amount detecting means for detecting the inclination amount of the sample mounting plate facing the ultrasonic transducer, the detected inclination amount, and the reflected wave signal averaged by the signal calculating means. A thickness calculating means for calculating the thickness of the object to be measured, and a sound speed measuring means for measuring the sound speed of the object to be measured based on the thickness of the object to be measured electrically corrected by the amount of inclination. Features.

また、請求項11に記載したこの発明に係る超音波顕微鏡を使用した音速像取得方法は、超音波振動子をパルス励起することで被測定物の性状を音速像として取得するパルス励起型超音波顕微鏡を使用した音速像取得方法であって、所定の間隔を有するパルス励起用のトリガー信号を生成するトリガー信号生成ステップと、トリガー信号が供給されて、トリガー信号よりも間隔の狭い励起パルスで超音波振動子を集中的に励起する集中パルス励起ステップと、この集中パルス励起期間中に得られる励起パルスからの反射波信号を平均化するステップと、超音波振動子と対峙する試料載置板の傾き量を検出するステップと、検出されたこの傾き量と、平均化した反射波信号とに基づいて被測定物の厚みを算出するステップと、反射波信号に基づいて算出された被測定物の厚みを傾き量によって電気的に補正するステップと、被測定物の厚み情報から当該被測定物の音速を算出するステップと、この音速から被測定物の音速像を生成するステップとを有することを特徴とする。
According to another aspect of the present invention, there is provided a sonic image acquisition method using an ultrasonic microscope according to the present invention. the speed-of-sound image acquisition method using a microscope, and a trigger signal generating step of generating a trigger signal for pulse excitation with a predetermined interval, and a trigger signal is supplied, the interval than the trigger signal narrow I励 electromotive pulse in a centralized pulsed excitation step of intensively excited ultrasonic transducer, comprising the steps of averaging the reflected wave signal from the obtained that excitation pulses during the intensive pulsed excitation period, the sample facing the ultrasonic vibrator A step of detecting the amount of inclination of the mounting plate, a step of calculating the thickness of the object to be measured based on the detected amount of inclination and the averaged reflected wave signal, and based on the reflected wave signal The step of electrically correcting the thickness of the measured object by the amount of inclination, the step of calculating the sound speed of the measured object from the thickness information of the measured object, and the sound speed image of the measured object from the sound speed And generating.

さらに、請求項17に記載したこの発明に係る超音波顕微鏡を使用した画像診断装置では、超音波振動子をパルス励起することで被測定物の性状を音速像として取得するパルス励起型超音波顕微鏡を使用した画像診断装置であって、所定の間隔を有するパルス励起用のトリガー信号を生成するトリガー信号生成部と、トリガー信号生成部で生成されたトリガー信号が供給されて、トリガー信号よりも間隔の狭い励起パルスで超音波振動子を集中的に励起する集中パルス励起機能を有したパルサー部と、集中パルス励起期間中に得られるパルサー部で励起された励起パルスからの反射波信号を平均化する信号演算手段と、超音波振動子と対峙する試料載置板の傾き量を検出する傾き量検出手段と、検出されたこの傾き量と、信号演算手段で平均化された反射波信号とに基づいて被測定物の厚みを算出する厚み算出手段と、反射波信号に基づいて算出された被測定物の厚みを傾き量によって電気的に補正する補正手段と、被測定物の厚み情報から当該被測定物の音速を算出する音速算出手段と、この音速から被測定物の音速像を生成する音速像生成手段とを有することを特徴とする。 Furthermore, in the diagnostic imaging apparatus using the ultrasonic microscope according to the present invention as set forth in claim 17, a pulse excitation type ultrasonic microscope which acquires the property of the object to be measured as a sound velocity image by pulse exciting the ultrasonic transducer. , A trigger signal generation unit that generates a trigger signal for pulse excitation having a predetermined interval, and a trigger signal generated by the trigger signal generation unit is supplied, and the interval is longer than the trigger signal. narrow I励 and pulser unit having a centralized pulsed excitation function to focus exciting the ultrasonic transducers in force pulse, the reflected wave signal from the excited excitation pulse pulser portion obtained during intensive pulsed excitation period The signal calculation means for averaging, the inclination amount detection means for detecting the inclination amount of the sample mounting plate facing the ultrasonic transducer, and the detected inclination amount are averaged by the signal calculation means. A thickness calculating means for calculating the thickness of the object to be measured based on the reflected wave signal, a correcting means for electrically correcting the thickness of the object to be measured calculated based on the reflected wave signal by the amount of inclination, It has a sound speed calculating means for calculating the sound speed of the object to be measured from the thickness information of the object to be measured, and a sound speed image generating means for generating a sound speed image of the object to be measured from the sound speed.

この発明では、被測定物の音速を測定する前に、試料載置板としてのプレパラートの平面方程式Zに用いられる係数a,b,cを算出する。平面方程式を算出するためには少なくともプレパラートの平面上の任意の3点をとる。3点以上多点測定ももちろん可能であるし、その分精度が向上するが、実施例では3点を例示する。   In this invention, before measuring the sound velocity of the object to be measured, the coefficients a, b, and c used in the plane equation Z of the preparation as the sample mounting plate are calculated. In order to calculate the plane equation, at least three arbitrary points on the plane of the preparation are taken. Of course, multipoint measurement of three or more points is possible, and the accuracy is improved accordingly, but three points are exemplified in the embodiment.

この3点は何れも試料(被測定物)が載置されていない面内の点が選ばれる。3点のうちの1点例えばPoを参照点とする。参照点Poでのz軸方向の値(傾き)Zoはゼロとする。残りの2点でのz軸方向の値Z1,Z2は、実際に試料に超音波を照射して得た反射波信号Srを用いて算出する。具体的には3つの測定点Po,P1,P2から得られた反射波信号Sro,Sr1,Sr2を周波数成分に分解し、その周波数成分比からz軸方向の値Z1,Z2を算出することで、最終的に係数a,b,cを求めることができる。   These three points are all points on the surface where the sample (object to be measured) is not placed. One of the three points, for example, Po is set as a reference point. The value (slope) Zo in the z-axis direction at the reference point Po is zero. The z-axis direction values Z1 and Z2 at the remaining two points are calculated using the reflected wave signal Sr obtained by actually irradiating the sample with ultrasonic waves. Specifically, the reflected wave signals Sro, Sr1, Sr2 obtained from the three measurement points Po, P1, P2 are decomposed into frequency components, and values Z1, Z2 in the z-axis direction are calculated from the frequency component ratio. Finally, the coefficients a, b, and c can be obtained.

係数a,b,cが判れば、試料が載置されたプレパラートの平面内の任意の点におけるz軸方向の値Zp、つまり傾き量を算出できる。このz軸方向の値Zpが試料の厚みdを算出するときの補正値として利用されることで、プレパラートが超音波出射方向に対して傾いていたとしても、この傾きを考慮した状態で測定点での試料の音速を正確に算出できる。つまり、電気的に傾きを補正した状態で試料である被測定物の音速を求めることができる。   If the coefficients a, b, and c are known, the value Zp in the z-axis direction at an arbitrary point in the plane of the preparation on which the sample is placed, that is, the amount of inclination can be calculated. This z-axis direction value Zp is used as a correction value when calculating the thickness d of the sample, so that even if the preparation is tilted with respect to the ultrasonic wave emission direction, the measurement point is taken into account while taking this tilt into account. The sound speed of the sample can be accurately calculated. That is, the sound speed of the object to be measured, which is a sample, can be obtained in a state where the inclination is electrically corrected.

したがって、被測定物を画像として取得して診断を行う場合にあっても、プレパラートに被測定物を載せた状態で、上述した係数a,b,cを求めてから実際の画像(音速像)を求める作業が行われる。生体組織の厚みを高精度に算出できれば、その測定点での音速も正確に算出できることになり、音速の値を測定点ごとに等音速線として描画すれば、生体組織の音速像が得られる。したがってこの音速像の精度も大幅に向上することになり、生体組織をより緻密に診断することが可能になる。   Accordingly, even when the object to be measured is acquired as an image for diagnosis, the actual image (sonic image) is obtained after the above-described coefficients a, b, and c are obtained with the object to be measured placed on the slide. Work to ask for is done. If the thickness of the living tissue can be calculated with high accuracy, the sound velocity at the measurement point can also be calculated accurately. If the value of the sound velocity is drawn as an isosonic line for each measurement point, a sound velocity image of the living tissue can be obtained. Therefore, the accuracy of the sound velocity image is also greatly improved, and it becomes possible to diagnose the living tissue more precisely.

この発明は、生体組織が載置される載置面の超音波の出射方向に対する傾きを信号処理によって補正できるようにしたものである。   According to the present invention, the inclination of the placement surface on which the living tissue is placed with respect to the emission direction of the ultrasonic wave can be corrected by signal processing.

これによれば、載置面の傾きを信号処理によって補正した状態で被測定物の厚みや音速を算出できるから、音速測定や画像診断のための前処理時間を従来よりも大幅に短縮できる。さらに、載置面の傾きを電気的に補正することで被測定物の音速を算出するものであるから、機械的調整のように調整ミスによる誤差がなくなり、その結果、音速測定や画像診断の都度、得られたデータにばらつきが殆どなくなり、算出された音速の精度が格段に向上する。したがって、被測定物の音速測定を高精度に、かつ極めて短時間に行うことができると共に、被測定物をより正確に反映した音速像を取得できる実益を有する。   According to this, since the thickness and sound speed of the object to be measured can be calculated in a state in which the inclination of the placement surface is corrected by signal processing, the preprocessing time for sound speed measurement and image diagnosis can be greatly shortened compared to the conventional case. Furthermore, since the sound speed of the object to be measured is calculated by electrically correcting the inclination of the mounting surface, there is no error due to misadjustment like mechanical adjustment. Each time, the obtained data has almost no variation, and the accuracy of the calculated sound speed is remarkably improved. Therefore, the sound velocity of the object to be measured can be measured with high accuracy and in a very short time, and the sound velocity image reflecting the object to be measured can be obtained more accurately.

この発明に係る音速測定では、被測定物の音速を極めて高い精度で測定できるので、この考えをこの発明に係る画像診断に応用することは極めて重要である。画像診断の対象となる生体組織の切片は、その厚みが10μm程度と非常に薄いから、数μm単位の僅かな傾きでも音速像算出結果には甚大な影響を及ぼすからである。   In the sound speed measurement according to the present invention, the sound speed of the object to be measured can be measured with extremely high accuracy. Therefore, it is extremely important to apply this idea to the image diagnosis according to the present invention. This is because a slice of a living tissue to be subjected to image diagnosis has a very thin thickness of about 10 μm, and even a slight inclination of several μm has a great influence on the sound velocity image calculation result.

続いて、この発明に係る超音波顕微鏡を使用した音速測定方法、音速測定装置、画像診断方法およびその画像診断装置の好ましい実施例を図面を参照して詳細に説明する。以下に示す実施例としては、被測定物として生体組織の切片を用いてその厚みと音速を測定することで、生体組織の音速像を得るものに適用した例を述べる。   Subsequently, preferred embodiments of a sound speed measuring method, a sound speed measuring apparatus, an image diagnostic method and an image diagnostic apparatus using the ultrasonic microscope according to the present invention will be described in detail with reference to the drawings. As an example described below, an example will be described in which a slice of a living tissue is used as an object to be measured and its thickness and sound speed are measured to obtain a sound velocity image of the living tissue.

被測定物は生体組織切片である。図17は傾き補正を電気的に行うための原理的説明に供する図である。
機械的調整によってプレパラート12の傾きを調整するのではなく、電気的な処理でこの機械的な傾きを補正するには、まずプレパラート12の平面方程式を解く必要がある。X−Yステージ16の傾きはプレパラート12の傾きに反映されるので、結局は生体組織切片13を載置しているプレパラート12の傾きを補正すればよいことになる。
The object to be measured is a biological tissue section. FIG. 17 is a diagram for explaining the principle for electrically performing tilt correction.
Rather than adjusting the inclination of the preparation 12 by mechanical adjustment, in order to correct this mechanical inclination by electrical processing, it is necessary to first solve the plane equation of the preparation 12. Since the inclination of the XY stage 16 is reflected in the inclination of the preparation 12, it is only necessary to correct the inclination of the preparation 12 on which the biological tissue section 13 is placed.

プレパラート12が傾いた状態での任意の点での平面方程式が算出できれば、その平面方程式によって表されるz軸の値そのものが傾きを含んだ値となっているので、生体組織切片13の厚みdをこの平面方程式を利用して算出するだけで、傾きを考慮した厚み算出となる。つまり、傾きによる影響を排除して厚みが算出されたことになる。   If the plane equation at an arbitrary point in the state where the preparation 12 is tilted can be calculated, the z-axis value itself represented by the plane equation is a value including the tilt. Is calculated using this plane equation, and the thickness is calculated in consideration of the inclination. That is, the thickness is calculated excluding the influence of the inclination.

そのため、図17のように生体組織切片13が載置されていないプレパラート12の平面12aにおける任意の3点Po,P1,P2を指定する。そのうちの1点例えばPo(xo,yo)を参照点とする。残りの2点P1,P2の座標を(x1,y1)、(x2,y2)とする。   Therefore, arbitrary three points Po, P1, and P2 on the plane 12a of the preparation 12 on which the biological tissue section 13 is not placed are designated as shown in FIG. One of them, for example, Po (xo, yo) is used as a reference point. The coordinates of the remaining two points P1 and P2 are (x1, y1) and (x2, y2).

平面方程式は周知のように、
Z=ax+by+c ・・・・・(1)
(ここに、a,b,cは係数である)
である。したがって3点Po,P1,P2の平面方程式Zo,Z1,Z2は
Zo=axo+byo+c ・・・・・(2)
Z1=ax1+by1+c ・・・・・(3)
Z2=ax2+by2+c ・・・・・(4)
参照点Poの傾きは、基準点であるためゼロ(Zo=0)とする。
(2)〜(4)式から、係数a,b,cは、
As is well known, the plane equation is
Z = ax + by + c (1)
(Where a, b, and c are coefficients)
It is. Therefore, the plane equations Zo, Z1, Z2 of the three points Po, P1, P2 are: Zo = axo + byo + c (2)
Z1 = ax1 + by1 + c (3)
Z2 = ax2 + by2 + c (4)
The inclination of the reference point Po is zero (Zo = 0) because it is a reference point.
From the equations (2) to (4), the coefficients a, b, and c are

Figure 0004521584
Figure 0004521584

となる。 It becomes.

したがって、生体組織切片13が載置されているプレパラート12の任意の測定点Px(xp、yp)での平面方程式は、上述した係数a,b,cから
Zp=axp+byp+c ・・・・・(6)
となる。このz軸方向の値Zp、すなわち参照点Poと測定点Pxとのz軸方向における距離差分Zpを傾き補正値として、その測定点Pxにおける厚みdpが算出される。
Accordingly, the plane equation at an arbitrary measurement point Px (xp, yp) of the preparation 12 on which the biological tissue section 13 is placed is expressed by the following equations from the coefficients a, b, c: Zp = axp + byp + c (6 )
It becomes. Using the value Zp in the z-axis direction, that is, the distance difference Zp in the z-axis direction between the reference point Po and the measurement point Px, the thickness dp at the measurement point Px is calculated.

上述した測定点P1,P2での平面方程式Z1,Z2は以下のようにして算出できる。まず、各点Po,P1,P2での反射波信号Sro、Sr1,Sr2を得る。   The plane equations Z1 and Z2 at the measurement points P1 and P2 described above can be calculated as follows. First, the reflected wave signals Sro, Sr1, Sr2 at the points Po, P1, P2 are obtained.

次に、これら反射波信号Sro、Sr1,Sr2の周波数成分を得るためフーリエ変換(例えば高速フーリエ変換FFT)処理を施す。それぞれのフーリエ変換出力Fo、F1,F2から参照点Poを基準にしたフーリエ変換出力の比(F1/Fo)および(F2/Fo)を求める。このフーリエ変換出力の比は規格化スペクトルとなる。   Next, in order to obtain frequency components of these reflected wave signals Sro, Sr1, and Sr2, a Fourier transform (for example, a fast Fourier transform FFT) process is performed. Ratios (F1 / Fo) and (F2 / Fo) of Fourier transform outputs based on the reference point Po are obtained from the respective Fourier transform outputs Fo, F1, and F2. The ratio of the Fourier transform output is a normalized spectrum.

次に、この規格化スペクトルに関連した複素平面(周波数に対する位相の関係)より、φ/2πf(fは超音波信号の中心周波数)を求める。このφ/2πfは、結局参照点Poと測定点P1との間の傾きによって生ずる時間差Δto1となる。同様に参照点Poと測定点P2との傾きによって生ずる時間差Δto2が求まる。   Next, φ / 2πf (f is the center frequency of the ultrasonic signal) is obtained from a complex plane (phase relationship with respect to frequency) related to the normalized spectrum. This φ / 2πf eventually becomes the time difference Δto1 caused by the inclination between the reference point Po and the measurement point P1. Similarly, a time difference Δto2 caused by the inclination between the reference point Po and the measurement point P2 is obtained.

したがって、参照点Poと測定点P1までのz軸方向の値が、すなわち平面方程式Z1であり、これは超音波振動子10とプレパラート12との間に介在された媒質(通常は水)の音速をCoとしたとき、
Z1=Co・Δto1 ・・・・(7)
として求めることができる。Coは既知の値(=1600m/s)である。同様に、測定点P2における平面方程式Z2は、
Z2=Co・Δto2 ・・・・(8)
となり、これらの値Z1,Z2および上述した参照値Zo(=0)から上述した係数a,b,cを算出できる(式(5))。算出されたこれら係数a,b,cは生体組織切片13の厚み測定時に参照されるために、メモリ手段(RAMなど)に保存される。
Accordingly, the value in the z-axis direction from the reference point Po to the measurement point P1, that is, the plane equation Z1, is the sound velocity of the medium (usually water) interposed between the ultrasonic transducer 10 and the preparation 12. When Co is
Z1 = Co · Δto1 (7)
Can be obtained as Co is a known value (= 1600 m / s). Similarly, the plane equation Z2 at the measurement point P2 is
Z2 = Co · Δto2 (8)
Thus, the above-described coefficients a, b, c can be calculated from these values Z1, Z2 and the above-described reference value Zo (= 0) (formula (5)). The calculated coefficients a, b, and c are stored in a memory means (such as a RAM) for reference when measuring the thickness of the biological tissue section 13.

上述したように、生体組織切片13に超音波信号を照射すると、生体組織切片13の表面と背面からの反射波が干渉した干渉反射波信号Soが得られる。そこで、プレパラート12の表面からの参照反射波信号Srと干渉反射波信号Soとをそれぞれフーリエ変換して、参照反射波信号Srと干渉反射波信号So(具体的には一次反射波信号Su1)と比較することで、規格化した強度スペクトルと同じく規格化した位相スペクトルを得ることができる。このスペクトル例を図18A,図18Bに示す。   As described above, when an ultrasonic signal is applied to the biological tissue section 13, an interference reflected wave signal So in which reflected waves from the front and back surfaces of the biological tissue section 13 interfere is obtained. Therefore, the reference reflected wave signal Sr and the interference reflected wave signal So from the surface of the preparation 12 are Fourier transformed, respectively, and the reference reflected wave signal Sr and the interference reflected wave signal So (specifically, the primary reflected wave signal Su1) and By comparison, a normalized phase spectrum can be obtained as well as a normalized intensity spectrum. Examples of this spectrum are shown in FIGS. 18A and 18B.

そして、信号強度の極小点または極大点の周波数をfm、そのときの位相をφmとすると、生体組織切片13の表面と背面からの反射波は極小点では逆位相、極大点では同位相となっている。すなわち極小点においては生体組織切片13の表面からの反射波は背面からの反射波より(2n−1)πだけ進んでおり、{φm+(2n−1)π}となる(nは自然数)。生体組織切片13の厚さをd、水の音速をCoとすれば、   If the frequency of the signal intensity minimum point or maximum point is fm, and the phase at that time is φm, the reflected waves from the front and back surfaces of the biological tissue section 13 are opposite in phase at the minimum point and in phase at the maximum point. ing. That is, at the minimum point, the reflected wave from the surface of the biological tissue section 13 is advanced by (2n−1) π from the reflected wave from the back surface, and becomes {φm + (2n−1) π} (n is a natural number). If the thickness of the biological tissue section 13 is d and the sound speed of water is Co,

Figure 0004521584
Figure 0004521584

極大点においてはこの位相差が2nπとなっていることから、 Since this phase difference is 2nπ at the maximum point,

Figure 0004521584
Figure 0004521584

が成立している。一方、2dの距離を生体組織切片13内の音速Cで通過した反射波と、参照音速Coで通過した反射波の位相差がφmであるから、 Is established. On the other hand, the phase difference between the reflected wave that has passed through the distance of 2d at the sound velocity C in the biological tissue slice 13 and the reflected wave that has passed at the reference sound velocity Co is φm.

Figure 0004521584
Figure 0004521584

(9)式と(11)式を連立させて解くと、
d={φm+(2n−1)π}Co/4πfm ・・・(12)
のようにして生体組織切片13の厚さdを求めることができる。また、(10)式と(11)式を連立させて解くと、
Solving equations (9) and (11) simultaneously,
d = {φm + (2n−1) π} Co / 4πfm (12)
Thus, the thickness d of the biological tissue section 13 can be obtained. Also, solving (10) and (11) together,

Figure 0004521584
Figure 0004521584

のように生体組織の音速Cが求まる。 Thus, the sound velocity C of the living tissue is obtained.

一方、生体組織切片13の厚みdと音速Cを演算するときの周波数faは、図18Aに示す規格化強度スペクトルの最小値fm(この場合98MHz)を採用する。またその最小値fmでの位相φmは、図18Bの規格化位相スペクトル曲線より92degとなる。   On the other hand, the minimum value fm (98 MHz in this case) of the normalized intensity spectrum shown in FIG. 18A is adopted as the frequency fa when calculating the thickness d and the sound velocity C of the biological tissue section 13. The phase φm at the minimum value fm is 92 deg from the normalized phase spectrum curve of FIG. 18B.

厚みdと音速Cを演算するときには、fm=faとして代入するが、そのときの位相はリファレンス点となるプレパラート12上の参照点Poからの傾きによって生じたz軸方向(超音波出射方向)の距離差分Zpに相当する位相だけずれていることになる。したがって、この距離差分Zpに相当する位相分を補正する必要がある。距離差分Zpは(6)式によって算出できる。   When calculating the thickness d and the sound velocity C, fm = fa is substituted, but the phase at that time is in the z-axis direction (ultrasonic emission direction) generated by the inclination from the reference point Po on the preparation 12 serving as the reference point. The phase is shifted by a phase corresponding to the distance difference Zp. Therefore, it is necessary to correct the phase corresponding to the distance difference Zp. The distance difference Zp can be calculated by equation (6).

そこで、この距離差分Zpだけ補正した位相φmを以下のように演算する。   Therefore, the phase φm corrected by the distance difference Zp is calculated as follows.

Figure 0004521584
Figure 0004521584

このようにして算出した位相φmと、極小値での周波数fm(=fa)を、(12)式と(13)式に代入することで、距離差分Zpを補正した最終的な厚みdと音速Cを得ることができる。   By substituting the phase φm calculated in this way and the frequency fm (= fa) at the minimum value into the equations (12) and (13), the final thickness d and the sound velocity corrected for the distance difference Zp. C can be obtained.

続いて、この発明に係る超音波顕微鏡を使用した音速測定装置をさらに応用した画像診断装置100の一例を図1以下を参照して説明する。   Next, an example of an image diagnostic apparatus 100 to which the sound velocity measuring apparatus using the ultrasonic microscope according to the present invention is further applied will be described with reference to FIG.

この画像診断装置100は大別して超音波顕微鏡1の本体部1Aと、第1の信号処理部1Bと第2の信号処理部1Cとで構成される。この画像診断装置100は、当然のことながら音速測定装置としても機能する。相違する点は、音速測定装置として使用する場合には、画像診断における最終工程である音速像の算出処理部57が含まれないことである。   The diagnostic imaging apparatus 100 is roughly composed of a main body 1A of the ultrasonic microscope 1, a first signal processing unit 1B, and a second signal processing unit 1C. Naturally, the diagnostic imaging apparatus 100 also functions as a sound velocity measuring apparatus. The difference is that, when used as a sound velocity measuring device, the sound velocity image calculation processing unit 57 which is the final step in image diagnosis is not included.

ここで、この画像診断装置100は、超音波振動子11と対峙する試料載置板12の傾き量を検出する傾き量検出手段と、検出されたこの傾き量と、超音波振動子11からの反射波信号Soとに基づいて被測定物13の厚みを算出する厚み算出手段と、反射波信号Soに基づいて算出された被測定物13の厚みを傾き量によって電気的に補正する補正手段と、被測定物13の厚み情報から当該被測定物13の音速を算出する音速算出手段と、この音速から被測定物13の音速像を生成する音速像生成手段とで構成されていることになる。   Here, the diagnostic imaging apparatus 100 includes an inclination amount detection unit that detects an inclination amount of the sample mounting plate 12 facing the ultrasonic transducer 11, the detected inclination amount, and the ultrasonic transducer 11 Thickness calculating means for calculating the thickness of the object to be measured 13 based on the reflected wave signal So, and correcting means for electrically correcting the thickness of the object to be measured 13 calculated based on the reflected wave signal So by the amount of inclination. The sound velocity calculating means for calculating the sound velocity of the object to be measured 13 from the thickness information of the object to be measured 13 and the sound velocity image generating means for generating the sound velocity image of the object to be measured 13 from this sound velocity. .

傾き量検出手段、厚み算出手段、傾き量補正手段及び音速像生成手段は、主として本体部1Aと第2の信号処理部1Cとによって構成される。   The inclination amount detection means, the thickness calculation means, the inclination amount correction means, and the sonic image generation means are mainly composed of the main body 1A and the second signal processing part 1C.

顕微鏡本体部1Aは、上述したように被測定物であるこの例では生体組織切片13を二次元的に動かすために二次元走査手段として使用されるX−Yステージ16が超音波振動子11と正対する。X−Yステージ16はそれぞれのステージ16X,16Yを駆動するモータ31X,33Yが設けられている。駆動モータ31X,33Yはステッピングモータやリニアモータを使用することができ、好ましくは高速かつ低振動のリニアモータが好適である。   In this example, which is the object to be measured as described above, the microscope main body 1A has an XY stage 16 used as a two-dimensional scanning means to move the biological tissue section 13 two-dimensionally with the ultrasonic transducer 11. Face up. The XY stage 16 is provided with motors 31X and 33Y for driving the stages 16X and 16Y. A stepping motor or a linear motor can be used as the drive motors 31X and 33Y, and a high-speed and low-vibration linear motor is preferable.

ステッピングモータ、特にリニアモータを使用することで、Xステージ16Xを連続走査(連続送り)することが可能となり、そして、Yステージ16Yを間欠送りとなるように制御することで、X−Yステージ16の高速走査が可能になる。   By using a stepping motor, particularly a linear motor, the X stage 16X can be continuously scanned (continuous feed), and by controlling the Y stage 16Y to be intermittent feed, the XY stage 16 can be controlled. High speed scanning becomes possible.

駆動モータ31X,33Yにはそれぞれこれらに対する駆動信号を生成するためのコントローラ30X,32Yが設けられ、一方のコントローラ30Xには第2の信号処理部1Cで生成された駆動制御信号(パルス信号)が供給される。   The drive motors 31X and 33Y are provided with controllers 30X and 32Y for generating drive signals for them, respectively, and the drive control signal (pulse signal) generated by the second signal processing unit 1C is provided in one controller 30X. Supplied.

この例ではXステージ16Xに関連してエンコーダ34が設けられ、このエンコーダ34によってXステージ16Xの走査位置が検出される。例えば、図13のようにプレパラート12を含めて300×300個の測定点(ピクセル)に分割したときには、1回のx方向における走査が300分割されることになるので、それぞれのピクセルの位置がこのエンコーダ34によって検出され、このエンコーダ出力に同期して駆動制御信号が生成される。   In this example, an encoder 34 is provided in association with the X stage 16X, and the scanning position of the X stage 16X is detected by the encoder 34. For example, when divided into 300 × 300 measurement points (pixels) including the preparation 12 as shown in FIG. 13, one scan in the x direction is divided into 300, so that the position of each pixel is A drive control signal is generated in synchronization with the encoder output detected by the encoder 34.

そのため、エンコーダ出力が第2の信号処理部1C内に設けられたCPU55に同期信号として供給されており、またCPU55にはパルス生成部56が関連され、したがってパルス生成部56ではエンコーダ出力に同期した駆動制御信号が生成される。   Therefore, the encoder output is supplied as a synchronization signal to the CPU 55 provided in the second signal processing unit 1C, and the pulse generation unit 56 is associated with the CPU 55. Therefore, the pulse generation unit 56 is synchronized with the encoder output. A drive control signal is generated.

Xステージ16X用のコントローラ30Xからはさらに駆動制御信号に同期したトリガー信号Stが生成されて、これがパルサー部36に供給される。パルサー部36ではこのトリガー信号Stに同期したタイミングで励起パルスが生成される。   A trigger signal St synchronized with the drive control signal is further generated from the controller 30X for the X stage 16X, and this is supplied to the pulsar section 36. In the pulsar unit 36, an excitation pulse is generated at a timing synchronized with the trigger signal St.

パルサー部36は超音波振動子11に対する励起パルスを生成するためのもので、この励起パルスは送受波分離部37を介して超音波振動子11に供給される。パルサー部36からの励起パルスSiは広帯域幅のパルス信号であって(図12参照)、この例では図11Aに示すように1サイクルの正弦波状をなすパルス信号が出力されるように構成されている。   The pulsar unit 36 is for generating an excitation pulse for the ultrasonic transducer 11, and this excitation pulse is supplied to the ultrasonic transducer 11 via the transmission / reception wave separation unit 37. The excitation pulse Si from the pulsar unit 36 is a broadband pulse signal (see FIG. 12). In this example, as shown in FIG. 11A, a pulse signal having a sine wave shape of one cycle is output. Yes.

パルサー部36は後述するように、トリガー信号Stによって間欠的に駆動されるが、この単位間欠周期の中で、さらに集中的に複数個の励起パルスが短時間に出力できるように構成されている。   As will be described later, the pulsar unit 36 is intermittently driven by the trigger signal St, and is configured so that a plurality of excitation pulses can be output more intensively in this unit intermittent period in a short time. .

超音波振動子11は送受波兼用の振動子が使用され、励起パルスによって励起されることで、ビーム状に絞られた超音波信号Siが出射される。同じ超音波振動子11で受波した反射波信号SrやSoのうち、生体組織切片13から得られる反射波信号Soは干渉波形である。これら反射波信号Sr,Soは送受波分離部37および受信部38を介して、信号演算部として機能する第1の信号処理部1Bに供給される。   As the ultrasonic transducer 11, a transducer for both transmission and reception is used, and when excited by an excitation pulse, an ultrasonic signal Si narrowed in a beam shape is emitted. Of the reflected wave signals Sr and So received by the same ultrasonic transducer 11, the reflected wave signal So obtained from the biological tissue section 13 is an interference waveform. These reflected wave signals Sr and So are supplied to the first signal processing unit 1B functioning as a signal calculation unit via the transmission / reception wave separation unit 37 and the reception unit 38.

第1の信号処理部1Bはゲート回路40とゲートパルス生成回路42を有する。ゲートパルス生成回路42には励起パルスに同期したパルス信号が供給される。そのため、パルサー部36の出力段には抵抗分圧部41が設けられ、抵抗分圧部41で分圧された励起パルスSiがゲートパルス生成回路42に供給され、励起パルスSiより所定時間遅延したゲートパルスSgが生成される。   The first signal processing unit 1B includes a gate circuit 40 and a gate pulse generation circuit 42. The gate pulse generation circuit 42 is supplied with a pulse signal synchronized with the excitation pulse. Therefore, a resistance voltage dividing unit 41 is provided at the output stage of the pulsar unit 36, and the excitation pulse Si divided by the resistance voltage dividing unit 41 is supplied to the gate pulse generation circuit 42 and delayed by a predetermined time from the excitation pulse Si. A gate pulse Sg is generated.

一方、ゲート回路40には参照反射波信号Srと干渉反射波信号Soが供給され、上述したゲートパルスSgで干渉反射波信号Soのうち所望の信号部分のみがゲートされる。   On the other hand, the reference reflected wave signal Sr and the interference reflected wave signal So are supplied to the gate circuit 40, and only the desired signal portion of the interference reflected wave signal So is gated by the gate pulse Sg described above.

ここで、干渉反射波信号Soは図11Cにも示すように超音波振動子11を励起する励起パルスSiに関連した反射波信号Si’とその反射波信号Suとが合成されたものである。超音波の反射波は超音波振動子11と生体組織切片13との間で反射が繰り返されることになるから、本来ではn次の反射波信号Suが多重信号として得られることになる。   Here, as shown in FIG. 11C, the interference reflected wave signal So is a combination of the reflected wave signal Si ′ related to the excitation pulse Si for exciting the ultrasonic transducer 11 and the reflected wave signal Su. Since the reflected ultrasonic wave is repeatedly reflected between the ultrasonic transducer 11 and the biological tissue section 13, the nth-order reflected wave signal Su is originally obtained as a multiplexed signal.

実施例では、そのうち最初に反射する反射波信号(一次反射波信号Su1という。)のみを抽出するため、励起パルスSiから所定時間だけ遅延したゲートパルスSgが使用されるものである。   In the embodiment, in order to extract only the reflected wave signal that is reflected first (referred to as the primary reflected wave signal Su1), the gate pulse Sg delayed by a predetermined time from the excitation pulse Si is used.

ゲートされた一次反射波信号Su1はA/D変換回路43に供給されてA/D変換処理されると共に、処理されたA/D変換出力がメモリ手段44に一時的に保存される。後述するように、超音波振動子11は所定周期で間欠的に励起されるが、この単位間欠周期の中で、さらに集中的に複数個の励起パルスが短時間に励起される。詳細は後述するとして、この例では7〜10回程度高速励起される。以下では8回高速励起した場合を示す。   The gated primary reflected wave signal Su1 is supplied to the A / D conversion circuit 43 and subjected to A / D conversion processing, and the processed A / D conversion output is temporarily stored in the memory means 44. As will be described later, the ultrasonic transducer 11 is intermittently excited in a predetermined cycle, and a plurality of excitation pulses are excited in a short time in a more concentrated manner in this unit intermittent cycle. As will be described in detail later, in this example, high-speed excitation is performed about 7 to 10 times. In the following, a case where high-speed excitation is performed eight times is shown.

そして、それぞれがメモリ手段44に格納されると共に、8個目の一次反射波信号Su(8)を取得した後、後段の平均化回路45でこれら一次反射波信号Su1(1)、Su2(2)、・・・Su8(8)を用いて平均化される。このような平均化処理を施すのは、1つの一次反射波信号Su1だけでは充分なS/N(つまりC/N)が得られないためであり、S/Nを改善することが主目的である。   Each of them is stored in the memory means 44, and after the eighth primary reflected wave signal Su (8) is acquired, the primary reflected wave signals Su1 (1) and Su2 (2) are obtained by the averaging circuit 45 in the subsequent stage. ),..., And averaged using Su8 (8). Such an averaging process is performed because a sufficient S / N (that is, C / N) cannot be obtained with only one primary reflected wave signal Su1, and the main purpose is to improve the S / N. is there.

また、上述したようにパルサー部38の出力段側に設けられた抵抗分圧部41で分圧された励起パルスSiを用いてゲートパルスSgを生成したのは、8回高速励起して8つの一次反射波信号Su1を得るときの、ゲートパルスSgのジッタを回避するためである。ゲートパルスSgにジッタがなければ、一次反射波信号Su1を平均化したときの同期ジッタ成分がゼロとなるため、安定した平均化出力を得ることができるからである。   In addition, as described above, the gate pulse Sg is generated using the excitation pulse Si divided by the resistance voltage dividing unit 41 provided on the output stage side of the pulsar unit 38. This is to avoid the jitter of the gate pulse Sg when the primary reflected wave signal Su1 is obtained. This is because if there is no jitter in the gate pulse Sg, the synchronization jitter component when the primary reflected wave signal Su1 is averaged becomes zero, and a stable averaged output can be obtained.

なお、このような目的で使用できる高速処理可能なゲート機能付きA/Dボードとしては汎用の高速A/Dボード(例えば販売先が横河製作所である型番741025などのA/Dボード)を使用できる。   Note that a general-purpose high-speed A / D board (for example, an A / D board such as model number 741025 sold by Yokogawa Seisakusho) is used as an A / D board with a gate function capable of high-speed processing that can be used for such purposes. it can.

平均化処理された一次反射波信号Suは第2の信号処理部1Cに供給される。第2の信号処理部1Cでは、音速像を得るための前処理として、平面方程式のZの算出処理や係数a,b,cを算出するための補正値算出処理が行われる。その他に音速の算出処理および音速像の生成処理が行われる。   The averaged primary reflected wave signal Su is supplied to the second signal processing unit 1C. In the second signal processing unit 1C, as preprocessing for obtaining a sound velocity image, Z calculation processing of a plane equation and correction value calculation processing for calculating coefficients a, b, and c are performed. In addition, sound speed calculation processing and sound speed image generation processing are performed.

そのため、この第2の信号処理部1Cでは、一次反射波信号Su1が供給される高速フーリエ変換処理部51、フーリエ変換出力を演算処理する演算部52、演算処理された規格化スペクトルから時間差Δtを算出するΔt算出部53およびこれらの処理を行うときのワーキング用メモリ手段54などが設けられる。   Therefore, in the second signal processing unit 1C, the fast Fourier transform processing unit 51 to which the primary reflected wave signal Su1 is supplied, the computing unit 52 that computes the Fourier transform output, and the time difference Δt from the computed normalized spectrum. A Δt calculating unit 53 for calculating, a working memory means 54 for performing these processes, and the like are provided.

さらに、時間差Δtのデータやフーリエ変換出力がそれぞれ供給される音速の算出と音速像生成処理部57が設けられる。この音速算出および音速像生成処理部57は画像描画部として機能することになる。時間差Δtは係数a,b,cを算出するときに供給されるデータであり、フーリエ変換出力は通常の音速像を生成するときに供給されるデータである。   Further, a sound speed calculation and sound speed image generation processing unit 57 to which data of the time difference Δt and Fourier transform output are respectively supplied are provided. The sound speed calculation and sound speed image generation processing unit 57 functions as an image drawing unit. The time difference Δt is data supplied when calculating the coefficients a, b, and c, and the Fourier transform output is data supplied when generating a normal sound velocity image.

そして、CPU55内のROM60(図2参照)にストアされているz値算出のための処理プログラムや音速算出用処理プログラムさらには音速像算出用処理プログラムを実行することによって、上述したz値算出処理などが行われることになる。   Then, by executing a processing program for calculating the z value, a processing program for calculating the sound speed, and a processing program for calculating the sound speed stored in the ROM 60 (see FIG. 2) in the CPU 55, the above-described z value calculating process is executed. Etc. will be performed.

このCPU55にはさらにコントローラ30Xに供給するためのパルス生成部56が関連され、CPU55からの同期パルスに同期してパルス生成部56が駆動される。CPU55にはさらに顕微鏡本体部1Aに設けられたエンコーダ34のエンコーダ出力が供給される。このエンコーダ出力に同期してCPU55から同期パルスが出力される。   The CPU 55 is further associated with a pulse generation unit 56 for supplying to the controller 30X, and the pulse generation unit 56 is driven in synchronization with the synchronization pulse from the CPU 55. The CPU 55 is further supplied with the encoder output of the encoder 34 provided in the microscope main body 1A. A synchronization pulse is output from the CPU 55 in synchronization with the encoder output.

したがって最終的にはエンコーダ34のエンコーダ出力に同期してトリガー信号Stが生成されるため、生体組織切片13の各ピクセルに同期してパルサー部36が励起されるので、確実にそれぞれのピクセルごとに超音波信号を照射できるようになる。   Therefore, since the trigger signal St is finally generated in synchronization with the encoder output of the encoder 34, the pulsar unit 36 is excited in synchronization with each pixel of the biological tissue section 13, so that each pixel is reliably detected. An ultrasonic signal can be irradiated.

第2の信号処理部1Cの各部構成の関係は、コンピュータを中心にすると、図2のようにも書き替えることができる。上述したワーキング用のメモリ手段54は図2のRAMが使用されることになる。インタフェース61を介して一次反射波信号(デジタル信号)Su1が取り込まれる。   The relationship between the components of the second signal processing unit 1C can be rewritten as shown in FIG. The working memory means 54 uses the RAM shown in FIG. The primary reflected wave signal (digital signal) Su1 is taken in via the interface 61.

画像描画部としても機能するこの音速像算出処理部57では、等音速線によって描画されることで、生体組織切片13の組織がミクロン単位で描画される。同一音速同士をつなぎ合わせることで得られる等音速線をモニタ(図示せず)に描画することで生体組織の可視像(音速像)が得られる。この音速像から生体組織を観察し、診断できる。   In the sonic image calculation processing unit 57 that also functions as an image drawing unit, the tissue of the biological tissue section 13 is drawn in units of microns by drawing with isosonic lines. A visible image (sonic image) of a living tissue can be obtained by drawing on a monitor (not shown) an isosonic line obtained by connecting the same sound velocities. A biological tissue can be observed and diagnosed from this sound velocity image.

等音速線で囲まれる領域を同じ色で表示することでカラー表示が可能になり、これによって変質した生体組織を明確に識別できるようになる。表示色にグラデーションをかけることで組織の境界遷移がスムーズになる。もちろん、一次反射波信号Su1の減衰度(強度)を基準にして描画すると生体像として減衰像が得られる。   By displaying the area surrounded by the isosonic line in the same color, color display becomes possible, and this makes it possible to clearly identify the altered biological tissue. By applying gradation to the display color, the boundary transition of the organization becomes smooth. Of course, when drawing is performed based on the attenuation (intensity) of the primary reflected wave signal Su1, an attenuation image is obtained as a living body image.

図2では、第2の信号処理部1C用としてCPU55を説明したが、装置全体としての制御を司るCPUを別に設ける他に、このCPU55を装置全体の制御を司るCPUとして兼用できることはもちろんである。   In FIG. 2, the CPU 55 has been described for the second signal processing unit 1 </ b> C. However, it is a matter of course that the CPU 55 can also be used as a CPU for controlling the entire apparatus, in addition to providing a separate CPU for controlling the entire apparatus. .

続いて、図1に示した画像診断装置100の動作を図3以下を参照して説明する。
図3は一次反射波信号Su1のS/Nの説明図である。上述したように生体組織切片13を高速で二次元走査して最終的には計測を開始してから、5分、就中2分以下で音速像を描画できるようにするには、x方向への走査を高速連続走査して1ラインの走査時間を高速化する必要がある。
Next, the operation of the diagnostic imaging apparatus 100 shown in FIG. 1 will be described with reference to FIG.
FIG. 3 is an explanatory diagram of S / N of the primary reflected wave signal Su1. As described above, in order to be able to draw a sound velocity image within 5 minutes, especially 2 minutes or less after the biological tissue section 13 is two-dimensionally scanned at high speed and finally starts measurement, in the x direction. It is necessary to increase the scanning time of one line by performing high-speed continuous scanning.

この例では、1ピクセルに対する計測時間が80μsec程度となるからトリガー信号Stの間隔も80μsecとなる(図3A)。   In this example, since the measurement time for one pixel is about 80 μsec, the interval of the trigger signal St is also 80 μsec (FIG. 3A).

トリガー信号Stに同期して励起パルスSiが得られる(図3B)。そしてその反射波は図3Cのようになる。励起パルスSiに関連した反射波信号Si’に対して大凡2μsec遅れて一次反射波信号Su1が、さらに2μsec遅れて二次反射波信号Su2が得られる。反射波信号Soの幅は100nsec程度である。   An excitation pulse Si is obtained in synchronization with the trigger signal St (FIG. 3B). The reflected wave is as shown in FIG. 3C. The primary reflected wave signal Su1 is obtained with a delay of about 2 μsec with respect to the reflected wave signal Si ′ related to the excitation pulse Si, and the secondary reflected wave signal Su2 is obtained with a delay of 2 μsec. The width of the reflected wave signal So is about 100 nsec.

最も単純な信号処理は、図3DのゲートパルスSgを用いて一次反射波信号Su1を抽出する処理であるが、これでは充分なS/Nが得られない。
S/Nを改善するには、複数回超音波振動子11を励起して同じ信号を得、これを平均化すればよい。しかし、そうするためには同じ走査位置で複数回の励起処理を行わなければならないので、次のピクセルに対する計測開始までに複数倍の時間がかかってしまう。これでは、高速画像処理を実現できない。
The simplest signal processing is processing for extracting the primary reflected wave signal Su1 using the gate pulse Sg of FIG. 3D, but this does not provide sufficient S / N.
In order to improve the S / N, the ultrasonic transducer 11 may be excited a plurality of times to obtain the same signal and averaged. However, in order to do so, excitation processing must be performed a plurality of times at the same scanning position, so that it takes a multiple of times before the start of measurement for the next pixel. With this, high-speed image processing cannot be realized.

そこで、超音波振動子11を間欠的に励起する時間間隔を保ったまま、さらにこの間欠励起時間内で集中的に複数個の励起パルスを短時間に励起することで、短時間に複数個の一次反射波信号を取得し、その平均化した信号を反射波信号として用いるようにする。具体例を図4を参照して説明する。   Therefore, while maintaining a time interval for intermittently exciting the ultrasonic transducer 11, a plurality of excitation pulses are intensively excited within a short period of time within the intermittent excitation time. The primary reflected wave signal is acquired, and the averaged signal is used as the reflected wave signal. A specific example will be described with reference to FIG.

上述したように励起パルスSiに対する反射波信号Si’から一次反射波信号Su1が得られるまでの時間は大凡2μsecであり、この一次反射波信号Su1が最もS/Nがよく、二次以降の反射波信号Su2、・・・は超音波振動子11と生体組織切片13との間の多重反射信号であるためS/Nが悪い。   As described above, the time until the primary reflected wave signal Su1 is obtained from the reflected wave signal Si ′ with respect to the excitation pulse Si is about 2 μsec, and this primary reflected wave signal Su1 has the best S / N, and the secondary and subsequent reflections. Since the wave signal Su2,... Is a multiple reflection signal between the ultrasonic transducer 11 and the biological tissue section 13, the S / N is poor.

したがって、図4Aのように高速画像処理を実現するため間欠励起時間は図3の場合と同じにする。次に、使用する反射波信号としては一次反射波信号Su1のみとする。この反射波信号の取捨を行うため、ほぼ2μsec間隔で複数回(2回以上)、好ましくは7〜10回、例えば8回に亘り連続的に超音波振動子11を励起する。すなわち、図4Aに示すようにトリガー信号Stに対して、パルサー部36では図4Bに示すような2μsec間隔の励起パルスが8個連続的に生成される。   Therefore, the intermittent excitation time is made the same as in FIG. 3 in order to realize high-speed image processing as shown in FIG. 4A. Next, the primary reflected wave signal Su1 is the only reflected wave signal to be used. In order to dispose of the reflected wave signal, the ultrasonic transducer 11 is excited continuously several times (2 times or more), preferably 7 to 10 times, for example, 8 times at intervals of about 2 μsec. That is, as shown in FIG. 4A, for the trigger signal St, the pulsar unit 36 continuously generates eight excitation pulses at intervals of 2 μsec as shown in FIG. 4B.

この2μsec間隔の励起パルスSiで超音波振動子11が励起される。その結果、図4Cに示すような反射波信号が得られる。   The ultrasonic transducer 11 is excited by the excitation pulses Si at intervals of 2 μsec. As a result, a reflected wave signal as shown in FIG. 4C is obtained.

その一方で、励起パルスSiに同期したゲートパルスSg(図4D)が生成される。この連続的なゲートパルスSgで一次反射波信号Su1のみをゲートすれば、8個の一次反射波信号Su1が得られる(図4E)。図4Eは便宜的に並べて図示したものである。実際にはこのような時系列で信号が並ぶものではないことは明らかである。   On the other hand, a gate pulse Sg (FIG. 4D) synchronized with the excitation pulse Si is generated. If only the primary reflected wave signal Su1 is gated by this continuous gate pulse Sg, eight primary reflected wave signals Su1 are obtained (FIG. 4E). FIG. 4E is shown side by side for convenience. It is clear that the signals are not actually arranged in such a time series.

これら8個の一次反射波信号Su1(1)〜Su8(8)を平均化処理して、最終的な一次反射波信号Suが生成される(図4F)。この処理によって一次反射波信号SuのS/Nが改善されると共に、高速画像処理が可能になる。   These eight primary reflected wave signals Su1 (1) to Su8 (8) are averaged to generate a final primary reflected wave signal Su (FIG. 4F). This process improves the S / N of the primary reflected wave signal Su and enables high-speed image processing.

続いて、この発明に係る画像診断方法を実現する処理例を図5および図6のフローチャートを参照して説明する。   Next, a processing example for realizing the image diagnostic method according to the present invention will be described with reference to the flowcharts of FIGS.

図5は、平面方程式における係数a,b,cを算出するための処理例を示すフローチャートであって、この算出処理プログラムは音速像を生成するための生成処理プログラムの一部としてプログラムされている場合と、単独にサブルーチンとしてプログラムされている場合がある。   FIG. 5 is a flowchart showing a processing example for calculating the coefficients a, b, and c in the plane equation, and this calculation processing program is programmed as a part of the generation processing program for generating the sound velocity image. Sometimes it is programmed as a subroutine alone.

以下はサブルーチンとしてプログラムしたときの処理例であり、この算出プログラムが起動されると、まずXYテーブル16を走査してプレパラート12上における任意の測定ポイント(測定点)Pi(i=0,1,2であって、Po,P1,P2)の座標データ(xi,yi){i=0,1,2であって、(x0,y0)、(x1,y1)、(x2,y2)}を取得する(ステップ71)。   The following is an example of processing when programmed as a subroutine. When this calculation program is started, first, the XY table 16 is scanned and an arbitrary measurement point (measurement point) Pi (i = 0, 1, 1) on the slide 12 is scanned. 2 and Po, P1, P2) coordinate data (xi, yi) {i = 0, 1, 2 (x0, y0), (x1, y1), (x2, y2)}. Obtain (step 71).

次に、これら測定ポイントPiでの反射波信号Sri(i=0,1,2であって、Sro,Sr1,Sr2)を取得すると共に(ステップ72)、得られた反射波信号Sriをそれぞれ高速フーリエ変換してフーリエ変換出力Fiを求める(ステップ73)。   Next, the reflected wave signals Sri (i = 0, 1, 2 and Sro, Sr1, Sr2) at these measurement points Pi are acquired (step 72), and the obtained reflected wave signals Sri are respectively transmitted at high speed. Fourier transform is performed to obtain a Fourier transform output Fi (step 73).

続いて、これらフーリエ変換出力Fiより参照フーリエ変換出力Froとの比(Fr1/FroおよびFr2/Fro)を演算して規格化スペクトルを求める(ステップ74)。   Subsequently, a normalized spectrum is obtained by calculating a ratio (Fr1 / Fro and Fr2 / Fro) with the reference Fourier transform output Fro from these Fourier transform outputs Fi (step 74).

このようにして求めた規格化スペクトルの周波数成分fと位相成分φから、
φ/2πfでの時間差Δti(i=1,2)を算出する(ステップ75)。
From the frequency component f and the phase component φ of the normalized spectrum thus obtained,
A time difference Δti (i = 1, 2) at φ / 2πf is calculated (step 75).

測定ポイントPi(i=1,2)における参照点Poからの時間差Δt1,Δt2と既知の音速Co(この例では水の音速)から、測定ポイントP1,P2での平面方程式Zi(i=1,2であって、Z1,Z2)を求める(ステップ76)。この平面方程式Ziが参照点Poに対する測定ポイントP1及びP2でのz値(傾き量)となる。   From the time differences Δt1, Δt2 from the reference point Po at the measurement point Pi (i = 1, 2) and the known sound velocity Co (the sound velocity of water in this example), the plane equation Zi (i = 1, 1) at the measurement points P1, P2. 2 and Z1 and Z2) are obtained (step 76). This plane equation Zi becomes the z value (inclination amount) at the measurement points P1 and P2 with respect to the reference point Po.

このような処理を経て算出した3つの平面方程式Zo,Z1,Z2より係数a,b,cが算出され、算出された係数a,b,cがメモリ手段54に保存される(ステップ77)。したがって、これら生体組織切片13の全ての測定点Pp(xp,yp)における平面方程式Zp、つまり参照点Poを基準にした測定ポイントPpでのz値Zpを算出できることになる。なお、この係数算出処理は音速像生成の前処理として被測定物である生体組織を観察する都度実行される。   Coefficients a, b and c are calculated from the three plane equations Zo, Z1 and Z2 calculated through such processing, and the calculated coefficients a, b and c are stored in the memory means 54 (step 77). Therefore, the plane equation Zp at all the measurement points Pp (xp, yp) of these biological tissue sections 13, that is, the z value Zp at the measurement point Pp based on the reference point Po can be calculated. This coefficient calculation process is executed every time a biological tissue that is an object to be measured is observed as a pre-process for generating a sonic image.

図6は傾き補正を含めた音速算出例を示すフローチャートであって、係数a,b,cはすでに保存されているものとする。
音速像生成処理プログラムが起動されると、CPU55からの指示に基づいてX−Yステージ16が駆動されてx方向への走査が開始されると共に、エンコーダ34からのエンコーダ出力に同期して駆動制御信号がコントローラ30Xに供給される。
FIG. 6 is a flowchart showing an example of sound speed calculation including inclination correction, and it is assumed that the coefficients a, b, and c are already stored.
When the sonic image generation processing program is started, the XY stage 16 is driven based on an instruction from the CPU 55 to start scanning in the x direction, and drive control is performed in synchronization with the encoder output from the encoder 34. A signal is supplied to the controller 30X.

この駆動制御信号に同期してトリガー信号Stがパルサー部36に供給され、超音波振動子11が励起される。この一連の処理によって得られたx方向における走査線上の測定点(ピクセル)Piの座標(xi,yi)(i=1〜300)を取得すると共に(ステップ81)、測定点Piでの反射波信号Soiを取得する(ステップ82)。反射波信号Soiとしては、ほぼ同一測定点Piからの得られた複数個、この例では8個の一次反射波信号Su1を平均した信号が利用される。   In synchronization with this drive control signal, a trigger signal St is supplied to the pulsar unit 36, and the ultrasonic transducer 11 is excited. The coordinates (xi, yi) (i = 1 to 300) of the measurement point (pixel) Pi on the scanning line in the x direction obtained by this series of processing are acquired (step 81) and the reflected wave at the measurement point Pi. The signal Soi is acquired (step 82). As the reflected wave signal Soi, a signal obtained by averaging a plurality of, in this example, eight primary reflected wave signals Su1 obtained from substantially the same measurement point Pi is used.

続いて、この反射波信号Soiより測定点Piでの平面方程式、換言すれば測定点Piでのz値(傾き量)Zpが算出されて、これが保存される(ステップ83)。   Subsequently, the plane equation at the measurement point Pi, in other words, the z value (inclination amount) Zp at the measurement point Pi is calculated from the reflected wave signal Soi, and stored (step 83).

その後、反射波信号Soiと、参照すべき反射波信号Srとの時間差、反射波信号Soiの強度、位相差などから測定点Piでの生体組織の厚みdiが算出される(ステップ84)。この厚み算出過程でz値Zpが参照される。つまり、このステップ84でプレパラート12の傾きを含めた状態での厚みdiが得られ、傾きによる影響を除去した厚みdiとなる。   Thereafter, the thickness di of the living tissue at the measurement point Pi is calculated from the time difference between the reflected wave signal Soi and the reflected wave signal Sr to be referred to, the intensity of the reflected wave signal Soi, the phase difference, and the like (step 84). The z value Zp is referred to in this thickness calculation process. That is, in step 84, the thickness di including the inclination of the preparation 12 is obtained, and the thickness di is obtained by removing the influence of the inclination.

次に、算出された厚みdiから測定点Piでの音速Ciが算出され、その値が描画プレーン(フレーム)として機能するメモリ手段54に保存される(ステップ85)。そしてこの音速Ciが等音速線のデータとして画面上に描画されると共に、描画情報としてさらにメモリ手段54の別のエリアに保存される(ステップ86)。   Next, the sound velocity Ci at the measurement point Pi is calculated from the calculated thickness di, and the value is stored in the memory means 54 functioning as a drawing plane (frame) (step 85). The sound speed Ci is drawn on the screen as isosonic line data, and is further saved as drawing information in another area of the memory means 54 (step 86).

このような音速算出処理が全てのピクセル(例えば300×300ピクセル)に対して実行されるので(ステップ87)、全てのピクセルに対し音速データを取得することで、生体組織の等音速線を描画できる。等音速線をカラー表示して区別することで、カラー音速線が得られる。等音速線の描画によって生体組織の音速像が得られる。また各測定点Piでの反射波信号Soiより反射強度(減衰量)を得ることで、生体組織の減衰像を得ることができる。   Since such a sound speed calculation process is executed for all pixels (for example, 300 × 300 pixels) (step 87), by obtaining sound speed data for all pixels, an isosonic line of biological tissue is drawn. it can. A color sound speed line can be obtained by distinguishing the isosonic speed line by color display. A sonic image of a living tissue can be obtained by drawing an isosonic line. Further, an attenuation image of the biological tissue can be obtained by obtaining the reflection intensity (attenuation amount) from the reflected wave signal Soi at each measurement point Pi.

上述した実施例では、平面方程式を得るのにプレパラート12に対して任意の3点を指定して行ったが、それ以上の点を用いて平面方程式を求め、これより平面方程式の係数を算出するようにしてもよい。多点(4点以上)であればあるほど、算出した係数値の精度が増す。   In the above-described embodiment, the plane equation is obtained by designating any three points with respect to the preparation 12, but the plane equation is obtained using more points and the coefficient of the plane equation is calculated from this. You may do it. The more points (four or more), the more accurate the calculated coefficient value.

さらに、プレパラート12や生体組織13自体の凹凸、反りなどを考慮して平面方程式を算出すれば、より一層正確な生体組織切片13の厚みdおよび音速を算出できるようになる。その他は実施例1と同じである。   Furthermore, if the plane equation is calculated in consideration of the unevenness and warpage of the preparation 12 and the living tissue 13 itself, the thickness d and sound speed of the living tissue section 13 can be calculated more accurately. Others are the same as in the first embodiment.

上述した実施例では、被測定物として生体組織切片を用いたときの音速測定を電気的に(信号的に)補正して求めると共に、算出された音速から生体組織の音速像を得る場合に適用したが、被測定物の音速測定及び音速像の生成としては生体組織に限らず、種々の被測定物が測定対象となり得ることは明らかである。その他は実施例1と同じである。   In the above-described embodiment, the sound velocity measurement when using a biological tissue section as the object to be measured is obtained by correcting electrically (in a signal manner), and the sound velocity image of the biological tissue is obtained from the calculated sound velocity. However, it is obvious that the measurement of the sound velocity of the object to be measured and the generation of the sound velocity image are not limited to the living tissue, and various objects to be measured can be the measurement object. Others are the same as in the first embodiment.

この発明は、超音波顕微鏡を使用して試料の音速を測定する音速測定装置、生体組織、その他の被測定物の表面状態を観察したり診断したりする画像診断装置に適用できる。   The present invention can be applied to a sound velocity measuring device that measures the sound velocity of a sample using an ultrasonic microscope, an image diagnostic device that observes and diagnoses the surface state of a biological tissue and other objects to be measured.

この発明に係る超音波顕微鏡を使用した画像診断装置の実施例を示す要部の系統図である。It is a systematic diagram of the principal part which shows the Example of the image diagnostic apparatus using the ultrasonic microscope which concerns on this invention. 第2の信号処理部の実施例を示す要部の系統図である。It is a systematic diagram of the principal part which shows the Example of a 2nd signal processing part. 超音波振動子の励起と一次反射波信号との関係を示す波形図である。It is a wave form diagram which shows the relationship between excitation of an ultrasonic transducer | vibrator and a primary reflected wave signal. 集中パルス励起を示す波形図である。It is a wave form diagram which shows concentrated pulse excitation. この発明に係る画像診断方法を実現する場合に必要な平面方程式の係数算出例を示すフローチャートである。It is a flowchart which shows the coefficient calculation example of a plane equation required when implement | achieving the image diagnostic method which concerns on this invention. プレパラートの傾き補正を信号処理によって処理する場合の処理例を示すフローチャートである。It is a flowchart which shows the process example in the case of processing the inclination correction of a preparation by signal processing. 超音波顕微鏡の概念図である。It is a conceptual diagram of an ultrasonic microscope. 被測定物と二次元走査手段との関係を示す図である。It is a figure which shows the relationship between a to-be-measured object and a two-dimensional scanning means. 連続走査例を示す図である。It is a figure which shows the example of continuous scanning. 超音波の反射波の説明図である。It is explanatory drawing of the reflected wave of an ultrasonic wave. そのときの励起パルスと反射波との関係を示す図である。It is a figure which shows the relationship between the excitation pulse and reflected wave at that time. 超音波の帯域を示す図である。It is a figure which shows the zone | band of an ultrasonic wave. 走査領域を示す図である。It is a figure which shows a scanning area | region. 超音波によって得られた画像例を示す図である。It is a figure which shows the example of an image obtained by the ultrasonic wave. プレパラートの傾きと生体組織との関係を示す図である。It is a figure which shows the relationship between the inclination of a preparation, and a biological tissue. 従来の傾き補正を機械的に行うときの説明図である。It is explanatory drawing when performing the conventional inclination correction mechanically. プレパラートの傾きを信号処理によって補正するときの説明図である。It is explanatory drawing when correcting the inclination of a preparation by signal processing. 規格化スペクトルの曲線図である。It is a curve figure of a normalized spectrum.

符号の説明Explanation of symbols

1・・・超音波顕微鏡
11・・・超音波振動子
12・・・プレパラート
13・・・被測定物(生体組織切片)
16・・・X−Yステージ
16X・・・Xステージ
16Y・・・Yステージ
17・・・媒質(水)
21・・・音速像
1A・・・顕微鏡本体部
1B・・・第1の信号処理部
1C・・・第2の信号処理部
36・・・パルサー部
40・・・ゲート回路
42・・・ゲートパルス生成回路
45・・・平均化回路
51・・・FFT処理部
55・・・CPU
57・・・音速像生成処理部
DESCRIPTION OF SYMBOLS 1 ... Ultrasonic microscope 11 ... Ultrasonic vibrator 12 ... Preparation 13 ... Object to be measured (living tissue section)
16 ... X-Y stage 16X ... X stage 16Y ... Y stage 17 ... medium (water)
21 ... sound velocity image 1A ... microscope main body 1B ... first signal processing unit 1C ... second signal processing unit 36 ... pulsar unit 40 ... gate circuit 42 ... gate Pulse generation circuit 45 ... averaging circuit 51 ... FFT processing unit 55 ... CPU
57 ... Sound velocity image generation processing unit

Claims (22)

超音波振動子をパルス励起することで被測定物の音速を測定するパルス励起型超音波顕微鏡を使用した音速測定方法であって、
所定の間隔を有するパルス励起用のトリガー信号を生成するトリガー信号生成ステップと、
上記トリガー信号が供給されて、上記トリガー信号よりも間隔の狭い励起パルスで上記超音波振動子を集中的に励起する集中パルス励起ステップと、
この集中パルス励起期間中に得られる上記励起パルスからの反射波信号を平均化するステップと、
上記超音波振動子と対峙する試料載置板の傾き量を検出するステップと、
検出されたこの傾き量と、平均化した上記反射波信号とに基づいて上記被測定物の厚みを算出するステップと、
上記傾き量によって電気的に補正された上記被測定物の厚みに基づいて上記被測定物の音速を算出するステップとを
有することを特徴とする超音波顕微鏡を使用した音速測定方法。
A sound velocity measurement method using a pulse excitation type ultrasonic microscope that measures the sound velocity of an object to be measured by pulse exciting an ultrasonic transducer,
A trigger signal generating step for generating a trigger signal for pulse excitation having a predetermined interval ;
And the trigger signal is supplied, and concentrated pulsed excitation step of intensively excite the ultrasonic transducer with narrow I励 cause pulse intervals than the trigger signal,
A step of averaging the reflected wave signal from Ki励 force pulse on obtained during the intensive pulsed excitation period,
Detecting the amount of inclination of the sample mounting plate facing the ultrasonic transducer;
Calculating the thickness of the object to be measured based on the detected amount of inclination and the averaged reflected wave signal;
And a step of calculating a sound speed of the object to be measured based on a thickness of the object to be electrically corrected by the tilt amount. A method of measuring a sound speed using an ultrasonic microscope.
上記傾き量を検出するステップは、上記試料載置面の平面方程式の係数を算出するステップと、
上記被測定物に対する任意の被測定点における傾き量が、上記平面方程式の係数に基づいて算出されるステップとを含む
ことを特徴とする請求項1記載の超音波顕微鏡を使用した音速測定方法。
The step of detecting the amount of inclination includes calculating a coefficient of a plane equation of the sample mounting surface;
The method for measuring the speed of sound using an ultrasonic microscope according to claim 1, further comprising a step of calculating an inclination amount at an arbitrary measurement point with respect to the measurement object based on a coefficient of the plane equation.
上記平面方程式の係数を算出するステップは、上記試料載置板の少なくとも3つの測定点からの反射波信号を用いて行う
ことを特徴とする請求項1記載の超音波顕微鏡を使用した音速測定方法。
2. The sound speed measuring method using an ultrasonic microscope according to claim 1, wherein the step of calculating the coefficient of the plane equation is performed using reflected wave signals from at least three measurement points of the sample mounting plate. .
上記測定点は、上記試料載置板の平面であって、上記被測定物の非載置面が使用される
ことを特徴とする請求項3記載の超音波顕微鏡を使用した音速測定方法。
4. The sound speed measuring method using an ultrasonic microscope according to claim 3, wherein the measurement point is a plane of the sample mounting plate and a non-mounting surface of the object to be measured is used.
記励起パルスは、広帯域のパルス信号である
ことを特徴とする請求項1記載の超音波顕微鏡を使用した音速測定方法。
Upper Ki励 causing pulses, speed-of-sound measurement method using the ultrasonic microscope according to claim 1, characterized in that the wide-band pulse signals.
超音波振動子をパルス励起することで被測定物の音速を測定するパルス励起型超音波顕微鏡を使用した音速測定装置あって、
所定の間隔を有するパルス励起用のトリガー信号を生成するトリガー信号生成部と、
上記トリガー信号生成部で生成された上記トリガー信号が供給されて、上記トリガー信号よりも間隔の狭い励起パルスで上記超音波振動子を集中的に励起する集中パルス励起機能を有したパルサー部と、
集中パルス励起期間中に得られる上記パルサー部で励起された上記励起パルスからの反射波信号を平均化する信号演算手段と、
上記超音波振動子と対峙する試料載置板の傾き量を検出する傾き量検出手段と、
検出されたこの傾き量と、上記信号演算手段で平均化された上記反射波信号とに基づいて上記被測定物の厚みを算出する厚み算出手段と、
上記傾き量によって電気的に補正された上記被測定物の厚みに基づいて、上記被測定物の音速を測定する音速測定手段と
で構成されたことを特徴とする超音波顕微鏡を使用した音速測定装置。
There is a sound velocity measuring device using a pulse excitation type ultrasonic microscope that measures the sound velocity of an object to be measured by pulse exciting an ultrasonic transducer,
A trigger signal generator for generating a trigger signal for pulse excitation having a predetermined interval ;
The trigger signal the trigger signal generated by the generator is supplied, pulser unit having a centralized pulsed excitation function of centrally excite the ultrasonic transducer with narrow I励 cause pulse intervals than the triggering signal When,
And signal calculation means for averaging the reflected wave signal from Ki励 electromotive pulse after being excited by the pulser unit obtained during intensive pulsed excitation period,
A tilt amount detecting means for detecting a tilt amount of the sample mounting plate facing the ultrasonic transducer;
Thickness calculating means for calculating the thickness of the object to be measured based on the detected amount of inclination and the reflected wave signal averaged by the signal calculating means;
Sound velocity measurement using an ultrasonic microscope characterized by comprising a sound velocity measuring means for measuring the sound velocity of the object to be measured based on the thickness of the object to be measured electrically corrected by the tilt amount. apparatus.
上記傾き量検出は、上記試料載置板の平面方程式を算出する算出手段を有し、
上記被測定物に対する任意の被測定点における傾き量が、上記平面方程式の係数に基づいて算出される
ことを特徴とする請求項6記載の超音波顕微鏡を使用した音速測定装置。
The inclination amount detection has a calculation means for calculating a plane equation of the sample mounting plate,
The sound velocity measuring apparatus using an ultrasonic microscope according to claim 6, wherein an inclination amount at an arbitrary measured point with respect to the measured object is calculated based on a coefficient of the plane equation.
上記平面方程式の算出手段では、上記試料載置板の少なくとも3つの測定点からの反射波信号を用いて上記平面方程式の係数が算出される
ことを特徴とする請求項6記載の超音波顕微鏡を使用した音速測定装置。
The ultrasonic microscope according to claim 6, wherein the plane equation calculation means calculates the coefficient of the plane equation using reflected wave signals from at least three measurement points of the sample mounting plate. Used sound speed measurement device.
上記測定点は、上記試料載置板の平面であって、上記被測定物の非載置面が使用される
ことを特徴とする請求項8記載の超音波顕微鏡を使用した音速測定装置。
9. The sound velocity measuring apparatus using an ultrasonic microscope according to claim 8, wherein the measurement point is a flat surface of the sample mounting plate and a non-mounting surface of the object to be measured is used.
記励起パルスは、広帯域のパルス信号である
ことを特徴とする請求項6記載の超音波顕微鏡を使用した音速測定装置。
Upper Ki励 electromotive pulse, acoustic velocity measuring apparatus using an acoustic microscope according to claim 6, characterized in that the wide-band pulse signals.
超音波振動子をパルス励起することで被測定物の性状を音速像として取得するパルス励起型超音波顕微鏡を使用した音速像取得方法であって、
所定の間隔を有するパルス励起用のトリガー信号を生成するトリガー信号生成ステップと、
上記トリガー信号が供給されて、上記トリガー信号よりも間隔の狭い励起パルスで上記超音波振動子を集中的に励起する集中パルス励起ステップと、
この集中パルス励起期間中に得られる上記励起パルスからの反射波信号を平均化するステップと、
上記超音波振動子と対峙する試料載置板の傾き量を検出するステップと、
検出されたこの傾き量と、平均化した上記反射波信号とに基づいて上記被測定物の厚みを算出するステップと、
上記反射波信号に基づいて算出された上記被測定物の厚みを上記傾き量によって電気的に補正するステップと、
上記被測定物の厚み情報から当該被測定物の音速を算出するステップと、
この音速から上記被測定物の音速像を生成するステップとを有する
ことを特徴とする超音波顕微鏡を使用した音速像取得方法。
A sound velocity image acquisition method using a pulse excitation type ultrasonic microscope that acquires a property of an object to be measured as a sound velocity image by pulse exciting an ultrasonic transducer,
A trigger signal generating step for generating a trigger signal for pulse excitation having a predetermined interval ;
And the trigger signal is supplied, and concentrated pulsed excitation step of intensively excite the ultrasonic transducer with narrow I励 cause pulse intervals than the trigger signal,
A step of averaging the reflected wave signal from Ki励 force pulse on obtained during the intensive pulsed excitation period,
Detecting the amount of inclination of the sample mounting plate facing the ultrasonic transducer;
Calculating the thickness of the object to be measured based on the detected amount of inclination and the averaged reflected wave signal;
Electrically correcting the thickness of the measurement object calculated based on the reflected wave signal by the amount of inclination;
Calculating the sound velocity of the device under test from the thickness information of the device under test;
Generating a sound speed image of the object to be measured from the sound speed, and a sound speed image acquiring method using an ultrasonic microscope.
上記傾き量を検出するステップは、上記試料載置面の平面方程式の係数を算出するステップと、
上記被測定物に対する任意の被測定点における傾き量が、上記平面方程式の係数に基づいて算出されるステップとを含む
ことを特徴とする請求項11記載の超音波顕微鏡を使用した音速像取得方法。
The step of detecting the amount of inclination includes calculating a coefficient of a plane equation of the sample mounting surface;
The method according to claim 11, further comprising: calculating an inclination amount at an arbitrary measurement point with respect to the measurement object based on a coefficient of the plane equation. .
上記平面方程式の係数を算出するステップは、上記試料載置板の少なくとも3つの測定点からの反射波信号を用いて行う
ことを特徴とする請求項11記載の超音波顕微鏡を使用した音速像取得方法。
The sound velocity image acquisition using an ultrasonic microscope according to claim 11, wherein the step of calculating the coefficient of the plane equation is performed using reflected wave signals from at least three measurement points of the sample mounting plate. Method.
上記測定点は、上記試料載置板の平面であって、上記被測定物の非載置面が使用される
ことを特徴とする請求項13記載の超音波顕微鏡を使用した音速像取得方法。
The method according to claim 13, wherein the measurement point is a plane of the sample mounting plate and a non-mounting surface of the object to be measured is used.
記励起パルスは、広帯域のパルス信号である
ことを特徴とする請求項11記載の超音波顕微鏡を使用した音速像取得方法。
Upper Ki励 electromotive pulse, acoustic velocity image acquisition method using an ultrasonic microscope according to claim 11, wherein it is a wide-band pulse signals.
上記被測定物は、生体組織切片である
ことを特徴とする請求項11記載の超音波顕微鏡を使用した音速像取得方法。
The method according to claim 11, wherein the object to be measured is a biological tissue section.
超音波振動子をパルス励起することで被測定物の性状を音速像として取得するパルス励起型超音波顕微鏡を使用した画像診断装置であって、
所定の間隔を有するパルス励起用のトリガー信号を生成するトリガー信号生成部と、
上記トリガー信号生成部で生成された上記トリガー信号が供給されて、上記トリガー信号よりも間隔の狭い励起パルスで上記超音波振動子を集中的に励起する集中パルス励起機能を有したパルサー部と、
集中パルス励起期間中に得られる上記パルサー部で励起された上記励起パルスからの反射波信号を平均化する信号演算手段と、
上記超音波振動子と対峙する試料載置板の傾き量を検出する傾き量検出手段と、
検出されたこの傾き量と、上記信号演算手段で平均化された上記反射波信号とに基づいて上記被測定物の厚みを算出する厚み算出手段と、
上記反射波信号に基づいて算出された上記被測定物の厚みを上記傾き量によって電気的に補正する補正手段と、
上記被測定物の厚み情報から当該被測定物の音速を算出する音速算出手段と、
この音速から上記被測定物の音速像を生成する音速像生成手段とを有する
ことを特徴とする超音波顕微鏡を使用した画像診断装置。
An image diagnostic apparatus using a pulse excitation type ultrasonic microscope that acquires a property of an object to be measured as a sound velocity image by pulse exciting an ultrasonic transducer,
A trigger signal generator for generating a trigger signal for pulse excitation having a predetermined interval ;
The trigger signal the trigger signal generated by the generator is supplied, pulser unit having a centralized pulsed excitation function of centrally excite the ultrasonic transducer with narrow I励 cause pulse intervals than the triggering signal When,
And signal calculation means for averaging the reflected wave signal from Ki励 electromotive pulse after being excited by the pulser unit obtained during intensive pulsed excitation period,
A tilt amount detecting means for detecting a tilt amount of the sample mounting plate facing the ultrasonic transducer;
Thickness calculating means for calculating the thickness of the object to be measured based on the detected amount of inclination and the reflected wave signal averaged by the signal calculating means;
Correction means for electrically correcting the thickness of the measurement object calculated based on the reflected wave signal by the amount of inclination;
A sound velocity calculating means for calculating the sound velocity of the object to be measured from the thickness information of the object to be measured;
An image diagnostic apparatus using an ultrasonic microscope, comprising: a sound velocity image generating means for generating a sound velocity image of the object to be measured from the sound velocity.
上記傾き量検出は、上記試料載置板の平面方程式を算出する算出手段を有し、
上記被測定物に対する任意の被測定点における傾き量が、上記平面方程式の係数に基づいて算出される
ことを特徴とする請求項17記載の超音波顕微鏡を使用した画像診断装置。
The inclination amount detection has a calculation means for calculating a plane equation of the sample mounting plate,
18. The diagnostic imaging apparatus using an ultrasonic microscope according to claim 17, wherein an inclination amount at an arbitrary measurement point with respect to the measurement object is calculated based on a coefficient of the plane equation.
上記平面方程式の算出手段では、上記試料載置板の少なくとも3つの測定点からの反射波信号を用いて上記平面方程式の係数が算出される
ことを特徴とする請求項17記載の超音波顕微鏡を使用した画像診断装置。
The ultrasonic microscope according to claim 17, wherein the plane equation calculating means calculates the coefficient of the plane equation using reflected wave signals from at least three measurement points of the sample mounting plate. Used diagnostic imaging equipment.
上記測定点は、上記試料載置板の平面であって、上記被測定物の非載置面が使用される
ことを特徴とする請求項19記載の超音波顕微鏡を使用した画像診断装置。
20. The diagnostic imaging apparatus using an ultrasonic microscope according to claim 19, wherein the measurement point is a plane of the sample mounting plate and a non-mounting surface of the object to be measured is used.
記励起パルスは、広帯域のパルス信号である
ことを特徴とする請求項17記載の超音波顕微鏡を使用した画像診断装置。
Upper Ki励 causing pulses, image diagnosis apparatus using an ultrasonic microscope according to claim 17, wherein it is a wide-band pulse signals.
上記被測定物は、生体組織切片である
ことを特徴とする請求項17記載の超音波顕微鏡を使用した画像診断装置。
18. The diagnostic imaging apparatus using an ultrasonic microscope according to claim 17, wherein the object to be measured is a biological tissue section.
JP2004105213A 2004-03-31 2004-03-31 Sonic velocity measuring method using ultrasonic microscope, sonic velocity measuring device, sonic velocity image acquiring method using ultrasonic microscope, and diagnostic imaging device thereof Expired - Fee Related JP4521584B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2004105213A JP4521584B2 (en) 2004-03-31 2004-03-31 Sonic velocity measuring method using ultrasonic microscope, sonic velocity measuring device, sonic velocity image acquiring method using ultrasonic microscope, and diagnostic imaging device thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2004105213A JP4521584B2 (en) 2004-03-31 2004-03-31 Sonic velocity measuring method using ultrasonic microscope, sonic velocity measuring device, sonic velocity image acquiring method using ultrasonic microscope, and diagnostic imaging device thereof

Publications (2)

Publication Number Publication Date
JP2005291827A JP2005291827A (en) 2005-10-20
JP4521584B2 true JP4521584B2 (en) 2010-08-11

Family

ID=35324928

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004105213A Expired - Fee Related JP4521584B2 (en) 2004-03-31 2004-03-31 Sonic velocity measuring method using ultrasonic microscope, sonic velocity measuring device, sonic velocity image acquiring method using ultrasonic microscope, and diagnostic imaging device thereof

Country Status (1)

Country Link
JP (1) JP4521584B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007171051A (en) * 2005-12-22 2007-07-05 Honda Electronic Co Ltd Sonic velocity measuring method and device, and ultrasonic image inspecting device
JP4787914B2 (en) * 2005-12-22 2011-10-05 本多電子株式会社 Sonic velocity measuring method, sonic velocity measuring apparatus, and ultrasonic image inspection apparatus
JP5017572B2 (en) * 2006-07-06 2012-09-05 本多電子株式会社 Ultrasonic image inspection method and ultrasonic image inspection apparatus
JP5473683B2 (en) * 2010-03-09 2014-04-16 三菱電機株式会社 Feature detection system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2624749B2 (en) * 1988-03-04 1997-06-25 日立建機株式会社 Ultrasonic flaw detector
JP2002014083A (en) * 2000-04-28 2002-01-18 Nkk Corp Method and device for judging solidification state of cast piece, and method for manufacturing continuously casting cast piece
JP2003050117A (en) * 2001-08-08 2003-02-21 Olympus Optical Co Ltd Method and apparatus for measuring thickness

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03262909A (en) * 1990-03-13 1991-11-22 Olympus Optical Co Ltd Ultrasonic measuring instrument

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2624749B2 (en) * 1988-03-04 1997-06-25 日立建機株式会社 Ultrasonic flaw detector
JP2002014083A (en) * 2000-04-28 2002-01-18 Nkk Corp Method and device for judging solidification state of cast piece, and method for manufacturing continuously casting cast piece
JP2003050117A (en) * 2001-08-08 2003-02-21 Olympus Optical Co Ltd Method and apparatus for measuring thickness

Also Published As

Publication number Publication date
JP2005291827A (en) 2005-10-20

Similar Documents

Publication Publication Date Title
US6336899B1 (en) Ultrasonic diagnosis apparatus
EP2881041A1 (en) Apparatus and method for ultrasonic diagnosis
US6464642B1 (en) Ultrasonic diagnosis apparatus
JP2006280532A (en) Ultrasonic diagnostic equipment
JP5917037B2 (en) Subject information acquisition apparatus and subject information acquisition method
KR20140020486A (en) Method and apparatus for analyzing elastography of tissue using ultrasound
CN110151270B (en) Vibration puncture device, frequency adjustment method and ultrasonic imaging system
WO2013153968A1 (en) Ultrasonic diagnosis device
JP2021514769A (en) Hybrid elastography methods, probes and devices for hybrid elastography
JP7052530B2 (en) Ultrasonic diagnostic equipment and ultrasonic signal processing method
JP2005319199A (en) Ultrasonic diagnostic apparatus
JP4521584B2 (en) Sonic velocity measuring method using ultrasonic microscope, sonic velocity measuring device, sonic velocity image acquiring method using ultrasonic microscope, and diagnostic imaging device thereof
JP4654335B2 (en) Ultrasonic image inspection method, ultrasonic image inspection apparatus
JP4830100B2 (en) Inspected object measuring method and inspected object measuring apparatus
JP2020178905A (en) Ultrasonic signal processing device, ultrasonic diagnostic device, ultrasonic signal processing method, and ultrasonic signal processing program
JP4521585B2 (en) Sonic velocity measuring method using ultrasonic microscope, sonic velocity measuring device, sonic velocity image acquiring method, and diagnostic imaging device
JP2015519120A (en) Method for imaging specular object and target anatomy in tissue using ultrasound and ultrasound imaging apparatus
JP4996141B2 (en) Ultrasonic diagnostic equipment
JP5317395B2 (en) Ultrasonic diagnostic apparatus and ultrasonic diagnostic image display method
JP4613269B2 (en) Sonic velocity measuring method and sonic velocity measuring device
CN210447159U (en) Vibration puncture device and ultrasonic imaging system
JP2008046097A (en) Acoustic velocity measuring method and acoustic velocity measuring instrument
JP4830099B2 (en) Sonic velocity measuring method and sonic velocity measuring device
JP4787914B2 (en) Sonic velocity measuring method, sonic velocity measuring apparatus, and ultrasonic image inspection apparatus
JP5124741B2 (en) Acoustic impedance measuring method and acoustic impedance measuring device

Legal Events

Date Code Title Description
RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20060607

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070220

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20070220

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20091022

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20091104

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20091228

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100126

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100329

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: 20100413

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20100507

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

Ref document number: 4521584

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130604

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130604

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313117

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees