JP2011045457A - 超音波診断装置 - Google Patents

超音波診断装置 Download PDF

Info

Publication number
JP2011045457A
JP2011045457A JP2009195004A JP2009195004A JP2011045457A JP 2011045457 A JP2011045457 A JP 2011045457A JP 2009195004 A JP2009195004 A JP 2009195004A JP 2009195004 A JP2009195004 A JP 2009195004A JP 2011045457 A JP2011045457 A JP 2011045457A
Authority
JP
Japan
Prior art keywords
unit
physical quantity
correlation
average
elastic image
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.)
Granted
Application number
JP2009195004A
Other languages
English (en)
Other versions
JP5484826B2 (ja
Inventor
Shunichiro Tanigawa
俊一郎 谷川
Hiroji Miyama
広二 見山
Haruji Funaya
晴二 船矢
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.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
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 GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Priority to JP2009195004A priority Critical patent/JP5484826B2/ja
Priority to CN201010273768.8A priority patent/CN101999907B/zh
Priority to US12/868,926 priority patent/US8684931B2/en
Publication of JP2011045457A publication Critical patent/JP2011045457A/ja
Application granted granted Critical
Publication of JP5484826B2 publication Critical patent/JP5484826B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • A61B8/463Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/485Diagnostic techniques involving measuring strain or elastic properties
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52042Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/5205Means for monitoring or calibrating
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52053Display arrangements
    • G01S7/52057Cathode ray tube displays
    • G01S7/52073Production of cursor lines, markers or indicia by electronic means

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Remote Sensing (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

【課題】弾性画像データに基づいて表示される弾性画像が、生体組織の弾性を正確に反映した弾性画像であるか否かを、従来とは異なる観点又は従来よりも幅広い観点から評価することができる超音波診断装置を提供することである。
【解決手段】生体組織における各部の弾性に関する物理量を算出する物理量算出部51と、前記物理量に基づいて、生体組織の弾性画像データを、超音波の送受信面における弾性画像作成領域について作成する弾性画像データ作成部52と、前記弾性画像作成領域における前記物理量の平均をフレーム毎に算出する物理量平均部8と、予め設定された前記物理量の平均値に対する前記物理量平均部8による算出値の比を算出する比算出部9と、比算出部9の算出結果を表示する表示部7と、を備えることを特徴とする。
【選択図】図2

Description

本発明は、超音波診断装置に関し、特に生体組織の硬さ又は軟らかさを表す弾性画像を表示する超音波診断装置に関する。
通常のBモード画像と、生体組織の硬さ又は軟らかさを表す弾性画像とを合成して表示させる超音波診断装置が、例えば特許文献1などに開示されている。この種の超音波診断装置において、弾性画像は次のようにして作成される。先ず、生体組織に対し、超音波プローブによって体表面からの圧迫とその弛緩を繰り返しながら超音波の送受信を行い、エコー信号を取得する。そして、得られたエコー信号に基づいて、生体組織の弾性に関する物理量を算出し、この物理量を色相情報に変換してカラーの弾性画像を作成する。ちなみに、生体組織の弾性に関する物理量としては、例えば生体組織の変形による変位(以下、単に「変位」と云う)などを算出している。
前記物理量の算出手法の一例についてもう少し説明すると、先ず時間的に異なる同一音線上の二つのエコー信号に、所定のデータ数分の幅を有する相関ウィンドウをそれぞれ設定し、この相関ウィンドウ間で相関演算を行なって前記物理量を算出する。例えば特許文献2では、相関ウィンドウ間で相関演算を行なうことによって、両エコー信号の波形のずれを算出し、この波形のずれを変位とみなしている。
ところで、例えば、圧迫とその弛緩の度合いが足りないなど、生体組織の変形が不十分な場合には、相関演算の算出値が生体組織の弾性の違いに応じた差となって現れないことがある。この場合、弾性画像が生体組織の弾性を正確に反映したものとならない。
一方、圧迫とその弛緩の度合いが過剰である場合には、生体組織に横ずれが生じることがある。このような場合に取得されたエコー信号には横ずれによるノイズが含まれ、相関演算における相関係数が低くなるおそれがある。また、圧迫とその弛緩の度合いが過剰であると、生体組織の変形が大きすぎ、二つのエコー信号に設定される相関ウィンドウのマッチングがとれずに相関係数が低くなるおそれがある。ここで、相関演算における相関係数が低くなると、得られる算出値が、生体組織の弾性を正確に反映した算出値を得ることができなくなる。従って、相関演算の相関係数が低くなるほど、生体組織の弾性を正確に反映した弾性画像を得ることができなくなる。
また、超音波の反射体が少ない領域や送信超音波が減衰によって到達しにくい生体組織の深部などにおいては、エコー信号の強度が不十分となる。このように信号強度が不十分なエコー信号についての相関演算の相関係数は低くなる。また、前記超音波プローブの圧迫とその弛緩の方向が超音波の音線方向と一致していない場合、上述の横ずれが生じるため、このような状態で取得されたエコー信号についての相関演算の相関係数も低くなる。従って、これらの場合にも、生体組織の弾性を正確に反映した弾性画像を得ることができない。
しかし、このように生体組織の弾性を正確に反映した弾性画像が得られていないとしても、その弾性画像を見た操作者が、生体組織の弾性を正確に反映した弾性画像であるか否かを判断することは困難である。
そこで、特許文献3には、操作者が生体組織の弾性を正確に反映した弾性画像であるか否かを容易に判断できるようにするために、相関演算で得られた弾性を表す値のフレーム毎の平均を算出し、その平均値について、「H」(high),「M」(medium),「L」(low)の評価結果を表示する手法が開示されている。
特開2005−118152号公報 特開2008−126079号公報 米国特許第6558324号明細書
しかし、特許文献3においては、生体組織の変形が不十分な状態で弾性画像が作成されたことを検出することを目的としているため、平均値が低くなるほど悪い評価(具体的には「L」)がなされる。一方で、平均値が高ければ「H」の評価となる。しかし、上述のように、生体組織への圧迫やその弛緩が過剰に行なわれた場合や、相関係数が低い相関演算で得られた算出値が含まれている場合であっても、生体組織の弾性を正確に反映した弾性画像にはならない。このため、特許文献3の手法では、生体組織の弾性を正確に反映した弾性画像であるか否かの判断をするには不十分な場合があり、判断の正確性に欠ける場合があった。
本発明が解決しようとする課題は、弾性画像データに基づいて表示される弾性画像が、生体組織の弾性を正確に反映した弾性画像であるか否かを、従来とは異なる観点又は従来よりも幅広い観点から評価することができる超音波診断装置を提供することである。
この発明は、前記課題を解決するためになされたもので、第1の観点の発明は、生体組織に対する超音波の送受信により得られた同一音線上の時間的に異なる二つのエコー信号に相関ウィンドウを設定し、該相関ウィンドウ間で相関演算を行なって生体組織における各部の弾性に関する物理量を算出する物理量算出部と、前記物理量に基づいて、生体組織の弾性画像データを、超音波の送受信面における弾性画像作成領域について作成する弾性画像データ作成部と、前記弾性画像作成領域における前記物理量の平均をフレーム毎に算出する物理量平均部と、該物理量平均部による算出値を、予め設定された前記物理量の平均値と比較する比較部と、該比較部による比較結果を報知する報知部と、を備えることを特徴とする超音波診断装置である。
第2の観点の発明によれば、第1の観点の発明において、前記物理量平均部は、所定の閾値以上の相関係数の相関演算が行なわれた相関ウィンドウについて得られた物理量の平均算出を行なうことを特徴とする超音波診断装置である。
第3の観点の発明は、第1又は2の観点の発明において、前記比較部は、前記比較結果として、予め設定された前記物理量の平均値に対する前記物理量平均部による算出値の比を算出することを特徴とする超音波診断装置である。
第4の観点の発明は、第1〜3のいずれか一の観点の発明において、前記報知部には、前記比較部によってフレーム毎に得られる比較結果が複数フレームにわたって平均化されて報知されることを特徴とする超音波診断装置である。
第5の観点の発明は、生体組織に対する超音波の送受信により得られた同一音線上の時間的に異なる二つのエコー信号に相関ウィンドウを設定し、該相関ウィンドウ間で相関演算を行なって生体組織における各部の弾性に関する物理量を算出する物理量算出部と、前記物理量に基づいて、生体組織の弾性画像データを、超音波の送受信面における弾性画像作成領域について作成する弾性画像データ作成部と、前記相関ウィンドウ間の相関演算における相関係数の前記弾性画像作成領域における平均をフレーム毎に算出する相関係数平均部と、該相関係数平均部による算出結果を報知する報知部と、を備えることを特徴とする超音波診断装置である。
第6の観点の発明は、第5の観点の発明において、前記報知部には、前記相関係数平均部によってフレーム毎に得られる算出結果が複数フレームにわたって平均化されて報知されることを特徴とする超音波診断装置である。
第7の観点の発明は、生体組織に対する超音波の送受信により得られた同一音線上の時間的に異なる二つのエコー信号に基づいて生体組織における各部の弾性に関する物理量を算出する物理量算出部と、前記物理量に基づいて、生体組織の弾性画像データを、超音波の送受信面における弾性画像作成領域について作成する弾性画像データ作成部と、所定の閾値以上の相関係数の相関演算が行なわれた相関ウィンドウについて得られた物理量の前記弾性画像作成領域における平均をフレーム毎に算出する物理量平均部と、予め設定された前記物理量の平均値に対する前記物理量平均部による算出値の比を算出する比算出部と、前記相関ウィンドウ間の相関演算における相関係数の前記弾性画像作成領域における平均をフレーム毎に算出する相関係数平均部と、前記比算出部の算出値と、前記相関係数平均部の算出値とを乗算する乗算部と、該乗算部による乗算結果を報知する報知部と、を備えることを特徴とする超音波診断装置である。
第8の観点の発明は、第7の観点の発明において、前記報知部には、前記乗算部によってフレーム毎に得られる乗算結果が複数フレームにわたって平均化されて報知されることを特徴とする超音波診断装置である。
第9の観点の発明は、第7,8の観点の発明において、前記乗算部は、前記比算出部の算出値と、前記相関係数平均部の算出値との重み付け演算を行なうことを特徴とする超音波診断装置である。
第10の観点の発明は、生体組織に対する超音波の送受信により得られた同一音線上の時間的に異なる二つのエコー信号に基づいて生体組織における各部の弾性に関する物理量を算出する物理量算出部と、前記物理量に基づいて、生体組織の弾性画像データを、超音波の送受信面における弾性画像作成領域について作成する弾性画像データ作成部と、所定の閾値以上の相関係数の相関演算が行なわれた相関ウィンドウについて得られた物理量の前記弾性画像作成領域における平均をフレーム毎に算出する物理量平均部と、予め設定された前記物理量の平均値に対する前記物理量平均部による算出値の比を算出する比算出部と、前記相関ウィンドウ間の相関演算における相関係数の前記弾性画像作成領域における平均をフレーム毎に算出する相関係数平均部と、前記比算出部の算出値と、前記相関係数平均部の算出値とを乗算する乗算部と、前記比算出部による算出結果、前記相関係数平均部による算出結果又は前記乗算部による乗算結果を切り替えて報知可能な報知部と、該報知部における切替の指示を入力する操作部と、を備えることを特徴とする超音波診断装置である。
本発明によれば、前記生体組織における各部の物理量の前記弾性画像作成領域における平均値と、予め設定された前記物理量の平均値とを比較した比較結果が表示されるので、操作者は、生体組織に対する圧迫とその弛緩の度合いが足りなかったり、また過剰であったりしないかどうかを容易に判断することができる。これにより、前記弾性画像データに基づいて表示される弾性画像が、生体組織の弾性を正確に反映した弾性画像か否かを従来よりも幅広い観点から評価することができる。
また、前記物理量平均部により、所定の閾値以上の相関係数の相関演算が行なわれた相関ウィンドウについて得られた物理量の平均算出を行なって得られた平均値は、エコー信号の強度が不十分な部分、生体組織の横ずれが生じている部分など、相関係数が低い部分の変位が除かれて得られた平均値である。従って、このような平均値と予め設定された前記物理量の平均値(理想値)との比較結果は、生体組織に対する圧迫とその弛緩が適切な強さで行なわれているか否かを示すものとなる。以上より、操作者は、前記比較結果の表示から、生体組織に対する圧迫とその弛緩が適切な強さで行なわれているか否かをより正確に把握することができる。
また、前記表示部に、前記比較部によってフレーム毎に得られる比較結果が複数フレームにわたって平均化されて表示されることにより、安定した比較結果を表示することができる。
また、他の発明によれば、前記相関ウィンドウ間の相関演算における相関係数の前記弾性画像作成領域における平均の算出結果が表示されるので、操作者は、前記弾性画像データに基づいて表示される弾性画像が、生体組織の弾性を正確に反映した弾性画像か否かを、従来とは異なる観点から評価することができる。
また、前記表示部に、前記相関係数平均部によってフレーム毎に得られる算出結果が複数フレームにわたって平均化されて表示されることにより、安定した算出結果を表示することができる。
また、他の発明によれば、所定の閾値以上の相関係数の相関演算が行なわれた相関ウィンドウについて得られた物理量の平均値を用いて前記比算出部によって算出された算出値と、前記相関係数平均部によって算出された算出値とを前記乗算部によって乗算し、この乗算結果が表示される。この乗算結果は、生体組織への圧迫とその弛緩の度合いの要素と、相関係数の要素とが加味されたものである。これにより、生体組織の弾性を正確に反映した弾性画像であるか否かを、従来よりも幅広い観点から評価することができる。
また、前記表示部に、前記乗算部によってフレーム毎に得られる乗算結果が複数フレームにわたって平均化されて表示されることにより、安定した乗算結果を表示することができる。
また、本発明によれば、前記比算出部による算出結果、前記相関係数平均部による算出結果及び前記乗算部による乗算結果を切り替えて表示することができるので、生体組織の弾性を正確に反映した弾性画像であるか否かを、従来よりも幅広い観点から評価することができる。
本発明に係る超音波診断装置の第一実施形態の概略構成を示すブロック図である。 図1に示す超音波診断装置における要部の構成を示すブロック図である。 Bモード画像データ及び弾性画像データの作成の説明図である。 図1に示す超音波診断装置における表示部の表示の一例を示す図である。 弾性画像データを作成する際における物理量の算出を説明するための図である。 比算出部で用いられる関数のグラフを示す図である。 表示部の表示の一例を示し、時間の経過とともにクオリティ表示が左から右へ流れるように表示されることを説明するための図である。 表示部の表示の一例を示し、時間の経過とともにクオリティ表示が左から右へ流れるように表示されることを説明するための図である。 表示部の表示の一例を示し、時間の経過とともにクオリティ表示が左から右へ流れるように表示されることを説明するための図である。 表示部の表示の一例を示し、クオリティ表示の他例を示す図である。 本発明に係る超音波診断装置の第二実施形態の概略構成を示すブロック図である。 図11に示す超音波診断装置における要部の構成を示すブロック図である。 本発明に係る超音波診断装置の第三実施形態の概略構成を示すブロック図である。 図13に示す超音波診断装置における要部の構成を示すブロック図である。 本発明に係る超音波診断装置の第四実施形態の概略構成を示すブロック図である。 図15に示す超音波診断装置における要部の構成を示すブロック図である。 クオリティ表示の他例を示す図である。 クオリティ表示の他例が表示された表示部を示す図である。
以下、本発明の実施形態について図面に基づいて詳細に説明する。
(第一実施形態)
先ず、第一実施形態について図1〜図10に基づいて説明する。図1に示す超音波診断装置1は、超音波プローブ2、送受信部3、Bモード画像処理部4、弾性画像処理部5、合成部6、表示部7、物理量平均部8、比算出部9、クオリティ表示作成部10を備え、さらに制御部11及び操作部12を備える。
前記超音波プローブ2は、生体組織に対して超音波を送信しそのエコーを受信する。この超音波プローブ2を生体組織の表面に当接させた状態で圧迫と弛緩を繰り返しながら超音波の送受信を行なって取得されたエコー信号に基づいて、後述のように弾性画像が作成される。
前記送受信部3は、前記超音波プローブ2を所定の走査条件で駆動させて音線毎の超音波の走査を行なう。また、前記超音波プローブ2で受信したエコー信号について、整相加算処理等の信号処理を行なう。前記送受信部3で信号処理されたエコー信号は、前記Bモード画像処理部4及び前記弾性画像処理部5に出力される。また、送受信部3から出力されたエコー信号は、ローデータ(Raw Data)として図示しない記憶部に記憶されてもよい。
ちなみに、前記送受信部3は、Bモード画像を作成するための走査と、弾性画像を作成するための走査とを別に行なう。弾性画像を作成するための走査としては、被検体における弾性画像を作成する領域において、同一音線上に二回の走査を行なう。
前記Bモード画像処理部4は、前記送受信部3から出力されたエコー信号に対し、対数圧縮処理、包絡線検波処理等のBモード処理を行い、Bモード画像データを作成する。ただし、前記Bモード画像処理部4は、前記記憶部にローデータとして記憶されたエコー信号に基づいてBモード画像データを作成してもよい。
前記弾性画像処理部5は、前記送受信部3から出力されたエコー信号に基づいて、弾性画像データを作成する。ただし、前記弾性画像処理部5は、前記Bモード画像処理部4と同様に、前記記憶部にローデータとして記憶されたエコー信号に基づいてBモード画像データを作成してもよい。
前記弾性画像処理部5についてもう少し詳しく説明すると、この弾性画像処理部5は、図2に示すように、物理量算出部51と弾性画像データ作成部52とを有している。前記物理量算出部51は、生体組織における各部の弾性に関する物理量として、前記超音波プローブ2による圧迫とその弛緩によって生じた生体組織における各部の変形による変位(以下、単に「変位」と云う)を算出する。前記物理量算出部51は、図3に示すように時間的に異なる二つのフレーム(i),(ii)に属する同一音線上における二つのエコー信号に基づいて変位を算出する。後述するように、前記エコー信号に相関ウィンドウW1,W2を設定し(図5参照)、これら相関ウィンドウW1,W2間で相関演算を行なって変位を算出する。そして、この変位に基づいて前記弾性画像データ作成部52によって一画素分の弾性画像データが作成される。前記物理量算出部51は本発明における物理量算出部の実施の形態の一例である。
前記弾性画像データ作成部52は、前記物理量算出部51によって算出された変位を色相情報に変換し、超音波の送受信面における弾性画像作成領域(本例では後述の関心領域R)についての弾性画像データを作成する。前記弾性画像データ作成部52は本発明における弾性画像データ作成部の実施の形態の一例である。
ちなみに、図3に示すように、異なる二つのフレーム(i),(ii)に属するエコー信号から一フレーム分の弾性画像データが作成される。一方で、Bモード画像データは、前記フレーム(i),(ii)のいずれかのエコー信号から作成される。
ここで、本例では、図4に示すように前記表示部7に表示されたBモード画像BG上に関心領域(ROI:Region Of Interest)Rが設定され、この関心領域Rについて前記弾性画像データが作成される。前記関心領域Rは本発明における弾性画像作成領域の実施の形態の一例である。ただし、本発明は、このように前記Bモード画像BGの一部について弾性画像を作成する場合に限られるものではなく、前記Bモード画像BGの全体について前記弾性画像データを作成してもよい。
前記Bモード画像処理部4で作成されたBモード画像データと、前記弾性画像処理部5で作成された弾性画像データは、前記合成部6で合成される。具体的には、この合成部6は、一フレーム分の前記Bモード画像データと前記弾性画像データとを加算処理し、前記表示部7に表示する一フレーム分の超音波画像データを作成する。そして、前記合成部6で得られた超音波画像データは、図4に示すように白黒のBモード画像BGとカラーの弾性画像EGとが合成された超音波画像Gとして前記表示部7に表示される。本例では、前記弾性画像EGは、前記関心領域R内に半透明で(背景のBモード画像が透けた状態で)表示される。
また、前記表示部7には、後述するクオリティ表示QGが表示される。前記表示部7は、本発明における報知部の実施の形態の一例である。
前記物理量平均部8は、一画素毎に算出された変位の前記関心領域Rにおける平均をフレーム毎に算出する。前記物理量平均部8の算出値を平均値XrAVとする。前記物理量平均部8は、本発明における物理量平均部の実施の形態の一例である。
前記比算出部9は、変位の平均の理想値XiAVに対する前記平均値XrAVの比Raを算出し、さらに後述するように(式1)の演算を行なう。前記比算出部9は、本発明における比較部及び比算出部の実施の形態の一例である。また、前記理想値XiAVは、本発明における予め設定された物理量の平均値の実施の形態の一例である。
ここで、前記理想値XiAVは、生体組織の弾性をより正確に反映した弾性画像を得ることができる強さで、超音波の送受信時に前記超音波プローブ2による生体組織への圧迫とその弛緩が行なわれた場合に、任意に設定される領域において得られる変位の平均値である。この理想値XiAVは、例えば腫瘍と同じ硬さの部分や正常組織と同じ硬さの部分などからなるファントム等を対象として実験を行ない、経験上得られる値である。また、この理想値XiAVは、操作者が前記操作部12において設定できるようになっていてもよいし、デフォルトとして装置に記憶されていてもよい。
前記クオリティ表示作成部10は、図4に示すように、超音波画像Gとともに前記表示部7に表示されるクオリティ表示QGを作成する。このクオリティ表示QGは、本例では横軸が時間、縦軸が、後述するようにフレーム毎に算出されるクオリティ値Qnを表すグラフgrからなる。前記クオリティ表示QGの作成については、後で詳述する。
前記制御部11は、CPU(Central Processing Unit)で構成され、図示しない記憶部に記憶された制御プログラムを読み出し、前記超音波診断装置1の各部における機能を実行させる。また、前記操作部9は、操作者が指示や情報を入力するためのキーボード及びポインティングデバイス(図示省略)などを含んで構成されている。
さて、本例の超音波診断装置1の作用について説明する。先ず、前記送受信部3は、前記超音波プローブ2から被検体の生体組織へ超音波を送信させ、そのエコー信号を取得する。このとき、前記超音波プローブ2により、被検体への圧迫とその弛緩を繰り返しながら超音波の送受信を行う。
そして、前記Bモード画像処理部4は、前記エコー信号に基づいてBモード画像データを作成する。また、前記弾性画像処理部5は、前記エコー信号に基づいて弾性画像データを作成する。前記Bモード画像データと前記弾性画像データは、前記合成部6で合成され、図3に示すようにBモード画像BGと弾性画像EGとが合成された超音波画像Gが前記表示部7に表示される。
また、前記表示部7には、前記超音波画像Gの下方に、前記クオリティ表示作成部10により作成されたクオリティ表示QGが表示される。
前記弾性画像処理部5における弾性画像データの作成及び前記クオリティ表示QGの作成について詳細に説明する。前記弾性画像データを作成するにあたり、前記物理量算出部51は、フレーム(i),(ii)に属するエコー信号のそれぞれに相関ウィンドウを設定する。具体的には、前記物理量算出部51は、図5に示すようにフレーム(i)に属するエコー信号に相関ウィンドウW1を設定し、フレーム(ii)に属するエコー信号に相関ウィンドウW2を設定する。そして、前記物理量算出部51は、前記相関ウィンドウW1,W2間で相関演算を行なって変位を算出する。
具体的に説明すると、図5において、前記フレーム(i),(ii)は、複数本の音線上において取得されたエコー信号からなる。図5では、前記フレーム(i)における複数本の音線の一部として、五本の音線L1a,L1b,L1c,L1d,L1eが示され、また前記フレーム(ii)において前記音線L1a〜L1eに対応する音線として、音線L2a,L2b,L2c,L2d,L2eが示されている。すなわち、前記音線L1a及び前記音線L2a、前記音線L1b及び前記音線L2b、前記音線L1c及び前記音線L2c、前記音線L1d及び前記音線L2d、前記音線L1e及び前記音線L2eは、異なる二つのフレームに属する同一音線に該当する。また、図5においてR(i),R(ii)は、前記関心領域Rに対応する領域を示している。
例えば、前記音線L1c上のエコー信号に、前記相関ウィンドウW1として相関ウィンドウW1cが設定され、前記音線L2c上のエコー信号に、前記相関ウィンドウW2として相関ウィンドウW2cが設定されたとする。前記物理量算出部51は、前記相関ウィンドウW1c,W2c間で相関演算を行ない、変位を算出する。前記物理量算出部51は、前記音線L1c,L2c上において、前記領域R(i),R(ii)の上端100から下端101まで相関ウィンドウW1c,W2cを順次設定し、変位を算出する。また、前記物理量算出部51は、前記領域R(i),R(ii)内の他の音線についても同様にして変位を算出する。そして、このようにして前記物理量算出部51によって変位が算出されると、この変位に基づいて前記弾性画像データ作成部52が弾性画像データを作成する。
次に、前記クオリティ表示GRの作成について説明する。このクオリティ表示GRの作成にあたり、先ず前記物理量平均部8が、前記関心領域R(前記領域R(i),R(ii))における変位の平均値XrAVを算出する。ちなみに、変位は負になることもあることから、前記平均値XrAVは負になることもあるものとする。次に、前記比算出部9が、XrAV/XiAVの演算を行ない、前記比Raを算出する。さらに、前記比算出部9は、前記比Raを次の(式1)に代入し、数値Yを得る。
Y=1.0−|log10|Ra||・・・(式1)
ここで、Yは、本発明において比較部による比較結果及び比較部の算出値の実施の形態の一例である。
ちなみに、この(式1)は、前記比Raを0から1までの範囲にするためのものであり、この(式1)で得られるYは、前記理想値XiAVに対する平均値XrAVの比と同等である。この(式1)で表される関数をグラフで表すと、図6に示すグラフとなる。この図6に示すように、0≦Y≦1となる。
また、0.1≦|Ra|≦10であるものとし、|Ra|がこの範囲を超えた場合、Yは零とする。
前記比算出部9の算出値Yは、前記クオリティ表示作成部10へ入力される。ここで、前記算出値Yはフレーム毎に算出される。前記クオリティ表示作成部10では、フレーム毎の前記算出値Yをクオリティ値Qnとしてプロットし、横軸が時間、縦軸が前記クオリティ値Qnを表すグラフgrからなるクオリティ表示QGを作成する。この時、前記クオリティ表示作成部10は、前記クオリティ値Qnの複数フレーム分の平均を算出し、この平均値をプロットしていってもよい。これにより、数値のばらつきのない安定したグラフgrを得ることができる。
0≦Y≦1であるため、0≦Qn≦1となる。クオリティ値Qnが1に近くなるほど、弾性画像のクオリティとしては良好であることを意味し、一方でクオリティ値Qnが0に近くなるほど、弾性画像のクオリティとしては悪くなることを意味する。ここで、弾性画像のクオリティが良好であるとは、生体組織の弾性をより正確に反映した弾性画像であることを意味し、一方で弾性画像のクオリティが悪いとは、生体組織の弾性を正確に反映した弾性画像ではないことを意味する。
クオリティ値Qnと弾性画像のクオリティとの関係についてより詳細に説明すると、図6のグラフから分かるように、前記平均値XrAVが前記理想値XiAVと等しい場合(すなわち、|Ra|が1)、YすなわちQnは1となる。従って、Qnが1、または1に近い値であれば、前記超音波プローブ2による生体組織に対する圧迫とその弛緩の度合いが適切であり、生体組織の弾性を正確に反映した弾性画像EGが得られていることになる。
一方で、前記平均値XrAVが前記理想値XiAVと離れた値になるほど(すなわち、|Ra|が1から離れた値になるほど)、Qnは零に近づく。ここで、前記平均値XrAVが前記理想値XiAVと離れた値になるということは、前記超音波プローブ2による生体組織に対する圧迫やその弛緩の度合いが足りない、または過剰であることを意味する。従って、Qnが零に近づくほど、生体組織に対する圧迫やその弛緩の度合いが足りないか、または過剰である結果、生体組織の弾性を正確に反映した弾性画像EGが得られていないことになる。
ちなみに、前記クオリティ値Qnが低いフレームについては、前記弾性画像EGの表示を行わないようにしてもよい。
前記クオリティ表示作成部10によって作成されたクオリティ表示QGは、前記合成部6において前記超音波画像Gと合成される。これにより、前記表示部7には前記超音波画像Gの下方に前記クオリティ表示QGが表示される。
前記クオリティ表示QGについてさらに詳細に説明すると、前記超音波画像Gが動画で表示される場合、前記クオリティ表示作成部10は、現在表示されている超音波画像Gにおけるクオリティ値Qnをフレーム毎にプロットすることにより、前記グラフgrを作成する。従って、前記表示部7において、前記グラフgrは、図7、図8、図9に示すように、時間の経過とともに左から右へ流れるように表示される。
ただし、前記超音波画像Gが、前記記憶部(図示省略)にローデータとして記憶されたエコー信号に基づいて作成される動画である場合、前記クオリティ表示作成部10は、再生の初めから終わりまでのグラフgrを作成し、これを前記表示部7に表示させてもよい。この場合、図10に示すように、前記クオリティ表示QGは、前記グラフgrのほかに、現在表示されている超音波画像Gがどの時間のフレームのものであるかを示す縦方向の線分bを含んでいてもよい。この線分bは、前記超音波画像Gの再生中に、左から右へ移動する(図中矢印の方向)。
ちなみに、前記線分bの長さは、前記クオリティ表示作成部10でフレーム毎に算出される前記クオリティ値Qnの最小値と最大値の間の長さになっている。
本例の超音波診断装置1によれば、前記理想値XiAVに対する前記平均値XrAVの比Raに基づいて算出される前記クオリティ値Qnの時間変化を表すグラフgrからなるクオリティ表示QGが表示されるので、操作者は、前記超音波プローブ2による生体組織に対する圧迫とその弛緩の度合いが足りなかったり、また過剰であったりしないかどうかを容易に判断することができる。これにより、生体組織の弾性を正確に反映した弾性画像であるか否かを従来よりも幅広い観点から評価することができる。
また、操作者は、前記グラフgrを見ることにより、クオリティ値Qnが高い所で前記超音波画像Gをフリーズし、この超音波画像Gを印刷等によって出力してもよい。これにより、生体組織の弾性をより正確に反映した超音波画像を印刷等によって出力することができる。さらに、前記超音波画像Gがリアルタイムで表示されている場合、操作者は、前記グラフgrを見ることによって前記超音波プローブ2による生体組織への圧迫とその弛緩の度合いを調節することもできる。
次に、第一実施形態の変形例について説明する。この変形例では、前記物理量平均部8は、相関係数C(0≦C≦1)が所定の閾値CTH以上である相関演算が行なわれた相関ウィンドウを選択してその変位の平均算出を行ない、平均値XrAV′を得る。そして、前記比算出部9が、前記平均値XrAV′を用いて前記比Raを算出し、また(式1)を用いてYを算出する。さらに、前記クオリティ表示作成部10が、前記算出値Yを用いて前記クオリティ表示QGを作成する。
前記平均値XrAV′は、エコー信号の強度が不十分な部分、生体組織の横ずれが生じている部分など、相関係数が低い部分の変位が除かれて得られた平均値である。従って、このような平均値XrAV′から得られたクオリティ値Qnは、前記超音波プローブ2による圧迫とその弛緩が適切な強さで行なわれているか否かを示すものとなる。以上より、操作者は、前記クオリティ表示QGから、前記超音波プローブ2による圧迫とその弛緩が適切な強さで行なわれているか否かをより正確に把握することができる。例えば、前記クオリティ値Qnが1から離れている場合、前記超音波プローブ2による圧迫とその弛緩が適切な強さで行なわれていないことを把握することができる。一方で、前記クオリティ値Qnが1或いは1に近い値であれば、操作者は前記超音波プローブ2による圧迫が適切な強さで行なわれていることを把握することができる。
また、仮に相関係数が低い相関演算で得られた変位を含めて前記平均値XrAVの算出を行った場合、前記超音波プローブ2による生体組織への圧迫とその弛緩の度合いが適切であっても、例えばエコー信号の強度が弱い場合は、前記平均値XrAVが小さくなり、前記クオリティ値Qnが1から離れてしまう。従って、この変形例のように、相関係数が低い部分の変位を除いて前記平均値XrAVの算出を行なうことにより、前記超音波プローブ2による生体組織への圧迫とその弛緩の度合いが適切であれば、常に前記クオリティ値Qnが1に近くなる。以上より、生体組織に対する圧迫とその弛緩の度合いが適切であるか否かをより正確に反映したクオリティ表示QGを表示させることができる。
(第二実施形態)
次に、第二実施形態について図11及び図12に基づいて説明する。なお、第一実施形態と同一の構成については同一の符号を付して説明を省略する。
本例の超音波診断装置20は、前記物理量平均部8及び比算出部9を備えておらず、代わりに相関係数平均部21を備えている。この相関係数平均部21は、本発明における相関係数平均部21の実施の形態の一例である。
本例の超音波診断装置20の作用について説明する。本例においては、前記クオリティ表示QGの作成方法が第一実施形態と異なっている。具体的に説明すると、前記相関係数平均部21は、前記クオリティ表示QGの作成にあたり、前記物理量算出部51によって行なわれた各相関演算における相関係数Cの関心領域R(領域R(i),R(ii))における平均値CAVをフレーム毎に算出する。本例では、この相関係数Cの平均値CAVをクオリティ値Qnとする。ここで、0≦C≦1であるので、本例においても、0≦Qn≦1である。相関演算における相関係数は、1に近づくほど生体組織の弾性をより正確に反映した変位を得ることができ、一方で零に近づくほど生体組織の弾性を正確に反映した変位を得ることができなくなる。従って、本例においても、Qnが1に近づくほど弾性画像EGのクオリティが良好になり、一方でQnが零に近づくほど弾性画像EGのクオリティが悪くなる。
そして、前記クオリティ表示作成部10は、前記平均値CAVを前記クオリティ値Qnとしてプロットし、前記グラフgrからなるクオリティ表示QGを作成する。この時、前記クオリティ表示作成部10は、第一実施形態と同様に、前記クオリティ値Qnの複数フレーム分の平均を算出し、この平均値をプロットしていってもよい。
本例においても、前記クオリティ表示QGを構成するグラフgrは、第一実施形態と同様に、時間の経過とともに左から右へ流れるように表示されてもよい。また、前記クオリティ表示QGは、前記グラフgrのほかに縦方向の線分bを含んでいてもよい。
本例の超音波診断装置20によれば、相関係数Cの平均値CAVであるクオリティ値Qnの時間変化を表すグラフgrからなるクオリティ表示QGが表示されるので、操作者は、表示されている弾性画像について、例えば生体組織に対する圧迫とその弛緩が過剰であったり、エコー信号の強度が不十分であったりすることなどに起因して相関係数が低い相関演算で得られた変位に基づいて作成された弾性画像データの画像であるか否かを把握することができる。これにより、生体組織の弾性画像を正確に反映した画像であるか否かを従来とは異なる観点から評価することができる。
(第三実施形態)
次に、第三実施形態について図13及び図14に基づいて説明する。なお、第一、第二実施形態と同一の構成については同一の符号を付して説明を省略する。
本例の超音波診断装置30は、前記物理量平均部8、前記比算出部9、前記クオリティ表示作成部10、前記相関係数平均部20等を備え、さらに乗算部31を備えている。前記乗算部31は、本発明における乗算部の実施の形態の一例である。
本例の超音波診断装置30における前記クオリティ表示QGの作成について説明する。前記物理量平均部8は、第一実施形態の変形例と同様に、相関係数Cが所定の閾値CTH以上である相関演算が行なわれた相関ウィンドウを選択してその変位の平均値XrAV′を算出し、また前記比算出部9が、前記平均値XrAV′を用いて前記比Raを算出し、前記(式1)からYを算出する。また、第二実施形態と同様に、前記相関係数平均部21が相関係数の平均値CAVを算出する。
そして、前記乗算部31は、前記比算出部9で得られた算出値Yと、前記相関係数平均部21で得られた相関係数の平均値CAVとを乗算し、乗算値Mを算出する。この乗算値Mはフレーム毎に算出される。
ここで、前記乗算部31は、前記算出値Yと前記相関係数の平均値CAVとを乗算する時に、重み付けをして乗算してもよい。
前記クオリティ表示作成部10は、前記乗算値Mを前記クオリティ値Qnとしてプロットし、前記クオリティ表示QGを作成する。この時、前記クオリティ表示作成部10は、第一、第二実施形態と同様に、前記クオリティ値Qnの複数フレーム分の平均を算出し、この平均値をプロットしていってもよい。
本例においても、前記クオリティ表示QGを構成するグラフgrは、第一、第二実施形態と同様に、時間の経過とともに左から右へ流れるように表示されてもよい。また、前記クオリティ表示QGは、前記グラフgrのほかに縦方向の線分bを含んでいてもよい。
ここで、0≦Y≦1、0≦CAV≦1であるので、0≦M≦1となる。従って、本例においても、0≦Qn≦1である。前記乗算値Mは、前記算出値Yと前記相関係数の平均値CAVとの乗算値であるため、乗算値M、すなわちQnが1に近づくほど弾性画像EGのクオリティが良好になり、一方でQnが零に近づくほど弾性画像EGのクオリティが悪くなる。
ここで、第一実施形態の変形例のように、所定の閾値CTH以上の相関係数の相関演算で得られた変位の平均値XrAV′から算出されたクオリティ値Qnを前記クオリティ表示QGとして表示すると、相関係数は弾性画像のクオリティの評価の要素として全く反映されないことになる。一方で、第二実施形態のように、相関係数Cの平均値CAVを前記クオリティ表示QGとして表示すると、前記超音波プローブ2による生体組織への圧迫とその弛緩の度合いが足りなかったとしても、相関係数Cとしては高くなるために、前記クオリティ値Qnとしては良好な値が表示されることがある。従って、本例では、前記平均値XrAV′を用いて算出された前記比Raを用いて得られる算出値Yと前記相関係数Cの平均値CAVとを乗算することにより、生体組織への圧迫とその弛緩の度合いの要素と、相関係数の要素とを加味したクオリティ表示QGを表示することができる。これにより、生体組織の弾性を正確に反映した弾性画像であるか否かを、従来よりも幅広い観点から評価することができる。
(第四実施形態)
次に、第四実施形態について図15及び図16に基づいて説明する。なお、第一〜第三実施形態と同一の構成については同一の符号を付して説明を省略する。
本例の超音波診断装置40は、第三実施形態の超音波診断装置30と基本的構成は同じであるが、前記比算出部9で得られる算出値Y、前記相関係数平均部21で得られる相関係数の平均値CAV及び前記乗算部31で得られる乗算値Mのいずれかがクオリティ値Qnとして選択され、前記クオリティ表示作成部10において、クオリティ表示QGが表示されるようになっている。
具体的には、前記操作部12において、前記算出値Y、前記相関係数の平均値CAV及び前記乗算値Mのいずれかを用いて前記クオリティ表示QGの表示を行なうかが操作者により入力される。そして、この操作部12における指示入力に基づいて、前記制御部11が制御を行なって、前記比算出部9による算出値Yの算出、前記相関係数平均部21による相関係数の平均値CAVの算出及び前記乗算部31による乗算値Mの算出のいずれかを行なわせ、得られた算出値に基づいて、前記クオリティ表示部10によるクオリティ表示QGの作成が行なわれる。これにより、操作者により選択された種類の値を用いて作成されたクオリティ表示QGが前記表示部7に表示される。
操作者は、前記操作部12の指示入力を行なって、前記算出値Y、前記相関係数の平均値CAV及び前記乗算値Mの中から一旦選択した値の種類を変更してもよい。このように、前記操作部12によって切替の指示入力を行なうことにより、新たに選択された種類の値を用いて作成されたクオリティ表示QGが前記表示部7に表示される。
本例の超音波診断装置40によれば、前記比算出部9で得られる算出値Yを用いて作成されたクオリティ表示QG、前記相関係数平均部21で得られる相関係数の平均値CAVを用いて作成されたクオリティ表示QG、前記乗算部31で得られる乗算値Mを用いて作成されたクオリティ表示QGを切り替えて表示させることができるので、生体組織の弾性を正確に反映した弾性画像であるか否かを、従来よりも幅広い観点から評価することができる。
以上、本発明を前記各実施形態によって説明したが、本発明はその主旨を変更しない範囲で種々変更実施可能なことはもちろんである。例えば、前記物理量算出部51は、生体組織の弾性に関する物理量として、生体組織の変形による変位の代わりに生体組織の歪みや弾性率を算出してもよい。
また、前記比算出部9では、前記比Raのみを算出し、(式1)の演算を行わなくてもよい。この場合、前記クオリティ表示作成部10は、前記比|Ra|を前記クオリティ値Qnとしてプロットしてなるグラフgrを、前記クオリティ表示QGとして作成してもよい。
前記比|Ra|を前記クオリティ値Qnとしてプロットして作成され、前記表示部7に表示されるクオリティ表示QGの一例を図17に示す。図17において、横軸は時間、縦軸は比|Ra|である。この図17に示すように、前記比|Ra|が1に近い所定の範囲に、帯状の部分Oを表示してもよい。この帯状の部分Oは、生体組織の弾性を正確に反映した弾性画像EGが得られる比|Ra|の範囲に設定される。このような帯状の部分Oを表示することにより、クオリティ表示QGがこの帯状の部分Oに入るように、操作者が前記超音波プローブ2による生体組織への圧迫と弛緩を行なえば、生体組織の弾性を正確に反映した弾性画像を得ることができる。
また、前記クオリティ表示QGは、グラフgrからなるものに限られず、例えば図18に示すように、バーBからなるものであってもよい。このバーBは、縦方向の長さが前記クオリティ値Qnの値(0≦Qn≦1)に相当し、クオリティ値Qnの変化とともに、縦方向に伸縮する。
また、バーBは、クオリティ値Qnに応じて縦方向に伸縮するものではなく、クオリティ値Qnに応じて色が変化するものであってもよい。
その他、クオリティ表示QGは、前記表示部7に数値で表示されてもよい。さらに、前記クオリティ値Qnをクオリティ表示QGとして表示するものに限られない。例えば、クオリティ値Qnを音として発するためのスピーカー(図示省略)を備えていてもよい。このスピーカーは、本発明における報知部の実施の形態の一例である。この場合には、クオリティ値Qnの高低を、音の高低で表すようにする。
1,20,30,40 超音波診断装置
7 表示部(報知部)
8 物理量平均部
9 比算出部(比較部)
12 操作部
21 相関係数平均部
31 乗算部
51 物理量算出部
52 弾性画像データ作成部

Claims (10)

  1. 生体組織に対する超音波の送受信により得られた同一音線上の時間的に異なる二つのエコー信号に相関ウィンドウを設定し、該相関ウィンドウ間で相関演算を行なって生体組織における各部の弾性に関する物理量を算出する物理量算出部と、
    前記物理量に基づいて、生体組織の弾性画像データを、超音波の送受信面における弾性画像作成領域について作成する弾性画像データ作成部と、
    前記弾性画像作成領域における前記物理量の平均をフレーム毎に算出する物理量平均部と、
    該物理量平均部による算出値を、予め設定された前記物理量の平均値と比較する比較部と、
    該比較部による比較結果を報知する報知部と、
    を備えることを特徴とする超音波診断装置。
  2. 前記物理量平均部は、所定の閾値以上の相関係数の相関演算が行なわれた相関ウィンドウについて得られた物理量の平均算出を行なうことを特徴とする請求項1に記載の超音波診断装置。
  3. 前記比較部は、前記比較結果として、予め設定された前記物理量の平均値に対する前記物理量平均部による算出値の比を算出することを特徴とする請求項1又は2に記載の超音波診断装置。
  4. 前記報知部には、前記比較部によってフレーム毎に得られる比較結果が複数フレームにわたって平均化されて報知されることを特徴とする請求項1〜3のいずれか一項に記載の超音波診断装置。
  5. 生体組織に対する超音波の送受信により得られた同一音線上の時間的に異なる二つのエコー信号に相関ウィンドウを設定し、該相関ウィンドウ間で相関演算を行なって生体組織における各部の弾性に関する物理量を算出する物理量算出部と、
    前記物理量に基づいて、生体組織の弾性画像データを、超音波の送受信面における弾性画像作成領域について作成する弾性画像データ作成部と、
    前記相関ウィンドウ間の相関演算における相関係数の前記弾性画像作成領域における平均をフレーム毎に算出する相関係数平均部と、
    該相関係数平均部による算出結果を報知する報知部と、
    を備えることを特徴とする超音波診断装置。
  6. 前記報知部には、前記相関係数平均部によってフレーム毎に得られる算出結果が複数フレームにわたって平均化されて報知されることを特徴とする請求項5に記載の超音波診断装置。
  7. 生体組織に対する超音波の送受信により得られた同一音線上の時間的に異なる二つのエコー信号に基づいて生体組織における各部の弾性に関する物理量を算出する物理量算出部と、
    前記物理量に基づいて、生体組織の弾性画像データを、超音波の送受信面における弾性画像作成領域について作成する弾性画像データ作成部と、
    所定の閾値以上の相関係数の相関演算が行なわれた相関ウィンドウについて得られた物理量の前記弾性画像作成領域における平均をフレーム毎に算出する物理量平均部と、
    予め設定された前記物理量の平均値に対する前記物理量平均部による算出値の比を算出する比算出部と、
    前記相関ウィンドウ間の相関演算における相関係数の前記弾性画像作成領域における平均をフレーム毎に算出する相関係数平均部と、
    前記比算出部の算出値と、前記相関係数平均部の算出値とを乗算する乗算部と、
    該乗算部による乗算結果を報知する報知部と、
    を備えることを特徴とする超音波診断装置。
  8. 前記報知部には、前記乗算部によってフレーム毎に得られる乗算結果が複数フレームにわたって平均化されて報知されることを特徴とする請求項7に記載の超音波診断装置。
  9. 前記乗算部は、前記比算出部の算出値と、前記相関係数平均部の算出値との重み付け演算を行なうことを特徴とする請求項7又は8に記載の超音波診断装置。
  10. 生体組織に対する超音波の送受信により得られた同一音線上の時間的に異なる二つのエコー信号に基づいて生体組織における各部の弾性に関する物理量を算出する物理量算出部と、
    前記物理量に基づいて、生体組織の弾性画像データを、超音波の送受信面における弾性画像作成領域について作成する弾性画像データ作成部と、
    所定の閾値以上の相関係数の相関演算が行なわれた相関ウィンドウについて得られた物理量の前記弾性画像作成領域における平均をフレーム毎に算出する物理量平均部と、
    予め設定された前記物理量の平均値に対する前記物理量平均部による算出値の比を算出する比算出部と、
    前記相関ウィンドウ間の相関演算における相関係数の前記弾性画像作成領域における平均をフレーム毎に算出する相関係数平均部と、
    前記比算出部の算出値と、前記相関係数平均部の算出値とを乗算する乗算部と、
    前記比算出部による算出結果、前記相関係数平均部による算出結果又は前記乗算部による乗算結果を切り替えて報知可能な報知部と、
    該報知部における切替の指示を入力する操作部と、
    を備えることを特徴とする超音波診断装置。
JP2009195004A 2009-08-26 2009-08-26 超音波診断装置 Active JP5484826B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2009195004A JP5484826B2 (ja) 2009-08-26 2009-08-26 超音波診断装置
CN201010273768.8A CN101999907B (zh) 2009-08-26 2010-08-25 超声诊断设备
US12/868,926 US8684931B2 (en) 2009-08-26 2010-08-26 Ultrasonic diagnostic apparatus for elasticity imaging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009195004A JP5484826B2 (ja) 2009-08-26 2009-08-26 超音波診断装置

Publications (2)

Publication Number Publication Date
JP2011045457A true JP2011045457A (ja) 2011-03-10
JP5484826B2 JP5484826B2 (ja) 2014-05-07

Family

ID=43625873

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009195004A Active JP5484826B2 (ja) 2009-08-26 2009-08-26 超音波診断装置

Country Status (3)

Country Link
US (1) US8684931B2 (ja)
JP (1) JP5484826B2 (ja)
CN (1) CN101999907B (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014530055A (ja) * 2011-09-27 2014-11-17 コーニンクレッカ フィリップス エヌ ヴェ 超音波エラストグラフィシステム及び方法

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1864612A4 (en) * 2005-03-30 2009-10-28 Hitachi Medical Corp ULTRASOUND DEVICE
FR2899336B1 (fr) 2006-03-29 2008-07-04 Super Sonic Imagine Procede et dispositif pour l'imagerie d'un milieu viscoelastique
EP2146640B1 (en) * 2007-05-16 2018-05-23 Super Sonic Imagine Method and device for measuring a mean value of visco-elasticity of a region of interest
WO2011010626A1 (ja) * 2009-07-24 2011-01-27 株式会社 日立メディコ 超音波診断装置、弾性画像の保存/再生方法、及び弾性画像の保存/再生プログラム
EP2620102A4 (en) * 2010-09-21 2016-12-07 Hitachi Ltd ULTRASONIC DIAGNOSIS DEVICE AND METHOD FOR DISPLAYING ULTRASONIC IMAGES
JP5294340B2 (ja) 2010-10-27 2013-09-18 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置
JP5209026B2 (ja) 2010-10-27 2013-06-12 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置
JP5373749B2 (ja) 2010-11-30 2013-12-18 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 医用画像表示装置及びその制御プログラム
CN102824193B (zh) * 2011-06-14 2016-05-18 深圳迈瑞生物医疗电子股份有限公司 一种弹性成像中的位移检测方法、装置及系统
JP5774448B2 (ja) 2011-10-28 2015-09-09 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置
JP5774560B2 (ja) * 2012-08-20 2015-09-09 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びその制御プログラム
JP5879230B2 (ja) * 2012-08-21 2016-03-08 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びその制御プログラム
JP5710566B2 (ja) * 2012-09-18 2015-04-30 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びその制御プログラム
CN103845081B (zh) * 2012-11-28 2018-04-10 深圳迈瑞生物医疗电子股份有限公司 超声弹性成像系统和方法、实时动态帧间处理方法
JP5863628B2 (ja) 2012-11-30 2016-02-16 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びその制御プログラム
US9332962B2 (en) * 2013-03-13 2016-05-10 Siemens Medical Solutions Usa, Inc. Ultrasound ARFI displacement imaging using an adaptive time instance
JP2015008733A (ja) * 2013-06-26 2015-01-19 ソニー株式会社 超音波処理装置および方法
WO2016061802A1 (zh) * 2014-10-24 2016-04-28 通用电气公司 一种在当前超声影像上显示感兴趣区域的方法及装置
EP3639752A1 (en) * 2018-10-03 2020-04-22 Canon Medical Systems Corporation Analyzing apparatus and analyzing program
JP6782747B2 (ja) * 2018-10-24 2020-11-11 ゼネラル・エレクトリック・カンパニイ 超音波装置及びその制御プログラム
FR3104736B1 (fr) * 2019-12-13 2022-12-09 Supersonic Imagine Procédé ultrasonore pour quantifier l’élasticité non linéaire par ondes de cisaillement d’un milieu, et dispositif pour mettre en œuvre ce procédé
CN113939232B (zh) * 2020-05-14 2024-06-28 深圳迈瑞生物医疗电子股份有限公司 一种超声成像设备及方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11128222A (ja) * 1997-11-04 1999-05-18 Aloka Co Ltd 超音波診断装置
JP2005118152A (ja) * 2003-10-14 2005-05-12 Hitachi Medical Corp 超音波診断装置
WO2005122907A1 (ja) * 2004-06-22 2005-12-29 Hitachi Medical Corporation 超音波診断装置および弾性画像表示方法
JP2006166955A (ja) * 2004-12-13 2006-06-29 Matsushita Electric Ind Co Ltd 超音波診断装置
JP3991282B2 (ja) * 2004-08-05 2007-10-17 株式会社日立メディコ 弾性像表示方法及び超音波診断装置
JP2008073144A (ja) * 2006-09-20 2008-04-03 Aloka Co Ltd 超音波診断装置
JP2008126079A (ja) * 2006-11-22 2008-06-05 General Electric Co <Ge> 組織の弾力特性を計測するための直接式歪み推定器
JP2008161546A (ja) * 2006-12-28 2008-07-17 Matsushita Electric Ind Co Ltd 超音波診断装置
JP2008259555A (ja) * 2007-04-10 2008-10-30 Hitachi Medical Corp 超音波診断装置

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5107837A (en) 1989-11-17 1992-04-28 Board Of Regents, University Of Texas Method and apparatus for measurement and imaging of tissue compressibility or compliance
SE504575C2 (sv) * 1994-10-06 1997-03-10 Lorentzen & Wettre Ab Anordning för ultraljudsmätning av elastiska egenskaper hos en pappersbana i rörelse
US6270459B1 (en) * 1998-05-26 2001-08-07 The Board Of Regents Of The University Of Texas System Method for estimating and imaging of transverse displacements, transverse strains and strain ratios
US6558324B1 (en) 2000-11-22 2003-05-06 Siemens Medical Solutions, Inc., Usa System and method for strain image display
JP4257696B2 (ja) 2002-10-28 2009-04-22 株式会社日立メディコ 超音波装置
JP4314035B2 (ja) 2003-01-15 2009-08-12 株式会社日立メディコ 超音波診断装置
US7914456B2 (en) 2003-05-30 2011-03-29 Hitachi Medical Corporation Ultrasonic probe and ultrasonic elasticity imaging device
JP5203605B2 (ja) * 2004-06-18 2013-06-05 株式会社日立メディコ 超音波診断装置
US7678051B2 (en) 2005-09-27 2010-03-16 Siemens Medical Solutions Usa, Inc. Panoramic elasticity ultrasound imaging
JP5356140B2 (ja) * 2009-07-22 2013-12-04 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びその制御プログラム
JP4999969B2 (ja) * 2010-07-13 2012-08-15 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びその制御プログラム

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11128222A (ja) * 1997-11-04 1999-05-18 Aloka Co Ltd 超音波診断装置
JP2005118152A (ja) * 2003-10-14 2005-05-12 Hitachi Medical Corp 超音波診断装置
WO2005122907A1 (ja) * 2004-06-22 2005-12-29 Hitachi Medical Corporation 超音波診断装置および弾性画像表示方法
JP3991282B2 (ja) * 2004-08-05 2007-10-17 株式会社日立メディコ 弾性像表示方法及び超音波診断装置
JP2006166955A (ja) * 2004-12-13 2006-06-29 Matsushita Electric Ind Co Ltd 超音波診断装置
JP2008073144A (ja) * 2006-09-20 2008-04-03 Aloka Co Ltd 超音波診断装置
JP2008126079A (ja) * 2006-11-22 2008-06-05 General Electric Co <Ge> 組織の弾力特性を計測するための直接式歪み推定器
JP2008161546A (ja) * 2006-12-28 2008-07-17 Matsushita Electric Ind Co Ltd 超音波診断装置
JP2008259555A (ja) * 2007-04-10 2008-10-30 Hitachi Medical Corp 超音波診断装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014530055A (ja) * 2011-09-27 2014-11-17 コーニンクレッカ フィリップス エヌ ヴェ 超音波エラストグラフィシステム及び方法

Also Published As

Publication number Publication date
CN101999907B (zh) 2015-09-23
CN101999907A (zh) 2011-04-06
US20110054314A1 (en) 2011-03-03
JP5484826B2 (ja) 2014-05-07
US8684931B2 (en) 2014-04-01

Similar Documents

Publication Publication Date Title
JP5484826B2 (ja) 超音波診断装置
JP4455003B2 (ja) 超音波診断装置
JP4999969B2 (ja) 超音波診断装置及びその制御プログラム
JP2007282932A (ja) 弾性画像生成方法及び超音波診断装置
KR101666324B1 (ko) 초음파 진단 장치
JP2010099378A (ja) 超音波診断装置
JP5280379B2 (ja) 超音波診断装置、超音波弾性情報処理方法及び超音波弾性情報処理プログラム
JP5535575B2 (ja) 超音波診断装置
JP5356140B2 (ja) 超音波診断装置及びその制御プログラム
JP5473527B2 (ja) 超音波診断装置
JP5415669B2 (ja) 超音波診断装置
JP5489178B2 (ja) 超音波診断装置
JP2012061075A (ja) 超音波診断装置及びその制御プログラム
JP5535574B2 (ja) 超音波診断装置
KR101329944B1 (ko) 초음파 진단 장치
JP5623609B2 (ja) 超音波診断装置
JP5485418B2 (ja) 超音波診断装置
JP2010082308A (ja) 超音波診断装置
JP5389497B2 (ja) 超音波診断装置およびその制御プログラム
JP5449896B2 (ja) 超音波診断装置
KR101574821B1 (ko) 초음파 진단 장치 및 그 제어 프로그램
JP5394693B2 (ja) 超音波診断装置
JP2010246818A (ja) 超音波診断装置及びその制御プログラム

Legal Events

Date Code Title Description
A625 Written request for application examination (by other person)

Free format text: JAPANESE INTERMEDIATE CODE: A625

Effective date: 20120530

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120620

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130729

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130731

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20131024

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140219

R150 Certificate of patent or registration of utility model

Ref document number: 5484826

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250