JP2013529098A - せん断波を使用する撮像方法および装置 - Google Patents
せん断波を使用する撮像方法および装置 Download PDFInfo
- Publication number
- JP2013529098A JP2013529098A JP2013505551A JP2013505551A JP2013529098A JP 2013529098 A JP2013529098 A JP 2013529098A JP 2013505551 A JP2013505551 A JP 2013505551A JP 2013505551 A JP2013505551 A JP 2013505551A JP 2013529098 A JP2013529098 A JP 2013529098A
- Authority
- JP
- Japan
- Prior art keywords
- excitation
- shear
- viscoelastic medium
- original image
- shear wave
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 56
- 230000005284 excitation Effects 0.000 claims abstract description 99
- 238000012935 Averaging Methods 0.000 claims abstract description 11
- YTAHJIFKAKIKAV-XNMGPUDCSA-N [(1R)-3-morpholin-4-yl-1-phenylpropyl] N-[(3S)-2-oxo-5-phenyl-1,3-dihydro-1,4-benzodiazepin-3-yl]carbamate Chemical compound O=C1[C@H](N=C(C2=C(N1)C=CC=C2)C1=CC=CC=C1)NC(O[C@H](CCN1CCOCC1)C1=CC=CC=C1)=O YTAHJIFKAKIKAV-XNMGPUDCSA-N 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 16
- 238000012512 characterization method Methods 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims 1
- 230000006835 compression Effects 0.000 description 13
- 238000007906 compression Methods 0.000 description 13
- 238000002604 ultrasonography Methods 0.000 description 13
- 239000000523 sample Substances 0.000 description 11
- 239000002245 particle Substances 0.000 description 7
- 238000002592 echocardiography Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- 241001465754 Metazoa Species 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012285 ultrasound imaging Methods 0.000 description 2
- 102000008186 Collagen Human genes 0.000 description 1
- 108010035532 Collagen Proteins 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 229920001436 collagen Polymers 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 210000001087 myotubule Anatomy 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
-
- G—PHYSICS
- 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/52046—Techniques for image enhancement involving transmitter or receiver
-
- 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/5206—Two-dimensional coordinated display of distance and direction; B-scan display
- G01S7/52065—Compound scan display, e.g. panoramic imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52071—Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52077—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 with means for elimination of unwanted signals, e.g. noise or interference
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
粘弾性媒質(2)を観測するためにせん断波を使用する撮像方法であって、その間に弾性せん断波が、撮像デバイスによって粘弾性媒質(2)中の異なる励起場所(Lj)でそれぞれ生成されるいくつかの連続励起ステップjであって、異なる励起場所は、最大距離Dmだけ互いに分離される、いくつかの連続励起ステップjと、時間tkでの粘弾性媒質の連続原画像Imj(tk)の組jが、せん断波の伝搬の間に決定される、各励起ステップjに対応する撮像ステップであって、原画像は、最大距離Dmよりも大きい分解能Rを有する、撮像ステップと、同じ相対時間tkに対応する原画像Imj(tk)が、平均画像Im'(tk)を決定するために平均化される、平均化ステップとを含む、撮像方法。
Description
本発明は、せん断波を使用する撮像方法および装置に関する。
より詳しくは、本発明は、粘弾性媒質を観測するためにせん断波を使用する撮像方法であって、
- その間に弾性せん断波が、撮像デバイスによって粘弾性媒質中の励起場所で生成される励起ステップjと、
- 粘弾性媒質の一組の連続画像が、前記粘弾性媒質中のせん断波の伝搬の間に前記撮像デバイスによって決定される撮像ステップとを含む、撮像方法に関する。
- その間に弾性せん断波が、撮像デバイスによって粘弾性媒質中の励起場所で生成される励起ステップjと、
- 粘弾性媒質の一組の連続画像が、前記粘弾性媒質中のせん断波の伝搬の間に前記撮像デバイスによって決定される撮像ステップとを含む、撮像方法に関する。
文書米国特許第7252004(B2)号は、そのような方法の例を述べ、そこではせん断波の伝搬の動画を構成する連続画像が、粘弾性媒質中の1つまたはいくつかの点でせん断波の少なくとも1つの伝搬パラメータを計算することによって、粘弾性媒質をマッピングするために使用される。
この周知の方法は、例えば癌領域または同様のものを見つけるための特に良好な結果を与えるけれども、そのような方法は、ある状況では、例えば高い減衰および回折を生じる非常に複雑な構造を有する生体媒質(例えば、せん断波が筋線維に対して横断方向に伝播されるときの動物の筋肉)では比較的低い信号対雑音比を有する画像をもたらすこともある。
本発明の特定の目的は、この欠点を軽減することである。
この目的のために、本発明は、粘弾性媒質を観測するためにせん断波を使用する撮像方法であって、
- その間に弾性せん断波が、撮像デバイスによって粘弾性媒質中の異なる励起場所でそれぞれ生成されるいくつかの連続励起ステップjであって、前記異なる励起場所は、最大距離Dmだけ互いに分離される、いくつかの連続励起ステップjと、
- 時間tkでの粘弾性媒質の連続原画像Imj(tk)の組jが、前記粘弾性媒質中のせん断波の伝搬の間に前記撮像デバイスによって決定される、各励起ステップjに対応する撮像ステップであって、前記原画像は、前記最大距離Dmに少なくとも等しい分解能Rを有し、tkは、対応するせん断波の生成から数えられる相対時間である、撮像ステップと、
- 同じ相対時間tkに対応する原画像Imj(tk)が、前記相対時間tkに対応する平均画像Im'(tk)を決定するために平均化される平均化ステップとを含む、撮像方法を提案する。
- その間に弾性せん断波が、撮像デバイスによって粘弾性媒質中の異なる励起場所でそれぞれ生成されるいくつかの連続励起ステップjであって、前記異なる励起場所は、最大距離Dmだけ互いに分離される、いくつかの連続励起ステップjと、
- 時間tkでの粘弾性媒質の連続原画像Imj(tk)の組jが、前記粘弾性媒質中のせん断波の伝搬の間に前記撮像デバイスによって決定される、各励起ステップjに対応する撮像ステップであって、前記原画像は、前記最大距離Dmに少なくとも等しい分解能Rを有し、tkは、対応するせん断波の生成から数えられる相対時間である、撮像ステップと、
- 同じ相対時間tkに対応する原画像Imj(tk)が、前記相対時間tkに対応する平均画像Im'(tk)を決定するために平均化される平均化ステップとを含む、撮像方法を提案する。
それ故に、本発明は、原画像での雑音の大部分を除去することを可能にし、より高い信号対雑音比を有する平均化画像を得ることを可能にする。従って、平均化画像は、せん断波伝搬パラメータのより良好なマッピングを得るために、かつオプションとして粘弾性媒質のあるレオロジー値および特に粘性効果を表すレオロジー値を得ることを可能にする前記伝搬パラメータのスペクトル(すなわち、周波数の関数としての前記伝搬パラメータの値)を得るためにより効果的に処理することができる。
本発明の方法の様々な実装形態では、次の処理の1つまたは複数に頼ることもまたオプションとして可能であり、
- 本方法は、撮像ステップと平均化ステップとの間に、原画像Imj(tk)が、異なる励起場所が原画像Imj(tk)で相互に対応して位置決めされるように空間的に再位置決めされる再位置決めステップをさらに含み、
- 前記原画像Imj(tk)は、異なる励起場所が原画像Imj(tk)で同じ座標を有するように各原画像の座標系をオフセットすることによって空間的に再位置決めされ、
- 連続原画像の前記組jの数Jは、少なくとも5であり、
- 各励起場所は、せいぜいR/5であり、好ましくはせいぜいR/10であるピッチだけ隣接励起場所から分離され、
- 励起場所は、励起ステップごとに撮像デバイスによって粘弾性媒質中で移動され、
- 粘弾性媒質は、励起ステップごとに励起場所が撮像デバイスに対して固定されたままである間に移動され、その結果励起場所は、励起ステップごとに粘弾性媒質中で移動され、
- 前記せん断励起は、前記撮像デバイスに属する超音波振動子のアレイによって粘弾性媒質中に放出される少なくとも1つの超音波によって生成され、
- 前記せん断励起は、外部の機械力を粘弾性媒質に各励起場所で加えることによって生成され、
- 各原画像Imj(tk)は、超音波撮像によってか、または MRIによって得られ、
- 本方法は、その間に時間にわたる前記平均画像Im'(tk)に基づいて、少なくとも1つのせん断波伝搬パラメータが、粘弾性媒質の少なくとも1つの点で計算される特徴付けステップを含み、
- 特徴付けステップの間に計算されるせん断波伝搬パラメータは、せん断波速度、せん断係数、ヤング率、せん断波減衰率、せん断弾性、せん断粘度、機械的緩和時間および異方性から選択され、
- 前記励起は、少なくとも500Hzの帯域幅を有し、前記せん断波伝搬パラメータは、前記帯域幅内で複数の周波数で計算される。
- 本方法は、撮像ステップと平均化ステップとの間に、原画像Imj(tk)が、異なる励起場所が原画像Imj(tk)で相互に対応して位置決めされるように空間的に再位置決めされる再位置決めステップをさらに含み、
- 前記原画像Imj(tk)は、異なる励起場所が原画像Imj(tk)で同じ座標を有するように各原画像の座標系をオフセットすることによって空間的に再位置決めされ、
- 連続原画像の前記組jの数Jは、少なくとも5であり、
- 各励起場所は、せいぜいR/5であり、好ましくはせいぜいR/10であるピッチだけ隣接励起場所から分離され、
- 励起場所は、励起ステップごとに撮像デバイスによって粘弾性媒質中で移動され、
- 粘弾性媒質は、励起ステップごとに励起場所が撮像デバイスに対して固定されたままである間に移動され、その結果励起場所は、励起ステップごとに粘弾性媒質中で移動され、
- 前記せん断励起は、前記撮像デバイスに属する超音波振動子のアレイによって粘弾性媒質中に放出される少なくとも1つの超音波によって生成され、
- 前記せん断励起は、外部の機械力を粘弾性媒質に各励起場所で加えることによって生成され、
- 各原画像Imj(tk)は、超音波撮像によってか、または MRIによって得られ、
- 本方法は、その間に時間にわたる前記平均画像Im'(tk)に基づいて、少なくとも1つのせん断波伝搬パラメータが、粘弾性媒質の少なくとも1つの点で計算される特徴付けステップを含み、
- 特徴付けステップの間に計算されるせん断波伝搬パラメータは、せん断波速度、せん断係数、ヤング率、せん断波減衰率、せん断弾性、せん断粘度、機械的緩和時間および異方性から選択され、
- 前記励起は、少なくとも500Hzの帯域幅を有し、前記せん断波伝搬パラメータは、前記帯域幅内で複数の周波数で計算される。
さらにその上、本発明はまた、粘弾性媒質を観測するためにせん断波を使用する先行クレームのいずれかによる方法を実装するための撮像装置も提供し、前記装置は、
- いくつかの連続励起ステップjの間に粘弾性媒質中の異なる励起場所でそれぞれ弾性せん断波を生成するように適合され、前記異なる励起場所は、最大距離Dmだけ互いに分離され、
- 各励起ステップjについてそれぞれ前記粘弾性媒質中のせん断波の伝搬の間に時間tkでの粘弾性媒質の連続原画像Imj(tk)の組jを決定するように適合され、前記原画像は、前記最大距離Dmに少なくとも等しい分解能Rを有し、tkは、対応するせん断波の生成から数えられる相対時間であり、
- 前記相対時間tkについて平均画像Im'(tk)を決定するために同じ相対時間tkに対応する原画像Imj(tk)を平均化するように適合される、少なくとも1つの電子中央ユニットを含む。
- いくつかの連続励起ステップjの間に粘弾性媒質中の異なる励起場所でそれぞれ弾性せん断波を生成するように適合され、前記異なる励起場所は、最大距離Dmだけ互いに分離され、
- 各励起ステップjについてそれぞれ前記粘弾性媒質中のせん断波の伝搬の間に時間tkでの粘弾性媒質の連続原画像Imj(tk)の組jを決定するように適合され、前記原画像は、前記最大距離Dmに少なくとも等しい分解能Rを有し、tkは、対応するせん断波の生成から数えられる相対時間であり、
- 前記相対時間tkについて平均画像Im'(tk)を決定するために同じ相対時間tkに対応する原画像Imj(tk)を平均化するように適合される、少なくとも1つの電子中央ユニットを含む。
本発明の装置の様々な実装形態では、次の処理の1つまたは複数に頼ることもまたオプションとして可能であり、
- 前記電子中央ユニットはさらに、異なる励起場所が原画像Imj(tk)で相互に対応して位置決めされるように、前記画像を平均化する前に前記原画像Imj(tk)を空間的に再位置決めするように適合され、
- 前記電子中央ユニットは、異なる励起場所が原画像Imj(tk)で同じ座標を有するように各原画像の座標系をオフセットすることによって前記原画像Imj(tk)を空間的に再位置決めするように適合され、
- 前記電子中央ユニットは、励起ステップごとに粘弾性媒質中の励起場所を移動させるように適合され、
- 前記電子中央ユニットは、時間にわたる連続平均画像に基づいて粘弾性媒質の少なくとも1つの点で少なくとも1つのせん断波伝搬パラメータを計算するように適合され、前記せん断波伝搬パラメータは、せん断波速度、せん断係数、ヤング率、せん断波減衰率、せん断弾性、せん断粘度、機械的緩和時間および異方性から選択され、
- 前記装置は、少なくとも500Hzの帯域幅を持つ前記励起を加えるように適合され、前記電子中央ユニットは、前記帯域幅での複数の周波数で前記せん断波伝搬パラメータを計算するように適合される。
- 前記電子中央ユニットはさらに、異なる励起場所が原画像Imj(tk)で相互に対応して位置決めされるように、前記画像を平均化する前に前記原画像Imj(tk)を空間的に再位置決めするように適合され、
- 前記電子中央ユニットは、異なる励起場所が原画像Imj(tk)で同じ座標を有するように各原画像の座標系をオフセットすることによって前記原画像Imj(tk)を空間的に再位置決めするように適合され、
- 前記電子中央ユニットは、励起ステップごとに粘弾性媒質中の励起場所を移動させるように適合され、
- 前記電子中央ユニットは、時間にわたる連続平均画像に基づいて粘弾性媒質の少なくとも1つの点で少なくとも1つのせん断波伝搬パラメータを計算するように適合され、前記せん断波伝搬パラメータは、せん断波速度、せん断係数、ヤング率、せん断波減衰率、せん断弾性、せん断粘度、機械的緩和時間および異方性から選択され、
- 前記装置は、少なくとも500Hzの帯域幅を持つ前記励起を加えるように適合され、前記電子中央ユニットは、前記帯域幅での複数の周波数で前記せん断波伝搬パラメータを計算するように適合される。
本発明の他の特徴および利点は、限定されない例として与えられ、添付の図面を参照するそれらの実施形態の次の説明から明らかになる。
様々な図では、類似の参照番号は、同一のまたは同様の要素を指定する。
図1で示される撮像デバイス1は、超音波を圧縮状態で拡散する粘弾性媒質2中の弾性せん断波の伝搬を調べるためであり、その粘弾性媒質は、例えば、
- 特に工業的応用での品質管理のために、不活性体によって、または
- 医学的応用での生体、例えば患者の身体の一部によって構成されることもある。
- 特に工業的応用での品質管理のために、不活性体によって、または
- 医学的応用での生体、例えば患者の身体の一部によって構成されることもある。
例として、これらの動きは、超音波圧縮波を媒質2中にそれの外表面3から送る役目を果たす、マイクロコンピュータ4(少なくともキーボード、その他などの入力インターフェース4b、およびスクリーン、その他などの出力インターフェース4aを含む)または任意の他の電子中央ユニットを用いて追跡される。これらの圧縮波は、媒質2に含有される拡散粒子5と相互作用し、その粒子は、超音波圧縮波を反射する。粒子5は、媒質2中の任意の不均一性によって構成されることもあり、特に医学的応用では、それらは、人体組織に存在するコラーゲンの粒子によって構成されることもある(超音波検査画像では、そのような粒子は、「スペックル」として周知の点を形成する)。
せん断波の伝搬を観測するために、観測媒質2の外表面3に接触して配置される超音波プローブ6が、使用される。このプローブは、軸Xに沿って超音波圧縮波パルスを供給し、そのパルスは、例えば0.5MHzから100MHzの範囲、好ましくは0.5MHzから15MHzの範囲にあり、例えば約4MHzである周波数を有する、超音波検査で一般に使用される種類である。
超音波プローブ6は、n個の超音波振動子T1、T2、．．．、Ti、．．．、Tnのアレイによって構成され、ここでnは、1以上の整数である。例として、プローブ6は、例えば軸Xに垂直である軸Yに沿って一直線になったn=128個の振動子を含む振動子の直線状ストリップの形であってもよい。しかしながら、問題になっているプローブは同等に、振動子の二次元アレイ(平面またはその他)の形である可能性もある。
振動子T1〜Tnは、おそらくは例えば柔軟なケーブルを介してプローブ6に接続される電子ラック7に含有される中央ユニットCPUを介して、マイクロコンピュータ4によって互いから独立して制御される。振動子T1〜Tnはそれ故に、
- 「平面」である超音波圧縮波(すなわち、その波面がX、Y平面で直線に進む波)、または媒質2中の全観測領域を照らす任意の他の種類の集束もしくは非集束波、例えばランダム音響信号を様々な振動子T1〜Tnによって放出させることによって生成される波か、
- さもなければ媒質2の1つもしくは複数の点に集束される超音波圧縮波を選択的に放出することができる。
- 「平面」である超音波圧縮波(すなわち、その波面がX、Y平面で直線に進む波)、または媒質2中の全観測領域を照らす任意の他の種類の集束もしくは非集束波、例えばランダム音響信号を様々な振動子T1〜Tnによって放出させることによって生成される波か、
- さもなければ媒質2の1つもしくは複数の点に集束される超音波圧縮波を選択的に放出することができる。
媒質2中のせん断波の伝搬を観測するために、いくつかのステップ、すなわち、
a)その間にマイクロコンピュータ4が、粘弾性媒質中で集束される少なくとも1つの超音波をプローブ6によって放出させることによって弾性せん断波を粘弾性媒質2中に生成する励起ステップ、
b)その間にせん断波の伝搬が、粘弾性媒質2中の多数の観測領域点で同時に観測される観測ステップ、
c)およびその間にマイクロコンピュータ4が、連続する瞬間tk(tkは、せん断波の生成から数えられる相対時間であってもよい)での連続伝搬画像を決定するためにサブステップb2)の間に粘弾性媒質2から受け取る音響信号を処理する撮像ステップが、連続して行われる。
a)その間にマイクロコンピュータ4が、粘弾性媒質中で集束される少なくとも1つの超音波をプローブ6によって放出させることによって弾性せん断波を粘弾性媒質2中に生成する励起ステップ、
b)その間にせん断波の伝搬が、粘弾性媒質2中の多数の観測領域点で同時に観測される観測ステップ、
c)およびその間にマイクロコンピュータ4が、連続する瞬間tk(tkは、せん断波の生成から数えられる相対時間であってもよい)での連続伝搬画像を決定するためにサブステップb2)の間に粘弾性媒質2から受け取る音響信号を処理する撮像ステップが、連続して行われる。
a)励起ステップ
励起ステップa)の間に、せん断励起は、超音波振動子のアレイ6によって粘弾性媒質2中に放出される少なくとも1つの集束超音波によって生じさせることができる。
励起ステップa)の間に、せん断励起は、超音波振動子のアレイ6によって粘弾性媒質2中に放出される少なくとも1つの集束超音波によって生じさせることができる。
励起ステップa)の間に放出される集束超音波は、0.5MHzから15MHzの範囲にある、例えば約4MHzに等しい周波数fの単色波であってもよく、それは、p/f秒の間放出され、ここでpは、50から5000の範囲にある(例えば、約500である)整数であり、fは、Hz単位で表される。矩形信号または同類の形のそのような励起は、少なくとも500Hzの比較的大きい帯域幅(例えば、0〜1000Hzの帯域幅)を持つせん断波を生成することができる。
励起ステップa)の間に放出される集束超音波は、生成されるようなせん断波が、所望の波形を提示しかつ媒質2中の所望の領域を照らすように、単一集束点または複数の集束点に集束されてもよい。集束点は、せん断波が生成される励起場所を構成する。例えば、励起場所Lは、いくつかの連続する超音波を直線に沿って素早く集束することによって得られるこの直線の形であってもよく、その場合には例えばBercoffらによって説明されるように[「Supersonic Shear Imaging: a New Technique for Soft Tissue Elasticity Mapping」、IEEE transactions on ultrasonics, ferroelectrics, and frequency control、Vol. 51、no. 4、2004年4月、396〜409頁]、平面であるせん断波を生成することが可能である。
b)観測ステップ
観測ステップは、次のサブステップを含んでもよく、
b1)マイクロコンピュータ4は、プローブ6に一連の超音波圧縮波衝撃を粘弾性媒質中に毎秒少なくとも500ショットの割合で放出させ、それは好ましくは、集束されないまたは軽度に集束されてもよく[ステップa)で放出される超音波の集束、および前記超音波のタイミングは、前記超音波の少なくともいくつかが、超音波放出の少なくともいくつかについて、観測領域を通るせん断波の伝搬の間に観測領域に達するように適合され]、
b2)マイクロコンピュータ4は、プローブ6に粘弾性媒質2から受け取る音響信号を実時間で検出させかつ記録させ、前記信号は、粘弾性媒質中の反射粒子5と相互作用する超音波圧縮波によって生成されるエコーを含み、これらのエコーは、粘弾性媒質の変位の連続画像に(直接的にまたは間接的に)対応する。
観測ステップは、次のサブステップを含んでもよく、
b1)マイクロコンピュータ4は、プローブ6に一連の超音波圧縮波衝撃を粘弾性媒質中に毎秒少なくとも500ショットの割合で放出させ、それは好ましくは、集束されないまたは軽度に集束されてもよく[ステップa)で放出される超音波の集束、および前記超音波のタイミングは、前記超音波の少なくともいくつかが、超音波放出の少なくともいくつかについて、観測領域を通るせん断波の伝搬の間に観測領域に達するように適合され]、
b2)マイクロコンピュータ4は、プローブ6に粘弾性媒質2から受け取る音響信号を実時間で検出させかつ記録させ、前記信号は、粘弾性媒質中の反射粒子5と相互作用する超音波圧縮波によって生成されるエコーを含み、これらのエコーは、粘弾性媒質の変位の連続画像に(直接的にまたは間接的に)対応する。
例えば1秒未満の間続くこともあるステップb1)の間に、例えば毎秒500から10,000ショットの範囲、好ましくは毎秒1000から5000ショットの範囲にある割合で非集束超音波圧縮波を放出することが可能である(この割合が、媒質2を通る圧縮波についての往復の移動時間によって、すなわち方向Xでの媒質2の厚さによって制限される場合、圧縮波によって生成されるエコーはすべて、新しい圧縮波が送られる前にプローブ6によって受け取られていることが必要である)。
各非集束超音波圧縮波は、せん断波の伝搬速度よりもはるかに速い伝搬速度(例えば、人体中で約1500メートル毎秒(m/s))で媒質2を通って伝播し、反射粒子5と相互作用し、それによって超音波検査の分野では名称「スペックル雑音」の下でそれ自体周知であるエコーまたは他の類似のかく乱を信号中に生成する。
スペックル雑音は、非集束超音波圧縮波の各ショットkの後、サブステップb2)の間に振動子T1〜Tnによって拾い上げられる。ショット番号kの後に各振動子Tiによってこのようにして拾い上げられるような信号si(tk)は、最初に高周波数(例えば、30MHzから100MHz)でサンプリングされ、ラック7の一部を形成しかつ前記振動子に接続されるサンプラーによって実時間でデジタル化(例えば、12ビットに)され、そのサンプラーは、それぞれE1、E2、．．．、Enと参照される。
このようにしてサンプリングされかつデジタル化されるような信号si(tk)は次いで、同様に実時間で、ラック7に属しかつ振動子Tiに特有のメモリMiに保存される。
例として、各メモリMiは、約128メガバイト(MB)の容量を提示し、ショットk=1からqにわたって連続して受け取る信号si(tk)をすべて含有し、ここでqは、超音波ショットの総数である。
c)撮像ステップ
撮像ステップc)は、せん断波の同じ伝搬に対応する信号si(tk)がすべて保存された後、例えば延期された時間にマイクロコンピュータ4によって行うことができ、中央ユニットCPUは、例えば米国特許第7252004(B2)号で説明されるように、サブステップc1)に対応する従来の経路形成(path-forming)ステップを使用して、これらの信号をラック7に属する加算回路Sによって再処理させる(さもなければ中央ユニットCPUは、この処理をそれ自身で行う、または実際その処理は、マイクロコンピュータ4で行われてもよい)。
撮像ステップc)は、せん断波の同じ伝搬に対応する信号si(tk)がすべて保存された後、例えば延期された時間にマイクロコンピュータ4によって行うことができ、中央ユニットCPUは、例えば米国特許第7252004(B2)号で説明されるように、サブステップc1)に対応する従来の経路形成(path-forming)ステップを使用して、これらの信号をラック7に属する加算回路Sによって再処理させる(さもなければ中央ユニットCPUは、この処理をそれ自身で行う、または実際その処理は、マイクロコンピュータ4で行われてもよい)。
これは、画像Im(tk)を、例えば今の場合には2D画像を生成し、各々は、観測領域での座標(x,y)の離散位置に対応するマトリクスの形であり、各々は、ショット番号kの後の観測領域の画像に対応する。
経路形成ステップの後、中央ユニットCPUは、ラック7の一部を形成する中央メモリMまたはマイクロコンピュータ4に画像Im(tk)を保存する。
各画像Im(tk)は、約1mm以下の分解能Rを有してもよい。
ステップa〜cの繰り返し
本発明によると、ステップa〜cは、例えば4と10との間に含まれてもよい回数Jにわたって繰り返される。励起ステップj=1からJの間に、励起場所Ljは、撮像デバイス1によって(振動子アレイ6によって伝送される超音波の集束を変えることによって、すなわち当技術分野で周知であるように、振動子T1〜Tnに加えられる遅延を変えることによって)繰り返しjごとにオフセットされる。連続する励起場所Ljは、せいぜい前記分解能Rに等しく、好ましくは前記分解能Rよりも小さい最大距離Dmだけ互いに分離される。特定の実施形態では、各励起場所Ljは、せいぜいR/5であり、好ましくはせいぜいR/10であるピッチだけ隣接励起場所から分離できる。例えば、人は、J=7の繰り返しを使用することができ、励起場所Ljは、各繰り返しでR/10のピッチだけオフセットでき、その結果最初の繰り返し(j=1)から最後(j=J)の間の全オフセットDmは、7R/10に等しい。
本発明によると、ステップa〜cは、例えば4と10との間に含まれてもよい回数Jにわたって繰り返される。励起ステップj=1からJの間に、励起場所Ljは、撮像デバイス1によって(振動子アレイ6によって伝送される超音波の集束を変えることによって、すなわち当技術分野で周知であるように、振動子T1〜Tnに加えられる遅延を変えることによって)繰り返しjごとにオフセットされる。連続する励起場所Ljは、せいぜい前記分解能Rに等しく、好ましくは前記分解能Rよりも小さい最大距離Dmだけ互いに分離される。特定の実施形態では、各励起場所Ljは、せいぜいR/5であり、好ましくはせいぜいR/10であるピッチだけ隣接励起場所から分離できる。例えば、人は、J=7の繰り返しを使用することができ、励起場所Ljは、各繰り返しでR/10のピッチだけオフセットでき、その結果最初の繰り返し(j=1)から最後(j=J)の間の全オフセットDmは、7R/10に等しい。
励起場所Ljが、例えば図2で示されるように軸Xに平行な直線であるときは、前記励起場所Ljは、繰り返しごとに軸Yに平行な前記ピッチだけオフセットできる。上記のオフセットは、励起場所Ljが励起ステップa)ごとに撮像デバイス1によって粘弾性媒質2中で移動されるという事実によって得ることができる。励起場所が直線であるときは、この線はまた、この線を小角度(数度)だけ傾けることによって繰り返しごとにオフセットでき、この場合には、1つの励起線の各点から別の励起線までの距離(前記距離はそのとき、前記点が属する励起線に垂直に測定される)は、前記線に沿って変化し、前記最大距離Dmは、1つの励起場所の1つの点からもう一方の励起場所までの前記距離の最大値である。
変形形態では、粘弾性媒質2は、アレイ6が固定されたままであり、励起場所が振動子アレイ6に対して固定されたままである間に励起ステップごとに移動され、その結果励起場所Ljは、励起ステップごとに粘弾性媒質2中で移動される。これは、例えば媒質2が、生きている人間または動物の心臓もしくは肝臓の一部であるときに当てはまることもある。そのような場合には、連続励起ステップa)は、時には媒質2の位置が基準位置に非常に近いところで実行されてもよく、その結果励起場所間の最大距離Dmは、分解能R未満のままである。
各励起ステップjの後、人は、連続する瞬間tkに上で述べられるような一組の画像Imj(tk)を得る。これらの画像Imj(tk)は、今後は原画像と呼ばれることになる。
d)再位置決めステップ
図3で示されるように、原画像Imj(tk)は、異なる励起場所(Lj)が様々な原画像Imj(tk)で相互に対応して位置決めされるように空間的に再位置決めされてもよい。この再位置決めは、例えば異なる励起場所(Lj)がすべての原画像Imj(tk)で同じ座標を有するように原画像Imj(tk)の座標系(x,y)[または3D撮像の場合には3D座標系(x,y,z)]をオフセットすることによって得ることができる。再位置決めステップは、各撮像ステップc)の後の各原画像についてか、またはステップa〜cのすべての繰り返しの後の同じ時間でのすべての原画像について実行できる。
図3で示されるように、原画像Imj(tk)は、異なる励起場所(Lj)が様々な原画像Imj(tk)で相互に対応して位置決めされるように空間的に再位置決めされてもよい。この再位置決めは、例えば異なる励起場所(Lj)がすべての原画像Imj(tk)で同じ座標を有するように原画像Imj(tk)の座標系(x,y)[または3D撮像の場合には3D座標系(x,y,z)]をオフセットすることによって得ることができる。再位置決めステップは、各撮像ステップc)の後の各原画像についてか、またはステップa〜cのすべての繰り返しの後の同じ時間でのすべての原画像について実行できる。
この再位置決めステップは、オプションとして省略されてもよい。粘弾性媒質2が、可動であり、励起場所が、上で説明されるように振動子アレイに対して固定されたままであるときは、この再位置決めステップは、好ましくは省略でき、反対に、励起場所Ljが、振動子アレイに対して励起ごとに移動されるときは、前記再位置決めは、好ましくは使用される。
e)平均化ステップ
同じ相対時間tkに対応する原画像Imj(tk)(もしあればそれらの再位置決め後の)は次いで、前記相対時間tkに対応する平均画像Im'(tk)を決定するために平均化される。平均は、原画像の画素値sj(tk)(x,y)の単純算術平均とすることができる。
同じ相対時間tkに対応する原画像Imj(tk)(もしあればそれらの再位置決め後の)は次いで、前記相対時間tkに対応する平均画像Im'(tk)を決定するために平均化される。平均は、原画像の画素値sj(tk)(x,y)の単純算術平均とすることができる。
図4で示されるように、平均化ステップは、図5で示されるそれのような原画像よりもはるかに良好な品質の平均画像を得ることを可能にする。特に、画像の信号対雑音比は、平均化ステップによって劇的に増大し、それは、次のステップでの粘弾性媒質2のレオロジー的特徴のより良好な特徴付けを可能にする。
押し付けている領域(pushing area)から2mmの横方向距離では、最大せん断波周波数は、200Hz(@-6dB)である。200Hzより上のすべての周波数成分は、組織によって減衰された。新しい平均化方法を使用するとき、最大せん断波周波数は、500Hz(@-6dB)であり、所与の場所について、従来の方法でよりも組織粘弾性についてはるかに多くの情報にアクセスできる。
f)特徴付けステップ
時間にわたる平均画像Im'(tk)に基づいて、少なくとも1つのせん断波伝搬パラメータが、粘弾性媒質の少なくとも1つの点で、好ましくは粘弾性媒質のマッピングを得るために完全な画像について、計算される。このせん断波伝搬パラメータは、例えばせん断波速度、せん断係数、ヤング率、せん断波減衰率、せん断弾性、せん断粘度、機械的緩和時間および異方性から選択できる。そのようなパラメータは、例えば米国特許第7252004(B2)号で説明されるように決定できる。
時間にわたる平均画像Im'(tk)に基づいて、少なくとも1つのせん断波伝搬パラメータが、粘弾性媒質の少なくとも1つの点で、好ましくは粘弾性媒質のマッピングを得るために完全な画像について、計算される。このせん断波伝搬パラメータは、例えばせん断波速度、せん断係数、ヤング率、せん断波減衰率、せん断弾性、せん断粘度、機械的緩和時間および異方性から選択できる。そのようなパラメータは、例えば米国特許第7252004(B2)号で説明されるように決定できる。
さらに、平均化画像の良好な信号対雑音比に起因して、前記パラメータの分光分析を実行することが可能である。例えば、図6で示されるように、前記パラメータが、せん断波の速度vであるときは、前記速度vをせん断波の帯域幅内で周波数fの関数として決定することが可能である。
測定されるパラメータが、せん断波減衰αであるときは、そのような減衰は、
α(f) = α0 + α1fa
などの法則に従って変化すると考えられてもよく、ここでfは、せん断波周波数であり、aは、生体媒質での機械的波について一般に0と2との間に含まれる、べき指数である。この場合には、本発明は、例えばWO-A-2009/007582で教示されるように、べき指数aを正確に評価し、粘弾性媒質中のyをマッピングすることを可能にする。
α(f) = α0 + α1fa
などの法則に従って変化すると考えられてもよく、ここでfは、せん断波周波数であり、aは、生体媒質での機械的波について一般に0と2との間に含まれる、べき指数である。この場合には、本発明は、例えばWO-A-2009/007582で教示されるように、べき指数aを正確に評価し、粘弾性媒質中のyをマッピングすることを可能にする。
変形形態
変形形態では、せん断波の励起は、例えばWO-A-00/55616で説明されるように、例えばマイクロコンピュータ4によって制御される振動器を通じて、外部の機械力を粘弾性媒質に各励起場所で加えることによって得られる可能性がある。
変形形態では、せん断波の励起は、例えばWO-A-00/55616で説明されるように、例えばマイクロコンピュータ4によって制御される振動器を通じて、外部の機械力を粘弾性媒質に各励起場所で加えることによって得られる可能性がある。
別の変形形態では、各原画像Imj(tk)は、超音波撮像によってよりもむしろMRIによって得られる。
1 撮像デバイス
2 粘弾性媒質
3 外表面
4 マイクロコンピュータ
4a 出力インターフェース
4b 入力インターフェース
5 拡散粒子
6 超音波プローブ、超音波振動子のアレイ
7 電子ラック
2 粘弾性媒質
3 外表面
4 マイクロコンピュータ
4a 出力インターフェース
4b 入力インターフェース
5 拡散粒子
6 超音波プローブ、超音波振動子のアレイ
7 電子ラック
Claims (19)
- 粘弾性媒質(2)を観測するためにせん断波を使用する撮像方法であって、
- いくつかの連続励起ステップjであって、該ステップの間に弾性せん断波が、撮像デバイス(1)によって前記粘弾性媒質(2)中の異なる励起場所(Lj)でそれぞれ生成され、前記異なる励起場所は、最大距離Dmだけ互いに分離される、いくつかの連続励起ステップjと、
- 時間tkでの前記粘弾性媒体の連続原画像Imj(tk)の組jが、前記粘弾性媒質中の前記せん断波の伝搬の間に前記撮像デバイスによって決定される、各励起ステップjに対応する撮像ステップであって、前記原画像は、前記最大距離Dmに少なくとも等しい分解能Rを有し、tkは、前記対応するせん断波の生成から数えられる相対時間である、撮像ステップと、
- 前記同じ相対時間tkに対応する原画像Imj(tk)が、前記相対時間tkに対応する平均画像Im'(tk)を決定するために平均化される平均化ステップとを含む、撮像方法。 - 前記撮像ステップと前記平均化ステップとの間に、前記原画像Imj(tk)が、前記異なる励起場所(Lj)が前記原画像Imj(tk)で相互に対応して位置決めされるように空間的に再位置決めされる再位置決めステップをさらに含む、請求項1に記載の撮像方法。
- 前記原画像Imj(tk)は、前記異なる励起場所(Lj)が前記原画像Imj(tk)で同じ座標を有するように各原画像の座標系をオフセットすることによって空間的に再位置決めされる、請求項2に記載の撮像方法。
- 連続原画像の前記組jの数Jは、少なくとも5である、請求項1から3のいずれか一項に記載の撮像方法。
- 各励起場所(Lj)は、せいぜいR/5であり、好ましくはせいぜいR/10であるピッチだけ隣接励起場所から分離される、請求項1から4のいずれか一項に記載の撮像方法。
- 前記励起場所(Lj)は、励起ステップごとに前記撮像デバイス(1)によって前記粘弾性媒質(2)中で移動される、請求項1から5のいずれか一項に記載の撮像方法。
- 前記粘弾性媒質(2)は、励起ステップごとに前記励起場所が前記撮像デバイス(1)に対して固定されたままである間に移動され、その結果前記励起場所(Lj)は、励起ステップごとに前記粘弾性媒質(2)中で移動される、請求項1から5のいずれか一項に記載の撮像方法。
- 前記せん断励起は、前記撮像デバイス(1)に属する超音波振動子のアレイ(6)によって前記粘弾性媒質(2)中に放出される少なくとも1つの超音波によって生成される、請求項1から7のいずれか一項に記載の撮像方法。
- 前記せん断励起は、外部の機械力を前記粘弾性媒質に各励起場所で加えることによって生成される、請求項1から7のいずれか一項に記載の撮像方法。
- 各原画像Imj(tk)は、超音波撮像によってか、またはMRIによって得られる、請求項1から9のいずれか一項に記載の撮像方法。
- その間に時間にわたる前記平均画像Im'(tk)に基づいて、少なくとも1つのせん断波伝搬パラメータが、前記粘弾性媒質の少なくとも1つの点で計算される特徴付けステップを含む、請求項1から10のいずれか一項に記載の撮像方法。
- 前記特徴付けステップの間に計算される前記せん断波伝搬パラメータは、せん断波速度、せん断係数、ヤング率、せん断波減衰率、せん断弾性、せん断粘度、機械的緩和時間および異方性から選択される、請求項11に記載の方法。
- 前記励起は、少なくとも500Hzの帯域幅を有し、前記せん断波伝搬パラメータは、前記帯域幅内の複数の周波数で計算される、請求項10または請求項11に記載の方法。
- 粘弾性媒質(2)を観測するためにせん断波を使用する請求項1から13のいずれかに記載の方法を実装するための撮像装置であって、前記装置は、
- いくつかの連続励起ステップjの間に前記粘弾性媒質(2)中の異なる励起場所(Lj)でそれぞれ弾性せん断波を生成するように適合され、前記異なる励起場所は、最大距離Dmだけ互いに分離され、
- 各励起ステップjについてそれぞれ前記粘弾性媒質中の前記せん断波の伝搬の間に時間tkでの前記粘弾性媒質の連続原画像Imj(tk)の組jを決定するように適合され、前記原画像は、前記最大距離Dmに少なくとも等しい分解能Rを有し、tkは、前記対応するせん断波の生成から数えられる相対時間であり、
- 前記相対時間tkについて平均画像Im'(tk)を決定するために前記同じ相対時間tkに対応する原画像Imj(tk)を平均化するように適合される
少なくとも1つの電子中央ユニット(4)を含む、撮像装置。 - 前記電子中央ユニット(4)はさらに、前記異なる励起場所(Lj)が前記原画像Imj(tk)で相互に対応して位置決めされるように、前記画像を平均化する前に前記原画像Imj(tk)を空間的に再位置決めするように適合される、請求項14に記載の装置。
- 前記電子中央ユニット(4)は、前記異なる励起場所(Lj)が前記原画像Imj(tk)で同じ座標を有するように各原画像の座標系をオフセットさせることによって前記原画像Imj(tk)を空間的に再位置決めするように適合される、請求項15に記載の装置。
- 前記電子中央ユニット(4)は、励起ステップごとに前記粘弾性媒質(2)中の前記励起場所(Lj)を移動させるように適合される、請求項14から16のいずれか一項に記載の装置。
- 前記電子中央ユニット(4)は、時間にわたる連続平均画像に基づいて前記粘弾性媒質(2)の少なくとも1つの点で少なくとも1つのせん断波伝搬パラメータを計算するように適合され、前記せん断波伝搬パラメータは、せん断波速度、せん断係数、ヤング率、せん断波減衰率、せん断弾性、せん断粘度、機械的緩和時間および異方性から選択される、請求項14から17のいずれか一項に記載の装置。
- 前記装置は、少なくとも500Hzの帯域幅を持つ前記励起を加えるように適合され、前記中央ユニット(4)は、前記帯域幅での複数の周波数で前記せん断波伝搬パラメータを計算するように適合される、請求項18に記載の装置。
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/IB2010/001282 WO2011132014A1 (en) | 2010-04-20 | 2010-04-20 | Imaging method and apparatus using shear waves |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2013529098A true JP2013529098A (ja) | 2013-07-18 |
| JP5760080B2 JP5760080B2 (ja) | 2015-08-05 |
Family
ID=43416784
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2013505551A Active JP5760080B2 (ja) | 2010-04-20 | 2010-04-20 | せん断波を使用する撮像方法および装置 |
Country Status (9)
| Country | Link |
|---|---|
| US (1) | US9140781B2 (ja) |
| EP (1) | EP2561380B1 (ja) |
| JP (1) | JP5760080B2 (ja) |
| KR (1) | KR101581369B1 (ja) |
| CN (1) | CN103026257B (ja) |
| CA (1) | CA2797262C (ja) |
| HK (1) | HK1178607A1 (ja) |
| IL (1) | IL222267A (ja) |
| WO (1) | WO2011132014A1 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20190028268A (ko) * | 2017-09-08 | 2019-03-18 | 삼성메디슨 주식회사 | 초음파 영상장치 및 그 제어방법 |
Families Citing this family (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2012364909B2 (en) | 2011-02-15 | 2016-11-24 | Hemosonics, Llc | Characterization of blood hemostasis and oxygen transport parameters |
| US20120253194A1 (en) * | 2011-03-30 | 2012-10-04 | Tadashi Tamura | Methods and apparatus for ultrasound imaging |
| EP2599444B1 (en) | 2011-12-01 | 2017-08-02 | Theraclion SA | Method and device for determining the elastic modulus of a biological tissue |
| CN103800038B (zh) * | 2012-11-12 | 2016-09-21 | 通用电气公司 | 改善的系统和装置以用于确定目标组织的机械特性 |
| CN103064047B (zh) * | 2012-12-25 | 2014-12-10 | 深圳先进技术研究院 | 磁共振弹性成像的弹性重建修正方法和系统 |
| PL2959312T3 (pl) * | 2013-02-19 | 2018-05-30 | Echosens | Metoda elastografii wieloimpulsowej |
| EP2958496B1 (en) * | 2013-02-25 | 2017-05-03 | Koninklijke Philips N.V. | Determination of the concentration distribution of sonically dispersive elements |
| KR20150014315A (ko) | 2013-07-29 | 2015-02-06 | 삼성전자주식회사 | 횡파를 이용한 초음파 진단 방법 및 장치 |
| US9332963B2 (en) | 2014-01-21 | 2016-05-10 | Siemens Medical Solutions Usa, Inc. | Swept focus for acoustic radiation force impulse |
| CN107205724B (zh) * | 2014-10-28 | 2020-10-30 | 超声成像 | 成像方法以及进行剪切波弹性成像的装置 |
| US11058397B2 (en) * | 2015-01-02 | 2021-07-13 | Esaote S.P.A. | Method for quantifying the elasticity of a material by ultrasounds |
| US9726647B2 (en) | 2015-03-17 | 2017-08-08 | Hemosonics, Llc | Determining mechanical properties via ultrasound-induced resonance |
| CN104751037B (zh) * | 2015-04-10 | 2018-06-12 | 无锡海斯凯尔医学技术有限公司 | 医疗检测设备的使用控制方法、系统和医疗检测设备 |
| US10376242B2 (en) | 2015-04-16 | 2019-08-13 | Siemens Medical Solutions Usa, Inc. | Quantitative viscoelastic ultrasound imaging |
| US10495159B2 (en) * | 2015-04-24 | 2019-12-03 | Magna Powertrain Inc. | Electronically-controlled selectable clutch assembly with clutch module having integrated safety switch |
| CN107970043B (zh) * | 2017-12-28 | 2021-01-19 | 深圳开立生物医疗科技股份有限公司 | 一种剪切波的检测方法及装置 |
| JP7284769B2 (ja) * | 2018-06-27 | 2023-05-31 | コーニンクレッカ フィリップス エヌ ヴェ | 解剖学的粘度の剪断波検出 、関連するデバイス、システム、及び方法 |
| CN111077227B (zh) * | 2019-12-18 | 2021-11-02 | 华南理工大学 | 一种超声阵列扫查反演方法、系统、存储介质及设备 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002538911A (ja) * | 1999-03-15 | 2002-11-19 | ソシエテ・デラストグラフィ・アンピュルショネル・プール・レ・システム・ドゥ・メジュール・ドゥ・レラスティシテ（エスウイエスエムウ） | 剪断波を用いたイメージング方法および装置 |
| JP2003530941A (ja) * | 2000-04-26 | 2003-10-21 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 剪断波パラメター評価のための超音波方法およびシステム |
| US7444875B1 (en) * | 2005-02-18 | 2008-11-04 | University Of Rochester | Real time visualization of shear wave propagation in soft materials with sonoelastography |
| WO2008139245A1 (en) * | 2007-05-16 | 2008-11-20 | Super Sonic Imagine | Method and device for measuring a mean value of visco-elasticity of a region of interest |
| JP2009534198A (ja) * | 2006-04-21 | 2009-09-24 | レニショウ パブリック リミテッド カンパニー | 誤差補正の方法 |
| WO2010012092A1 (en) * | 2008-07-30 | 2010-02-04 | Centre Hospitalier De L'universite De Montreal | A system and method for detection, characterization and imaging of heterogeneity using shear wave induced resonance |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003235839A (ja) | 2002-02-18 | 2003-08-26 | Matsushita Electric Ind Co Ltd | 超音波診断装置 |
| FR2844058B1 (fr) | 2002-09-02 | 2004-11-12 | Centre Nat Rech Scient | Procede et dispositif d'imagerie utilisant des ondes de cisaillement |
| US7175599B2 (en) * | 2003-04-17 | 2007-02-13 | Brigham And Women's Hospital, Inc. | Shear mode diagnostic ultrasound |
| CN100562293C (zh) * | 2003-04-17 | 2009-11-25 | 布赖汉姆妇女医院 | 剪切式诊断超声波系统 |
| FR2913875B1 (fr) | 2007-03-21 | 2009-08-07 | Echosens Sa | Dispositif pour mesurer des proprietes viscoelastiques de tissus biologiques et procede utilisant ce dispositif |
| FR2917831B1 (fr) | 2007-06-25 | 2009-10-30 | Super Sonic Imagine Sa | Procede de caracterisation rheologique d'un milieu viscoelastique |
| US8187187B2 (en) * | 2008-07-16 | 2012-05-29 | Siemens Medical Solutions Usa, Inc. | Shear wave imaging |
| CN201370595Y (zh) | 2009-02-20 | 2009-12-30 | 刘丹宁 | 一种肝脏多维超声弹性检测装置 |
| US8147410B2 (en) * | 2009-03-23 | 2012-04-03 | The Hong Kong Polytechnic University | Method and apparatus for ultrasound imaging and elasticity measurement |
| CN101690677B (zh) | 2009-10-13 | 2011-09-21 | 西安交通大学 | 聚焦超声治疗组合阵元相控阵及多焦点剪切波成像系统 |
| US8726734B1 (en) * | 2010-09-15 | 2014-05-20 | Sonowise, Inc. | Shear wave generation system and methods for ultrasound imaging |
-
2010
- 2010-04-20 EP EP20100759969 patent/EP2561380B1/en active Active
- 2010-04-20 US US13/642,478 patent/US9140781B2/en active Active
- 2010-04-20 WO PCT/IB2010/001282 patent/WO2011132014A1/en active Application Filing
- 2010-04-20 KR KR1020127029983A patent/KR101581369B1/ko active IP Right Grant
- 2010-04-20 CN CN201080067697.4A patent/CN103026257B/zh active Active
- 2010-04-20 JP JP2013505551A patent/JP5760080B2/ja active Active
- 2010-04-20 CA CA2797262A patent/CA2797262C/en not_active Expired - Fee Related
-
2012
- 2012-10-09 IL IL222267A patent/IL222267A/en active IP Right Grant
-
2013
- 2013-05-09 HK HK13105567.7A patent/HK1178607A1/xx not_active IP Right Cessation
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002538911A (ja) * | 1999-03-15 | 2002-11-19 | ソシエテ・デラストグラフィ・アンピュルショネル・プール・レ・システム・ドゥ・メジュール・ドゥ・レラスティシテ（エスウイエスエムウ） | 剪断波を用いたイメージング方法および装置 |
| JP2003530941A (ja) * | 2000-04-26 | 2003-10-21 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 剪断波パラメター評価のための超音波方法およびシステム |
| US7444875B1 (en) * | 2005-02-18 | 2008-11-04 | University Of Rochester | Real time visualization of shear wave propagation in soft materials with sonoelastography |
| JP2009534198A (ja) * | 2006-04-21 | 2009-09-24 | レニショウ パブリック リミテッド カンパニー | 誤差補正の方法 |
| WO2008139245A1 (en) * | 2007-05-16 | 2008-11-20 | Super Sonic Imagine | Method and device for measuring a mean value of visco-elasticity of a region of interest |
| WO2010012092A1 (en) * | 2008-07-30 | 2010-02-04 | Centre Hospitalier De L'universite De Montreal | A system and method for detection, characterization and imaging of heterogeneity using shear wave induced resonance |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20190028268A (ko) * | 2017-09-08 | 2019-03-18 | 삼성메디슨 주식회사 | 초음파 영상장치 및 그 제어방법 |
| KR102627723B1 (ko) | 2017-09-08 | 2024-01-23 | 삼성메디슨 주식회사 | 초음파 영상장치 및 그 제어방법 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103026257A (zh) | 2013-04-03 |
| US9140781B2 (en) | 2015-09-22 |
| CA2797262C (en) | 2016-12-13 |
| KR101581369B1 (ko) | 2015-12-30 |
| KR20130057435A (ko) | 2013-05-31 |
| IL222267A0 (en) | 2012-12-31 |
| IL222267A (en) | 2016-08-31 |
| CA2797262A1 (en) | 2011-10-27 |
| EP2561380B1 (en) | 2014-02-26 |
| JP5760080B2 (ja) | 2015-08-05 |
| CN103026257B (zh) | 2015-05-20 |
| WO2011132014A1 (en) | 2011-10-27 |
| HK1178607A1 (en) | 2013-09-13 |
| EP2561380A1 (en) | 2013-02-27 |
| US20130031981A1 (en) | 2013-02-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5760080B2 (ja) | せん断波を使用する撮像方法および装置 | |
| US11672509B2 (en) | Shear wave elastrography method and apparatus for imaging an anisotropic medium | |
| JP6053860B2 (ja) | 医療用超音波イメージングシステムにおいて剪断波情報を求めるためのコンピュータプログラム | |
| KR101501479B1 (ko) | 초음파 최적화 방법과 그를 위한 초음파 의료 장치 | |
| JP5530685B2 (ja) | 剛性が変化した領域を検出するためのシステム及び方法 | |
| CN109077754B (zh) | 一种测量组织力学特性参数的方法及设备 | |
| CN103505243A (zh) | 测量超声波的声吸收或衰减 | |
| CN102724917A (zh) | 通过传播剪切波来测量哺乳动物软组织物理参数的方法和装置 | |
| JP2012081269A5 (ja) | ||
| CA2218682A1 (en) | Method and apparatus for osteoporosis diagnosis | |
| Scorza et al. | A comparative study on the influence of probe placement on quality assurance measurements in B-mode ultrasound by means of ultrasound phantoms | |
| JP2008136880A (ja) | 超音波診断システム、歪み分布表示方法及び弾性係数分布表示方法 | |
| CN109758091A (zh) | 一种超声成像方法及装置 | |
| JP2012200443A (ja) | 超音波診断装置および超音波画像生成方法 | |
| JP5836241B2 (ja) | 超音波検査装置、超音波検査装置の信号処理方法およびプログラム | |
| JP2023540954A (ja) | 媒体の超音波特性評価用の方法及びシステム | |
| Fiori et al. | A novel method for the gain conversion factor estimation in quality assessment of ultrasound diagnostic systems | |
| WO2020113381A1 (zh) | 一种超声成像方法及装置 | |
| Kofler et al. | Optimization of medical ultrasound transducers with simulation and genetic algorithms | |
| JP2012249925A (ja) | 超音波画像生成方法および超音波画像診断装置 | |
| Wüst et al. | Towards a Model-based High-resolution Ultrasonic Measurement System for Non-destructive Testing | |
| JP2014124231A (ja) | 超音波検査装置、超音波検査装置の信号処理方法およびプログラム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20140307 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140714 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20141003 |
|
| A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20141010 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20141212 |
|
| 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: 20150511 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20150608 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5760080 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 |
|
| S633 | Written request for registration of reclamation of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313633 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313115 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |