JP6274517B2 - 超音波診断装置及び超音波画像処理プログラム - Google Patents

超音波診断装置及び超音波画像処理プログラム Download PDF

Info

Publication number
JP6274517B2
JP6274517B2 JP2014048258A JP2014048258A JP6274517B2 JP 6274517 B2 JP6274517 B2 JP 6274517B2 JP 2014048258 A JP2014048258 A JP 2014048258A JP 2014048258 A JP2014048258 A JP 2014048258A JP 6274517 B2 JP6274517 B2 JP 6274517B2
Authority
JP
Japan
Prior art keywords
blood flow
flow information
unit
change
image data
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
JP2014048258A
Other languages
English (en)
Other versions
JP2014198240A (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.)
Canon Medical Systems Corp
Original Assignee
Canon Medical Systems 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 Canon Medical Systems Corp filed Critical Canon Medical Systems Corp
Priority to JP2014048258A priority Critical patent/JP6274517B2/ja
Publication of JP2014198240A publication Critical patent/JP2014198240A/ja
Application granted granted Critical
Publication of JP6274517B2 publication Critical patent/JP6274517B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06—Measuring blood flow
    • 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
    • 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/488—Diagnostic techniques involving Doppler signals
    • 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/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • A61B8/5276—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts due to motion
    • 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/8979—Combined Doppler and pulse-echo imaging systems
    • G01S15/8981—Discriminating between fixed and moving objects or between objects moving at different speeds, e.g. wall clutter filter
    • 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/8979—Combined Doppler and pulse-echo imaging systems
    • G01S15/8988—Colour Doppler 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
    • 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
    • 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/8934—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a dynamic transducer configuration
    • G01S15/8945—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a dynamic transducer configuration using transducers mounted for linear mechanical movement
    • 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/52074—Composite displays, e.g. split-screen displays; Combination of multiple images or of images and alphanumeric tabular information

Description

本発明の実施形態は、超音波診断装置及び超音波画像処理プログラムに関する。
超音波診断装置は、超音波プローブを用いて被検体に超音波を送信してその反射波を受信することにより、被検体の生体情報を取得する装置である。
超音波診断装置には、断層画像（Ｂモード画像）とともに血流画像（ＣＦＭ画像：Ｃｏｌｏｒ Ｆｌｏｗ Ｍａｐｐｉｎｇ画像）を一定の時間間隔毎に生成し、動画情報を生成する装置がある。
血流画像は、ドプラ信号に基づいて生成される。血流画像において、ドプラ信号に基づいて求められた血流情報と色情報とが対応づけられ、色情報により２次元血流像が表示される。それにより、血流は血流画像において可視化される。
特開２０１１−２４８８９号公報
ドプラ信号の中には、臓器の拍動や蠕動などの被検体組織の動きによって生じたクラッタが含まれる場合がある。超音波診断中において、被検体組織の動きが大きいときはドプラ信号に含まれるクラッタが多く、被検体組織の動きが小さいときはドプラ信号に含まれるクラッタが少ない。ドプラに信号含まれたクラッタがそのまま血流画像に描出される場合、クラッタは血流画像のうち部分的な領域に描出される。クラッタが血流画像に描出された像は血流を表す像ではないので、クラッタが描出された像は、誤診の誘因となる場合があった。また、ドプラ信号の中には、超音波プローブの動きによって生じたモーションアーチファクトが含まれる場合がある。超音波プローブの動きは、例えばユーザの手ぶれによって生じる。超音波診断中において、超音波プローブの動きが大きいときは、ドプラ信号に含まれるモーションアーチファクトが多く、超音波プローブの動きが小さいときはドプラ信号に含まれるモーションアーチファクトが少ない。ドプラ信号に含まれたモーションアーチファクトがそのまま血流画像に描出される場合、モーションアーチファクトは血流画像の広範な領域に描出される。このようなモーションアーチファクトの像は、断層画像や血流を表す像をユーザが視認する上で邪魔になっていた。言い換えると、ユーザにとってモーションアーチファクトが描出された像は目障りであった。以上のことから、超音波診断装置について、クラッタやモーションアーチファクトを低減させる要望があった。
本発明が解決しようとする課題は、クラッタやモーションアーチファクトを低減することができる超音波診断装置及び超音波画像処理プログラムを提供することである。
実施形態の超音波診断装置は、断層画像データ生成部と、血流情報生成部と、変化検出部と、条件変更部とを有する。断層画像データ生成部は、超音波を送信された被検体からの反射波に基づく受信信号を受け、被検体の断層画像データを複数の時間について逐次求める。血流情報生成部は、受信信号を受け、被検体の血流情報に係る速度閾値を含む信号処理条件に基づく処理を受信信号に施して血流情報を求める。変化検出部は、時間が互いに異なる複数の断層画像データを受け、受けた複数の断層画像データの時間軸方向における変化の大きさを検出する。条件変更部は、変化検出部が検出した時間軸方向における変化の大きさに基づいて信号処理条件を変更する。血流情報生成部は、条件変更部が信号処理条件を変更したとき、変更された信号処理条件に基づく処理を受信信号に施して血流情報を求める。
また、実施形態の超音波画像処理プログラムは、超音波を送信された被検体からの反射波に基づく受信信号を処理するコンピュータに、被検体の断層画像データを複数の時間について逐次求めさせ、前記被検体の血流情報に係る速度閾値を含む信号処理条件に基づく処理を受信信号に施して血流情報を求めさせ、時間が互いに異なる複数の断層画像データの時間軸方向における変化の大きさを検出させ、検出された時間軸方向における変化の大きさに基づいて信号処理条件を変更させた後、変更された信号処理条件に基づいて血流情報を求めさせる。
実施形態の超音波診断装置の構成を示すブロック図。 実施形態の超音波診断装置の構成を示すブロック図。 実施形態の超音波診断装置の概略を示す模式図。 実施形態の超音波診断装置の動作を示すフローチャート。 実施形態の超音波診断装置により生成された超音波画像。 実施形態の超音波診断装置により生成された超音波画像に対する比較画像。 実施形態の変形例に係る超音波診断装置の構成を示すブロック図。 実施形態の変形例に係る超音波診断装置の動作を示すフローチャート。 実施形態の変形例に係る超音波診断装置の構成を示すブロック図。
［構成］
図１及び図２は、この実施形態の超音波診断装置１の構成を示すブロック図である。超音波診断装置１は、超音波プローブ１０と、送信部１１と、受信部１２と、断層画像データ生成部１３と、血流情報生成部１４と、変化検出部１５と、条件変更部１６と、画像生成部１７と、表示制御部１８と、表示部１９と、システム制御部２０と、操作部２１とを有する。
（超音波プローブ１０）
超音波プローブ１０には、複数の超音波振動子が走査方向に１列に配置された１次元アレイプローブ、又は、複数の超音波振動子が２次元的に配置された２次元アレイプローブが用いられる。また、走査方向に１列に配置された複数の超音波振動子を、走査方向に直交する揺動方向に揺動させる機械式１次元アレイプローブを用いてもよい。超音波プローブ１０は被検体に超音波を送信し、被検体からの反射波をエコー信号として受信する。
（送信部１１）
送信部１１は、超音波プローブ１０に電気信号を供給して超音波を発生させる。送信部１１は、図示しない送信遅延回路及びパルサ回路を有する。送信遅延回路は、超音波の送信時に遅延を掛けて送信フォーカスを実施する。パルサ回路は、各超音波振動子に対応した経路（チャンネル）の数に応じたパルサを備え、遅延が掛けられた送信タイミングで駆動パルスを発生し、超音波プローブ１０の各超音波振動子に供給する。
（受信部１２）
受信部１２は、超音波プローブ１０が受信したエコー信号を受信する。受信部１２は、超音波プローブ１０が受信したエコー信号を受信し、そのエコー信号に対して遅延処理を行うことにより、アナログのエコー信号を整相された（受信ビームフォームされた）デジタルのデータに変換する。
受信部１２は、例えば図示しないプリアンプ回路と、Ａ／Ｄ変換器と、受信遅延回路と、加算器とを有する。プリアンプ回路は、超音波プローブ１０の各超音波振動子から出力されるエコー信号を受信チャンネルごとに増幅する。Ａ／Ｄ変換器は、増幅されたエコー信号をデジタル信号に変換する。受信遅延回路は、デジタル信号に変換されたエコー信号に、受信指向性を決定するために必要な遅延時間を与える。加算器は、遅延時間が与えられたエコー信号を整相加算する。その整相加算によって、受信指向性に応じた方向からの反射成分が強調される。受信部１２は、加算した信号を直交検波し、断層画像データ生成部１３と血流情報生成部１４とに受信信号として出力する。
（断層画像データ生成部１３）
断層画像データ生成部１３は、超音波を送信された被検体からの反射波に基づく受信信号を受信部１２から受け、被検体の断層画像データを複数の時間について逐次求める。換言すると、断層画像データ生成部１３は、受信部１２から受けた受信信号に基づいて、一定の時間間隔（フレームレート）ごとに断層画像データを求める。ある時間について求められた断層画像データは、求められた時間についての静止画（フレーム）に対応する。断層画像データには、後述する、ログ圧縮部１３０からの出力、エッジ強調部１３１からの出力、断層平滑化処理部１３２からの出力、及び断層パーシスタンス部１３３からの出力のいずれかが含まれる。断層画像データ生成部１３は、ログ圧縮部１３０と、エッジ強調部１３１と、断層平滑化処理部１３２と、断層パーシスタンス部１３３とを有する。ログ圧縮部１３０は、対数変換による圧縮処理を受信信号に施し、エッジ強調部１３１と記憶部１５０とへ出力する。エッジ強調部１３１は、ログ圧縮部１３０からの信号のエッジを強調し、断層平滑化処理部１３２へ出力する。断層平滑化処理部１３２は、エッジ強調部１３１からの信号に平滑化フィルタリング処理を施し、断層パーシスタンス部１３３へ出力する。断層パーシスタンス部１３３は、断層平滑化処理部１３２から過去に受けた信号と断層平滑化処理部１３２から新たに受けた信号とを重み付け加算する。断層パーシスタンス部１３３は、重み付け加算した信号を画像生成部１７へ出力する。
（血流情報生成部１４）
血流情報生成部１４は、受信信号を受信部１２から受け、指定された信号処理条件に基づく処理を受信信号に施して被検体の血流情報を求める。血流情報生成部は、受信信号を受け、被検体の血流情報に係る速度閾値を含む信号処理条件に基づく処理を受信信号に施して血流情報を求める。信号処理条件は、後述する、周波数特性、速度閾値、分散閾値、下限パワー閾値、若しくは上限パワー閾値、又はこれらの組み合わせを含む。血流情報は、被検体における血流の速度、分散、若しくはパワー、又はこれらの組み合わせを含む。また、血流情報生成部１４は、一定の時間間隔ごとに血流情報を求める。ある時間について求められた血流情報は、求められた時間についての静止画（フレーム）に対応する。また、血流情報生成部１４は、条件変更部１６が信号処理条件を変更したとき、変更された信号処理条件に基づく処理を受信信号に施して血流情報を求める。信号処理条件は、ユーザによって指定されてもよく、プリセットされてもよい。血流情報生成部１４は、周波数フィルタ１４０、自己相関部１４１、算出部１４２、ブランク処理部１４３、血流平滑化処理部１４４、及び血流パーシスタンス部１４５を有する。
周波数フィルタ１４０は、指定された周波数特性に基づくフィルタ処理を受信部１２からの受信信号に施す。また、周波数フィルタ１４０は、条件変更部１６が周波数特性を変更したとき、変更された周波数特性に基づくフィルタ処理を受信信号に施す。それにより、周波数フィルタ１４０は、受信信号に含まれるクラッタを低減する。周波数特性には、例えば、バタワース、ベッセル、及びチェビシェフのフィルタタイプ、また、各フィルタタイプにおける遮断周波数、通過帯域、フィルタ次数が含まれる。周波数フィルタ１４０は、フィルタ処理を施した信号を自己相関部１４１へ出力する。
自己相関部１４１は、周波数フィルタ１４０からの信号について自己相関演算を行う。自己相関部１４１は自己相関演算を行った信号を算出部１４２へ出力する。
算出部１４２は、自己相関部１４１から信号を受け、被検体における血流の速度を算出する。算出部１４２は、算出した速度をブランク処理部１４３へ出力する。また、算出部１４２は、被検体における血流の分散を算出する。算出部１４２は、算出した分散をブランク処理部１４３へ出力する。また、算出部１４２は、被検体における血流のパワーを算出する。算出部１４２は、算出したパワーをブランク処理部１４３へ出力する。
ブランク処理部１４３は、算出部１４２により算出された速度について、血流情報として求めるための境目の値である速度閾値を記憶し、算出部１４２から受けた速度のうち速度閾値以上の速度を血流平滑化処理部１４４へ出力する。すなわち、ブランク処理部１４３は、算出部１４２から受けた速度のうち速度閾値未満の速度をクラッタから算出されたものとして取り除き、そして、取り除かれなかった速度を血流平滑化処理部１４４へ出力する。なお、ブランク処理部１４３は、条件変更部１６が速度閾値を変更したとき、変更された速度閾値以上の速度を血流平滑化処理部１４４へ出力する。
また、ブランク処理部１４３は、算出部１４２により算出された分散について、血流情報として求めるための境目の値である分散閾値を記憶し、算出部１４２から受けた分散のうち分散閾値以上の分散を血流平滑化処理部１４４へ出力する。すなわち、ブランク処理部１４３は、算出部１４２から受けた分散のうち分散閾値未満の分散をクラッタから算出されたものとして取り除き、そして、取り除かれなかった分散を血流平滑化処理部１４４へ出力する。なお、ブランク処理部１４３は、条件変更部１６が分散閾値を変更したとき、変更された分散閾値以上の分散を血流平滑化処理部１４４へ出力する。
また、ブランク処理部１４３は、算出部１４２により算出されたパワーについて、血流情報として求めるための境目の値である下限パワー閾値と上限パワー閾値とを記憶し、算出部１４２から受けたパワーのうち下限パワー閾値以上かつ上限パワー以下のパワーを血流平滑化処理部１４４へ出力する。すなわち、ブランク処理部１４３は、算出部１４２から受けたパワーのうち下限パワー閾値未満のパワーと上限パワー閾値超のパワーとをクラッタから算出されたものとして取り除き、そして、取り除かれなかったパワーを血流平滑化処理部１４４へ出力する。なお、ブランク処理部１４３は、条件変更部１６が下限パワー閾値と上限パワー閾値とを変更したとき、変更された下限パワー閾値以上かつ上限パワー以下のパワーを血流平滑化処理部１４４へ出力する。
血流平滑化処理部１４４は、ブランク処理部１４３による出力を受け平滑化フィルタリング処理を施し、血流パーシスタンス部１４５へ出力する。なお、この実施形態では、算出部１４２による出力をブランク処理部１４３が受け、ブランク処理部１４３による出力を血流平滑化処理部１４４が受ける構成を説明するが、算出部１４２による出力を血流平滑化処理部１４４が受け、ブランク処理部１４３による出力を血流平滑化処理部１４４が受ける構成でもよい。
血流パーシスタンス部１４５は、血流平滑化処理部１４４から過去に受けた信号と血流平滑化処理部１４４から新たに受けた信号とを重み付け加算する。血流パーシスタンス部１４５は、重み付け加算した信号を画像生成部１７へ出力する。
（変化検出部１５）
変化検出部１５は、時間が互いに異なる複数の断層画像データを断層画像データ生成部１３から受け、受けた複数の断層画像データの時間軸方向における変化の大きさを検出する。変化検出部１５は、最新の時間について求められた断層画像データである最新フレームと、最新の時間より過去の時間について求められた断層画像データである過去フレームとを記憶部１５０から読み出す。変化検出部１５は、最新フレームにおける信号の振幅の平均値である最新平均値と、過去フレームにおける信号の振幅の平均値である過去平均値とを算出する。変化検出部１５は、最新平均値と過去平均値との差分を算出し、この差分を断層画像データの時間軸方向における変化の大きさとして検出する。このとき、算出された差分が大きいほど、検出される変化の大きさは大きい。変化検出部１５は、検出した変化の大きさを条件変更部１６へ出力する。なお、変化検出部１５は、記憶部１５０が過去フレームを記憶していないとき、断層画像データの時間軸方向における変化の大きさを検出しない。
ここで、断層画像データの時間軸方向における変化と被検体組織の動き及び超音波プローブ１０の動きとの関係について説明する。断層画像データの時間軸方向における変化とは、過去フレームに描出された被検体の組織形状（画像の内容）に対する最新フレームに描出された被検体の組織形状（画像の内容）の差異のことである。例えば、被検体組織が動いてクラッタが発生したとき、動いた組織に相当する部分画像に差異が生じる。また、超音波プローブ１０が動いてモーションアーチファクトが発生したとき、画像に描出される組織が平行移動や回転移動したように描出される。超音波診断装置１が所定のフレームレートで断層画像を取得しているとき、被検体組織又は超音波プローブ１０が動くと、最新フレームに描出される被検体の組織形状と過去フレームに描出される被検体の組織形状（画像の内容）とに差異が生じる。また、被検体組織又は超音波プローブ１０が大きく動くほど、断層画像に描出される被検体の組織形状の変化が大きい。変化検出部１５はこの変化の大きさを検出し、条件変更部１６へ出力する。
また、変化検出部１５は、最新フレームと過去フレームとを複数の領域に分割して最新平均値と過去平均値とを算出し、分割された領域毎に最新平均値と過去平均値との差分を算出することによって断層画像データの時間軸方向における変化の大きさを検出してもよい。
また、変化検出部１５は、最新フレームと過去フレームとに指定された部分領域（関心領域）について、最新フレームの信号の平均値と過去フレームの信号の平均値とを算出し、これら平均値の差分を算出して断層画像データの時間軸方向における変化の大きさを検出してもよい。このとき、部分領域はユーザによって指定されてもよく、プリセットされてもよい。
また、変化検出部１５は、相互相関解析によって最新フレームと過去フレームとの類似度を求め、この類似度を断層画像データの時間軸方向における変化の大きさとして検出してもよい。このとき、求められた類似度が小さいほど、検出される変化の大きさは大きい。
なお、変化検出部１５は、最新フレームと、時間が互いに異なる複数の過去フレームとを記憶部１５０から読み出してもよい。このとき、変化検出部１５は、３以上のフレームに基づいて断層画像データの時間軸方向における変化の大きさを検出する。例えば、変化検出部１５は、２つの時間の過去フレームから過去の断層画像データの時間軸方向における変化の大きさを求め、また、最新フレームと最新フレームの時間に最も近い時間の過去フレームとから最新の断層画像データの時間軸方向における変化の大きさを求める。そして変化検出部１５は、過去の断層画像データの時間軸方向における変化の大きさに重み付けをして最新の断層画像データの時間軸方向における変化の大きさに加算又は乗算することによって断層画像データの時間軸方向における変化の大きさを検出する。また、例えば、変化検出部１５は、複数の過去フレームの信号を外挿し、最新のフレームと同じ時間についてのフレームを予測する。この予測されたフレームを予測フレームとする。変化検出部１５は、予測フレームと記憶部１５０から読み出した最新フレームとの差分を算出する。変化検出部１５は、この差分に重み付けをして最新の断層画像データの時間軸方向における変化の大きさに加算又は乗算することによって断層画像データの時間軸方向における変化の大きさを検出してもよい。
また、変化検出部１５は、記憶部１５０を有する。記憶部１５０は、ログ圧縮部１３０から受けた断層画像データを記憶する。ここで、記憶部１５０は、変化検出部１５が断層画像データの時間軸方向における変化の大きさの検出に用いるフレーム数の断層画像データを少なくとも記憶する。例えば、変化検出部１５が最新フレームと１つの過去フレームに基づいて断層画像データの時間軸方向における変化の大きさを検出するとき、記憶部１５０は最新フレームと１つの過去フレームとの２つの断層画像データを記憶する。記憶部１５０は、変化検出部１５が断層画像データの時間軸方向における変化の大きさを検出したとき、記憶していた１つの過去フレームを削除し、記憶していた最新フレームを過去フレームとして記憶し、新たな断層画像データを最新フレームとして記憶する。同様に、変化検出部１５がｎ枚のフレームに基づいて断層画像データの時間軸方向における変化の大きさを検出するとき、記憶部１５０は、最新フレームとｎ−１枚の過去フレームを記憶する。記憶部１５０は、変化検出部１５が断層画像データの時間軸方向における変化の大きさを検出したとき、ｎ−１枚の過去フレームのうち最も過去の時間の過去フレームを削除し、新たな断層画像データを最新フレームとして記憶する。また、記憶部１５０は、超音波診断中のフリーズ動作、超音波の送信条件の変更動作が行われたとき、当該する制御信号をシステム制御部２０から受け、記憶していた断層画像データをすべて削除してもよい。なお、記憶部１５０は、ログ圧縮部１３０ではなく、受信部１２、エッジ強調部１３１、断層平滑化処理部１３２、又は断層パーシスタンス部１３３による出力を受けて記憶してもよい。
（条件変更部１６）
条件変更部１６は、変化検出部１５が検出した断層画像データの時間軸方向における変化の大きさに基づいて信号処理条件を変更する。例えば、条件変更部１６は、信号処理条件のうち周波数フィルタ１４０の周波数特性を変更する。条件変更部１６は、断層画像データの時間軸方向における変化の大きさと周波数特性とを対応付けて予め記憶する。条件変更部１６は、変化検出部１５による断層画像データの時間軸方向における変化の大きさに対応付けられた周波数特性を読み出し、周波数フィルタ１４０の周波数特性を読み出した周波数特性に変更する。
また、条件変更部１６は、例えば、信号処理条件のうち速度閾値を変更する。このとき、条件変更部１６は、変化検出部１５が検出した断層画像データの時間軸方向における変化が大きいほど、速度閾値を高くする。また、条件変更部１６は、例えば、信号処理条件のうち分散閾値を変更する。このとき、条件変更部１６は、変化検出部１５が検出した断層画像データの時間軸方向における変化が大きいほど、分散閾値を低くする。また、条件変更部１６は、例えば、信号処理条件のうち下限パワー閾値と上限パワー閾値とを変更する。このとき、条件変更部１６は、変化検出部１５が検出した断層画像データの時間軸方向における変化が大きいほど、下限パワー閾値を高くし、上限パワー閾値を低くする。例えば、条件変更部１６は、断層画像データの時間軸方向における変化の大きさと、速度閾値、分散閾値、下限パワー閾値、及び上限パワー閾値とを対応付けて記憶する。条件変更部１６は、変化検出部１５から受けた断層画像データの時間軸方向における変化の大きさに対応付けられた速度閾値、分散閾値、下限パワー閾値、及び上限パワー閾値を読み出し、ブランク処理部１４３へ出力することによって、ブランク処理部１４３における速度閾値、分散閾値、下限パワー閾値、及び上限パワー閾値を変更する。図３は、条件変更部１６が対応付けて記憶した断層画像データの時間軸方向における変化の大きさと、速度閾値Ｖ、分散閾値Ｔ、下限パワー閾値Ｐ１、及び上限パワー閾値Ｐ２との例を示す模式図である。図３において、断層画像データの時間軸方向における変化の大きさは大きさＡ、大きさＢ、大きさＣ、大きさＤ、大きさＥの順に大きいとする。大きさＣと、速度閾値Ｖの値である「α」、分散閾値Ｔの値である「β」、下限パワー閾値Ｐ１の値である「γ」、及び上限パワー閾値Ｐ２の値である「δ」が対応付けられている。断層画像データの時間軸方向における変化の他の大きさ（大きさＡ、大きさＢ、大きさＤ、大きさＥ）には、図３にて示される係数が「α」、「β」、「γ」、及び「δ」のそれぞれに乗算した値が対応付けられる。「α」、「β」、「γ」、及び「δ」それぞれの値は、例えば、「α＝０．１」、「β＝０．９」、「γ＝０．８」、「δ＝０．１」などの値がユーザによって指定されてもよく、自動的にプリセットされてもよい。
（画像生成部１７）
画像生成部１７は、断層画像データ生成部１３により出力された断層画像データと血流情報生成部１４により出力された血流情報とに基づいて超音波画像データを生成する。画像生成部１７は、例えばＤＳＣ（Ｄｉｇｉｔａｌ Ｓｃａｎ Ｃｏｎｖｅｒｔｅｒ：デジタルスキャンコンバータ）を有する。画像生成部１７は、走査線の信号列で表される断層画像データと血流情報とを直交座標系で表される画像データに変換する（スキャンコンバージョン処理）。画像生成部１７が生成した超音波画像データは、断層画像データ（Ｂモード画像データ）を表す断層画像（Ｂモード画像）と血流情報（カラードプラ情報）を表す血流画像（カラードプラ画像）とが重ねられた画像を表すデータである。画像生成部１７は、生成した超音波画像データを表示制御部１８へ出力する。
（表示制御部１８）
表示制御部１８は、超音波画像データを画像生成部１７から受け、超音波画像データに基づく超音波画像を表示部１９に表示させる。
表示部１９は、超音波画像を表示する。表示部１９は、例えば、ＣＲＴ（Ｃａｔｈｏｄｅ Ｒａｙ Ｔｕｂｅ）やＬＣＤ（Ｌｉｑｕｉｄ Ｃｒｙｓｔａｌ Ｄｉｓｐｌａｙ）などの表示デバイスで構成される。表示部１９は、必ずしも超音波診断装置１の一体として備えられる必要はなく、一般的なインターフェイスを介して表示制御部１８によって制御され、超音波画像を表示する構成でもよい。
（システム制御部２０）
システム制御部２０は、超音波診断装置１の各部を制御する。システム制御部２０は、例えば、記憶装置と処理装置とを含んで構成される。記憶装置には、超音波診断装置１の各部の機能を実行するためのコンピュータプログラムが記憶されている。処理装置は、これらコンピュータプログラムを実行することで、上記機能を実現する。
（操作部２１）
操作部２１は、ユーザによる操作を受けて、この操作の内容に応じた信号や情報を装置各部に入力する。操作部２１は、例えば、キーボード、マウス、タッチパネルなどによって構成される。また、操作部２１は、必ずしも超音波診断装置１の一体として備えられる必要はなく、一般的なインターフェイスを介して信号や情報を装置各部に入力する構成でもよい。
［動作］
図４は、この実施形態の超音波診断装置１の動作を示すフローチャートである。
（Ｓ０１）
断層画像データ生成部１３は、超音波を送信された被検体からの反射波に基づく受信信号を受信部１２から受け、被検体の断層画像データを求める。このとき、ログ圧縮部１３０は、対数変換による圧縮処理を受信信号に施し、エッジ強調部１３１と記憶部１５０とへ出力する。エッジ強調部１３１は、ログ圧縮部１３０からの信号のエッジを強調し、断層平滑化処理部１３２へ出力する。断層平滑化処理部１３２は、エッジ強調部１３１からの信号に平滑化フィルタリング処理を施し、断層パーシスタンス部１３３へ出力する。断層パーシスタンス部１３３は、断層平滑化処理部１３２から過去に受けた信号と断層平滑化処理部１３２から新たに受けた信号とを重み付け加算する。断層パーシスタンス部１３３は、重み付け加算した信号を画像生成部１７へ出力する。
（Ｓ０２）
記憶部１５０は、ログ圧縮部１３０から受けた断層画像データを記憶する。
（Ｓ０３）
記憶部１５０が、ログ圧縮部１３０から受けた最新の断層画像データ（最新フレーム）とともに該断層画像データより過去についての断層画像データ（過去フレーム）を記憶しているとき（ステップＳ０３；ＹＥＳ）、ステップＳ０４へ進む。記憶部１５０が、ログ圧縮部１３０から受けた最新の断層画像データ（最新フレーム）とともに該断層画像データより過去についての断層画像データ（過去フレーム）を記憶していないとき（ステップＳ０３；ＮＯ）、ステップＳ０７へ進む。
（Ｓ０４）
変化検出部１５は、最新の時間について求められた断層画像データである最新フレームと、最新の時間より過去の時間について求められた断層画像データである過去フレームとを記憶部１５０から読み出す。変化検出部１５は、読み出した断層画像データに基づいて断層画像データの時間軸方向における変化の大きさを検出する。変化検出部１５は、検出した変化の大きさを条件変更部１６へ出力する。
（Ｓ０５）
記憶部１５０は、記憶していた過去フレームを削除し、記憶していた最新フレームを過去フレームとして記憶する。
（Ｓ０６）
条件変更部１６は、変化検出部１５による断層画像データの時間軸方向における変化の大きさに基づいて血流情報生成部１４の信号処理条件を変更する。信号処理条件には、周波数フィルタ１４０の周波数特性、ブランク処理部１４３の速度閾値Ｖ、分散閾値Ｔ、若しくは下限パワー閾値Ｐ１及び上限パワー閾値Ｐ２、又はこれらの組み合わせが含まれる。このとき、条件変更部１６は、変化検出部１５が検出した断層画像データの時間軸方向における変化が大きいほど、速度閾値Ｖを高くし、分散閾値Ｔを低くし、下限パワー閾値Ｐ１を高くし、上限パワー閾値Ｐ２を低くする。
（Ｓ０７）
血流情報生成部１４は、受信信号を受信部１２から受け、信号処理条件に基づく処理を受信信号に施して被検体の血流情報を求める。このとき、周波数フィルタ１４０は、周波数特性に基づくフィルタ処理を受信部１２からの受信信号に施す。自己相関部１４１は、周波数フィルタ１４０からの信号について自己相関演算を行う。算出部１４２は、被検体における血流の速度、分散、若しくはパワー、又はこれらの組み合わせを算出する。ブランク処理部１４３は、算出部１４２から受けた速度のうち速度閾値以上の速度を血流平滑化処理部１４４へ出力する。また、ブランク処理部１４３は、算出部１４２から受けた分散のうち分散閾値以上の分散を血流平滑化処理部１４４へ出力する。また、ブランク処理部１４３は、算出部１４２から受けたパワーのうち下限パワー閾値以上かつ上限パワー以下のパワーを血流平滑化処理部１４４へ出力する。血流平滑化処理部１４４は、ブランク処理部１４３による出力を受け平滑化フィルタリング処理を施し、血流パーシスタンス部１４５へ出力する。血流パーシスタンス部１４５は、血流平滑化処理部１４４から過去に受けた信号と血流平滑化処理部１４４から新たに受けた信号とを重み付け加算する。血流パーシスタンス部１４５は、重み付け加算した信号を画像生成部１７へ出力する。
（Ｓ０８）
画像生成部１７は、断層画像データ生成部１３により出力された断層画像データと血流情報生成部１４により出力された血流情報とに基づいて超音波画像データを生成する。
（Ｓ０９）
表示制御部１８は、超音波画像データを画像生成部１７から受け、超音波画像データに基づく超音波画像を表示部１９に表示させる。
（Ｓ１０）
超音波診断を続行するとき（ステップＳ１０；ＹＥＳ）、ステップＳ０１へ戻る。超音波診断を続行しないとき（ステップＳ１０；ＮＯ）、動作を終了する。
［効果］
この実施形態の超音波診断装置１の効果について説明する。
実施形態の超音波診断装置１により表示された超音波画像を図５Ａに、この実施形態の変化検出部１５と条件変更部１６との機能を停止した超音波診断装置による超音波画像を図５Ｂに示す。説明のため、図５Ａ及び図５Ｂともに、白黒反転させている。図５Ａと図５Ｂとにおいて、断層画像ＢＲと血流画像ＣＤとが重ねられている。図５Ｂにおいて、破線ＢＫで囲まれた領域に描出された像はクラッタが描出された像ＣＬである。図５Ａでは、血流画像ＣＤ内にクラッタが描出された像ＣＬが描出されず、画像の視認性が向上している。
実施形態の超音波診断装置１は、断層画像データ生成部１３と、血流情報生成部１４と、変化検出部１５と、条件変更部１６とを有する。断層画像データ生成部１３は、超音波を送信された被検体からの反射波に基づく受信信号を受け、被検体の断層画像データを複数の時間について逐次求める。血流情報生成部１４は、受信信号を受け、被検体の血流情報に係る速度閾値を含む信号処理条件に基づく処理を受信信号に施して血流情報を求める。変化検出部１５は、時間が互いに異なる複数の断層画像データを受け、受けた複数の断層画像データの時間軸方向における変化の大きさを検出する。条件変更部１６は、変化検出部１５が検出した時間軸方向における変化の大きさに基づいて信号処理条件を変更する。血流情報生成部１４は、条件変更部１６が信号処理条件を変更したとき、変更された信号処理条件に基づく処理を受信信号に施して血流情報を求める。このように、断層画像データの時系列的な変化を検出し、その変化の大きさに基づいて血流情報の信号処理条件を変更する。それにより、血流情報生成部１４は、被検体組織又は超音波プローブ１０の動きの大きさに応じた信号処理条件に基づいてクラッタやモーションアーチファクトを低減し、血流情報を求めることができる。従って、クラッタやモーションアーチファクトを低減することができる超音波診断装置を提供することができる。
〈第１の変形例〉
［構成］
図６は、この変形例の超音波診断装置１の構成を示すブロック図である。この変形例の超音波診断装置１は、前述の実施形態の超音波診断装置１に比べ、主に血流情報生成部１４の構成が異なる。なお、この変形例において、前述の実施形態と同じ構成については同一符号を付してその説明を省略し、異なる構成について主に説明する。
血流情報生成部１４は、血流情報生成部１４は、周波数フィルタ１４０、自己相関部１４１、算出部１４２、ブランク処理部１４３、血流平滑化処理部１４４、及び血流パーシスタンス部１４５に加え、クラッタ推定部１４６を有する。
クラッタ推定部１４６は、血流情報に含まれるクラッタを推定する。換言すると、クラッタ推定部１４６は、算出部１４２が算出した血流の速度、分散、及びパワーに基づいて、周波数フィルタ１４０を通過した信号にどの程度クラッタが含まれているかを推定する。例えば、クラッタ検出部は、血流の速度、分散、及びパワーそれぞれについて基準値を信号処理条件として記憶し、それら基準値のそれぞれと、算出部１４２から受けた血流の速度、分散、及びパワーについての算出値とを比較し、基準値と算出部１４２との差分が大きいほどクラッタが多いと推定する。なお、クラッタ推定部１４６は、基準値と算出部１４２との差分とクラッタの多さとの相関を示す表データ又は函数を記憶し、表データを参照又は函数を用いた演算を行うことにより、クラッタの多さを推定する。例えば、クラッタ推定部１４６は、推定されたクラッタが多いほど、通過帯域が狭くなるように周波数フィルタ１４０のフィルタ特性を変更する。
条件変更部１６は、変化検出部１５が検出した断層画像データの時間軸方向における変化の大きさを受け、受けた断層画像データの時間軸方向における変化の大きさに基づいてクラッタ推定部１４６の信号処理条件である基準値を変更する。条件変更部１６は、断層画像データの時間軸方向における変化が大きいほど、クラッタ推定部１４６がクラッタを多く推定するようにクラッタ推定部１４６の基準値を変更する。
［動作］
図７は、この変形例の超音波診断装置１の動作を示すフローチャートである。
（Ｓ２１）
断層画像データ生成部１３は、超音波を送信された被検体からの反射波に基づく受信信号を受信部１２から受け、被検体の断層画像データを求める。このとき、ログ圧縮部１３０は、対数変換による圧縮処理を受信信号に施し、エッジ強調部１３１と記憶部１５０とへ出力する。エッジ強調部１３１は、ログ圧縮部１３０からの信号のエッジを強調し、断層平滑化処理部１３２へ出力する。断層平滑化処理部１３２は、エッジ強調部１３１からの信号に平滑化フィルタリング処理を施し、断層パーシスタンス部１３３へ出力する。断層パーシスタンス部１３３は、断層平滑化処理部１３２から過去に受けた信号と断層平滑化処理部１３２から新たに受けた信号とを重み付け加算する。断層パーシスタンス部１３３は、重み付け加算した信号を画像生成部１７へ出力する。
（Ｓ２２）
記憶部１５０は、ログ圧縮部１３０から受けた断層画像データを記憶する。
（Ｓ２３）
記憶部１５０が、ログ圧縮部１３０から受けた最新の断層画像データ（最新フレーム）とともに該断層画像データより過去についての断層画像データ（過去フレーム）を記憶しているとき（ステップＳ２３；ＹＥＳ）、ステップＳ２４へ進む。記憶部１５０が、ログ圧縮部１３０から受けた最新の断層画像データ（最新フレーム）とともに該断層画像データより過去についての断層画像データ（過去フレーム）を記憶していないとき（ステップＳ２３；ＮＯ）、ステップＳ２７へ進む。
（Ｓ２４）
変化検出部１５は、最新の時間について求められた断層画像データである最新フレームと、最新の時間より過去の時間について求められた断層画像データである過去フレームとを記憶部１５０から読み出す。変化検出部１５は、読み出した断層画像データに基づいて断層画像データの時間軸方向における変化の大きさを検出する。変化検出部１５は、検出した断層画像データの時間軸方向における変化の大きさを条件変更部１６へ出力する。
（Ｓ２５）
記憶部１５０は、記憶していた過去フレームを削除し、記憶していた最新フレームを過去フレームとして記憶する。
（Ｓ２６）
条件変更部１６は、変化検出部１５による断層画像データの時間軸方向における変化の大きさに基づいて血流情報生成部１４の信号処理条件を変更する。信号処理条件には、周波数フィルタ１４０の周波数特性、ブランク処理部１４３の速度閾値、分散閾値、若しくは下限パワー閾値及び上限パワー閾値、又はこれらの組み合わせが含まれる。このとき、条件変更部１６は、変化検出部１５が検出した変化の大きさが大きいほど、速度閾値を高くし、分散閾値を低くし、下限パワー閾値を高くし、上限パワー閾値を低くする。また、条件変更部１６は、変化検出部１５が検出した断層画像データの時間軸方向における変化の大きさを受け、受けた断層画像データの時間軸方向における変化の大きさの大きさに基づいてクラッタ推定部１４６の信号処理条件である基準値を変更する。
（Ｓ２７）
血流情報生成部１４は、受信信号を受信部１２から受け、信号処理条件に基づく処理を受信信号に施して被検体の血流情報を求める。このとき、周波数フィルタ１４０は、周波数特性に基づくフィルタ処理を受信部１２からの受信信号に施す。自己相関部１４１は、周波数フィルタ１４０からの信号について自己相関演算を行う。算出部１４２は、被検体における血流の速度、分散、若しくはパワー、又はこれらの組み合わせを算出する。ブランク処理部１４３は、算出部１４２から受けた速度のうち速度閾値以上の速度を血流平滑化処理部１４４へ出力する。また、ブランク処理部１４３は、算出部１４２から受けた分散のうち分散閾値以上の分散を血流平滑化処理部１４４へ出力する。また、ブランク処理部１４３は、算出部１４２から受けたパワーのうち下限パワー閾値以上かつ上限パワー以下のパワーを血流平滑化処理部１４４へ出力する。血流平滑化処理部１４４は、ブランク処理部１４３による出力を受け平滑化フィルタリング処理を施し、血流パーシスタンス部１４５へ出力する。血流パーシスタンス部１４５は、血流平滑化処理部１４４から過去に受けた信号と血流平滑化処理部１４４から新たに受けた信号とを重み付け加算する。血流パーシスタンス部１４５は、重み付け加算した信号を画像生成部１７へ出力する。
（Ｓ２８）
クラッタ推定部１４６は、算出部１４２が算出した血流の速度、分散、及びパワーに基づいて、周波数フィルタ１４０を通過した信号にどの程度クラッタが含まれているかを推定し、周波数フィルタ１４０のフィルタ特性を変更する。
（Ｓ２９）
画像生成部１７は、断層画像データ生成部１３により出力された断層画像データと血流情報生成部１４により出力された血流情報とに基づいて超音波画像データを生成する。
（Ｓ３０）
表示制御部１８は、超音波画像データを画像生成部１７から受け、超音波画像データに基づく超音波画像を表示部１９に表示させる。
（Ｓ３１）
超音波診断を続行するとき（ステップＳ３１；ＹＥＳ）、ステップＳ０１へ戻る。超音波診断を続行しないとき（ステップＳ３１；ＮＯ）、動作を終了する。
［効果］
この変形例の超音波診断装置１の効果について説明する。
この変形例の超音波診断装置１は、断層画像データ生成部１３と、血流情報生成部１４と、変化検出部１５と、条件変更部１６とを有する。断層画像データ生成部１３は、超音波を送信された被検体からの反射波に基づく受信信号を受け、被検体の断層画像データを複数の時間について逐次求める。血流情報生成部１４は、受信信号を受け、被検体の血流情報に係る速度閾値を含む信号処理条件に基づく処理を受信信号に施して血流情報を求める。変化検出部１５は、時間が互いに異なる複数の断層画像データを受け、受けた複数の断層画像データの時間軸方向における変化の大きさを検出する。条件変更部１６は、変化検出部１５が検出した時間軸方向における変化の大きさに基づいて信号処理条件を変更する。血流情報生成部１４は、条件変更部１６が信号処理条件を変更したとき、変更された信号処理条件に基づく処理を受信信号に施して血流情報を求める。さらに、血流情報生成部１４は、クラッタ推定部１４６を有する。クラッタ推定部１４６は、算出された血流の速度、分散、及びパワーに基づいて、周波数フィルタ１４０を通過した信号にどの程度クラッタが含まれているかを推定する。条件変更部１６は、変化検出部１５が検出した断層画像データの時間軸方向における変化の大きさを受け、受けた断層画像データの時間軸方向における変化の大きさに基づいてクラッタ推定部１４６の信号処理条件である基準値を変更する。このように、超音波診断装置１は、周波数フィルタ１４０を通過した信号に含まれるクラッタを推定する基準値を断層画像データの時間軸方向における変化の大きさに対応して推定し、推定されたクラッタに基づいて周波数フィルタ１４０のフィルタ特性を変更する。それにより、被検体組織の動きの大きさによって異なるクラッタの多さを逐次推定し、周波数フィルタ１４０にフィードバックすることができる。従って、クラッタやモーションアーチファクトを低減することができる超音波診断装置を提供することができる。
〈第２の変形例〉
この変形例の超音波診断装置１は、クラッタやモーションアーチファクトに基づく像を、血流情報を表す像よりも弱く表示可能な超音波画像を生成する。この変形例の超音波診断装置１は、前述の実施形態の超音波診断装置１に比べ、主に、血流情報生成部１４の構成が異なる。なお、この変形例において、前述の実施形態と同じ構成については同一符号を付してその説明を省略し、異なる構成について主に説明する。
図８は、この変形例の血流情報生成部１４の構成を表すブロック図である。この血流情報生成部１４は、信号処理条件に基づいて、受信信号が血流情報を表すか否か判別する。この血流情報生成部１４は、ブランク処理部１４３に換えて、判別処理部１４７を有する。判別処理部１４７は、算出部１４２により算出された血流情報が被検体の血流動態を表す情報か否かを判別する。
判別処理部１４７は、速度閾値を記憶する。判別処理部１４７は、算出部１４２から受けた速度と記憶した速度閾値とを参照して、該速度が速度閾値以上の速度か、または速度閾値未満の速度かの判別を行う。速度閾値以上と判別された速度は、被検体の血流動態を表す情報に相当する。判別処理部１４７は、判別結果を表す付帯情報を当該速度の速度情報に付帯して、血流平滑化処理部１４４へ出力する。このとき、判別処理部１４７は、速度閾値未満の速度を取り除かず、速度閾値未満である判定結果を表す付帯情報を付帯して、当該速度の速度情報を血流平滑化処理部１４４へ出力する。なお、判別処理部１４７は、条件変更部１６が速度閾値を変更したとき、変更された速度閾値を参照して、この判別を行う。
また、判別処理部１４７は、分散閾値を記憶する。判別処理部１４７は、算出部１４２から受けた分散と記憶した分散閾値とを参照して、該分散が分散閾値以上の分散か、または分散閾値未満の分散かの判別を行う。分散閾値以上と判別された分散は、被検体の血流動態を表す情報に相当する。判別処理部１４７は、判別結果を表す付帯情報を当該分散に付帯して、当該分散を血流平滑化処理部１４４へ出力する。このとき、判別処理部１４７は、分散閾値未満の分散を取り除かず、分散閾値未満である判定結果を表す付帯情報を付帯して、当該分散を血流平滑化処理部１４４へ出力する。なお、判別処理部１４７は、条件変更部１６が分散閾値を変更したとき、変更された速度閾値を参照して、この判別を行う。
また、判別処理部１４７は、下限パワー閾値及び上限パワー閾値を記憶する。判別処理部１４７は、算出部１４２から受けたパワーと記憶した下限パワー閾値及び上限パワー閾値とを参照して、該パワーが下限パワー閾値以上かつ上限パワー閾値以下であるか否かを判別する。下限パワー閾値以上かつ上限パワー閾値以下と判別されたパワーは、被検体の血流動態を表す情報に相当する。判別処理部１４７は、判別結果を表す付帯情報を当該パワーに付帯して、当該パワーを血流平滑化処理部１４４へ出力する。このとき、判別処理部１４７は、下限パワー閾値以上かつ上限パワー閾値以下でないパワーを取り除かず、下限パワー閾値以上かつ上限パワー閾値以下でない判定結果を表す付帯情報を付帯して、当該パワーを血流平滑化処理部１４４へ出力する。なお、判別処理部１４７は、条件変更部１６が下限パワー閾値及び上限パワー閾値を変更したとき、変更された下限パワー閾値及び上限パワー閾値を参照して、この判別を行う。
判別部がこのような各種判別処理を行うことにより、各種判別結果を表す付帯情報を含む血流情報（カラードプラ情報）が画像生成部１７へ出力される。
変化検出部１５は、検出した変化の大きさを条件変更部１６へ出力するとともに、画像生成部１７へも出力する。画像生成部１７は、判別処理部１４７による判別結果に基づいて、血流情報を表すと判別された受信信号に基づく像よりも、血流情報を表さないと判別された受信信号に基づく像を弱く表示可能な超音波画像データを生成する。画像生成部１７は、断層画像データ生成部１３による断層画像データと血流情報生成部１４による血流情報と変化検出部１５による変化の大きさとに基づいて超音波画像データを生成する。画像生成部１７は、血流情報のうち、被検体の血流動態を表すと判別された血流情報と被検体の血流動態を表さないと判別された血流情報とで描出態様が異なる超音波画像データを生成する。
例えば、画像生成部１７は、血流情報の付帯情報を参照して判別処理部１４７による判別結果を特定する。画像生成部１７は、被検体の血流動態を表すと判別された血流情報を表す画素の輝度より、被検体の血流動態を表さないと判別された血流情報を表す画素の輝度が小さくなるように超音波画像データを生成する。
このとき、画像生成部１７は、変化検出部１５による変化の大きさが大きいほど、被検体の血流動態を表さないと判別された血流情報を表す画素の輝度を小さくする。輝度を小さくすることは、クラッタやモーションアーチファクトに基づく像の表示出力を弱くすることに相当する。例えば、画像生成部１７は、変化検出部１５による変化の大きさに応じて、被検体の血流動態を表さないと判別された血流情報を表す画素の輝度を段階的に小さくする。なお、画像生成部１７は、変化検出部１５による変化の大きさに応じて、被検体の血流動態を表さないと判別された血流情報を表す画素の輝度を無段階に小さくしてもよい。
画像生成部１７は、このように生成した超音波画像データを表示制御部１８へ出力する。表示制御部１８は、この超音波画像データに基づく超音波画像を表示部１９に表示させる。それにより、被検体の血流動態を表さないと判別された血流情報を表す画素は、変化検出部１５による変化の大きさが大きいほど、被検体の血流動態を表すと判別された血流情報を表す画素よりも暗い色で表示される。したがって、ユーザは、超音波画像を視認しながら、暗い色で表示された画素の部分をクラッタやモーションアーチファクトにより描出された部分であると確認することができる。
なお、画像生成部１７は、被検体の血流動態を表さないと判別された血流情報を表す画素が透かし表示されるように、超音波画像データを生成してもよい。例えば、画像生成部１７は、被検体の血流動態を表さないと判別された血流情報を表す画素と、該画素と同じ座標である断層画像データの画素とのＲＧＢ（Ｒｅｄ−Ｇｒｅｅｎ−Ｂｌｕｅ）信号を加算する。
このとき、画像生成部１７は、変化検出部１５による変化の大きさが大きいほど、被検体の血流動態を表さないと判別された血流情報を表す画素のＲＧＢ信号の加算比率を小さくし、該画素と同じ座標である断層画像データの画素のＲＧＢ信号の加算比率を大きくする。被検体の血流動態を表さないと判別された血流情報を表す画素のＲＧＢ信号の加算比率を小さくし、該画素と同じ座標である断層画像データの画素のＲＧＢ信号の加算比率を大きくすることは、クラッタやモーションアーチファクトに基づく像の表示出力を弱くすることに相当する。例えば、画像生成部１７は、変化検出部１５による変化の大きさに応じて、被検体の血流動態を表さないと判別された血流情報を表す画素のＲＧＢ信号の加算比率を段階的に小さくし、該画素と同じ座標である断層画像データの画素のＲＧＢ信号の加算比率を段階的に大きくする。なお、画像生成部１７は、変化検出部１５による変化の大きさに応じて、被検体の血流動態を表さないと判別された血流情報を表す画素のＲＧＢ信号の加算比率を無段階に小さくし、該画素と同じ座標である断層画像データの画素のＲＧＢ信号の加算比率を無段階に大きくしてもよい。なお、画像生成部１７は、被検体の血流動態を表すと判別された血流情報を表す画素については、この加算を行わない。画像生成部１７は、被検体の血流動態を表すと判別された血流情報を表す画素について、断層画像にからドプラ画像が重ねられるように超音波画像データを生成する。
このように生成された超音波画像データに基づく超音波画像が表示されることにより、被検体の血流動態を表さないと判別された血流情報を表す画素部分は、変化検出部１５による変化の大きさが大きいほど、断層画像の画素のＲＧＢ信号が強い比率で表示され、カラードプラ画像のＲＧＢ信号が弱い比率で表示され、カラードプラ画像が透かし表示される。被検体の血流動態を表さないと判別された血流情報を表す画素部分は、この透かし表示がされず、断層画像にカラードプラ画像が重ねられて表示さる。したがって、ユーザは、超音波画像を視認しながら、カラードプラ画像が透かし表示された画素の部分をクラッタやモーションアーチファクトにより描出された部分であると確認することができる。
上記のいくつかの実施形態を実現するための超音波画像処理プログラムを、コンピュータによって読み取り可能な任意の記録媒体に記憶させることができる。この記録媒体としては、たとえば、半導体メモリ、光ディスク、光磁気ディスク、磁気記憶媒体などを用いることが可能である。また、インターネットやＬＡＮ等のネットワークを通じてこのプログラムを送受信することも可能である。
この発明の実施形態を説明したが、上記の実施形態は例として提示したものであり、発明の範囲を限定することを意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
１ 超音波診断装置
１０ 超音波プローブ
１１ 送信部
１２ 受信部
１３ 断層画像データ生成部
１４ 血流情報生成部
１５ 変化検出部
１６ 条件変更部
１７ 画像生成部
１８ 表示制御部
１９ 表示部
２０ システム制御部
２１ 操作部
１３０ ログ圧縮部
１３１ エッジ強調部
１３２ 断層平滑化処理部
１３３ 断層パーシスタンス部
１４０ 周波数フィルタ
１４１ 自己相関部
１４２ 算出部
１４３ ブランク処理部
１４４ 血流平滑化処理部
１４５ 血流パーシスタンス部
１４６ クラッタ推定部
１４７ 判別処理部
１５０ 記憶部

Claims (11)

  1. 超音波を送信された被検体からの反射波に基づく受信信号を受け、前記被検体の断層画像データを複数の時間について逐次求める断層画像データ生成部と、
    前記受信信号を受け、前記被検体の血流情報に係る速度閾値を含む信号処理条件に基づく処理を前記受信信号に施して前記血流情報を求める血流情報生成部と、
    時間が互いに異なる複数の前記断層画像データを受け、受けた複数の前記断層画像データの時間軸方向における変化の大きさを検出する変化検出部と、
    前記変化検出部が検出した時間軸方向における変化の大きさに基づいて前記信号処理条件を変更する条件変更部と
    を有し、
    前記血流情報生成部は、前記条件変更部が前記信号処理条件を変更したとき、変更された前記信号処理条件に基づく処理を前記受信信号に施して前記血流情報を求める
    ことを特徴とする超音波診断装置。
  2. 前記血流情報生成部は、前記被検体における血流の速度を算出し、
    前記速度閾値は、算出された前記速度について、前記血流情報生成部が前記血流情報として求めるための境目の値であり、
    前記血流情報生成部は、算出した前記速度のうち、前記速度閾値以上の速度を前記血流情報として求め、
    前記条件変更部は、前記速度閾値を変更し、
    前記血流情報生成部は、前記条件変更部が前記速度閾値を変更したとき、変更された前記速度閾値以上の速度を前記血流情報として求める
    ことを特徴とする請求項１に記載の超音波診断装置。
  3. 前記条件変更部は、前記変化検出部が検出した時間軸方向における変化の大きさが大きいほど、前記速度閾値を高くすることを特徴とする請求項２に記載の超音波診断装置。
  4. 前記血流情報生成部は、前記被検体における血流の分散を算出し、
    前記信号処理条件は、算出された前記分散について、前記血流情報生成部が前記血流情報として求めるための境目の値である分散閾値を含み、
    前記血流情報生成部は、算出した前記分散のうち、前記分散閾値以上の分散を前記血流情報として求め、
    前記条件変更部は、前記分散閾値を変更し、
    前記血流情報生成部は、前記条件変更部が前記分散閾値を変更したとき、変更された前記分散閾値以上の分散を前記血流情報として求める
    ことを特徴とする請求項１に記載の超音波診断装置。
  5. 前記条件変更部は、前記変化検出部が検出した時間軸方向における変化の大きさが大きいほど、前記分散閾値を低くすることを特徴とする請求項４に記載の超音波診断装置。
  6. 前記血流情報生成部は、前記被検体における血流のパワーを算出し、
    前記信号処理条件は、算出された前記パワーについて、前記血流情報生成部が前記血流情報として求めるための境目の値である下限パワー閾値と上限パワー閾値とを含み、
    前記血流情報生成部は、算出した前記パワーのうち、前記下限パワー閾値以上かつ前記上限パワー閾値以下のパワーを前記血流情報として求め、
    前記条件変更部は、前記下限パワー閾値と前記上限パワー閾値とを変更し、
    前記血流情報生成部は、前記条件変更部が前記下限パワー閾値と前記上限パワー閾値とを変更したとき、変更された前記下限パワー閾値以上かつ前記上限パワー閾値以下のパワーを前記血流情報として求める
    ことを特徴とする請求項１に記載の超音波診断装置。
  7. 前記条件変更部は、前記変化検出部が検出した時間軸方向における変化の大きさが大きいほど、前記下限パワー閾値を高くし、前記上限パワー閾値を低くすることを特徴とする請求項６に記載の超音波診断装置。
  8. 前記信号処理条件は、周波数フィルタの周波数特性を含み、
    前記血流情報生成部は、指定された前記周波数特性に基づくフィルタ処理を前記受信信号に施す前記周波数フィルタを有し、
    前記条件変更部は、前記信号処理条件に含まれる前記周波数特性を変更し、
    前記血流情報生成部は、前記条件変更部が前記周波数特性を変更したとき、変更された前記周波数特性に基づく処理を前記受信信号に施して前記血流情報を求める
    ことを特徴とする請求項１に記載の超音波診断装置。
  9. 血流情報生成部は、予め記憶した基準値と前記血流情報とに基づいて、前記血流情報に含まれるクラッタを推定し、推定したクラッタに基づいて前記周波数特性を変更するクラッタ推定部を有し、
    前記条件変更部は、前記変化検出部が検出した時間軸方向における変化の大きさに基づいて、前記基準値を変更する
    ことを特徴とする請求項８に記載の超音波診断装置。
  10. 前記血流情報生成部は、前記信号処理条件に基づいて、前記受信信号が前記血流情報を表すか否か判別する判別処理部を有し、
    前記判別処理部による判別結果に基づいて、前記血流情報を表すと判別された前記受信信号に基づく像よりも、前記血流情報を表さないと判別された前記受信信号に基づく像を弱く表示可能な超音波画像データを生成する画像生成部をさらに有する
    ことを特徴とする請求項１に記載の超音波診断装置。
  11. 超音波を送信された被検体からの反射波に基づく受信信号を処理するコンピュータに、
    前記被検体の断層画像データを複数の時間について逐次求めさせ、
    前記被検体の血流情報に係る速度閾値を含む信号処理条件に基づく処理を前記受信信号に施して前記血流情報を求めさせ、
    時間が互いに異なる複数の前記断層画像データの時間軸方向における変化の大きさを検出させ、
    検出された時間軸方向における変化の大きさに基づいて前記信号処理条件を変更させた後、変更された前記信号処理条件に基づいて前記血流情報を求めさせる
    ことを特徴とする超音波画像処理プログラム。
JP2014048258A 2013-03-12 2014-03-12 超音波診断装置及び超音波画像処理プログラム Active JP6274517B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2014048258A JP6274517B2 (ja) 2013-03-12 2014-03-12 超音波診断装置及び超音波画像処理プログラム

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2013049745 2013-03-12
JP2013049745 2013-03-12
JP2014048258A JP6274517B2 (ja) 2013-03-12 2014-03-12 超音波診断装置及び超音波画像処理プログラム

Publications (2)

Publication Number Publication Date
JP2014198240A JP2014198240A (ja) 2014-10-23
JP6274517B2 true JP6274517B2 (ja) 2018-02-07

Family

ID=51536825

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014048258A Active JP6274517B2 (ja) 2013-03-12 2014-03-12 超音波診断装置及び超音波画像処理プログラム

Country Status (4)

Country Link
US (1) US20150359507A1 (ja)
JP (1) JP6274517B2 (ja)
CN (1) CN105073020B (ja)
WO (1) WO2014142174A1 (ja)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2989986B1 (en) * 2014-09-01 2019-12-18 Samsung Medison Co., Ltd. Ultrasound diagnosis apparatus and method of operating the same
CN106028947B (zh) * 2015-05-07 2019-11-12 深圳迈瑞生物医疗电子股份有限公司 血流成像方法及系统
JP7052313B2 (ja) * 2017-11-20 2022-04-12 コニカミノルタ株式会社 超音波信号処理装置、超音波診断装置、および、超音波信号処理方法
EP3820375A1 (en) * 2018-07-11 2021-05-19 Koninklijke Philips N.V. Ultrasound imaging system with pixel extrapolation image enhancement
US20200121296A1 (en) * 2018-10-22 2020-04-23 EchoNous, Inc. Motion artifact suppression in ultrasound color flow imaging
CN109363722B (zh) * 2018-11-29 2021-08-06 武汉中旗生物医疗电子有限公司 彩色血流成像中运动伪像的抑制方法及设备
US11141138B2 (en) 2019-05-28 2021-10-12 Siemens Medical Solutions Usa, Inc. Kalman filtering for flash artifact suppression in ultrasound imaging
JP7343342B2 (ja) 2019-09-25 2023-09-12 キヤノンメディカルシステムズ株式会社 超音波診断装置、及び画像処理装置
JP7469877B2 (ja) 2019-12-27 2024-04-17 キヤノンメディカルシステムズ株式会社 超音波診断装置、医用画像処理装置、および医用画像処理プログラム
JP2021159276A (ja) * 2020-03-31 2021-10-11 キヤノンメディカルシステムズ株式会社 超音波診断装置
JP7440328B2 (ja) 2020-04-03 2024-02-28 キヤノンメディカルシステムズ株式会社 超音波診断装置、およびプログラム

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5170792A (en) * 1989-11-27 1992-12-15 Acoustic Imaging Technologies Corporation Adaptive tissue velocity compensation for ultrasonic Doppler imaging
JP3029315B2 (ja) * 1990-03-30 2000-04-04 株式会社東芝 超音波血流イメージング装置
DE4134724C2 (de) * 1990-10-24 1995-11-16 Hitachi Medical Corp Einrichtung zur farbigen Strömungsaufzeichnung mit Ultraschall
US5445156A (en) * 1994-11-02 1995-08-29 General Electric Company Method for adaptively filtering doppler signals using a complex time domain filter
US5487389A (en) * 1994-12-29 1996-01-30 Siemens Medical Systems, Inc. Ultrasonic Doppler imager having an adaptive tissue rejection filter with enhanced tissue motion sensitivity
JPH1075955A (ja) * 1996-07-11 1998-03-24 Fujitsu Ltd 超音波診断装置
US6520915B1 (en) * 2000-01-28 2003-02-18 U-Systems, Inc. Ultrasound imaging system with intrinsic doppler capability
JP2003250802A (ja) * 2002-03-04 2003-09-09 Fukuda Denshi Co Ltd 超音波診断装置
US20090209861A1 (en) * 2005-02-17 2009-08-20 Panasonic Corporation Ultrasonic Doppler Blood Flow Meter
JP5593116B2 (ja) * 2009-04-30 2014-09-17 三星メディソン株式会社 適応クラッターフィルタリングを行う超音波システムおよび方法
JP5366678B2 (ja) * 2009-06-25 2013-12-11 株式会社東芝 ３次元超音波診断装置及びプログラム
CN102429684B (zh) * 2010-09-28 2013-10-09 深圳迈瑞生物医疗电子股份有限公司 一种多普勒彩色血流成像方法和装置

Also Published As

Publication number Publication date
US20150359507A1 (en) 2015-12-17
CN105073020B (zh) 2017-06-20
JP2014198240A (ja) 2014-10-23
CN105073020A (zh) 2015-11-18
WO2014142174A1 (ja) 2014-09-18

Similar Documents

Publication Publication Date Title
JP6274517B2 (ja) 超音波診断装置及び超音波画像処理プログラム
KR100908252B1 (ko) 영상 처리 시스템 및 방법
JP5645628B2 (ja) 超音波診断装置
JP6218400B2 (ja) 超音波診断装置及び超音波診断装置の制御プログラム
US10893848B2 (en) Ultrasound diagnosis apparatus and image processing apparatus
WO2011036891A1 (ja) 超音波診断装置
US20240023940A1 (en) Ultrasound diagnostic apparatus
KR20150117685A (ko) 초음파 이미지 디스플레이 장치 및 초음파 이미지를 디스플레이하기 위한 방법
JP4971080B2 (ja) 超音波診断装置
US10667792B2 (en) Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus and ultrasonic diagnostic apparatus control method
JP2008154891A (ja) カラードップラ装置および超音波撮像装置
JP4808373B2 (ja) Ｂモード画像のバンディング抑制に関連する応用のための方法及び装置
US9839412B2 (en) Ultrasonic image display apparatus and control program for controlling the same
US11844651B2 (en) Analyzing apparatus and analyzing method using distribution information
KR100998443B1 (ko) 스캔 변환을 고려하여 초음파 데이터를 처리하는 초음파시스템 및 방법
US20210338207A1 (en) Image analyzing apparatus
US20220330920A1 (en) Ultrasonic diagnostic apparatus and medical image processing apparatus
US11786209B2 (en) Ultrasonic diagnostic apparatus and image processing apparatus
JP2017192516A (ja) 超音波診断装置及びプログラム
JP2009261436A (ja) 超音波診断装置
JP2019097794A (ja) 超音波診断装置、医用画像処理装置及びそのプログラム
US20210228177A1 (en) Ultrasonic diagnostic apparatus, learning apparatus, and image processing method
JP2002034987A (ja) Ｂモード画像生成方法および超音波診断装置
JP6253901B2 (ja) 超音波診断装置
JP7469877B2 (ja) 超音波診断装置、医用画像処理装置、および医用画像処理プログラム

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20160527

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20170201

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20171227

R150 Certificate of patent or registration of utility model

Ref document number: 6274517

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150