JP6963677B2 - 超音波診断装置および超音波診断装置の制御方法 - Google Patents
超音波診断装置および超音波診断装置の制御方法 Download PDFInfo
- Publication number
- JP6963677B2 JP6963677B2 JP2020509206A JP2020509206A JP6963677B2 JP 6963677 B2 JP6963677 B2 JP 6963677B2 JP 2020509206 A JP2020509206 A JP 2020509206A JP 2020509206 A JP2020509206 A JP 2020509206A JP 6963677 B2 JP6963677 B2 JP 6963677B2
- Authority
- JP
- Japan
- Prior art keywords
- signal
- ultrasonic
- bubble
- likeness
- unit
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating 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/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/54—Control of the diagnostic device
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
- A61B8/4488—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer the transducer being a phased array
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/461—Displaying means of special interest
- A61B8/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/481—Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating 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/22—Details, e.g. general constructional or apparatus details
- G01N29/24—Probes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating 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/34—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating 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/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8959—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes
- G01S15/8963—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes using pulse inversion
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details 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/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52038—Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
- G01S7/52039—Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target exploiting the non-linear response of a contrast enhancer, e.g. a contrast agent
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details 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/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52071—Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details 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/52085—Details related to the ultrasound signal acquisition, e.g. scan sequences
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
- A61B8/14—Echo-tomography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/024—Mixtures
- G01N2291/02475—Tissue characterisation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/044—Internal reflections (echoes), e.g. on walls or defects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/10—Number of transducers
- G01N2291/106—Number of transducers one or more transducer arrays
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- General Physics & Mathematics (AREA)
- Public Health (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Computer Networks & Wireless Communication (AREA)
- Acoustics & Sound (AREA)
- Hematology (AREA)
- Signal Processing (AREA)
- Nonlinear Science (AREA)
- Gynecology & Obstetrics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Description
本発明は、超音波診断装置および超音波診断装置の制御方法に係り、特に、ハーモニックイメージング法により超音波画像の生成を行う超音波診断装置および超音波診断装置の制御方法に関する。
医用超音波診断装置において、被検体に造影剤を導入して診断を行なう場合等には、特許文献1および2に開示されるように、造影剤の非線形性を利用し、振動子アレイにより受信した超音波エコーから非線形成分を抽出して画像化を行う、いわゆるハーモニックイメージング法が知られている。ハーモニックイメージング法を用いることにより、被検体の組織と造影剤のバブルとのコントラストが高い画像を生成することができる。
ハーモニックイメージング法において、超音波エコーから非線形成分を抽出する方法として、例えば、同一の走査線上に、互いに位相を反転させた第1の超音波パルスと第2の超音波パルスを順次送信し、第1の超音波パルスによる受信信号と第2の超音波パルスによる受信信号を加算するパルスインバージョン法がある。このパルスインバージョン法により、超音波エコーに基づく受信信号から、第1の超音波パルスおよび第2の超音波パルスを形成する基本波の帯域を有する基本波信号を除去して、造影剤のバブルに起因する非線形信号を抽出することができる。
しかしながら、通常、パルスインバージョン法を用いたとしても、組織の境界等からの超音波エコーに起因する信号は、造影剤のバブルに起因する非線形信号と同程度に残存することが多く、非線形信号との区別をすることが難しいという問題があった。従来の技術では、造影剤のバブルが血管内を移動することに伴って生じる、信号の輝度の時間変化を観察することにより、組織の境界等からの超音波エコーに起因する信号と造影剤のバブルに起因する非線形信号とを区別することが多く、取得された信号のうち、どの信号が造影剤のバブルに起因する非線形信号であるかを判断するために、時間と手間を要していた。
また、被検体内において拍動等の組織の動きがある場合には、パルスインバージョン法により相殺されずに残存する基本波信号が多くなり、基本波信号と造影剤のバブルに起因する非線形信号を区別することが難しいことがあるという問題があった。
本発明は、このような従来の問題点を解消するためになされたものであり、被検体の組織に起因する信号と造影剤のバブルに起因する信号を容易に且つ短時間に区別することができる超音波診断装置および超音波診断装置の制御方法を提供することを目的とする。
上記目的を達成するために、本発明の超音波診断装置は、振動子アレイと、振動子アレイから被検体内に互いに位相を反転させた第1の超音波パルスと第2の超音波パルスの組を同一の走査線上に少なくとも2回以上のN回送信させる送信部と、被検体内において発生した超音波エコーを受けた振動子アレイから出力される信号により受信信号を取得する受信部と、受信部により取得された受信信号に対して定められた帯域で直交検波を行うことにより第1の超音波パルスに対応するIQ信号列と第2の超音波パルスに対応するIQ信号列を取得する直交検波部と、直交検波部により取得されたIQ信号列を用いて、第1の超音波パルスに対応するIQ信号と第2の超音波パルスに対応するIQ信号を加算することにより基本波成分が除去された画像信号を取得するパルスインバージョン加算部と、
直交検波部により取得されたIQ信号列に基づいてバブル信号らしさを算出するバブル信号らしさ算出部と、バブル信号らしさ算出部により算出されたバブル信号らしさとパルスインバージョン加算部により取得された画像信号とに基づいて超音波画像を生成する画像生成部とを備えたことを特徴とする。
直交検波部により取得されたIQ信号列に基づいてバブル信号らしさを算出するバブル信号らしさ算出部と、バブル信号らしさ算出部により算出されたバブル信号らしさとパルスインバージョン加算部により取得された画像信号とに基づいて超音波画像を生成する画像生成部とを備えたことを特徴とする。
バブル信号らしさ算出部は、直交検波部により取得されたIQ信号列から自己相関を求め、求められた前記自己相関に基づいてバブル信号らしさを算出することができる。
もしくは、バブル信号らしさ算出部は、直交検波部により取得されたIQ信号列から位相差の分散値を算出し、算出された位相差の分散値を用いてバブル信号らしさを算出することもできる。
もしくは、バブル信号らしさ算出部は、直交検波部により取得されたIQ信号列から振幅の分散値を算出し、算出された振幅の分散値を用いてバブル信号らしさを算出することができる。
もしくは、バブル信号らしさ算出部は、直交検波部により取得されたIQ信号列から位相差の分散値を算出し、算出された位相差の分散値を用いてバブル信号らしさを算出することもできる。
もしくは、バブル信号らしさ算出部は、直交検波部により取得されたIQ信号列から振幅の分散値を算出し、算出された振幅の分散値を用いてバブル信号らしさを算出することができる。
パルスインバージョン加算部により取得された画像信号から非線形信号のパワーおよび速度の少なくとも一方を算出する非線形信号情報算出部をさらに備えることができる。
この際に、画像生成部は、非線形信号情報算出部により算出された非線形信号のパワーおよび速度の少なくとも一方に基づいて超音波画像を生成することができる。
さらに、画像生成部は、非線形信号情報算出部により算出された非線形信号のパワーおよび速度の少なくとも一方とバブル信号らしさ算出部により算出されたバブル信号らしさとに基づいたカラーマップにより超音波画像を生成することもできる。
また、超音波画像を表示する表示部をさらに備えることができる。
この際に、画像生成部は、非線形信号情報算出部により算出された非線形信号のパワーおよび速度の少なくとも一方に基づいて超音波画像を生成することができる。
さらに、画像生成部は、非線形信号情報算出部により算出された非線形信号のパワーおよび速度の少なくとも一方とバブル信号らしさ算出部により算出されたバブル信号らしさとに基づいたカラーマップにより超音波画像を生成することもできる。
また、超音波画像を表示する表示部をさらに備えることができる。
本発明の超音波診断装置の制御方法は、振動子アレイから被検体内に互いに位相を反転させた第1の超音波パルスと第2の超音波パルスの組を同一の走査線上に少なくとも2回以上のN回送信させ、被検体内において発生した超音波エコーを受けた振動子アレイから出力される信号により受信信号を取得し、取得された受信信号に対して定められた帯域で直交検波を行うことにより第1の超音波パルスに対応するIQ信号列と第2の超音波パルスに対応するIQ信号列を取得し、取得されたIQ信号列を用いて第1の超音波パルスに対応するIQ信号と第2の超音波パルスに対応するIQ信号を加算することにより基本波成分が除去された画像信号を取得し、取得されたIQ信号列に基づいてバブル信号らしさを算出し、算出されたバブル信号らしさと、取得された画像信号とに基づいて超音波画像を生成することを特徴とする。
本発明によれば、直交検波部により取得されたIQ信号列のランダム性を表す指標をバブル信号らしさとして算出するバブル信号らしさ算出部を備えるため、被検体の組織に起因する信号と造影剤のバブルに起因する信号を容易に且つ短時間に区別することができる。
以下、この発明の実施の形態を添付図面に基づいて説明する。また、以下においては、被検体に造影剤が導入されているものとする。
実施の形態
図1に、本発明の実施の形態に係る超音波診断装置1の構成を示す。図1に示すように、超音波診断装置1は、振動子アレイ2を備えており、振動子アレイ2に送信部3および受信部4がそれぞれ接続されている。受信部4には、直交検波部5が接続され、直交検波部5に、バブル信号らしさ算出部6およびパルスインバージョン加算部7が接続されている。また、パルスインバージョン加算部7には、非線形信号情報算出部8、画像生成部9、表示制御部10および表示部11が順次接続されている。さらに、バブル信号らしさ算出部6は、画像生成部9に接続されている。
実施の形態
図1に、本発明の実施の形態に係る超音波診断装置1の構成を示す。図1に示すように、超音波診断装置1は、振動子アレイ2を備えており、振動子アレイ2に送信部3および受信部4がそれぞれ接続されている。受信部4には、直交検波部5が接続され、直交検波部5に、バブル信号らしさ算出部6およびパルスインバージョン加算部7が接続されている。また、パルスインバージョン加算部7には、非線形信号情報算出部8、画像生成部9、表示制御部10および表示部11が順次接続されている。さらに、バブル信号らしさ算出部6は、画像生成部9に接続されている。
さらに、送信部3、受信部4、直交検波部5、バブル信号らしさ算出部6、パルスインバージョン加算部7、非線形信号情報算出部8、画像生成部9および表示制御部10に、装置制御部12が接続されており、装置制御部12に、格納部13および操作部14が接続されている。装置制御部12と格納部13は、互いに双方向の情報の受け渡しが可能に接続されている。
また、送信部3、受信部4、直交検波部5、バブル信号らしさ算出部6、位相補正部7、パルスインバージョン加算部7、非線形信号情報算出部8、画像生成部9、表示制御部10および装置制御部12により、プロセッサ15が構成されている。
また、送信部3、受信部4、直交検波部5、バブル信号らしさ算出部6、位相補正部7、パルスインバージョン加算部7、非線形信号情報算出部8、画像生成部9、表示制御部10および装置制御部12により、プロセッサ15が構成されている。
図1に示す超音波診断装置1の振動子アレイ2は、1次元または2次元に配列された複数の振動子を有している。これらの振動子は、それぞれ送信部3から供給される駆動信号に従って超音波を送信すると共に、被検体からの超音波エコーを受信して、超音波エコーに基づく信号を出力する。各振動子は、例えば、PZT(Lead Zirconate Titanate:チタン酸ジルコン酸鉛)に代表される圧電セラミック、PVDF(Poly Vinylidene Di Fluoride:ポリフッ化ビニリデン)に代表される高分子圧電素子およびPMN−PT(Lead Magnesium Niobate-Lead Titanate:マグネシウムニオブ酸鉛−チタン酸鉛固溶体)に代表される圧電単結晶等からなる圧電体の両端に電極を形成することにより構成される。
プロセッサ15の送信部3は、例えば、複数のパルス発生器を含んでおり、装置制御部12からの制御信号に応じて選択された送信遅延パターンに基づいて、振動子アレイ2の複数の振動子から送信される超音波が超音波ビームを形成するように、それぞれの駆動電圧を、遅延量を調節して複数の振動子に供給する。このように、振動子アレイ2の複数の振動子の電極にパルス状の駆動電圧が印加されると、圧電体が伸縮し、それぞれの振動子からパルス状の超音波が発生して、それらの超音波の合成波から、パルス状の超音波ビームすなわち超音波パルスが形成される。送信部3は、このようにして、互いに位相を反転させた第1の超音波パルスおよび第2の超音波パルスを振動子アレイ2から順次発生させ、振動子アレイ2を介して、第1の超音波パルスと第2の超音波パルスの組を同一の走査線に沿って被検体内に複数回送信する。
被検体内に送信された第1の超音波パルスおよび第2の超音波パルスは、例えば、被検体の部位等の対象において反射され、いわゆる超音波エコーとして振動子アレイ2に向かって被検体内を伝搬する。このように振動子アレイ2に向かって伝搬する超音波エコーは、振動子アレイ2を構成するそれぞれの振動子により受信される。この際に、振動子アレイ2を構成するそれぞれの振動子は、伝搬する超音波エコーを受信することにより伸縮して電気信号を発生させ、これらの電気信号を受信部4に出力する。
プロセッサ15の受信部4は、装置制御部12からの制御信号に従って、振動子アレイ2から出力される信号の処理を行う。図2に示すように、受信部4は、増幅部16、AD(Analog Digital)変換部17およびビームフォーマ18が直列接続された構成を有している。
受信部4の増幅部16は、振動子アレイ2を構成するそれぞれの振動子から入力された信号を増幅し、増幅した信号をAD変換部17に送信する。AD変換部17は、増幅部16から送信された信号をデジタルデータに変換し、これらのデータをビームフォーマ18に送信する。ビームフォーマ18は、装置制御部12からの制御信号に応じて選択された受信遅延パターンに基づいて設定される音速または音速の分布に従い、AD変換部17により変換された各データに対してそれぞれの遅延を与えて加算することにより、いわゆる受信フォーカス処理を行う。この受信フォーカス処理により、AD変換部17により変換された各データが整相加算され且つ超音波エコーの焦点が絞り込まれた受信信号が取得される。
ここで、被検体内を伝搬する超音波エコーには、第1の超音波パルスおよび第2の超音波パルスを形成する基本波の帯域を有する基本波成分と、被検体の組織の動きに起因する2次高調波成分と、被検体に導入された造影剤のバブルの微小な振動に起因する非線形成分が含まれている。そのため、受信部4により取得された受信信号には、例えば図3および図4に示すように、超音波エコーの基本波成分に基づく基本波信号E1またはE4、2次高調波成分に基づく2次高調波信号E2および造影剤のバブルによる非線形成分に基づくバブル信号E3を含んでいる。図3に示す例では、正の位相を有する超音波パルスに基づく受信信号が示されており、基本波信号E1は、正の値を有している。一方、図4に示す例では、負の位相を有する超音波パルスに基づく受信信号が示されており、基本波信号E4は、負の値を有している。
プロセッサ15の直交検波部5は、受信部4により取得された受信信号に参照周波数のキャリア信号を混合することにより、受信信号を直交検波して複素データであるIQ信号に変換し、第1の超音波パルスに対応するIQ信号列と第2の超音波パルスに対応するIQ信号列を取得する。この際に、直交検波部5は、被検体に導入されている造影剤の検出精度を向上させるため、図3および図4に示すように、基本波信号E1の信号強度および2次高調波信号E2の信号強度と比較して、造影剤のバブルに起因する非線形な信号からなるバブル信号E3の信号強度が相対的に大きくなる周波数を含むように、直交検波の帯域FBを設定することが望ましい。さらに、基本波信号E1の信号強度と、2次高調波信号E2およびバブル信号E3を含む非線形信号の信号強度とが、互いに比較的近い値を有するように、図3および図4に示すように、バブル信号E3が相対的に大きくなる周波数に加えて、基本波信号E1またはE4の周波数帯域の一部を含むように、直交検波の帯域FBを設定することが好ましい。
プロセッサ15のバブル信号らしさ算出部6は、直交検波部5により取得されたIQ信号列に基づいて算出された自己相関または分散値を用いた指標を、被検体に導入された造影剤のバブルに基づくバブル信号らしさとして算出する。バブル信号らしさ算出部6によるバブル信号らしさの算出については、後に詳しく説明する。
プロセッサ15のパルスインバージョン加算部7は、直交検波部5により取得されたIQ信号列を用いて、第1の超音波パルスに対応するIQ信号と、第2の超音波パルスに対応するIQ信号とを加算することにより、図3および図4に示すような基本波信号E1およびE4が除去された加算信号を取得する。
プロセッサ15の非線形信号情報算出部8は、パルスインバージョン加算部7により取得された加算信号から、被検体の組織に起因する2次高調波信号E2および被検体に導入された造影剤のバブルに起因するバブル信号E3を含む非線形信号のパワーおよび速度の少なくとも一方を非線形信号情報として算出する。
プロセッサ15の画像生成部9は、非線形信号情報算出部8により算出された非線形信号のパワーおよび速度の少なくとも一方と、バブル信号らしさ算出部6により算出されたバブル信号らしさとに基づいて超音波画像を生成する。
プロセッサ15の表示制御部10は、装置制御部12の制御の下、画像生成部9により生成された超音波画像等に所定の処理を施して、表示部11に超音波画像等を表示させる。
プロセッサ15の画像生成部9は、非線形信号情報算出部8により算出された非線形信号のパワーおよび速度の少なくとも一方と、バブル信号らしさ算出部6により算出されたバブル信号らしさとに基づいて超音波画像を生成する。
プロセッサ15の表示制御部10は、装置制御部12の制御の下、画像生成部9により生成された超音波画像等に所定の処理を施して、表示部11に超音波画像等を表示させる。
超音波診断装置1の表示部11は、表示制御部10の制御の下、画像等を表示するものであり、例えば、LCD(Liquid Crystal Display:液晶ディスプレイ)等のディスプレイ装置を含む。
超音波診断装置1の操作部14は、ユーザが入力操作を行うためのものであり、キーボード、マウス、トラックボール、タッチパッドおよびタッチパネル等を備えて構成することができる。
超音波診断装置1の操作部14は、ユーザが入力操作を行うためのものであり、キーボード、マウス、トラックボール、タッチパッドおよびタッチパネル等を備えて構成することができる。
格納部13は、超音波診断装置1の動作プログラム等を格納するもので、フラッシュメモリ、HDD(Hard Disc Drive:ハードディスクドライブ)、SSD(Solid State Drive:ソリッドステートドライブ)、FD(Flexible Disc:フレキシブルディスク)、MOディスク(Magneto-Optical disc:光磁気ディスク)、MT(Magnetic Tape:磁気テープ)、RAM(Random Access Memory:ランダムアクセスメモリ)、CD(Compact Disc:コンパクトディスク)、DVD(Digital Versatile Disc:デジタルバーサタイルディスク)、SDカード(Secure Digital card:セキュアデジタルカード)、USBメモリ(Universal Serial Bus memory:ユニバーサルシリアルバスメモリ)等の記録メディア、またはサーバ等を用いることができる。
なお、送信部3、受信部4、直交検波部5、バブル信号らしさ算出部6、パルスインバージョン加算部7、非線形信号情報算出部8、画像生成部9、表示制御部10および装置制御部12を有するプロセッサ15は、CPU(Central Processing Unit:中央処理装置)、および、CPUに各種の処理を行わせるための制御プログラムから構成されるが、FPGA(Field Programmable Gate Array:フィードプログラマブルゲートアレイ)、DSP(Digital Signal Processor:デジタルシグナルプロセッサ)、ASIC(Application Specific Integrated Circuit:アプリケーションスペシフィックインテグレイテッドサーキット)、GPU(Graphics Processing Unit:グラフィックスプロセッシングユニット)、その他のIC(Integrated Circuit:集積回路)を用いて構成されてもよい。また、これらの送信部3、受信部4、直交検波部5、バブル信号らしさ算出部6、パルスインバージョン加算部7、非線形信号情報算出部8、画像生成部9、表示制御部10および装置制御部12を部分的にあるいは全体的に1つのCPU等に統合させて構成することもできる。
次に、図5に示すフローチャートを用いて、実施の形態における超音波診断装置1の動作を詳細に説明する。実施の形態において、超音波診断装置1は、同一の走査線上に互いに位相を反転させた第1の超音波パルスと第2の超音波パルスを順次送信し、第1の超音波パルスによる受信信号と第2の超音波パルスによる受信信号を加算するパルスインバージョン法を用いて、超音波画像を生成する。
まず、ステップS1において、送信部3は、互いに位相を反転させた第1の超音波パルスと第2の超音波パルスを、振動子アレイ2を介して同一の走査線上に複数回送信する。この際に、送信部3は、第1の超音波パルスと第2の超音波パルスの組を同一の走査線上にN回送信した後に、次の走査線上に第1の超音波パルスと第2の超音波パルスの組をN回送信する。ここで、Nは2以上の整数である。例えば、送信部3は、図6に示すように、各走査線L1、L2、L3、L4、L5上において、第1の超音波パルスFPと第2の超音波パルスSPを、交互に4回ずつ送信している。また、図6に示す例においては、時系列に隣り合う第1の超音波パルスFP同士の時間間隔PRT1と、時系列に隣り合う第2の超音波パルスSP同士の時間間隔PRT1は、互いに同一である。
ステップS2において、受信部4は、ステップS1で被検体内に送信された第1の超音波パルスFPと第2の超音波パルスSPに基づいて被検体内で発生した超音波エコーを受信した振動子アレイ2から出力される信号により、受信信号を取得する。
続くステップS3において、直交検波部5は、ステップS2で取得された受信信号に対して定められた帯域FBで直交検波を行うことにより、第1の超音波パルスFPに対応するIQ信号列と第2の超音波パルスに対応するIQ信号列を取得する。この際に、直交検波部5は、例えば、図3および図4に示すように、基本波信号E1の信号強度および2次高調信号E2の信号強度と比較して、造影剤のバブルに起因する非線形な信号からなるバブル信号E3の信号強度が相対的に大きくなるように、且つ、基本波信号E1の信号強度と、2次高調波信号E2およびバブル信号E3を含む非線形信号の信号強度とが互いに近い値を有するように、直交検波の帯域FBを設定し、この帯域FBにおいて直交検波を実行する。
また、第1の超音波パルスFPに対応するIQ信号列と第2の超音波パルスSPに対応するIQ信号列は、互いに極性の異なる位相を有している。例えば、第1の超音波パルスFPが正の位相を有し、第2の超音波パルスSPが負の位相を有している場合に、図7に示すように、第1の超音波パルスFPに対応するIQ信号P1、P3、P5、P7、P9、P11を含むIQ信号列C1は、正の位相を有し、第2の超音波パルスSPに対応するIQ信号P2、P4、P6、P8、P10、P12を含むIQ信号列C2は、負の位相を有する。
続くステップS4において、バブル信号らしさ算出部6は、ステップS3で取得されたIQ信号列C1およびC2を用いて、バブル信号らしさを表す指標を算出する。例えば、バブル信号らしさ算出部6は、ステップS3で得られたIQ信号により下記式(1)に示すように自己相関VTを算出し、下記式(2)に示すようにIQ信号のパワーPTを算出し、さらに、算出された自己相関VTおよびパワーPTを用いて下記式(3)に示すように、バブル信号らしさとして、分散値VSを算出する。
VT=[Σ(Pj+1・P* j)]/(2n−1) (j=1,2,・・・,2n−1)・・・(1)
PT=(Σ|Pk|2)/(2n) (k=1,2,・・・,2n)・・・(2)
VS=1−(|VT|/PT) ・・・(3)
VT=[Σ(Pj+1・P* j)]/(2n−1) (j=1,2,・・・,2n−1)・・・(1)
PT=(Σ|Pk|2)/(2n) (k=1,2,・・・,2n)・・・(2)
VS=1−(|VT|/PT) ・・・(3)
ここで、数式(1)および(2)におけるnは自然数である。また、IQ信号の自己相関とは、時系列に異なる2つのIQ信号のうち、時系列において後のIQ信号と、時系列において前のIQ信号の複素共役との積により計算されるものである。
また、バブル信号らしさの指標値は、信号の大小によらない規格化された値であることが望ましい。規格化の方法については特に限定しないが、ここでは最も単純な方法として、PTで|VT|を除し、バブル信号らしさの指標値が0から1の範囲に値を有するように、指標値を規格化している。後述する他の指標値の例についても、同様の考えに基づき規格化を行っている。
このようにして算出された分散値VSが、バブル信号らしさを表す指標値として用いられることができる理由について説明する。
また、バブル信号らしさの指標値は、信号の大小によらない規格化された値であることが望ましい。規格化の方法については特に限定しないが、ここでは最も単純な方法として、PTで|VT|を除し、バブル信号らしさの指標値が0から1の範囲に値を有するように、指標値を規格化している。後述する他の指標値の例についても、同様の考えに基づき規格化を行っている。
このようにして算出された分散値VSが、バブル信号らしさを表す指標値として用いられることができる理由について説明する。
まず、ステップS3で得られたIQ信号列C1およびC2において、図3および図4に示すような2次高調波信号E2と造影剤のバブルに起因するバブル信号E3を含む非線形信号の影響が比較的小さく、基本波信号E1およびE4が支配的である場合には、IQ信号列C1およびC2は、例えば、図8に示すようなIQ信号により構成される。ここで、図8には、超音波エコーにおいて造影剤のバブルの影響が無く、被検体の組織に起因する成分が支配的な場合のIQ信号が示されており、IQ信号の例として、3つのIQ信号P1、P2、P3が示されている。また、実際には、IQ信号P1とP3は互いに異なる大きさおよび位相を有する信号であるが、説明のために、IQ信号P1およびP3は互いに等しいとする。
図8に示すように、IQ信号P1およびP3は、例えば、基本波信号E1に対応する基本波ベクトルG1と、被検体の組織に起因する2次高調波信号E2等の高次の高調波信号に対応する非線形信号ベクトルH1の和で表され、IQ信号P2は、基本波信号E4に対応する基本波ベクトルG2と高次の高調波信号に対応する非線形信号ベクトルH2との和で表される。ここで、高次の高調波信号の信号強度は、基本波信号E1の信号強度に対して、例えば10分の1以下のオーダーを有する等、非常に小さく、非線形信号ベクトルH1およびH2の大きさは、基本波ベクトルG1およびG2の大きさと比較して非常に小さい。
そのため、IQ信号P1およびP3は、概ね基本波ベクトルG1に等しく、IQ信号P2は、概ね基本波ベクトルG2に等しい。このように、ステップS3で取得されたIQ信号列C1およびC2において基本波信号E1およびE4が支配的である場合には、IQ信号列C1およびC2のランダム性が低い。また、この際に、IQ信号P1およびP3は、IQ信号P2の位相からIQ信号P1の位相を減じた位相差DA1と、IQ信号P3の位相からIQ信号P2の位相を減じた位相差DA2は、概ね180度近傍の値となる。
そのため、例えば、下記式(4)に示すように、数式(1)中に現れる時系列に隣り合うIQ信号同士の自己相関を自己相関ベクトルVTjとすると、自己相関ベクトルVTjは、図9に示すように、概ね1方向に向かって延びるベクトルとなる。ここで、図9においては、IQ信号P1およびP2に基づいて算出された自己相関ベクトルVT1と、IQ信号P2およびP3に基づいて算出された自己相関ベクトルVT2が示されており、自己相関ベクトルVT1およびVT2は、それぞれ180度近傍の位相DA1およびDA2を有している。
VTj=Pj+1・P* j (j=1,2,・・・,2n−1)・・・(4)
VTj=Pj+1・P* j (j=1,2,・・・,2n−1)・・・(4)
このような場合には、数式(1)に示される自己相関VTの絶対値と数式(2)に示されるパワーPTの絶対値との比が、概ね1に等しくなるため、数式(3)に示される分散VSは、ゼロの近傍の値となる。
次に、ステップS3で得られたIQ信号列C1およびC2において、造影剤のバブルに起因する非線形な信号からなるバブル信号E3の影響が比較的大きい場合には、IQ信号列C1およびC2は、例えば、図10に示すようなIQ信号により構成される。ここで、造影剤のバブルに起因するバブル信号E3の信号強度は、組織に起因する2次高調波信号E2の信号強度と比べて比較的強いため、図3および図4に示すような帯域FBにおいて直交検波がなされることにより、例えば、基本波信号E1および2次高調波信号E2の信号強度に対して同程度もしくはそれ以上の信号強度を有するバブル信号E3を含む、IQ信号が得られる。このようなIQ信号の例として、図10には、IQ信号P1、P2、P3が示されている。図8と同様に、実際には、IQ信号P1とP3は互いに異なる大きさおよび位相を有する信号であるが、説明のために、IQ信号P1およびP3は互いに等しいとする。
図10に示すように、IQ信号P1およびP3は、例えば、基本波信号E1に対応する基本波ベクトルG1と、2次高調波信号E2およびバブル信号E3を含む非線形信号に対応する非線形信号ベクトルH1の和で表され、IQ信号P2は、基本波信号E4に対応する基本波ベクトルG2と2次高調波信号E2およびバブル信号E3を含む非線形信号に対応する非線形信号ベクトルH2との和で表される。図10に示す例では、基本波ベクトルG1の大きさと非線形信号ベクトルH1の大きさとの比、および、基本波ベクトルG2の大きさと非線形信号ベクトルH2の大きさとの比は、概ね1に等しい。
このような場合には、IQ信号P1およびP3は、基本波ベクトルG1から大きくずれた信号であり、IQ信号P2は、基本波ベクトルG2から大きくずれた信号である。このように、ステップS3で取得されたIQ信号列C1およびC2において基本波信号E1およびE4が支配的ではなく、非線形信号が比較的大きい場合には、IQ信号のランダム性が高い。また、この際に、IQ信号P2の位相からIQ信号P1の位相を減じた位相差DA1と、IQ信号P3の位相からIQ信号P1の位相を減じた位相差DA2は、例えば、90度近傍の値または270度近傍の値等の、180度から離れた値となる。
そのため、数式(4)に示す自己相関ベクトルVTjは、例えば図11に示すように、概ね反対方向に向かって延びるベクトルとなる。ここで、図11においては、IQ信号P1およびP2に基づいて算出された自己相関ベクトルVT1と、IQ信号P2およびP3に基づいて算出された自己相関ベクトルVT2が示されており、自己相関ベクトルVT1は、90度近傍の位相を有しており、VT2は、270度近傍の位相を有している。
このような場合には、数式(1)に示される自己相関VTの計算において、数式(4)に示される自己相関ベクトルVTjが概ね相殺されるため、数式(3)に示される分散VSは、1の近傍の値となる。
したがって、分散VSの値が1に近づくほど、ステップS3で取得されたIQ信号列C1およびC2のランダム性が高くなり、すなわち、IQ信号列C1およびC2において造影剤のバブルに起因するバブル信号E3の影響が大きくなるため、IQ信号列C1およびC2が造影剤のバブルに起因する信号らしいことがわかる。また、分散VSの値がゼロに近づくほど、IQ信号列C1およびC2のランダム性が低くなり、すなわち、IQ信号列C1およびC2において造影剤のバブルに起因するバブル信号E3の影響が小さくなり、基本波信号E1およびE4の影響が大きくなるため、IQ信号列C1およびC2が被検体の組織に起因する信号らしいことがわかる。
したがって、分散VSの値が1に近づくほど、ステップS3で取得されたIQ信号列C1およびC2のランダム性が高くなり、すなわち、IQ信号列C1およびC2において造影剤のバブルに起因するバブル信号E3の影響が大きくなるため、IQ信号列C1およびC2が造影剤のバブルに起因する信号らしいことがわかる。また、分散VSの値がゼロに近づくほど、IQ信号列C1およびC2のランダム性が低くなり、すなわち、IQ信号列C1およびC2において造影剤のバブルに起因するバブル信号E3の影響が小さくなり、基本波信号E1およびE4の影響が大きくなるため、IQ信号列C1およびC2が被検体の組織に起因する信号らしいことがわかる。
ここで、自己相関VTは、時間的に互いに隣接するIQ信号対の一方と、他方の複素共役との積により求められるが、これらのIQ信号に含まれるバブル信号E3の信号強度が基本波信号E1および2次高調波信号E2の信号に対して同程度もしくはそれ以上である場合には、このような積の計算において、複素共役となる、第1の超音波パルスFPからの超音波エコーに対応する第1のIQ信号と第2の超音波パルスSPからの超音波エコーに対応する第2のIQ信号とが交互に入れ替わることにより、算出される位相差DAjが180度から交互に離れた値を有するという性質が得られている。このような場合には、つまり、自己相関VTを算出する際に、積が計算される第1のIQ信号と第2のIQ信号の順番が周期的に入れ替われば同様の性質が得られるため、IQ信号Pj+3とPjおよびPj+5とPj等から算出される自己相関でも、同様の結果を得ることができる。
このようにして、ステップS4でバブル信号らしさが算出されると、続くステップS5において、パルスインバージョン加算部7は、ステップS3で取得されたIQ信号列C1およびC2を用いて、第1の超音波パルスFPに対応するIQ信号と第2の超音波パルスSPに対応するIQ信号を加算することにより基本波信号E1およびE4が除去された加算信号を取得する。
例えば、パルスインバージョン加算部7は、図12に示すように、互いに時系列に隣り合い且つ互いに極性の異なるIQ信号を加算することにより、加算信号a1,a2,a3,a4,a5,a6およびb1,b2,b3,b4,b5を算出する。この際に、パルスインバージョン加算部7は、例えば、下記式(5)を用いて加算信号a1,a2,a3,a4,a5,a6を算出し、下記式(6)を用いて加算信号b1,b2,b3,b4,b5を算出する。
am=P2m−1+P2m (m=1,2,・・・,n)・・・(5)
bq=P2q+P2q+1 (q=1,2,・・・,n−1)・・・(6)
am=P2m−1+P2m (m=1,2,・・・,n)・・・(5)
bq=P2q+P2q+1 (q=1,2,・・・,n−1)・・・(6)
続くステップS6において、非線形信号情報算出部8は、ステップS5で算出された加算信号amおよびbqを用いて、被検体の組織に起因する2次高調波信号E2および被検体に導入された造影剤のバブルに起因する非線形な信号からなるバブル信号E3を含む非線形信号のパワーおよび速度ベクトルの少なくとも一方を算出する。例えば、非線形信号情報算出部8は、下記式(7)を用いてバブル信号E3のパワーPBを算出し、下記式(8)を用いてバブル信号E3の速度VBを算出することができる。
PB=[Σ|am|2+Σ|bq|2]/(2n−1)
(m=1,2,・・・,n、 q=1,2,・・・,n−1)・・・(7)
VB=[Σ(ar+1・a* r)+Σ(bt+1・b* t)]/(2n−3)
(r=1,2,・・・,n−2、 t=1,2,・・・,n−3)・・・(8)
PB=[Σ|am|2+Σ|bq|2]/(2n−1)
(m=1,2,・・・,n、 q=1,2,・・・,n−1)・・・(7)
VB=[Σ(ar+1・a* r)+Σ(bt+1・b* t)]/(2n−3)
(r=1,2,・・・,n−2、 t=1,2,・・・,n−3)・・・(8)
続くステップS7において、画像生成部9は、ステップS6で算出された非線形信号のパワーPBおよび速度VBのうち少なくとも一方と、ステップS4で算出されたバブル信号らしさに基づいて超音波画像を生成し、生成した超音波画像を表示部11に表示する。例えば、画像生成部9は、図13に示すように、ステップS6で算出された非線形信号のパワーPBの値を明度変化により表し、バブル信号らしさの指標値を色Bの彩度変化により表す、いわゆるカラーマップにより超音波画像を生成し、生成した超音波画像を表示部11に表示することができる。図13に示す例では、パワーPBの値が大きくなるほど明度が大きくなり、バブル信号らしさの指標値が大きくなるほど色Bの彩度が高くなっている。
このようにして、本発明の実施の形態に係る超音波診断装置1の動作が終了する。
このようにして、本発明の実施の形態に係る超音波診断装置1の動作が終了する。
以上から、本発明の実施の形態1に係る超音波診断装置1によれば、バブル信号らしさ算出部6が、直交検波部5により取得されたIQ信号列C1およびC2から求めた自己相関VTに基づいてバブル信号らしさを算出するため、被検体の組織に起因する信号と、被検体に導入された造影剤のバブルに起因する信号とを容易に且つ短時間に区別することができる。
さらに、算出されたバブル信号らしさの指標値と、非線形信号のパワーPBおよび速度VBのうち一方とに基づいて超音波画像を生成し、この超音波画像を表示部11に表示するため、ユーザが、被検体に導入された造影剤のバブルに起因する信号を容易に把握することができる。
さらに、算出されたバブル信号らしさの指標値と、非線形信号のパワーPBおよび速度VBのうち一方とに基づいて超音波画像を生成し、この超音波画像を表示部11に表示するため、ユーザが、被検体に導入された造影剤のバブルに起因する信号を容易に把握することができる。
なお、実施の形態において、バブル信号らしさ算出部6は、数式(3)に示す分散VSをバブル信号らしさの指標値として算出しているが、IQ信号列に基づいて算出された自己相関VTを用いた指標値であれば、これに限定されない。例えば、バブル信号らしさ算出部6は、下記式(9)に示すように、数式(3)中に現れる|VT|/PTを用いて分散値VXを計算することにより、分散値VXをバブル信号らしさとして算出することができる。この場合には、分散値VXが1に近づくほど、直交検波部5により取得されたIQ信号列C1およびC2が被検体の組織に起因する信号らしいと判断することができ、分散値VXがゼロに近づくほど、IQ信号列C1およびC2が造影剤のバブルに起因する信号らしいと判断することができる。
VX=|VT|/PT ・・・(9)
VX=|VT|/PT ・・・(9)
また、例えば、バブル信号らしさ算出部6は、直交検波部5により取得されたIQ信号列C1およびC2において、互いに時系列に隣り合うIQ信号の位相差の分散値を、バブル信号らしさとして算出することもできる。例えば、下記式(10)に示すように、IQ信号Pkの実部をXk、虚部をYkとすると、互いに時系列に隣り合うIQ信号の位相差DAjは、下記式(11)で表される。ここで、下記式(10)におけるiは、虚数単位を表す。
Pk=Xk+iYk (k=1,2,・・・,2n)・・・(10)
DAj=tan−1[(Yj+1Xj−Xj+1Yj)/(Xj+1Xj+Yj+1Yj)] (j=1,2,・・・,2n−1)・・・(11)
Pk=Xk+iYk (k=1,2,・・・,2n)・・・(10)
DAj=tan−1[(Yj+1Xj−Xj+1Yj)/(Xj+1Xj+Yj+1Yj)] (j=1,2,・・・,2n−1)・・・(11)
さらに、バブル信号らしさ算出部6は、下記式(12)に示すように、数式(11)により算出された位相差DAjの分散値VS1をバブル信号らしさの指標値として算出することができる。ここで、下記式(12)におけるE(DAj)は、位相差DAjの算術平均である。図9に示すように、非線形な信号からなるバブル信号E3の影響が小さい場合には、それぞれのDAjは概ね互いに等しくなるため、VS1はゼロに近づく。一方、図11に示すように、非線形な信号からなるバブル信号E3の影響が大きい場合には、DAjとDAj+1は互いに離れた値となるため、VS1はゼロよりも有意に大きな値となる。よって、分散値VS1が1に近づくほど、直交検波部5により取得されたIQ信号列C1およびC2が造影剤のバブルに起因する信号らしいと判断することができ、分散値VS1がゼロに近づくほど、IQ信号列C1およびC2が被検体の組織に起因する信号らしいと判断することができる。また、自己相関VTと同様の考え方により、Pj+3とPjおよびPj+5とPj等から算出される位相差でも同様の結果を得ることができる。
VS1=1−[E(DAj)2/{ΣDAj 2/(2n−1)}] (j=1,2,・・・、2n−1)・・・(12)
VS1=1−[E(DAj)2/{ΣDAj 2/(2n−1)}] (j=1,2,・・・、2n−1)・・・(12)
また、例えば、バブル信号らしさ算出部6は、下記式(13)に示すように、数式(12)中に現れる分散値VX1を、バブル信号らしさとして算出することもできる。この場合には、分散値VX1が1に近づくほど、直交検波部5により取得されたIQ信号列C1およびC2が被検体の組織に起因する信号らしいと判断することができ、分散値VX1がゼロに近づくほど、造影剤のバブルに起因する信号らしいと判断することができる。
VX1=[E(DAj)2/{ΣDAj 2/(2n−1)}] (j=1,2,・・・、2n−1)・・・(13)
VX1=[E(DAj)2/{ΣDAj 2/(2n−1)}] (j=1,2,・・・、2n−1)・・・(13)
また、例えば、バブル信号らしさ算出部6は、下記式(14)に示すように、直交検波部5により取得さえたIQ信号列C1およびC2の振幅の分散値VS2を、バブル信号らしさの指標値として算出することができる。ここで、下記式(14)におけるE(|Pk|)は、IQ信号Pkの絶対値すなわち振幅の算術平均である。図8に示すように、非線形な信号からなるバブル信号E3の影響が小さい場合には、それぞれの|Pk|は概ね互いに等しくなるため、VS2はゼロに近づく。一方、図10に示すように、非線形な信号からなるバブル信号E3の影響が大きい場合には、|Pk|と|Pk+1|は互いに離れた値となるため、VS2はゼロよりも有意に大きな値となる。よって、分散値VS2が1に近づくほど、直交検波部5により取得されたIQ信号列C1およびC2が造影剤のバブルに起因する信号らしいと判断することができ、分散値VS2がゼロに近づくほど、IQ信号列C1およびC2が被検体の組織に起因する信号らしいと判断することができる。
VS2=1−[E(|Pk|)2/{Σ|Pk|2/(2n−1)}] (k=1,2,・・・,2n)・・・(14)
VS2=1−[E(|Pk|)2/{Σ|Pk|2/(2n−1)}] (k=1,2,・・・,2n)・・・(14)
また、例えば、バブル信号らしさ算出部6は、下記式(15)に示すように、数式(14)中に現れる分散値VX2を、バブル信号らしさの指標値として算出することもできる。この場合には、分散値VX2が1に近づくほど、直交検波部5により取得されたIQ信号列C1およびC2が被検体の組織に起因する信号らしいと判断することができ、分散値VX2がゼロに近づくほど、造影剤のバブルに起因する信号らしいと判断することができる。
VX2=[E(|Pk|)2/{Σ|Pk|2/(2n−1)}] (k=1,2,・・・,2n)・・・(15)
VX2=[E(|Pk|)2/{Σ|Pk|2/(2n−1)}] (k=1,2,・・・,2n)・・・(15)
また、実施の形態において画像生成部9により生成される超音波画像の表示例として、図13に示すような、バブル信号E3のパワーPBの値とバブル信号らしさとに基づく超音波画像を示しているが、画像生成部9により生成される超音波画像は、これに限定されない。例えば、画像生成部9は、図13に示す超音波画像の表示例において、バブル信号らしさの指標値が一定の値よりも大きい箇所のみ色Bを用いて表示し、バブル信号らしさの指標値が一定の値よりも小さい箇所については色Bを用いた表示を行わないで、パワーPBの値をグレースケールにより表示することができる。
また、画像生成部9は、例えば、図14に示すように、被検体の組織に起因する2次高調波信号E2および被検体に導入された造影剤のバブルに起因する非線形な信号からなるバブル信号E3を含む非線形信号の速度VBの位相の極性に応じて色B1および色B2の一方を選択し、パワーPBの値を明度変化により表し、バブル信号らしさの指標値を色B1および色B2の彩度変化により表した超音波画像を生成し、生成した超音波画像を表示部11に表示することができる。図14に示す例では、パワーPBの値が大きくなるほど明度が大きくなり、速度VBの位相が正の領域においては、バブル信号らしさの指標値が大きくなるほど色B1の彩度が高くなり、速度VBの位相が負の領域においては、バブル信号らしさの指標値が大きくなるほど色B2の彩度が高くなっている。
また、画像生成部9は、例えば、非線形信号の速度VBの位相の値とバブル信号らしさの指標値とに基づいて超音波画像を生成し、生成した超音波画像を表示部11に表示することもできる。例えば、この際に、画像生成部9は、図14に示す例と同様に、速度VBの位相の極性に応じて色B1および色B2の一方を選択し、速度VBの位相の絶対値を明度変化により表し、バブル信号らしさの指標値を色B1および色B2の彩度変化により表した超音波画像を生成することができる。
また、画像生成部9は、非線形信号のパワーPBおよび速度VBのうち少なくとも一方に基づいて超音波画像を生成し、生成した超音波画像を表示部11に表示することもできる。例えば、画像生成部9は、図示しないが、パワーPBが大きいほど明度が大きくなるようなグレースケールにより超音波画像を生成することができる。
また、画像生成部9は、非線形信号のパワーPBおよび速度VBのうち少なくとも一方に基づいて超音波画像を生成する場合に、バブル信号らしさの指標値を超音波画像に重畳して、または、並べて表示部11に表示することができる。この際に、例えば、画像生成部9は、図示しないが、操作部14を介してユーザにより指定された超音波画像上の位置に対応するバブル信号らしさの指標値を、表示部11に表示することができる。また、この際に、例えば、操作部14を介してユーザが操作することができるカーソルを表示部11に表示させておき、このカーソルにより、操作部14を介してユーザが超音波画像上の位置を指定することができる。
なお、超音波診断装置1に表示部11とは異なるディスプレイを設け、このディスプレイにバブル信号らしさの指標値を表示することもできる。
なお、超音波診断装置1に表示部11とは異なるディスプレイを設け、このディスプレイにバブル信号らしさの指標値を表示することもできる。
また、図示しないが、超音波画像装置1にBモード画像を生成するためのBモード処理部を設けることにより、被検体の断層画像を表すBモード画像上に、非線形信号のパワーPBおよび速度VB、ならびに、バブル信号らしさを画像化して表示部11に重畳表示させることができる。また、被検体の断層画像を表すBモード画像に並べて、非線形信号のパワーPBおよび速度VB、ならびに、バブル信号らしさを画像化して表示部11に表示させることもできる。
また、実施の形態において、パルスインバージョン加算部7は、時系列に隣り合う第1の超音波パルスFPに対応するIQ信号と、第2の超音波パルスSPに対応するIQ信号を加算しているが、時系列に隣り合わない任意の組み合わせにより、第1の超音波パルスFPに対応するIQ信号と、第2の超音波パルスSPに対応するIQ信号とを加算することもできる。しかしながら、時系列に隣り合う第1の超音波パルスFPに対応するIQ信号と、第2の超音波パルスSPに対応するIQ信号を加算する方が、被検体の組織の動きの影響が少ないため、好ましい。
また、パルスインバージョン加算部7は、数式(5)を用いて加算信号amを算出し、数式(6)を用いて加算信号bqを算出しているが、加算信号amおよびbqのうち、一方のみを算出することもできる。
1 超音波診断装置、2 振動子アレイ、3 送信部、4 受信部、5 直交検波部、6 バブル信号らしさ算出部、7 パルスインバージョン加算部、8 非線形信号情報算出部、9 画像生成部、10 表示制御部、11 表示部、12 装置制御部、13 格納部、14 操作部、15 プロセッサ、16 増幅部、17 AD変換部、18 ビームフォーマ、a1,a2,a3,a4,a5,a6,b1,b2,b3,b4,b5 加算信号、B,B1,B2 色、C1,C2 IQ信号列、DA1,DA2 位相差、E1,E4 基本波信号、E2 2次高調波信号、E3 バブル信号、FB 帯域、FP 第1の超音波パルス、G1,G2 基本波ベクトル、H1,H2 非線形信号ベクトル、VT1,VT2 自己相関ベクトル、L1,L2,L3,L4,L5 走査線、P1,P2,P3,P4,P5,P6,P7,P8,P9,P10,P11,P12 IQ信号、PRT1 時間間隔、SP 第2の超音波パルス。
Claims (9)
- 振動子アレイと、
前記振動子アレイから被検体内に互いに位相を反転させた第1の超音波パルスと第2の超音波パルスの組を同一の走査線上に少なくとも2回以上のN回送信させる送信部と、
前記被検体内において発生した超音波エコーを受けた前記振動子アレイから出力される信号により受信信号を取得する受信部と、
前記受信部により取得された前記受信信号に対して定められた帯域で直交検波を行うことにより前記第1の超音波パルスに対応するIQ信号列と前記第2の超音波パルスに対応するIQ信号列を取得する直交検波部と、
前記直交検波部により取得された前記IQ信号列を用いて、前記第1の超音波パルスに対応するIQ信号と前記第2の超音波パルスに対応するIQ信号を加算することにより基本波成分が除去された画像信号を取得するパルスインバージョン加算部と、
前記直交検波部により取得された前記IQ信号列に基づいてバブル信号らしさを算出するバブル信号らしさ算出部と、
前記バブル信号らしさ算出部により算出された前記バブル信号らしさと前記パルスインバージョン加算部により取得された前記画像信号とに基づいて超音波画像を生成する画像生成部と
を備えた超音波診断装置。 - 前記バブル信号らしさ算出部は、前記直交検波部により取得された前記IQ信号列から自己相関を求め、求められた前記自己相関に基づいて前記バブル信号らしさを算出する請求項1に記載の超音波診断装置。
- 前記バブル信号らしさ算出部は、前記直交検波部により取得された前記IQ信号列から位相差の分散値を算出し、算出された前記位相差の分散値を用いて前記バブル信号らしさを算出する請求項1に記載の超音波診断装置。
- 前記バブル信号らしさ算出部は、前記直交検波部により取得された前記IQ信号列から振幅の分散値を算出し、算出された前記振幅の分散値を用いて前記バブル信号らしさを算出する請求項1に記載の超音波診断装置。
- 前記パルスインバージョン加算部により取得された前記画像信号から非線形信号のパワーおよび速度の少なくとも一方を算出する非線形信号情報算出部をさらに備える請求項1〜4のいずれか一項に記載の超音波診断装置。
- 前記画像生成部は、前記非線形信号情報算出部により算出された前記非線形信号のパワーおよび速度の少なくとも一方に基づいて前記超音波画像を生成する請求項5に記載の超音波診断装置。
- 前記画像生成部は、前記非線形信号情報算出部により算出された前記非線形信号のパワーおよび速度の少なくとも一方と前記バブル信号らしさ算出部により算出された前記バブル信号らしさとに基づいたカラーマップにより前記超音波画像を生成する請求項6に記載の超音波診断装置。
- 前記超音波画像を表示する表示部をさらに備える請求項1〜7のいずれか一項に記載の超音波診断装置。
- 振動子アレイから被検体内に互いに位相を反転させた第1の超音波パルスと第2の超音波パルスの組を同一の走査線上に少なくとも2回以上のN回送信させ、
前記被検体内において発生した超音波エコーを受けた前記振動子アレイから出力される信号により受信信号を取得し、
取得された前記受信信号に対して定められた帯域で直交検波を行うことにより前記第1の超音波パルスに対応するIQ信号列と前記第2の超音波パルスに対応するIQ信号列を取得し、
取得された前記IQ信号列を用いて前記第1の超音波パルスに対応するIQ信号と前記第2の超音波パルスに対応するIQ信号を加算することにより基本波成分が除去された画像信号を取得し、
取得された前記IQ信号列に基づいてバブル信号らしさを算出し、
算出された前記バブル信号らしさと、取得された前記画像信号とに基づいて超音波画像を生成する
超音波診断装置の制御方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018067682 | 2018-03-30 | ||
JP2018067682 | 2018-03-30 | ||
PCT/JP2019/013237 WO2019189386A1 (ja) | 2018-03-30 | 2019-03-27 | 超音波診断装置および超音波診断装置の制御方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2019189386A1 JPWO2019189386A1 (ja) | 2021-02-18 |
JP6963677B2 true JP6963677B2 (ja) | 2021-11-10 |
Family
ID=68062073
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020509206A Active JP6963677B2 (ja) | 2018-03-30 | 2019-03-27 | 超音波診断装置および超音波診断装置の制御方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US11969296B2 (ja) |
EP (1) | EP3777699A4 (ja) |
JP (1) | JP6963677B2 (ja) |
WO (1) | WO2019189386A1 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP4312793A1 (en) * | 2021-03-31 | 2024-02-07 | Exo Imaging Inc. | Processing circuitry, system and method to test pixels in an ultrasonic imaging device |
EP4173572A1 (en) | 2021-10-29 | 2023-05-03 | SuperSonic Imagine | Method and system for determining a physical characteristic of a medium |
Family Cites Families (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ATE285711T1 (de) * | 1995-10-10 | 2005-01-15 | Advanced Tech Lab | Ultraschall-bilderzeugung zur diagnostik mittels kontrastmitteln |
US6537222B1 (en) * | 1997-08-26 | 2003-03-25 | Koninklijke Philips Electronics N.V. | Methods for the detection of contrast agents in ultrasonic imaging |
JP4117383B2 (ja) | 2000-03-22 | 2008-07-16 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波撮影装置 |
JP4567842B2 (ja) * | 2000-04-10 | 2010-10-20 | 株式会社東芝 | 超音波診断装置 |
JP4723747B2 (ja) | 2001-04-09 | 2011-07-13 | 株式会社東芝 | 超音波診断装置 |
JP4945040B2 (ja) * | 2001-09-28 | 2012-06-06 | 株式会社東芝 | 超音波診断装置 |
JP3959257B2 (ja) * | 2001-11-08 | 2007-08-15 | 株式会社東芝 | 超音波診断装置 |
JP3908555B2 (ja) | 2002-02-08 | 2007-04-25 | 株式会社東芝 | 超音波診断装置 |
JP2005081073A (ja) * | 2003-09-11 | 2005-03-31 | Toshiba Corp | 超音波診断装置 |
JP4583068B2 (ja) * | 2004-05-11 | 2010-11-17 | 株式会社東芝 | 超音波診断装置 |
WO2006003556A1 (en) * | 2004-06-30 | 2006-01-12 | Koninklijke Philips Electronics, N.V. | Nonlinear ultrasonic diagnostic imaging using intermodulation product signals |
JP4738061B2 (ja) * | 2005-05-30 | 2011-08-03 | 株式会社東芝 | 超音波診断装置及び超音波診断プログラム |
US8315693B2 (en) * | 2006-02-28 | 2012-11-20 | Physio-Control, Inc. | Electrocardiogram monitoring |
US10203274B2 (en) * | 2006-10-06 | 2019-02-12 | California Institute Of Technology | Optical focusing inside scattering media with time-reversed ultrasound microbubble encoded (TRUME) light |
US20100056924A1 (en) * | 2006-11-20 | 2010-03-04 | Koninklijke Philips Electronics N.V. | Control and display of ultrasonic microbubble cavitation |
JP5081251B2 (ja) * | 2006-12-08 | 2012-11-28 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 結合された治療を計画するシステム、方法、コンピュータ読み取り可能な記録媒体及び使用 |
JP4928989B2 (ja) * | 2007-03-07 | 2012-05-09 | 株式会社東芝 | 超音波診断装置および超音波診断装置の制御プログラム |
US7713209B2 (en) * | 2007-05-21 | 2010-05-11 | Siemens Medical Solutions Usa, Inc. | Targeted contrast agent imaging with medical diagnostic ultrasound |
EP2419020A4 (en) * | 2009-04-17 | 2014-01-15 | Visualsonics Inc | METHOD FOR NONLINEAR IMAGING OF ULTRASONIC CONTRASTING AT HIGH FREQUENCIES |
JP6139186B2 (ja) | 2013-03-11 | 2017-05-31 | 東芝メディカルシステムズ株式会社 | 超音波診断装置、画像処理装置及び画像処理プログラム |
US9743909B1 (en) * | 2013-05-15 | 2017-08-29 | University Of Washington Through Its Center For Commercialization | Imaging bubbles in a medium |
JP6426747B2 (ja) * | 2014-01-23 | 2018-11-21 | スーパー・ソニック・イマジン | 媒体内の正確な位置における物理的特性を決定するための方法、媒体の画像を決定するための方法、及び当該方法を実装する装置 |
EP3097538B1 (en) * | 2014-01-23 | 2018-09-26 | Koninklijke Philips N.V. | Evaluation of carotid plaque using contrast enhanced ultrasonic imaging |
JP6547612B2 (ja) * | 2015-12-14 | 2019-07-24 | コニカミノルタ株式会社 | 画像処理装置、画像処理方法、および、画像処理装置を備える超音波診断装置 |
JP6677042B2 (ja) * | 2016-03-25 | 2020-04-08 | コニカミノルタ株式会社 | 画像処理装置、超音波診断装置及び画像処理プログラム |
WO2018209193A1 (en) * | 2017-05-11 | 2018-11-15 | Verathon Inc. | Probability map-based ultrasound scanning |
US20190154822A1 (en) * | 2017-11-21 | 2019-05-23 | The Charles Stark Draper Laboratory, Inc. | System and method for imaging and localization of contrast-enhanced features in the presence of accumulating contrast agent in a body |
CN111902074A (zh) * | 2018-03-13 | 2020-11-06 | 博信生物科技股份有限公司 | 用于灵敏分子分析的组合物和方法 |
-
2019
- 2019-03-27 JP JP2020509206A patent/JP6963677B2/ja active Active
- 2019-03-27 WO PCT/JP2019/013237 patent/WO2019189386A1/ja active Application Filing
- 2019-03-27 EP EP19774254.7A patent/EP3777699A4/en active Pending
-
2020
- 2020-09-04 US US17/013,051 patent/US11969296B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
EP3777699A1 (en) | 2021-02-17 |
EP3777699A4 (en) | 2021-05-26 |
JPWO2019189386A1 (ja) | 2021-02-18 |
WO2019189386A1 (ja) | 2019-10-03 |
US11969296B2 (en) | 2024-04-30 |
US20200397412A1 (en) | 2020-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3462584B2 (ja) | 超音波診断装置 | |
JP6963677B2 (ja) | 超音波診断装置および超音波診断装置の制御方法 | |
JPWO2006121034A1 (ja) | 超音波診断装置 | |
US11331080B2 (en) | Ultrasound diagnostic apparatus that generates an ultrasound image using a harmonic imaging method and method of controlling ultrasound diagnostic apparatus that generates an ultrasound image using a harmonic imaging | |
JP2018029788A (ja) | 超音波診断装置 | |
US10617395B2 (en) | Ultrasound diagnostic apparatus and doppler waveform image generating method | |
KR102245671B1 (ko) | 음향 방사력-기반 초음파 이미징에서의 적응식 클러터 필터링 | |
JP5148194B2 (ja) | 超音波診断装置 | |
JP6494784B2 (ja) | 超音波診断装置および超音波診断装置の制御方法 | |
JP5869411B2 (ja) | 超音波診断装置および超音波画像生成方法 | |
WO2021020035A1 (ja) | 超音波診断装置および超音波診断装置の制御方法 | |
JP2015128554A (ja) | 超音波診断装置 | |
JP4912982B2 (ja) | 超音波診断装置 | |
US11751853B2 (en) | Ultrasound diagnostic apparatus and ultrasound image generating method | |
TW201232476A (en) | Detection system and signal processing method thereof | |
JPH09299369A (ja) | 超音波連続波ドプラ血流計 | |
JP7095177B2 (ja) | 超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサ | |
JP6434642B2 (ja) | 音速算出システムおよび音速算出方法 | |
JPWO2018181127A1 (ja) | 超音波診断装置および超音波診断装置の制御方法 | |
WO2017122411A1 (ja) | 超音波診断装置および音速定量化方法 | |
JP6536357B2 (ja) | 超音波画像診断装置 | |
JP2011218049A (ja) | 超音波診断装置 | |
JP2011212066A (ja) | 環境音速取得方法および装置 | |
JP6364084B2 (ja) | 音響波診断装置およびその制御方法 | |
JP5299427B2 (ja) | 超音波診断装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200909 |
|
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: 20210928 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20211015 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6963677 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |