JP6365243B2 - 超音波測定装置、超音波画像装置及び超音波測定方法 - Google Patents

超音波測定装置、超音波画像装置及び超音波測定方法 Download PDF

Info

Publication number
JP6365243B2
JP6365243B2 JP2014222474A JP2014222474A JP6365243B2 JP 6365243 B2 JP6365243 B2 JP 6365243B2 JP 2014222474 A JP2014222474 A JP 2014222474A JP 2014222474 A JP2014222474 A JP 2014222474A JP 6365243 B2 JP6365243 B2 JP 6365243B2
Authority
JP
Japan
Prior art keywords
wave
ultrasonic
base
signal
waves
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2014222474A
Other languages
English (en)
Other versions
JP2016086958A (ja
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.)
Seiko Epson Corp
Original Assignee
Seiko Epson Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Seiko Epson Corp filed Critical Seiko Epson Corp
Priority to JP2014222474A priority Critical patent/JP6365243B2/ja
Priority to US14/928,076 priority patent/US10849598B2/en
Publication of JP2016086958A publication Critical patent/JP2016086958A/ja
Application granted granted Critical
Publication of JP6365243B2 publication Critical patent/JP6365243B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8959Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes
    • G01S15/8963Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes using pulse inversion
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8977Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52038Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4405Device being mounted on a trolley
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4427Device being portable or laptop-like
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Nonlinear Science (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Theoretical Computer Science (AREA)

Description

本発明は、超音波測定装置、超音波画像装置及び超音波測定方法等に関係する。
被検体である人体の内部を検査するために用いる装置として、対象物に向けて超音波を出射し、対象物内部における音響インピーダンスの異なる界面からの反射波を受信する超音波測定装置が注目されている。さらに、超音波測定装置は、内臓脂肪の測定や血流量の測定など、被検体の表層の画像診断にも応用されている。
このような超音波測定装置を用いて画像診断を行う場合には、超音波エコーの画像処理の高分解能化を図る必要があり、例えばハーモニックイメージング(ハーモニックイメージング法)などが利用されている。
ここで、ハーモニックイメージングにおいては、超音波エコーのハーモニック成分を抽出する必要があるが、そのためのハーモニック成分抽出方法としては、フィルター法や、特許文献1等に記載される位相反転法がある。特許文献1では、3次以上の高次の高調波を用いて位相反転法を行う超音波画像装置が開示されている。
特開2002−360569号公報
従来のBモード画像生成処理では、距離分解能と方位分解能がほぼ同じであったため、距離分解能の向上は課題とはならなかった。しかし、ハーモニックイメージングや適応型ビームフォーミングを適応すると、距離分解能よりも方位分解能が高くなり、生成された画像における分解能の異方性が新たに生じることとなった。例えば、前述した特許文献1でも、方位分解能に対して距離分解能が相対的に低下している。そのため、距離分解能の向上が必要である。
本発明の幾つかの態様によれば、超音波による対象物の測定結果の方位分解能だけでなく、距離分解能も向上させることができる超音波測定装置、超音波画像装置及び超音波測定方法等を提供することができる。
本発明の一態様は、対象物に対して超音波を送信する処理を行う送信処理部と、送信した前記超音波に対する超音波エコーの受信処理を行う受信処理部と、前記受信処理部からの受信信号に対して処理を行う処理部と、を含み、前記処理部は、前記送信処理部により送信された送信パルス信号に対応する受信信号に対して、前記受信信号を構成する複数の第1基底波の結合係数特定処理を行い、前記結合係数特定処理により特定された複数の結合係数と、前記第1基底波よりも波数が少ない第2基底波とに基づいて、前記受信信号から再構成信号への変換処理を行う超音波測定装置に関係する。
本発明の一態様では、送信された送信パルス信号に対応する受信信号に対して、受信信号を構成する複数の第1基底波の結合係数特定処理を行い、特定された複数の結合係数と、第1基底波よりも波数が少ない第2基底波とに基づいて、受信信号から再構成信号への変換処理を行う。よって、超音波による対象物の測定結果の方位分解能だけでなく、距離分解能も向上させることが可能となる。
また、本発明の一態様では、前記処理部は、前記受信信号に対応する高調波に対し、前記第1基底波の前記結合係数特定処理を行った後に、前記変換処理を行って、前記再構成信号として、前記第2基底波による再構成波を生成してもよい。
これにより、第1基底波よりも波数が少ない第2基底波から構成される再構成波に、受信波を変換し、対象物の測定結果における距離分解能を向上させること等が可能になる。
また、本発明の一態様では、前記複数の第1基底波は、M個(Mは2以上の整数)の第1基底波であり、前記M個の第1基底波のうちの第i(iは、1≦i≦Mの整数)の第1基底波は、第iの測定点に配置された第iの点散乱体からの前記超音波の受信信号に対応する高調波であり、前記M個の第1基底波のうちの第(i+1)の第1基底波は、前記超音波の送信点から、前記第iの測定点よりも遠い位置である第(i+1)の測定点に配置された第(i+1)の点散乱体からの超音波の受信信号に対応する高調波であってもよい。
これにより、設定した各測定点の間隔に対応する距離分解能で、受信波から第1基底波成分を抽出すること等が可能になる。
また、本発明の一態様では、前記第1基底波は、前記受信信号から抽出可能な高調波であってもよい。
これにより、受信信号を複数の第1基底波に分解すること等が可能になる。
また、本発明の一態様では、前記送信処理部は、互いに位相を反転させた2つのパルス信号を前記対象物に送信し、前記処理部は、送信された前記2つのパルス信号に対応する2つの受信信号に基づいて減算処理を行って、1つの差分信号を求め、求めた前記差分信号に対して第1フィルター処理を行って、高調波成分を抽出し、抽出された前記高調波成分に基づいて、所与の測定点に配置された点散乱体からの反射波成分に対応する高調波を、前記第1基底波として求めてもよい。
これにより、対象物内の所与の測定点の点散乱体からの反射波成分に対応する第1基底波を特定すること等が可能になる。
また、本発明の一態様では、前記複数の第1基底波は、M個(Mは2以上の整数)の第1基底波であり、前記M個の第1基底波のうちの第i(iは、1≦i≦Mの整数)の第1基底波は、前記送信パルス信号のパルス幅又は前記受信信号のパルス幅に対応する位相差よりも短い位相差で、第(i+1)の第1基底波と相互に位相がずれていてもよい。
これにより、送信パルス信号のパルス幅又は受信信号のパルス幅に対応する位相差よりも短い距離の距離分解能で、対象物を測定すること等が可能になる。
また、本発明の一態様では、前記処理部は、複数の第2基底波に基づいて、前記変換処理を行い、前記複数の第2基底波は、N個(Nは2以上の整数)の第2基底波であり、前記N個の第2基底波のうちの第j(jは、1≦j≦Nの整数)の第2基底波は、第jの測定点に配置された第jの点散乱体からの前記超音波の受信信号に対応する高調波であり、前記N個の第2基底波のうちの第(j+1)の第2基底波は、前記超音波の送信点から、前記第jの測定点よりも遠い位置である第(j+1)の測定点に配置された第(j+1)の点散乱体からの超音波の受信信号に対応する高調波であってもよい。
これにより、超音波による対象物の測定結果による距離分解能を、設定した各測定点の間隔まで向上させること等が可能になる。
また、本発明の一態様では、前記第2基底波は、前記受信信号から抽出可能な基本波に対して、時間成分の圧縮処理を行って得られてもよい。
これにより、減算処理、フィルター処理及び時間成分の圧縮処理等の簡易な処理で、第2基底波を特定すること等が可能になる。
また、本発明の一態様では、前記送信処理部は、互いに位相を反転させた2つのパルス信号を前記対象物に送信し、前記処理部は、送信された前記2つのパルス信号に対応する2つの受信信号に基づいて減算処理を行って、1つの差分信号を求め、求めた前記差分信号に対して第2フィルター処理を行って、基本波成分を抽出し、抽出された前記基本波成分に基づいて、所与の測定点に配置された点散乱体からの反射波成分に対応する基本波を求め、求めた前記基本波に対して時間成分の圧縮処理を行って、前記第2基底波を求めてもよい。
これにより、対象物内の所与の測定点の点散乱体からの反射波成分に対応する第1基底波の波長を短くした第2基底波を特定すること等が可能になる。
また、本発明の一態様では、前記第2基底波は、前記第1基底波と位相差が同じで、波数が少なくてもよい。
これにより、受信波を構成する基底波のパルス幅を短くすること等が可能になる。
また、本発明の一態様では、前記処理部は、前記変換処理後に得られた前記再構成信号に対して、包絡線検波処理を行ってもよい。
これにより、ユーザーが測定結果を判別しやすい波形データを表示部に表示すること等が可能になる。
また、本発明の一態様では、前記処理部は、前記結合係数特定処理として、前記受信信号のデコンボリューション処理を行ってもよい。
これにより、受信波を構成する第1基底波の結合係数を特定すること等が可能になる。
また、本発明の一態様では、前記処理部は、前記再構成信号の前記変換処理として、前記第2基底波のコンボリューション処理を行ってもよい。
これにより、受信波又は受信波から抽出可能な高調波を構成する第1基底波の波数を削減した第2基底波によって、再構成波を生成すること等が可能になる。
また、本発明の他の態様は、超音波測定装置と、前記再構成信号に基づいて生成された表示用画像データを表示する表示部と、を含む超音波画像装置に関係する。
また、本発明の他の態様は、対象物に対して超音波を送信する処理を行い、送信した前記超音波に対する超音波エコーの受信処理を行い、送信された送信パルス信号に対応する受信信号に対して、前記受信信号を構成する複数の第1基底波の結合係数特定処理を行い、前記結合係数特定処理により特定された複数の結合係数と、前記第1基底波よりも波数が少ない第2基底波とに基づいて、前記受信信号から再構成信号への変換処理を行う超音波測定方法に関係する。
フィルター法の説明図。 図2(A)〜図2(C)は、位相反転法の説明図。 図3(A)、図3(B)は、位相反転法とフィルター処理を併用する処理の説明図。 図4(A)、図4(B)は、距離分解能とパルス幅の関係についての説明図。 本実施形態のシステム構成例。 本実施形態の超音波画像装置の詳細なシステム構成例。 図7(A)〜図7(C)は、超音波測定装置の具体的な機器構成の一例。 図8(A)〜図8(D)は、再構成波の生成処理の説明図。 本実施形態の全体の処理の流れを説明するフローチャート。 再構成波の生成処理の流れを説明するフローチャート。 第1基底波及び第2基底波の生成処理の流れを説明するフローチャート。 図12(A)〜図12(C)は、第2基底波の生成処理の詳細な説明図。 図13(A)〜図13(C)は、点散乱体と第1基底波及び第2基底波の対応関係の説明図。 図14(A)、図14(B)は、測定結果の説明図。 図15(A)〜図15(C)は、超音波トランスデューサー素子の構成例。 超音波トランスデューサーデバイスの構成例。 図17(A)、図17(B)は、各チャンネルに対応して設けられる超音波トランスデューサー素子群の構成例。
以下、本実施形態について説明する。なお、以下に説明する本実施形態は、特許請求の範囲に記載された本発明の内容を不当に限定するものではない。また、本実施形態で説明される構成の全てが、本発明の必須構成要件であるとは限らない。
1.概要
被検体である人体の内部を検査するために用いる装置として、対象物に向けて超音波を出射し、対象物内部における音響インピーダンスの異なる界面からの反射波を受信する超音波測定装置が知られている。さらに、超音波測定装置の応用例としては、内臓脂肪の測定や血流量の測定など、被検体の表層の画像診断を行うポケット型超音波ビューアなどがあり、ヘルスケア分野への展開が期待されている。
上記のように、超音波測定装置を用いて画像診断を行う場合には、超音波エコーの画像処理の高分解能化を図る必要がある。そして、高分解能化を実現する為の画像処理技術としては、ハーモニックイメージング(ハーモニックイメージング法)がある。
ハーモニックイメージングとは、後述するハーモニック成分を映像化する手法のことをいう。ここで、媒質中を伝搬する超音波(粗密波)の速度は、音圧の高い部分は速く、低い部分では遅くなるという性質がある。したがって、単純な正弦波であっても伝搬過程で徐々に歪みが生じて波形が変化し、基本波には含まれなかった基本周波数の整数倍の高調波成分(これをハーモニック成分又は非線形成分とも言う)が含まれるようになる。このような非線形効果は、超音波の音圧の2乗に比例して大きくなり、また伝搬距離に比例して蓄積する。
そして、ハーモニックイメージングは、超音波が組織を伝搬する時に組織自身から発生するハーモニック成分を映像化するティッシュハーモニックイメージングと、超音波造形剤の微小気泡が共振、崩壊する時に発生するハーモニック成分を映像化する造影ハーモニックイメージングの二つに大別される。本実施形態では、ティッシュハーモニックイメージングを用いる。
また、ハーモニックイメージングには2つの利点がある。まず、ハーモニック成分の振幅は送信超音波の振幅の2乗に比例するという特徴があることから、ハーモニック成分の振幅は、音圧の高い送信ビーム中央では強いが、ビーム中央から端になるほど急激に弱くなる。これにより、ハーモニックイメージングでは、非線形効果の生じる範囲はビーム中央に制限され、結果的に他の手法に比べて方位分解能が向上する。これが第1の利点である。
また、超音波画像に乗る主なノイズとしては、多重反射によるノイズとサイドローブによるノイズがある。ここで、反射した超音波エコーは音圧が低く、ハーモニック成分自体が発生しない。そのため、多重反射によるノイズが低減される。さらに、サイドローブは音圧が低く、サイドローブでもハーモニック成分自体が発生しない。そのため、サイドローブによるノイズも低減される。このように、ハーモニックイメージングでは、多重反射によるノイズも、サイドローブによるノイズも低減することができる。これが第2の利点である。
本実施形態では、ハーモニックイメージングの中でも、3次高調波成分を映像化する3次ハーモニックイメージングを行う。3次ハーモニックイメージングは、2次高調波成分を映像化する手法に対して、ビーム幅が細くなるため、更に方位分解能を向上することができる。
ここで、3次ハーモニックイメージングにおいては、超音波エコーの3次高調波成分を抽出する必要があるが、そのための抽出方法としては、フィルター法と位相反転法とがある。
まず、フィルター法とは、周波数フィルター(ハイパスフィルター)により基本波成分及び2次高調波成分と、3次高調波成分とを分離し、3次高調波成分だけを抽出し、映像化する手法である。例えば、フィルター法を説明する図として、基本波帯域の中心周波数がfであり、2次高調波帯域の中心周波数が2fであり、3次高調波帯域の中心周波数が3fである受信信号を、縦軸を信号強度、横軸を周波数とした図1のグラフに示す。実際には図1に示すように、受信する基本波成分及び2次高調波成分と、3次高調波成分は、ある帯域幅を有しているため、2次高調波成分と3次高調波成分は重複し、両者を分離できなくなり、画像劣化の要因となる。この重複を少なくするためには、パルス幅を長くする必要性がある。しかし、パルス幅が長くなると距離分解能が低下する。
一方で、位相反転法は、フィルター法の欠点を改善するために開発された手法である。この手法は、同一方向に続けて2回の超音波の送信を行う。図2(A)に示すように、2回目の送信波は、1回目の送信波に対して位相が180度異なる。
そして、生体や造影剤から反射して戻ってくる受信波は、その非線形な伝播特性によりハーモニック成分を含むため、歪んだ波形となる。各回の送信波に対する受信波を、基本波、2次高調波及び3次高調波に分解して図示すると、図2(B)のようになる。図2(B)に示すように、この2回の受信波の間には、送信波を1回目と2回目で反転させているために基本波成分及び奇数次高調波成分(3次高調波成分)は反転しているが、偶数次高調波成分(2次高調波成分)は反転していないという関係がある。つまり、この2回の送信波に対する2回の受信波は、基本波成分及び奇数次高調波成分は互いに位相反転しているが、偶数次高調波成分は同相となる。
そのため、2回の受信波を減算すると、図2(C)に示すように、2次高調波成分は除去され、基本波成分と3次高調波成分は振幅が2倍になって残る。従って、基本波成分と3次高調波成分を抽出することが可能となる。
さらに、位相反転法により抽出した基本波成分、奇数次高調波成分、または、偶数次高調波成分から、目的とするN次高調波成分のみを抽出する場合は(Nは2以上の整数)、前述したフィルター法と組合せる必要がある。本実施形態では、図3(A)に示す、位相反転法により抽出した基本波成分と3次高調波成分から、周波数フィルター(ハイパスフィルター又はバンドパスフィルター)により基本波成分と3次高調波成分とを分離し、図3(B)に示すように3次高調波成分だけを抽出して映像化する。
このように、位相反転法とフィルター法を併用することにより、従来の基本波成分のみからBモード画像を生成する手法と比べて、サイドローブや多重反射によるアーチファクトが少なく、方位分解能を向上させた高品質なBモード画像を生成することができる。
しかし、上記のような方法でBモード画像を生成した場合でも、距離分解能は向上しない。距離分解能Δxは、パルス幅によって決まり、下式(1)により求められる。なお、下式(1)において、nは波数であり、λは波長である。
例えば、図4(A)に示す送信波に対する受信波では、基本波成分PS1の波長λ1に比べて、3次高調波成分PS3の波長λ3が1/3になるが、3次高調波成分PS3の波数n3が、基本波成分PS1の波数n1の3倍となる。そのため、3次高調波成分を用いて画像を生成しても、距離分解能Δxは基本波成分を用いる場合と変わらない。図4(B)に示すように、波長も波数も基本波成分よりも小さい波を用いることができれば、距離分解能を向上させることができる。
そこで、以下で説明する本実施形態では、超音波の受信信号に対して、受信信号を構成する複数の第1基底波の結合係数特定処理を行い、特定された複数の結合係数と、第1基底波よりも波数が少ない第2基底波とに基づいて、受信信号から再構成信号への変換処理を行う。つまり、本実施形態では、後述する図8(A)に示すような受信波Xを再構成して、図8(D)に示す再構成波X’を生成する。この再構成波X’は、図8(C)に示す第2基底波により構成される。第2基底波は、元の受信波Xを構成する図8(B)の第1基底波と、波長が同じで、波数が少ない。そのため、距離分解能を向上させることができる。
2.システム構成例
次に、本実施形態の超音波測定装置の構成例を図5に示す。超音波測定装置100は、送信処理部110と、受信処理部120と、処理部130とを含む。
さらに、本実施形態の超音波画像装置の具体的な構成例を図6に示す。超音波画像装置は、超音波測定装置100と、超音波プローブ200と、表示部300と、を含む。また、図6に示す超音波測定装置100は、送信処理部110と、受信処理部120と、処理部130と、送受信切替スイッチ140と、DSC(Digital Scan Convertor)150と、制御回路160とを含む。
なお、超音波測定装置100及びこれを含む超音波画像装置は、図5及び図6の構成に限定されず、これらの一部の構成要素を省略したり、他の構成要素を追加したりするなどの種々の変形実施が可能である。また、本実施形態の超音波測定装置100及びこれを含む超音波画像装置の一部又は全部の機能は、通信により接続されたサーバーにより実現されてもよい。
次に各部で行われる処理について説明する。
超音波プローブ200は、超音波トランスデューサーデバイスを含む。
そして、超音波トランスデューサーデバイスは、走査面に沿って対象物をスキャンしながら、対象物に対して超音波ビームを送信すると共に、超音波ビームによる超音波エコーを受信する。圧電素子を用いるタイプを例にとれば、超音波トランスデューサーデバイスは、複数の超音波トランスデューサー素子(超音波素子アレイ)と、複数の開口がアレイ状に配置された基板とを有する。そして、超音波トランスデューサー素子としては、薄手の圧電素子と金属板(振動膜)を貼り合わせたモノモルフ(ユニモルフ)構造を用いたものを用いる。超音波トランスデューサー素子(振動素子)は、電気的な振動を機械的な振動に変換するものであるが、この場合には、圧電素子が面内で伸び縮みすると貼り合わせた金属板(振動膜)の寸法はそのままであるため反りが生じる。
また、超音波トランスデューサーデバイスでは、近隣に配置された数個の超音波トランスデューサー素子で一つのチャンネルを構成し、1回に複数のチャンネルを駆動しながら、超音波ビームを順次移動させるものであってもよい。
なお、超音波トランスデューサーデバイスとしては、圧電素子(薄膜圧電素子)を用いるタイプのトランスデューサーを採用できるが、本実施形態はこれに限定されない。例えばc‐MUT(Capacitive Micro-machined Ultrasonic Transducers)などの容量性素子を用いるタイプのトランスデューサーを採用してもよいし、バルクタイプのトランスデューサーを採用してもよい。超音波トランスデューサー素子及び超音波トランスデューサーデバイスのさらに詳細な説明については、後述する。
次に、送信処理部110は、対象物に対して超音波を送信する処理を行う。また、例えば図6に示す送信処理部110は、送信パルス発生器111と、送信遅延回路113とを含む。
具体的に、送信パルス発生器111は、送信パルス電圧を印加させ、超音波プローブ200を駆動させる。
また、送信遅延回路113は、送波ビームをフォーカシングする。そのために、送信遅延回路113は、送信パルス電圧の印加タイミングに関して、チャンネル間で時間差を与え、複数の振動素子から発生した超音波を集束させる。このように、遅延時間を変化させることにより、焦点距離を任意に変化させることが可能である。
また、送受信切替スイッチ140は、超音波の送受信の切り替え処理を行う。送受信切替スイッチ140は、送信時の振幅パルスが受信回路に入力されないように保護し、受信時の信号を受信回路に通す。
一方で、受信処理部120は、送信した超音波に対する超音波エコーの受信処理を行う。また、例えば図6に示す受信処理部120は、受信遅延回路121と、フィルター回路123と、メモリー125とを含む。
受信遅延回路121は、受波ビームをフォーカシングする。ある反射体からの反射波は球面上に広がるため、受信遅延回路121は、各振動子に到達する時間が同じになるように遅延時間を与え、遅延時間を考慮して反射波を加算する。
そして、フィルター回路123は、受信信号に対して帯域通過フィルターによりフィルター処理を行い、雑音を除去する。
また、メモリー125は、フィルター回路123から出力された受信信号を記憶するもので、その機能はRAM等のメモリーやHDDなどにより実現できる。
また、処理部130は、受信処理部120からの受信信号に対して処理を行う。例えば図6に示す処理部130は、ハーモニック処理部131と、再構成波生成部132と、検波処理部133と、対数変換処理部135と、ゲイン・ダイナミックレンジ調整部137と、STC(Sensitivity Time Control)139とを含む。
具体的に、ハーモニック処理部131は、前述したように、ハーモニック成分(主に3次高調波成分)の抽出処理を行う。
再構成波生成部132は、後に詳述するように、抽出されたハーモニック成分(主に3次高調波成分)に基づいて、受信信号を再構成信号へ変換する変換処理を行う。
そして、検波処理部133は、絶対値(整流)処理後に、低域通過フィルターをかけて、非変調信号を抽出する。
さらに、対数変換処理部135は、Log圧縮を行い、受信信号の信号強度の最大部分と最小部分を同時に確認しやすいように、表現形式を変換する。
そして、ゲイン・ダイナミックレンジ調整部137は、信号強度及び関心領域を調整する。具体的に、ゲイン調整処理では、Log圧縮後の入力信号に対して、直流成分を加える。また、ダイナミックレンジ調整処理では、Log圧縮後の入力信号に対して、任意の数を乗算する。
また、STC139は、深さに応じて増幅度(明るさ)を補正し、画面全体で一様な明るさの画像を取得する。
なお、処理部130の機能は、各種プロセッサー(CPU等)、ASIC(ゲートアレイ等)などのハードウェアや、プログラムなどにより実現できる。
そして、DSC150は、Bモード画像データに走査変換処理を行う。例えば、DSC150は、バイリニアなどの補間処理により、ライン信号を画像信号に変換する。
また、制御回路160は、送信パルス発生器111と、送信遅延回路113と、受信遅延回路121と、送受信切替スイッチ140と、ハーモニック処理部131の制御を行う。
また、表示部300は、再構成信号に基づいて生成された表示用画像データを表示する。表示部300は、例えば液晶ディスプレイや有機ELディスプレイ、電子ペーパーなどにより実現できる。
ここで、本実施形態の超音波画像装置(広義には電子機器)の具体的な機器構成の例を図7(A)〜図7(C)に示す。図7(A)はハンディタイプの超音波画像装置の例であり、図7(B)は据置タイプの超音波画像装置の例である。図7(C)は超音波プローブ200が本体に内蔵された一体型の超音波画像装置の例である。
図7(A)、図7(B)の超音波画像装置は、超音波プローブ200と超音波測定装置100を含み、超音波プローブ200と超音波測定装置100はケーブル210により接続される。超音波プローブ200の先端部分には、プローブヘッド220が設けられており、超音波測定装置本体101には、画像を表示する表示部300が設けられている。図7(C)では、表示部300を有する超音波画像装置に超音波プローブ200が内蔵されている。図7(C)の場合、超音波画像装置は、例えばスマートフォンなどの汎用の携帯情報端末により実現できる。
3.処理の詳細
3.1. 再構成波の生成処理
本実施形態の処理部130は、送信処理部110により送信された送信パルス信号に対応する受信信号に対して、受信信号を構成する複数の第1基底波の結合係数特定処理を行い、結合係数特定処理により特定された複数の結合係数と、第1基底波よりも波数が少ない第2基底波とに基づいて、受信信号から再構成信号への変換処理を行う。
これにより、生成した再構成信号に基づいて、対象物内の測定を行うことができる。この再構成信号は、第1基底波よりも波数が少ない第2基底波から構成されている。なお、第1基底波及び第2基底波の波数は、必ずしも自然数である必要はなく、例えば1.5波等であってもよい。
ここで、再構成信号(再構成波)とは、受信信号(受信波)に含まれる対象物内の各点散乱体からの各反射信号成分(反射波成分)を、元の反射信号成分よりもパルス幅の短い信号(波形)にそれぞれ置き換えて、置き換えた後の信号(波形)を元の反射信号成分の受信タイミングと同じタイミングで重ね合わせ直した信号である。再構成信号は、例えば後述する図8(D)のような再構成波である。つまり、前述した第1基底波が、各点散乱体からの各反射信号成分(反射波成分)に対応し、第2基底波が、第1基底波を置換する、元の反射信号成分よりもパルス幅の短い信号(波形)に対応する。
そして、第1基底波は、受信信号のうち、対象物内の所与の深さに存在する点散乱体からの反射波成分に対応する波であり、例えば後述する図8(B)に示すような波である。後述するように、受信信号に第1基底波成分が含まれているか否かを判断することにより、その第1基底波に対応する点散乱体が対象物に含まれているか否かを判断することができる。また、その第1基底波に対応する点散乱体が対象物に含まれていると判断できる場合には、受信信号に含まれる第1基底波成分の信号強度(反射強度)により、点散乱体の反射特性等を特定することができる。この反射強度を特定する処理が、結合係数特定処理である。なお、第2基底波については後述する。
このように、本実施形態では、受信波の各反射波成分をパルス幅の短い第2基底波に置換して、受信波を再構成するため、再構成信号に基づく対象物の測定結果における距離分解能を向上させることができる。よって、超音波による対象物の測定結果の方位分解能だけでなく、距離分解能も向上させることが可能となる。
具体的には、処理部130は、例えば図8(A)に示すように、受信信号に対応する高調波(受信波X)に対し、図8(B)に示すように第1基底波(S)の結合係数特定処理を行った後に、図8(C)に示すように受信波Xを構成する第1基底波(S)を第2基底波(S’)へ変換する変換処理を行って、図8(D)に示すように再構成信号として、第2基底波(S’)による再構成波X’を生成する。なお、図8(A)に示す高調波Xは、元の受信波に対してフィルター処理を行って抽出された高調波等であってもよい。
これにより、第1基底波よりも波数が少ない第2基底波から構成される再構成波に、受信波を変換し、対象物の測定結果における距離分解能を向上させること等が可能になる。
ここで、第1基底波は、受信信号から抽出可能な高調波である。これにより、受信信号を複数の第1基底波に分解すること等が可能になる。
例えば、図8(B)の例では、第1基底波は、基底関数Sにより表される波である。本実施形態では、第1基底波は1つではなく、複数の波である。例えば図8(B)の例では、複数の第1基底波は、M個の第1基底波であるものとし、Mは2以上の整数であり、iは、1≦i≦Mの整数であるものとする。
そして、M個の第1基底波のうちの第iの第1基底波は、送信パルス信号のパルス幅又は受信信号のパルス幅に対応する位相差よりも短い位相差で、第(i+1)の第1基底波と相互に位相がずれている。
これにより、送信パルス信号のパルス幅又は受信信号のパルス幅に対応する位相差よりも短い距離の距離分解能で、対象物を測定すること等が可能になる。
さらに、受信波Xを構成する複数の第1基底波の結合係数特定処理を行うと、図8(B)に示すように、各第1基底波Sに対して各結合係数aが求められる。この結合係数aは、対応する第1基底波Sが受信波Xにどの程度の割合で含まれているかを決定付ける値である。つまり、受信波Xは下式(2)に示すように、各第1基底波と各結合係数の積の和により表される。
また、第2基底波は、第1基底波と位相差が同じで、波数が少ない波である。なお、この位相差とは、第iの第1基底波と第(i+1)の第1基底波の間の位相差のことである。
そして、各第2基底波S’に対しては、受信波Xを複数の第1基底波Sに分解する際(結合係数特定処理)に求めた各結合係数aを対応付ける。よって、受信波Xは下式(3)により表される。
これにより、受信波を構成する基底波のパルス幅を短くすること等が可能になる。
次に、図9のフローチャートを用いて、本実施形態の処理の流れについて説明する。
まず、走査線番号nの初期値を1に設定する(S101)。
次に、送信パルス発生器111が、位相0°のパルス電圧を生成する(S102)。
そして、送信遅延回路113が送信フォーカス制御を行い(S103)、超音波プローブ200が、生成されたパルス電圧に対応する超音波ビームを対象物に対して出射する(S104)。さらに、超音波プローブ200は、出射した超音波ビームが対象物に反射し、返ってきた超音波エコーを受信する(S104)。
これに対して、受信遅延回路121は受信フォーカス制御を行い(S105)、フィルター回路123が、受信フォーカス制御後の受信信号に対してBPF(バンドパスフィルター)処理を行い(S106)、BPF処理後の受信信号をメモリー125に保存する(S107)。
次に、送信パルス発生器111が、位相180°のパルス電圧を生成する(S108)。そして、この位相反転後のパルスについて、前述したステップS102〜ステップS107の処理と同様の処理を行う(S109〜S113)。
その後に、全ての走査線についてステップS102〜S113の処理を行ったか否かを判断する(S114)。具体的には、現在の走査線番号nが全走査線数Nよりも小さいか否かを判定する。
全ての走査線についてステップS102〜S113の処理を行っていないと判断した場合、すなわち、現在の走査線番号nが全走査線数Nよりも小さいと判定した場合には、現在の走査線番号nに1を加算して(S115)、再度ステップS102〜S114の処理を行う。
一方、ステップS115において、全ての走査線についてステップS102〜S114の処理を行ったと判断した場合、すなわち、現在の走査線番号nが全走査線数Nと等しいと判定した場合には、ハーモニック処理部131が高調波成分(ハーモニック成分)の抽出処理を行う(S116)。具体的に本抽出処理では、前述した図2(A)〜図2(C)に示すように、位相0°の送信波に対応する受信波と、位相180°の送信波に対応する受信波の減算処理を行って、基本波と3次高調波を抽出する。その後に、図3(A)及び図3(B)に示すように、抽出した基本波と3次高調波に対して、ハイパスフィルター処理を行って、3次高調波のみを抽出する。
次に、再構成波生成部132が、抽出された3次高調波に基づいて再構成波を生成する(S117)。ここで、本実施形態における再構成波の生成処理の流れを図10のフローチャートに示す。
まず、再構成波生成部132が、不図示のメモリーから基底関数(第1基底波及び第2基底波)を読み込む(S201)。次に、再構成波生成部132が、周波数フィルター処理(BPF)を行い、メモリーから読み出した第1基底波に基づいて、前述した図8(B)に示すように、3次高調波における第1基底波成分を抽出する(S202)。
そして、再構成波生成部132は、抽出結果に基づいて、3次高調波を構成する第1基底波の結合係数の特定処理を行う(S203)。
具体的には、再構成波生成部132(処理部130)は、結合係数特定処理として、受信信号のデコンボリューション(deconvolution)処理を行う。
これにより、受信波を構成する第1基底波の結合係数を特定すること等が可能になる。
その後に、再構成波生成部132は、各第1基底波を位相が同じ各第2基底波に置き換えて、特定した結合係数と第2基底波に基づいて、再構成波を生成する(S204)。言い換えれば、受信信号を再構成信号へ変換する変換処理を行う。
具体的に、処理部130は、再構成信号の変換処理として、第2基底波のコンボリューション(convolution)処理を行う。
これにより、受信波又は受信波から抽出可能な高調波を構成する第1基底波の波数を削減した第2基底波によって、再構成波を生成すること等が可能になる。
そして、検波処理部133が、生成した再構成波に対して、絶対値(整流)処理後に、低域通過フィルターをかけて、非変調信号を抽出し(S118)、対数変換処理部135が、対数変換処理を行う(S119)。
そして、ゲイン・ダイナミックレンジ調整部137が、信号強度及び関心領域を調整し(S120)、STC139が、深さに応じて増幅度(明るさ)を補正する(S121)。
さらに、DSC150が、走査変換処理を行って、Bモード画像データ(表示用画像データ)を生成し(S122)、表示部300が生成された表示用画像データを表示し(S123)、処理を終了する。
3.2. 第1基底波及び第2基底波の生成処理
図9及び図10のフローチャートの処理で用いた第1基底波及び第2基底波は、これらの処理を行う前に(前処理で)予め生成して、不図示のメモリーに記憶しておく必要がある。以下では、第1基底波及び第2基底波の生成処理の流れについて、図11のフローチャートを用いて説明する。
まず、超音波ファントムのワイヤー初期位置Pを設定する(S301)。ここでは、初期位置Pを、超音波測定装置の超音波プローブの測定面に最も近い位置Pminに設定する。またワイヤー位置とは、点散乱体の位置のことである。
そして、送信処理部110が、互いに位相を反転させた2つのパルス信号を対象物(超音波ファントム)に送信し、受信処理部120が、送信された2つのパルス信号に対応する2つの受信信号(RFデータ)をそれぞれ受信する(S302)。
そして、処理部130は、前述した図2(B)及び図2(C)に示すように、送信された2つのパルス信号に対応する2つの受信信号に基づいて減算処理を行って、1つの差分信号を求め、前述した図3(A)及び図3(B)に示すように、求めた差分信号に対して第1フィルター処理を行って、高調波成分(3次高調波成分)を抽出する。ここで、第1フィルター処理とは、例えばハイパスフィルター処理やバンドパスフィルター処理等である。
さらに、処理部130は、抽出された高調波成分に基づいて、所与の測定点に配置された点散乱体からの反射波成分に対応する高調波を、第1基底波として求める(S303)。
前述した図8(A)及び図8(B)の例では、受信波(高調波)Xが、対象物内の様々な点散乱体からの反射波成分を含んでおり、予め特定されている第1基底波と比較しなければ、受信波Xに含まれる所与の測定点からの第1基底波成分を特定することができない。これに対して、本処理では第1フィルター処理後の高調波成分が、超音波ファントム内の所与の測定点に配置された点散乱体からの反射波成分だけを含んでいる。この高調波成分には、ノイズ等も含まれるが、容易にノイズ成分と切り分け可能である。
そのため、対象物内の所与の測定点の点散乱体からの反射波成分に対応する第1基底波を特定すること等が可能になる。そして、処理部130は、特定した第1基底波をメモリーに保存する(S304)。
さらに、処理部130は、図12(A)に示す減算処理により求めた差分信号に対して第2フィルター処理を行って、図12(B)に示すように基本波成分を抽出する(S305)。そして、処理部130は、抽出された基本波成分に基づいて、所与の測定点に配置された点散乱体からの反射波成分に対応する基本波を求め、求めた基本波に対して、図12(C)に示すように時間成分の圧縮処理を行って、第2基底波を求める(S306)。時間成分を1/yに圧縮する処理を行った場合には、図12(C)の第2基底波の波長λ3は、図12(B)の基本波の波長λ1の1/y(yは正の数)になる。
このように、対象物内の所与の測定点の点散乱体からの反射波成分に対応する第1基底波の波長を短くした第2基底波を特定すること等が可能になる。
また、第2基底波は、受信信号から抽出可能な基本波に対して、時間成分の圧縮処理を行って得られる。
これにより、減算処理、フィルター処理及び時間成分の圧縮処理等の簡易な処理で、第2基底波を特定すること等が可能になる。そして、処理部130は、特定した第2基底波をメモリーに保存する(S307)。
その後に、処理部130は、ワイヤー位置Pが測定範囲の最大値Pmaxよりも大きいか否かを判断し(S308)、ワイヤー位置Pが測定範囲の最大値Pmax以下であると判断した場合には、ワイヤー位置Pを下式(4)に基づいて更新して(S309)、ステップS302に戻る。なお、式(4)において、Kは所与の定数であり、λは波長である。
一方、処理部130は、ワイヤー位置Pが測定範囲の最大値Pmaxよりも大きいと判断した場合には、処理を終了する。
以上をまとめると、例えば図13(A)及び図13(B)に示すように、M個の第1基底波のうちの第iの第1基底波sは、第iの測定点に配置された第iの点散乱体SPiからの超音波の受信信号に対応する高調波である。なお、Mは2以上の整数であり、iは、1≦i≦Mの整数である。
さらに、M個の第1基底波のうちの第(i+1)の第1基底波s(i+1)は、超音波の送信点TPから、第iの測定点よりも遠い位置である第(i+1)の測定点に配置された第(i+1)の点散乱体SP(i+1)からの超音波の受信信号に対応する高調波である。
これにより、設定した各測定点の間隔に対応する距離分解能で、受信波から第1基底波成分を抽出すること等が可能になる。
そして、例えば図13(A)及び図13(C)に示すように、N個の第2基底波のうちの第jの第2基底波s’は、第jの測定点に配置された第jの点散乱体SPjからの超音波の受信信号に対応する高調波である。なお、Nは2以上の整数であり、jは、1≦j≦Nの整数である。本例では、M=Nであるが、M≠Nであってもよい。
さらに、N個の第2基底波のうちの第(j+1)の第2基底波s’(j+1)は、超音波の送信点TPから、第jの測定点よりも遠い位置である第(j+1)の測定点に配置された第(j+1)の点散乱体SP(j+1)からの超音波の受信信号に対応する高調波である。
これにより、超音波による対象物の測定結果による距離分解能を、設定した各測定点の間隔まで向上させること等が可能になる。
ただし、本実施形態の第1基底波及び第2基底波の生成処理は、以上で説明した処理に限定されない。例えば、各第2基底波を、対応する各第1基底波の波数を削減することにより生成してもよい。さらに、第1基底波及び第2基底波の生成処理を行わずに、あらかじめメモリーに記憶された第1基底波及び第2基底波を用いても良い。なお、上記の超音波プローブ及び測定点は、シミュレーション上のものであってもよい。
3.3. 測定結果
図14(A)及び図14(B)に、本実施形態の測定結果の一例を示す。本例では、図14(A)に示すように、波長λ=0.44mmで、波数1、周波数3.5MHzの超音波パルスを対象物に送信する例を示す。
図14(A)の例では、図面の上方向を対象物の深さ方向とし、対象物内に3つの点散乱体(SP1〜SP3)が存在するものとする。そして、超音波トランスデューサー素子と点散乱体SP1は30mm離れており、点散乱体SP1と点散乱体SP2は(4/10)λ離れており、点散乱体SP2と点散乱体SP3は(5/10)λ離れているものとする。さらに、点散乱体SP1からの反射強度は0.5、点散乱体SP2からの反射強度は1.0、点散乱体SP3からの反射強度は0.7であるものとする。
この場合に、位相反転法とフィルター法のみを用いてBモード画像を生成すると、同図の画像BIM1が生成される。この画像BIM1では、各点散乱体からの反射波が互いに重複しているため、画像全体の色味の変化が少なく、3つの点散乱体の詳細な位置を特定することは困難である。
一方、前述した本実施形態による再構成波に基づいて、Bモード画像を生成すると、同図の画像BIM2が生成される。この画像BIM2では、色層L1が点散乱体SP1からの反射に対応しており、色層L2が点散乱体SP2からの反射に対応しており、色層L3が点散乱体SP3からの反射に対応していることがはっきりと分かる。つまり、位相反転法とフィルター法のみを用いる手法に比べて、距離分解能が向上している。
また、処理部130は、変換処理後に得られた再構成信号に対して、包絡線検波処理を行ってもよい。
例えば、図14(B)に示すように、位相反転法とフィルター法による処理を受信信号RSに行った後の波形に対して、包絡線検波処理を行うと、波形AS1が得られる。波形AS1では、大きな2つの山が確認できるが、ここから対象物に3つの点散乱体が含まれていると判断することは困難である。
これに対し、受信信号RSに本実施形態の処理を行って、波形ARSを得て、波形ARSに対して、包絡線検波処理を行うと、波形AS2が得られる。波形AS2では、3つの山が確認でき、これらが点散乱体からの反射に対応する山であると容易に判断することができる。また、波形AS1に比べて、波形AS2の山の位置が実際の対象物における点散乱体の分布と近くなっている。
このように、ユーザーが測定結果を判別しやすい波形データを表示部300に表示すること等が可能になる。
4.超音波トランスデューサー素子
図15(A)〜図15(C)に、超音波トランスデューサーデバイスの超音波トランスデューサー素子10の構成例を示す。この超音波トランスデューサー素子10は、振動膜(メンブレン、支持部材)50と圧電素子部とを有する。圧電素子部は、第1電極層(下部電極)21、圧電体層(圧電体膜)30、第2電極層(上部電極)22を有する。
図15(A)は、基板(シリコン基板)60に形成された超音波トランスデューサー素子10の、素子形成面側の基板60に垂直な方向から見た平面図である。図15(B)は、図15(A)のA−A’に沿った断面を示す断面図である。図15(C)は、図15(A)のB−B’に沿った断面を示す断面図である。
第1電極層21は、振動膜50の上層に例えば金属薄膜で形成される。この第1電極層21は、図15(A)に示すように素子形成領域の外側へ延長され、隣接する超音波トランスデューサー素子10に接続される配線であってもよい。
圧電体層30は、例えばPZT(ジルコン酸チタン酸鉛)薄膜により形成され、第1電極層21の少なくとも一部を覆うように設けられる。なお、圧電体層30の材料は、PZTに限定されるものではなく、例えばチタン酸鉛(PbTiO3)、ジルコン酸鉛(PbZrO3)、チタン酸鉛ランタン((Pb、La)TiO3)などを用いてもよい。
第2電極層22は、例えば金属薄膜で形成され、圧電体層30の少なくとも一部を覆うように設けられる。この第2電極層22は、図15(A)に示すように素子形成領域の外側へ延長され、隣接する超音波トランスデューサー素子10に接続される配線であってもよい。
振動膜(メンブレン)50は、例えばSiO2薄膜とZrO2薄膜との2層構造により開口40を塞ぐように設けられる。この振動膜50は、圧電体層30及び第1、第2電極層21、22を支持すると共に、圧電体層30の伸縮に従って振動し、超音波を発生させることができる。
開口40は、基板60(シリコン基板)の裏面(素子が形成されない面)側から反応性イオンエッチング(RIE)等によりエッチングすることで形成される。この開口40の開口部45のサイズによって超音波の共振周波数が決定され、その超音波は圧電体層30側(図15(A)において紙面奥から手前方向)に放射される。
超音波トランスデューサー素子10の下部電極(第1電極)は、第1電極層21により形成され、上部電極(第2電極)は、第2電極層22により形成される。具体的には、第1電極層21のうちの圧電体層30に覆われた部分が下部電極を形成し、第2電極層22のうちの圧電体層30を覆う部分が上部電極を形成する。即ち、圧電体層30は、下部電極と上部電極に挟まれて設けられる。
5.超音波トランスデューサーデバイス
図16に、超音波トランスデューサーデバイス(素子チップ)の構成例を示す。本構成例の超音波トランスデューサーデバイスは、複数の超音波トランスデューサー素子群UG1〜UG64、駆動電極線DL1〜DL64(広義には第1〜第nの駆動電極線。nは2以上の整数)、コモン電極線CL1〜CL8(広義には第1〜第mのコモン電極線。mは2以上の整数)を含む。なお、駆動電極線の本数(n)やコモン電極線の本数(m)は、図16に示す本数には限定されない。
複数の超音波トランスデューサー素子群UG1〜UG64は、第2の方向D2(スキャン方向)に沿って64列に配置される。UG1〜UG64の各超音波トランスデューサー素子群は、第1の方向D1(スライス方向)に沿って配置される複数の超音波トランスデューサー素子を有する。
図17(A)に、超音波トランスデューサー素子群UG(UG1〜UG64)の例を示す。図17(A)では、超音波トランスデューサー素子群UGは第1〜第4の素子列により構成される。第1の素子列は、第1の方向D1に沿って配置される超音波トランスデューサー素子UE11〜UE18により構成され、第2の素子列は、第1の方向D1に沿って配置される超音波トランスデューサー素子UE21〜UE28により構成される。第3の素子列(UE31〜UE38)、第4の素子列(UE41〜UE48)も同様である。これらの第1〜第4の素子列には、駆動電極線DL(DL1〜DL64)が共通接続される。また、第1〜第4の素子列の超音波トランスデューサー素子にはコモン電極線CL1〜CL8が接続される。
そして図17(A)の超音波トランスデューサー素子群UGが、超音波トランスデューサーデバイスの1チャンネルを構成する。即ち、駆動電極線DLが1チャンネルの駆動電極線に相当し、送信回路からの1チャンネルの送信信号は駆動電極線DLに入力される。また駆動電極線DLからの1チャンネルの受信信号は駆動電極線DLから出力される。なお、1チャンネルを構成する素子列数は図17(A)のような4列には限定されず、4列よりも少なくてもよいし、4列よりも多くてもよい。例えば図17(B)に示すように、素子列数は1列であってもよい。
図16に示すように、駆動電極線DL1〜DL64(第1〜第nの駆動電極線)は、第1の方向D1に沿って配線される。駆動電極線DL1〜DL64のうちの第j(jは1≦j≦nである整数)の駆動電極線DLj(第jのチャンネル)は、第jの超音波トランスデューサー素子群UGjの超音波トランスデューサー素子が有する第1の電極(例えば下部電極)に接続される。
超音波を出射する送信期間には、送信信号VT1〜VT64が駆動電極線DL1〜DL64を介して超音波トランスデューサー素子に供給される。また、超音波エコー信号を受信する受信期間には、超音波トランスデューサー素子からの受信信号VR1〜VR64が駆動電極線DL1〜DL64を介して出力される。
コモン電極線CL1〜CL8(第1〜第mのコモン電極線)は、第2の方向D2に沿って配線される。超音波トランスデューサー素子が有する第2の電極は、コモン電極線CL1〜CL8のうちのいずれかに接続される。具体的には、例えば図16に示すように、コモン電極線CL1〜CL8のうちの第i(iは1≦i≦mである整数)のコモン電極線CLiは、第i行に配置される超音波トランスデューサー素子が有する第2の電極(例えば上部電極)に接続される。
コモン電極線CL1〜CL8には、コモン電圧VCOMが供給される。このコモン電圧VCOMは一定の直流電圧であればよく、0V、即ちグランド電位(接地電位)でなくてもよい。
そして送信期間では、送信信号電圧とコモン電圧との差の電圧が超音波トランスデューサー素子に印加され、所定の周波数の超音波が放射される。
なお、超音波トランスデューサー素子の配置は、図16に示すマトリックス配置に限定されず、いわゆる千鳥配置等であってもよい。
また図17(A)〜図17(B)では、1つの超音波トランスデューサー素子が送信素子及び受信素子の両方に兼用される場合について示したが、本実施形態はこれに限定されない。例えば送信素子用の超音波トランスデューサー素子、受信素子用の超音波トランスデューサー素子を別々に設けて、アレイ状に配置してもよい。
なお、本実施形態の超音波測定装置及び超音波画像装置等は、その処理の一部または大部分をプログラムにより実現してもよい。この場合には、CPU等のプロセッサーがプログラムを実行することで、本実施形態の超音波測定装置及び超音波画像装置等が実現される。具体的には、非一時的な情報記憶装置に記憶されたプログラムが読み出され、読み出されたプログラムをCPU等のプロセッサーが実行する。ここで、情報記憶装置(コンピューターにより読み取り可能な装置)は、プログラムやデータなどを格納するものであり、その機能は、光ディスク(DVD、CD等)、HDD(ハードディスクドライブ)、或いはメモリー(カード型メモリー、ROM等)などにより実現できる。そして、CPU等のプロセッサーは、情報記憶装置に格納されるプログラム(データ)に基づいて本実施形態の種々の処理を行う。即ち、情報記憶装置には、本実施形態の各部としてコンピューター(操作部、処理部、記憶部、出力部を備える装置)を機能させるためのプログラム(各部の処理をコンピューターに実行させるためのプログラム)が記憶される。
また、本実施形態の超音波測定装置及び超音波画像装置等は、プロセッサーとメモリーを含んでも良い。ここでのプロセッサーは、例えばCPU(Central Processing Unit)であってもよい。ただし、プロセッサーはCPUに限定されるものではなく、GPU(Graphics Processing Unit)、或いはDSP(Digital Signal Processor)等、各種プロセッサーを用いることが可能である。また、プロセッサーはASIC(Application Specific Integrated Circuit)によるハードウェア回路でもよい。また、メモリーはコンピューターにより読み取り可能な命令を格納するものであり、当該命令がプロセッサーにより実行されることで、本実施形態に係る超音波測定装置及び超音波画像装置等の各部が実現されることになる。ここでのメモリーは、SRAM(Static Random Access Memory)、DRAM(Dynamic Random Access Memory)などの半導体メモリーであってもよいし、レジスターやハードディスク等でもよい。また、ここでの命令は、プログラムを構成する命令セットの命令でもよいし、プロセッサーのハードウェア回路に対して操作を指示する命令であってもよい。
以上のように本実施形態について詳細に説明したが、本発明の新規事項および効果から実体的に逸脱しない多くの変形が可能であることは当業者には容易に理解できるであろう。従って、このような変形例はすべて本発明の範囲に含まれるものとする。例えば、明細書又は図面において、少なくとも一度、より広義または同義な異なる用語と共に記載された用語は、明細書又は図面のいかなる箇所においても、その異なる用語に置き換えることができる。また、超音波測定装置、超音波画像装置及び超音波測定方法の構成、動作も本実施形態で説明したものに限定されず、種々の変形実施が可能である。
10 超音波トランスデューサー素子、21 第1電極層、22 第2電極層、
30 圧電体層、40 開口、45 開口部、50 振動膜、60 基板、
100 超音波測定装置、101 超音波測定装置本体、110 送信処理部、
111 送信パルス発生器、113 送信遅延回路、120 受信処理部、
121 受信遅延回路、123 フィルター回路、125 メモリー、130 処理部、
131 ハーモニック処理部、132 再構成波生成部、133 検波処理部、
135 対数変換処理部、137 ゲイン・ダイナミックレンジ調整部、
139 STC、140 送受信切替スイッチ、150 DSC、160 制御回路、
200 超音波プローブ、210 ケーブル、220 プローブヘッド、300 表示部

Claims (15)

  1. 対象物に対して超音波を送信する処理を行う送信処理部と、
    送信した前記超音波に対する超音波エコーの受信処理を行う受信処理部と、
    前記受信処理部からの受信信号に対して処理を行う処理部と、
    を含み、
    前記処理部は、
    前記送信処理部により送信された送信パルス信号に対応する受信信号に対して、前記受信信号を構成する複数の第1基底波の結合係数特定処理を行い、
    前記結合係数特定処理により特定された複数の結合係数と、前記第1基底波よりも波数が少ない第2基底波とに基づいて、前記受信信号から再構成信号への変換処理を行うことを特徴とする超音波測定装置。
  2. 請求項1において、
    前記処理部は、
    前記受信信号に対応する高調波に対し、前記第1基底波の前記結合係数特定処理を行った後に、前記変換処理を行って、前記再構成信号として、前記第2基底波による再構成波を生成することを特徴とする超音波測定装置。
  3. 請求項1又は2において、
    前記複数の第1基底波は、
    M個(Mは2以上の整数)の第1基底波であり、
    前記M個の第1基底波のうちの第i(iは、1≦i≦Mの整数)の第1基底波は、
    第iの測定点に配置された第iの点散乱体からの前記超音波の受信信号に対応する高調波であり、
    前記M個の第1基底波のうちの第(i+1)の第1基底波は、
    前記超音波の送信点から、前記第iの測定点よりも遠い位置である第(i+1)の測定点に配置された第(i+1)の点散乱体からの超音波の受信信号に対応する高調波であることを特徴とする超音波測定装置。
  4. 請求項1乃至3のいずれかにおいて、
    前記第1基底波は、
    前記受信信号から抽出可能な高調波であることを特徴とする超音波測定装置。
  5. 請求項1乃至4のいずれかにおいて、
    前記送信処理部は、
    互いに位相を反転させた2つのパルス信号を前記対象物に送信し、
    前記処理部は、
    送信された前記2つのパルス信号に対応する2つの受信信号に基づいて減算処理を行って、1つの差分信号を求め、
    求めた前記差分信号に対して第1フィルター処理を行って、高調波成分を抽出し、
    抽出された前記高調波成分に基づいて、所与の測定点に配置された点散乱体からの反射波成分に対応する高調波を、前記第1基底波として求めることを特徴とする超音波測定装置。
  6. 請求項1乃至5のいずれかにおいて、
    前記複数の第1基底波は、
    M個(Mは2以上の整数)の第1基底波であり、
    前記M個の第1基底波のうちの第i(iは、1≦i≦Mの整数)の第1基底波は、
    前記送信パルス信号のパルス幅又は前記受信信号のパルス幅に対応する位相差よりも短い位相差で、第(i+1)の第1基底波と相互に位相がずれていることを特徴とする超音波測定装置。
  7. 請求項1乃至6のいずれかにおいて、
    前記処理部は、
    複数の第2基底波に基づいて、前記変換処理を行い、
    前記複数の第2基底波は、
    N個(Nは2以上の整数)の第2基底波であり、
    前記N個の第2基底波のうちの第j(jは、1≦j≦Nの整数)の第2基底波は、
    第jの測定点に配置された第jの点散乱体からの前記超音波の受信信号に対応する高調波であり、
    前記N個の第2基底波のうちの第(j+1)の第2基底波は、
    前記超音波の送信点から、前記第jの測定点よりも遠い位置である第(j+1)の測定点に配置された第(j+1)の点散乱体からの超音波の受信信号に対応する高調波であることを特徴とする超音波測定装置。
  8. 請求項1乃至7のいずれかにおいて、
    前記第2基底波は、
    前記受信信号から抽出可能な基本波に対して、時間成分の圧縮処理を行って得られることを特徴とする超音波測定装置。
  9. 請求項1乃至8のいずれかにおいて、
    前記送信処理部は、
    互いに位相を反転させた2つのパルス信号を前記対象物に送信し、
    前記処理部は、
    送信された前記2つのパルス信号に対応する2つの受信信号に基づいて減算処理を行って、1つの差分信号を求め、
    求めた前記差分信号に対して第2フィルター処理を行って、基本波成分を抽出し、
    抽出された前記基本波成分に基づいて、所与の測定点に配置された点散乱体からの反射波成分に対応する基本波を求め、
    求めた前記基本波に対して時間成分の圧縮処理を行って、前記第2基底波を求めることを特徴とする超音波測定装置。
  10. 請求項1乃至9のいずれかにおいて、
    前記第2基底波は、
    前記第1基底波と位相差が同じで、波数が少ないことを特徴とする超音波測定装置。
  11. 請求項1乃至10のいずれかにおいて、
    前記処理部は、
    前記変換処理後に得られた前記再構成信号に対して、包絡線検波処理を行うことを特徴とする超音波測定装置。
  12. 請求項1乃至11のいずれかにおいて、
    前記処理部は、
    前記結合係数特定処理として、前記受信信号のデコンボリューション処理を行うことを特徴とする超音波測定装置。
  13. 請求項1乃至12のいずれかにおいて、
    前記処理部は、
    前記再構成信号の前記変換処理として、前記第2基底波のコンボリューション処理を行うことを特徴とする超音波測定装置。
  14. 請求項1乃至13のいずれかに記載の超音波測定装置と、
    前記再構成信号に基づいて生成された表示用画像データを表示する表示部と、
    を含むことを特徴とする超音波画像装置。
  15. 対象物に対して超音波を送信する処理を行い、
    送信した前記超音波に対する超音波エコーの受信処理を行い、
    送信された送信パルス信号に対応する受信信号に対して、前記受信信号を構成する複数の第1基底波の結合係数特定処理を行い、
    前記結合係数特定処理により特定された複数の結合係数と、前記第1基底波よりも波数が少ない第2基底波とに基づいて、前記受信信号から再構成信号への変換処理を行うことを特徴とする超音波測定方法。
JP2014222474A 2014-10-31 2014-10-31 超音波測定装置、超音波画像装置及び超音波測定方法 Active JP6365243B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2014222474A JP6365243B2 (ja) 2014-10-31 2014-10-31 超音波測定装置、超音波画像装置及び超音波測定方法
US14/928,076 US10849598B2 (en) 2014-10-31 2015-10-30 Ultrasonic measurement apparatus, ultrasonic imaging apparatus, and ultrasonic measurement method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2014222474A JP6365243B2 (ja) 2014-10-31 2014-10-31 超音波測定装置、超音波画像装置及び超音波測定方法

Publications (2)

Publication Number Publication Date
JP2016086958A JP2016086958A (ja) 2016-05-23
JP6365243B2 true JP6365243B2 (ja) 2018-08-01

Family

ID=55851357

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014222474A Active JP6365243B2 (ja) 2014-10-31 2014-10-31 超音波測定装置、超音波画像装置及び超音波測定方法

Country Status (2)

Country Link
US (1) US10849598B2 (ja)
JP (1) JP6365243B2 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7317332B2 (ja) * 2017-10-19 2023-07-31 ザインエレクトロニクス株式会社 送信装置および送受信システム
JP6939414B2 (ja) * 2017-10-26 2021-09-22 セイコーエプソン株式会社 超音波デバイス、及び超音波測定装置
JP6922651B2 (ja) * 2017-10-26 2021-08-18 セイコーエプソン株式会社 超音波デバイス、及び超音波測定装置
JP7298186B2 (ja) * 2019-02-26 2023-06-27 セイコーエプソン株式会社 超音波計測装置、及び超音波計測方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3575252B2 (ja) * 1997-11-13 2004-10-13 Jfeスチール株式会社 電磁超音波計測方法および装置
JP4744727B2 (ja) 2001-06-06 2011-08-10 株式会社東芝 超音波画像診断装置
US7729201B2 (en) * 2003-03-25 2010-06-01 Chester Wildey Method and apparatus for spread spectrum distance and velocity profile measurement
JP2006166957A (ja) * 2004-12-13 2006-06-29 Matsushita Electric Ind Co Ltd 超音波診断装置
EP1949298A4 (en) * 2005-09-19 2010-04-07 Univ Virginia SYSTEM AND METHOD FOR ADAPTIVE BEAM FORMATION FOR RECONSTITUTING IMAGE AND / OR LOCATING A TARGET / SOURCE
JP4627556B2 (ja) * 2008-08-08 2011-02-09 アロカ株式会社 超音波診断装置
JP2011010793A (ja) 2009-07-01 2011-01-20 Konica Minolta Medical & Graphic Inc 超音波画像診断装置

Also Published As

Publication number Publication date
JP2016086958A (ja) 2016-05-23
US20160120513A1 (en) 2016-05-05
US10849598B2 (en) 2020-12-01

Similar Documents

Publication Publication Date Title
US9867594B2 (en) Ultrasonic measurement apparatus, ultrasonic image apparatus, and ultrasonic measurement method
JP5389722B2 (ja) 超音波診断装置及びその作動方法
JP6111870B2 (ja) 超音波測定装置、超音波画像装置及び超音波画像の処理方法
JP6365243B2 (ja) 超音波測定装置、超音波画像装置及び超音波測定方法
JP2015016144A (ja) 超音波測定装置、超音波画像装置及び超音波測定方法
US20180299537A1 (en) Beamforming apparatus, beamforming method, and ultrasonic imaging apparatus
JP6398614B2 (ja) 超音波測定装置、超音波診断装置及び超音波測定方法
US20120245467A1 (en) Ultrasound diagnostic apparatus and ultrasound image producing method
JP6402590B2 (ja) 超音波測定装置、超音波画像装置及び超音波測定方法
JP2014188149A (ja) 超音波測定装置、超音波画像装置及び超音波測定方法
JP7005206B2 (ja) 超音波診断装置及び超音波イメージングプログラム
JP5637960B2 (ja) 超音波診断装置および超音波画像生成方法
JP6089741B2 (ja) 超音波測定装置、超音波画像装置及び超音波測定方法
JP6186737B2 (ja) 超音波測定装置、超音波画像装置及び超音波測定方法
JP6237007B2 (ja) 超音波測定装置、超音波画像装置及び超音波測定方法
JP7273519B2 (ja) 超音波診断装置、及び学習プログラム
JP5917718B2 (ja) 超音波診断装置
JP5851345B2 (ja) 超音波診断装置およびデータ処理方法
JP2005245479A (ja) 超音波診断装置
CN106175823A (zh) 产生谐波图像的超声诊断装置和产生超声图像的方法
JP2014226186A (ja) 超音波測定装置及び超音波測定装置の制御方法
JP5681755B2 (ja) 超音波診断装置及びその作動方法
US20140073920A1 (en) Ultrasound diagnostic apparatus and ultrasound image producing method
JP2014180490A (ja) 超音波測定装置、超音波画像装置及び超音波測定方法

Legal Events

Date Code Title Description
RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20160623

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20171003

TRDD Decision of grant or rejection written
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180531

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20180605

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20180618

R150 Certificate of patent or registration of utility model

Ref document number: 6365243

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150