JP2021000400A - 超音波診断装置、パノラマ画像生成方法及びプログラム - Google Patents

超音波診断装置、パノラマ画像生成方法及びプログラム Download PDF

Info

Publication number
JP2021000400A
JP2021000400A JP2019116842A JP2019116842A JP2021000400A JP 2021000400 A JP2021000400 A JP 2021000400A JP 2019116842 A JP2019116842 A JP 2019116842A JP 2019116842 A JP2019116842 A JP 2019116842A JP 2021000400 A JP2021000400 A JP 2021000400A
Authority
JP
Japan
Prior art keywords
image data
ultrasonic
unit
evaluation result
ultrasonic image
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2019116842A
Other languages
English (en)
Other versions
JP7298331B2 (ja
Inventor
亮汰 鳥井
Ryota Torii
亮汰 鳥井
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Konica Minolta Inc
Original Assignee
Konica Minolta Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Konica Minolta Inc filed Critical Konica Minolta Inc
Priority to JP2019116842A priority Critical patent/JP7298331B2/ja
Priority to US16/881,608 priority patent/US11721426B2/en
Publication of JP2021000400A publication Critical patent/JP2021000400A/ja
Application granted granted Critical
Publication of JP7298331B2 publication Critical patent/JP7298331B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13—Tomography
    • A61B8/14—Echo-tomography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/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/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A61B8/5246—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • A61B8/5276—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts due to motion
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/60—Editing figures and text; Combining figures or text
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
    • 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/48—Diagnostic techniques
    • A61B8/488—Diagnostic techniques involving Doppler signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image

Abstract

【課題】診断に有用なパノラマ画像データを得ることである。また、診断に有用なパノラマ画像データを効率よく得ることである。【解決手段】超音波診断装置１００は、移動されている超音波探触子２から得た受信信号に基づき超音波画像データを生成する画像生成部１４と、画像生成部１４で生成され画像メモリー部１５ａに記憶された複数の超音波画像データの合成の適否に関する指標を評価して評価結果を生成し、生成された評価結果に応じて、超音波画像データを選択して画像処理部１５に合成させてパノラマ画像データを生成する制御部１８と、を備える。【選択図】図２

Description

本発明は、超音波診断装置、パノラマ画像生成方法及びプログラムに関する。
超音波診断は、超音波探触子を患者の被検体の体表又は体腔内から当てるという簡単な操作で心臓や胎児の様子が超音波画像として得られ、かつ安全性が高いため繰り返して検査を行うことができる。このような超音波診断を行うために用いられる超音波診断装置が知られている。超音波画像データは、圧電素子を有する超音波探触子から超音波が被検体に送信され、反射した超音波を超音波探触子が受信し、その受信した信号に様々な処理を行うことで得られる。
また、超音波探触子を開口方向（方位方向（ラテラル方向）、走査方向）に移動させることによって得られた位置が隣り合う複数の超音波画像データを連続的に合成し、広い視野の超音波画像データ（パノラマ画像データ）を得るパノラマ画像取得機能を有する超音波診断装置が知られている。しかし、超音波探触子の縦方向の移動量が大きいと、観察対象の歪み量の違いにより、滑らかなパノラマ画像データが生成できないおそれがある。このため、複数の超音波画像の垂直方向のずれ幅を検出し、そのずれ幅を補正してパノラマ画像データを生成する超音波診断装置が知られている（特許文献１参照）。
また、複数の超音波画像の移動量に基づいて重畳画像（パノラマ画像データ）の信頼性を評価し、評価結果を表示する超音波診断装置が知られている（特許文献２参照）。
特開２００３−３８４８７号公報 特開２０１５−１４４６２３号公報
しかし、特許文献１の超音波診断装置では、垂直（縦）方向のずれ幅（移動量）を検出、補正するが、平行移動による補正のため合成画像中に歪みや不連続な箇所が発生し、診断に有用でないパノラマ画像データを得る可能性がある。また、特許文献２の超音波診断装置では、重畳画像の信頼性を移動量に基づいて評価し、ユーザーに提示するが、信頼性の低い箇所における合成画像は歪みなどを含み、診断に利用できない可能性があり、目的とする観察部位付近の信頼性が低い場合、超音波画像の再取得が必要となり、無駄な時間がかかってしまうという問題が生じる。
本発明の課題は、診断に有用なパノラマ画像データを得ることである。また、本発明の他の課題は、診断に有用なパノラマ画像データを効率よく得ることである。
上記課題を解決するため、請求項１に記載の発明の超音波診断装置は、
移動されている超音波探触子から得た受信信号に基づき超音波画像データを生成する画像生成部と、
前記画像生成部で生成された複数の超音波画像データの合成の適否に関する指標を評価して評価結果を生成する評価部と、
前記生成された評価結果に応じて、超音波画像データを選択して合成してパノラマ画像データを生成する合成部と、を備える。
請求項２に記載の発明は、請求項１に記載の超音波診断装置において、
前記合成部は、前記生成された評価結果が高い超音波画像データを選択する。
請求項３に記載の発明は、請求項１又は２に記載の超音波診断装置において、
前記複数の超音波画像データにおける縦方向の移動量及び前記超音波探触子の移動方向である横方向の移動量を算出する算出部を備え、
前記合成部は、前記算出された縦方向の移動量及び前記横方向の移動量に基づいて、前記選択された複数の超音波画像データを合成する。
請求項４に記載の発明は、請求項３に記載の超音波診断装置において、
前記縦方向は、前記超音波探触子の移動方向と略直交する方向である。
請求項５に記載の発明は、請求項３又は４に記載の超音波診断装置において、
前記算出部は、前記複数の超音波画像データにおける回転量を算出し、
前記合成部は、前記算出された回転量に基づいて、前記選択された複数の超音波画像データを合成する。
請求項６に記載の発明は、請求項１から５のいずれか一項に記載の超音波診断装置において、
前記指標は、前記複数の超音波画像データにおける被検体に対する前記超音波探触子の接触面に垂直な縦方向の移動量であり、
前記評価部は、前記縦方向の移動量が第１の閾値より大きい場合に前記評価結果を低くする。
請求項７に記載の発明は、請求項１から６のいずれか一項に記載の超音波診断装置において、
前記指標は、前記複数の超音波画像データにおける被検体に対する前記超音波探触子の接触面に平行な横方向の移動量に対する縦方向の移動量の移動量比であり、
前記評価部は、前記移動量比が第２の閾値より大きい場合に前記評価結果を低くする。
請求項８に記載の発明は、請求項１から７のいずれか一項に記載の超音波診断装置において、
前記指標は、前記複数の超音波画像データに設定される関心領域内の上側及び下側に配置された小関心領域の類似度の差分であり、
前記評価部は、前記小関心領域の類似度の差分が第３の閾値より大きい場合に前記評価結果を低くする。
請求項９に記載の発明は、請求項８に記載の超音波診断装置において、
前記複数の超音波画像データの各々における前記上側に配置された小関心領域は、略同一の超音波照射領域に対応する。
請求項１０に記載の発明は、請求項１から９のいずれか一項に記載の超音波診断装置において、
前記合成部は、前記複数の超音波画像データが記憶されている記憶部から、前記生成された評価結果が低い超音波画像データを削除する。
請求項１１に記載の発明は、請求項１から１０のいずれか一項に記載の超音波診断装置において、
前記合成部は、前記複数の超音波画像データの評価結果が連続して低い場合に、当該評価結果が連続して低い複数の超音波画像データの合成を中断する。
請求項１２に記載の発明は、請求項１１に記載の超音波診断装置において、
前記合成部は、前記合成が中断されたパノラマ画像データと、当該合成の中断を示す中断情報と、を表示部に表示させる。
請求項１３に記載の発明は、請求項１から１０のいずれか一項に記載の超音波診断装置において、
前記合成部は、前記複数の超音波画像データの評価結果が連続して低い場合に、当該評価結果が連続して低い複数の超音波画像データのうちの最も評価結果が高い超音波画像データを選択する。
請求項１４に記載の発明は、請求項１から１３のいずれか一項に記載の超音波診断装置において、
前記評価結果に基づき合成係数を決定する決定部を備え、
前記合成部は、前記決定された合成係数に基づいて、前記選択された複数の超音波画像データを合成する。
請求項１５に記載の発明は、請求項１４に記載の超音波診断装置において、
前記合成部は、前記生成されたパノラマ画像データと、当該パノラマ画像データの評価結果を示す評価結果情報と、を表示部に表示させる。
請求項１６に記載の発明のパノラマ画像生成方法は、
移動されている超音波探触子から得た受信信号に基づき超音波画像データを生成する画像生成工程と、
前記画像生成工程で生成された複数の超音波画像データの合成の適否に関する指標を評価して評価結果を生成する評価工程と、
前記生成された評価結果に応じて、超音波画像データを選択して合成してパノラマ画像データを生成する合成工程と、を含む。
請求項１５に記載の発明のプログラムは、
コンピューターを、
移動されている超音波探触子から得た受信信号に基づき超音波画像データを生成する画像生成部、
前記画像生成部で生成された複数の超音波画像データの合成の適否に関する指標を評価して評価結果を生成する評価部、
前記生成された評価結果に応じて、超音波画像データを選択して合成してパノラマ画像データを生成する合成部、
として機能させる。
本発明によれば、診断に有用なパノラマ画像データを得ることができる。また、本発明の別の効果としては、診断に有用なパノラマ画像データを効率よく得ることができる。
本発明の実施の形態の超音波診断装置の外観図である。 超音波診断装置の機能構成を示すブロック図である。 第１のパノラマ画像生成処理を示すフローチャートである。 （ａ）は、超音波画像と次フレームの超音波画像とにおけるＲＯＩを示す図である。（ｂ）は、超音波画像に対する次フレームの超音波画像の回転量を示す図である。 第２のパノラマ画像生成処理を示すフローチャートである。 第３のパノラマ画像生成処理を示すフローチャートである。 超音波画像と次フレームの超音波画像とにおける類似度が高い位置のＲＯＩを示す図である。 超音波画像と次フレームの超音波画像とにおけるＲＯＩと小ＲＯＩとを示す図である。 第４のパノラマ画像生成処理を示すフローチャートである。 （ａ）は、第１のパノラマ画像表示画面を示す図である。（ｂ）は、第２のパノラマ画像表示画面を示す図である。 第５のパノラマ画像生成処理を示すフローチャートである。 （ａ）は、第３のパノラマ画像表示画面を示す図である。（ｂ）は、第４のパノラマ画像表示画面を示す図である。
添付図面を参照して本発明に係る実施の形態、第１〜第４の変形例を順に詳細に説明する。なお、本発明は、図示例に限定されるものではない。
（実施の形態）
図１〜図４を参照して、本発明に係る実施の形態を説明する。まず、図１及び図２を参照して、本実施の形態の装置構成を説明する。図１は、実施の形態の超音波診断装置１００の外観図である。図２は、超音波診断装置１００の機能構成を示すブロック図である。
図１及び図２に示すように、本実施の形態の超音波診断装置１００は、超音波診断装置本体１と、超音波探触子２と、を備える。超音波探触子２は、図示しない患者の生体などの被検体に対して超音波（送信超音波）を送信するとともに、この被検体で反射した超音波の反射波（反射超音波：エコー）を受信する。超音波診断装置本体１は、ケーブル３を介して、超音波探触子２と接続され、超音波探触子２に電気信号の駆動信号を送信することによって超音波探触子２に被検体に対して送信超音波を送信させるとともに、超音波探触子２にて受信した被検体内からの反射超音波に応じて超音波探触子２で生成された電気信号である受信信号に基づいて被検体内の内部状態を超音波画像として画像化する。
超音波探触子２は、圧電素子からなる振動子を備えており、この振動子は、例えば、方位方向（走査方向）に一次元アレイ状に複数配列されている。本実施の形態では、例えば、１９２個の振動子を備えた超音波探触子２を用いている。なお、振動子は、二次元アレイ状に配列されたものであってもよい。また、振動子の個数は、任意に設定することができる。また、本実施の形態では、超音波探触子２について、リニア走査方式の電子スキャンプローブを採用したが、電子走査方式あるいは機械走査方式の何れを採用してもよく、また、リニア走査方式、セクタ走査方式あるいはコンベックス走査方式の何れの方式を採用することもできる。
図２に示すように、超音波診断装置本体１は、例えば、操作入力部１１と、送信部１２と、受信部１３と、画像生成部１４と、合成部としての画像処理部１５と、ＤＳＣ（Digital Scan Converter）１０６と、表示部１７と、評価部、合成部、算出部、決定部としての制御部１８と、記憶部１９と、を備える。
操作入力部１１は、例えば、診断開始を指示するコマンドや被検体の個人情報などのデータの入力などを行うための各種スイッチ、各種キー（ハードキー）、トラックボール、マウス、キーボードなどを備える操作卓（コンソール）であり、医師、技師などのユーザーからの操作入力に応じた操作信号を制御部１８に出力する。
送信部１２は、制御部１８の制御に従って、超音波探触子２にケーブル３を介して電気信号である駆動信号を供給して超音波探触子２に送信超音波を発生させる回路である。また、送信部１２は、例えば、クロック発生回路、遅延回路、パルス発生回路を備えている。クロック発生回路は、駆動信号の送信タイミングや送信周波数を決定するクロック信号を発生させる回路である。遅延回路は、駆動信号の送信タイミングを振動子毎に対応した個別経路毎に遅延時間を設定し、設定された遅延時間だけ駆動信号の送信を遅延させて送信超音波によって構成される送信ビームの集束を行うための回路である。パルス発生回路は、所定の周期で駆動信号としてのパルス信号を発生させるための回路である。上述のように構成された送信部１２は、例えば、超音波探触子２に配列された複数（例えば、１９２個）の振動子のうちの連続する一部（例えば、６４個）を駆動して送信超音波を発生させる。そして、送信部１２は、送信超音波を発生させる毎に駆動する振動子を方位方向にずらすことで走査（スキャン）を行う。
受信部１３は、制御部１８の制御に従って、超音波探触子２からケーブル３を介して電気信号である受信信号を受信する回路である。受信部１３は、例えば、増幅器、Ａ／Ｄ変換回路、整相加算回路を備えている。増幅器は、受信信号を、振動子毎に対応した個別経路毎に、予め設定された増幅率で増幅させるための回路である。Ａ／Ｄ変換回路は、増幅された受信信号をＡ／Ｄ変換するための回路である。整相加算回路は、Ａ／Ｄ変換された受信信号に対して、振動子毎に対応した個別経路毎に遅延時間を与えて時相を整え、これらを加算（整相加算）して音線データを生成するための回路である。
画像生成部１４は、制御部１８の制御に従って、受信部１３からの音線データに対して包絡線検波処理やログ圧縮などを実施し、ダイナミックレンジやゲインの調整を行って輝度変換することにより、Ｂモード画像データを生成する。すなわち、Ｂモード画像データは、受信信号の強さを輝度によって表したものである。画像生成部１４は、Ｂモード画像データの他、ドプラ法による画像データが生成できるものであってもよい。
画像処理部１５は、ＤＲＡＭ（Dynamic Random Access Memory）などの半導体メモリーによって構成された画像メモリー部１５ａ及びパノラマ画像メモリー部１５ｂを有する。画像メモリー部１５ａは、画像生成部１４から出力されたＢモード画像データをフレーム単位で記憶する。画像処理部１５は、制御部１８の制御に従って、画像生成部１４から出力されたＢモード画像データをフレーム単位で画像メモリー部１５ａに記憶し、フレームごとのＢモード画像データに適宜画像処理を施す。フレーム単位での画像データを超音波画像データ、あるいはフレーム画像データということがある。
画像処理部１５は、画像処理の一つとして、超音波探触子２が開口方向（方位方向、走査方向）に移動されたときに、位置が隣り合う複数のフレームのＢモード画像データを連続的に合成して広い視野のパノラマ画像データを生成する機能を有する。画像処理部１５は、例えば、平行移動、回転、拡大・縮小を考慮したヘルマート変換や、剪断変形まで考慮したアフィン変換を用いてパノラマ画像データへの合成を行う。ヘルマート変換は、
次式（１）、（２）などで表される。
Figure 2021000400
 ただし、ｘ’：変換後の水平座標、ｙ’：変換後の垂直座標、ｘ：変換前の水平座標、ｙ：変換前の垂直座標、Ｔ：変換行列、θ：回転量、Δｘ：水平平行移動量、Δｙ：垂直平行移動量、Ｒ：拡大・縮小率、である。
アフィン変換は、式（１）、次式（３）などで表される。
Figure 2021000400
 ただし、ｘ’：変換後の水平座標、ｙ’：変換後の垂直座標、ｘ：変換前の水平座標、ｙ：変換前の垂直座標、Ｔ：変換行列、ａ，ｂ，ｃ，ｄ：変換係数、である。
パノラマ画像メモリー部１５ｂは、画像処理部１５により生成されたパノラマ画像データを記憶する。画像メモリー部１５ａに記憶されたフレーム画像データ、パノラマ画像メモリー部１５ｂに記憶されたパノラマ画像データは、制御部１８の制御に従って、ＤＳＣ１６に出力される。
ＤＳＣ１６は、制御部１８の制御に従って、画像処理部１５より受信したフレーム画像データ、パノラマ画像データに座標変換などを施して表示部１７用の画像信号に変換し、表示部１７に出力する。
表示部１７は、ＬＣＤ（Liquid Crystal Display）、ＣＲＴ（Cathode-Ray Tube）ディスプレイ、有機ＥＬ（Electronic Luminescence）ディスプレイ、無機ＥＬティスプレイ及びプラズマディスプレイなどの表示装置が適用可能である。表示部１７は、ＤＳＣ１６から出力された画像信号に従って表示画面上に画像の表示を行う。なお、超音波診断装置１００は、表示部１７の画面上に設けられユーザーのタッチ入力を受け付けるタッチパネルを有する構成としてもよい。
制御部１８は、例えば、ＣＰＵ（Central Processing Unit）、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）を備えて構成され、ＲＯＭに記憶されているシステムプログラムなどの各種処理プログラムを読み出してＲＡＭに展開し、展開したプログラムに従って超音波診断装置１００の各部の動作を集中制御する。ＲＯＭは、半導体などの不揮発メモリーなどにより構成され、超音波診断装置１００に対応するシステムプログラム及び該システムプログラム上で実行可能な、例えば、後述する第１のパノラマ画像生成処理を実行するための第１のパノラマ画像生成プログラムなどの各種処理プログラムや、ガンマテーブルなどの各種データなどを記憶する。これらのプログラムは、コンピューターが読み取り可能なプログラムコードの形態で格納され、ＣＰＵは、当該プログラムコードに従った動作を逐次実行する。ＲＡＭは、ＣＰＵにより実行される各種プログラム及びこれらプログラムに係るデータを一時的に記憶するワークエリアを形成する。
記憶部１９は、例えば、ＨＤＤ（Hard Disk Drive）、ＳＳＤ（Solid State Drive）などの大容量記録媒体によって構成されており、患者情報に紐づいて保存されるフレーム画像データ、パノラマ画像データなどを記憶する。
超音波診断装置１００が備える各部について、各々の機能ブロックの一部又は全部の機能は、集積回路などのハードウェア回路として実現することができる。集積回路とは、例えばＬＳＩ（Large Scale Integration）であり、ＬＳＩは集積度の違いにより、ＩＣ（Integrated Circuit）、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩと呼称されることもある。また、集積回路化の手法はＬＳＩに限るものではなく、専用回路又は汎用プロセッサーで実現してもよいし、ＦＰＧＡ（Field Programmable Gate Array）やＬＳＩ内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサーを利用してもよい。また、各々の機能ブロックの一部又は全部の機能をソフトウェアにより実行するようにしてもよい。この場合、このソフトウェアは一つ又はそれ以上のＲＯＭなどの記憶媒体、光ディスク、又はハードディスクなどに記憶されており、このソフトウェアが演算処理器により実行される。
つぎに、図３、図４を参照して、本実施の形態の超音波診断装置１００の動作を説明する。図３は、第１のパノラマ画像生成処理を示すフローチャートである。図４（ａ）は、超音波画像３０Ａと次フレームの超音波画像３０ＢとにおけるＲＯＩ（Region Of Interest；関心領域）３１Ａ，３１Ｂを示す図である。図４（ｂ）は、超音波画像３０Ａに対する次フレームの超音波画像３０Ｂの回転量θを示す図である。
本実施の形態では、複数フレームの断層画像データとしてのＢモード画像データを用いてパノラマ画像データを生成する処理を説明するものとするが、これに限定するものではなく、カラードプラモードなど、他の断層画像データを用いてパノラマ画像データを生成する処理としてもよく、他の変形例でも同様である。
予め、超音波診断装置１００において、操作入力部１１を介する、医師、技師などの操作からの指示入力により、パノラマ画像表示モードに切り替えられているものとする。そして、ユーザーにより、患者の被検体に超音波探触子２が当てられて開口方向に移動され、制御部１８の制御で、送信部１２、超音波探触子２、受信部１３、画像生成部１４、画像処理部１５により、パノラマ画像データ生成用の連続する複数フレームのＢモード画像データ（超音波画像データ）が生成され、画像メモリー部１５ａに記憶され、これらの処理が、操作入力部１１を介してユーザーからパノラマ画像データ生成用の超音波画像データの取得終了の指示入力があるまで続けられるものとする。
そして、超音波診断装置１００において、パノラマ画像データ生成用の複数フレームの超音波画像データが画像メモリー部１５ａに記憶されたことをトリガーとして、制御部１８は、ＲＯＭに記憶された第１のパノラマ画像生成プログラムに従い、第１のパノラマ画像生成処理を実行する。
まず、制御部１８は、パノラマ画像データ生成用の時間的に最も古い１フレームの超音波画像データを画像メモリー部１５ａから読み出して取得する（ステップＳ１１）。そして、制御部１８は、ステップＳ１１で取得された１フレームの超音波画像データ又は削除されていなく直前に比較済（比較元）の超音波画像データの時間的に次の１フレームの超音波画像データを画像メモリー部１５ａから読み出して取得し、ステップＳ１１で取得された１フレームの超音波画像データ又は削除されていない比較済の１フレーム前の超音波画像データと、取得した超音波画像データとに、予め大きさが設定されたＲＯＩを設定する（ステップＳ１２）。例えば、図４（ａ）に示すように、時間的に前の超音波画像データの超音波画像３０Ａと、次フレームの超音波画像データの超音波画像３０Ｂとに、ＲＯＩ３１Ａ，３１Ｂが設定される。ＲＯＩ３１Ａ，３１Ｂの形状及び大きさは同じである。超音波画像３０Ａ中のＲＯＩ３１Ａは、位置を固定して設定される。超音波画像３０Ｂ中のＲＯＩ３１Ｂは、縦方向及び横方向に移動される。縦方向とは、パノラマ画像データ生成における超音波探触子２の移動方向と略直交する方向である。横方向とは、パノラマ画像データ生成における超音波探触子２の移動方向である。
そして、制御部１８は、ステップＳ１２で設定された２つの超音波画像データのＲＯＩ内の輝度値の類似度を算出する（ステップＳ１３）。そして、制御部１８は、ステップＳ１３で算出されたＲＯＩ内の輝度値の類似度に基づいて、比較元の超音波画像データに対するステップＳ１２で取得された次フレームの超音波画像データの縦方向の移動量、横方向の移動量及び回転量を算出する（ステップＳ１４）。例えば、図４に示す超音波画像３０Ａ中のＲＯＩ３１Ａに類似度が高い超音波画像３０Ｂ中のＲＯＩ３１Ｂの位置情報が取得され、ＲＯＩ３１Ａに対するＲＯＩ３１Ｂの相対位置情報から、被検体に対する超音波探触子２の接触面に垂直な縦方向の移動量及び横方向の移動量が算出される。また、図４（ｂ）に示すように、超音波画像３０ＡのＲＯＩ３１Ａを回転して、ＲＯＩ３１Ａに対する超音波画像３０ＢのＲＯＩ３１Ｂの類似度が高くなる回転量θが算出される。
そして、制御部１８は、ステップＳ１４で算出されたパノラマ画像データ生成用の合成の適否に関する指標としての縦方向の移動量が、予め設定された所定の縦方向の移動量用の閾値より大きいか否かを判別する（ステップＳ１５）。つまり、ステップＳ１５は、指標の評価に該当し、指標の値が大きいほど合成に適せず評価結果（信頼度）が低くなり、評価結果が閾値より低いか否かの判別に対応している。また、評価結果（信頼度）とは、超音波他像データがパノラマ画像データ生成用の合成に適する（信頼できる）ほど大きな値を示す。縦方向の移動量が閾値以下である場合（評価結果が閾値以上である場合）（ステップＳ１５；ＮＯ）、制御部１８は、比較した新しい方の１フレームの超音波画像データを画像メモリー部１５ａに保持する（ステップＳ１６）。
そして、制御部１８は、画像処理部１５により、ステップＳ１４で算出された縦方向の移動量、横方向の移動量及び回転量を用いて、パノラマ画像メモリー部１５ｂに記憶されているパノラマ画像データ（パノラマ画像データが記憶されていない場合はステップＳ１１で取得された１フレームの超音波画像データ）と、ステップＳ１２で取得された１フレームの超音波画像データと、を合成して新たなパノラマ画像データを生成させ、パノラマ画像メモリー部１５ｂに記憶（古いパノラマ画像データを更新）する（ステップＳ１７）。
そして、制御部１８は、未比較の次フレームの超音波画像データが画像メモリー部１５ａに記憶されているか否かを判別する（ステップＳ１８）。縦方向の移動量が閾値より大きい場合（評価結果が閾値より低い場合）（ステップＳ１５；ＹＥＳ）、制御部１８は、ステップＳ１２で取得された１フレームの超音波画像データを画像メモリー部１５ａから削除し（ステップＳ１９）、ステップＳ１８に移行される。
次フレームの超音波画像データが記憶されている場合（ステップＳ１８；ＹＥＳ）、ステップＳ１２に移行される。次フレームの超音波画像データが記憶されていない場合（ステップＳ１８；ＮＯ）、第１のパノラマ画像生成処理が終了する。第１のパノラマ画像生成処理の終了後、パノラマ画像メモリー部１５ｂに記憶されているパノラマ画像データは、制御部１８により、ＤＳＣ１６を介して表示部１７に表示され、操作入力部１１を介するユーザーからの保存指示入力に応じて記憶部１９に記憶される。
以上、本実施の形態によれば、超音波診断装置１００は、移動されている超音波探触子２から得た受信信号から超音波画像データを生成する画像生成部１４と、画像生成部１４で生成され画像メモリー部１５ａに記憶された複数の超音波画像データの合成の適否に関する指標を評価して評価結果を生成し、生成された評価結果に応じて、超音波画像データを選択して画像処理部１５に合成させてパノラマ画像データを生成する制御部１８と、を備える。制御部１８は、生成された評価結果が高い超音波画像データを選択する。
このため、指標の評価結果が高い超音波画像データを選択して合成に使用でき、指標の評価結果が低い超音波画像データを合成に使用しないことで、診断に有用なパノラマ画像データを得ることができるとともに、パノラマ画像データの再取得リスクを下げることができる。また、診断に有用なパノラマ画像データを効率よく得ることができる。
また、制御部１８は、複数の超音波画像データにおける縦方向の移動量及び超音波探触子２の移動方向である横方向の移動量を算出し、算出された縦方向の移動量及び横方向の移動量に基づいて、選択された複数の超音波画像データを合成する。縦方向は、超音波探触子２の移動方向と略直交する方向である。このため、複数の超音波画像データの合成の精度を上げることができる。
また、制御部１８は、複数の超音波画像データにおける回転量を算出し、算出された回転量に基づいて、選択された複数の超音波画像データを合成する。このため、回転量も合成に用いることで、湾曲した被検体をパノラマ撮像できる。
また、指標は、複数の超音波画像データにおける被検体に対する超音波探触子２の接触面に垂直な縦方向の移動量である。制御部１８は、縦方向の移動量が、縦方向の移動量用の閾値より大きい場合に評価結果を低くする。このため、パノラマ画像データ生成用の合成の適否を縦方向の移動量に応じて適切に評価できる。
また、制御部１８は、複数の超音波画像データが記憶されている画像メモリー部１５ａから、生成された評価結果が低い超音波画像データを削除する。このため、合成に不要な超音波画像データを破棄することで、画像メモリー部１５ａの容量を節約できる。
（第１の変形例）
図５を参照して、上記実施の形態の第１の変形例を説明する。図５は、第２のパノラマ画像生成処理を示すフローチャートである。
上記実施の形態では、指標として、縦方向の移動量を用いてパノラマ画像データ生成の合成用に適する超音波画像データを選択する構成であったが、本変形例では、指標として、横方向の移動量に対する縦方向の移動量である移動量比を用いてパノラマ画像データ生成の合成用に適する超音波画像データを選択する構成である。
本変形例では、上記実施の形態と同様に、装置構成として、超音波診断装置１００を用いる。ただし、制御部１８のＲＯＭには、第１のパノラマ画像生成プログラムに代えて、後述する第２のパノラマ画像生成処理を実行するための第２のパノラマ画像生成プログラムが記憶されているものとする。
つぎに、図５を参照して、本変形例の超音波診断装置１００の動作を説明する。上記実施の形態と同様に、超音波画像データの取得終了の指示入力があるまで、パノラマ画像データ生成用の被検体の複数フレームの超音波画像データが取得されて、画像メモリー部１５ａに記憶される。
超音波診断装置１００において、パノラマ画像データ生成用の複数フレームの超音波画像データが画像メモリー部１５ａに記憶されたことをトリガーとして、制御部１８は、ＲＯＭに記憶された第２のパノラマ画像生成プログラムに従い、第２のパノラマ画像生成処理を実行する。
図５に示すように、ステップＳ２１〜Ｓ２３は、図３の第１のパノラマ画像生成処理のステップＳ１１〜Ｓ１３と同様である。そして、制御部１８は、ステップＳ１３で算出されたＲＯＩ内の輝度値の類似度に基づいて、比較元の超音波画像データに対するステップＳ１２で取得された次フレームの超音波画像データの、縦方向の移動量と、横方向の移動量と、当該横方向の移動量に対する当該縦方向の移動量としての移動量比と、回転量とを算出する（ステップＳ２４）。
そして、制御部１８は、ステップＳ２４で算出された指標としての移動量比が、予め設定された所定の移動量比用の閾値より大きいか否かを判別する（ステップＳ２５）。つまり、ステップＳ２５は、指標の評価に該当し、指標の値が大きいほど評価結果（信頼度）が低くなり、評価結果が閾値より低いか否かの判別に対応している。移動量比が閾値以下である場合（評価結果が閾値以上である場合）（ステップＳ２５；ＮＯ）、ステップＳ２６が実行される。移動量比が閾値より大きい場合（評価結果が閾値より低い場合）（ステップＳ２５；ＹＥＳ）、ステップＳ２９が実行される。ステップＳ２６〜Ｓ２９は、図３の第１のパノラマ画像生成処理のステップＳ１６〜Ｓ１９と同様である。
以上、本変形例によれば、指標は、複数の超音波画像データにおける被検体に対する超音波探触子２の接触面に平行な横方向の移動量に対する縦方向の移動量の移動量比である。制御部１８は、移動量比が移動量比用の閾値より大きい場合に評価結果を低くする。このため、縦方向の移動量に横方向の移動量を加味して移動量比とすることで、過度に評価結果が低くならないようにできる。なお、移動量比だけでなく、縦方向、横方向それぞれの方向ベクトルのなす角度を指標にしてもよい。
（第２の変形例）
図６〜図８を参照して、上記実施の形態の第２の変形例を説明する。図６は、第３のパノラマ画像生成処理を示すフローチャートである。図７は、超音波画像３０Ａと次フレームの超音波画像３０Ｂとにおける類似度が高い位置のＲＯＩ３１Ａ，３１Ｂを示す図である。図８は、超音波画像３０Ａと次フレームの超音波画像３０ＢとにおけるＲＯＩ３１Ａ，３１Ｂと小ＲＯＩ３１１Ａ，３１２Ａ，３１１Ｂ，３１２Ｂとを示す図である。
上記実施の形態では、１つのＲＯＩを用いて得た縦方向の移動量を用いてパノラマ画像データ生成の合成用に適する超音波画像データを選択する構成であったが、本変形例では、ＲＯＩとそのＲＯＩ中の２つの小ＲＯＩとを用いてパノラマ画像データ生成の合成用に適する超音波画像データを選択する構成である。
本変形例では、上記実施の形態と同様に、装置構成として、超音波診断装置１００を用いる。ただし、制御部１８のＲＯＭには、第１のパノラマ画像生成プログラムに代えて、後述する第３のパノラマ画像生成処理を実行するための第３のパノラマ画像生成プログラムが記憶されているものとする。
つぎに、図６〜図８を参照して、本変形例の超音波診断装置１００の動作を説明する。上記実施の形態と同様に、超音波画像データの取得終了の指示入力があるまで、パノラマ画像データ生成用の被検体の複数フレームの超音波画像データが取得されて、画像メモリー部１５ａに記憶される。
超音波診断装置１００において、パノラマ画像データ生成用の複数フレームの超音波画像データが画像メモリー部１５ａに記憶されたことをトリガーとして、制御部１８は、ＲＯＭに記憶された第３のパノラマ画像生成プログラムに従い、第３のパノラマ画像生成処理を実行する。
図６に示すように、ステップＳ３１は、図３の第１のパノラマ画像生成処理のステップＳ１１と同様である。そして、制御部１８は、ステップＳ３１で取得された１フレームの超音波画像データ又は削除されていなく直前に比較済（比較元）の超音波画像データの時間的に次の１フレームの超音波画像データを画像メモリー部１５ａから取得し、ステップＳ３１で取得された１フレームの超音波画像データ又は削除されていない比較済の１フレーム前の超音波画像データと、取得した超音波画像データとに、予め大きさが設定されたＲＯＩと当該ＲＯＩ内の上側及び下側の２つの小ＲＯＩとを設定する（ステップＳ３２）。
例えば、図７に示すように、上記実施の形態のように、時間的に前の超音波画像データの超音波画像３０Ａと、次フレームの超音波画像データの超音波画像３０Ｂとに、ＲＯＩ３１Ａ，３１Ｂを設定するものとする。ＲＯＩ３１Ａ，３１Ｂ内には、被検体の上部（被検体の表面側、超音波探触子２側）の構造物ＳＴ１と下部（被検体の深部側）の構造物ＳＴ２とが含まれるものとする。また、超音波画像３０Ｂは、超音波探触子２の移動時に、被検体に圧迫を加えてしまった場合の超音波画像とする。また、ＲＯＩ３１Ｂは、ＲＯＩ３１Ａと輝度値の類似度が高い位置に配置されている。構造物ＳＴ１は、被検体の表面側にあるので、圧迫によって歪みやすい。構造物ＳＴ２は、構造物ＳＴ１よりも深い位置にあるので、圧迫によって歪みにくい（上方向への平行移動のみ起こっている）。このように、ＲＯＩが１つであると、当該ＲＯＩ内の上部又は下部の構造物が歪みを生じたままパノラマ画像データに合成されてしまう。
そこで、図８に示すように、超音波画像３０ＡのＲＯＩ３１Ａ内に小ＲＯＩ３１１Ａ（上部），３１２Ａ（下部）を設定する。小ＲＯＩ３１１Ａ，３１２Ａには、それぞれ、構造物ＳＴ１，ＳＴ２が含まれている。ＲＯＩ３１Ａ内の小ＲＯＩ３１１Ａ，３１２Ａの位置を固定し（例えば、ＲＯＩ３１Ａの上辺から下方にｚ１，ｚ２［ｍｍ］の位置に小ＲＯＩ３１１Ａ，３１２Ａ配置）、ＲＯＩ３１Ｂ内の小ＲＯＩ３１１Ｂ，３１２Ｂの位置も同様とする。小ＲＯＩ３１１Ｂ，３１２Ｂにも、それぞれ、構造物ＳＴ１，ＳＴ２が含まれている。超音波画像３０Ａ，３０Ｂの各々における上側に配置された小ＲＯＩ３１１Ａ，３１１Ｂは、略同一の超音波照射領域に対応する。
すると、ＲＯＩ３１Ａに対するＲＯＩ３１Ｂの移動により縦方向の移動量、横方向の移動量が算出できるとともに、小ＲＯＩ３１１Ａに対する小ＲＯＩ３１１Ｂの類似度と小ＲＯＩ３１２Ａに対する小ＲＯＩ３１２Ｂの類似度との差により、超音波画像３０Ｂの超音波画像データがパノラマ画像データ生成の合成用に適するか否かが判別できる。類似度との差が所定の閾値より大きいと、例えば小ＲＯＩ３１１Ａに対する小ＲＯＩ３１１Ｂの類似度が低くなり圧迫状態であるおそれがあり、超音波画像３０Ｂがパノラマ画像データ生成の合成用に適しない。ステップＳ３２では、このため、比較元のフレームの超音波画像データに固定の位置のＲＯＩ及びその内の上側及び下側の２つの小ＲＯＩが配置され、次フレームの超音波画像データに移動可能な位置のＲＯＩ及びその内の上側及び下側の２つの小ＲＯＩが配置される。
図６に戻り、制御部１８は、ステップＳ３２で設定された２つの超音波画像データのＲＯＩ、小ＲＯＩ内の輝度値の類似度を算出する（ステップＳ３３）。そして、制御部１８は、ステップＳ３３で算出されたＲＯＩ内の輝度値の類似度に基づいて、比較元の超音波画像データに対するステップＳ３２で取得された次フレームの超音波画像データの縦方向の移動量、横方向の移動量及び回転量を算出する（ステップＳ４４）。
ステップＳ３５，Ｓ３６，Ｓ３７，Ｓ３８，Ｓ４０は、図３の第１のパノラマ画像生成処理のステップＳ１５，Ｓ１６，Ｓ１７，Ｓ１８，Ｓ１９と同様である。縦方向の移動量が閾値より大きい場合（ステップＳ３５；ＹＥＳ）、制御部１８は、ステップＳ３３で算出された指標としてのＲＯＩ内の上側及び下側の２つの小ＲＯＩの類似度の差が予め設定された所定の類似度の差用の閾値より大きいか否かを判別する（ステップＳ３９）。例えば、図８の小ＲＯＩ３１１Ａに対する小ＲＯＩ３１１Ｂの類似度と、小ＲＯＩ３１２Ａに対する小ＲＯＩ３１２Ｂの類似度と、の差である。つまり、ステップＳ３９は、指標の評価に該当し、指標の値が大きいほど評価結果（信頼度）が低くなり、評価結果が閾値より低いか否かの判別に対応している。
上下の２つの小ＲＯＩの類似度の差が閾値以下の場合（評価結果が閾値以上である場合）（ステップＳ３９；ＮＯ）、ステップＳ３６に移行される。上下の２つの小ＲＯＩの類似度の差が閾値より大きい場合（評価結果が閾値より低い場合）（ステップＳ３９；ＹＥＳ）、ステップＳ４０に移行される。
以上、本変形例によれば、指標は、複数の超音波画像データに設定されるＲＯＩ内の上側及び下側に配置された小ＲＯＩの類似度の差分である。制御部１８は、小ＲＯＩの類似度の差分が類似度の差分用の閾値より大きい場合に評価結果を低くする。複数の超音波画像データの各々における上側に配置された小ＲＯＩは、略同一の超音波照射領域に対応する。このため、上側と下側との小ＲＯＩでそれぞれ類似度を算出することで、複数の超音波画像データの評価の精度を向上できる（過度に評価結果が低くならないようにすることができる）。
（第３の変形例）
図９、図１０を参照して、上記実施の形態の第３の変形例を説明する。図９は、第４のパノラマ画像生成処理を示すフローチャートである。図１０（ａ）は、パノラマ画像表示画面４０Ａを示す図である。図１０（ｂ）は、パノラマ画像表示画面４０Ｂを示す図である。
上記実施の形態では、指標としての縦方向の移動量が閾値以下の超音波画像データを用いて超音波画像データを生成する構成であったが、本変形例では、超音波画像データの指標の評価結果が連続して低い（信頼度が小さい）場合にパノラマ画像データ生成を自動的に中断する構成である。
本変形例では、上記実施の形態と同様に、装置構成として、超音波診断装置１００を用いる。ただし、制御部１８のＲＯＭには、第１のパノラマ画像生成プログラムに代えて、後述する第４のパノラマ画像生成処理を実行するための第４のパノラマ画像生成プログラムが記憶されているものとする。
つぎに、図９、図１０を参照して、本変形例の超音波診断装置１００の動作を説明する。上記実施の形態と同様に、超音波画像データの取得終了の指示入力があるまで、パノラマ画像データ生成用の被検体の複数フレームの超音波画像データが取得されて、画像メモリー部１５ａに記憶される。
超音波診断装置１００において、パノラマ画像データ生成用の複数フレームの超音波画像データが画像メモリー部１５ａに記憶されたことをトリガーとして、制御部１８は、ＲＯＭに記憶された第４のパノラマ画像生成プログラムに従い、第４のパノラマ画像生成処理を実行する。
図９に示すように、ステップＳ５１，Ｓ５２は、図３の第１のパノラマ画像生成処理のステップＳ１１〜Ｓ１４と同様とする。そして、制御部１８は、ステップＳ５２で算出された指標としての縦方向の移動量を評価し評価結果を算出する（ステップＳ５３）。例えば、図３の第１のパノラマ画像生成処理のステップＳ１５と同様に、縦方向の移動量が、移動量用の閾値より大きいか否かなどにより、値が高くなるほど信頼度（評価）が高いことを示す評価結果が生成される。
そして、制御部１８は、ステップＳ５３で生成された評価結果が予め設定された所定の評価結果用の閾値より小さいか否かを判別し、２回連続して評価結果が当該閾値より低いか否かを判別する（ステップＳ５４）。連続して評価結果が閾値より低くない場合（ステップＳ５４；ＮＯ）、制御部１８は、次フレームの超音波画像データが画像メモリー部１５ａに記憶されていなく、超音波画像データの取得を終了するか否かを判別する（ステップＳ５５）。超音波画像データの取得を終了しない場合（ステップＳ５５；ＮＯ）、ステップＳ５２に移行される。
超音波画像データの取得を終了する場合（ステップＳ５５；ＹＥＳ）、制御部１８は、画像処理部１５により、ステップＳ５２で算出された各フレームの縦方向の移動量及び横方向の移動量と回転量とを用いて、画像メモリー部１５ａに記憶されている全ての超音波画像データを合成して撮像した全フレームからなるパノラマ画像データを生成してパノラマ画像メモリー部１５ｂに記憶する（ステップＳ５６）。そして、制御部１８は、ステップＳ５６で生成されたパノラマ画像データを、ＤＳＣ１６を介して表示部１７に表示し（ステップＳ５７）、第４のパノラマ画像生成処理を終了する。
連続して評価結果が閾値より低い場合（ステップＳ５４；ＹＥＳ）、制御部１８は、画像処理部１５により、ステップＳ５２で算出された各フレームの縦方向の移動量及び横方向の移動量と回転量とを用いて、画像メモリー部１５ａに記憶されているステップＳ１１，Ｓ１２で取得済の超音波画像データを合成して途中までのフレームからなるパノラマ画像データを生成してパノラマ画像メモリー部１５ｂに記憶する（ステップＳ５８）。合成に用いる超音波画像データは、評価結果が連続して閾値より小さくなったフレームを含まないものとする。そして、制御部１８は、ステップＳ５８で生成されたパノラマ画像データと、パノラマ画像データ生成が途中までになった中断メッセージ欄とを、ＤＳＣ１６を介して表示部１７に表示し（ステップＳ５９）、第４のパノラマ画像生成処理を終了する。
ステップＳ５９では、例えば、図１０（ａ）に示すパノラマ画像表示画面４０Ａが表示される。パノラマ画像表示画面４０Ａは、パノラマ画像４１と、中断メッセージ欄４２Ａと、を有する。パノラマ画像４１は、例えばリニア走査方式の超音波探触子２が円形状に移動された場合の、評価結果が連続して低くなるまでのパノラマ画像データに基づくパノラマ画像である。中断メッセージ欄４２Ａは、パノラマ画像４１に重畳され、中断メッセージ及び閉じるボタンを有し、操作入力部１１を介するユーザーからの当該閉じるボタンのクリック入力により表示消去が可能である。また、ステップＳ５９では、図１０（ｂ）に示すパノラマ画像表示画面４０Ｂが表示されてもよい。パノラマ画像表示画面４０Ｂは、パノラマ画像４１と、中断メッセージ欄４２Ｂと、を有する。中断メッセージ欄４２Ｂは、パノラマ画像４１に重畳されず、中断メッセージを有する。中断メッセージ欄４２Ａ，４２Ｂのメッセージは、合成の中断の旨と、その理由とを含むものとする。
以上、本変形例によれば、制御部１８は、複数の超音波画像データの評価結果が連続して低い場合に、当該評価結果が連続して低い超音波画像データの合成を中断する。このため、評価結果が連続して低くなる前の診断に有用な超音波画像データを確実に合成してパノラマ画像データを生成できるとともに、診断に有用でない超音波画像データが合成されることを自動的に防ぐことができる。
また、制御部１８は、合成が中断されたパノラマ画像データと、合成の中断を示す中断メッセージ部と、を表示部１７に表示させる。このため、診断に有用でない超音波画像をユーザーに提示するのを防ぐことができ、ユーザーが中断メッセージ部により合成を中断した旨及びその理由を容易に視認できる。
（第４の変形例）
図１１、図１２を参照して、上記実施の形態の第４の変形例を説明する。図１１は、第５のパノラマ画像生成処理を示すフローチャートである。図１２（ａ）は、パノラマ画像表示画面５０Ａを示す図である。図１２（ｂ）は、パノラマ画像表示画面５０Ｂを示す図である。
上記実施の形態では、パノラマ画像データのみを生成し表示する構成であったが、本変形例では、超音波画像データの指標の評価結果情報をパノラマ画像データとともに表示する構成である。
本変形例では、上記実施の形態と同様に、装置構成として、超音波診断装置１００を用いる。ただし、制御部１８のＲＯＭには、第１のパノラマ画像生成プログラムに代えて、後述する第５のパノラマ画像生成処理を実行するための第５のパノラマ画像生成プログラムが記憶されているものとする。
つぎに、図１１を参照して、本変形例の超音波診断装置１００の動作を説明する。上記実施の形態と同様に、超音波画像データの取得終了の指示入力があるまで、パノラマ画像データ生成用の被検体の複数フレームの超音波画像データが取得されて、画像メモリー部１５ａに記憶される。
超音波診断装置１００において、パノラマ画像データ生成用の複数フレームの超音波画像データが画像メモリー部１５ａに記憶されたことをトリガーとして、制御部１８は、ＲＯＭに記憶された第５のパノラマ画像生成プログラムに従い、第５のパノラマ画像生成処理を実行する。
図１１に示すように、ステップＳ６１〜Ｓ６４は、図９の第３のパノラマ画像生成処理のステップＳ５１〜Ｓ５４と同様である。連続して評価結果が閾値より低くない場合（ステップＳ６４；ＮＯ）、ステップＳ６５に移行される。ステップＳ６５は、図６の第３のパノラマ画像生成処理のステップＳ５５と同様である。連続して評価結果が閾値より低い場合（ステップＳ６４；ＹＥＳ）、制御部１８は、評価結果が閾値より小さい連続した複数フレームのうち、ステップＳ６３で生成された評価結果が最も高いフレームの超音波画像データを画像メモリー部１５ａに保持し、当該複数フレームのうち、保持した以外のフレームの超音波画像データを画像メモリー部１５ａから削除し（ステップＳ６６）、ステップＳ６５に移行する。
超音波画像データの取得を終了する場合（ステップＳ６５；ＹＥＳ）、制御部１８は、ステップＳ６３で生成された評価結果に基づき合成係数を決定する（ステップＳ６７）。例えば、評価結果が高いほど、信頼度が高いため合成係数（合成の重み付け量）が高くなるように決定される。そして、制御部１８は、ステップＳ６２で算出された縦方向の移動量及び縦方向の移動量と回転量と、ステップＳ６７で決定された合成係数とを用いて、画像メモリー部１５ａに記憶されている超音波画像データを合成してパノラマ画像データを生成してパノラマ画像メモリー部１５ｂに記憶する（ステップＳ６８）。
ステップＳ６８でのパノラマ画像データの画素値の計算方法（重みづけ方法）は、次式（４）、（５）を用いる。
合成元の画像データにパノラマ画像データ（既に合成された画像データ）が存在する場合に、パノラマ画像データの画素値は、
Figure 2021000400
 で得られる。ただし、Ｏ：合成後の画素値、Ｉ：合成元の画像の画素値、Ｉ’：合成する画像の画素値、α：合成元の画像の信頼度（評価結果）、β：合成する画像の信頼度、である。
合成元の画像データにパノラマ画像データが存在しない場合に、パノラマ画像データの画素値は、
Ｏ（ｘ’、ｙ’）＝Ｉ’（ｘ’、ｙ’） …（５）
で得られる。
そして、制御部１８は、ステップＳ６８で生成されたパノラマ画像データと、パノラマ画像データ生成の評価結果（信頼度）を示す評価結果情報とを、ＤＳＣ１６を介して表示部１７に表示し（ステップＳ６９）、第５のパノラマ画像生成処理を終了する。
ステップＳ６９では、例えば、図１２（ａ）に示すパノラマ画像表示画面５０Ａが表示される。パノラマ画像表示画面５０Ａは、パノラマ画像５１Ａと、評価結果情報として評価結果バー５２Ａと、を有する。パノラマ画像５１Ａは、例えばリニア走査方式の超音波探触子２が円形状に移動された場合のパノラマ画像データに基づくパノラマ画像である。評価結果バー５２Ａは、パノラマ画像５１Ａの所定深さに沿って配置され、パノラマ画像５１Ａの横方向（超音波探触子２の移動方向）の各位置に対応する評価結果（信頼度）の値を白（信頼度が最も高い）⇔黒（信頼度が最も低い）の間で多段階にグレー表示するバーである。また、図１２（ｂ）に示すパノラマ画像表示画面５０Ｂが表示されることとしてもよい。パノラマ画像表示画面５０Ｂは、パノラマ画像５１Ｂと、評価結果バー５２Ｂと、を有する。パノラマ画像５１Ａは、例えばリニア走査方式の超音波探触子２が被検体に沿って略平行に移動された場合のパノラマ画像データに基づくパノラマ画像である。評価結果バー５２Ｂは、パノラマ画像５１Ｂの横方向（超音波探触子２の移動方向）の各位置に対応する評価結果（信頼度）の値をグレー表示するバーである。ただし、評価結果バー５２Ａ，５２Ｂは、白黒の濃淡表示としてのグレー表示に限定されるものではなく、他の色の濃淡表示や、カラーマッピングなど、他の表現方法で表現されてもよい。
以上、本変形例によれば、制御部１８は、複数の超音波画像データの評価結果が連続して低い場合に、当該評価結果が連続して低い複数の超音波画像データのうちの最も評価結果が高い超音波画像データを選択する。このため、できるだけ評価結果の高い超音波画像データを利用することで、診断に有用な超音波画像データを合成してパノラマ画像データを生成できる。
また、制御部１８は、評価結果に基づき合成係数を決定し、決定された合成係数に基づいて、選択された複数の超音波画像データを合成する。このため、自動的な合成中断を避けることでユーザーの操作を阻害することを防ぐことができ、合成係数を用いた合成により、診断に有用でない超音波画像データが大きく含まれるパノラマ画像データを生成するリスクを低減できる。
また、制御部１８は、生成されたパノラマ画像データと、パノラマ画像データの評価結果を示す評価結果バーと、を表示部１７に表示させる。このため、パノラマ画像中の診断に有用でない可能性のある箇所をユーザーが視認しやすくすることができる。
なお、上記実施の形態及び変形例における記述は、本発明に係る好適な超音波診断装置、パノラマ画像生成方法及びプログラムの一例であり、これに限定されるものではない。上記実施の形態、各変形例の少なくとも２つを適宜組み合わせる構成としてもよい。
また、以上の実施の形態及び変形例における超音波診断装置１００を構成する各部の細部構成及び細部動作に関して本発明の趣旨を逸脱することのない範囲で適宜変更可能である。
１００ 超音波診断装置
１ 超音波診断装置本体
１１ 操作入力部
１２ 送信部
１３ 受信部
１４ 画像生成部
１５ 画像処理部
１５ａ 画像メモリー部
１５ｂ パノラマ画像メモリー部
１６ ＤＳＣ
１７ 表示部
１８ 制御部
１９ 記憶部
２ 超音波探触子
２ａ 振動子
３ ケーブル

Claims (17)

  1. 移動されている超音波探触子から得た受信信号に基づき超音波画像データを生成する画像生成部と、
    前記画像生成部で生成された複数の超音波画像データの合成の適否に関する指標を評価して評価結果を生成する評価部と、
    前記生成された評価結果に応じて、超音波画像データを選択して合成してパノラマ画像データを生成する合成部と、を備える超音波診断装置。
  2. 前記合成部は、前記生成された評価結果が高い超音波画像データを選択する請求項１に記載の超音波診断装置。
  3. 前記複数の超音波画像データにおける縦方向の移動量及び前記超音波探触子の移動方向である横方向の移動量を算出する算出部を備え、
    前記合成部は、前記縦方向の移動量及び前記横方向の移動量に基づいて、前記選択された複数の超音波画像データを合成する請求項１又は２に記載の超音波診断装置。
  4. 前記縦方向は、前記超音波探触子の移動方向と略直交する方向である請求項３に記載の超音波診断装置。
  5. 前記算出部は、前記複数の超音波画像データにおける回転量を算出し、
    前記合成部は、前記算出された回転量に基づいて、前記選択された複数の超音波画像データを合成する請求項３又は４に記載の超音波診断装置。
  6. 前記指標は、前記複数の超音波画像データにおける被検体に対する前記超音波探触子の接触面に垂直な縦方向の移動量であり、
    前記評価部は、前記縦方向の移動量が第１の閾値より大きい場合に前記評価結果を低くする請求項１から５のいずれか一項に記載の超音波診断装置。
  7. 前記指標は、前記複数の超音波画像データにおける被検体に対する前記超音波探触子の接触面に平行な横方向の移動量に対する縦方向の移動量の移動量比であり、
    前記評価部は、前記移動量比が第２の閾値より大きい場合に前記評価結果を低くする請求項１から６のいずれか一項に記載の超音波診断装置。
  8. 前記指標は、前記複数の超音波画像データに設定される関心領域内の上側及び下側に配置された小関心領域の類似度の差分であり、
    前記評価部は、前記小関心領域の類似度の差分が第３の閾値より大きい場合に前記評価結果を低くする請求項１から７のいずれか一項に記載の超音波診断装置。
  9. 前記複数の超音波画像データの各々における前記上側に配置された小関心領域は、略同一の超音波照射領域に対応する請求項８に記載の超音波診断装置。
  10. 前記合成部は、前記複数の超音波画像データが記憶されている記憶部から、前記生成された評価結果が低い超音波画像データを削除する請求項１から９のいずれか一項に記載の超音波診断装置。
  11. 前記合成部は、前記複数の超音波画像データの評価結果が連続して低い場合に、当該評価結果が連続して低い複数の超音波画像データの合成を中断する請求項１から１０のいずれか一項に記載の超音波診断装置。
  12. 前記合成部は、前記合成が中断されたパノラマ画像データと、当該合成の中断を示す中断情報と、を表示部に表示させる請求項１１に記載の超音波診断装置。
  13. 前記合成部は、前記複数の超音波画像データの評価結果が連続して低い場合に、当該評価結果が連続して低い複数の超音波画像データのうちの最も評価結果が高い超音波画像データを選択する請求項１から１０のいずれか一項に記載の超音波診断装置。
  14. 前記評価結果に基づき合成係数を決定する決定部を備え、
    前記合成部は、前記決定された合成係数に基づいて、前記選択された複数の超音波画像データを合成する請求項１から１３のいずれか一項に記載の超音波診断装置。
  15. 前記合成部は、前記生成されたパノラマ画像データと、当該パノラマ画像データの評価結果を示す評価結果情報と、を表示部に表示させる請求項１４に記載の超音波診断装置。
  16. 移動されている超音波探触子から得た受信信号に基づき超音波画像データを生成する画像生成工程と、
    前記画像生成工程で生成された複数の超音波画像データの合成の適否に関する指標を評価して評価結果を生成する評価工程と、
    前記生成された評価結果に応じて、超音波画像データを選択して合成してパノラマ画像データを生成する合成工程と、を含むパノラマ画像生成方法。
  17. コンピューターを、
    移動されている超音波探触子から得た受信信号に基づき超音波画像データを生成する画像生成部、
    前記画像生成部で生成された複数の超音波画像データの合成の適否に関する指標を評価して評価結果を生成する評価部、
    前記生成された評価結果に応じて、超音波画像データを選択して合成してパノラマ画像データを生成する合成部、
    として機能させるためのプログラム。
JP2019116842A 2019-06-25 2019-06-25 超音波診断装置、パノラマ画像生成方法及びプログラム Active JP7298331B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2019116842A JP7298331B2 (ja) 2019-06-25 2019-06-25 超音波診断装置、パノラマ画像生成方法及びプログラム
US16/881,608 US11721426B2 (en) 2019-06-25 2020-05-22 Ultrasound diagnostic apparatus, panoramic image generation method, and recording medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019116842A JP7298331B2 (ja) 2019-06-25 2019-06-25 超音波診断装置、パノラマ画像生成方法及びプログラム

Publications (2)

Publication Number Publication Date
JP2021000400A true JP2021000400A (ja) 2021-01-07
JP7298331B2 JP7298331B2 (ja) 2023-06-27

Family

ID=73994530

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019116842A Active JP7298331B2 (ja) 2019-06-25 2019-06-25 超音波診断装置、パノラマ画像生成方法及びプログラム

Country Status (2)

Country Link
US (1) US11721426B2 (ja)
JP (1) JP7298331B2 (ja)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009029096A (ja) * 2007-07-30 2009-02-12 Fuji Kosho Kk 会員番号付きクレジットカードクレジットカード機能付き会員カード
JP2013075079A (ja) * 2011-09-30 2013-04-25 Fujifilm Corp 医用画像処理装置および方法ならびにプログラム
JP2015039521A (ja) * 2013-08-22 2015-03-02 セイコーエプソン株式会社 画像処理システム、超音波測定システム及び画像処理方法
JP2015144623A (ja) * 2012-05-14 2015-08-13 日立アロカメディカル株式会社 超音波診断装置および画像評価表示方法
JP2018068495A (ja) * 2016-10-26 2018-05-10 株式会社日立製作所 超音波画像処理装置及びプログラム

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5782766A (en) * 1995-03-31 1998-07-21 Siemens Medical Systems, Inc. Method and apparatus for generating and displaying panoramic ultrasound images
JP4785105B2 (ja) 2001-08-03 2011-10-05 株式会社日立メディコ 超音波画像処理装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009029096A (ja) * 2007-07-30 2009-02-12 Fuji Kosho Kk 会員番号付きクレジットカードクレジットカード機能付き会員カード
JP2013075079A (ja) * 2011-09-30 2013-04-25 Fujifilm Corp 医用画像処理装置および方法ならびにプログラム
JP2015144623A (ja) * 2012-05-14 2015-08-13 日立アロカメディカル株式会社 超音波診断装置および画像評価表示方法
JP2015039521A (ja) * 2013-08-22 2015-03-02 セイコーエプソン株式会社 画像処理システム、超音波測定システム及び画像処理方法
JP2018068495A (ja) * 2016-10-26 2018-05-10 株式会社日立製作所 超音波画像処理装置及びプログラム

Also Published As

Publication number Publication date
US20200411165A1 (en) 2020-12-31
US11721426B2 (en) 2023-08-08
JP7298331B2 (ja) 2023-06-27

Similar Documents

Publication Publication Date Title
CN103251429B (zh) 超声波成像装置
JP6881629B2 (ja) 超音波診断装置及び超音波診断装置の作動方法
JP2005511235A (ja) 超音波撮像システム及び方法
US10617394B2 (en) Ultrasound diagnostic apparatus, method of producing ultrasound image
JP6744141B2 (ja) 超音波診断装置及び画像処理装置
US10788459B2 (en) Ultrasound diagnostic apparatus, ultrasound image generation method, and recording medium
CN110575198B (zh) 解析装置及解析方法
JP7298331B2 (ja) 超音波診断装置、パノラマ画像生成方法及びプログラム
JP5869958B2 (ja) 超音波信号処理装置および超音波信号処理方法
JP5810631B2 (ja) 超音波診断装置
CN111839583A (zh) 超声波诊断装置以及分析装置
JP5337446B2 (ja) 超音波画像診断装置、画像処理装置及び超音波画像診断支援プログラム
JP2009039475A (ja) 超音波診断装置
JP4369206B2 (ja) 超音波画像生成装置
JP6834738B2 (ja) 超音波診断装置
JP2008212522A (ja) 超音波診断装置
US20220358643A1 (en) Medical image processing apparatus, ultrasonic diagnosis apparatus, and method
JP2012143358A (ja) 超音波診断装置及びプログラム
JP2002165795A (ja) 関心領域設定方法および装置並びに超音波撮影装置
JP6733339B2 (ja) 超音波画像診断装置
US20210128108A1 (en) Loosely coupled probe position and view in ultrasound imaging
JP6707363B2 (ja) 超音波診断装置
JP2020074996A (ja) 超音波診断装置及び超音波画像調整方法
JP3917929B2 (ja) 超音波診断装置
JP2005137422A (ja) 超音波画像生成装置および超音波画像生成方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20211224

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20221121

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20221213

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20230209

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230216

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230529

R150 Certificate of patent or registration of utility model

Ref document number: 7298331

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150