JP6604272B2 - 超音波診断装置及び超音波信号処理方法 - Google Patents

超音波診断装置及び超音波信号処理方法 Download PDF

Info

Publication number
JP6604272B2
JP6604272B2 JP2016112802A JP2016112802A JP6604272B2 JP 6604272 B2 JP6604272 B2 JP 6604272B2 JP 2016112802 A JP2016112802 A JP 2016112802A JP 2016112802 A JP2016112802 A JP 2016112802A JP 6604272 B2 JP6604272 B2 JP 6604272B2
Authority
JP
Japan
Prior art keywords
wave
detection
transmission
shear wave
push pulse
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
JP2016112802A
Other languages
English (en)
Other versions
JP2017217147A (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.)
Konica Minolta Inc
Original Assignee
Konica Minolta Inc
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 Konica Minolta Inc filed Critical Konica Minolta Inc
Priority to JP2016112802A priority Critical patent/JP6604272B2/ja
Priority to US15/612,295 priority patent/US11399805B2/en
Publication of JP2017217147A publication Critical patent/JP2017217147A/ja
Application granted granted Critical
Publication of JP6604272B2 publication Critical patent/JP6604272B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/52019—Details of transmitters
    • G01S7/5202—Details of transmitters for pulse systems
    • G01S7/52022—Details of transmitters for pulse systems using a sequence of pulses, at least one pulse manipulating the transmissivity or reflexivity of the medium
    • 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
    • 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/485—Diagnostic techniques involving measuring strain or elastic properties
    • 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
    • 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/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • 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
    • 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/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • 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/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • 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/467—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B8/469—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means for selection of a region of interest

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biophysics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Animal Behavior & Ethology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Acoustics & Sound (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

本発明は、超音波診断装置における画像化処理に関し、特に、せん断波を用いた組織の硬さ測定に関する。
近年、被検体内の組織の硬さを評価する機能を有する超音波診断装置が普及してきている。特に、被検体内の着目領域（ＲＯＩ；Ｒｅｇｉｏｎ ｏｆ Ｉｎｔｅｒｅｓｔ）にせん断波（ｓｈｅａｒ ｗａｖｅ）を生じさせ、着目領域内の組織の変位を時系列に取得することで、せん断波の伝播速度を計測する方法が広く用いられている。せん断波の伝播速度は組織の弾性率に応じて変化するため、組織の絶対的な硬さ（例えば、弾性率）を評価することができるからである。せん断波を生じさせる手法としては、例えば、ＡＲＦＩ（Ａｃｏｕｓｔｉｃ Ｒａｄｉａｔｉｏｎ Ｆｏｒｃｅ Ｉｍｐｕｌｓｅ）と呼ばれる、プッシュパルスを焦点に集中させることで、超音波であるプッシュパルスの音圧により、焦点上にある被検体内の組織に変位を発生させる手法が用いられる。
特開２０１４−１１３３２３号公報
一方で、せん断波の伝播速度を測定するには、着目領域内の組織の変位の時間変化を観測するため、プッシュパルスに続いて観測のための検出波（トラッキングパルス）の送受信を高速度に繰り返す必要がある。そのため、特許文献１に開示されているように、高頻度の超音波の送受信によって超音波探触子が発熱しやすくなり、また、検出波の受信で得た受信信号が多くなるため変位検出のための演算量が増加するという課題がある。
一方で、せん断波の伝播速度の測定精度は検出波の送信頻度、すなわち、フレームレートに連動して高くなる。また、せん断波の観測時間は着目領域の広さに依存する。そのため、単純に検出波の送受信頻度を低下させたり、送受信回数を削減したりすると、せん断波の伝播速度の精度が低下する恐れがある。
本開示は、上記開示を鑑みてなされたものであり、検出波として平面波を用いる場合に、せん断波の伝播速度の精度低下を抑止しながら、プッシュパルスに続く一連の検出波の送受信において、平均送受信頻度を低下させることを目的とする。
本開示の一態様に係る超音波診断装置は、超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して特定部位にある組織を物理的に押圧した後検出波の送受信を被検体内に対して繰り返し行うことで、押圧された特定部位の組織を振動源とするせん断波の伝播速度を検出する超音波診断装置であって、プッシュパルスを送信するプッシュパルス送信部と、前記プッシュパルスに続き平面波を検出波として被検体内の着目領域へ複数回送信し、検出波に対応する被検体からの反射検出波を受信して複数の受信信号を時系列に生成する検出波送受信部と、前記反射検出波の受信時刻のそれぞれにおける、前記プッシュパルスに起因するせん断波による前記着目領域内の組織の変位を前記複数の受信信号からそれぞれ検出する変位検出部と、前記変位検出部が検出した変位に基づいて前記着目領域内におけるせん断波の伝播速度を検出するせん断波解析部とを備え、前記検出波送受信部は、プッシュパルスに続く一連の検出波の送信において、連続する２つの検出波の送信間隔が、第１の送信間隔と、前記第１の送信間隔より長い第２の送信間隔とを少なくとも含むように、検出波を送信することを特徴とする。
本開示によれば、上記構成により、検出波の送信間隔を調整することで、せん断波の伝播速度の精度低下を抑止しながら検出波の送受信頻度を低下させることができる。
実施の形態１に係る超音波診断装置１のブロック図である。 実施の形態１に係る超音波診断装置１の動作を示すフローチャートである。 各実施の形態に係る検出波送受信の動作を示すフローチャートである。 （ａ）は検出波の送信プロファイルの一例を示す模式図であり、（ｂ）および（ｃ）は送信プロファイルとせん断波の波面との関係を示す模式図である。 せん断波の発生と伝播の様子を示す模式図である。 実施の形態１に係るせん断波の伝播解析の動作を示す模式図である。 実施の形態１に係る弾性率測定の処理を示す模式図である。 （ａ）〜（ｃ）は、送信プロファイルとせん断波の波面との関係を示す模式図である。 実施の形態２に係る超音波診断装置の動作を示すフローチャートである。 実施の形態２に係る検出波の送信プロファイル決定方法を示す模式図である。 実施の形態３に係る超音波診断装置の動作を示すフローチャートである。 実施の形態３に係るせん断波の伝播解析結果の合成動作を示す模式図である。 実施の形態３の変形例に係る超音波診断装置の動作を示すフローチャートである。 実施の形態３の変形例に係る検出波の送信プロファイル決定方法を示す模式図である。
≪発明を実施するための形態に至った経緯≫
発明者は、せん断波を用いて組織の硬さの測定を行う超音波診断装置において、測定精度の低下を抑止しながら平面波たる検出波の送受信頻度を低下させるために各種の検討を行った。なお、本開示では、平面波を検出波として用いる場合について検討を行っている。これは、平面波を検出波として用いる場合、１送信で着目領域全域に対して超音波の送受信が行えるため、検出波の送信頻度を変位検出のためのＢモード画像のフレームレートより高くする必要がないためである。
上述したように、検出波の送受信頻度が高くなると、超音波探触子が発熱する。特に、検出波の送受信頻度が高い状態を長時間維持すると、超音波探触子が発熱しやすい。また、検出波として平面波を用いる場合には、焦点波を検出波として用いる場合と比べ、超音波送信に用いる振動子の数が多く、より発熱量が大きい。その結果、超音波探触子の発熱によって短時間に超音波探触子の温度上限に達する可能性があり、超音波診断装置の連続動作の妨げになっている。この問題を解決するためには超音波探触子が発熱量を減らす、および／または発熱時間を短縮・低頻度化する必要がある。したがって、検出波の送受信頻度を削減し、および／または、検出波の送受信頻度が高い状態が維持される時間を短縮することが好ましい。また、送受信回数が増加すると、取得される受信信号が比例して増加するため、演算量が膨大となる。演算量が増加すると超音波診断装置がユーザ（検査者）に結果を提示するための所要時間が長くなり、結果のリアルタイム性が低下し、使い勝手の低下に結びつくこととなる。
一方で、検出波の送受信頻度を単純に低下させると、せん断波の伝播解析において時間解像度が低下する。結果として、せん断波の速度が時間方向および空間方向に平均化され、小さな硬い異物を発見しづらくなる、硬い組織と周辺組織との界面が不明瞭になる、など、測定精度の低下に結びつくこととなる。また、検出波の送受信時間を単純に短縮すると、着目領域の範囲が狭小化することとなる。これは、せん断波の速度が組織の硬さに依存するため、観測時間を短縮すれば、当然に観測時間中にせん断波が移動する範囲が狭くなるためである。
そこで、発明者は、上記課題に鑑み、せん断波の速度検出精度と着目領域の面積とに影響を与えないように検出波の送受信頻度を低下させる技術について実施の形態に係る超音波診断装置に相当するに至ったものである。
以下、実施の形態に係る超音波診断装置について図面を用いて詳細に説明する。
≪実施の形態≫
実施の形態に係る超音波診断装置１のブロック図を図１に示す。超音波診断装置１は、制御部１１、せん断波励起部１２、超音波信号取得部１３、変位検出部１４、伝播解析部１５、断層画像記憶部１６、変位量記憶部１７、波面画像記憶部１８、弾性画像記憶部１９を備える。また、制御部１１には、超音波探触子２と表示部３とがそれぞれ接続可能に構成されている。図１は超音波診断装置１に超音波探触子２、表示部３が接続された状態を示している。
超音波探触子２は、例えば、一次元方向に配列された複数の振動子（不図示）を有する。各振動子は、例えば、ＰＺＴ（チタン酸ジルコン酸鉛）により構成される。超音波探触子２は、せん断波励起部１２で生成された電気信号（以下、「ＡＲＦＩ駆動信号」と呼ぶ）、または、超音波信号取得部１３で生成された電気信号（以下、「検出駆動信号」と呼ぶ）を制御部１１から受け取り、超音波に変換する。超音波探触子２は、超音波探触子２の振動子側外表面を被検体の皮膚表面など表面に接触させた状態で、ＡＲＦＩ駆動信号又は検出駆動信号から変換され、複数の振動子から発せられる複数の超音波からなる超音波ビームを被検体内の測定対象に向けて送信する。そして、超音波探触子２は、検出駆動信号に基づく送信検出波に対する、測定対象からの複数の反射検出波を受信し、複数の振動子によりこれら反射超音波をそれぞれ電気信号（以下、「素子受信信号」と呼ぶ）に変換し、制御部１１を介して素子受信信号を超音波信号取得部１３に供給する。なお、ここではせん断波励起部１２及び超音波信号取得部１３を別構成として説明しているが、超音波信号取得部１３の検出駆動信号を生成する構成と同一の構成により、ＡＲＦＩ駆動信号を生成してもよい。
せん断波励起部１２は、プッシュパルスを超音波探触子２に送出させるための電気信号であるＡＲＦＩ駆動信号を生成する。プッシュパルスとは、被検体内にせん断波を発生させるため、被検体内の組織に変位をもたらすためのパルス状の超音波である。具体的には、被検体内の着目領域内のある１点を焦点とする、波数が後述する送信超音波より多い超音波である。したがって、ＡＲＦＩ駆動信号は、超音波探触子２を構成する各振動素子から送出される超音波が焦点に届くように生成される、振動素子ごとにタイミングの異なるパルス状の電気信号である。
超音波信号取得部１３は、送信検出波を超音波探触子２に送出させるための電気信号である検出駆動信号を生成する。検出駆動信号は、超音波探触子２を構成する各振動素子から送出される検出波の波面が、進行方向に対して直交する平面となるように生成される電気信号である。具体的には、すべての振動子の駆動タイミングが同一となるように、または、振動子列の一方の端から他方の端まで固定ピッチで段階的に動作タイミングをずらした、検出駆動信号を生成する。これにより、送信検出波は着目領域全体にいきわたるように送信される。また、超音波信号取得部１３は、反射検出波に基づく素子受信信号に整相加算を行って、着目領域全体の音響線信号を生成する。超音波信号取得部１３は、生成した音響線信号を、制御部１１を介して断層画像記憶部１６に出力する。超音波信号取得部１３は、所定のプロファイルに従って検出駆動信号を制御部１１に出力することで、所定のタイミングで検出波を繰り返し送信させる。詳細は後述する。
変位検出部１４は、変位検出の対象となる１つの断層画像に係る複数の音響線信号（以下、「断層画像信号」と呼ぶ）と、基準となる１つの断層画像に係る複数の音響線信号（以下、「基準断層画像信号」と呼ぶ）とを、制御部１１を介して断層画像記憶部１６から取得する。基準断層画像信号とは、断層画像信号からせん断波による変位を抽出するために用いるものであり、具体的には、プッシュパルス送出前に着目領域を撮像した断層画像信号である。そして、変位検出部１４は、断層画像信号と基準断層画像信号との差分から、断層画像信号の各画素の変位を検出し、変位を各画素の座標と関連付けて変位画像を生成する。変位検出部１４は、生成した変位画像を、制御部１１を介して変位量記憶部１７に出力する。
伝播解析部１５は、変位画像を、制御部１１を介して変位量記憶部１７から取得する。伝播解析部１５は、変位画像から、変位画像を取得した各時刻における、せん断波の波面の位置、進行方向および速度を検出し、波面画像を生成する。伝播解析部１５は、せん断波の波面の位置、進行方向および速度から、変位画像の各画素に対応する被検体組織の弾性率を算出し、弾性画像を生成する。伝播解析部１５は、生成した波面画像を波面画像記憶部１８に、弾性画像を弾性画像記憶部１９に、制御部１１を介してそれぞれ出力する。
制御部１１は、上述したような各構成要素の制御に加え、伝播解析部１５が生成した弾性画像を表示部３に出力する。
断層画像記憶部１６、変位量記憶部１７、波面画像記憶部１８、弾性画像記憶部１９は、それぞれ、断層画像、変位画像、波面画像、弾性画像データを記憶する。断層画像記憶部１６、変位量記憶部１７、波面画像記憶部１８、弾性画像記憶部１９のそれぞれは、例えば、ＲＡＭ、フラッシュメモリ、ハードディスク、光学ディスクなどの記憶媒体で実現される。なお、断層画像記憶部１６、変位量記憶部１７、波面画像記憶部１８、弾性画像記憶部１９のうち２以上を単一の記憶媒体で実現してもよい。また、断層画像記憶部１６、変位量記憶部１７、波面画像記憶部１８、弾性画像記憶部１９のうち１以上は、超音波診断装置１の外部に構成され、ＵＳＢやｅＳＡＴＡなどのインターフェースを介して超音波診断装置１と接続されてもよいし、超音波診断装置１からネットワークを介してアクセス可能に構成された資源、例えば、ファイルサーバやＮＡＳ（Ｎｅｔｗｏｒｋ Ａｔｔａｃｈｅｄ Ｓｔｏｒａｇｅ）であってもよい。
制御部１１、せん断波励起部１２、超音波信号取得部１３、変位検出部１４、伝播解析部１５のそれぞれは、例えば、ＦＰＧＡ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｇｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）などのハードウェアにより実現される。なお、これらの一部または全部は、単一のＦＰＧＡ、または、ＡＳＩＣで実現されてもよい。また、これらは、それぞれ個別に、または、２以上を１まとめとして、メモリと、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＧＰＵ（Ｇｒａｐｈｉｃ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）などのプログラマブルデバイスとソフトウェアで実現されてもよい。
＜動作＞
実施の形態１に係る超音波診断装置１の動作について説明する。図２は、超音波診断装置１全体の動作を示すフローチャートである。
まず、制御部１１が、着目領域を設定する（ステップＳ１０）。着目領域を設定する方法は、例えば、表示部３に断層画像記憶部１６に記録されている最新の断層画像を表示し、タッチパネル、マウス、トラックボールなどの入力部（図示しない）を通して検査者に着目領域を指定させる。なお、着目領域の設定方法はこの場合に限られず、例えば、断層画像の全域を着目領域としてもよいし、あるいは、断層画像の中央部分を含む一定範囲を着目領域としてもよい。また、着目領域を設定する際に、改めて断層画像を取得してもよい。
次に、被検体に超音波を送受信し、取得した受信信号を記憶する（ステップＳ２０）。具体的には、次のような動作となる。まず、以下のように送信イベントを行う。最初に、超音波信号取得部１３が、パルス状の送信信号を生成する。次に、超音波信号取得部１３は、送信信号に対して、超音波探触子２の素子毎の遅延時間を設定する送信ビームフォーミングを行い、超音波探触子２の各素子に対応する複数の検出駆動信号を生成する。超音波探触子２の各振動子は対応する検出駆動信号を超音波に変換することで、超音波ビームが被検体内に送出される。次に、超音波探触子２の各振動子は、被検体内から反射された反射超音波を取得して素子受信信号に変換する。超音波信号取得部１３は、素子受信信号に整相加算を行い、音響線信号を生成する。制御部１１は、１回の送信イベントで生成した１つの断層画像を構成する複数の音響線信号を超音波信号取得部１３から取得し、基準断層画像信号として、断層画像記憶部１６に記憶する。
次に、プッシュパルスの送信終了からの経過時間が長くなると検出波の送信間隔が長くなるように、検出波の送信プロファイルを決定する（ステップＳ３０）。本実施の形態では、図４（ａ）に示すような所定の送信プロファイルを用いる。図４（ａ）は、本実施の形態に係る検出波の送信プロファイルを示すタイムチャートである。ここでは、プッシュパルス４００の送信が完了した後、第１期間Ｐ１においては、検出波の送信間隔をＩ１とする。すなわち、第１の検出波４１１−１の送信後、Ｉ１の間隔を置いて第２の検出波４１１−２が送信される。最後に、ｎ番目の検出波４１１−ｎが送信される。次に、第１期間Ｐ１に続く第２期間Ｐ２においては、検出波の送信間隔を、Ｉ１より長いＩ２とする。すなわち、検出波４１１−ｎの送信後、Ｉ２の間隔を置いて検出波４１２−１が送信され、さらにＩ２の間隔を置いて検出波４１２−２が送信される。最後に、検出波４１２−ｍが送信される。次に、第２期間Ｐ２に続く第３期間Ｐ３においては、検出波の送信間隔を、Ｉ２より長いＩ３とする。すなわち、検出波４１２−ｍの送信後、Ｉ３の間隔を置いて検出波４１３−１が送信される。つまり、プッシュパルス４００の送信時刻と検出波の送信時刻との間の時間が長くなるに従い、段階的に送信間隔が大きくなる送信プロファイルである。より具体的には、例えば、第１期間Ｐ１は３．８ｍｓであり、送信間隔Ｉ１は１００μｓである（つまり、送信回数は３８回である）。また、第２期間Ｐ２は３．８ｍｓであり、送信間隔Ｉ２は２００μｓである（送信回数は１９回である）。また、第３期間Ｐ３は１１．２ｍｓであり、送信間隔Ｉ３は４００μｓである（送信回数は２８回である）。また、第４期間Ｐ４は１２ｍｓであり、送信間隔Ｉ４は８００μｓである（送信回数は１５回である）。つまり、１００回の送信を３０．８ｍｓかけて行う。
続いて、プッシュパルスを送信する（ステップＳ４０）。具体的には、せん断波励起部１２が、パルス状のＡＲＦＩ信号を生成する。次に、せん断波励起部１２は、ＡＲＦＩ信号に対して、超音波探触子２の素子毎の遅延時間を設定する送信ビームフォーミングを行い、超音波探触子２の各素子に対応する複数のＡＲＦＩ駆動信号を生成する。超音波探触子２の各振動子は対応する送信駆動信号を超音波に変換することで、プッシュパルスが被検体内に送出される。
ここで、プッシュパルスによるせん断波の生成について、図５（ａ）〜（ｅ）の模式図を用いて説明する。図５（ａ）は、着目領域に対応した被検体内領域の、プッシュパルス印加前における組織を示した模式図である。図５（ａ）〜（ｅ）において、個々の“○”は、着目領域における被検体内の組織の一部を、破線の交点は、負荷がない場合の組織”○“の中心位置を、それぞれ示している。ここで、超音波探触子２を皮膚表面１００に密接させた状態で焦点１０１に対してプッシュパルスを印加すると、図５（ｂ）の模式図に示すように、焦点１０１に位置していた組織１３２が、プッシュパルスの進行方向に押されて移動する。また、組織１３２からプッシュパルスの進行方向側にある組織１３３は、組織１３２に押されてプッシュパルスの進行方向に移動する。次に、プッシュパルスの送信が終了すると、組織１３２、１３３が元の位置に復元しようとするので、図５（ｃ）の模式図に示すように、組織１３１〜１３３がプッシュパルスの進行方向に沿った振動を開始する。すると、図５（ｄ）の模式図に示すように、振動が組織１３１〜１３３に隣接する、組織１２１〜１２３および組織１４１〜１４３に伝播する。さらに、図５（ｅ）の模式図に示すように、振動がさらに組織１１１〜１１３および組織１５１〜１５３に伝播する。したがって、被検体内において、振動が振動の方向と直交する向きに伝播する。すなわち、せん断波がプッシュパルスの印加場所に発生し、被検体内を伝播する。
図２に戻って説明を続ける。次に、着目領域に検出波を複数回送受信し、取得した複数の超音波信号を記憶する（ステップＳ５０）。以下、図３を用いて詳細に説明する。図３は、ステップＳ５０の詳細を示すフローチャートである。まず、検出波送受信部１３は、第１タイマーを起動した後（ステップＳ５１）、第１タイマーの値が、１回目の検出波の送信タイミングになるまで待機する（ステップＳ５２）。次に、期間カウンタｊを１に初期化した後（ステップＳ５３）、第２タイマーを起動して（ステップＳ５４）、検出波を着目領域に送信する（ステップＳ５５）。ステップＳ５５においては、ステップＳ２０と同様、超音波信号取得部１３は、パルス状の送信信号を生成し、送信信号に対して超音波探触子２の素子毎の遅延時間を設定する送信ビームフォーミングを行い、超音波探触子２の各素子に対応する複数の検出駆動信号を生成する。超音波探触子２の各振動子は対応する検出駆動信号を超音波に変換することで、検出波が被検体内に送出される。次に、超音波探触子２の各振動子は、被検体内から反射された反射超音波を取得して素子受信信号に変換する。超音波信号取得部１３は、素子受信信号を取得して記憶する（ステップＳ５６）。次に、超音波信号取得部１３は、第１タイマーの値が、第１期間〜第ｊ期間の長さ（ここでは、第１期間Ｐ１の長さ）を超えているかどうかを判定する（ステップＳ５７）。１回目の検出波の送信後最初のステップＳ５７ではＮｏとなるので、超音波信号取得部１３は、第２タイマーの値が送信間隔Ｉｊ（ここでは、送信間隔Ｉ１）となるまで待機した後（ステップＳ５８）、第２タイマーをゼロにリセットして（ステップＳ５４）、検出波を着目領域に送信する（ステップＳ５５）。つまり、ステップＳ５４からＳ５８のループによって検出波が繰り返し送信間隔Ｉ１で送信される。ここで、第１期間Ｐ１が完了すると、最後の検出波の送信後、ステップＳ５７においてＹｅｓとなるので、第２期間Ｐ２が存在すれば（ステップＳ５９でＹｅｓ）、ｊを２にインクリメントし（ステップＳ６１）、第２タイマーの値が送信間隔Ｉｊ（ここでは、送信間隔Ｉ２）となるまで待機した後（ステップＳ５８）、第２タイマーをゼロにリセットして（ステップＳ５４）、検出波を着目領域に送信する（ステップＳ５５）。これにより、第２期間Ｐ２の最初の検出波は第１期間Ｐ１の最後の検出波より送信間隔Ｉ２だけ遅れて送信される。以下、第１期間Ｐ１と同様の処理により、第２期間Ｐ２の間、検出波が繰り返し送信間隔Ｉ２で送信される。ここで、第２期間Ｐ２が完了すると、同様に、最後の検出波の送信後、ステップＳ５７においてＹｅｓとなるので、第３期間Ｐ３が存在すれば（ステップＳ５９でＹｅｓ）、ｊを３にインクリメントし（ステップＳ６１）、第２タイマーの値が送信間隔Ｉｊ（ここでは、送信間隔Ｉ３）となるまで待機した後（ステップＳ５８）、第２タイマーをゼロにリセットして（ステップＳ５４）、検出波を着目領域に送信する（ステップＳ５５）。この繰り返しにより、最後の第ｊ期間まで、検出波が送信プロファイル通りに送信される。最後に、超音波信号取得部１３は、検出波ごとに、素子受信信号に整相加算を行い、音響線信号を生成する（ステップＳ６２）。制御部１１は、検出波ごとに生成した複数の音響線信号を超音波信号取得部１３から取得し、断層画像信号として、断層画像記憶部１６に記憶する。
次に、各画素の変位を検出する（ステップＳ６０）。具体的には、まず、変位検出部１４が、ステップＳ２０で断層画像記憶部１６に記憶された基準断層画像信号を取得する。次に、変位検出部１４は、ステップＳ５０で断層画像記憶部１８に記憶された各断層画像信号に対し、基準断層画像信号との差分から、当該断層画像信号に係る反射検出波を受信した時刻における、各画素の変位を検出する。より具体的には、例えば、断層画像信号と基準断層画像信号との相関処理を行うことで、断層画像信号上の画素が基準断層画像信号上のいずれの画素に対応するのかを検索し、その座標差を、断層画像信号上の画素に対応する変位として特定する。なお、変位の検出方法は相関処理に限られず、例えば、パターンマッチングなど、２つの断層画像信号間の動き量を検出する任意の技術を用いてよい。パターンマッチングの例としては、例えば、断層画像信号を８ピクセル×８ピクセルなどの所定の大きさの領域に分割し、各領域と基準断層画像信号とをパターンマッチングすることで、断層画像信号の各画素の変位を検出することができる。パターンマッチングの方法としては、例えば、各領域と基準断層画像信号内の同サイズの基準領域との間で、対応する画素毎に輝度値の差分を算出してその絶対値の合計値を算出し、その合計値が最も小さくなる領域と基準領域との組み合わせについて、領域と基準領域とが同一の領域であるものとし、領域の基準点（例えば、左上の角）と基準領域の基準点との距離を変位として検出する。なお、所定の大きさの領域は他のサイズであってもよいし、輝度値の差分の絶対値の合計値に替えて、例えば、輝度値の差分の２乗の合計値を用いてもよい。また、相関処理やパターンマッチングにより変位を検出する際、対応する画素の座標差ではなく、対応する画素のｙ座標の差（深さの差）を変位の大きさとして用いてもよい。これは、せん断波の伝播方向が原則として素子列方向（ｘ方向）であるため、せん断波による変位は伝播方向とは直交する方向であり、原則として深さ方向（ｙ方向）であるからである。以上の処理により、各断層画像信号の各画素に対応する被検体の組織が、プッシュパルスまたはせん断波によってどれだけ動いたかが変位として算出される。変位検出部１４は、１の断層画像に係る各画素の変位を当該画素の座標と対応付けることで変位画像を生成し、生成した変位画像を変位量記憶部１７に出力する。
次に、せん断波の伝播解析を行う（ステップＳ７０）。具体的には、各変位画像からせん断波の波面を抽出して波面画像を生成する。この波面画像より、波面の位置、振幅、進行方向および速度を容易に検出することができる。波面画像の生成は、例えば、変位領域の抽出、細線化処理、空間フィルタリング、時間フィルタリングの手順によって行う。
図６を用いて具体的な処理を説明する。図６（ａ）は、変位画像の一例を示している。図５と同じく、図中の“○”は着目領域における被検体内の組織の一部を示しており、プッシュパルスを印加する前の位置は破線の交点である。伝播解析部１５は、ｙ座標ごとに変位量δを座標ｘの関数として、動的閾値を用いることで変位量δが大きい領域を抽出する。また、ｘ座標ごとに変位量δを座標ｙの関数として、動的閾値を用いて、ある閾値を超える領域を変位量δが大きい領域として抽出する。動的閾値とは、対象領域内について信号解析又は画像解析を行って閾値を決定することである。閾値は一定値ではなく、対象領域の信号の幅や最大値などによって異なる値となる。図６（ａ）に、ｙ＝ｙ1の直線２１０上における変位量をプロットしたグラフ２１１と、ｘ＝ｘ1の直線２２０上における変位量をプロットしたグラフ２２１とを示す。これにより、例えば、変位量δが閾値より大きな変位領域２３０が抽出できる。
次に、伝播解析部１５は、変位領域に細線化処理をおこなって波面を抽出する。図６（ｂ）の模式図に示している変位領域２４０、２５０は、それぞれ、ステップＳ５２において変位領域として抽出された領域である。伝播解析部１５は、例えば、Ｈｉｌｄｉｔｃｈの細線化アルゴリズムを用いて、波面を抽出する。例えば、図６（ｂ）の模式図において、変位領域２４０から波面２４１が、変位領域２５０から波面２５１が、それぞれ抽出される。なお、細線化のアルゴリズムはＨｉｌｄｉｔｃｈに限らず、任意の細線化アルゴリズムを用いてよい。また、各変位領域に対して、変位量δが閾値以下の座標を変位領域から取り除く処理を、変位領域が幅１ピクセルの線になるまで、閾値を大きくしながら繰り返し行ってもよい。
次に、伝播解析部１５は、細線化処理後の波面画像データに対して空間フィルタリングを行い、長さが短い波面を除去する。例えば、ステップＳ５３で抽出した各波面の長さを検出し、全ての波面の長さの平均値の１／２よりも長さが短い波面を、ノイズとして削除する。具体的には、図６（ｃ）の波面画像に示すように、波面２６１〜２６４の長さの平均値を算出し、それよりも短い波面２６３、２６４を、ノイズとして消去する。これにより、誤検出された波面を消去できる。
伝播解析部１５は、変位領域の抽出、細線化処理、空間フィルタリングの動作を、全ての変位画像に対して行う。これにより、変位画像に対して１対１で波面画像データが生成される。
最後に、伝播解析部１５は、複数の波面画像データに対して時間フィルタリングを行い、伝播していない波面を除去する。具体的には、時間的に連続する２以上の波面画像において、波面位置の時間変化を検出し、速度が異常である波面をノイズとして除去する。伝播解析部１５は、例えば、時刻ｔ＝ｔ1の波面画像２７０、時刻ｔ＝ｔ1＋Δｔの波面画像２８０、時刻ｔ＝ｔ1＋２Δｔの波面画像２９０との間で、波面位置の時間変化を検出する。例えば、波面２７１に対して、波面画像２８０のうち、波面２７１と同じ位置を中心に、波面と垂直な向き（図７においてはｘ軸方向）にΔｔの間にせん断波が移動しうる領域２７６で、波面２７１との相関処理を行う。このとき、波面２７１のｘ軸の正方向（図の右側）と負方向（図の左側）の双方を含む範囲内で相関処理を行う。これは、透過波と反射波の両方を検出するためである。これにより、波面２７１の移動先が波面画像２８０内の波面２８１であると検出し、時間Δｔにおける波面２７１の移動距離を算出する。同様に、波面２７２、２７３のそれぞれについて、波面画像２８０において当該波面と同じ位置を中心に、波面と垂直な向きにΔｔの間にせん断波が移動しうる領域で相関処理を行う。これにより、波面２７２が波面２８３の位置に、波面２７３が波面２８２の位置に、それぞれ移動したことを検出する。波面画像２８０と波面画像２９０との間でも同様の処理を行い、波面２８１が波面２９１の位置に、波面２８２が波面２９２の位置に、波面２８３が波面２９３の位置に、それぞれ移動したことを検出する。ここで、波面２７３、波面２８２、波面２９２で示される１の波面については、他の波面と比べて移動距離が著しく小さい（伝播速度が著しく遅い）。このような波面は誤検知である可能性が高いので、ノイズとして消去する。これにより、図６（ｅ）の波面画像３００に示すように、波面３０１、３０２が検出できる。
伝播解析部１５は、生成した時刻ごとの波面画像データと、波面の対応情報とを用いて、波面の位置と速度を算出する。ここで、波面の対応情報とは、同一の波面が各波面画像のどの波面に対応するかを示した情報であり、例えば、図６（ｄ）において、波面２７２が波面２８３の位置に移動したことが検出された場合、波面２８３と波面２７２とが同一の波面であるという情報である。以下、図７を用いて波面の速度算出について説明する。図７（ａ）は、ある時刻ｔ1における波面画像と、時刻ｔ2（ｔ1＜ｔ2）における波面画像を１つの波面画像３１０として合成したものである。ここで、時刻ｔ1における波面３１１と、時刻ｔ2における波面３１２とが同一の波面であるとする対応情報が存在するものとする。伝播解析部１５は、対応情報から、波面３１１上の座標（ｘt1、ｙt1）に対応する波面３１２上の座標（ｘt2、ｙt2）を検出する。これにより、時刻ｔ1に座標（ｘt1、ｙt1）を通過したせん断波が、時刻ｔ2に座標（ｘt2、ｙt2）に到達していると推定できる。したがって、座標（ｘt1、ｙt1）を通過したせん断波の速度ｖ（ｘt1、ｙt1）は、座標（ｘt1、ｙt1）と座標（ｘt2、ｙt2）との間の距離ｄを所要時間Δｔ＝ｔ2−ｔ1で割った値と推定できる。すなわち、ｖ（ｘt、ｙt）＝ｄ／Δｔ＝√｛（ｘt2−ｘt1）2＋（ｙt2−ｙt1）2｝／Δｔとなる。伝播解析部１５は、全ての波面に対して上述の処理を行い、波面が通過した全座標についてせん断波の速度を取得し、各画素とせん断波の速度とを対応付けて速度分布図を作成してこれを保持する。
図２に戻って説明を続ける。最後に、弾性画像を生成して表示する（ステップＳ８０）。具体的には、まず、伝播解析部１５が、速度分布図の各画素についてせん断波の速度から弾性率を算出し、各画素と弾性率とを対応付けて弾性画像を生成する。座標（ｘt、ｙt）における弾性率Ｅ（ｘt、ｙt）は、当該座標におけるせん断波の速度ｖ（ｘt、ｙt）を用いて、以下のように算出できる。
Ｅ（ｘt、ｙt）＝２（１＋γ）ρ・ｖ（ｘt、ｙt）2
ここで、γは座標（ｘt、ｙt）における組織のポアソン比、ρは密度である。簡易的には、例えば、γ＝０．５、ρ＝１ｇ／ｃｍ3として、以下のように算出してよい。
Ｅ（ｘt、ｙt）≒３・ｖ（ｘt、ｙt）2
各画素と弾性率との対応付けは、例えば、色情報をマッピングすることで行う。これにより、例えば、図７（ｂ）に示すように、弾性率が一定値以上の座標は赤、弾性率が一定値未満の座標は緑、弾性率が取得できなかった座標は黒、というように色分けした弾性画像３２０を生成する。分類は二値化に限られず、所定の段階で分類及び色分けを行ってよい。図７（ｂ）において、領域３２２は弾性率が一定値以上の領域であり、内包物３２１に対応する。なお、図７（ｂ）では説明のために内包物３２１を明示しているが、内包物３２１は実際の弾性画像上には直接現れない。伝播解析部１５は、生成した弾性画像を制御部１１に出力し、制御部１１は弾性画像を弾性画像記憶部１９に出力する。制御部１１は、弾性画像と超音波画像とを表示する。具体的には、制御部１１は、ステップＳ８０で生成された弾性画像と、ステップＳ２０で取得した基準断層画像信号とに対して画面表示用の画像データとなるよう幾何変換を行い、幾何変換後の弾性画像と超音波画像とを表示部３に出力する。
＜検出波の送信プロファイルがせん断波の伝播解析にもたらす効果＞
以下、検出波の送信プロファイルがせん断波の伝播解析にもたらす効果について説明する。
まず、せん断波の伝播速度が速い場合について、図４（ｂ）の模式図を用いて説明する。図４（ｂ）は、第１期間Ｐ１の間にせん断波が着目領域外に伝播した場合について説明する模式図である。
せん断波の伝播速度が速い場合、プッシュパルス焦点４５０で発生したせん断波は、第１期間Ｐ１の間にせん断波が着目領域の外周に達し、外に伝播する。この場合、第２期間Ｐ２以降には波面が検出されないため、着目領域内には、第１期間Ｐ１の波面群４６０、すなわち波面４６１〜４６８のみが検出されることとなる。したがって、連続する２つの検出波の送信間隔は、常に最も短い送信間隔であるＩ１である。そのため、せん断波の伝播解析において、時間解像度が低下する場所がなく、せん断波の速度の精度低下を抑止することができる。
一方、せん断波の伝播速度が遅い場合について、図４（ｃ）の模式図を用いて説明する。図４（ｃ）は、第３期間Ｐ３の間にせん断波が着目領域外に伝播した場合について説明する模式図である。
せん断波の伝播速度が遅い場合、プッシュパルス焦点４７０で発生したせん断波は、第２期間Ｐ２が終了しても着目領域内を伝播している。この場合、第１期間Ｐ１の波面群４８０、すなわち波面４８１〜４８８と、第２期間Ｐ２の波面群４９０、すなわち波面４９１〜４９４と、第３期間Ｐ３の波面群５００、すなわち波面５０１が検出されることとなる。ここで、せん断波の進行速度が遅いことから、第２期間Ｐ２に係る送信間隔Ｉ２におけるせん断波の波面の移動距離が過大となることによる空間解像度および時間解像度の低下が発生せず、せん断波の伝播解析の精度は低下しない。すなわち、第１期間Ｐ１と第２期間Ｐ２の間においては、せん断波の伝播解析には十分に短い間隔で検出波が送信されており、せん断波の速度の精度低下は発生しない。一方、第３期間Ｐ３に取得された波面５０１では、送信間隔Ｉ３におけるせん断波の波面の移動距離が大きくなり、空間解像度および時間解像度が低下する。しかしながら、第３期間Ｐ３ではプッシュパルスの送信時刻より時間が経過しているため、せん断波のエネルギーが減衰しており、せん断波による変位量の絶対値が小さい。そのため、第３期間Ｐ３では、仮に送信間隔Ｉ１で検出波を送信しても、せん断波の伝播解析結果を精度良く得ることができない。したがって、第３期間Ｐ３では、送信間隔Ｉ３がせん断波の伝播解析のために必要な間隔より長いとしても、精度の高くない領域の精度がさらに低下するに過ぎず、その影響は小さい。
また、着目領域内の一部においてせん断波の伝播速度が異なる場合について、図８を用いて説明する。図８（ａ）は、深さ方向（ｙ方向）に一定の帯状の領域８１２におけるせん断波の伝播速度が、周囲の領域８１１、８１３より速い場合について示した模式図である。この場合、せん断波の伝播速度は、領域８１２内では速く、領域８１１、８１３では遅い。そのため、領域８１２の全域に対して、第１期間Ｐ１に取得された波面８０１〜８０８により、伝播解析を行うことができる。したがって、せん断波の伝播速度の速い領域８１２については、せん断波の伝播解析において、時間解像度が低下せず、せん断波の速度の精度低下を抑止することができる。一方、領域８１１、８１３では、第１期間Ｐ１に取得された波面８０１〜８０８に加え、第２期間Ｐ２に取得された波面８２１を用いて伝播解析を行うこととなる。しかしながら、領域８１１、８１３では、せん断波の伝播速度が速くないため、せん断波の伝播解析を精度良く行うには十分に短い間隔で検出波が送信されており、せん断波の速度の精度低下は発生しない。
一方、図８（ｂ）は、硬い組織８４０（せん断波の伝播速度の速い領域）がプッシュパルスの焦点８３０に近い位置に存在する場合を示した模式図である。この場合、硬い組織８４０内をせん断波が第１期間Ｐ１の間に透過する。したがって、硬い組織８４０の内部について、第１期間Ｐ１の波面群８５０、すなわち波面８５１〜８５８を用いて伝播解析を行うことができる。そのため、硬い組織８４０の内部を最も短い送信間隔であるＩ１を用いて取得した波面を用いて解析でき、せん断波の速度の精度低下を抑止することができる。また、硬い組織８４０の外部では、せん断波の伝播速度が速くないため、せん断波の伝播解析を精度良く行うためには送信間隔Ｉ２でも十分短い。そのため、第１期間Ｐ１の波面群８５０、すなわち波面８５１〜８５８と、第２期間Ｐ２の波面群８６０、すなわち波面８６１〜８６３との両方を用いて精度よく伝播解析を行うことができる。
一方、図８（ｃ）は、硬い組織８８０がプッシュパルスの焦点８７０より遠い位置に存在する場合を示した模式図である。この場合、せん断波の伝播速度は、硬い組織８８０の外部では遅いため、せん断波の波面が第１期間Ｐ１のうちに硬い組織８８０に到達せず、第２期間Ｐ２にはいってから、せん断波の波面が硬い組織８８０に到達している。ここで、硬い組織８８０の外部ではせん断波の伝播速度が速くないため、せん断波の伝播解析を精度良く行うためには送信間隔Ｉ２でも十分短い。そのため、第１期間Ｐ１の波面群８９０、すなわち波面８９１〜８９８と、第２期間Ｐ２の波面群９００、すなわち波面９０１〜９０３との両方を用いて伝播解析を行うことができ、せん断波の速度の精度低下は発生しない。一方で、硬い組織８８０については、せん断波の伝播速度が速いのに対し、送信間隔Ｉ２で取得された波面群９００を用いて伝播解析を行うこととなるため、伝播速度の解析精度が低下する。しかしながら、プッシュパルスの焦点と硬い領域との間の距離が大きい場合、せん断波のエネルギーは硬い領域に到達する前に大きく減衰している。さらに、硬い領域の表面ではせん断波の反射が発生する。したがって、波面９０３は反射波波面９０３−２と透過波波面９０３−１に分離しており、さらに、せん断波のエネルギーが小さいため、透過波波面９０３−２における変位量は非常に小さい。したがって、硬い組織８８０のせん断波の伝播解析は、検出波の送信間隔とは無関係に困難であり、検出波を送信間隔Ｉ２で送信した場合の精度低下の影響は無視することができる。
実施の形態に係る検出波の送信タイミングによれば、速度の速いせん断波を観測すべき、プッシュパルスの送信時刻と近い期間においては検出波の送信間隔を短くし、速度の遅いせん断波を観測すべき、プッシュパルスの送信時刻から遠い期間においては検出波の送信間隔を長くする。したがって、せん断波の伝播速度にかかわらず、本実施の形態に係る検出波の送信間隔によれば、せん断波の速度の精度低下を抑止することができる。また、検出波の送信間隔は第１期間Ｐ１を除いて間隔Ｉ１より長くなるため、超音波探触子の発熱を抑止することができる。また、時間当たりの検出波の送信回数は常に間隔Ｉ１で送信する場合より少なくなるため、受信信号の量を削減して演算量を削減することができる。さらに、演算量の削減により、プッシュパルスを送信してから弾性画像を生成するまでの時間を短縮することが可能となり、超音波診断装置のリアルタイム性を向上しユーザビリティの向上に奏功する。
＜小括＞
実施の形態に係る超音波診断装置の構成により、せん断波速度の精度低下を抑止しながら、検出波の送受信頻度を低下させることができる。
≪実施の形態２≫
実施の形態１では、予め用意した検出波の送信プロファイルを用いることで、検出波の送受信頻度を低下させる場合について説明した。
これに対し、実施の形態２では、事前取得した組織の硬さ情報を基に、検出波の送受信プロファイルを最適化することに特徴がある。
＜動作＞
実施の形態２に係る超音波診断装置の動作について説明する。図９は実施の形態２に係る超音波診断装置の動作を示すフローチャートである。なお、図２と同じ動作については同じステップ番号を付し説明を省略する。
実施の形態２に係る超音波診断装置は、ステップＳ１４０において、伝播解析部がステップＳ７０で行った伝播解析の結果を弾性画像記憶部に記憶する。さらに、超音波信号取得部は、ステップＳ１２０において、同一の着目領域を用いて以前に行われたステップＳ１４０で記憶された伝播解析の結果を取得し、ステップＳ１３０において、伝播解析の結果からせん断波の速度を推定し、せん断波の速度が速い期間では検出波の送信間隔を短く、せん断波の速度が遅い期間では検出波の送信間隔を長くなるように、検出波の送信タイミングを決定する。
以下、ステップＳ１３０について詳細に説明する。
まず、超音波信号取得部は、弾性画像記憶部から、すでに取得済みの、着目領域に係る弾性画像を取得する。次に、弾性画像に基づいて、検出波の送信タイミングを決定する。図１０は、実施の形態２に係る検出タイミングの決定方法を示す模式図である。図１０において、弾性画像５３０は、着目領域の弾性画像である。ここで、内包物５４２、５４３は周辺組織より硬い。なお、弾性画像を取得した際のプッシュパルスの焦点は焦点５４１である。速度グラフ５２０は、プッシュパルスの焦点５４１を含む一定の深さ上の直線５３１における、せん断波の速度を示すグラフである。速度グラフ５２０に示すように、せん断波の速度は、内包物５４２に対応する区間５２２、内包物５４３に対応する区間５２４では速く、これに対し、その他の区間５２１、５２３、５２５では遅い。そこで、せん断波の速度は弾性画像５３０と同一であるものと推定し、せん断波の速度が速い期間については検出波を時間的に密に送信し、それ以外の期間については検出波の送信間隔を広くする。すなわち、生成される検出波の送信プロファイル５１０は、以下のようになる。せん断波の速度の速い区間５２２に対応する期間Ｐ１１、区間５２４に対応する期間Ｐ１２では、検出波の送信間隔を短くする。これに対し、せん断波速度の遅い区間５２１に対応する期間Ｐ２１、区間５２３に対応する期間Ｐ２２、区間５２５に対応する期間Ｐ２３では、検出波の送信間隔を長くする。より具体的には、せん断波の推定速度が所定の閾値（例えば、５ｍ／ｓ）以上であれば、検出波の送信間隔を１００μｓとし、閾値未満であれば、２００μｓとする。
なお、さらに、せん断波が着目領域の端に達する時刻ｔ＝ｔ1より後の時刻については、検出波を送信しない、としてもよい。このようにすることで、せん断波が着目領域外に伝播した後における不必要な検出波の送受信を抑止することができ、検出波の送受信回数をさらに削減することができる。
＜補足＞
上述の例では、着目領域に係る弾性画像と新たに行う測定との間でプッシュパルスの焦点位置が同一である場合について説明した。しかしながら、例えば、着目領域に係る弾性画像と新たに行う測定との間でプッシュパルスの焦点位置が同一でなくてもよい。この場合、以下のような処理を行うことができる。
例えば、着目領域に係る弾性画像と新たに行う測定との間でプッシュパルスの焦点位置について、深さ（ｙ軸方向）のみが同一で素子列方向（ｘ軸方向）が異なる場合、速度グラフにおいて新たなプッシュパルスの焦点位置におけるｔ座標を、検出波送信プロファイルのｔ＝０（プッシュパルスの送信直後）とする。具体的には、せん断波の進行方向が右側（ｙ軸の正の方向）で、新たなプッシュパルスの焦点位置が元のプッシュパルスの焦点より右側（ｙ軸の正の方向）であるとした場合、せん断波が新たなプッシュパルスの焦点位置を通過する時刻における時刻を検出波送信プロファイルのｔ＝０（プッシュパルスの送信直後）とするように、検出波送信プロファイル上のｔ＝０をｔ軸の正の方向に移動させる。同様に、プッシュパルスの焦点位置がせん断波の進行方向とは逆方向に移動する場合には、せん断波が新たなプッシュパルスの焦点位置を通過する時刻における時刻を検出波送信プロファイルのｔ＝０とするように、検出波送信プロファイル上のｔ＝０をｔ軸の負の方向に移動させる。
また、着目領域に係る弾性画像と新たに行う測定との間でプッシュパルスの焦点位置について、深さが異なる場合は、新たなプッシュパルスの焦点位置の深さにおけるせん断波の速度を用いることができる。なお、弾性画像に係るプッシュパルスの焦点位置の深さにおけるせん断波の速度を、新たなプッシュパルスの焦点位置の深さにおけるせん断波の速度の代用としてもよい。
また、せん断波の速度はプッシュパルスの焦点位置の深さに限らず、例えば、着目領域の中心点を含む深さ、皮膚表面から所定の深さ、などの特定の深さにおけるせん断波の伝播速度を用いてもよい。例えば、上述の例において、直線５３１を、深さ方向に着目領域の中心となる深さに設けてもよい。また、弾性画像５３０に単一の深さの直線を複数引き、複数の直線それぞれにおけるせん断波の伝播速度を検出し、それらを統合してせん断波の伝播速度を推定するとしてもよい。単一の深さの直線としては、例えば、着目領域を深さ方向に１：３に分割する直線、１：１に分割する直線、３：１に分割する直線、の組み合わせを用いることができる。統合する手段としては、例えば、せん断波の素子列方向（ｘ方向）における同一の座標に対し、複数のｙ座標におけるせん断波の伝播速度の平均値を用いるとしてもよいし、平均値に代えて、中央値、最大値等他の代表値を用いるとしてもよい。
なお、上述の説明では、弾性画像に基づいてせん断波の速度を算出したが、伝播解析の結果に基づいてせん断波の速度を算出してもよい。さらに、例えば、プッシュパルスの送信後最初に変位の検出された時刻を波面の到達とみなす方法により、せん断波の速度を算出してもよい。
なお、着目領域におけるせん断波の速度を取得していない場合、超音波信号取得部は、実施の形態１に係る方法で、検出波の送信プロファイルを設定することができる。
＜小括＞
以上のように検出波の送信プロファイルを設定することで、せん断波の速度の速い期間については検出波を時間的に密に送信し、せん断波の速度の精度低下を抑止することができる。また、せん断波の速度の遅い期間では検出波の送信間隔が長くなるため、超音波探触子の発熱を抑止することができる。また、常に検出波を時間的に密に送信する場合と比較して検出波の送信回数も少なくなるため、受信信号の量を削減して演算量を削減することができる。さらに、演算量の削減および不必要な検出波の送受信の抑止により、プッシュパルスを送信してから弾性画像を生成するまでの時間を短縮することが可能となり、超音波診断装置のリアルタイム性を向上しユーザビリティの向上に奏功する。
≪実施の形態３≫
実施の形態１、２では、せん断波の速度の算出において、プッシュパルスの送信を１度だけ行う場合について説明した。
これに対し、実施の形態３では、せん断波の速度算出をプッシュパルスの複数送信で行うことに特徴がある。
＜動作＞
実施の形態３に係る超音波診断装置の動作について説明する。図１１は実施の形態３に係る超音波診断装置の動作を示すフローチャートである。なお、図２、図９と同じ動作については同じステップ番号を付し説明を省略する。
着目領域の設定（ステップＳ１０）から検出波の送信プロファイルの決定（ステップＳ１３０）については、実施の形態２と同一であるので説明を省略する。
次に、プッシュパルスの送信と後に続く検出波の送受信について説明する（ステップＳ２１０〜２５０）。ここでは、プッシュパルスの焦点を４つ用い、各プッシュパルスの送出後に検出波の送受信を行ったのち伝播解析を行う。４つの焦点位置は、ここでは、着目領域を振動子の列方向に４分割した小領域の１つの中央に設定する。そして、各プッシュパルスの送信後、ステップＳ１３０で決定した送信プロファイルに従って検出波の送受信を行い、変位の検出とせん断波の伝播解析を行う。
以下、せん断波の伝播解析結果の合成（ステップＳ２７０）について、図１２の模式図を用いて説明する。図１２（ａ−１）〜（ａ−４）は、それぞれ、第１〜第４のプッシュパルスに対応する、着目領域とせん断波の速度分布、および、プッシュパルスの焦点との位置関係を示している。例えば、図１２（ａ−１）における速度分布図６１０では、着目領域６１１の左端近くにプッシュパルスの焦点位置６１２が存在し、せん断波の伝播速度が高い領域６１３が検出されている。同様に、図１２（ａ−２）における速度分布図６２０では、着目領域６２１の左寄りにプッシュパルスの焦点位置６２２が存在し、せん断波の伝播速度が高い領域６２３が検出されている。また、図１２（ａ−３）における速度分布図６３０では、着目領域６３１の右寄りにプッシュパルスの焦点位置６３２が存在し、せん断波の伝播速度が高い領域４３３が検出されている。同様に、図１２（ａ−４）における速度分布図６４０では、着目領域６４１の右端近くにプッシュパルスの焦点位置６４２が存在し、せん断波の伝播速度が高い領域６４３が検出されている。なお、せん断波の伝播速度が高い領域６１３、６２３、６３３、６４３は、実際には１つの組織に対応したものであるが、せん断波の減衰により、プッシュパルスの焦点位置より近い部分はその他の領域との境界が明確である半面、プッシュパルスの焦点位置より遠い部分についてはその他の領域との境界が不明確となることがある。これらを統合し、図１２（ｂ）に示すような１つの速度分布図６５０を作成する。具体的には、各座標について、速度分布図６１０、６２０、６３０、６４０のそれぞれからせん断波の速度を取得し、代表値を算出する。代表値の算出方法としては、例えば、重みづけ平均、最大値などを用いてもよいし、または、無効なデータ（例えば、速度が取得できておらず値がない、他の速度分布図から取得した速度のいずれとも差異が大きい、など）を除いた平均値を用いるとしてもよい。重みづけ平均を用いる場合、重みづけ係数ａiとしては、プッシュパルスの焦点位置に近づくほど大きく、プッシュパルスの焦点位置から遠ざかるほど小さくなる値を用いることができる。これは、プッシュパルスの焦点距離から近いほどせん断波のエネルギーが大きく、せん断波速度の精度が高いことが期待できるからである。重みづけ係数ａiの一例として、図１２（ｃ）の係数６６１に示すように、プッシュパルスの焦点位置（ｘ座標がｘf）とｘ座標の差が一定以上であれば０としてもよいし、図１２（ｃ）の係数６６２〜６６４に示すように、プッシュパルスの焦点位置（ｘ座標がｘf）とｘ座標の差が小さいほど大きくなるような任意の関数を使用してよい。このようにすることで、せん断波の伝播速度が高い領域６５３の全貌が確認できるような速度分布図６５０を生成することができる。
せん断波の伝播解析結果の合成後、伝播解析部は、合成後の速度分布図６５０を保持し（ステップＳ１４０）、弾性画像を生成し制御部がこれを表示部に表示する（ステップＳ８０）。
＜小括＞
上記構成により、着目領域内の各観測点と最近接のプッシュパルスの焦点との間の距離を短縮することができ、精度の高いせん断波の速度測定を行うことが可能となる。さらに、せん断波の速度の速い期間についてはせん断波を時間的に密に送信し、せん断波の速度の精度低下を抑止することができる。また、せん断波の速度の遅い期間では検出波の送信間隔が長くなるため、超音波探触子の発熱を抑止することができる。また、検出波の送信回数も少なくなるため、受信信号の量を削減して演算量を削減することができる。さらに、演算量の削減により、プッシュパルスを送信してから弾性画像を生成するまでの時間を短縮することが可能となり、超音波診断装置のリアルタイム性を向上しユーザビリティの向上に奏功する。
≪変形例≫
実施の形態３では、複数のプッシュパルスに係るせん断波の伝播解析結果を合成し、次回以降の測定において、合成後のせん断波の伝播解析結果を用いて検出波の送受信プロファイルを最適化する場合について説明した。
本変形例ではこれに対し、複数のプッシュパルスを用いた測定を行ってその結果を合成する場合に、２回目以降のプッシュパルスの送信後の検出波の送信において、１回目以降のプッシュパルスに係るせん断波の伝播解析結果を用いる点に特徴がある。
＜動作＞
変形例に係る超音波診断装置の動作について説明する。図１３は変形例に係る超音波診断装置の動作を示すフローチャートである。なお、図２、図９、図１１と同じ動作については同じステップ番号を付し説明を省略する。
本実施例では実施の形態３と異なり、プッシュパルスの送信ごとに、ステップＳ１４０で記憶された伝播解析の結果を取得し（ステップＳ１２０）、検出波の送信タイミングを決定する（ステップＳ３３０）。また、伝播解析部は、１回のプッシュパルスに係る伝播解析が完了するごとに、伝播解析結果を保持する（Ｓ１４０）。したがって、ｉ≧２におけるステップＳ１２０では、ｉ＝ｉ−１におけるステップＳ１４０で記憶されたせん断波の伝播解析結果が取得されることになる。そして、すべてのプッシュパルスに係る伝播解析の完了後、せん断波の伝播解析結果の合成（ステップＳ２７０）と弾性画像の生成および表示（ステップＳ８０）を行う。
以下、ステップＳ３３０について詳細に説明する。
まず、超音波信号取得部は、弾性画像記憶部から、直前に送出したプッシュパルスに係る、着目領域に係る弾性画像を取得する。次に、弾性画像に基づいて、検出波の送信タイミングを決定する。図１４は、変形例に係る検出タイミングの決定方法を示す模式図である。図１４の弾性画像７００は、第１回目のプッシュパルスに係る着目領域の弾性画像である。ここで、内包物７０３は周辺組織より硬く、内包物７０４は内包物７０３よりさらに硬い。なお、第１回目のプッシュパルスの焦点は焦点７０１である。速度グラフ７１０は、焦点７０１および第２回目のプッシュパルスの焦点７０２を含む一定の深さ上の直線７１０における、せん断波の速度を示すグラフである。速度グラフ７１０に示すように、せん断波の速度は、内包物７０３に対応する区間７２２では速く、内包物７０４に対応する区間７１５ではさらに速く、その他の区間７１１、７１４、７１６では遅い。そこで、せん断波の速度が速い期間については検出波を時間的に密に送信し、それ以外の期間については検出波の送信間隔を広くする。ここで、焦点の位置が異なるため、次の方法によりせん断波速度の予測を行う。まず、プッシュパルスの送信終了タイミングにおいてｔ＝０となるように、せん断波が焦点７０２に到達する時刻ｔ2を、検出波の送信プロファイル７２０におけるｔ＝０（プッシュパルスの送信終了時）とする。次に、焦点の位置が異なることにより、せん断波の進行向きが逆転する領域について、データの補正を行う。すなわち、焦点７０１と焦点７０２の間の領域では、焦点７０１を発生源とするせん断波が右向き（ｘ軸の正の方向）に伝播するところ、焦点７０２を発生源とするせん断波は左向き（ｘ軸の負の方向）に伝播することになる。そこで、当該領域について、せん断波の伝播速度の大きさ（絶対値）は進行方向に依存しないと仮定する。すなわち、焦点７０２から見て焦点７０１側の領域（焦点７０２より左側の領域）に対応する区間７１１、７１２、７１３について、速度グラフ７１０における点７２２を中心とした点対称となるように折り返し、それらを区間７１１’、７１２’、７１３’とする。そして、せん断波の速度の速い区間に対応する期間では検出波の送信間隔を短くし、せん断波速度の遅い区間に対応する期間では、検出波の送信間隔を長くする。したがって、区間７１５および区間７１２’のうち少なくとも一方に対応する期間Ｐ４１では検出波の送信間隔を短くし、それ以外の期間Ｐ３１、Ｐ３２では検出波の送信間隔を長くする。また、折り返し区間７１１’、７１２’、７１３’を含めてせん断波が着目領域内に存在しなくなる時刻ｔ＝ｔ3以降については、検出波を送信しない。
＜小括＞
上記構成により、実施の形態３の構成に加えて、プッシュパルスを複数回送信する場合に、着目領域に関するせん断波の速度を事前に得ていなくても、２回目以降のプッシュパルスの送信について、検出波の送受信プロファイルを最適化することが可能となる。
≪実施の形態に係るその他の変形例≫
（１）実施の形態１では、第１期間Ｐ１は３．８ｍｓ、送信間隔Ｉ１は１００μｓ、第２期間Ｐ２は３．８ｍｓ、送信間隔Ｉ２は２００μｓ、第３期間Ｐ３は１１．２ｍｓ、送信間隔Ｉ３は４００μｓ、第４期間Ｐ４は１２ｍｓ、送信間隔Ｉ４は８００μｓである場合について説明したが、本発明は必ずしもこの場合に限られない。具体的には、Ｉ１＜Ｉ２＜Ｉ３＜Ｉ４であればよく、また、Ｐ１、Ｐ２、Ｐ３、Ｐ４は任意に設定してよい。例えば、Ｉ１：Ｉ２：Ｉ３：Ｉ４＝１：２：３：４としてもよいし、Ｉ１：Ｉ２：Ｉ３：Ｉ４＝１：３：９：２７としてもよい。なお、Ｉ１は最も速い伝播速度ｖ1を有するせん断波の伝播解析を精度良く行える間隔であることが好ましく、Ｐ１は伝播速度ｖ1のせん断波が発生してから着目領域外に伝播する時間であることが好ましい。また、Ｉ２は次に速い伝播速度ｖ2を有するせん断波の伝播解析を精度良く行える間隔であることが好ましく、Ｐ１＋Ｐ２が伝播速度ｖ2のせん断波が発生してから着目領域外に伝播する時間であることが好ましい。同様に、Ｉ４は最も遅い伝播速度ｖ4を有するせん断波の伝播解析を精度良く行える間隔であることが好ましく、Ｐ１＋Ｐ２＋Ｐ３＋Ｐ４が伝播速度ｖ4のせん断波が発生してから着目領域外に伝播する時間であることが好ましい。
なお、実施の形態及び上述の変形例では期間および送信間隔を４種類ずつ設定するものとしたが、２種類や３種類であってもよいし、５種類以上を用いてもよい。
また、着目領域の大きさや観測すべき部位に応じて、プロファイルを調整、または、あらかじめ用意した複数のプロファイルから１つを選択する、としてもよい。例えば、着目領域の素子列方向の幅が広い場合には、送信間隔Ｉ２〜Ｉ４をそれぞれ長くする、などの手法で、プッシュパルスの送信回数を削減することができる。また、例えば、平均的な弾性率が大きい着目領域については、送信間隔Ｉ１〜Ｉ４をそれぞれ短くし、平均的な弾性率が小さい着目領域については、送信間隔Ｉ１〜Ｉ４をそれぞれ長くする、としてもよい。
（２）各実施の形態及び変形例では、検出波の送信プロファイルを、送信間隔と期間との組み合わせで示しているが、例えば、検出波の送信プロファイルを、送信間隔と送信回数の組み合わせで示してもよい。例えば、第１期間Ｐ１は３．８ｍｓ、送信間隔Ｉ１は１００μｓ、第２期間Ｐ２は３．８ｍｓ、送信間隔Ｉ２は２００μｓ、第３期間Ｐ３は１１．２ｍｓ、送信間隔Ｉ３は４００μｓ、第４期間Ｐ４は１２ｍｓ、送信間隔Ｉ４は８００μｓである場合について、第１送信回数Ｃ１は３８回、送信間隔Ｉ１は１００μｓ、第２送信回数Ｃ２は１９回、送信間隔Ｉ２は２００μｓ、第３送信回数Ｃ３は２８回、送信間隔Ｉ３は４００μｓ、第４送信回数Ｃ４は１５回、送信間隔Ｉ４は８００μｓとして示してもよい。この２つの送信プロファイルは本質的には同一であり、表現を変えているに過ぎないからである。なお、検出波の送信プロファイルを、送信間隔と送信回数との組み合わせとして示してもよいのは勿論である。
（３）各実施の形態及び変形例では、検出波の送信プロファイルは、期間ごとに送信間隔を定め、段階的に送信間隔を変更しているが、送信間隔を連続的に変化させてもよい。例えば、実施の形態１の場合において、送信間隔を、プッシュパルスの送信から直前の検出波送信までの時間に依存して増加する関数、例えば、一次関数として定めてもよい。同様に、例えば、実施の形態２の場合において、送信間隔を、せん断波の推定速度が高いほど小さくなる関数、例えば、反比例として定めてもよい。
（４）実施の形態２、３および変形例では、せん断波の波面の推定位置が着目領域外である時刻には検出波を送信しないこととしたが、さらに、以下のような処理を行ってもよい。例えば、検出波の総送信回数が閾値より小さくなった場合には、第１期間Ｐ１を延長する、または、送信間隔Ｉ１を短縮するなどの処理を行ってもよい。具体的には、検出波の総送信回数が８０回に満たない場合には、総送信回数が１００となるように、第１期間Ｐ１を延長してもよい。このようにすることで、検出波の総送信回数が少なく超音波探触子の発熱を考慮する必要がない場合、せん断波の伝播解析の精度をさらに高めることが可能となる。
（５）各実施の形態及び変形例では、せん断波の伝播解析を、変位領域の抽出、細線化処理、空間フィルタリング、時間フィルタリングの手順によって行う場合について説明したが、各場所の変位のピーク時刻検出、時間フィルタリング、空間フィルタリングの手順で行ってもよい。
また、せん断波の伝播解析を、上述したような手法ではなく、単に、着目領域内の各場所において、プッシュパルスの送信後最初に変位が観測された時刻を、当該場所にせん断波の波面が到達したとみなして伝播解析を行ってもよい。
（６）実施の形態３の変形例では、ｉ＝１のプッシュパルスに係る検出波の送信プロファイルについては実施の形態１と同様の方法で決定するとしたが、実施の形態３と同様の方法で検出波の送信プロファイルを決定してもよい。また、ｉ≧２のプッシュパルスに係る検出波の送信プロファイルについては直前のプッシュパルスに係るせん断波の伝播解析結果を用いるとしたが、例えば、２以上前のプッシュパルスに係るせん断波の伝播解析結果を用いるとしてもよい。
（７）実施の形態および各変形例では、超音波診断装置は表示部３と接続される構成であるとしたが、本発明は必ずしもこの場合に限られない。例えば、超音波診断装置１は表示部３を内蔵しているとしてもよい。または、超音波診断装置１は表示部３と接続されず、伝播解析部１５が生成し弾性画像記憶部１９に記憶されている弾性画像を他の記憶媒体に保存、あるいは、ネットワークを通じて他の装置に出力するとしてもよい。
また、同様に、超音波診断装置は超音波探触子２を内蔵しているとしてもよいし、あるいは、超音波探触子２が超音波信号取得部１３を備え、超音波信号取得部１３を有さない超音波診断装置が超音波探触子２から音響線信号を取得するとしてもよい。
（８）実施の形態および各変形例に係る超音波診断装置は、その構成要素の全部又は一部を、１チップ又は複数チップの集積回路で実現してもよいし、コンピュータのプログラムで実現してもよいし、その他どのような形態で実施してもよい。例えば、伝播解析部と評価部とを１チップで実現してもよいし、超音波信号取得部のみを１チップで実現し、変位検出部等を別のチップで実現してもよい。
集積回路で実現する場合、典型的には、ＬＳＩ（Ｌａｒｇｅ Ｓｃａｌｅ Ｉｎｔｅｇｒａｔｉｏｎ）として実現される。ここでは、ＬＳＩとしたが、集積度の違いにより、ＩＣ、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩと呼称されることもある。
また、集積回路化の手法はＬＳＩに限るものではなく、専用回路、又は汎用プロセッサで実現してもよい。ＬＳＩ製造後に、プログラムすることが可能なＦＰＧＡ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）や、ＬＳＩ内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサを利用してもよい。
さらには、半導体技術の進歩、又は派生する別技術により、ＬＳＩに置き換わる集積回路化の技術が登場すれば、当然その技術を用いて機能ブロックの集積化を行ってもよい。
また、各実施の形態および各変形例に係る超音波診断装置は、記憶媒体に書き込まれたプログラムと、プログラムを読み込んで実行するコンピュータとで実現されてもよい。記憶媒体は、メモリカード、ＣＤ−ＲＯＭなどいかなる記録媒体であってもよい。また、本発明に係る超音波診断装置は、ネットワークを経由してダウンロードされるプログラムと、プログラムをネットワークからダウンロードして実行するコンピュータとで実現されてもよい。
（９）以上で説明した実施の形態は、いずれも本発明の好ましい一具体例を示すものである。実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置及び接続形態、工程、工程の順序などは一例であり、本発明を限定する主旨ではない。また、実施の形態における構成要素のうち、本発明の最上位概念を示す独立請求項に記載されていない工程については、より好ましい形態を構成する任意の構成要素として説明される。
また、発明の理解の容易のため、上記各実施の形態で挙げた各図の構成要素の縮尺は実際のものと異なる場合がある。また本発明は上記各実施の形態の記載によって限定されるものではなく、本発明の要旨を逸脱しない範囲において適宜変更可能である。
さらに、超音波診断装置においては基板上に回路部品、リード線等の部材も存在するが、電気的配線、電気回路について当該技術分野における通常の知識に基づいて様々な態様を実施可能であり、本発明の説明として直接的には無関係のため、説明を省略している。尚、上記示した各図は模式図であり、必ずしも厳密に図示したものではない。
≪まとめ≫
（１）実施の形態に係る超音波診断装置は、超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して特定部位にある組織を物理的に押圧した後検出波の送受信を被検体内に対して繰り返し行うことで、押圧された特定部位の組織を振動源とするせん断波の伝播速度を検出する超音波診断装置であって、プッシュパルスを送信するプッシュパルス送信部と、前記プッシュパルスに続き平面波を検出波として被検体内の着目領域へ複数回送信し、検出波に対応する被検体からの反射検出波を受信して複数の受信信号を時系列に生成する検出波送受信部と、前記反射検出波の受信時刻のそれぞれにおける、前記プッシュパルスに起因するせん断波による前記着目領域内の組織の変位を前記複数の受信信号からそれぞれ検出する変位検出部と、前記変位検出部が検出した変位に基づいて前記着目領域内におけるせん断波の伝播速度を検出するせん断波解析部とを備え、前記検出波送受信部は、プッシュパルスに続く一連の検出波の送信において、連続する２つの検出波の送信間隔が、第１の送信間隔と、前記第１の送信間隔より長い第２の送信間隔とを少なくとも含むように、検出波を送信することを特徴とする。
また、実施の形態に係る超音波信号処理装置は、超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して特定部位にある組織を物理的に押圧した後検出波の送受信を被検体内に対して繰り返し行うことで、押圧された特定部位の組織を振動源とするせん断波の伝播速度を検出する超音波信号処理方法であって、プッシュパルスを送信し、前記プッシュパルスに続き平面波を検出波として被検体内の着目領域へ複数回送信し、検出波に対応する被検体からの反射検出波を受信して複数の受信信号を時系列に生成し、前記反射検出波の受信時刻のそれぞれにおける、前記プッシュパルスに起因するせん断波による前記着目領域内の組織の変位を前記複数の受信信号からそれぞれ検出し、検出した前記変位に基づいて前記着目領域内におけるせん断波の伝播速度を検出する処理を含み、プッシュパルスに続く一連の検出波の送信において、連続する２つの検出波の送信間隔が、第１の送信間隔と、前記第１の送信間隔より長い第２の送信間隔とを少なくとも含むように、検出波を送信することを特徴とする。
本開示によれば、上記構成により、検出波の送信間隔を調整することで、せん断波の伝播速度の精度低下を抑止しながら検出波の送受信頻度を低下させることができる。
（２）また、上記（１）の超音波診断装置は、前記検出波送受信部は、前記プッシュパルスの送信時刻からの経過時間が所定時間以下である場合は前記第１の送信間隔で検出波を送信し、前記プッシュパルスの送信時刻からの経過時間が所定時間を超えた場合は前記第２の送信間隔で検出波を送信する、としてもよい。
（３）また、上記（１）の超音波診断装置は、前記検出波送受信部は、前記プッシュパルスに続いてすでに送信した検出波の回数が所定回数以下である場合は前記第１の送信間隔で検出波を送信し、前記プッシュパルスに続いてすでに送信した検出波の回数が所定回数を超えた場合は前記第２の送信間隔で検出波を送信する、としてもよい。
上記（２）または（３）の構成により、速度の速いせん断波を観測すべき、プッシュパルスの送信時刻と近い期間においては検出波の送信間隔を短くし、速度の遅いせん断波を観測すべき、プッシュパルスの送信時刻から遠い期間においては検出波の送信間隔を長くする。したがって、せん断波の伝播速度にかかわらず、せん断波の速度の精度低下を抑止することができる。
（４）また、上記（１）の超音波診断装置は、前記せん断波解析部は、前記着目領域内におけるせん断波の伝播速度を保持し、前記検出波送受信部は、前記プッシュパルスの送信前において前記せん断波解析部が保持しているせん断波の伝播速度からせん断波の伝播速度を推定し、前記プッシュパルスの送信後の期間のうち、推定されたせん断波の伝播速度が所定速度以上である期間においては前記第１の送信間隔で、所定速度未満である期間においては前記第２の送信間隔で、それぞれ検出波を送信する、としてもよい。
上記構成により、せん断波の速度の速い期間については検出波を時間的に密に送信し、せん断波の速度の精度低下を抑止することができる。また、せん断波の速度の遅い期間では検出波の送信間隔が長くなるため、検出波の送信頻度を低下させることができる。
（５）また、上記（４）の超音波診断装置は、前記プッシュパルス送信部は、複数のプッシュパルスを送信し、前記検出波送受信部は、プッシュパルスのそれぞれに続いて一連の検出波の送受信を行い、前記変位検出部は、プッシュパルスごとに、プッシュパルスに続く一連の検出波の送受信で得た前記複数の受信信号に基づいて前記着目領域内の組織の変位を検出し、前記せん断波解析部は、プッシュパルスごとに、変位に基づいて前記着目領域内のせん断波の伝播速度を検出した後、検出した伝播速度を合成することで前記着目領域内のせん断波の合成伝播速度を検出し、前記せん断波の合成伝播速度を保持する、としてもよい。
上記構成により、せん断波の伝播速度の合成によってせん断波の伝播速度の解析結果の精度を高めることができる。また、精度の高いせん断波の伝播速度の解析結果に基づいて、検出波の送信方法を最適化することが可能となる。
（６）また、上記（４）の超音波診断装置は、前記プッシュパルス送信部は、複数のプッシュパルスを送信し、前記検出波送受信部は、プッシュパルスのそれぞれに続いて一連の検出波の送受信を行い、前記変位検出部は、プッシュパルスごとに、プッシュパルスに続く一連の検出波の送受信で得た前記複数の受信信号に基づいて前記着目領域内の組織の変位を検出し、前記せん断波解析部は、プッシュパルスごとに、変位に基づいて前記着目領域内のせん断波の伝播速度を検出した後、前記せん断波の伝播速度を保持し、さらに、すべてのプッシュパルスの送信後、検出した伝播速度を合成することで前記着目領域内のせん断波の合成伝播速度を検出する、としてもよい。
上記構成により、せん断波の伝播速度の合成によってせん断波の伝播速度の解析結果の精度を高めることができる。また、２回目以降のプッシュパルスの送信に続く検出波の送信方法を１回目以降のプッシュパルスに起因するせん断波の解析結果に基づいて決定できる。したがって、同一の着目領域に対する合成伝播速度がなくても検出波の送信方法を最適化することが可能となる。
（７）また、上記（４）〜（６）の超音波診断装置は、前記検出波送受信部は、前記プッシュパルスの送信前において前記せん断波解析部が保持しているせん断波の伝播速度からせん断波の位置を推定し、せん断波の推定位置が前記着目領域の外部である場合に検出波を送信しない、としてもよい。
上記構成により、不必要な検出波の送受信を抑止することができ、せん断波の伝播速度の精度に影響を与えることなく検出波の送受信回数を削減できる。
（８）また、上記（４）〜（７）の超音波診断装置は、前記検出波送受信部は、せん断波の速度および／または位置を推定する場合において、保持しているせん断波の伝播速度のうち、前記プッシュパルスが集中する焦点の深さに関する情報を用いる、としてもよい。
上記構成により、せん断波の速度または位置を効率よく推定することができる。
本開示に係る超音波診断装置、および、超音波信号処理方法は、超音波を用いた組織の硬さの測定に有用である。そのため、組織の硬さの測定精度を向上させることが可能となり、医療診断機器等において高い利用可能性を持つ。
１ 超音波診断装置
２ 超音波探触子
３ 表示部
１１ 制御部
１２ せん断波励起部
１３ 超音波信号取得部
１４ 変位検出部
１５ 伝播解析部
１６ 断層画像記憶部
１７ 変位量記憶部
１８ シーケンス保持部
１９ 弾性画像記憶部

Claims (9)

  1. 超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して特定部位にある組織を物理的に押圧した後検出波の送受信を被検体内に対して繰り返し行うことで、押圧された特定部位の組織を振動源とするせん断波の伝播速度を検出する超音波診断装置であって、
    プッシュパルスを送信するプッシュパルス送信部と、
    前記プッシュパルスに続き平面波を検出波として被検体内の着目領域へ複数回送信し、検出波に対応する被検体からの反射検出波を受信して複数の受信信号を時系列に生成する検出波送受信部と、
    前記反射検出波の受信時刻のそれぞれにおける、前記プッシュパルスに起因するせん断波による前記着目領域内の組織の変位を前記複数の受信信号からそれぞれ検出する変位検出部と、
    前記変位検出部が検出した変位に基づいて前記着目領域内におけるせん断波の伝播速度を検出するせん断波解析部と
    を備え、
    前記検出波送受信部は、プッシュパルスに続く一連の検出波の送信において、連続する２つの検出波の送信間隔が、第１の送信間隔と、前記第１の送信間隔より長い第２の送信間隔とを少なくとも含むように、検出波を送信する
    ことを特徴とする超音波診断装置。
  2. 前記検出波送受信部は、前記プッシュパルスの送信時刻からの経過時間が所定時間以下である場合は前記第１の送信間隔で検出波を送信し、前記プッシュパルスの送信時刻からの経過時間が所定時間を超えた場合は前記第２の送信間隔で検出波を送信する
    請求項１に記載の超音波診断装置。
  3. 前記検出波送受信部は、前記プッシュパルスに続いてすでに送信した検出波の回数が所定回数以下である場合は前記第１の送信間隔で検出波を送信し、前記プッシュパルスに続いてすでに送信した検出波の回数が所定回数を超えた場合は前記第２の送信間隔で検出波を送信する
    請求項１に記載の超音波診断装置。
  4. 前記せん断波解析部は、前記着目領域内におけるせん断波の伝播速度を保持し、
    前記検出波送受信部は、前記プッシュパルスの送信前において前記せん断波解析部が保持しているせん断波の伝播速度からせん断波の伝播速度を推定し、前記プッシュパルスの送信後の期間のうち、推定されたせん断波の伝播速度が所定速度以上である期間においては前記第１の送信間隔で、所定速度未満である期間においては前記第２の送信間隔で、それぞれ検出波を送信する
    請求項１に記載の超音波診断装置。
  5. 前記プッシュパルス送信部は、複数のプッシュパルスを送信し、
    前記検出波送受信部は、プッシュパルスのそれぞれに続いて一連の検出波の送受信を行い、
    前記変位検出部は、プッシュパルスごとに、プッシュパルスに続く一連の検出波の送受信で得た前記複数の受信信号に基づいて前記着目領域内の組織の変位を検出し、
    前記せん断波解析部は、プッシュパルスごとに、変位に基づいて前記着目領域内のせん断波の伝播速度を検出した後、検出した伝播速度を合成することで前記着目領域内のせん断波の合成伝播速度を検出し、前記せん断波の合成伝播速度を保持する
    ことを特徴とする請求項４に記載の超音波診断装置。
  6. 前記プッシュパルス送信部は、複数のプッシュパルスを送信し、
    前記検出波送受信部は、プッシュパルスのそれぞれに続いて一連の検出波の送受信を行い、
    前記変位検出部は、プッシュパルスごとに、プッシュパルスに続く一連の検出波の送受信で得た前記複数の受信信号に基づいて前記着目領域内の組織の変位を検出し、
    前記せん断波解析部は、プッシュパルスごとに、変位に基づいて前記着目領域内のせん断波の伝播速度を検出した後、前記せん断波の伝播速度を保持し、さらに、すべてのプッシュパルスの送信後、検出した伝播速度を合成することで前記着目領域内のせん断波の合成伝播速度を検出する
    ことを特徴とする請求項４に記載の超音波診断装置。
  7. 前記検出波送受信部は、前記プッシュパルスの送信前において前記せん断波解析部が保持しているせん断波の伝播速度からせん断波の位置を推定し、せん断波の推定位置が前記着目領域の外部である場合に検出波を送信しない
    ことを特徴とする請求項４から６のいずれか１項に記載の超音波診断装置。
  8. 前記検出波送受信部は、せん断波の速度および／または位置を推定する場合において、保持しているせん断波の伝播速度のうち、前記プッシュパルスが集中する焦点の深さに関する情報を用いる
    ことを特徴とする請求項４から７のいずれか１項に記載の超音波診断装置。
  9. 超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して特定部位にある組織を物理的に押圧した後検出波の送受信を被検体内に対して繰り返し行うことで、押圧された特定部位の組織を振動源とするせん断波の伝播速度を検出する超音波信号処理方法であって、
    プッシュパルスを送信し、
    前記プッシュパルスに続き平面波を検出波として被検体内の着目領域へ複数回送信し、検出波に対応する被検体からの反射検出波を受信して複数の受信信号を時系列に生成し、
    前記反射検出波の受信時刻のそれぞれにおける、前記プッシュパルスに起因するせん断波による前記着目領域内の組織の変位を前記複数の受信信号からそれぞれ検出し、
    検出した前記変位に基づいて前記着目領域内におけるせん断波の伝播速度を検出する
    処理を含み、
    プッシュパルスに続く一連の検出波の送信において、連続する２つの検出波の送信間隔が、第１の送信間隔と、前記第１の送信間隔より長い第２の送信間隔とを少なくとも含むように、検出波を送信する
    ことを特徴とする超音波信号処理方法。
JP2016112802A 2016-06-06 2016-06-06 超音波診断装置及び超音波信号処理方法 Active JP6604272B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2016112802A JP6604272B2 (ja) 2016-06-06 2016-06-06 超音波診断装置及び超音波信号処理方法
US15/612,295 US11399805B2 (en) 2016-06-06 2017-06-02 Ultrasound diagnostic device and ultrasound signal processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016112802A JP6604272B2 (ja) 2016-06-06 2016-06-06 超音波診断装置及び超音波信号処理方法

Publications (2)

Publication Number Publication Date
JP2017217147A JP2017217147A (ja) 2017-12-14
JP6604272B2 true JP6604272B2 (ja) 2019-11-13

Family

ID=60482495

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016112802A Active JP6604272B2 (ja) 2016-06-06 2016-06-06 超音波診断装置及び超音波信号処理方法

Country Status (2)

Country Link
US (1) US11399805B2 (ja)
JP (1) JP6604272B2 (ja)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11129598B2 (en) * 2017-12-13 2021-09-28 Siemens Medical Solutions Usa, Inc. Calibration for ARFI imaging
US11576654B2 (en) * 2017-12-21 2023-02-14 Samsung Medison Co., Ltd. Ultrasound diagnosis apparatus for measuring and displaying elasticity of object and method of operating the same
WO2019192970A1 (en) * 2018-04-02 2019-10-10 Koninklijke Philips N.V. Ultrasonic shear wave imaging with improved accuracy and reliability
CN111050661B (zh) * 2018-05-15 2023-08-11 深圳迈瑞生物医疗电子股份有限公司 一种剪切波弹性测量方法及剪切波弹性成像***
CN108852416B (zh) 2018-05-17 2021-03-30 深圳开立生物医疗科技股份有限公司 一种剪切波传播速度的确定方法及装置
JP7354632B2 (ja) * 2019-07-12 2023-10-03 コニカミノルタ株式会社 超音波診断装置、および、超音波診断装置の制御方法
JP7304230B2 (ja) * 2019-07-26 2023-07-06 富士フイルムヘルスケア株式会社 超音波撮像装置

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6053860B2 (ja) * 1978-10-31 1985-11-27 ウエスト電気株式会社 電子閃光装置
US7547283B2 (en) * 2000-11-28 2009-06-16 Physiosonics, Inc. Methods for determining intracranial pressure non-invasively
US7004906B1 (en) * 2004-07-26 2006-02-28 Siemens Medical Solutions Usa, Inc. Contrast agent imaging with agent specific ultrasound detection
US7887484B2 (en) * 2005-08-12 2011-02-15 Siemens Medical Solutions Usa, Inc. Automatic velocity scale identification for medical diagnostic ultrasound
US8123692B2 (en) * 2006-12-06 2012-02-28 General Electric Company Apparatus, system, and method for adaptively controlling a frame interval between ultrasound scanning frames for an ultrasound elasticity imaging scan
US9364194B2 (en) * 2008-09-18 2016-06-14 General Electric Company Systems and methods for detecting regions of altered stiffness
US9554770B2 (en) * 2008-09-29 2017-01-31 Siemens Medical Solutions Usa, Inc. High pulse repetition frequency for detection of tissue mechanical property with ultrasound
EP2474266A4 (en) * 2009-09-04 2014-11-05 Hitachi Medical Corp ULTRASOUND DIAGNOSTIC TOOL
US8961418B2 (en) * 2010-10-06 2015-02-24 Siemens Medical Solutions Usa, Inc. Solving for shear wave information in medical ultrasound imaging
KR101929198B1 (ko) * 2011-02-25 2018-12-14 메이오 파운데이션 포 메디칼 에쥬케이션 앤드 리써치 비집속식 초음파에 의한 초음파 바이브로메트리
US9237878B2 (en) * 2011-04-22 2016-01-19 Mayo Foundation For Medical Education And Research Generation and assessment of shear waves in elasticity imaging
US8801614B2 (en) * 2012-02-10 2014-08-12 Siemens Medical Solutions Usa, Inc. On-axis shear wave characterization with ultrasound
WO2013153857A1 (ja) * 2012-04-13 2013-10-17 日立アロカメディカル株式会社 超音波診断装置及び軌跡表示方法
JP6257892B2 (ja) * 2012-12-10 2018-01-10 東芝メディカルシステムズ株式会社 超音波診断装置及び制御プログラム
WO2015029651A1 (ja) * 2013-08-26 2015-03-05 日立アロカメディカル株式会社 超音波診断装置および弾性評価方法
JP2015128554A (ja) * 2014-01-09 2015-07-16 日立アロカメディカル株式会社 超音波診断装置
JP6386853B2 (ja) * 2014-09-26 2018-09-05 株式会社日立製作所 超音波診断装置
JP6384340B2 (ja) * 2015-01-28 2018-09-05 コニカミノルタ株式会社 超音波診断装置
US11096671B2 (en) * 2015-09-10 2021-08-24 Siemens Medical Solutions Usa, Inc. Sparkle artifact detection in ultrasound color flow

Also Published As

Publication number Publication date
JP2017217147A (ja) 2017-12-14
US20170347990A1 (en) 2017-12-07
US11399805B2 (en) 2022-08-02

Similar Documents

Publication Publication Date Title
JP6604272B2 (ja) 超音波診断装置及び超音波信号処理方法
JP6390516B2 (ja) 超音波診断装置及び超音波診断装置の制御方法
JP6053860B2 (ja) 医療用超音波イメージングシステムにおいて剪断波情報を求めるためのコンピュータプログラム
JP6741012B2 (ja) 超音波診断装置、及び超音波信号処理方法
US20170112471A1 (en) Ultrasound diagnostic apparatus and ultrasound signal processing method
JP2010526626A (ja) 関心領域の粘弾性の平均値を測定するための方法および装置
JP6384340B2 (ja) 超音波診断装置
JP6601320B2 (ja) 超音波診断装置、及び超音波診断装置の制御方法
JP2012081269A5 (ja)
CN110974296B (zh) 一种获取剪切波弹性图像的方法、装置及超声成像***
JP6259953B2 (ja) せん断波を用いる測定に関する超音波システム及び方法
EP3429476B1 (en) Shear wave group velocity estimation using spatiotemporal peaks and amplitude thresholding
CN109259801B (zh) 一种剪切波弹性成像方法及装置
JP6321162B2 (ja) 繊維状媒体のマッピングの方法及びデバイス
JP7354632B2 (ja) 超音波診断装置、および、超音波診断装置の制御方法
JP6672809B2 (ja) 超音波診断装置、及び超音波信号処理方法
JP7027924B2 (ja) 超音波診断装置、及び超音波診断装置の制御方法
JP5491671B2 (ja) 関心領域の粘弾性の平均値を測定するための方法
JP2016087140A (ja) 超音波診断装置
JP7347445B2 (ja) 超音波信号処理装置、超音波診断装置、および、超音波信号処理方法
JP7302651B2 (ja) 超音波信号処理装置、超音波診断装置、超音波信号処理方法、およびプログラム
JP7010082B2 (ja) 超音波診断装置、及び超音波診断装置の制御方法
JP2020199162A (ja) 超音波診断装置、および、超音波診断装置の制御方法
Wüst et al. Towards a Model-based High-resolution Ultrasonic Measurement System for Non-destructive Testing

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20181225

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

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20190913

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190930

R150 Certificate of patent or registration of utility model

Ref document number: 6604272

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150