JP4847665B2 - 医療診断用超音波イメージングシステム及び被送信超音波ビームの音響出力パラメータを決定する方法 - Google Patents
医療診断用超音波イメージングシステム及び被送信超音波ビームの音響出力パラメータを決定する方法 Download PDFInfo
- Publication number
- JP4847665B2 JP4847665B2 JP2001557472A JP2001557472A JP4847665B2 JP 4847665 B2 JP4847665 B2 JP 4847665B2 JP 2001557472 A JP2001557472 A JP 2001557472A JP 2001557472 A JP2001557472 A JP 2001557472A JP 4847665 B2 JP4847665 B2 JP 4847665B2
- Authority
- JP
- Japan
- Prior art keywords
- ultrasound
- output parameter
- ultrasound image
- imaging system
- region
- 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.)
- Expired - Lifetime
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/42—Details of probe positioning or probe attachment to the patient
- A61B8/4272—Details of probe positioning or probe attachment to the patient involving the acoustic interface between the transducer and the tissue
- A61B8/4281—Details of probe positioning or probe attachment to the patient involving the acoustic interface between the transducer and the tissue characterised by sound-transmitting media or devices for coupling the transducer to the tissue
-
- 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/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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/56—Details of data transmission or power supply
-
- 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
-
- 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/5205—Means for monitoring or calibrating
-
- 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
-
- 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/52076—Luminous indicators
-
- 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/52079—Constructional features
- G01S7/52084—Constructional features related to particular user interfaces
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Human Computer Interaction (AREA)
- Acoustics & Sound (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Description
(背景技術)
連邦薬品局(FDA)は、患者に進入する超音波ビームのピークレアフラクショナル圧力(peak rarefractional pressure)を指定レベル以下にすることを要求している。この要求に合致させるために、イメージング場内のピーク音響圧力に関係する機械的インデックスを表示するようになっている医療診断用超音波イメージングシステムが多い。表示された機械的インデックスは、コントラストイメージング検査を準備して遂行するためにも使用することができる。コントラスト剤の非線形レスポンス(調波または破壊)は、部分的に、超音波の音響圧力に依存する。もし非線形レスポンスを望まなければ、超音波システムのユーザは、例えば送信電力を低下させて表示された機械的インデックスをあるレベルまで低下させ、コントラスト剤の望ましくないレスポンスを最小にすることができる。しかしながら、表示された機械的インデックスが、イメージング場内のコントラスト剤が存在している位置には無関係になる可能性がある。従って、表示された機械的インデックスを使用することは、関連圧力の粗な測度でしかない場合が多く、イメージング状態をそれ程最適ではなくしてしまう可能性がある。例えば、コントラスト剤からの非線形レスポンスを回避するために、送信電力を必要なレベルより低いレベルまで低下させることはできるが、それによって、イメージ品質を望ましくなく犠牲にすることになる。また、現在のアプローチを用いる場合、イメージング手順を最適化するために、患者内へコントラスト剤を複数回注入する必要があり得る。更に、2つの関心領域からのコントラスト剤のレスポンスを比較する際に、表示された機械的インデックスに対応付けられた空間的曖昧さが誤差をもたらす恐れがある。
【0002】
従って、上述した欠陥を解消する医療診断用超音波イメージングシステムに対する要望が存在している。
【0003】
(発明の概要)
本発明は、特許請求の範囲によって限定されるものであり、以下の説明はこれらの特許請求の範囲を何等制限するものではない。
【0004】
序説として以下に説明する好ましい実施の形態は、医療診断用超音波イメージングシステム及び送信される(以下、被送信という)超音波ビームの音響出力パラメータを決定する方法を提供する。一つの好ましい実施の形態では、超音波システムは、ユーザが選択した領域における被送信超音波ビームの音響出力パラメータを決定する。別の好ましい実施の形態では、超音波システムは、超音波イメージングシステムの動作パラメータを自動的に調整することによって、選択された領域における被送信超音波ビームの指定された音響出力パラメータを達成する。更に別の好ましい実施の形態では、被送信超音波ビームのピーク音響出力パラメータを含まない領域が超音波イメージ内に選択される。次いでシステムは、その領域における被送信超音波ビームの音響出力パラメータを決定し、決定された音響出力パラメータの指示を供給する。
【0005】
以下に、添付図面を参照して好ましい実施の形態を説明する。
【0006】
(実施の形態)
図1は、現在では好ましい実施の形態の医療診断用超音波イメージングシステム100及びトランスデューサ105のブロック図である。超音波システム100はどのような適当なイメージングモード(例えば、Bモードイメージング、ドップラーイメージング、組織調波(ハーモニック)イメージング、コントラスト剤調波イメージング等)でも使用することができ、またトランスデューサ105はどのような型(例えば、1D、1.5D、平・凹、単要素、フェーズドアレイ等)であっても差し支えない。トランスデューサ105は、送信ビームフォーマー110及び受信ビームフォーマー115に結合されている。本明細書に使用する用語“…に結合されている”とは、直接的に結合されていること、または1つまたはそれ以上の成分を通して間接的に結合されていることを意味する。
【0007】
各ビームフォーマー110、115はプロセッサ120に結合されており、プロセッサ120はユーザインタフェース125及びディスプレイ130に結合されている。用語“プロセッサ”は、広義に、以下に説明する好ましい実施の形態を実現するために使用することができる超音波システム100の適切なハードウェア及び/またはソフトウェア成分を意味している。どのような適切なハードウェア(アナログまたはデジタル)またはソフトウェアも使用できること、及び以下に説明する実施の形態はハードウェアを用いることなく実現できることを理解されたい。更に、プロセッサ120は、簡易化のために図1には示されていない超音波システム100(付随プロセッサを含む)の他のプロセッサと分離することも、またはそれらと(部分的に、または全体的に)組合わせることも可能である。
【0008】
動作を説明する。プロセッサ120は、送信ビームフォーマー110からトランスデューサ105へ電圧を印加せしめてそれを振動させ、人の組織(即ち、患者の身体)のようなオブジェクト140内へ超音波ビーム135を放出させる。身体から反射した超音波エネルギはトランスデューサ105上に衝突し、トランスデューサ105から得られた電圧は受信ビームフォーマー115によって受信される。プロセッサ120はこれらの検知された電圧を処理し、反射信号に関連する超音波イメージを作り、ディスプレイ130上にイメージを表示する。典型的には、超音波イメージを生成するために、幾つかの超音波ビームが使用される。ユーザインタフェース125は、例えば送信、受信、及び表示動作に使用されるパラメータを調整するために使用することができる。超音波イメージングシステム100が、付加的な成分を含むことができることを理解されたい。
【0009】
超音波システム100は、以下に説明するように、被送信超音波ビームの音響出力パラメータの決定及び/または較正に関する1つまたはそれ以上の動作を遂行するように動作する。被送信超音波ビームの“音響出力パラメータ”なる用語は、広義に、医療診断用超音波イメージングシステムのトランスデューサから放出される超音波ビームの音響出力パラメータをカバーすることを意図している。音響出力パラメータが熱的及び/または機械的音響出力であり、音響出力パラメータはコントラスト剤を変更するように影響を与えることができる(例えば、音響出力パラメータがコントラスト剤に非線形レスポンスを生じさせることができる)ことが好ましい。熱的音響出力インデックスは、限定するものではないが、音響パワー、音響エネルギ、熱的インデックス(TI)、骨熱的インデックス(TIB)、頭骨熱的インデックス(TIC)、柔組織熱的インデックス(TIS)、及びパルス強度積分(PII)を含む。機械的音響出力のインデックスは、限定するものではないが、圧力(圧縮またはレアフラクショナル,rarefractional)、瞬時空間ピーク時間平均(ISPTA)、及び普通は送信周波数の平方根で除した被送信超音波ビームのピークレアフラクショナル圧力として定義される機械的インデックスを含む。
【0010】
超音波システム100は、図2のフローチャートに示すように、超音波イメージ内のユーザによって選択された領域における被送信超音波ビームの音響出力パラメータを決定するための方法に使用することができる。図3−5は、この方法の説明を援助するための超音波イメージである。先ず、超音波ビームがトランスデューサ105から送信され(ステップ210)、次いで超音波システム100は図3に示す超音波イメージ300を生成する(ステップ215)。次に、超音波システム100は、超音波イメージ内のある領域の選択をユーザから受信する(ステップ220)。ここに使用する“領域”は、その領域を1つの線、または任意のまたは予め定められた形状によって限定される単一の点であることも、または複数の点であることもできる。例えば、図3の超音波イメージ300においては、領域はキャリパー305によって指示されている点であり、一方図4及び5の超音波イメージ400においては、領域405は五角形によって囲まれた複数の点である。また“領域”は、ある線内の、または任意の形状または予め定められた形状によって囲まれた点の部分集合ということもできる。例えば“領域”は、ある線上の点、またはある形状によって囲まれた1つの点、幾つかの点、または全ての点であることができる。ユーザは、例えばユーザインタフェース125(例えば、トラックボール、マウス、キーボード、タッチパッド、タッチスクリーン、音声認識装置等)と対話してディスプレイ130上に表示された超音波イメージ上にカーソルその他の視覚標識を位置決めし、ある領域を選択することができる。もしカーソルの下に位置するイメージされつつある組織が運動していれば、カーソルは、その組織内のオリジナルの位置を追跡するようにフレーム毎に自動的に再位置決めすることができる。
【0011】
超音波システム100はユーザが選択した領域を受けた後に、ユーザが選択し領域における被送信超音波ビームの音響出力パラメータを決定する(ステップ225)。上述したように、“領域”は、ある線によって限定された、または任意のまたは予め定められた形状によって囲まれた点の部分集合であることができる。例えば、図4において、音響出力パラメータが決定される“領域”は、五角形405によって囲まれた点の1つ、幾つか、または全てであることができる。音響出力パラメータは、測定、計算、推定、予測、または他の何等かの適当な方法によって“決定”することができる。2つの文献、即ち、American Institute of Ultrasound in Medicine and National Electrical Manufacturers Associationから刊行された“Standard for Real-Time Display of Thermal and Mechanical Acoustic Output Indices on Diagnostic Ultrasound Equipment”Revision 1 (1998)、及び“Acoustic Output Measurement Standard for Diagnostic Ultrasound Equipment”(1998) に、送信される音場のパラメータを決定するために使用することができる適当な方法が記載されているので参照されたい。被送信超音波ビームの音響出力パラメータを決定するために使用されるハードウェア及び/またはソフトウェアは、従来の超音波システムにおいて機械的インデックスを決定するために使用されているものと同一であることができる。異なる点は、音響出力パラメータを、固定されたピーク圧力位置ではなく、ユーザが選択した領域内で決定することである。この決定動作は、捕捉されたフレームまたはクリップ上での取得中に、または取得後に行うことができる。取得後に決定する場合には、超音波システム100または外部解析及び定量化システムの何れかによる取得後決定を援助するために、超音波システムの動作状態を記録しておくことが好ましい。
【0012】
音響出力パラメータを決定した後に、超音波システム100は、決定した音響出力パラメータの指示を供給することができる(ステップ230)。この指示は、限定するものではないが、視覚、聴覚、または触覚指示を含むどのような適当な形状であることもできる。この指示は、超音波システム自体に供給する(例えば、表示されたイメージ上に、またはシステム自体に視覚指示を(例えば、LEDを介して)供給する)ことも、または超音波システムの外部デバイスによって供給する(例えば、その指示を外部スピーカーによって発声させる)こともできる。勿論、決定された音響出力パラメータがあるしきい値に等しいか、それより上にあるか、またはそれより下にあることを、音、光、または触覚のような標識を用いて指示を供給する他の方法を使用することができる。
【0013】
図3−5のディスプレイは、決定された音響出力パラメータを指示することができる種々の形状の幾つかを示している。これらの図においては、音響出力パラメータは、被送信超音波ビームの圧力である。この特定の音響出力パラメータの使用は単に説明の目的からに過ぎず、本発明の範囲を限定する意図は毛頭ない。決定された圧力と機械的インデックスとの間を有用に対応付けるために、決定された圧力を、送信周波数の平方根で除すことが好ましい。これらの図では、この量をコントラスト剤機械的インデックス(CA MI)と称している。図3においては、ユーザが選択した領域は単一の点であり、決定されたCA MIがイメージと共に表示されている。図3に示すように、機械的インデックス(MI)も表示されている。これらの好ましい実施の形態を実現するためには、機械的インデックスは必ずしも表示する必要はないが、国の規則に従うために機械的インデックスを表示することが好ましいであろう。勿論、複数の音響出力パラメータを決定して表示することができる。
【0014】
図3においては、ユーザが選択した領域は単一の点であり、供給される指示はその単一の点に対応付けられた単一の値である。図4及び5に示すように、もしユーザが選択した領域が複数の点であれば、供給される指示はその領域内の点の1つ、幾つか、または全てのための指示であることができる。図4及び5においては、ユーザが選択した領域は、五角形405によって囲まれた複数の点からなっている。これらの複数の点に関して決定された音響出力パラメータの分布(例えば、最大、最小、平均、モード、分散)を決定し、指示することができる。例えば、図4には、決定された音響出力パラメータの平均値が表示されており(CA MI MEAN)、また図4及び5には最大値(CA MI MAX)が表示されている。供給される指示は、ユーザが選択した領域において決定された音響出力パラメータの範囲を指示するために、一次元、二次元、または三次元の等圧線、またはマップ(グレースケール、またはカラー)の形状であることもできる。例えば、ユーザが選択した領域内の点または面積を、決定された音響出力パラメータがあるしきい値に等しいか、それより上にあるか、またはそれより下にあるものとして指示することができる。図5では等圧線表現が使用されており、五角形405は輪郭線410によって満たされている。図5には、最大輪郭の値(CA MI MAX=1.0)及び輪郭間隔(CA MI Interval=0.1)が指示されている。
【0015】
これらの好ましい実施の形態は幾つかの長所を伴っている。例えば、音響出力パラメータはユーザが選択した領域において決定されるので、従来の超音波システムにおける機械的インデックスの決定とは異なり、パラメータが決定される時にイメージ内の何処にも空間的曖昧さは存在しない。またこれらの好ましい実施の形態は、音響出力パラメータの空間的非均一性の記述を供給することができ、トランスデューサ、モード、周波数、及びイメージング応用全体にまたがる音響出力パラメータの矛盾のない最適化を可能にする。
【0016】
これらの好ましい実施の形態は、コントラスト剤イメージング応用においても付加的な長所を有している。コントラスト剤イメージングにおいては、被送信超音波ビームの圧力(または、他の音響出力パラメータ)を、コントラスト剤の望ましくないレスポンスを最小にする値に制限することが好ましいことが多い。これらの好ましい実施の形態では、ユーザは、コントラスト剤が存在しているか、または存在するであろう超音波イメージ内の領域(例えば、心臓の心室)を選択することができ、またその領域における被送信ビームの圧力を監視することができる。ユーザは、このフィードバックを用いて所望のレベルを達成するように超音波システムの動作パラメータを調整し、それによってコントラスト剤レスポンスを最適化し、非均一音響出力パラメータによってもたらされる効果を除去乃至は最小化することができる。ここに使用した“超音波システムの動作パラメータ”とは、広義に、音響出力パラメータに影響を与えるように調整できる何等かの動作パラメータのことである。動作パラメータは、限定するものではないが、アポジゼーション、送信開口(アパーチャ)内の要素数、焦点距離、送信電圧、及び超音波パルスの持続時間を含む。
【0017】
これらの好ましい実施の形態は、コントラストイメージデータの定量化を改善することに加えて、スタディ間のコントラストイメージデータの解析を援助することができ、決定された音響出力パラメータを記録することによってコントラストイメージング検査の再現性を改善する。更に、コントラスト剤を注入する前に音響出力パラメータを決定することによって、ユーザは、コントラストプロトコルを開始する前に所望の音響出力パラメータが達成されることを確認できる。これにより検査はより効率的になり、検査に必要とされるよりも多くのコントラスト剤を患者に注入することが回避される。薬品運搬容器を超音波によって破壊する薬品送給分野におけるような他のイメージング応用においても、同じような長所を得ることができる。例えば、好ましい実施の形態は、薬品送給及び投与のレートを決定することができる。
【0018】
ユーザが超音波システムの動作パラメータを調整することに加えて、またはその代替として、超音波システム自体は動作パラメータを自動的に調整することができる。例えば、ユーザは、コントラスト剤イメージングを最適にするために、1つまたはそれ以上の点のための目標音響出力パラメータ(圧力のような)を領域内に指定、またはプリセットすることができる。音響出力パラメータを決定した後に超音波システムは、システムの動作パラメータを自動的に調整することによって指定された音響出力パラメータを達成することができる(ステップ245)。
【0019】
これらの好ましい実施の形態には、幾つかの代替を使用することができる。一代替実施の形態では、ユーザが領域を超音波イメージ内に手動で選択する代わりに、またはそれに加えて、超音波イメージングシステムが領域を自動的に選択することができる。例えば、超音波イメージングシステムは、省略時に、またはユーザが指定した点、イメージ深さ、または方位送信焦点に基づいて領域を選択することができる。
【0020】
別の代替実施の形態は、超音波ビームの減衰に関する。超音波信号は伝播経路に沿って減衰するので、決定された音響出力パラメータは正確ではないことがあり得る。より正確な決定を得るために、超音波システム100は、被送信超音波ビームの音響減衰を決定し(ステップ235)、決定された音響減衰に関して決定済み音響出力パラメータを調整/較正して実際のイメージング状態を補償することが好ましい(ステップ235)。送信された超音波ビームの音響減衰を決定するために、トランスデューサ105からユーザが選択した領域までの音響線に沿って1回または複数回発射してこの音響線に沿うデータを取得し、この音響線に沿う減衰係数の推定を決定することができる。この線に沿って発射されるパルス毎に動作パラメータ(周波数及び帯域幅のような)を変化させることによって、この線に沿う平均減衰を推定することができる。減衰を計算する1つの適当な技術については、例えば、1981年のProc. IEEE Symposium on Sonics and Ultrasonics, #81CH1689-9, pp. 607-611に所載のH.E. Melton, Jr. 及び D.J. Skorton の論文“Rational-Gain-Compensation for Attenuation in Cardiac Imaging”を参照されたい。この技術は、受信され、ビーム形成された信号のエコー輝度を解析することによって、組織または血液の何れかを表すものとして超音波イメージの領域を識別する。組織は高強度のエコーを有し、血液は実質的にエコーを有していない。減衰は、血液及び組織のための典型的なパラメータを使用して計算される。
【0021】
上述した音響減衰較正とは別に、またはそれに加えて使用することができる別の較正技術は、音響出力パラメータの生体内測定に関する。計算された音響出力パラメータは、実際の測定、及び均質(理想的)イメージング状態であるとする特定の超音波システム動作状態のパワー管理モデルに基づくことができる。臨床状態においては不均質性が導入されることが多く、これらの不均質性のために実際の音響出力パラメータと、決定された音響出力パラメータとが大幅に異なるようになる。音響出力パラメータの生体内測定は、音響出力パラメータの関数として異なるレベルの非線形レスポンス(例えば、破壊)を有する1つまたはそれ以上の薬剤からなるコントラスト剤の集合(ポピュレーション)を使用して行うことができる。音響出力パラメータを増加させるように動作パラメータを変化させると、非線形レスポンスの最初の始まりが、関心点(手動で、または自動的に選択された領域)における最初の非線形しきい値が達成される時点を指示する。その後の非線形レスポンスは、集合内に、変化する非線形しきい値を有する付加的なコントラスト剤を用いて観測することができる。類似状態下のイメージング場内で達成される音響出力レベルを予測するのに使用される音響モデルを較正するために、単一の、または複数のしきい値レベルを使用することができる。
【0022】
別の好ましい実施の形態においては、ある超音波イメージ内のある領域において被送信超音波ビームの音響出力パラメータが決定される。次いで、決定された音響出力パラメータの指示が、その領域が超音波イメージ内の何処に位置するかの指示と共に供給される。超音波システムは、例えば機械的インデックスを決定して表示するのに加えて、表示された超音波イメージ上に視覚標識(例えば、ドット)をも表示し、表示された機械的インデックスに対応付けられたイメージ内の位置を指示することができる。別の例として、イメージ上に視覚標識を表示する代わりに、またはそれに加えて、機械的インデックスに対応付けられた範囲を表示することができる。この好ましい実施の形態を用いると、ユーザは、超音波イメージ内の機械的インデックス(または、他の音響出力パラメータ)の空間的位置を知る。もし機械的インデックスがユーザの関心位置に(例えば、コントラスト剤の位置に)出現すれば、ユーザはさらなる動作を要求しないであろう。しかしながら、もし機械的インデックスの位置がユーザの関心位置になければ、上述した好ましい実施の形態を使用して別の位置における音響出力パラメータを決定することができる。
【0023】
上述した好ましい実施の形態の何れかの種々の面の何れかは、単独で、または組合わせて使用できることに注目することが重要である。例えば、図2のフローチャート200には示されていないが、ステップ230、235、240、及び245は、一緒に、または別々に遂行することも、または結局は遂行しないこともできる。別の例として、超音波システムは、システムの動作パラメータを自動的に調整することによって、自動的に選択された領域(ユーザが選択した領域の代わりに、またはそれに加えて)内で指定された音響出力パラメータを達成するように動作パラメータを自動的に調整することができる。
【0024】
更に、別の代替実施の形態では、ピーク音響出力パラメータ(例えば、機械的インデックス)の位置に対応していない領域を超音波イメージ内に選択する。この領域は、超音波システムによって自動的に選択することも(例えば、その領域が、指定されたイメージ深さに、または送信された超音波ビームの送信焦点に自動的に選択される時のように)、またはユーザが選択することもできる。次いで、超音波システムは、その領域のための音響出力パラメータの指示を決定し、供給する。上述したように、これらの好ましい実施の形態の種々の面の何れかを、この代替実施の形態と共に使用することができる。
【0025】
以上の詳細な説明は、本発明が取り得る、そして本発明を限定するものではない選択された形状を例示することを意図している。本発明の範囲は、特許請求の範囲によってのみ限定されることを意図するものである。
【0026】
以下は、送信された超音波ビームの種々の音響及び熱的パラメータを決定する好ましい方法である。
【0027】
PMSモデルパラメータの決定:
(以前は:“pms.calc”計算された電力管理パラメータ)
(Stefan Schmitz:96/01/02)
(改訂:96/01/02 新Cファクタ対Vppアルゴリズム)
(改訂:97/10/15 帯域幅による音響変換重み)
(移動:98/04/07 to/people/pmsw/dev/spcalc; put under SCCS wa
s in interleaf; made sofutlink there:
System5→Acuson→adia→Scanning-specs→
Other-Features→mps calc)
【0028】
このドキュメントは、最初にパワー測定(“spmeas”)の、次いでパワー計算(“spcalc”)の結果のパラメータを記述する。これらはイメージングスペック内に入力される。これらにより、マシンは、ODSモデルの種々の音響及び熱的パラメータのための値をオンラインで計算することができる。
(注:表面における温度上昇のためのパラメータΔTは含まれない)
5つまでの独立パラメータにそれぞれ従属する9つのモデルパラメータが存在する。以下の2つの表は、
1.独立パラメータ(実際には、それらの超集合)
2.9つのモデルパラメータ
を要約している。
【0029】
標準単一線測定パラメータ:
all=(frq、ang、tx、loc、apo、fnu、cyc、vol)
frq=出力搬送波周波数
ang=線角度
tx =方位送信焦点
loc=測定位置(距離)
apo=アポジゼーション型
fnu=f数
cyc=サイクル内のパルス長(g数)
vol=送信p−p電圧
【0030】
IS及びそれらの従属になる計算されたパワー管理パラメータの表
【0031】
1+2)音響変換: a avg ( freq, cyc )及び a max ( freq, cyc ):
a avg(freq, cyc)[1]=E pa ac[W] / E pa el[W] または
W pa ac[W] / W pa el[W] (何れか選択可能)
の測定(ang、freq、tx、loc、apod、f#、cycについて)
の平均(ang、tx、loc、apod、f#、voltの)
=パルス平均音響変換
a max(freq, cyc)[1]=W min ac[W] / W min el[W]
の測定(ang、freq、tx、loc、apod、f#、cycについて)
の平均(ang、tx、loc、apod、f#、voltの)
E pa ac[μJ] =測定されたパルスエネルギ[μJ]
=0.01 * PII [μJ/cm^2] * x eebw[mm] * y eebw[mm]
x eebw[mm]、y eebw[mm]は、エネルギ等価ビーム幅である。
W pa ac[W]=E pa ac[μJ] / meas-τ[μs]
meas-τ =測定されたパルス長(立上がり時間アルゴリズム)
E pa el[W] =W pa el[W] * elec-τ[μs]
E pa el[W] =Sequoiaによって計算される電気パルス平均電力
elec-τ =Sequoiaによって計算されたパルス長(立上がり時間アルゴリズム)
W min ac[W]=測定されたパルス最大電力Pr
W min ac =(Pr * Pr / ρ * c)* x eebw * y eebw
最大Prレアフラクション圧力:Mpa
水の密度ρ =1g/cm3
水中の音速c=1.5246 mm/μs
W min acは測定によって計算されるので、それは計算のための入力内にある。
測定量は組織減衰されていないことに注目されたい。もし実際の測定が減衰されていれば、対応する減衰されていない値を先ず計算しなければならない。
【0032】
Acq からの要望:
W pa el[W]:電気パルス平均電力
τ el:電気パルス持続時間
E pa el=W pa el *τ el
W min el[W]:電気的最小ピーク電力(負へ移行するピーク)
【0033】
処理
→最初に個別の変調周波数及びサイクルカウント毎に両音響変換が計算される。
平均パルス及び最大パルスのための分離した値が存在する。
→所望の出力サンプルグリッドが選択される。
1.出力搬送波周波数の場合、グリッドは個別の変調周波数(usm mod freq)
に等しい。
2.サイクルカウントの場合、グリッドは個別のサイクルカウント(usm cyc
)に等しい。
【0034】
→各グリッド点(サンプルグリッド周波数及びサイクルの別々の対)毎に、全て
の(拒絶されなかった)測定の重みつき合計平均を以下のように測定する。
a avg/max(gridFreq i, gridCycle j)
=Σmeasurementsa avg/max * weight/Σmeasurementsweight
ここに、weight=exp(−ac weight exp fact * distance2)
但し、distanceは以下の2つの選択の一方によって与えられる。
もしspcaleオプションac weight exp bandwidthが0ででなければ、
distance2=
ac weight exp fact freq * [(gridFreq i−measFreq)/avg-measFreq]2
+ac weight exp fact cycl *[(1/gridCycle j)−(1/measCycl e j)]2
それ以外ならば(ac weight exp bandwidthが0である時;古いアプローチ)
distance2=
ac weight exp fact freq * [(gridFreq i−measFreq)/1 MHz]2
+ac weight exp fact cycl *[(gridCycle j−measCycle j)/1 cycle]2ここに、ac weight by bandwidth、ac weight exp fact、ac weight exp fact freq、及びac weight exp fact cyclはspcalc入力オプション(数)である。
avg-measFreqは、(拒絶されなかった)measFreqの平均である。
【0035】
注;重みつき合計は、グリッド再サンプリング、特に縁における“グリッド孔
”(これら孔を充填する一定の補外及び次元に関する固定された次数を原
因として悪い結果をもたらす)の問題を回避する。
より新しい重み付けアプローチ(オプションac weight by bandwidthが0
でない時の帯域幅による)は、長いパルス長を、測定によって指示された
ように遙かに近く移動させる。これは1997年に導入された。
【0036】
オプション(テンプレートファイル: spcalc.options 内のコメント参照):
→測定τ/電気的τの比が所望範囲外にある測定を拒絶する。この範囲の低端及
び高端は入力オプションである。(省略時0及び999、即ちターンオフ)
→a avgは、パルス平均電力またはパルスエネルギを使用して計算することがで
きる。但し、エネルギ=パルス平均電力 * パルス幅である。ユーザは、どち
らを使用するかを選択することができる。
→a avg及びa maxは、それらの間の平均として計算することができ、その場合そ
れらは同一である(及び、平均に等しい)。後方互換性の場合。
→a avg及びa maxは、異なるサイクルで同一値を有する周波数だけの関数として
計算することができる。後方互換性の場合。
→ac weight by bandwidth、ac weight exp fact、ac weight exp fact freq、及びac weight exp fact cycl(上記参照)。
【0037】
2)迎角焦点:
z el(freq, angle)[mm]=最小迎角帯域幅位置の測定(周波数、角度の)
処理:
→迎角焦点位置が法定距離限度のマージン内にある測定だけを受入れる。
→出力グリッド周波数及び角度を、個別の変調周波数及び個別のusline角度に等
しくセット。
→(個別の、及び格納された)出力搬送波周波数及び角度のデータをビンに入れ
る。周波数ビンが(必ずしも)出力グリッド周波数に一致しないことに注目さ
れたい。
→各ビン毎の平均を計算。
→各出力グリッド角度毎に、補間及び一定補外(これは周波数内の孔を充填する
)によって周波数を出力グリッド周波数に再サンプル。
→角度内のどの残余の孔をも充填(補間または一定補外によって)。
【0038】
3+4 ) 焦点における迎角ビーム幅 (L0y) 及び実効開口ファクタ ( eファクタ ) :
L0y(freq, angle)[mm]=meas(freq, angle, tx, loc について)からの計算さ
れた焦点における迎角帯域幅の平均(tx,locの)
e-factor(freq, angle)
[1]=meas(freq, angle, tx, loc について)
からの計算された迎角実効開口ファクタの平均
(tx,locの)
z1 =meas locの距離
zfy =迎角焦点
Ly =距離(測定された)におけるee帯域幅
LSy =表面(計算された)におけるee帯域幅=a0y * e-factor
L0y =焦点(計算された)におけるee帯域幅
a0y =迎角開口
e-factor =表面ファクタにおける迎角帯域幅=LSy / a0y
u=[1−{z1[mm]/zfy[mm}]とすれば、
Ly(z1)[mm]=L0y[mm] * [1+{( LSy[mm] / L0y[mm] )2−1} * u2 ]1/2
=[ L0y[mm]2 * ( LSy[mm]2 − L0y[mm]2 ) * u2 ]1/2
L0yについて解けば:
L0y[mm]=[( Ly(z1)[mm] 2 − LSy[mm]2 ]* u2 )/(1− u2)]1/2
L0yのための相対平均合計重みは:
w=[ L0y * ( 1− u2 ) / Ly(z) 2 ] 2
eファクタ=LSy/a0yについて解けば:
e-factor=(1/a0y ) * [{( Ly(z1)[mm] 2 * L0y[mm] 2) /u2 }
+L0y[mm] 2 ]1/2
eファクタのための相対平均合計重みは
w=[ e-factor * a0y2 * u2/ Ly(z) 2 ] 2
【0039】
注:
L0y を計算する場合:
迎角焦点に関しては、測定距離がソフトリミットに近過ぎる測定を無視する。
Ly>abs{( z−zfy )/zfy)}* Lsy( z<2*zfyの時)が正根より下の数である場合に限る。もしそうでなければ、ビーム幅が狭過ぎて公式に当て嵌めることはできない。我々は単にこれらの場合を無視するであろう。もし本当に焦点に接近して測定すれば、これは発生しないであろう。
上記公式を使用してL0yがLSyより大きいことが分かった場合には、その測定を無視する。しかしながら、モデル計算では、これらの場合(焦点におけるビーム幅が表面におけるよりも大きい場合)のために線形公式を使用する。その公式は、
L0y[mm]=LSy[mm]+(Ly(z1) [mm]−LSy[mm])*(zfy[mm] / z1[mm])
である。
送信及び測定位置を平均する場合、重みつき合計平均を遂行する。但し、相対重みは推定される誤差の逆平方に比例する。
L0y内の誤差={ d-L0y/d-Ly(z) } * const * Ly(z)
=const * Ly(z)2/L0y * ( 1−u2 )
相対重み:
w={ L0y * ( 1−u2 )/Ly(z)2 }2
注:たとえ何を測定したとしても、関数はLy(0または2*zfy)=LSyを通過し
なければならないので、重みはz1=0またはz1=2*zfyの場合に0になる。
【0040】
eファクタを計算する場合:
迎角焦点に関しては、測定範囲がソフトリミットに近過ぎる測定を無視する。
eファクタを計算する場合、それは複素数または無限大になり得る(もしu2=Ly(z1) 2<(1−u2) * L0y2ならば)。その場合には、その測定は無視する。
送信及び測定位置を平均する場合、重みつき合計平均を遂行する。但し、相対重みは推定される誤差の自乗に逆比例する。
eファクタ内の誤差={ d-e-factor/d-Ly(z) } * const * Ly(z)
=const * { Ly(z)2/e-factor * a0y2 * u2 }
eファクタのための相対重み:
w={ e-factor * a0y2 * u2/Ly(z)2 }2
【0041】
処理:
1.eファクタの開始値を使用してL0yを計算する。この処理は、(直接的な平均の代わりに)eファクタ及び重みの加算のために上記公式を用いて重みつき合計平均を使用することを除き、迎角焦点の場合と同じである。
2.次いで、最終出力グリッドへサンプルされた結果的なL0yを使用してeファクタを計算する。この処理は、1.の場合と同じである。
3.次いで、1.及び2.を繰り返す。但し、繰り返しカウントは入力オプションである。
【0042】
オプション(テンプレートファイル: spcalc.options 内のコメント参照):
→eファクタ開始(省略時0.666667)
→繰り返し数(省略時1)
(各ステップを1繰り返しと数える、従ってもしn=1ならばL0yを適合させるだけであり、もしn=2ならばL0y及びeファクタを1回適合させ、もしn=3ならばbファクタを2回適合させるが、eファクタは1回だけ適合させる等々である。)
− − − − − − − − − − − − − − − − − −
【0043】
6+7)焦点における方位ビーム幅及び実効開口ファクタ(b及びeファクタ)b (freq, angle, apod){1}=bファクタ及びeファクタのmeas(angle,freq,
tx,loc,apod,f#について)の平均(tx,loc,f-numの)
これは、迎角処理に極めて類似している! 従って “y”を“x”に置換し、最早焦点zfxを使用せずに名目的な送信焦点であるとすること以外は、上記公式を適用する。またアポジゼーション型内に付加的な従属性が存在する。これは、ビンニング、再サンプリング、出力グリッド、及びロジスティックスだけにそのような影響を与える。
唯一の“実際の”差は、それが迎角のためのものであって、bファクタのものではないので、モデルパラメータが“焦点におけるビーム幅”(Lox)ではないことである。
【0044】
L0x =bファクタ * λ * f数
λ =波長=周波数/速度
速度 =1.5246mm/μs;水中の音速
f数 =焦点/開口
a0x =開口
zfx =焦点
従って、
b=L0x[mm] * LSx[mm] * freq{MHz}/speed{mm / μs} *zfx[mm]
【0045】
注:前記spcalcバージョンでは実効f#を使用したが、ここでは上記公式内に真のf#を使用する!
bファクタのための平均・加算・重みはL0yの場合とは僅かに、丁度スケーリングファクタだけ異なる。ここでの“ストーリー”は、次のようである。
焦点から遠くの測定の場合(特に、焦点の2倍付近の深さにおける場合)、測定の小さな変動がbファクタに極めて大きい変動をもたらし得る。
従って、合計内の各項に、bの予測される測定誤差の逆平方だけ重み付けする。
【0046】
bの予測される測定誤差は:
const * derivative of b-factor w.r.t Lx(z) * error in Lx(z)
また、Lx(z)内の誤差はLx(z)に比例するものとして
d-b/d-Lx=b *(Lx/L02) * (1/1−u2) 但し、u=1−(z1/zfx)、
予測される誤差=const * b * (Lx2/L02) * (1/1−u2)
相対重み:
w={( L02/Lx2 ) * ( 1−u2 )}2
注:たとえ何を測定したとしても、関数はLx(0または2 * zfx)=LSxを通過しなければならないので、重みはz1=0またはz1=2 *zfxの場合に0になる。
【0047】
処理:
→出力搬送波周波数、角度、アポジゼーションの型の(格納された)個別の値の各組合わせ毎に、初期(または、最新の)eファクタを使用してbファクタの重み付き合計平均を計算する。
注:Lx(z)>u * LSxである場合に限ってbファクタは実の解を有している。
それ以外は、測定条件を無視する。
→結果を最終出力グリッドへ再サンプルする(最初に周波数を、次いで角度を再サンプルするが、アポジゼーション型は再サンプルしない)。
注:アポジゼーション型次元に関しては再サンプル“しない”。
しかしながら、我々は再サンプルしている(他の2つの次元を充填する)ので、何等かの孤立した孔はそれ以上存在し得ない。唯一の孔は、固定されたアポジゼーションの型に関して意味を有する完全な“平面”であり得、周波数または角度に関しては良好な値は存在しない。この場合、充填することは確かに誤りであろう(実際、これはむしろ病的である)。
【0048】
→今度は、eファクタのための上記ステップを繰り返す(それのための、及びその重みを加え合わせるための公式は、L0x=b-factor * λ *f数とする以外は迎角のための公式と同じである)。bファクタ
(最終出力グリッドへ再サンプルされ、次いで所望の値に戻すように補間されている)
を先に計算したように使用する。
→bファクタに関して(新しいeファクタを与えて)、
次いでeファクタに関して(新しいbファクタを与えて)繰り返す。
【0049】
オプション(テンプレートファイル: spcalc.options 内のコメント参照):
これらは、迎角の場合と同一である:
→eファクタ開始(省略時0.666667)
→繰り返し数(省略時1)
− − − − − − − − − − − − − − − − − −
【0050】
8) C ii : Iappa / Ispta /Cファクタ:
C ii ( freq, ang, tx, apod, volt){1}
=C iiの測定(全てについて)の最大(位置、f数、サイクル、sppa/spt
aの)
C iiには、3つの選択:C isppa、C ispta、及びC piiがある。1つだけが出力内にある(名目的に、C isppaと呼ぶ)。spcalc入力オプションによって、各測定毎に、これらの3つ、対としての最大、またはこれら3つの総合最大の何れかを選択することができる。標準及び省略時は以下の通りである:
C isppa final=max (C ispta, C pii )
C isppa, C ispta, C pii:
C isppa=Isppa meas[W/cm2]/Isppa calc[W/cm2]
及び
C ispta=Ispta meas[W/cm2]/Ispta calc[W/cm2] RES→∞の場合
但し、RES→∞はres/panボックスがビーム幅よりも遙かに大きい
限界内にあることを意味する。
及び
C pii=PII meas[μJ/cm2]/PII calc[μJ /cm2]
【0051】
ここに、
Isppa calc[W/cm2]=max overz{100 * W ac pa[W] * 10-0.003 * freq[MHZ] * z[mm]
/A ee beam(z)[mm2]
PII calc[μJ /cm2]=Isppa calc[W/cm2] * el τ[μs]
実際には、“max overz”とは、z=0、z=z elev focus、及びz=z focusにおいて評価し、これらの最大を採ることを意味する。
注:W ac pa[W] =a(freq)[1] * W el pa[W]
A ee beam[mm2]=Lx(z)[mm] * Ly(z)[mm]
【0052】
C ispta :
resボックスサイズが未知であるか、またはresボックス、線間隔あけ組合わせの全てについて上述した計算をすることは望ましくないから、RES→∞の限界が使用される。また、限界RES→0はPIIまたはIsppaケースによって処理されるから、以下の大きいRESボックス限界が有利である。
RES→∞の場合:
Ispta meas=(PII(z) * Lx meas(z) * NF/FT * RES) * 10-0.003 * f meas * z 注:Lx * NF/RESはBOF(または、BOF * FSC)である。
RES→∞の場合:
Ispta calc=(a * E pa el * NF/RES * Ly mod(z) * FT ) * 10-0.003 * f mod * z
但し、Lx,y:方位または迎角内のビーム幅(測定された、またはモデル化され
たの何れか)
NF :走査中の発射数
FT :フレーム時
a :モデルの音響変換
E pa el[W]=Sequoiaによって計算された電気パルスエネルギ
=パルス電力×τ
【0053】
C isptaは、オプションとして次の2つの方法で計算することができる。
1)この比は全ての測定zについて評価され、その最大が採られる。
C ispta=(PII(z) * Lx meas(z) Ly mod(z) / a * E pa el)
*(10-0.003 * f meas * z / 10-0.003 * f mod * z)
zのこの値は、全ての量、測定、またはモデルにおける測定範囲である。
要するに:
C ispta=meas-term(z)/calc-term(z)の測定されたzについての最大で
ある。
2)zは、それがC isppa及びC piiのためのものであるのと同様に処理される
。meas-termは測定範囲が“最大化”され、calc-termはz=0及び焦点が
“最大化”されている。
C ispta=meas-term(z-meas)の測定されたzの最大/
calc-term(z-calc)のz=0及び焦点の最大
この第2の方法は、C isppa及びC piiに類似している。
【0054】
*このセクションは置換された:
*そこで、C(V)をVpp2の関数として適合させ、以下を求める。
* C 0 ( freq, ang, tx, apod )[1]
* V 0 ( freq, ang, tx, apod )[V]
* p+ ( freq, ang, tx, apod )[1]
* p− ( freq, ang, tx, apod )[1]
*ここに、
* C(V)=C 0[ 1−P+{( V / V 0 )−1}2 ] V>V 0の場合
* C(V)=C 0[ 1−P−{( V / V 0 )−1}2 ] V<V 0の場合
*C 0及びV 0を見出す:
* C 0=V[ C(V) ]についての最大
* V 0=C(V)を最大にするV、またはC(V 0)=C 0
*P+/−を見出す:
* P+=V>V 0[( V 02 / C 0) * { C 0−C / ( V− V 0 ) 2}]についての最小
* P−=V<V 0[( V 02 / C 0) * { C 0−C / ( V− V 0 ) 2}]についての最小
* P+/−は、もしV>V 0またはV<V 0をそれぞれ有するC(V)が存在しなければ、
*0にセットされよう。
*置換されたセクションの終わり
【0055】
****より多く削除されたスタッフ
*電圧グリッドのためのCファクタ適合:
* 周波数、角度、送信の固定された各セット毎に、固定された電圧グリッド値V
* nのためのCnのセットを以下のように計算する(但し、n=1‥N)。
* 1)C * V2はV2内で区分的に線形である。ここに、区分はこれらのVnによ
* って区切られる。
* 2)全ての測定された点C * V2は、区分的には線形である曲線上に、または
* その下にある(決して曲線の上にはない)。
* 3)C * V2は、原点を通る、即ちV=0の場合には0であるとする。
【0056】
*アルゴリズム:
* 1) 先ずV間隔:0<V<V1で開始し、測定点がその曲線上に存在しないよ
* うなその曲線が取り得る最小勾配を見出す。
* これは、第1グリッド点V1における値C * V2をもたらす。
* オプショナルなグリッド重畳ファクタを有する「次の」間隔の部分を含ま
* せるために使用される間隔を拡張することができる。従って、より大きい
* 間隔内に存在する全ての測定点が検討される。それにも拘わらずこの結果
* は現在のグリッドのエンドポイントのために計算され、この手順が繰り返
* される。従って、グリッド重畳を用いずに、次のものを使用する。
* i番目の間隔の場合:
* V ( i−1 ) 2 < V < V i2
* グリッド重畳ファクタ(GOF)を用いる:
* V ( i−1 ) 2 < V < V i2+GOF * [ V ( i+1 ) 2− V i2 ]
* 再度、V iのみにおけるC*V2を計算するために、それでも、現在の間隔
* のために得られた最大勾配を使用する(グリッド重畳ファクタを用いて、
* または用いずに)。
* 2)全てのグリッド電圧間隔が遂行されるまで、次に高い電圧間隔について
* 上記述ステップを繰り返す。
*********削除されたスタッフの終わり
【0057】
C ファクタ電圧グリッド計算及びグリッドへの適合:
これは、1996年9月26日に変更されたアルゴリズムである(古いバージョンに関しては上述したので参照されたい)。
1.周波数、角度、送信焦点、及びアポジゼーション型の各独特のセット毎にVppの関数として全ての“生”Cファクタを計算し、このセットを“ケース”と呼ぶ。次いでxの関数としてyを計算する(yはC-factor *Vpp2であり、xはVppである。最終適合はy対xの区分的線形曲線であろう。
2.各ケース毎に、x-y対及び点(0,0)の部分集合からなるその凸のハルを計算する。
ハルは、先ず点(0,0)をそれへ加算することによって計算する。
次いで、最大勾配を有するx-y対を見出し、それをハルへ加算する。
次に、x>x0及びy>y0を有する生x-y対を見出す
(x0−y0は最後のハル点であり、これも最大x及びy値を有する点である)。
これはx0−y0から計算された最大勾配を有している。即ち全てのxi>x0及びyi>y0について、
最大勾配(xi,yi)=(yi−y0)/(xi−x0)
換言すれば、次のハル点は(xi,yi)であり(但し、xi>x0及びyi>y0である)、勾配(xi,yi)はxj>x0及びyj>y0を有する全ての点の最大である。
xi>x0及びyi>y0を有する生の点をそれ以上見出すことができなくなるまで、繰り返してハル点を加算する。
【0058】
3.初期xグリッド値を、将に2つの値、即ち全てのハルのケースについて0及び最大xにセットする。
4.各ハルのケースを、以下のようにしてグリッドに適合させる。点(0,0)を適合に加算する。
次いで、各xグリッド値毎に、先行xグリッド値からの最大勾配w.r.tを計算し、その先行グリッド点(x,y)からの勾配を有する線と、現在のグリッド線(xグリッドを通る垂直線)との交わりを計算する。これが、新しいグリッド点である。
次いで、全てのグリッド点がカバーされるまで繰り返す。換言すれば、次のyグリッド=最後のyグリッド+最大勾配*( 次のxグリッド−最後のxグリッド)但し、全てのxi>最後のxグリッド、及びyi>最後のyグリッドについて、最大勾配=(yi−最後のyグリッド)/(xi−最後のxグリッド)の最大である。
もし、最大勾配が0または負であるか、またはもしx>最後のxグリッド点を有するハル点が存在しなければ、最後のハル点から、以下のように補外する。
y(xグリッド)=(最大yハル/最大xハル)*xグリッド
注:最大xを有するハル点は、最大yをも有している。
【0059】
5.全てのケースに関して、各ヌル点における、及び各グリッド点におけるハルと適合との間の誤差を計算する。補外は上述した
(“一定Cファクタ補外”)通りである。
もし最大絶対誤差が入力許容差よりも大きければ、別のグリッド点を加える。そのグリッド点は、最大誤差を有する点
(1つのケースに関して、ハル上の点、またはグリッド上の点の何れか)を見出すことによって選択される。そのケースに関して、次に低いハル点を見出し
(それが既にグリッド内にない限り、次に低いものを選択する)、その点をグリッドへ追加する。
6.最大絶対誤差が許容差より小さいか等しくなるまで、ステップ4及び5を繰り返す。
【0060】
7.次に、再度グリッド点を除去することを繰り返し試みる。全てのケースに関して各グリッド点自体を除去し、最大絶対誤差を計算する。もしこの誤差が許容差より小さいか等しく留まっていれば、その点を除去する。さらなる縮小を達成できるまで、低から高x値へ、次いでその逆にグリッドを通過する。
8.全てのケースに関して、適合Cファクタは適合曲線のy/x比である。Vppグリッドは、xグリッド値の平方根からなる。
注:上述した“全てのケース”とは、Ispta及びMIの両者(及び、全ての周波数、角度、焦点、及びアポジゼーション型)に関していることを意味している。
【0061】
処理:
−変調周波数、角度、方位送信焦点、及びアポジゼーション型の個別の各値毎に、測定された、及びモデルパラメータ(組織減衰された)を計算し、比を形成する。これが、電圧の関数としてのCファクタである。
Cファクタ、電圧(Vpp)、及び測定の出力搬送波周波数を記録する。ここで、τ拒絶基準を適用し、測定τ/電気的τの比がオプション入力範囲外にある測定を無視する(音響変換の場合のように)。
−各ケース毎に、上述したようにグリッド及び適合Cファクタを計算する。
−3つの各独立パラメータ(角度、焦点、アポジゼーション型)毎に、(平均された)出力搬送波周波数によってCn値を分類し、そのリストを(個別の;分類された)変調周波数に再サンプルする。
−これにより、焦点及び角度を再サンプルすることによって(但し、アポジゼーション型は再サンプルされない。上述した方位ビーム幅ファクタに関する処理の注を参照されたい)
孔が充填されている規則的な矩形の4つの次元の表が得られる。
【0062】
オプション(テンプレートファイル: spcalc.options 内のコメント参照):
−低及び高τ拒絶比(音響変換参照)
−グリッド計算及び適合のための最大絶対誤差(%)
−C isptaのための範囲処理(前述した:測定範囲における比を計算して最大を採るか、「または」、測定範囲のmeas-termの最大を採ってz=0及び焦点のcalc-termの最大を採り、比を形成するの何れかを参照)
−省略時値によって、spc val.datファイルの代わりに、spc vsub.datファイル
及びspc vcfa.datファイルを作る。
− − − − − − − − − − − − − − − − − −
【0063】
9) c mi : MI c-factor :
c mi( freq, ang, tx, apod)
[1]=( MI meas / MI calc )2のmeas(全てについての)の最大(loc,f#,cycの)
MI meas[1]=Pr(z)[MPa]/Mpa*sqrt(10-0.003*freq[MHz]*z[mm])/sqrt(freq[MHz]/MHz
MI calc=max overz[W min a[W]*ρ[g/cm3]*c[mm/μs]*10-0.003*freq*z/
A ee beam(z)[mm2]*freq[MHz]]1/2
実際には、“max overz”とは、z=0、z=z elev focus、及びz=z focusにおいて評価し、これらの最大を採ることを意味する。
注:Pr2=(W min ac/A ee beam)* ρ *c
ρ=水の密度 1 g/cm3
c=水中の音速 1.5246 mm/μs
W mim ac=a(freq) * W mim el
A ee beam=Lx(z) * Ly(z)
電圧グリッド値Vn(n=1‥N)のためのCn値の計算は、C isppaに関する計算と同様である。
【0064】
処理:
C isppaについてと同様である。
【0065】
オプション(テンプレートファイル: spcalc.options 内のコメント参照):
C isppaに適用されるオプションの3つがC miにも適用される。
−低及び高τ拒絶比(音響変換参照)
−出力グリッド電圧のセット
−電圧グリッド重畳ファクタ(C ispta参照)
−省略時値によって、spc val.datファイルの代わりに、spc vsub.datファイル
及びspc vcfa.datファイルを作る。
− − − − − − − − − − − − − − − − − −
【0066】
電気的コンダクタンスの再スケーリング:
usmマシンは、3つの項(平均パルス電力、及び最小及び最大パルス電力)に3つの項、即ち、電圧項、アポジゼーション合計項、及び電気的コンダクタンスを乗ずることによって電力項を計算する。この電気的コンダクタンスは、マシンへの電力管理によって与えられる表である。この表は、電気的コンダクタンス(単位:μモー)対出力搬送波周波数の対からなる。オンライン動作中に、それは、現在の出力搬送波周波数のための電気的コンダクタンス値を発生させるために補間される。
spcalcは、測定が完了した後に、電気的コンダクタンスを変化させる能力を与える。それは、電気的コンダクタンス対出力搬送波周波数の新しい表を採り、この新しいコンダクタンスに従って電力項を再スケールする。勿論、これはモデルパラメータ、即ち音響変換及びCファクタを変化させることになろう。
同じ電気的コンダクタンス表を、spcalcへの入力として使用されたマシン内へロードすることが重要である。
【0067】
ここに公式が存在する:
old conductance[μモー] 測定中にマシンから読出す
old W pa el[W] 測定中にマシンから読出す
old W min el[W] 測定中にマシンから読出す
new conductance[μモー] 新しい値:分離した入力表からspcalcへ補間
電力値(音響変換、Cファクタ)を使用する上式の全てにおいて、次のことを行う。
W pa el=old W pa el * new conductance/old conductance
W min el=old W min el * new conductance/old conductance
【0068】
電力測定中に必要な Acq によって計算されるパラメータの概要:
信号ライン送信パラメータ、Plus
a0x[mm]: 送信方位開口
W pa el[W] 電気パルス平均エネルギ
W min el[W] 電気的最小ピークエネルギ(負に移行するピーク)
V[V] 送信電圧
conductance[μモー]:
elec-τ[μs] :
【0069】
電力計算の出力( spcalc ):
出力は2回生成される。即ち、1回はZのためにフォーマットされ、1回はSのためにフォーマットされる。プログラムは、出力フォーマットを説明するテキストファイルをも発生する。ターンオンすることができるデバッグ出力ファイルも存在する。これらは、“説明ファイル”内にも記述される。
【0070】
Zのための出力:
注:これは、現在も未だにそのモデルの初期バージョンと後方互換的なフォーマットである。幾つかのオプション(spcalcオプション内の)を正しくセットしなければならない(オプション入力ファイルを参照されたい)。また、方位bファクタは、後方互換性のために1.5で除される。
出力はmglフォーマットであり、各PMSパラメータ毎の表と共に単一のファイルにされる。各表は、パラメータ値及び独立パラメータのためのグリッドサンプル値を記録している。
1)音響変換:2Dテーブル
ac-conv a(freq)
2)迎角焦点及び最小ビーム幅:3Dテーブル
z el L0y(freq, angle) freq angle
3)方位ビーム幅パラメータb:3Dテーブル
b(freq, angle) / 1.5 freq angle
4)Isppa C-factor
C freq angle tx apod volt
5)MI C-factor
C freq angle tx apod volt
【0071】
Sのための出力:
1)Sは、6つの全ての次元においてオンライン補間を行い、従ってイメージングスペックは、正確にpms計算手順を格納する。即ち、イメージングスペックが作られる時には、補間または再サンプリングは存在しない。
2)省略時には、3つの表を充填する3つのデータファイルを生成する。これらのファイル(表)は、spc grid.dat(pms meas grid)、spc vsub.dat(pms meas sub values)、及びspc vcfa.dat(pms meas c factor)である。後方互換性に関して、2つのデータファイル、即ち、pms meas grid表を充填するためのspc grid.dat及び絶対pms meas values表を充填するためのspc val.datを生成することができる。pms meas gridは6つの独立変数(周波数、角度、焦点、アポジゼーション型、サイクル、電圧)を含み、他の表は9つの従属変数を含む。
3)pms meas sub values表は、独立変数、即ち周波数、角度、アポジゼーション型、及びサイクルカウントに従属する。pms meas c factorは、周波数、角度、焦点、アポジゼーション型、及び電圧に従属する。データを2つの表に分割することによって、ISデータベースのためのストレージを大幅に節約することができる。pms meas values表は、全てのデータを単一の表内へ管理する。これは、spcalcオプション“make spc val.=1”を通して未だに支援されている。省略時には、2つの新しい表、即ちpms meas sub values及びpmsmeas c factorのためのデータが代わりに作られる。
【0072】
4)2つの表pms meas sub values及びpms meas c factorは、以下の構造によって収集ソフトウェア内に表される。struct PmsMeasSubValuesIsm{Float acoustPulseAvgConv;Float acoustPulseMaxConv;Float elevFocusMm;Float elevMinBwidthMm;Float elevEffAperFactor;Float azimBFactor;Float aximEffAperFactor;};andstruct PmsMeasCFactorIsm //PMS data C factor values{Float isppaCFactor;Float miCFactor;};( both from the AcqIsmDataStruct.h file )
【0073】
“9”のテーブルは、各パラメータ毎の列を有し、その6次元マトリックスは列ベクトルとして線形の態様で格納される。実際には、これは、これに似た構造のアレイを意味する。
struct PmsMeasValuesIsm
//PMS data values{Float acoustPulseAvgConv;Float acoustPulseMaxConv;Float elevFocusMm;Float elevMinBwidthMm;Float elevEffAperFactor;Float azimBFactor;Float aximEffAperFactor;Float isppaCFactor;Float miCFactor;};従って、例えば、もし3つの周波数、2つの角度、5つの焦点、4つのアポジゼーション、6つのサイクル、及び7つの電圧を“悉く”測定すれは、表は、3×2×5×4×6×7=5040行を有することになる。
【0074】
5)6つの各独立パラメータ毎のサンプルグリッド値も、単一の表内にパックされる。それは、各パラメータ毎に1つずつの計6つの列を有している。
行の数は、4つのパラメータの何れかのための値の最大数である。
従って値リストが行数よりも短いパラメータの場合には、列ベクトルにターミネータ値を詰める。
struct PmsMeasGridIsm //PMS grid definition
{
Float outputCarrierFreqMhz;
Float usLineAngleDeg;
Float xmtFocusDepthMm;
Int baseApodId;
Float cycleCount;
Float xmtVppV;
};
従って、上例(3周波数、2角度、5焦点、4アポジゼーション、6サイクル、7電圧)においては、この表は7行を有することになる。
周波数、角度、アポジゼーション、及びサイクルに関しては、少なくとも4、5、2、3、及び1行がターミネータ値(多分、−1000000)を有するようになる。
【0075】
注:
1)spmeasは角度を測定することを試みる。
もし、迎角焦点を見出すことができなければ(それが浅過ぎるために)、プログラムはそのように“言う”が、進行はしない。
もし特定の送信焦点のピークを見出すことができなければ、プログラムは次の送信焦点へ移る。極めて大きい角度の場合には、
多分送信焦点はデータを有していないであろう。これらの孔は、出力グリッドへ再サンプルすることによって充填されよう(既述)。
【0076】
出力搬送波周波数のジッターの取扱いの要約説明:
これは、種々の“処理”セクションにおいて説明したものの要約である。
1)個別の変調周波数のリストを見出す。これらの値は、出力(サンプル)グリッド出力搬送波周波数としても役立つ。
2)個別の角度、送信焦点、アポジゼーション型、及びサイクルのリストを見出す。
3)出力搬送波周波数及びサイクルの個別の各出力グリッド値対毎に;
全ての音響変換に関して重み付き合計平均を計算する。
出力グリッド点に“より近い”測定は、“より遠い”測定よりも大きい重みを有している。
それでも、全ての測定は単一のグリッド点として平均される。
4)個別の出力搬送波周波数及び個別の角度毎に、
“z el”及び“L0y”及びeファクタf ocを計算して平均し、
これら3つを出力グリッド周波数値のリストへ戻して補間する。
5)個別の出力搬送波周波数及び個別の角度及びアポジゼーション型毎に、方位bファクタ及びeファクタ及びf ocを計算して平均し、これら3つを出力グリッド周波数値のリストへ戻して補間する。
6)個別の変調周波数、角度、送信焦点、アポジゼーション毎に、C(V)を計算し、f ocを平均する。
7)C(V)から、Cn(Vn)を決定する。
8)それらを出力搬送波周波数の出力リストへ戻して補間する。
【0077】
オンラインPMSモデル計算:
(以前は:“pms.online”音響及び熱的パラメータのオンライン計算)
(Stefan Schmitz:94/01/31)
(改訂:97/05/06)
(改訂:97/10/16 rel 2.5のためのz=0についてIspta変化)
(移動:98/04/07 to/people/pmsw/dev/spcalc; put under SCCS was in
interleaf; made sofutlink there:
System5→Acuson→adia→Scanning-specs→
Other-Features→mps calc)
【0078】
このドキュメントは、Sequoiaが動作している時のSequoiaによる音響及び熱的パラメータの計算を記述する。これらのパラメータは、FDA及び内部規則のためのものである。幾つかのパラメータは調整され、他のパラメータは単に表示されるだけである。
2つのグループ、即ち、単一の発射に適用される短時間一定パラメータ(STC)、及び長時間(フルフレーム)にわたって平均される長時間一定パラメータ(LTC)が存在する。
【0079】
STC:
1) Isppa t
2) MI
LTC:
4) Ispta t
5) TIS, TIB, TIC
6) dT skin, dT air
IEC1157のためのパラメータ
ヘッドルーム
注:次元は[ ]で示す。例えば、[mW/cm2]
[1]は無次元であることを意味する。
− − − − − − − − − − − − − − − − − −
【0080】
1)Isppa t[W/cm 2 ]
Isppa t[W/cm2]=100 * max overz[ c sppa[1] * ( W.3 pa acoust[W]/
A ee beam[mm2] ) ]
A ee beam(z, freq, ang, zfx, a0x, apod)[mm2]=Lx[mm] * Ly[mm]
Lx(z, freq, ang, zfx, a0x, apod)[mm]
=L0x[mm] * [1+{(LSx[mm]/L0x[mm]) 2−1} * {1−(z[mm]/zfx[mm])}2]1/2
LSx=a0x * azimBwidthAtSurface-fraction
L0x[mm]=(cH20[mm/μs]/freq[MHz]) * (zfx[mm]/a0x[mm]) (=bλf数)
Ly(z, freq, ang)[mm]
=L0y[mm] * [1+{(LSy[mm]/L0y[mm]) 2−1} * {1−(z[mm]/zfy[mm])}2]1/2
LSy=a0y * elevBwidthAtSurface-fraction
c sppa[1] * V2=V2内に区分的線形としてモデル化され、イメージングスペック
からの幾つかの値対(Vn、Cn)(n=1‥N)によって記述される。
これは、Cが
C(V) * V2=α * Cn * Vn2 + (1−α) * Cn+1 * Vn+12
のように補間されることを意味する。
ここに、Vn < V ≦ Vn+1
α=(Vn+12− V2)/(Vn+12− Vn2)
または
C(V) =β * Cn + (1−β) * Cn+1
β=(V〜 2−Vn2)/(Vn+12− Vn2)
ここに、V〜=Vn+1 * Vn/V
W.3 pa acoust[W]
=aAvg(freq,cyc)[1] * W0 pa elec[W] * 10-0.003 * freq[MHz] * z[mm]
W0 pa elec[W]
=10-6 * g(freq)[μモー] * Asq unsc sum(apod,f#)[1] * V pulse rms2 [V2]
【0081】
【0082】
2)MI[1] :
機械的インデックス
【0083】
3)Ispta t(Sequoiaソフトウェアバージョン1.61で開始、2.5に変化)
Ispta t[mW/cm2]=100 * max overz[ c sppa[1] * ( W.3 rms acoust[mW]
/A ispta(z)[mm2] ) ]
W.3 rms acoust[mW]
=aAvg(freq)[1] * W0 rms elec[mW] * 10-0.003 * freq[MHz] * z[mm]W0 rms elec[mW]
=10-3 * g(freq)[μモー] * Asq scan sum[1] * V total rms2 [V2]表面におけるA ispta (z=0):
A ispta(z=0)=(Lx(z=0)[mm] + Res-box-size(z=0)[mm]) * Ly(z=0)[mm]
(注:z>0とは異なる表面におけるIsptaを計算するためにレリーズ2.5を導入
した。z=0の場合、先の1.61公式に戻した。)
【0084】
【0085】
z>0の場合の A ispta の導出
Ispta=c * (PII.3 0 * FSC * BOF/ FT)
W.3=sum(PII.3) * Lx * Ly/FT
sum(PII)=PII 0 * FSC * num lines
Ispta=c * (PII.3 0 * FSC * num lines * Lx * Ly / FT * A ispta)
=c * (PII.3 0 * FSC * BOF/ FT)
ここに、
sum(PII)=走査の全ての縁の合計
FT =フレーム時
PII.3 0 =中心線(または、Isptaが採られる線)のPII.3
FSC =フローサンプルカウント(リバーブ、再発射、“動作促進処置”を含
む各走査線の発射の合計数)
BOF =ビーム重畳ファクタ
≧BOF =num-lines * Lx * Ly/A ispta
A ispta =num-lines * Lx * Ly/BOF
【0086】
幾つかの走査モードの組合わせ:
Ispta tは、多分分離した限度を用いて各モード毎に別々に管理される。しかし、合計Isptaはモード全体の和である。
【0087】
# 以下に、今では時代遅れのモデルを説明する:
【0088】
5)TI[1]
熱的インデックス
大開口:A aper>1cm2
小開口:A aper≦1cm2
【0089】
組合わせ TI :
N =フレーム内のN番目の送信成分
SUM =Nの合計またはNの適切な部分集合
MAX =2つの項の最大
X =S,C,またはBの1つ
とする。
一般的には:
TIX=MAX[ SUM( TIX N surface ),SUM( TIX N depth ) ]
詳述すれば:
TIS=MAX [SUM(TIS N scanned)+SUM(TIS N unscanned small aper),
SUM(TIS N unscanned large aper) ]
TISF≒SUM(TIS N unscanned large aper)(定義:特別な追加)
TIB=MAX [SUM(TIB N scanned), SUM(TIB N unscanned) ]
TIC=SUM(TIC N)
または、換言すれば:
TIS=MAX [SUM(TI A)+SUM(TI C), SUM(TI B)]
TIB=MAX [SUM(TI A), SUM(TI D)]
TIC=SUM(TI E)
TISF≒SUM(TI B)
【0090】
公式(表面):
式A:TI A=W01[mW] * freq[MHz]/210 * mW * MHz
(TISまたはTIB;走査された;表面)
式C:TI C=W0[mW] * freq[MHz]/210 * mW * MHz
(TIS;走査されていない;表面;小開口)
式E:TI E=W0[mW] * (10mm)/40 * Deq[mm] * mW
注:10mm(オリジナル 1cm)
(TIC;走査された、またはされていない;表面)
【0091】
公式(深さ):
式B:TI B=(freq[MHz]/210 * mW * MHz) *
max z>z bp[min(W.3(z)[mW], Ispta t(z)[mW/cm2] * 1cm2)]
(TIS;走査されていない;深さ;大開口)
【0092】
最大−最小決定:
最小(W, I)のz>z bpの最大に関しては3つのケースだけを検討すればよい。
ケース1:W.3(z bp)≦Ita.3(z bp) 従って最大−最小=W.3(z bp)
ケース2:W.3(z bp)>Ita.3(z bp)及びz>z bpに関して少なくとも1つのクロ
スオーバーが存在する。
z xoをそれらの最小であるものとすれば、
最大−最小=max[ W(z xo), Ita.3(z bp<z<z xo) ]
ケース3:W.3(z bp)>Ita.3(z bp)及びz>z bpに関してクロスオーバーは存在
しない。従って、
最大−最小=Ita.3(z)のz>z bp の最大
【0093】
要約:
式D:TI D=次の内の最小値
50mW-1 * [ W.3(zB.3)[mW] * Ispta.3(zB.3) [mW/cm2 ] * cm2,
4.4mW-1 * W.3(zB.3)[mW] ]1/2
(TIB;走査されていない;深さ)
ここに、
W0[mW]=W0 rms acoust[mW]=a(freq)[1] * W rms elec[mW]
(上記参照):平均音響パワー
W01[mW]=W01 rms acoust[mW]=a(freq)[1] * W01 rms elec[mW]
活動開口の中心の1cmを通る平均音響パワー
W01 rms elec[mW]=1000-1 * g(freq)[μモー] * A01[1] * V total rms2[V2]
W.3 [mW]=W.3 rms acoust [mW]=W0[mW] * 10-0.003 * freq[MHz] * z[mm]
Ispta t[mW/cm2] (上記参照)
A01[1]:1cm長開口(中心 1cm)について平方された1cmの走査され・合計さ
れたアポジゼーション値;成分のフレームのN番目の発射の場合中心1c
mで打ち切られたアポジゼーションプロファイルを平方し、合計する;
次いでこの成分のフレームを構成している全ての縁について平均する;
ACQによって計算
Deq[mm]=[ (4 / PI) * A aper[mm] ]1/2
z bp[mm]=1.5 * Deq[mm]
A aper[mm2]:活動開口面積
A aper[mm2]=a0y[mm] * a0x scanned 12db[mm]
a0x scanned 12db[mm]:走査された活動開口のサイズ、
即ち、V rms≧−12dB max(V rms)を有する全ての
要素からなる開口
ACQによって計算
もしA aper>1cm2=100mm2ならば、大開口
もしA aper≦1cm2=100mm2ならば、小開口
zB.3深さ:W.3(z) * Ispta t(z)を最大にするか、または同様に、0.6でデレートされたIsptaを最大にする
【0094】
成分の合計及びV依存性:
表面及びある深さにおけるTIを互いに加え合わせ、2つの合計の最大を求める。Wから計算された全てのTI項についてはV2依存性を使用し、Isptaまたは平方根(Ispta) から計算された全てのTI項についてはIspta Cファクタまたはその平方根を使用する。
N=N番目の成分とし、
Xバー=S, B, またはCとすれば、
各X=S, B, C毎に、
TIX=max (TIX surface, TIX depth )
注:電圧でスケーリング:全てのWは V2に比例し、一方IsptaはIsppa-C-factor 挙動を使用する。
− − − − − − − − − − − − − − − − − −
【0095】
6)Δ T : xdcr 表面における温度上昇
dT skin/air[C]
=1000-1k skin/air[C/W] * AperFct(SFSA)[1] * W0 rms elec[mW]
k skin/air:ヘッド容量[C/W];熱測定から;
gel/air及びgel/skinの場合はISから
AperFct:開口関数
AperFct=min( AperFct max, AperFct slope * SFSA+AperFct intercept )
ISからのAperFct max, AperFct slope, AperFct intercept:熱測定から
SFSA:静止断片表面積
活動走査された開口のサイズの測度
SFSA=成分に関する全ての発射のスパンされた開口/max-xdcr-aperture
以下のように近似する。
SFSA=a0x scanned 12db[mm] / max a0x aperture[mm]
「スパンされた」:スパンされた開口のための成分カウントについてあ
るフレーム中の何等かの発射の部分である全ての要素のことをいう。
W0 rms elec[mW]
=1000-1 * g(freq)[μモー] * Asq scan sum[1] * <V2> rms[V2]
注:A128に関して現れる走査ファクタ(SF)は、ACQによるW0 rms計算に
おいて斟酌される。
− − − − − − − − − − − − − − − − − −
【0096】
IEC1157 のためのパラメータ:
最大アンデレーテッドレアフラクショナル圧力: Pr.0 :
Mpaで表された最大アンデレーテッドレアフラクショナル圧力。
これは、次の2つの相違点を除いて、本質的にMIと同一である。
−周波数による除算は行わない。
−組織のデレーションは行わない。
【0097】
Ispta.0 :
アンデレーテッド空間ピーク時間平均強度
組織のデレーションを行わないIsptaと同一。モード及び順次焦点について加え合わせる。
W0 :
合計(アンデレーテッド)平均音響パワー
【0098】
Iob :
Iob=W0/max az aper * max el aper * max(sfsa)
W0:アンデレーテッド平均音響パワー(モード及び順次焦点について合計)
max az aper:最大方位開口
max el aper:最大迎角開口
sfsa: 静止断片的表面積(モード及び順次焦点について合計)
max(sfsa): 全モード及び順次焦点の最大
− − − − − − − − − − − − − − − − − −
【0099】
ヘッドルーム:
xdcr変化性を斟酌するために、音響及び熱的パラメータ(MI、Isppa、Ispta、TI、ΔT)をあるdBヘッドルームだけ上昇させることができる。このヘッドルームは周波数依存であることができ、イメージングスペックから得られる。(サンプルグリッド)周波数当たり2つのヘッドルーム、即ち音響パラメータのためのヘッドルーム、及び熱的パラメータのためのヘッドルームが存在する。
オンラインコードは、最後の値を繰り返し(一定補外)、次いで得られたdB値を使用して音響変換及びヘッド容量をスケールすることによって、周波数内を(dBで)補間するか、または補外する。
acoustDB:現出力搬送波周波数に補間された音響ヘッドルーム(dB)
thermDB:現出力搬送波周波数に補間された熱的ヘッドルーム(dB)
そこで
aAvg HR =aAvg * 10acoustDB/10
aMax HR =aMax * 10acoustDB/10
k skin HR=k skin * 10thermDB/10
k air HR =k air * 10thermDB/10
これらの変換及び特定熱ファクタは上式内に使用されている。
− − − − − − − − − − − − − − − − − −
【0100】
【0101】
【0102】
【0103】
Vpp[V] (Vmax=Vpp / 2) 正規化されたアポジゼーションのためのパルス
p−p
Vpp=Vpp ref * absolute scale
ここに、
Vpp refはCanのモデルによって計算される。
(Vmax またはVppは:
<V2> min(pulse shape(正規化されたアポジゼーションのための送
信電圧の自乗の最小)に置換)
従って:
V total rms2=V pulse rms2 * τ * FSC * num lines/frame time
ここに、
V pulse rmsは上記参照
τはCanのモデルにより立上がり時間として計算されたパルス持続時間FSC
* num linesは、フレーム内のモードのための発射の合計数を意味する。
【0104】
アポジゼーション合計の概要:
Aij=j番目の走査線についてのi番目の要素のアポジゼーション
Xij=j番目の走査線についてのi番目の要素の位置
Ni=要素の数
Nj=走査線の数
と定義すれば、
a0x=│X i1 j0− X i2 j0│のi1、i2の最大
j0=走査面積の中心線のインデックス
a0x scanned 12db=│X i1 j0− X i2 j0│のi1、i2の最大
i1及びi2は、V rms≧10-(12/20) * max overiVrms(iを有する要素の何
れかの対を表す
Asq unsc sum=ΣNi i=0│A i j0│2
j0=走査面積の中心線のインデックス
Asq scan sum=ΣNj j=0ΣNi i=0│A i j│2
A01=1/Nj ΣNj j=0ΣNi i=0 trunc 1cm(│A i j│2)
ここに、trunc 1cm(Aij2)は、
もしi番目が開口に沿う1cmの外側にあれば =0、
もしi番目が開口に沿う1cmの内側にあれば =Aij2
注:合計ではなく、線の数の平均を選択した。これは、Vrmsは何れかの線が発射
される時間にわたることを意味する。
Vrms total=Vrms pulse * τ-ppulse-length * FSC * num lines
/frame-time
但し、FSC * num lines は、(このモードの場合)フレーム当たりの合計発
射数を実際に意味する。
たとえA01の場合に合計をNj(線の合計数)によって除すとしても、それ
は合計に寄与するものではない。
【0105】
電力パラメータの概要:
W0 min el[W]
=10-6 * g(freq)[μモー] * Asq unsc sum[1] * <V2> min [V2]
W0 pa el[W]
=10-6 * g(freq)[μモー] * Asq unsc sum[1] * <V2> PA [V2]
W0 rms el[mW]
=10-3 * g(freq)[μモー] * Asq scan sum[1] * <V2> rms [V2]
W01 rms el[mW]
=10-3 * g(freq)[μモー] * A01 [1] * <V2> rms [V2]
− − − − − − − − − − − − − − − − − −
【0106】
実際の管理:
1)短時間定数(MI、 Isppa t )
これらは、モード特定である。2つの相で管理され、先ず、Vの“人工”値(これは、他のVの値をどのように“スケール”したかを告げるパラメータを含む)のための幾つかの値を推定する。
Isppa t(V)=Isppa t(V0) * V2 /V02 * [1−P+−* (1− V2 /V02)]
従って、Isppa t(V)は5つのパラメータによって決定される。
Isppa t(V)=fct(Isppa t(V0),V0,P+,P−,V)
MI(V) =fct(MI V0,V0,P+,P−,V)
従って、限界におけるIsppa tを与えるVを計算する。
2)長時間定数(TI、DT、 Ispta t )
a) Ispta :
3つのステップで管理される。
1)各モードは、ISからのモード最大値に制限される。
2)合計Isptaは、各モードをステップ1の後のその値に比例するように縮小させることによってISからの合計限度(FDA限度)に制限される。
3)各モードは、現在ユーザが選択したIspta限度に制限される。
b)TI及びDT
全てのモードが全フレームにわたって計算され、次いで加算公式に従って合計される。
c)管理
先ず、各モード毎に最大短時間電圧(音響的に、または熱的にの何れかに制限されている)を計算する。
この電圧に基づいて、各モードはIsptaを計算し、Ispta限度に合致させるためにもし必要ならば、電圧を縮小する。
TI及びDTが計算され、Vを各モード毎に比例的に縮小させてDTを限度以下にする。
電気的限度を計算し、もし必要ならば、Vを再び比例的に縮小させる。
最後に、もしユーザが、選択したMIまたはIsptaの下限を有しているか、または選択したdB縮小として有していれば、更に縮小させる。
【0107】
質問/ノート:
−注:ピークアポジゼーションは不要
−注:MIは、z PII.3 maxではなく、z MI maxにおいて評価される。
−注:ρ* c:密度/速度は水中値である。
−注:正規化されたアポジゼーションのためのV2
−注:Asq scan summedは、フレーム内の縁にわたって平均すべきである。
−注:res-box-sizeは、セクタ/ベクトル/湾曲した線形のためのアークである
−注:Ispta tのためのc sppa
−注:Ispta tは、ビーム面積及びc sppaのために中心線値を使用している。
−注:現在ではIspta tを追加しており、zについて最適化は行わない。
規則を確認するために、または規則は何を言っているかについて、オプションを使用することができる。
−注:たとえV I及びApod Iの両方がフレーム内のi番目の発射と共に変化しても、
W rmsの空間及び時間へのファクタリングは合理的である。
−注:A01は線発射にわたって平均され、線発射当たりの各要素にわたって合計される
−注:TIのV依存性:V2、cファクタ、または(cファクタ)1/2の何れか
−注:CWはフォールアウトすべきである。
−注:たとえブロード線がpan/res boxの外側にあっても、ブロードサイド線上の単一の線パラメータを根拠とする。
【0108】
?AuxCW?
単一線発射パラメータ(送信成分当たり)
Float outputCarrierFreqMHz; //Canのモデル
Float usLineAngleDeg;
Float xmtFocusMm;
Id baseApodTypeId;
Float xmtVppV;
Float xmtApertureMm; //湾曲した線形のアーク長
Float xdcrCurvatureRadiusMm;//線形スタックの場合無限大
または
Float xdcrAperChordMm; //線形スタックの場合=xmtApertureMm
//湾曲した線形スタックの場合=chord
//または代わりに曲率半径を使用
Float elevApertureMm;
Float refPulseVppV; //Canのモデル
Float refPulseRmsV; //Canのモデル
Float refPulseTauUSec; //Canのモデル
Float electConductanceUMho;
Float sumApodSquares; //単一線の場合:要素の合計
【0109】
走査パラメータ(送信成分当たり)
Float frameTimeUSec; //1つだけ(成分当たりではなく)
Float numFirings;
Float linearResBoxSizeMm;
Float angleResBoxSizeDeg; //セクタ、ベクトル、湾曲線形の場合
Float resRadiusOffsetMm; //セクタ、ベクトル、湾曲線形の場合
Float sumSumApodSquares;
//要素について合計された発射を合計したApod2
Float trunc1cmSumSumApodSq;
//要素について合計された発射を合計したApod2
//最大で送信開口の1cmで打ち切り
Float aperScanned12dbMm;
//最大−12db以内の走査されたapert.w.Vrms
【0110】
近似モード:
1)負の電圧振幅は、包絡線p−p電圧の1/2によって近似される。
2)aperScanned12dbは、SFSA * max-apertureによって近似される。
【0111】
CWモード:
CWは自然にモデルに適合し、それは別々に測定されない。
しかしながら、幾つかの入力パラメータは、幾つかの関係が保存される限り任意である。
1)RefPulseRmsV=1/4 * √2 * Vpp
2)τ* Num-pulses-per-frame=frame-time * mode-frame-time-fraction
但し、τ、Num..及びframe-timeは任意である。
【図面の簡単な説明】
【図1】 現在では好ましい実施の形態の医療診断用超音波イメージングシステムのブロック図である。
【図2】 現在では好ましい実施の形態の超音波イメージ内のユーザによって選択された領域における被送信超音波ビームの音響出力パラメータを決定する方法のフローチャートである。
【図3】 図2の好適な方法を例示する超音波イメージの例示であり、ユーザが選択した領域は単一の点である。
【図4】 図2の好適な方法を例示する超音波イメージの例示であり、ユーザが選択した領域は複数の点である。
【図5】 図2の好適な方法を例示する超音波イメージの例示であり、等圧線表現を使用して、ユーザが選択した領域の決定された音響出力パラメータを示す。
Claims (8)
- 超音波イメージを生成するように動作し、且つ超音波ビームを送信するように動作するトランスデューサ、および、プロセッサを備えている医療診断用超音波イメージングシステムと共に使用され、前記超音波イメージ内のユーザが選択した領域における被送信超音波ビームの音響出力パラメータを決定するためのプロセッサの作動方法において、
前記プロセッサが下記ステップを実行する、
(a)医療診断用超音波イメージングシステムのトランスデューサから超音波ビームを送信するステップと、
(b)前記医療診断用超音波イメージングシステムを用いて超音波イメージを生成するステップと、
(c)前記超音波イメージ内の領域の選択をユーザから受けるステップと、
(d)前記ユーザが選択した領域における被送信超音波ビームの音響出力パラメータを決定するステップと、を含み、
前記ユーザが選択した領域は複数の点からなり、前記ステップ(d)は、前記複数の各点毎にそれぞれの音響出力パラメータを決定することからなり、
更に、(e)前記複数の点のための前記決定された音響出力パラメータの分布を決定するステップ、
を含むことを特徴とする方法。 - (f)前記複数の点のための前記決定された音響出力パラメータの前記決定された分布の指示を供給するステップ、
を更に含むことを特徴とする請求項1に記載の方法。 - 超音波イメージを生成するように動作し、且つ超音波ビームを送信するように動作するトランスデューサ、および、プロセッサを備えている医療診断用超音波イメージングシステムと共に使用され、前記超音波イメージ内の選択された領域において被送信超音波ビームの指定された音響出力パラメータを達成するためのプロセッサの作動方法において、
前記プロセッサが下記ステップを実行する、
(a)医療診断用超音波イメージングシステムのトランスデューサから超音波ビームを送信するステップと、
(b)前記医療診断用超音波イメージングシステムを用いて超音波イメージを生成するステップと、
(c)前記超音波イメージ内にある領域を選択するステップと、
(d)前記選択された領域における被送信超音波ビームの音響出力パラメータを決定するステップと、を含み、
前記選択された領域は、複数の点からなり、前記ステップ(d)は、前記複数の各点毎にそれぞれの音響出力パラメータを決定することからなり、
更に、(e)前記医療診断用超音波イメージングシステムの動作パラメータを自動的に調整することによって、前記選択された領域において被送信超音波ビームの指定された音響出力パラメータを達成するステップと、
を含むことを特徴とする方法。 - 超音波イメージを生成するように動作し、且つ超音波ビームを送信するように動作するトランスデューサ、および、プロセッサを備えている医療診断用超音波イメージングシステムと共に使用され、前記超音波イメージ内の選択された領域における被送信超音波ビームの音響出力パラメータの指示を供給するためのプロセッサの作動方法において、
前記プロセッサが下記ステップを実行する、
(a)医療診断用超音波イメージングシステムのトランスデューサから超音波ビームを送信するステップと、
(b)前記医療診断用超音波イメージングシステムを用いて超音波イメージを生成するステップと、
(c)前記被送信超音波ビームのピーク音響出力パラメータを含む領域とは異なる領域を前記超音波イメージ内に選択するステップと、
(d)前記選択された領域における被送信超音波ビームの音響出力パラメータを決定するステップと、を含み、
前記選択した領域は複数の点からなり、前記ステップ(d)は、前記複数の各点毎にそれぞれの音響出力パラメータを決定することからなり、
更に、(e)前記決定された音響出力パラメータの指示を供給するステップと、
(f)前記複数の点のための前記決定された音響出力パラメータの分布を決定するステップと、
(g)前記複数の点のための前記決定された音響出力パラメータの前記決定された分布の指示を供給するステップと、
を含むことを特徴とする方法。 - 超音波イメージを生成するように動作し、且つ超音波ビームを送信するように動作するトランスデューサ、および、プロセッサを備えている医療診断用超音波イメージングシステムと共に使用され、前記超音波イメージ内の選択された領域における被送信超音波ビームの音響出力パラメータの指示を供給するためのプロセッサの作動方法において、
前記プロセッサが下記ステップを実行する、
(a)医療診断用超音波イメージングシステムのトランスデューサから超音波ビームを送信するステップと、
(b)前記医療診断用超音波イメージングシステムを用いて超音波イメージを生成するステップと、
(c)前記被送信超音波ビームのピーク音響出力パラメータを含む領域とは異なる領域を前記超音波イメージ内に選択するステップと、
(d)前記選択された領域における被送信超音波ビームの音響出力パラメータを決定するステップと、
前記選択された領域は複数の点からなり、前記ステップ(d)は、前記複数の各点毎にそれぞれの音響出力パラメータを決定することからなり、
(e)前記決定された音響出力パラメータの指示を供給するステップと、
(f)前記被送信超音波ビームの音響減衰を決定するステップと、
(g)前記決定された音響出力パラメータを、前記決定された音響減衰に関して補正するステップと、
を含むことを特徴とする方法。 - 超音波イメージを生成するように動作し、且つ超音波ビームを送信するように動作するトランスデューサ、および、プロセッサを備えている医療診断用超音波イメージングシステムと共に使用され、前記超音波イメージ内の選択された領域における被送信超音波ビームの音響出力パラメータの指示を供給するためのプロセッサの作動方法において、
前記プロセッサが下記ステップを実行する、
(a)医療診断用超音波イメージングシステムのトランスデューサから超音波ビームを送信するステップと、
(b)前記医療診断用超音波イメージングシステムを用いて超音波イメージを生成するステップと、
(c)前記被送信超音波ビームのピーク音響出力パラメータを含む領域とは異なる領域を前記超音波イメージ内に選択するステップと、
(d)前記選択された領域における被送信超音波ビームの音響出力パラメータを決定するステップと、
前記選択された領域は複数の点からなり、前記ステップ(d)は、前記複数の各点毎にそれぞれの音響出力パラメータを決定することからなり、
(e)前記決定された音響出力パラメータの指示を供給するステップと、を含み、
前記音響出力パラメータは音響モデルを使用して決定され、前記音響モデルを、異なる非線形レスポンスレベルを有するコントラスト剤の生体内測定によって較正するステップを更に含むことを特徴とする方法。 - 超音波イメージを生成するように動作し、且つ超音波ビームを送信するように動作するトランスデューサ、および、プロセッサを備えている医療診断用超音波イメージングシステムと共に使用され、前記超音波イメージ内の被送信超音波ビームの音響出力パラメータが決定される領域の位置の指示を供給するためのプロセッサの作動方法において、
前記プロセッサが下記ステップを実行する、
(a)医療診断用超音波イメージングシステムのトランスデューサから超音波ビームを送信するステップと、
(b)前記医療診断用超音波イメージングシステムを用いて超音波イメージを生成するステップと、
(c)前記超音波イメージ全体より小さい前記超音波イメージ内の領域における前記被送信超音波ビームの音響出力パラメータを決定するステップと、
前記領域は複数の点からなり、前記ステップ(c)は、前記複数の各点毎にそれぞれの音響出力パラメータを決定することからなり、
(d)前記決定された音響出力パラメータの指示を供給するステップと、
(e)前記超音波イメージ内の前記領域の位置の指示を供給するステップと、を含み、
前記ステップ(e)は、前記領域の範囲の指示を供給することからなる、
ことを特徴とする方法。 - 前記音響出力パラメータは、薬品運搬容器に影響を及ぼすように作動することを特徴とする請求項1乃至7の何れか1項に記載の方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/501,443 US6413218B1 (en) | 2000-02-10 | 2000-02-10 | Medical diagnostic ultrasound imaging system and method for determining an acoustic output parameter of a transmitted ultrasonic beam |
US09/501,443 | 2000-02-10 | ||
PCT/US2001/003534 WO2001058357A1 (en) | 2000-02-10 | 2001-02-02 | Medical diagnostic ultrasound imaging system and method for determining an acoustic output parameter of a transmitted ultrasonic beam |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2004514461A JP2004514461A (ja) | 2004-05-20 |
JP2004514461A5 JP2004514461A5 (ja) | 2005-01-20 |
JP4847665B2 true JP4847665B2 (ja) | 2011-12-28 |
Family
ID=23993595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001557472A Expired - Lifetime JP4847665B2 (ja) | 2000-02-10 | 2001-02-02 | 医療診断用超音波イメージングシステム及び被送信超音波ビームの音響出力パラメータを決定する方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US6413218B1 (ja) |
JP (1) | JP4847665B2 (ja) |
AU (1) | AU2001233275A1 (ja) |
DE (1) | DE10195485B3 (ja) |
WO (1) | WO2001058357A1 (ja) |
Families Citing this family (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6714667B1 (en) * | 2000-07-14 | 2004-03-30 | Koninklijke Philips Electronics N.V. | User interface for imaging system |
JP2002113006A (ja) * | 2000-10-10 | 2002-04-16 | Toshiba Medical System Co Ltd | 診断機器設定方法および診断装置 |
US7127401B2 (en) * | 2001-03-12 | 2006-10-24 | Ge Medical Systems Global Technology Company, Llc | Remote control of a medical device using speech recognition and foot controls |
JP4373699B2 (ja) * | 2003-04-28 | 2009-11-25 | 株式会社東芝 | 超音波診断装置 |
US20060030780A1 (en) * | 2004-08-03 | 2006-02-09 | Jean-Francois Gelly | System and method providing controllable attenuation of an ultrasound probe |
US20060030779A1 (en) * | 2004-08-09 | 2006-02-09 | Siemens Medical Solutions Usa, Inc. | Intensity level control for contrast agent imaging |
US8043219B2 (en) * | 2004-09-17 | 2011-10-25 | Siemens Medical Solutions Usa, Inc. | Automated power level for contrast agent imaging |
US20060074315A1 (en) * | 2004-10-04 | 2006-04-06 | Jianming Liang | Medical diagnostic ultrasound characterization of cardiac motion |
US20060074312A1 (en) * | 2004-10-06 | 2006-04-06 | Bogdan Georgescu | Medical diagnostic ultrasound signal extraction |
US7780602B2 (en) * | 2004-12-27 | 2010-08-24 | General Electric Company | Method and system for controlling an ultrasound system |
US7921717B2 (en) * | 2005-01-03 | 2011-04-12 | Siemens Medical Solutions Usa, Inc. | Ultrasonic imaging system |
US20060247544A1 (en) * | 2005-02-03 | 2006-11-02 | Maleeha Qazi | Characterization of cardiac motion with spatial relationship |
JP2007301181A (ja) * | 2006-05-11 | 2007-11-22 | Ge Medical Systems Global Technology Co Llc | 超音波診断装置および画像表示方法 |
JP4465018B2 (ja) * | 2008-09-09 | 2010-05-19 | オリンパスメディカルシステムズ株式会社 | 超音波診断装置 |
US20120165665A1 (en) * | 2010-12-22 | 2012-06-28 | Samsung Medison Co., Ltd. | Method for providing mechanical index map and/or pressure map based on depth value and diagnostic ultrasound system using the method |
US9629610B2 (en) * | 2010-12-22 | 2017-04-25 | Samsung Medison Co., Ltd. | Method of operating ultrasound diagnosis apparatus for providing map of interest index and ultrasound diagnosis apparatus using the method |
EP2468191B1 (en) * | 2010-12-22 | 2019-01-30 | Samsung Medison Co., Ltd. | Ultrasound diagnosis apparatus for providing map of interest index |
KR101219465B1 (ko) * | 2011-11-30 | 2013-01-09 | 삼성메디슨 주식회사 | 관심 지표 맵을 이용한 초음파 영상 보정 방법 및 초음파 영상 보정 장치 |
US9168027B2 (en) | 2013-02-22 | 2015-10-27 | Siemens Medical Solutions Usa, Inc. | Adaptive acoustic pressure estimation in medical ultrasound |
KR101501479B1 (ko) * | 2013-05-09 | 2015-03-11 | 알피니언메디칼시스템 주식회사 | 초음파 최적화 방법과 그를 위한 초음파 의료 장치 |
TWI485420B (zh) * | 2013-09-27 | 2015-05-21 | Univ Nat Taiwan | 超音波影像補償方法 |
EP3220829B1 (en) | 2014-11-18 | 2022-03-09 | C. R. Bard, Inc. | Ultrasound imaging system having automatic image presentation |
EP3220828B1 (en) | 2014-11-18 | 2021-12-22 | C.R. Bard, Inc. | Ultrasound imaging system having automatic image presentation |
KR20180085247A (ko) | 2017-01-18 | 2018-07-26 | 삼성메디슨 주식회사 | 초음파 영상 장치 및 초음파 영상 표시 방법 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08196537A (ja) * | 1995-01-31 | 1996-08-06 | Toshiba Corp | 超音波診断装置 |
JPH08280674A (ja) * | 1995-04-14 | 1996-10-29 | Toshiba Corp | 超音波診断装置 |
JPH1199152A (ja) * | 1997-09-26 | 1999-04-13 | Ge Yokogawa Medical Systems Ltd | 超音波診断装置 |
JPH11155858A (ja) * | 1997-11-26 | 1999-06-15 | Toshiba Corp | 超音波診断装置および超音波診断方法 |
JPH11318901A (ja) * | 1998-05-11 | 1999-11-24 | Toshiba Corp | 超音波診断装置 |
JP2000198521A (ja) * | 1998-12-28 | 2000-07-18 | Tiger Kawashima Co Ltd | 揚穀計量機 |
Family Cites Families (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5313947A (en) * | 1985-02-08 | 1994-05-24 | University Patents, Inc. | CW and pulsed doppler diagnostic system |
US5113706A (en) * | 1990-07-03 | 1992-05-19 | Hewlett-Packard Company | Ultrasound system with dynamic transmit focus |
US5357962A (en) * | 1992-01-27 | 1994-10-25 | Sri International | Ultrasonic imaging system and method wtih focusing correction |
US5331964A (en) * | 1993-05-14 | 1994-07-26 | Duke University | Ultrasonic phased array imaging system with high speed adaptive processing using selected elements |
US5509413A (en) * | 1993-08-11 | 1996-04-23 | Kabushiki Kaisha Toshiba | Ultrasonic diagnostic apparatus |
JP3059042B2 (ja) * | 1994-02-22 | 2000-07-04 | フクダ電子株式会社 | 超音波診断装置 |
US5551434A (en) * | 1994-06-22 | 1996-09-03 | Kabushiki Kaisha Toshiba | Ultrasonic imaging diagnosis apparatus |
US5482046A (en) * | 1994-11-23 | 1996-01-09 | General Electric Company | Acoustic power control technique |
US5634465A (en) | 1995-06-09 | 1997-06-03 | Advanced Technology Laboratories, Inc. | Continuous display of cardiac blood flow information |
US5709209A (en) * | 1996-03-29 | 1998-01-20 | Siemens Medical Systems, Inc. | Ultrasound signal processing system |
US5882306A (en) * | 1997-04-11 | 1999-03-16 | Acuson Corporation | Ultrasound imaging methods and systems |
JP3723663B2 (ja) * | 1997-07-15 | 2005-12-07 | フクダ電子株式会社 | 超音波診断装置 |
US6077225A (en) * | 1998-01-23 | 2000-06-20 | Hewlett-Packard Company | Ultrasound method for enhancing image presentation when contrast agents are used |
US5891040A (en) * | 1998-02-18 | 1999-04-06 | Hewlett-Packard Company | Method for maintaining a constant velocity to color map in an ultrasound flow imaging system |
JP4116143B2 (ja) * | 1998-04-10 | 2008-07-09 | 株式会社東芝 | 超音波診断装置 |
US6200267B1 (en) * | 1998-05-13 | 2001-03-13 | Thomas Burke | High-speed ultrasound image improvement using an optical correlator |
US6036643A (en) | 1998-05-14 | 2000-03-14 | Advanced Technology Laboratories, Inc. | Ultrasonic harmonic doppler imaging |
US6123670A (en) * | 1998-12-15 | 2000-09-26 | General Electric Company | Ultrasound imaging with optimal image quality in region of interest |
US6080107A (en) * | 1999-01-26 | 2000-06-27 | Hewlett-Packard Company | Methods for the use of contrast agents in ultrasonic imaging |
US6217516B1 (en) * | 1999-11-09 | 2001-04-17 | Agilent Technologies, Inc. | System and method for configuring the locus of focal points of ultrasound beams |
US6258033B1 (en) * | 1999-11-30 | 2001-07-10 | Agilent Technologies, Inc. | Ultrasound method employing echoes from a region of interest to enable quantization of backscatter signals |
US6210335B1 (en) * | 1999-12-08 | 2001-04-03 | General Electric Company | Acoustic flash to increase penetration |
-
2000
- 2000-02-10 US US09/501,443 patent/US6413218B1/en not_active Expired - Lifetime
-
2001
- 2001-02-02 AU AU2001233275A patent/AU2001233275A1/en not_active Abandoned
- 2001-02-02 DE DE10195485T patent/DE10195485B3/de not_active Expired - Lifetime
- 2001-02-02 JP JP2001557472A patent/JP4847665B2/ja not_active Expired - Lifetime
- 2001-02-02 WO PCT/US2001/003534 patent/WO2001058357A1/en active Application Filing
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08196537A (ja) * | 1995-01-31 | 1996-08-06 | Toshiba Corp | 超音波診断装置 |
JPH08280674A (ja) * | 1995-04-14 | 1996-10-29 | Toshiba Corp | 超音波診断装置 |
JPH1199152A (ja) * | 1997-09-26 | 1999-04-13 | Ge Yokogawa Medical Systems Ltd | 超音波診断装置 |
JPH11155858A (ja) * | 1997-11-26 | 1999-06-15 | Toshiba Corp | 超音波診断装置および超音波診断方法 |
JPH11318901A (ja) * | 1998-05-11 | 1999-11-24 | Toshiba Corp | 超音波診断装置 |
JP2000198521A (ja) * | 1998-12-28 | 2000-07-18 | Tiger Kawashima Co Ltd | 揚穀計量機 |
Also Published As
Publication number | Publication date |
---|---|
JP2004514461A (ja) | 2004-05-20 |
US6413218B1 (en) | 2002-07-02 |
AU2001233275A1 (en) | 2001-08-20 |
WO2001058357A1 (en) | 2001-08-16 |
DE10195485B3 (de) | 2013-06-27 |
DE10195485T1 (de) | 2003-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4847665B2 (ja) | 医療診断用超音波イメージングシステム及び被送信超音波ビームの音響出力パラメータを決定する方法 | |
US20230243966A1 (en) | Imaging methods and apparatuses for performing shear wave elastography imaging | |
EP1551305B1 (en) | Instantaneous ultrasonic measurement of bladder volume | |
US8002704B2 (en) | Method and system for determining contact along a surface of an ultrasound probe | |
Hinkelman et al. | Measurement and correction of ultrasonic pulse distortion produced by the human breast | |
US5235984A (en) | On-line acoustic densitometry tool for use with an ultrasonic imaging system | |
JP5185921B2 (ja) | 超音波画像診断のための方法および装置 | |
JP5760080B2 (ja) | せん断波を使用する撮像方法および装置 | |
US9168027B2 (en) | Adaptive acoustic pressure estimation in medical ultrasound | |
KR101501479B1 (ko) | 초음파 최적화 방법과 그를 위한 초음파 의료 장치 | |
US6858008B2 (en) | Automatic ultrasound transmit power setting method and system | |
US4511984A (en) | Ultrasound diagnostic apparatus | |
US20230404537A1 (en) | Ultrasound medical imaging with optimized speed of sound based on fat fraction | |
US11986356B2 (en) | Reduction of reverberation artifacts in ultrasound images and associated devices, systems, and methods | |
Bouakaz et al. | Harmonic ultrasonic field of medical phased arrays: Simulations and measurements | |
KR20170045985A (ko) | 초음파 영상장치 및 그 제어방법 | |
US7481769B2 (en) | Ultrasonic diagnosing apparatus | |
JP2022550020A (ja) | 心臓系をモデル化するための方法及びシステム | |
JP3474233B2 (ja) | 超音波診断装置 | |
US20220211342A1 (en) | Method Of Performing Automated Measurements Over Multiple Cardiac Cycles | |
Hemmsen | Image processing in medical ultrasound | |
Dolan | H EW LETT. PAC KAFI DJOU FI NAL | |
CN112087971A (zh) | 超声成像系统和方法 | |
Nakagawa et al. | Ultrasonic nonlinear parameter CT using nonlinear interaction of sound waves | |
JPH04183449A (ja) | 皮脂厚表示装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080129 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100715 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20101015 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20101026 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110117 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20110217 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20110324 |
|
RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20110401 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110509 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110615 |
|
A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20110628 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110817 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110823 |
|
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: 20110916 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20111014 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20141021 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4847665 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 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 |
|
EXPY | Cancellation because of completion of term |