JP2016097256A - 超音波画像処理装置 - Google Patents
超音波画像処理装置 Download PDFInfo
- Publication number
- JP2016097256A JP2016097256A JP2014239090A JP2014239090A JP2016097256A JP 2016097256 A JP2016097256 A JP 2016097256A JP 2014239090 A JP2014239090 A JP 2014239090A JP 2014239090 A JP2014239090 A JP 2014239090A JP 2016097256 A JP2016097256 A JP 2016097256A
- Authority
- JP
- Japan
- Prior art keywords
- frame
- measurement
- time phase
- cycle
- unit
- 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.)
- Pending
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
【課題】検査者に格別な負担を要せずに、計測に用いるための好適なフレームセットをフレーム列の中から選択する。【解決手段】フレーム列解析部32は、計測対象組織である胎児の心臓から取得されたフレーム列に含まれる各フレームの輝度情報に基づいて、心臓の心拍周期を示す心拍周期情報を取得する。フレーム特定部34は、心拍周期情報に基づいてフレーム列の中から拡張末期フレーム及び収縮末期フレームからなるフレームセットを特定する。表示処理部22は、特定されたフレームセットを並べて表示部24に表示させる。最初に表示されたフレームセットが計測用として適切でない場合、表示処理部22は、ユーザの操作に応じて、心拍周期単位でフレームセットの切り替え表示を行う。【選択図】図1
Description
本発明は超音波画像処理装置に関し、特に、超音波の送受波により取得したフレーム列の表示処理に関する。
超音波診断装置は、被検体に対して超音波を送受波し、これにより得られたデータに基づいて超音波画像を形成する装置である。例えばBモードにおいては、対象組織に対して超音波が送受波されることでフレームデータ列が取得される。取得されたフレームデータ列は、画像メモリ(シネメモリ)に一時記憶されるとともに、スキャンコンバータなどにより処理されることでBモード動画像としてリアルタイムに表示される。計測、印刷等の必要から、ユーザによりフリーズ操作が行われると、動画像のリアルタイム表示が停止し、フリーズ操作時点において表示されていた画像が静止画として表示される。その後、必要に応じて、画像メモリ上に記憶されたフレームデータ列が再生される。ユーザによりストア操作がされると、画像メモリに一時記憶されているフレームデータ列又は特定のフレームデータが不揮発性の記憶装置に記憶される。なお、本明細書においては、「フレーム」の用語は、フレームデータの簡略表現であり、その概念には、超音波ビームが電子走査されることで得られたフレームデータ、及びスキャンコンバータによる処理後のフレームデータが含まれる。
超音波画像は各種計測に用いられる。計測は、フリーズ操作によって表示された静止画像、若しくはシネメモリに記憶されたフレーム列を再生して得られる静止画像又は動画像を用いて行われる。超音波診断装置からフレーム列を受け取った情報処理装置において計測が実行される場合もある。計測の種類としては、Bモード画像上において所定の2点間の距離を測定し、当該測定値から対象組織内の所定の2点間の距離を算出する距離計測、Bモード画像上の所定範囲の面積から対象組織の断面積を計測する面積計測などがある。
例えば心臓のように、周期的に運動する組織について計測を行う場合、運動周期内の複数の時相において計測を行い、得られた複数の計測結果の変化を見る場合がある。例えば、FS計測と呼ばれる計測は、心臓の左室内径の短縮率(FS値)を計測するものであり、心臓の拡張末期における左室内径と心臓の収縮末期における左室内径とを算出し、これらの計測結果に基づいてFS値が算出される。
上記のような計測では、同じ心拍内における拡張末期のフレームと収縮末期のフレームからなるフレームセットが特定され、それらが同時に画面上に表示される。すなわち、予め1又は複数の運動周期に亘って対象組織に超音波を送受波してフレーム列が取得され、取得されたフレーム列の中から計測に用いる、同じ心拍内の複数の所定時相に対応した複数のフレーム(フレームセット)が事後的に選択されるのが一般的である。
例えば上述のFS計測においては、心拍周期の拡張末期と収縮末期における左室内径を計測する必要がある。つまり、フレーム列の中から拡張末期において取得されたフレームと収縮末期において取得されたフレームとを適切に選択しなければならない。
従来、計測に用いるフレームセットの選択は、医師などのユーザにより手動で行われていた。例えば、FS計測においては、被検体の心臓に対して超音波を複数心拍に亘り送受波して得られたフレーム列を1枚1枚表示させながら、ユーザが目視により拡張末期において取得されたフレームと収縮末期において取得されたフレームを選択していた。
上記のような目視判断による選択方法は、検査者の負担が大きく、時間がかかってしまうという問題が指摘できる。また、このような選択方法は、ユーザの主観に依存するという面も指摘できる。つまり、同じフレーム列であっても、ユーザに応じて異なるフレームセットが選択されてしまうという問題が指摘できる。
本発明の目的は、検査者に格別な負担を要せずに、計測に用いるための好適なフレームセットをフレーム列の中から選択することにある。
本発明に係る超音波画像処理装置は、周期的に運動する対象組織に超音波を送受波することで得られたフレーム列を解析することにより、前記対象組織の運動周期を示す周期情報を取得する周期情報取得手段と、前記周期情報に基づいて、前記フレーム列の中から、前記運動周期中の第1時相に対応する第1時相フレーム及び前記運動周期中の第2時相に対応する第2時相フレームを含むフレームセットを特定するフレームセット特定手段と、
前記フレームセットを並べて表示部に表示させる表示制御手段と、を備えることを特徴とする。
前記フレームセットを並べて表示部に表示させる表示制御手段と、を備えることを特徴とする。
上記構成によれば、周期情報取得手段により、取得されたフレーム列が解析されることで周期情報が取得される。対象組織が周期的に運動しているため、対象組織に対して超音波を送受波して得られるフレーム列においても対象組織の運動周期に応じた周期的な変化が現れる。例えば、当該フレーム列に含まれる各フレームの輝度情報(例えばフレーム内の特定領域の平均輝度)は、対象組織の運動に応じて変化することから、周期情報取得手段は、当該輝度の変化に基づいて周期情報を取得することができる。当該周期情報の取得方法は、例えば心電計が利用できない胎児の心臓についてFS計測を行うなど、対象組織の周期情報が直接取得できない場合に特に有効である。
フレームセット特定手段により、周期情報に基づいて複数の特定時相に対応した複数のフレームからなるフレームセットがフレーム列から選択され、表示制御手段によりそれが表示される。つまり、実行しようとしている計測に応じて、計測対象候補となる所定時相関係にある複数のフレームが同時表示される。よって、それらを個別観察あるいは対比観察して、セット単位で内容の良否を一括して評価することができるから、その評価を的確に行うことができ、また、計測対象を特定するまでの時間を短縮できる。ユーザは、計測用のフレームセットを選択するために、フレーム列に含まれる各フレームを1つ1つ目視する必要がなくなる。つまり、ユーザは、フレーム列の中から計測用のフレームセットをより容易に選択することが可能になる。上記構成によれば、実行しようとしている計測に応じて計測対象候補の関連付け(ペアリング、グルーピング)が自動的に行われるから、その点でユーザの負担が大幅に軽減され、計測効率を高められる。
望ましくは、前記表示制御手段は、前記フレーム列を用いた計測の開始操作時点から過去に直近の運動周期に対応するフレームセットを最初に表示させる。ユーザは、フレーム列の取得中にリアルタイム表示される動画像を目視して、計測に適した画像が撮れたと思った瞬間にフリーズ操作を行って計測を開始するのが一般的である。すなわち、計測の開始操作時点から直近の運動周期に対応するフレームセットが計測に用いるフレームセットとして適切なものである可能性が高い。したがって、計測の開始操作時点から直近の運動周期に対応する前記フレームセットを最初に表示することで、より計測に適している可能性の高いフレームセットを優先的に表示することができる。これにより、ユーザの計測用フレームセットの選択に係る手間あるいは時間をより低減させることができる。
望ましくは、前記表示制御手段は、ユーザの指示に従って、複数の運動周期に対応する複数のフレームセットを切り替え表示させる。最初に表示されたフレームセットが計測用のフレームセットとして適切でないとユーザが判断した場合、ユーザは他のフレームセットを選択する必要がある。このとき、候補となるのは、現在表示されているフレームセットに対応する運動周期以外の運動周期に対応するフレームセットである。したがって、表示制御手段は、ユーザ操作に応じて、各運動周期に対応する複数のフレームセットを切り替え表示させる。つまり、計測に使わない複数のフレームをジャンプして、計測対象候補だけが順次表示される。これにより、最初に表示されたフレームセットが計測用として適切でないフレームセットである場合であって他のフレームセットを選択する場合に、次なる計測用候補のフレームセットの選択にかかるユーザの手間あるいは時間をより低減させることができる。
望ましくは、前記対象組織は胎児の心臓であり、前記運動周期は心拍周期である。また、望ましくは、前記第1時相は、前記心拍周期における拡張末期であり、前記第2時相は、前記心拍周期における収縮末期である。
本発明によれば、検査者に格別な負担を要せずに、計測に用いるための好適なフレームセットをフレーム列の中から選択することができる。
以下、本発明の実施形態について説明する。本実施形態では、超音波診断装置10において計測対象組織である胎児の心臓の断層画像(Bモード画像)が形成され、これを用いてFS計測が実行される。
図1は、本実施形態に係る超音波画像処理装置としての超音波診断装置10の構成概略図である。超音波診断装置10は、一般に病院などの医療機関に設置される医療上の機器である。
プローブ12は、被検体に対して超音波の送受波を行う超音波探触子である。プローブ12は複数の振動子からなる振動子アレイを有している。複数の振動子には、送信部14からの複数の送信信号が供給される。これにより超音波ビームが形成される。振動子アレイを構成する複数の振動素子は送受波領域からの反射エコーを受信し、これにより複数の受信信号が受信部16へ出力される。本実施形態では、プローブ12は母体の腹部に当接され、母体の腹壁を介して胎児の心臓に対して超音波の送受波が行われる。
送信部14は、制御部26からの制御に応じて、プローブ12が有する複数の振動子に対して、それらを励振する複数の送信信号を送る。入力部30においてユーザによりフリーズ操作が行われると、送信部14は、プローブ12への送信信号の供給を一時停止する。ユーザによりフリーズ解除操作が行われると再度送信信号の供給を行う。
受信部16は、プローブ12の複数の振動子から出力された複数の受信信号を整相加算処理してビームデータを形成する。超音波ビームの電子走査に伴って、複数のビームデータが順次出力される。1回の電子走査当たり、1フレーム分の複数のビームデータが得られる。各ビームデータは、深さ方向に並ぶ複数の反射エコー信号により構成される。電子走査が繰り返し行われることにより、受信フレーム列が生成される。
画像形成部18は、例えばデジタルスキャンコンバータ(DSC)であり、受信部16からの受信フレーム列を表示フレーム列に変換するものである。具体的には、画像形成部18は座標変換機能及び補間処理機能などを備えている。本実施形態では、画像形成部18は、表示部24にBモード画像として表示される表示フレーム列を形成する。なお、後述の画像メモリ20が画像形成部18の前段に設けられている場合には、画像形成部18は、画像メモリに記憶された受信フレーム列を表示フレーム列に変換してもよい。
画像形成部18の前段又は後段に設けられる揮発性の画像メモリ20には、受信部16が出力する受信フレーム列又は画像形成部18で処理された表示フレーム列が一時的に記憶される。画像メモリ20は例えばリングバッファ構造を有しており、時系列順に入力される各フレームを記憶する。すなわち、画像メモリ20には、現時点から過去一定期間において取得されたフレーム列が記憶される。入力部30においてユーザによりストア操作が行われると、画像メモリ20に記憶されているフレーム列が不揮発性の画像記憶部(不図示)に転送される。本実施形態では画像メモリ20は画像形成部18の後段に設けられており、画像メモリには画像形成部18による処理後の表示フレーム列が記憶される。
表示処理部22は、画像形成部18が形成した表示フレーム列を処理し、表示部24に超音波画像としてのBモード画像を表示させる。また、後述のフレーム特定部34が特定したフレームセットに対応する複数のBモード画像を表示部24に表示させる。また、後述の計測部36が出力する計測結果情報(例えば計測値データやグラフデータなど)を処理し表示部24に計測値やグラフを表示させる。超音波画像と計測結果を組み合わせるなどして表示するようにしてもよい。表示部24は、例えば液晶ディスプレイなどで構成される。
制御部26は、例えばCPU又はマイクロプロセッサなどのハードウェア、あるいはCPU又はマイクロプロセッサの動作を規定するソフトウェア(プログラム)との協働により実現することができる。制御部26は、超音波診断装置10の各構成要素の制御を行う。
記憶部28は、不揮発性であり、例えばハードディスク、ROM、又はRAMなどで構成される。また、記憶部28には、超音波診断装置10の各構成要素を動作させるためのプログラムが記憶される。さらに、記憶部28には、後述の計測部36の計測結果などが記憶されてもよい。
入力部30は、ボタン、スイッチ、ダイアル式つまみ、及びトラックボールなどを含み、ユーザの指示を制御部26へ入力するためのものである。ユーザは、入力部30を用いて、フリーズ操作、フリーズ解除操作、ストア操作、計測開始操作、あるいは計測点入力操作などの各種操作を行うことができる。
次に、計測のために機能する構成について説明する。その計測は、リアルタイム動作後のフリーズ状態で実行されるものである。フレーム列解析部32は、図示された構成例では、記憶部28に記憶されたフレーム列を解析する。これにより、周期的に運動する対象組織の運動周期を示す周期情報を取得する。つまり、フレーム列解析部32は周期情報取得手段としての機能を有する。本実施形態では、対象組織が胎児の心臓であるため、フレーム列解析部32は、周期情報として胎児の心拍周期情報を取得する。
具体的には、フレーム列解析部32は、対象組織である胎児の心臓の拍動に応じて生じる、フレーム列に含まれる各フレームの輝度情報(輝度値)の変化に基づいて心拍周期情報を取得する。輝度情報としては、各フレームが有する各画素の輝度値の和あるいは平均を用いることができる。各フレームから得られた輝度情報は、胎児の心拍に応じた周期性を有するため、これに基づいて胎児の心拍周期情報を得ることができる。また、心拍周期情報を精度よく取得するために、フレーム内の所定の領域における輝度値を輝度情報として用いてもよい。フレーム列解析部32の処理の詳細については、図2及び図3を用いて後述する。
なお、周期情報の取得方法は、各フレームの輝度値の変化を解析する手法に限られない。例えば、フレーム列に含まれる各フレームについて対象組織の輪郭抽出処理を行い、抽出された輪郭の形状あるいは大きさなどに基づいて周期情報を取得する手法などであってもよい。また、計測対象が成人の心臓である場合などは、心電計で直接心拍情報を取得するようにしてもよい。
フレーム特定部34は、フレーム列解析部32が取得した心拍周期情報に基づいて、フレーム列を構成する複数のフレームの中から、対象組織の心拍周期における複数の時相において取得された複数のフレーム(フレームセット)を特定する。つまり、フレーム特定部34は、フレームセット特定手段としての機能を有する。フレームセットに含まれる各フレームは、同心拍周期内から選択される。具体的には、計測内容に応じて、要請される時相関係を満たす複数のフレームがフレームセットとして特定される。
特定時相とは、計測に必要なフレームが取得される時相である。したがって、特定時相は計測内容に応じて異なる。フレーム特定部34が、どの時相のフレームを特定するかは、予めユーザに指定されるか、あるいは計測内容に応じて自動的に決定されてよい。本実施形態では、FS計測を行うため、フレーム特定部34は、心拍周期における拡張末期において取得されたフレーム及び収縮末期において取得されたフレームを含むフレームセットを特定する。あるいは、これに加え、他の時相におけるフレーム(例えば拡張末期と収縮末期の中間時点において取得されたフレーム)などをフレームセットに含めるようにしてもよい。フレーム特定部34の処理の詳細については、図4を用いて後述する。
計測部36は、表示処理部22により表示部24に表示される超音波画像上において各種計測を行う。計測部36が実行可能な計測としては、対象組織の所定の2点間の距離を計測する距離計測、対象組織の所定領域の面積を計測する面積計測などがある。超音波画像上の計測位置(計測点)はユーザによって指定される。計測部36は、フレーム特定部34が特定したフレームセットにおいてこれらの計測を行い、複数の計測値を得る。そして、複数の計測値に基づいて、対象組織の運動機能を評価する値を算出する。例えば、上述のFS計測の他、心臓の拡張末期における左室容積と収縮末期における左室容積とに基づいて算出される左室の収縮機能の指標となるEF(Ejection Fraction)値などが算出される。EF値の算出においては、拡張末期及び収縮末期において取得されたフレームにおいて上述の面積計測により左室面積を測定し、それに基づいてそれぞれの左室容積が算出される。また、右室自由壁側の三尖弁輪の移動距離を計測するTAPSE(Tricuspid Annular Plane Systolic Excursion)計測なども実行可能である。TAPSE計測においては、拡張末期において取得されたフレームにおける三尖弁輪の位置と収縮末期において取得されたフレームにおける三尖弁輪の位置との距離に基づいて目的の距離が計測される。
対象組織は周期的な運動を行うため、フレーム列内の各フレームにおいて対象組織の注目部位の位置がそれぞれ異なる。トラッキング処理部38は、フレームごとに変化する注目部位に追従して、当該注目部位上に設定された計測点の位置をフレームごとに変化させる(すなわち計測点の位置を注目部位上に維持する)トラッキング処理を行う。トラッキングの手法としては、例えばパターンマッチング技術など、既知の技術を利用することができる。
本実施形態では、超音波画像処理装置として超音波診断装置10が用いられているが、超音波画像処理装置としては、超音波診断装置10以外のコンピュータ、例えばパーソナルコンピュータ(PC)を用いることが可能である。この場合は、当該コンピュータは、超音波診断装置で形成されたフレーム列を受け取って上記の処理を行う。
なお、図1に示す各構成(符号を付した各部)のうち、送信部14、受信部16、画像形成部18、表示処理部22、フレーム列解析部32、フレーム特定部34、計測部36、及びトラッキング処理部38の各部は、例えば電気電子回路やプロセッサ等のハードウェアを利用して実現することができ、その実現において必要に応じてメモリ等のデバイスが利用されてもよい。また、上記各部に対応した機能が、CPUやプロセッサやメモリなどのハードウェアと、CPUやプロセッサの動作を規定するソフトウェア(プログラム)との協働により実現されてもよい。
以下、図1を参照しながら、図2及び図3を用いて、フレーム列解析部32の処理の詳細を説明する。
図2は、胎児の心臓の4腔断面が示されたBモード画像が示されている。図2(a)には、拡張末期において取得されたフレームから形成される拡張末期画像50が示され、図2(b)には収縮末期において取得されたフレームから形成される収縮末期画像56が示されている。なお、以下の説明においては、簡単のためにBモード画像を用いて説明するが、実際にはフレーム列解析部32の処理はフレームに対して実施される。
図2(a)及び図2(b)を比較して分かるように、心臓の拡張・収縮に応じてBモード画像に示される心臓断面も拡張・収縮している。具体的には、図2(a)の拡張末期画像50においては、心臓断面52は心拍周期中において最大限に拡張し、心臓壁が伸び、各腔の面積が大きくなっている。一方、図2(b)の収縮末期画像56においては、心臓断面58は心拍周期中において最大限に収縮し、心臓壁が縮み、各腔の面積が小さくなっている。なお、本実施形態においては、心臓壁の部分(図2において斜線部分)が高輝度部分となっている。
上述の通り、フレーム列解析部32は、各フレームの輝度値の変化に基づいて心拍周期情報を取得する。1つのフレームの全領域の輝度値のフレームごとの変化に基づいて心拍周期情報を取得することも可能であるが、本実施形態では、より精度の高い心拍周期情報を取得するために、各フレームを複数の領域に分割し、各領域のうち輝度値変化の周期が最も安定した領域を特定し、当該領域から取得された輝度値のフレームごとの変化に基づいて心拍周期情報を取得する。
まず、フレーム列解析部32は、フレーム列に含まれる各フレームを16個(縦横4等分)の領域に分割する(図2参照)。そして、1つのフレームの16個の各領域について、領域に含まれる画素の輝度値の和(平均値でもよい)を算出する。つまり、1つのフレームについて16個の輝度値が算出される。そして、当該処理をフレーム列に含まれる全フレームについて行う。これにより、各領域に対応した16個の輝度値列(輝度値波形)が形成される。
フレーム列解析部32は、形成された各輝度値波形に対してノイズ除去のためにフィルタ(ローパスフィルタ)処理を行い、16個のフィルタ後波形を生成する。さらに、各フィルタ後波形について2次微分処理を行い、16個の2次微分波形を生成する。図3には、輝度値波形70、フィルタ後波形72、及び2次微分波形74の例が示されている。図3の例は、16個に分割された領域の1つの領域についての各波形が示されている。図3において横軸はフレーム番号(フレーム列を構成する各フレームに対して取得順に付した番号)を示す。つまり、横軸は時間軸と見ることもできる。縦軸は、輝度値波形70及びフィルタ後波形72については下端を0とする輝度値を表し、2次微分波形74については中央を0とする2次微分値を示す。
以下、1つの領域についての処理に注目して説明する。フレーム列解析部32は、フィルタ後波形に基づいて基準正弦波形を生成する。具体的には、まず、1つのフィルタ後波形から周期を抽出する。ここで、フィルタ後波形の周期が不安定な場合には多数の周期が抽出されることになる。そして、抽出された周期の平均値(AvgHR)を算出する。本実施形態では、フィルタ後波形から得られる周期のうち、1σ分の平均値をAvgHRとしている(すなわち極端な値を除外して算出している)。算出されたAvgHRを周期とする正弦波を基準正弦波形とする。
そして、フレーム列解析部32は、基準正弦波形形と当該領域に対応する2次微分波形との相関値を算出する。相関値の算出にあたっては、まず式(1)により相互相関波形を算出する。
ここで、g(x)は2次微分波形であり、FPSはフレームレート[Hz]を示す。また、xはフレーム番号を示し、yは位相調整量(何フレーム分位相をずらすか)を示す。Nはフレーム列に含まれるフレームの数である。
式(1)により算出される相互相関波形は、基準正弦波形と2次微分波形g(x)との位相が同位相の場合に最大値をとり、基準正弦波形と2次微分波形g(x)との位相が逆位相の場合に最小値をとる周期関数となる。そこで、フレーム列解析部32は、式(2)により相互相関波形の実効値RMSを算出し、算出された値を相関値とする。
算出された相関値が大きいということは、基準正弦波形と2次微分波形との類似度が高いということであり、これは、当該領域における輝度値波形の周期が安定していることを意味する。逆に言えば、輝度値波形の周期が不安定である領域については、基準正弦波形と2次微分波形との類似度が低く、相関値が低く算出される。
フレーム列解析部32は、上述の方法により、16個の領域全てについて相関値を算出する。そして、算出された複数の相関値の中から最も大きい値の相関値を特定する。そして、特定された最大相関値が算出された領域におけるフィルタ後波形あるいは2次微分波形を心拍周期情報(心拍周期波形)とする。本実施形態においては、図2に示す領域54について算出された相関値が最大であり、領域54におけるフィルタ後波形を心拍周期波形として利用する。図2を参照して分かるように、領域54においては、拡張末期において領域54の各画素の輝度値の和が最小となり、収縮末期において輝度値の和が最大となっている。したがって、心拍周期波形において極大値をとる時点が収縮末期を示し、極小値をとる時点が拡張末期を示している。
以下、図4を参照しながらフレーム特定部34の処理について説明する。図4には、フレーム列解析部32が取得した心拍周期波形80と、それに対応するフレーム列が示されている。フレーム特定部34は、心拍周期波形に基づいて、フレーム列の中から拡張末期において取得されたフレーム(拡張末期フレーム)及び収縮末期において取得されたフレーム(収縮末期フレーム)を特定する。具体的には、心拍周期波形80が極小値となる時点において取得されたフレームを拡張末期フレームとして特定し、心拍周期波形80が極大値となる時点において取得されたフレームを収縮末期フレームとして特定する。図4の例では、拡張末期フレーム82a及び82cが、収縮末期フレーム82b及び82dが特定される。
1組の拡張末期フレームと収縮末期フレームをフレームセットとして定義する。具体的には、1つの拡張末期フレームと、当該拡張末期フレームが取得された後最初に取得された収縮末期フレームとが1つのフレームセットとして定義する。図4の例では、拡張末期フレーム82a及び収縮末期フレーム82bからなるフレームセット84a、並びに拡張末期フレーム82c及び収縮末期フレーム82dからなるフレームセット84bが定義される。
以下、図5及び図6を参照しながら、フレームセットの表示態様及び計測について説明する。
図5に示される通り、表示処理部22は、取得されたフレーム列の中から選択された、拡張末期フレームと収縮末期フレーム(フレームセット)に対応する2つのBモード画像を並べて表示部24に表示する。これにより、ユーザはフレーム列を1枚1枚目視して所望の時相で得られたフレームセットを探し出す必要がなくなる。また、フレームセットに対応する2つのBモード画像が並べて表示されることでユーザは表示切り替え処理などを要することなく一目で計測に用いる画像として適切であるか否かを判断することができる。
また、表示処理部22は、フレーム特定部34が特定した複数のフレームセットのうち、計測開始時点あるいはフリーズ時点から最も近い心拍周期に対応するフレームセットを最初に表示させる。図4を参照すると、本実施形態では、時系列順に拡張末期フレーム82c、収縮末期フレーム82d、拡張末期フレーム82a、収縮末期フレーム82bの順で取得され、そこでフレームの取得が終了し(フリーズ操作され)計測が開始される。したがって、本実施形態では、フリーズ操作の直前の心拍周期に対応するフレームセット84aに対応する2つのBモード画像が最初に表示される。具体的には、拡張末期フレーム82aに対応するBモード画像86a及びと収縮末期フレーム82bに対応するBモード画像86bが最初に表示される(図5参照)。
ユーザは、フレームを取得しながら表示部に表示されるリアルタイム動画像を目視し、計測に用いるために適切な画像が撮れたと思った瞬間にフリーズ操作を行うのが一般的である。したがって、フレームセット84aが計測に用いるフレームセットとして適切である可能性が高い。そのため、フリーズ操作の直前の心拍周期に対応するフレームセットを最初に表示させることは、計測対象として適切である可能性の高いフレームセットを最初に表示させることを意味する。これにより、ユーザは、迅速に計測対象のフレームセットを選択することが可能になり、また、フレームセットの表示切り替えの手間の低減させることができる。
最初に表示されたフレームセットが計測に用いるのに適切なフレームセットではないとユーザが判断する場合も考えられる。したがって、表示処理部22は、ユーザの操作に応じて他の心拍周期に対応するフレームセットを切り替え表示させる。具体的には、ユーザが例えばダイアル式つまみを1目盛回転させた場合に、フレームセット84aとは異なる心拍周期において取得されたフレームセット84b(図4参照)に含まれる拡張末期フレーム82c及び収縮末期フレーム82dを並べて表示する。もちろん、フレームセットが3以上特定されている場合には、ユーザの操作に応じて次々とフレームセットが切り替わるようにしてもよい。これにより、最初に表示されたフレームセットが計測に用いるのに適切なフレームセットではない場合に、ユーザが他のフレームセットを検索するのにかかる手間及び時間を低減させることができる。
図6は、FS計測を実施した結果を示す図である。ユーザによりFS計測に用いるためのフレームセットが選択されると、当該フレームセットを用いてFS計測が行われる。FS値は、式(3)で算出される。
ここで、LVDdは拡張末期における左室内径であり、LVDsは収縮末期における左室内径である。したがって、FS計測のために、拡張末期フレーム82aに対応するBモード画像86a及び収縮末期フレーム82bに対応するBモード画像86b上のそれぞれにおいて左室内径が計測される。
ユーザは、Bモード画像86a上において、表示されるカーソルを動かし、左室内径計測の始点90a及び終点90bを指定する。すると、計測部36により始点90aと終点90b間の距離が計測される。左室内径は、Bモード画像86a上における距離とBモード画像86bの表示拡大率に基づいて算出される。同様にして、Bモード画像86b上において左室内径が算出され、2つの左室内径に基づいてFS値が算出される。
2つのBモード画像上における左室内径計測の始点及び終点は、ユーザにより手動で指定されてもよいが、本実施形態においては、トラッキング処理部38が、一方のBモード画像上において指定された始点及び終点のそれぞれについてフレーム列に含まれる各フレーム対してトラッキング処理を行うことで、他方のBモード画像における始点及び終点を自動的に指定する。
例えば、Bモード画像86aにおいて左室内径計測の始点90a及び終点90bがユーザにより指定された場合、トラッキング処理部38は、取得されたフレーム列において拡張末期フレーム82aから収縮末期フレーム82bに向けて(図4参照)始点90a及び終点90bのそれぞれについてトラッキングを行う。具体的には、拡張末期フレーム82aと、拡張末期フレーム82aの次に取得されたフレーム(以下「次フレーム」と記載)との間でパターンマッチング処理などを行う。当該処理により、拡張末期フレーム82の始点及び終点が指定された所定部位が次フレームにおいてどの位置へ移動したかを把握することができる。そして、次フレームにおける当該所定部位の位置に始点及び終点を指定する。この処理を収縮末期フレーム82bまで繰り返すことで、トラッキング処理部38は、収縮末期フレーム82bにおける始点90a及び終点90bの位置を指定する。
一方のBモード画像における始点及び終点がトラッキング処理により指定された場合には、指定された始点及び終点の位置が適切な位置であるか否かをユーザが確認する必要がある。そこで、表示処理部22は、2つのBモード画像のいずれか一方が表示された領域において、拡張末期フレームから収縮末期フレームまでの各フレームを動画像として表示する。当該動画像において、トラッキング処理により指定された各フレームにおける始点及び終点の位置が示される。したがって、当該動画像を再生すると、始点及び終点の位置が刻々と変化する様子が把握できる。ユーザは、当該動画像を確認することでトラッキングが正しく行われたか否かを容易に確認することができる。
2つのBモード画像において距離計測の始点及び終点が指定されると、計測部36によりFS値が算出される。図6に示されるように、算出されたFS値94は表示部24上に表示される。
また、取得されたフレーム列全体においてトラッキング処理を行い、フレーム列に含まれる各フレームにおいて始点及び終点を指定し、各フレームについて左室内径の計測を行ってもよい。そして、図6に示されるように、フレーム列における左室内径の推移を推移グラフ96として表示部24に表示してもよい。推移グラフ96の横軸はフレーム番号(時間)を示し、縦軸は左室内径を示している。画面上部に現在表示されているフレームセットの推移グラフ96における位置を強調表示するのも好適である。図6においては、推移グラフ96の一部が太線96aで表示されており、これが現在表示されているフレームセットを示している。推移グラフ96を表示することにより、ユーザはフレーム列における左室内径の推移を容易に把握することができる。
以下、図1を参照しつつ図7に示されるフローチャートを用いて、本実施形態の処理の流れを説明する。
ステップS10において、プローブ12から計測対象である胎児の心臓に対して一定期間(複数心拍に亘って)超音波の送受波が行われる。これによりフレーム列が取得される。一定期間超音波の送受波が行われるとユーザによりフリーズ操作が行われ、フレームの取得が停止される。取得されたフレーム列は画像メモリ20に記憶される。
ステップS12において、フレーム列解析部32は、取得されたフレーム列に含まれる各フレームの輝度値に基づいて胎児の心臓の心拍周期を示す心拍周期情報を取得する。
ステップS14において、フレーム特定部34は、ステップS12で心拍周期情報に基づいて、フレーム列の中から拡張末期において取得された拡張末期フレームと収縮末期において取得された収縮末期フレーム(フレームセット)を特定する。フレーム列は複数心拍に亘って取得されたフレームであるため、フレームセットは複数特定される。
ステップS16において、表示処理部22は、特定されたフレームセットに含まれる拡張末期フレーム及び収縮末期フレームに対応する2つのBモード画像を並べて表示部24に表示させる。ここで最初に表示されるフレームセットは、計測開始操作(フリーズ操作)から直近の心拍周期に対応するフレームセットである。
ステップS18において、制御部26は、ユーザからフレームセットの表示切り替え指示があったか否かを判定する。ユーザからフレームセットの表示切り替え指示があったということは、最初に表示されたフレームセットが計測対象として適切でないことを意味する。したがって、この場合、ステップS20において、表示処理部22は、他のフレームセットを表示させる。ステップS18及びS20の処理は、ユーザから表示切り替え指示がなくなるまで、すなわち計測に用いるフレームセットが確定されるまで繰り返し実行される。
ステップS22において、制御部26は、表示されたフレームセット、すなわち2つのBモード画像のうちいずれか一方において左室内径計測の始点及び終点が指定されたか否かを判定する。
Bモード画像のうちいずれか一方において距離計測の始点及び終点が指定された場合、ステップS24において、トラッキング処理部38は、一方のBモード画像(本例では拡張末期のBモード画像とする)において指定された始点及び終点についてトラッキングを行い、収縮末期のBモード画像上において左室内径計測の始点及び終点を指定する。
ステップS26において、表示処理部22は、拡張末期フレームから収縮末期フレームまでの各フレームを動画像としてループ再生する。このとき、トラッキング処理により拡張末期フレームから収縮末期フレームまでの各フレームについて指定された左室内径計測の始点及び終点が共に表示される。
ステップS28において、制御部26は、ユーザからFS値算出実行指示があったか否かを判定する。ユーザから当該指示がないということは、トラッキング処理により指定された始点及び終点の位置が適切でないことを意味するため、ステップS30において、収縮末期のBモード画像上における始点及び終点が修正される。当該修正はユーザにより行われてもよいし、トラッキング処理の条件などを変更して再度トラッキングを行うようにしてもよい。
ユーザからFS値算出実行指示がされると、ステップS32において、計測部36は、拡張末期及び収縮末期のBモード画像上においてそれぞれ左室内径を計測し、それに基づいてFS値を算出する。算出されたFS値は表示処理部22により表示部24に表示される。
以上説明した通り、本実施形態によれば、フレーム列の中から計測に用いるためのフレームセットが選択されそれらが表示される。このため、ユーザは格別の手間をかけることなく、フレーム列から計測に用いるフレームセットを選択することができる。また、最初に表示されたフレームセットが計測用として適切でなかった場合でも、計測用候補となるフレームセットが次々と切り替え表示されることから、ユーザは容易に計測用フレームセットを選択することができる。
なお、上述においては、胎児の心臓についてのFS計測が実施される例を用いて本発明の実施形態を説明したが、もちろん本発明の実施形態としてはこれに限られない。本発明は、周期的に運動する対象組織についての計測あって、対象組織の運動周期の中から特定時相の複数のフレームを用いる計測を実施する超音波画像処理装置に対して有効に適用することができる。
10 超音波診断装置、12 プローブ、14 送信部、16 受信部、18 画像形成部、20 画像メモリ、22 表示処理部、24 表示部、26 制御部、28 記憶部、30 入力部、32 フレーム列解析部、34 フレーム特定部、36 計測部、38 トラッキング処理部。
Claims (6)
- 周期的に運動する対象組織に超音波を送受波することで得られたフレーム列を解析することにより、前記対象組織の運動周期を示す周期情報を取得する周期情報取得手段と、
前記周期情報に基づいて、前記フレーム列の中から、前記運動周期中の第1時相に対応する第1時相フレーム及び前記運動周期中の第2時相に対応する第2時相フレームを含むフレームセットを特定するフレームセット特定手段と、
前記フレームセットを並べて表示部に表示させる表示制御手段と、
を備えることを特徴とする超音波画像処理装置。 - 前記表示制御手段は、前記フレーム列を用いた計測の開始操作時点から過去に直近の運動周期に対応するフレームセットを最初に表示させる、
ことを特徴とする、請求項1に記載の超音波画像処理装置。 - 前記表示制御手段は、ユーザの指示に従って、複数の運動周期に対応する複数のフレームセットを切り替え表示させる、
ことを特徴とする、請求項2に記載の超音波画像処理装置。 - 前記対象組織は胎児の心臓であり、
前記運動周期は心拍周期である、
ことを特徴とする、請求項1から3のいずれか1項に記載の超音波画像処理装置。 - 前記第1時相は、前記心拍周期における拡張末期であり、
前記第2時相は、前記心拍周期における収縮末期である、
ことを特徴とする、請求項4に記載の超音波画像処理装置。 - コンピュータを、
周期的に運動する対象組織に超音波を送受波することで得られたフレーム列を解析することにより、前記対象組織の運動周期を示す周期情報を取得する周期情報取得手段と、
前記周期情報に基づいて、前記フレーム列の中から、前記運動周期中の第1時相に対応する第1時相フレーム及び前記運動周期中の第2時相に対応する第2時相フレームを含むフレームセットを特定するフレームセット特定手段と、
前記フレームセットを並べて表示部に表示させる表示制御手段と、
として機能させることを特徴とする超音波画像処理プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014239090A JP2016097256A (ja) | 2014-11-26 | 2014-11-26 | 超音波画像処理装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014239090A JP2016097256A (ja) | 2014-11-26 | 2014-11-26 | 超音波画像処理装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2016097256A true JP2016097256A (ja) | 2016-05-30 |
Family
ID=56075692
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2014239090A Pending JP2016097256A (ja) | 2014-11-26 | 2014-11-26 | 超音波画像処理装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2016097256A (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019078054A1 (ja) * | 2017-10-17 | 2019-04-25 | 富士フイルム株式会社 | 音響波診断装置および音響波診断装置の制御方法 |
JP2019198389A (ja) * | 2018-05-14 | 2019-11-21 | キヤノンメディカルシステムズ株式会社 | 超音波診断装置、医用画像診断装置、医用画像処理装置及び医用画像処理プログラム |
CN110543912A (zh) * | 2019-09-02 | 2019-12-06 | 李肯立 | 自动获取胎儿关键切面超声视频中心动周期视频的方法 |
-
2014
- 2014-11-26 JP JP2014239090A patent/JP2016097256A/ja active Pending
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019078054A1 (ja) * | 2017-10-17 | 2019-04-25 | 富士フイルム株式会社 | 音響波診断装置および音響波診断装置の制御方法 |
CN111225619A (zh) * | 2017-10-17 | 2020-06-02 | 富士胶片株式会社 | 声波诊断装置及声波诊断装置的控制方法 |
EP3698722A4 (en) * | 2017-10-17 | 2020-12-23 | FUJIFILM Corporation | SOUND WAVE DIAGNOSTIC DEVICE AND METHOD OF CONTROLLING A SOUND WAVE DIAGNOSTIC DEVICE |
US11170519B2 (en) | 2017-10-17 | 2021-11-09 | Fujifilm Corporation | Acoustic wave diagnostic apparatus and control method of acoustic wave diagnostic apparatus |
CN111225619B (zh) * | 2017-10-17 | 2023-04-04 | 富士胶片株式会社 | 声波诊断装置及声波诊断装置的控制方法 |
US11636616B2 (en) | 2017-10-17 | 2023-04-25 | Fujifilm Corporation | Acoustic wave diagnostic apparatus and control method of acoustic wave diagnostic apparatus |
JP2019198389A (ja) * | 2018-05-14 | 2019-11-21 | キヤノンメディカルシステムズ株式会社 | 超音波診断装置、医用画像診断装置、医用画像処理装置及び医用画像処理プログラム |
CN110477952A (zh) * | 2018-05-14 | 2019-11-22 | 佳能医疗***株式会社 | 超声波诊断装置、医用图像诊断装置及存储介质 |
JP7136588B2 (ja) | 2018-05-14 | 2022-09-13 | キヤノンメディカルシステムズ株式会社 | 超音波診断装置、医用画像診断装置、医用画像処理装置及び医用画像処理プログラム |
CN110543912A (zh) * | 2019-09-02 | 2019-12-06 | 李肯立 | 自动获取胎儿关键切面超声视频中心动周期视频的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5414157B2 (ja) | 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム | |
US8526696B2 (en) | Medical imaging apparatus, medical image processing apparatus, ultrasonic imaging apparatus, ultrasonic image processing apparatus and method of processing medical images | |
JP4676334B2 (ja) | 生体信号モニタ装置 | |
JPWO2008023618A1 (ja) | 超音波診断装置 | |
JP2005342006A (ja) | 超音波診断装置、超音波画像処理装置、及び超音波信号処理プログラム | |
JP5292440B2 (ja) | 超音波診断装置 | |
JP2008080106A (ja) | 超音波診断装置および超音波診断装置のデータ処理プログラム | |
JP5820128B2 (ja) | 超音波診断装置 | |
JPWO2013105568A1 (ja) | 画像診断装置及び画像表示方法 | |
US9579084B2 (en) | Ultrasound diagnostic apparatus | |
JP5771306B1 (ja) | 超音波診断装置 | |
JP5842039B1 (ja) | 超音波画像処理装置 | |
JP2007143606A (ja) | 超音波診断装置 | |
JP2010269018A (ja) | 超音波診断装置 | |
JP2016097256A (ja) | 超音波画像処理装置 | |
EP2659839A1 (en) | Ultrasonic diagnosis device and image processing method | |
EP2716231B1 (en) | Ultrasonic diagnostic device | |
US20100234730A1 (en) | Image processing device, ultrasonic imaging apparatus including the same, and image processing method | |
US20220313214A1 (en) | Ultrasonic diagnostic apparatus, image processing apparatus, and image processing method | |
JP5784085B2 (ja) | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム | |
JP5346990B2 (ja) | 超音波診断装置 | |
JP2019000151A (ja) | 超音波診断装置及び関心領域設定方法 | |
JP4745455B2 (ja) | 超音波診断装置、超音波画像処理装置、及び超音波信号処理プログラム | |
JP2002330967A (ja) | 超音波診断装置 | |
JP2019000150A (ja) | 超音波診断装置及びサンプルゲート設定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20160518 |