JP2018157968A - 放射線画像処理装置及び放射線画像撮影システム - Google Patents

放射線画像処理装置及び放射線画像撮影システム Download PDF

Info

Publication number
JP2018157968A
JP2018157968A JP2017056991A JP2017056991A JP2018157968A JP 2018157968 A JP2018157968 A JP 2018157968A JP 2017056991 A JP2017056991 A JP 2017056991A JP 2017056991 A JP2017056991 A JP 2017056991A JP 2018157968 A JP2018157968 A JP 2018157968A
Authority
JP
Japan
Prior art keywords
image
value
analysis
radiation
processing apparatus
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
JP2017056991A
Other languages
English (en)
Other versions
JP6805918B2 (ja
Inventor
志行 金子
Shiko Kaneko
志行 金子
敬子 板屋
Keiko Itaya
敬子 板屋
慎介 勝原
Shinsuke Katsuhara
慎介 勝原
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 JP2017056991A priority Critical patent/JP6805918B2/ja
Priority to US15/926,328 priority patent/US10878560B2/en
Publication of JP2018157968A publication Critical patent/JP2018157968A/ja
Application granted granted Critical
Publication of JP6805918B2 publication Critical patent/JP6805918B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46Arrangements for interfacing with the operator or the patient
    • A61B6/461Displaying means of special interest
    • A61B6/463Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46Arrangements for interfacing with the operator or the patient
    • A61B6/461Displaying means of special interest
    • A61B6/464Displaying means of special interest involving a plurality of displays
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/486Diagnostic techniques involving generating temporal series of image data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/503Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/504Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/507Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5223Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data generating planar views from image data, e.g. extracting a coronal view from a 3D image
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Optics & Photonics (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Dentistry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Human Computer Interaction (AREA)
  • Theoretical Computer Science (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Pulmonology (AREA)
  • Vascular Medicine (AREA)
  • Quality & Reliability (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

【課題】医療用のモダリティーで撮影された動画から、ポイントとなる情報だけを容易に読影できるようにする。【解決手段】放射線画像処理装置1は、複数のフレーム画像からなる放射線動画の動画データを取得する画像取得手段と、画像取得手段が取得した動画データに所定の解析処理を加え、複数の解析フレーム画像からなる解析結果動画を生成する解析手段と、解析手段が生成した解析結果動画を構成する複数の解析フレーム画像の中から複数の特定解析フレーム画像を選択する選択手段と、選択手段によって選択された複数の特定解析フレーム画像における、共通する座標上にそれぞれ位置する各共通座標画素の信号値に基づいて一の演算信号値を導出する演算処理手段と、演算処理手段が座標ごとに生成した複数の演算信号値を用いて前記表示部に演算結果画像を表示する画像表示手段と、を備える。【選択図】図2

Description

本発明は、放射線画像処理装置、及びこの放射線画像処理装置を備えた放射線画像撮影システムに関する。
医療用のモダリティー(例えば、CR(Computed Radiography)やCT(Computed Tomography)等)を用いて撮影した動画は、表示装置に表示されることで医者の診断に供される。その際、撮影された動画をそのまま表示するだけでなく、様々な解析処理を施してから表示することがある。
従来、動画を構成する複数のフレーム画像の一覧を表示部に表示し、解析に用いるフレーム画像を一覧の中から一括選択できるようにする技術(特許文献1参照)や、画像を即時表示、参照表示あるいは転送する際に、それぞれに最適な画像処理を行う技術(特許文献2参照)が知られている。
特開2012−110397号公報 特開2012−110466号公報
しかしながら、特許文献1,2に記載されたような従来の技術では、ユーザーが元の動画と解析を加えた動画、あるいは異なる種類の解析を施した動画同士を見比べながら読影することを考えたとき、各動画のそれぞれの関心のある領域を交互に視点を変えながら読影する必要が出てくる。各動画は並行して再生が進むので、一方の動画で詳しく読影したい領域を見つけ、他方の動画の対応する領域を読影するのに視点を変えた時には、本来見たいフレーム画像よりも後のフレーム画像が表示されることになってしまう。このため、一方の動画から他方の動画へ移る際に、フレーム画像を戻して読影しなければならない場合があった。
つまり、動画は情報量が多い上に時間の経過に伴って変化するため、読影(診断)が容易でなく、ユーザーに過度の負担がかかっていた。
本発明は、上記の問題点に鑑みてなされたものであり、医療用のモダリティーで撮影された動画から、ポイントとなる情報だけを容易に読影できるようにすることを課題とする。
前記の問題を解決するために、本発明は、
複数のフレーム画像からなる放射線動画の動画データを取得する画像取得手段と、
前記画像取得手段が取得した動画データに所定の解析処理を加え、複数の解析フレーム画像からなる解析結果動画を生成する解析手段と、
前記解析手段が生成した解析結果動画を構成する複数の解析フレーム画像の中から複数の特定解析フレーム画像を選択する選択手段と、
前記選択手段によって選択された複数の特定解析フレーム画像における、共通する座標上にそれぞれ位置する各共通座標画素の信号値に基づいて一の演算信号値を算出する演算処理手段と、
前記演算処理手段が座標ごとに生成した複数の演算信号値を用いて前記表示部に演算結果画像を表示する画像表示手段と、を備えることを特徴とする。
本発明によれば、医療用のモダリティーで撮影された動画から、ポイントとなる情報だけを容易に読影できるようにすることができる。
本発明の第1〜第3実施形態に係る放射線画像撮影システムの全体構成を示すブロック図である。 第1実施形態に係る放射線画像処理装置が行う画像処理の概念図である。 同実施形態に係る放射線画像処理装置が表示する画像処理結果の一例である。 同実施形態に係る放射線画像処理装置が行う画像処理の概念図である。 同実施形態に係る放射線画像処理装置が表示する画像処理結果の一例である。 第2実施形態に係る放射線画像処理装置が行う画像選択の概念図である。 同実施形態に係る放射線画像処理装置が行う画像選択の概念図である。 同実施形態に係る放射線画像処理装置が行う画像選択の概念図である。 同実施形態に係る放射線画像処理装置が行う画像選択の概念図である。 本発明の第3実施形態に係る放射線画像処理装置が行う画像選択の概念図である。 同実施形態に係る放射線画像処理装置が表示する画像処理結果の一例である。 同実施形態に係る放射線画像処理装置が表示する画像処理結果の一例である。 第1〜第3実施形態の変形例に係る放射線画像処理装置が行う画像処理の概念図である。 同変形例に係る放射線画像処理装置が表示する画像処理結果の一例である。 同変形例に係る放射線画像処理装置が行う画像処理の概念図である。 同変形例に係る放射線画像処理装置が行う画像処理の概念図である。 同変形例に係る放射線画像処理装置が行う画像処理の概念図である。 同変形例に係る放射線画像処理装置が表示する画像処理結果の一例である。
以下、図面を参照して本発明の実施の形態を詳細に説明する。ただし、発明の範囲は、図示例に限定されない。
<第1実施形態>
〔放射線画像撮影システムの構成〕
まず、本発明の第1実施形態に係る放射線画像撮影システム100の構成について説明する。図1は本実施形態に係る放射線画像撮影システム100の機能的構成を表すブロック図である。
本実施形態の放射線画像撮影システム100は、放射線画像処理装置1や、医療用のモダリティー2等で構成されている。
また、放射線画像撮影システム100には、必要に応じて、図示しないコンソールやサーバー(例えば、医療用画像管理システム(Picture Archiving and Communication System:PACS)を構成するもの)等が接続される。
放射線画像撮影システム100を構成する各装置は、DICOM(Digital Image and Communications in Medicine)規格に準じており、各装置間の通信は、DICOMに則って行われる。
モダリティー2は。放射線照射装置21と、放射線画像撮影装置22と、を備えている。
放射線照射装置21は、図示を省略するが、放射線を生成可能な回転陽極や回転陽極に電子ビームを照射するフィラメント等を有する放射線源や、設定された管電圧や管電流、照射時間(mAs値)等に応じた線量の放射線を放射線源から照射させるジェネレーター等を備えている。
そして、ユーザーの操作に基づいて、放射線画像撮影装置22へ放射線を照射するようになっている。
放射線画像撮影装置22は、図示を省略するが、放射線を受けることで線量に応じた電荷を蓄積する複数の放射線検出素子が二次元状(マトリクス状)に配列された基板や、各放射線検出素子に蓄積された電荷を画像データとして読み出す読み出し回路、外部装置と通信したり画像データを送信したりするための通信部等を備えている。
そして、放射線画像撮影装置22は、放射線照射装置2から放射線の照射を受け、画像データを読み出すと、その画像データを、通信部を介して直ちに外部へ送信するようになっている。
このような放射線照射装置2及び放射線画像撮影装置22を備えたモダリティー2は、放射線照射装置21による放射線の照射と、放射線画像撮影装置22による電荷の蓄積・画像データの読み出しと、を連続的に繰り返すことで、複数のフレーム画像からなる放射線動画の動画データを生成することが可能となっている。
なお、放射線照射装置21と放射線画像撮影装置22は、CTのように一体化されたものでもよい。
〔放射線画像撮影システムの構成〕
次に、上記放射線画像撮影システム100を構成する放射線画像処理装置1の構成について説明する。
放射線画像処理装置1は、コンピューター又は専用の制御装置として構成されており、図1に示したように、制御部11、通信部12、記憶部13、操作部14、表示部15、を備えて構成され、各部はバス16により接続されている。
制御部11は、CPU、RAM等により構成される。制御部11のCPUは、操作部14の操作に応じて、記憶部13に記憶されているシステムプログラムや、各種処理プログラムを読み出してRAM内に展開し、展開されたプログラムに従って、後述する画像解析処理を始めとする各種処理を実行したり、表示部15の表示内容を制御したりする等、放射線画像処理装置1各部の動作を集中制御する。
通信部12は、LANアダプタやモデムやTA等を備え、通信ネットワークNTに接続された放射線画像撮影装置22やコンソール、サーバー等をはじめとする各装置との間のデータ送受信を制御する。
記憶部13は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部13は、制御部11で各種処理(例えば、後述する画像解析処理や画像演算処理)を実行するためのプログラム(画像解析処理プログラムや画像演算処理プログラム等)や、プログラムにより処理の実行に必要なパラメータ、処理結果、あるいは、各種処理を施す前の元動画や静止画、各種処理を施した後の動画や静止画のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納されている。
操作部14は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードや、マウス等のポインティングデバイス等によりユーザーが操作可能に構成されており、キーボードに対するキー操作あるいはマウス操作により入力された指示信号を制御部11に出力する。
なお、操作部14は、表示部15の表示画面に備えられたタッチパネルで構成されていても良く、この場合、タッチパネルを介して入力された指示信号を制御部11に出力する。
表示部15は、LCDやCRT等のモニタにより構成され、制御部11から入力される表示信号の指示に従って、各種画像や、操作部14からの入力指示、データ等を表示することが可能となっている。
〔放射線画像処理装置の動作〕
次に、本実施形態に係る放射線画像処理装置1の動作について説明する。図2,4は放射線画像処理装置1における画像処理の概念図であり、図3,5は処理結果の表示例である。
なお、本実施形態に係る放射線画像処理装置1は、静止画像の処理や表示も可能であるが、ここでは、主に動画を処理・表示する場合について説明する。また、各図には、肺野領域の画像を例示したが、本実施形態に係る放射線画像処理装置1は、肺野以外の部位を撮影した動画についても処理・表示することが可能である。
放射線画像処理装置1の制御部11は、操作部14を介してユーザーから処理を施す前の元の放射線動画を表示する旨の指示を受けると、複数のフレーム画像からなる放射線動画の元動画データ(ローデータ等)の取得(通信部12を介してモダリティー2、コンソールもしくはサーバーから受信、又は記憶部13から読み込み)を行い、準備した元動画データに基づく元動画又は図2(a)に示したような元動画を構成する複数の元フレーム画像Iの一覧(画像群))を表示部15に表示することが可能となっている。すなわち、制御部11は、本発明における画像取得手段をなす。なお、図2(a)では、元フレーム画像を符号Iで示したが、以下の説明では、元動画であるか元フレーム画像であるかにかかわらず、元動画データに基づいて表示部15に表示されるものをIで示すこととする。
また、制御部11は、操作部14を介してユーザーから元動画データに動態解析を加える旨の指示を受けると、元動画データ(ローデータ等)の取得を行い、記憶部13に記憶されている画像解析処理プログラムとの協働により元動画データに各種解析処理(換気解析や血流解析等)を施し、その解析結果である解析動画データを生成することが可能となっている。
そして、制御部11は、生成した解析動画データに基づく解析結果動画又は図2(b)に示したような解析結果動画を構成する複数の解析フレーム画像Iの一覧)を、自動で又はユーザーからの指示に基づいて表示部15に表示することが可能となっている。すなわち、制御部11は、本発明における解析手段をなす。なお、図2(b)では、解析フレーム画像を符号Iで示したが、以下の説明では、解析結果動画であるか解析フレーム画像であるかにかかわらず、解析動画データに基づいて表示部15に表示されるものをIで示すこととする。
また、制御部11は、操作部14を介してユーザーから解析動画データに演算処理を加える旨の指示を受けると、解析動画データの取得を行い、記憶部13に記憶されている画像演算処理プログラムとの協働により解析動画データに各種演算処理を施し、その演算結果である演算画像データを生成することが可能となっている。
なお、元動画データに演算処理を施して演算画像データを生成することも可能である。
演算処理は、表示部15に表示された、解析結果動画を構成する複数の解析フレーム画像Iの一覧の中から複数の解析フレーム画像I(解析結果動画を構成する全ての解析フレーム画像Iであってもよいし、その一部であってもよい。)が選択されることにより開始される。以下、選択された解析結果フレーム画像Iを特定解析フレーム画像Iと称する。
本実施形態において、特定解析フレーム画像Iの選択は、ユーザーが操作部14に対して行った操作に基づいて行われるようになっている。具体的には、例えば図2(b)に示したように、表示部15に表示された解析結果動画を構成する解析フレーム画像Iの一覧の中から、ユーザーが処理対象とする特定解析フレーム画像Iを所定領域(例えば後述するサマリー表示領域)へドラッグ及びドロップしたり、解析フレーム画像Iの端部に設けられたチェックボックスBにチェックCを入れたりすることで選択される。すなわち、制御部11は、本発明における選択手段をなす。
なお、選択する解析フレーム画像Iは連続している必要は無く、間を飛ばして任意の解析フレーム画像Iを選択することも可能である。
演算処理では、選択された複数の特定解析フレーム画像Iにおける共通する座標上にそれぞれ位置する各画素の信号値に基づいて一の演算信号値を導出する。そして、この処理を各座標((1,1)から(m、n)まで)について繰り返し行い、座標ごとに生成した複数の演算信号値を用いて演算結果画像Iとなる演算画像データを生成する。すなわち、制御部11は、本発明における演算処理手段をなす。以下、各特定解析フレーム画像Iにおける共通する座標上にそれぞれ位置する各画素のことを共通座標画素と称する。
演算信号値の導出の仕方としては、例えば、最大値(MIP)処理、最小値(MinP)処理、Pk−Pk値処理、標準偏差処理、最頻値処理、積算値処理、平均値処理、最大微分値処理、最小微分値処理等がある。制御部11は、自動で又はユーザーからの指示に基づいて、上記の処理の中から何れかを設定して演算信号値を導出する。
最大値処理は、複数の特定解析フレーム画像Iにおける、各共通座標画素の信号値の中から最大値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値(換気や血流)の最大値が分かるので、換気や血流の有無を一目で把握することができる。具体的には、無気肺や肺塞栓領域において、換気機能や血流機能がどこまで回復しているのかを把握しやすくなる。また、処理コストを低く抑えることができる。
最小値処理は、複数の特定解析フレーム画像Iにおける、各共通座標画素の信号値の中から最小値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値(換気や血流)の最小値が分かるので、換気機能や血流機能がどこまで低下しているのかを把握しやすくなる。また、処理コストを低く抑えることができる。
最頻値処理は、複数の特定解析フレーム画像Iにおける、各共通座標画素の信号値の中から最頻値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値(換気や血流)の最頻値が分かるので、特異的なノイズまたはアーチファクトの影響を低減させた上で、換気や血流の有無を一目で把握できるようになる。例えば複数呼吸分、複数心拍分の解析結果に対して最頻値処理を実施することで、呼吸や心拍のばらつき影響を低減させることができる。
Pk−Pk値処理は、複数の特定解析フレーム画像Iにおける、各共通座標画素の信号値の中から最大値と最小値を抽出してその差分または比を算出し、特定解析フレーム画像の数に対する最大値と最小値との差の割合を算出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値(換気や血流)の変動量が分かるので、換気や血流の有無を一目で把握できるようになる。例えば、無気肺や肺塞栓領域の把握を行いやすくなる(信号の変動に強い、ロバストな結果を出すことができる)。
標準偏差処理は、複数の特定解析フレーム画像Iにおける、各共通座標画素の信号値の標準偏差を算出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値(換気や血流)の変動量が(複雑な計算アルゴリズム不要で)分かるので、換気や血流の有無を一目で把握できるようになる。例えば、無気肺や肺塞栓領域の把握を行いやすくなる。また、処理コストを低く抑えることができる。
積算値処理は、複数の特定解析フレーム画像Iにおける、各共通座標画素の信号値の積算値を算出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値(換気や血流)のある一定時間分の積算結果が分かるので、ノイズまたはアーチファクトの影響を低減させることができ、換気や血流の有無の把握が行いやすくなる。また、複数呼吸分、複数心拍分の解析結果を積算することで、呼吸や心拍のばらつき影響を低減させることができる。また、ある一定時間の間に計測した値を積算した値、または画像を最終結果として出力する他モダリティー(例えばシンチグラフィ)と、類似した物理量を導出できる為、積算値処理を行うことで上記のような他モダリティーの情報と比較が行いやすくなる。その結果、診断の精度が向上する。更に、画素値は線量と関係が有る場合がある為、総照射線量(総被曝量)の推定が行いやすくなる為、被検者の被ばく管理などに有用である。
平均値処理は、複数の特定解析フレーム画像Iにおける、各共通座標画素の信号値の平均値を算出して演算信号値とする処理である。この処理により得られた解析画像からは、ノイズまたはアーチファクトの影響を低減させた上で、換気や血流の有無を一目で把握できるようになる。例えば複数呼吸分、複数心拍分の解析結果に対して平均値処理を実施することで、呼吸や心拍のばらつき影響を低減させることができる。
最大微分値処理は、複数の特定解析フレーム画像Iの前後する2枚のフレーム画像における各共通座標画素の信号値の差分(または比)を算出し、得られた各差分(または各比)の中の最大値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、手軽に時間変化速度の最大値あるいは、時間変化加速度の最大値を把握することができる。例えば呼気時のフレームを特定解析フレームとするとこで、息の吐きやすさを評価し易くなるため、閉塞性疾患や拘束性疾患など、動きや加速度に変化が生じる疾患領域を検出しやすくなる。
最小微分値処理は、複数の特定解析フレーム画像Iの前後する2枚のフレーム画像における各共通座標画素の信号値の差分または比を算出し、得られた各差分または各比の中の最小値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、手軽に時間変化速度の最小値あるいは、時間変化加速度の最小値を把握することができる。例えば呼気時のフレームを特定解析フレームとするとこで、息の吐きやすさを評価し易くなるため、閉塞性疾患や拘束性疾患など、動きや加速度に変化が生じる疾患領域を検出しやすくなる。
なお、制御部11は、一度に異なる複数種類の解析結果動画に演算処理を加える旨の指示を受ける(サマリー表示領域15Aにドラッグ及びドロップされる)と、解析結果動画ごとに演算処理を行うことが可能となっている。
また、制御部11は、上記演算処理を施すことで生成した演算画像データに基づく演算結果画像Iを、自動で又はユーザーからの指示に基づいて表示部15に表示することが可能となっている。すなわち、制御部11は、本発明における画像表示手段をなす。
また、本実施形態の放射線画像処理装置1は、表示部15の表示領域の少なくとも一部にサマリー表示領域15Aを形成することが可能となっている。サマリー表示領域15Aは、主に、上述した演算結果画像Iを表示するための領域である。
サマリー表示領域15Aは、常時表示されるものであってもよいし、必要に応じて表示・非表示が切り替えられるものであってもよい。また、サマリー表示領域15Aは、表示部15の一部に表示されるものであっても、全体に表示されるものであってもよい。
なお、表示部15とは異なる他の表示装置を備え、そこによりサマリー表示領域15Aを表示するようにすることも可能である。
サマリー表示領域15Aにおいては、演算結果画像I等を様々な表示態様で表示することが可能となっている。本実施形態においては、単独表示態様、並列表示態様、重畳表示態様、切替表示態様の4種類の表示態様がある。
また、演算結果画像Iを、サマリー表示領域15Aだけではなく、表示部15に表示するようにしてもよい。
単独表示態様は、例えば図2(c),(d)又は(e)に示したように、一のサマリー表示領域15Aに一の演算結果画像Iを表示する態様である。
なお、単独表示態様のサマリー表示領域15Aを、同時に複数表示するようにしてもよい。
重畳表示態様は、図2(f)や図4(c)に示したように、演算結果画像Iを、元動画Iや、元フレーム画像I0、解析結果動画I、解析フレーム画像I、他の演算結果画像Iに重畳させて表示する表示態様である。因みに、図2(f)は、演算結果画像Iを元動画又は元フレーム画像Iに重畳した例、図4(c)は異なる解析処理を施した二種類の解析画像データにそれぞれ同じ演算処理(MIP処理)を施して得られた二枚の演算結果画像I同士を重畳した例である。
なお、3種類以上の画像を重畳表示することも可能である。
以下、重畳表示態様で表示された画像を重畳画像Iと称することがある。
並列表示態様は、図3や図4(b),(d)に示したように、演算結果画像Iを、元動画Iや、元フレーム画像I0、解析結果動画I、解析フレーム画像I、他の演算結果画像I、重畳画像Iと並べて表示する表示態様である。因みに、図3は特定解析フレーム画像Iの選択枚数だけを変えて同一の演算処理(MIP処理)を施した複数の演算結果画像Iを並べた例、図4(b)は異なる解析処理を施した二種類の解析動画データにそれぞれ同じ演算処理(MIP処理)を施して得られた二枚の演算結果画像Iを並べた例、図4(d)は、異なる複数の演算結果画像I、これらの重畳画像I及び元動画又は元フレーム画像Iを並べた例である。
また、上述したように、異なる複数種類の演算処理を行った場合には、各演算結果を並列表示することが可能である。
切替表示態様は、図5に示したように、演算結果画像I2、元動画I、元フレーム画像I、解析結果動画I、解析フレーム画像I、他の演算結果画像I及び重畳画像Iの中から一の画像を所定時間表示した後、当該画像と入れ替えで当該画像とは異なる画像を表示する表示態様である。因みに、図5はある演算結果画像I、他の演算結果画像I、これらの重畳画像I、元動画又は元フレーム画像Iの順で変更した場合の例である。
なお、並列表示態様と切り替え表示態様を組み合わせることも可能である。すなわち、一のサマリー表示領域15Aに複数の画像を切り替えながら表示することも可能である。
また、上記表示態様の説明では、サマリー表示領域15Aに表示される画像には、演算結果画像Iが必ず含まれるものとしたが、含まない形で表示することも可能である。すなわち、元動画又は元フレーム画像Iと解析結果動画又は解析フレーム画像Iを重畳表示、並列表示、切替表示することも可能である。
上述したように、本実施形態に係る放射線画像処理装置1は、複数の元フレーム画像からなる元動画(放射線動画)の動画データを取得する画像取得手段と、画像取得手段が取得した動画データに所定の解析処理を加え、複数の解析フレーム画像からなる解析結果動画を生成する解析手段と、解析手段が生成した解析結果動画を構成する複数の解析フレーム画像の中から複数の特定解析フレーム画像を選択する選択手段と、選択手段によって選択された複数の特定解析フレーム画像における、共通する座標上にそれぞれ位置する各共通座標画素の信号値に基づいて一の演算信号値を導出する演算処理手段と、前記演算処理手段が座標ごとに生成した複数の演算信号値を用いて前記表示部に演算結果画像を表示する画像表示手段と、を備えたものとなっている。
サマリー表示領域15Aにおいて、上述したような動画又は静止画像の表示を行うことにより、ユーザーにとって関心の高い演算結果画像Iを直感的に得ることができる。その結果、ユーザー自身が診断に適切な表示方法を容易に探ることができ、読影精度が向上する。
また、表示部15に表示される動画のどこを注力して読影するべきかユーザーが理解し易くなり、効率的に読影できるようになるとともに、読影洩れを防止することができる。その結果、診断の精度を向上させることができる。
また、任意のフレーム画像を選択可能としているので、アーチファクト(ノイズ)が混入していると推定されるフレーム画像を容易に取り除くことができ、解析結果の精度が向上する。
すなわち、動画から、ポイントとなる情報だけを容易に読影することができ、ユーザーの負担が低減する。
また、図3に示したように、フレーム画像の選択枚数だけを変えて同一の演算処理(MIP処理)を施した複数の演算結果画像Iを並列表示すれば、演算処理の妥当性、すなわち、選択フレームの範囲を広くしすぎたのか、もう少し広げられそうか等をユーザーに直感的にすばやく把握させることができる。
また、図4に示したように、異なる種類の複数の動画(フレーム画像の一覧)からそれぞれ選択したフレーム画像に対して演算処理した結果を並列表示したり、複数の動画からそれぞれ選択したフレーム画像に対して演算処理した結果を重畳表示したりすることで、診断に有用な「V/Q(換気血流比)画像」をユーザーの手で容易に作成することが可能となる。
<第2実施形態>
次に、本発明の第2実施形態に係る放射線画像処理装置について説明する。図6〜9は本実施形態に係る放射線画像処理装置が行う画像選択の概念図である。なお、図6には、肺野領域の画像を例示したが、本実施形態に係る放射線画像処理装置1は、肺野以外にも、心臓領域、肺動脈、大動脈又は肺気管支等、所定周期で同じ動作を繰り返す部位の動画であれば扱うことが可能である。
第1実施形態では、演算処理を施す対象となる特定解析フレーム画像Iをユーザーが選択するようにしていたが、本実施形態に係る放射線画像処理装置1Aは、予め設定された選択パターンに基づいて特定解析フレーム画像Iを自動的に選択するようになっている。
このため、本実施形態に係る放射線画像処理装置1Aは、実行する制御(記憶部13Aの記憶内容)が第1実施形態と異なっている。なお、記憶部13Aの記憶内容を除く他の構成は第1実施形態と同様である。
すなわち、放射線画像処理装置1Aは、記憶部13Aに、自動選択の範囲を規定する選択パターンを複数格納しており、ユーザーが操作部14に対して行った操作に基づいて、予め用意された複数の選択パターンの中から何れかの選択パターンを設定する。そして、ユーザーが特に注意して読影したい関心領域(例えば図6(a)で符号Rで示したような領域)を選択すると、設定した選択パターンに基づいて前記特定解析フレーム画像を自動的に選択するようになっている。
本実施形態に係る放射線画像処理装置1Aにおいて設定可能な選択パターンには、大きく分けて(1)極大値(極小値)−極大値(極小値)、(2)極大値(極小値)−極小値(極大値)、(3)極大値(極小値)とその前後所定範囲、の3つがある。
なお、このうちのいずれかのパターンのみが予め設定されている(ユーザーによる設定はできない)ようにしてもよいし、このうちの2パターンの中からユーザーが設定できるものとしてもよい。更に、ここで例示した以外のパターンを含むようにしてもよい。
選択パターン(1)が選択されると、放射線画像処理装置1は、少なくとも肺野、心臓領域、肺動脈、大動脈、あるいは肺気管支の一部が含まれる関心領域R(Region Of Interest;ROI)内の複数の画素信号値の平均値の時間方向の変化を数値化する(図6(b)に示したようにグラフ化する)。そして、図6(b),(c)に示したように、ある極大値(極小値)から、当該極大値(極小値)と隣り合う極大値(極小値)までの間のフレーム画像(一周期分の動作を表す動画;選択範囲r)、又はある極大値(極小値)から、当該極大値と隣り合わない極大値(極小値)までの間の解析フレーム画像(複数周期分の動作を表す動画;選択範囲r)を特定解析フレーム画像として自動選択する。
なお、rとrのどちらの選択範囲とするかは、ユーザーが選択パターンを設定する際に併せて設定するよう構成してもよいし、制御部11が自動的に判断するよう構成してもよい。
診断対象が図6(a)に示したような肺野である場合、この選択範囲r,rは、吸気位から再度吸気位となるまでの動きを表すことになる。一方、診断対象が心臓である場合には、心収縮期から再度心収縮期となるまでの動きを表すことになる。
肺機能を把握する上では、肺の一呼吸分の動作や心臓の一心拍分の動作等、少なくとも一周期分の情報を読影することが重要である。このため、ここで挙げたような解析フレーム画像の選択パターンを設定することで、呼吸周期や心拍周期に関係する肺機能(例えば、肺への空気の出入り量(1回換気量)、肺野内への血液の出入り量(心拍出量、1回拍出量))を正確に把握することができる。
特に、複数周期に亘る選択範囲rを設定することで、呼吸や心拍のバラつきを低減することが可能である。
選択パターン(2)が選択されると、放射線画像処理装置1は、関心領域R内の平均値の時間方向の変化を数値化(グラフ化)する。そして、図7(a)に示したように、ある極大値から、当該極大値の後(グラフでは右側)に隣り合う極小値までの間の解析フレーム画像(選択範囲r)、又は図7(b)に示したように、ある極小値から、当該極小値の後(グラフでは右側)に隣り合う極大値までの間の解析フレーム画像(選択範囲r)を特定解析フレーム画像として自動選択する。
なお、図7に示したように、一のグラフに対して選択範囲r又はrを複数周期分設定することが可能である。
なお、rとrのどちらの選択範囲とするかは、ユーザーが解析フレーム画像の選択パターンを設定する際に併せて設定できるようにしてもよいし、制御部11が自動的に判断するようにしてもよい。
診断対象が肺野である場合には、選択範囲rは呼気中、選択範囲は吸気中の動きを表すことになる。一方、診断対象が心臓である場合には、選択範囲r又はrは心臓収縮中又は心臓膨張中の動きを表すことになる。
複数の呼気中(吸気中)あるいは心臓膨張中(心臓収縮中)の解析フレーム画像の結果に対して演算処理を加えることによって、呼吸や心拍のバラつきなどのノイズ影響を低減させることができ、血液が肺野に広がっていく心収縮中の様子、心臓に血液が戻っていく心膨張中の様子、吸気で肺が広がっていく様子、呼気で肺が縮んでいく様子等を精度よく把握することができる。
また、例えば閉塞性疾患患者の息の吐きにくさの評価、拘束性疾患や肥満の方の息の吸い難さの評価、心臓の動き方が不連続になる不整脈などの評価にも有用である。
選択パターン(3)が選択されると、放射線画像処理装置1は、関心領域R内内の平均値の時間方向の変化を数値化(グラフ化)する。そして、図8(a)に示したように、ある極大値となる時間tから所定時間tだけ前に撮影された解析フレーム画像から、当該極大値となる時間から所定時間tだけ後に撮影された解析フレーム画像まで(t−t〜t+t:選択範囲r)、又は図8(b)に示したように、ある極小値となる時間tから所定時間tだけ前に撮影された解析フレーム画像から、当該極小値となる時間から所定時間tだけ後に撮影された解析フレーム画像まで(t−t〜t+t:選択範囲r)を特定解析フレーム画像として自動選択する。所定時間t,tは一周期(呼吸周期、心拍周期)の1/4以下とするのが好ましく、より好ましくは1/8以下とする。
なお、図8には、所定時間tとtが同じ長さである場合を例示したが、tとtは異なっていてもよい。
また、時間ではなく、極大からの解析フレーム画像数で範囲を規定するようにしてもよい。
更に、図8に示したように、一のグラフに対して選択範囲r又はrを複数周期分設定することが可能である。
呼吸周期(心拍周期)は患者ごとに、スパイロメトリやモストグラフ、心電図等で計測した値を用いても良いし、各国、身長や性別、年齢などで定められた平均的な値を用いても良い。
診断対象が肺野である場合には、極大及びその周辺は吸気位辺り、極小及びその周辺は呼気位辺りの動きを表すことになる。一方、診断対象が心臓である場合には、心臓収縮期辺り又は心臓膨張期辺りの動きを表すことになる。
複数の呼気位(吸気位)又は心臓膨張期(心臓収縮期)付近の解析フレーム画像の結果に対して演算処理を加えることによって、血液が肺野に広がり切った心収縮期の様子、心臓に血液が戻り切った心膨張期の様子、呼吸により肺が広がりきった様子、肺が縮み切った様子等を精度よく把握することができる。
心臓から肺に血が広がる速度は速い為、肺野内のどこに血液が流れて行ったかを把握するのは通常の解析結果(動画)を見ているだけでは容易ではない。一方、全ての解析フレーム画像に対して演算処理を加えるとアーチファクトも多く付加されるデメリットがある。このため、血液が肺野に広がり切る前後の解析フレーム画像に対して演算処理を加えることで、アーチファクトの少ない肺血流分布画像を得ることが可能となる。つまり、選択パターン(3)によって得られた特定解析フレーム画像は、肺血流に関する疾患の把握に有用である。例えば、肺塞栓など局所的な肺血流機能欠損領域の把握や、心機能が低下する心不全や心肥大の評価に有用である。一方で、(3)によって得られた特定解析フレーム画像は同様な理由でアーチファクトの少ない換気機能画像となる為、局所的な肺換気機能欠損領域の把握や、横隔膜がどこまで広がっているか等、肺の広がり方の評価にも有用である。
なお、特定解析フレーム画像を複数周期分選択する場合には、サマリー表示領域15Aに各選択範囲の演算結果画像を並列表示するようにしてもよい。
こうすることで、どの呼吸(どの心拍)が最も正常/異常かを判断し易くなる。
なお、選択パターン(3)においては、所定時間t,tのいずれか一方を0としてもよい。すなわち、図9(a)に示したように、ある極大値となる時間tから所定時間tだけ後に撮影された解析フレーム画像まで(極大〜t+t:選択範囲r)、ある極小値となる時間tから所定時間tだけ後に撮影された解析フレーム画像まで(極小〜t+t:選択範囲r)、ある極大値となる時間tから所定時間tだけ前に撮影された解析フレーム画像から当該極大まで(t−t〜極大:選択範囲r)、又はある極小値となる時間tから所定時間tだけ前に撮影された解析フレーム画像から当該極小まで(t−t〜極小:選択範囲r10)を特定解析フレーム画像として自動選択するようにしてもよい。
なお、本実施形態では、関心領域R内の画素信号値の平均値の時間方向の変化を数値化(グラフ化)し、そこから特定解析フレーム画像を選択するようにしたが、他の統計値(最大値、最小値、最頻値、積算値、標準偏差等)を用いて特定解析フレーム画像を選択するようにしてもよい。
また、横隔膜移動量の時間変化や、肺野面積の時間変化、肋間距離の変化量の時間変化、肺野領域の水平方向の広がりの時間変化、心臓面積の時間変化、大動脈弓の位置の時間変化等の形態情報の時間方向変化を元に、特定解析フレーム画像を選択するようにしても良い。
撮影対象が、例えば肺野や心臓領域にて脂肪や筋肉などが不均一に分布している体型の患者(肥満の方や女性等に多い)であったり、撮影中に体が上下左右に大きく動いてしまう(体動が大きい)患者であったりする場合、画素値を用いた時間変化情報が正確に呼吸状態や心拍状態を反映していない場合がある。しかし、形態情報の時間方向の変化を用いることで、このような患者に対しても特定解析フレーム画像を精度よく選択することができるようになる。
<第3実施形態>
次に、本発明の第3実施形態に係る放射線画像処理装置1Bについて説明する。図10は放射線画像処理装置1Bにおける画像処理の概念図であり、図11,12は処理結果の表示例である。
なお、図10には、肺野領域の画像を例示したが、本実施形態に係る放射線画像処理装置1Bは、肺野以外の部位を撮影した動画についても処理・表示することが可能である。
第1,第2実施形態では、サマリー表示領域15Aに表示する内容を画像のみとしたが、本実施形態に係る放射線画像処理装置1Bは、所定の定量値をサマリー表示領域15Aに表示するようになっている。
このため、本実施形態に係る放射線画像処理装置1Bは、実行する制御(記憶部13Bの記憶内容)が第1,第2実施形態と異なっている。なお、記憶部13Bの記憶内容を除く他の構成は第1,第2実施形態と同様である。
放射線画像処理装置1Bの制御部11は、複数の演算信号値を算出した後、演算結果画像の関心領域R内に位置することになる複数の画素の信号値に基づいて定量値Qを導出するようになっている。すなわち、制御部11は、本発明における定量値算出手段をなす。なお、関心領域Rは、本発明における特定領域である。以下、関心領域R内に位置することになる画素を領域内画素と称する。なお、関心領域Rは複数設定されていても良く、その場合は、各々の関心領域Rに対応した定量値Qを各々個別に導出しても良いし、複数設定された各々の関心領域Rをまとめて一つとみなして、一つの定量値Qを導出しても良い。
定量値Qの導出の仕方としては、例えば、最大値処理、最小値処理、Pk−Pk値処理、標準偏差処理、平均値処理、最頻値処理、積算値処理等がある。
最大(小)値処理は、演算結果画像の関心領域R内に位置することになる複数の領域内画素の各信号値の中から最大(小)値を抽出して定量値Qとする処理である。この処理により得られた定量値Qを用いれば、経過観察(換気、血流の悪化・良化等)や、左右肺機能の比較等を行うことができる。特に、低い計算コストで、肺換気機能や肺血流機能、心機能に関する特徴量の最大(小)値を把握できるので、例えば、治療前後などの経過観察(肺活量の変化や血圧の変化を見る)に有用である。
最頻値処理は、複数の領域内画素の各信号値の中から最頻値を抽出して定量値Qとする処理である。この処理により得られた定量値Qを用いれば、経過観察や左右肺機能の比較等を行うことができる。特に、演算結果画像上に特異的に発生する画欠等のアーチファクトの影響を低減できる為、関心領域R内の肺換気機能や肺血流機能、心機能を精度よく把握することができる。
Pk−Pk値処理は、複数の領域内画素の各信号値の中から最大値と最小値を抽出してその差分を算出し、当該領域内画素の数に対する最大値と最小値との差の割合を算出して定量値Qとする処理である。この処理により得られた定量値Qを用いれば、経過観察や左右肺機能の比較等を行うことができる。特に、肺換気機能や肺血流機能、心機能に関する特徴量の変化量を把握できるので、関心領域R内における上記各機能の特徴量のバラつき(局所的に肺換気機能や肺血流機能が異常を示す疾患)を把握しやすくなる。
標準偏差処理は、複数の領域内画素の各信号値の標準偏差を算出して定量値Qとする処理である。この処理により得られた定量値Qを用いれば、経過観察や左右肺機能の比較等を行うことができる。特に、低い計算コストで、肺換気機能や肺血流機能、心機能に関する特徴量の変化量を把握できるので、関心領域R内の上記各機能のバラつきを把握しやすくなる。
平均値処理は、複数の領域内画素の各信号値の平均値を算出して定量値Qとする処理である。この処理により得られた定量値Qを用いれば、経過観察や左右肺機能比較等を行うことができる。特に、低い計算コストで、関心領域R内の肺換気機能や肺血流機能、心機能を把握することができる。
積算値処理は、複数の領域内画素の各信号値の積算値を算出して定量値Qとする処理である。この処理により得られた定量値Qを用いれば、経過観察や左右肺機能比較等を行うことができる。特に、関心領域Rの取り方にわずかな変化があっても算出される定量値Qに変化が少ない為、測定誤差が生じにくくなる。また、関心領域R内のトータルでの肺機能、心機能を把握しやすくなる。
なお、定量値Qは、関心領域R内の領域内画素の信号値に基づいて導出するのではなく、演算結果画像の全画素のうち、信号値が所定閾値以上、所定閾値以下又は所定範囲内にある画素の個数(面積)としてもよい。
このようにすれば、例えば、肺全体(全体肺面積)に対して、正常又は異常な肺領域が全体の何割あるか等の情報を容易に把握でき、肺の局所的な機能分布(どこに肺機能の低下、欠損があるか)を評価しやすくなる。
本実施形態において、関心領域Rとして設定できるのは、全画像領域、肺野内、心臓領域内、肺動脈内、大動脈内、肺気管支内である。
全画像領域を関心領域Rとした定量値Qを用いれば、計算コストをかけずに、手軽に大局的な肺換気機能、肺血流機能、心機能を把握することができる。
肺野内を関心領域Rとした定量値Qを用いれば、肺換気機能(呼吸に伴う濃度の変化、横隔膜の動き、肋骨の動き、鎖骨の動き、肩甲骨の動きなど)、肺血流機能(心拍に伴う濃度の変化、血管の動き)を把握することができる。
心臓領域内を関心領域Rとした定量値Qを用いれば、心機能(心壁の移動量(移動速度、加速度)、心拍に伴う濃度の変化量(濃度変化速度、濃度変化加速度))を把握することができる。
肺動脈内を関心領域Rとした定量値Qを用いれば、肺血流機能、心機能を把握することができる。
大動脈内を関心領域Rとした定量値Qを用いれば、体循環の血流機能(心拍に伴う濃度の変化、血管の動き)、心機能を把握することができる。
肺気管支内を関心領域Rとした定量値Qを用いれば、肺換気機能の把握、気道の肥厚の程度、気道収縮の程度、無気肺/胸水/腫瘤、のう胞/肺気腫の有無などを把握することができる。
なお、定量値Qの表示は、図10(a)に示したように単独で行ってもよいし、図10(b),(c)に示したように、画像と併せて行うようにしてもよい。
画像I,Iと併せて表示する場合には、図10(b)に示したように並列表示する他、切替表示、重畳表示することも可能である。
また、図10(c)に示したように、画像I,Iと関心領域Rを重畳表示し、更にその中に定量値Qを表示するようにしてもよい。このようにすれば、左右の肺機能を直感的に把握することが可能となる。
また、制御部11は、種類の異なる複数のフレーム画像群からそれぞれフレーム画像を選択(サマリー表示領域15Aにドラッグ及びドロップ等)した場合は、各動画についての演算結果画像Iを生成するだけでなく、各動画についての定量値Qも算出する。
このとき、換気解析動画(解析結果動画)について算出した定量値Qと、血流解析動画(解析結果動画)について算出した定量値Qを用いた特徴量(比や差分など)を算出し、図11に示したようにV/Q比(換気血流比)を示すV/Q画像Vを表示するようにしても良い。
このようにすれば、ユーザーが正常異常の判断を行いやすくなるだけでなく、経時的な比較が可能となる。
また、図12(a)に示したように、左右の肺ごとに定量値Qを算出し、さらに各定量値Qの左右肺比率を算出して定量値Qとともに表示するようにしても良い。こうすることで、左右肺の肺機能の不均等がわかり易くなり、正常と異常の把握が行いやすくなる。
また、心臓領域を診断対象とする場合には、図12(b)に示したように、心室、心房、右心系、左心系ごとに定量値Q及び比率を算出し、表示するようにしてもよい。こうすることで、心機能の不均衡がわかり易くなり、正常と異常の把握が行いやすくなる。
また、本実施形態において算出した定量値Qについては、演算結果画像Iを保存する際、画像のヘッダーに保存するようにしても良い。このようにすれば、演算結果画像Iを再度読み込む際、演算結果画像Iのファイル名一覧でヘッダー情報の有無を見ることにより、必要な場合のみ定量値Qを表示させることができる。
また、画像とは別データ(テキストやCSVファイルなど)として保存しても良い。そうすることで、定量値Qのみを再度見直すことが容易になるし、不要に画像を毎回開く手間が減り、経過観察等の作業効率が向上する。
(変形例1)
なお、上記実施形態においては、サマリー表示領域15Aに選択した複数の解析フレーム画像に基づく演算結果画像や定量値Qを表示するようにしたが、図13に示したように、単に選択(サマリー表示領域15Aにドラッグ及びドロップ)した複数のフレーム画像のみを表示するようにしてもよい。
なお、本変形例においても、上記実施形態と同様、選択するフレーム画像は必ずしも連続のフレーム画像である必要はない。また、サマリー表示領域15Aへの表示は、フレーム画像の選択後自動的に行ってもよいし、ユーザーの操作(例えばOKボタンをクリックする等)に基づいて行うようにしてもよい。
また、異なる複数種類のフレーム画像の一覧からそれぞれフレーム画像を選択した場合は、図14(a)に示したように重畳画像Iの表示(ここでは、元フレーム画像と解析フレーム画像の重畳表示)の他、図14(b)に示した並列表示(ここでは、解析フレーム画像I、元フレーム画像I及び重畳画像Iの並列表示)、図14(c)に示した切替表示(ここでは、解析フレーム画像I、元フレーム画像I及び重畳画像Iの切替表示)を行う。
なお、異なる複数種類の動画の中からそれぞれフレーム画像を選択する際、それぞれのフレーム画像数は同一であっても異なっていても良い。
本変形例のようにすれば、ユーザーにとって関心の高いフレーム画像のみがサマリー表示領域15Aに表示される為、ユーザーは読影に注力しやすくなり、読影精度が向上する。
特に、1呼吸分など呼吸周期の整数倍に相当するフレーム画像のみをサマリー表示領域15Aに表示することで、呼吸が正しく行われているかを評価しやすくなる。
また、1心拍分など、心拍周期の整数倍に相当するフレーム画像のみをサマリー表示領域15Aに表示することで、心臓が正常に動いているかを評価しやすくなる。
また、吸気始めあるいは、呼気始めに相当するフレーム画像のみをサマリー表示領域15Aに表示することで、呼吸器疾患症例に対しての、息の吸い難さや息の吐きにくさを評価しやすくなる。
また、心膨張あるいは心収縮に相当するフレーム画像のみをサマリー表示領域15Aに表示することで、心疾患症例に対しての、血液の拍入のし易さ、血液の拍出のし易さ、不整脈の有無や心肥大の有無を評価しやすくなる。
また、肺機能が正常あるいは異常と判断したフレーム画像群をサマリー表示領域15Aに表示させたうえで、表示部15に元々表示されているフレーム画像の一覧と比較することで、表示部15に元々表示されているフレーム画像群が総じて正常なのか異常なのかを判断しやすくなる。すなわち、画像表示部に元々表示されているフレーム画像群の読影効率が向上する。
なお、図15(a)に示したように、換気解析結果と、血流解析結果をサマリー表示領域15Aに重畳表示することで、容易に換気機能と血流機能を比較する画像(V/Q画像)と同等の画像を得ることができる。
このようにすれば、ユーザーが選択したフレーム画像のみを用いて換気機能と血流機能の比較が行える為、例えばアーチファクトの影響をあまり受けていないV/Q画像を得ることができ、診断精度が向上する。
また、図15(b)に示したように、換気解析結果と、血流解析結果、これらの重畳画像及び元画像の並列表示や、切替表示も可能である。
(変形例2)
また、上記実施形態では、表示部15に表示されたフレーム画像の一覧の中から、サマリー表示領域15Aに解析結果動画等を表示するためのフレーム画像を選択するようにしていたが、図16に示したように、表示部15において動画を再生し、動画から必要な画像を直接ドラックしサマリー表示領域15Aにドロップするようにしても良い。
また、図17に示したように、各フレーム画像にチェックボックスBを設け、チェックCを入れたフレーム画像がサマリー表示領域15Aに表示されるようにしても良い。このようにすれば、サマリー表示領域15Aに表示されたフレーム画像又は動画を見ながら、チェックCを入れるフレーム画像を選べるため、ユーザーが選択したフレーム画像が妥当かどうか(診断にふさわしいフレーム画像かどうか)判断を素早く行うことが可能となる。
(変形例3)
また、図18に示したように、選択されサマリー表示領域15Aに表示されるフレーム画像については、選択されていないフレーム画像や動画とは異なるルックアップテーブル(Lookup table;LUT)を適用した画像Iを表示してもよい。具体的には、濃度、コントラスト、カラーマップの変更し表示する。
このようにすれば、ユーザーが関心の高い画像が、ユーザーの関心の高い色づけでサマリー表示領域15Aに表示される為、ユーザーは読影に注力しやすくなり、読影精度が向上する。
1,1A,1B 放射線画像処理装置
11 制御部(画像取得手段、解析手段、選択手段、演算処理手段、画像表示手段)
12 通信部
13,13A,13B 記憶部
14 操作部
15 表示部
15A サマリー表示領域
16 バス
元動画、元フレーム画像
解析結果動画、解析フレーム画像
演算結果画像
重畳画像
R 関心領域
〜r10 (フレーム画像の)選択範囲
t 時間
,t 所定時間

Claims (10)

  1. 画像を表示可能な表示部と、
    複数のフレーム画像からなる放射線動画の動画データを取得する画像取得手段と、
    前記画像取得手段が取得した動画データに所定の解析処理を加え、複数の解析フレーム画像からなる解析結果動画を生成する解析手段と、
    前記解析手段が生成した解析結果動画を構成する複数の解析フレーム画像の中から複数の特定解析フレーム画像を選択する選択手段と、
    前記選択手段によって選択された複数の特定解析フレーム画像における、共通する座標上にそれぞれ位置する各共通座標画素の信号値に基づいて一の演算信号値を導出する演算処理手段と、
    前記演算処理手段が座標ごとに生成した複数の演算信号値を用いて前記表示部に演算結果画像を表示する画像表示手段と、を備えることを特徴とする放射線画像処理装置。
  2. 前記画像表示手段は、
    前記表示部の表示領域の少なくとも一部にサマリー表示領域を形成し、
    前記サマリー表示領域に前記演算結果画像を表示させることが可能であることを特徴とする請求項1に記載の放射線画像処理装置。
  3. ユーザーが操作可能な操作部を有し、
    前記選択手段は、ユーザーが前記操作部に対して行った操作に基づいて前記特定解析フレーム画像を選択することを特徴とする請求項1又は2に記載の放射線画像処理装置。
  4. 前記選択手段は、予め設定された選択パターンに基づいて前記特定解析フレーム画像を自動的に選択することを特徴とする請求項1又は2に記載の放射線画像処理装置。
  5. ユーザーが操作可能な操作部を有し、
    前記選択手段は、
    ユーザーが前記操作部に対して行った操作に基づいて、予め用意された複数の選択パターンの中から何れかの選択パターンを設定し、
    設定した選択パターンに基づいて前記特定解析フレーム画像を自動的に選択することを特徴とする請求項4に記載の放射線画像処理装置。
  6. 前記演算処理手段は、前記複数の演算信号値を算出した後、前記演算結果画像の特定領域内に位置することになる複数の領域内画素の信号値に基づいて所定の定量値を導出する定量値算出手段を備え、
    前記画像表示手段は、前記定量値算出手段が算出した定量値を前記サマリー表示領域15Aに表示させることが可能であることを特徴とする請求項1から5のいずれか一項に記載の放射線画像処理装置。
  7. 前記演算処理手段は、前記複数の特定解析フレーム画像における、前記各共通座標画素の信号値の中から最大値、最小値もしくは最頻値を抽出する、当該共通座標画素の信号値における、選択した特定解析フレーム画像の数に対する最大値と最小値との差の割合、標準偏差、平均値もしくは積算値を算出する、又は前後する2枚の特定解析フレーム画像における各共通座標画素の信号値の差分を算出し、得られた各差分の中の最大値もしくは最小値を抽出することにより、前記一の演算信号値を導出することを特徴とする請求項1から6のいずれか一項に記載の放射線画像処理装置。
  8. 前記定量値算出手段は、前記複数の領域内画素の信号値の中から最大値、最小値もしくは最頻値を抽出する、又は当該複数の領域内画素の信号値における、当該領域内画素の数に対する最大値と最小値との差の割合、標準偏差、平均値もしくは積算値を算出することにより、前記定量値を導出することを特徴とする請求項6又は7に記載の放射線画像処理装置。
  9. 前記特定領域は、肺野内、心臓領域内、肺動脈内、大動脈内又は肺気管支内であることを特徴とする請求項8に記載の放射線画像処理装置。
  10. 放射線を照射する放射線照射装置と、
    前記放射線照射装置から放射線を受けることで放射線画像の画像データを生成する放射線画像撮影装置と、を備え、
    前記放射線照射装置による放射線の照射と、放射線画像撮影装置による画像データの読出しを連続的に繰り返すことにより、放射線動画の動画データを生成することが可能な放射線画像撮影システムであって、
    前記放射線画像撮影装置と直接又は間接的に通信可能に接続された請求項1から9のいずれか一項に記載の放射線画像処理装置と、を備えることを特徴とする放射線画像撮影システム。
JP2017056991A 2017-03-23 2017-03-23 放射線画像処理装置及び放射線画像撮影システム Active JP6805918B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2017056991A JP6805918B2 (ja) 2017-03-23 2017-03-23 放射線画像処理装置及び放射線画像撮影システム
US15/926,328 US10878560B2 (en) 2017-03-23 2018-03-20 Radiation image processing apparatus and radiation image capturing system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017056991A JP6805918B2 (ja) 2017-03-23 2017-03-23 放射線画像処理装置及び放射線画像撮影システム

Publications (2)

Publication Number Publication Date
JP2018157968A true JP2018157968A (ja) 2018-10-11
JP6805918B2 JP6805918B2 (ja) 2020-12-23

Family

ID=63581191

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017056991A Active JP6805918B2 (ja) 2017-03-23 2017-03-23 放射線画像処理装置及び放射線画像撮影システム

Country Status (2)

Country Link
US (1) US10878560B2 (ja)
JP (1) JP6805918B2 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020080949A (ja) * 2018-11-16 2020-06-04 コニカミノルタ株式会社 画像処理装置及びプログラム
JP2020151232A (ja) * 2019-03-20 2020-09-24 コニカミノルタ株式会社 動態解析装置、動態解析システム、予測率算出方法及びプログラム
JP2020168176A (ja) * 2019-04-03 2020-10-15 コニカミノルタ株式会社 画像解析装置、画像解析システム及びプログラム
JP2020171427A (ja) * 2019-04-09 2020-10-22 コニカミノルタ株式会社 動態画像解析装置、動態画像解析システム及びプログラム
JP2021146013A (ja) * 2020-03-23 2021-09-27 コニカミノルタ株式会社 プログラム、動態解析装置及び動態解析システム

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6805918B2 (ja) * 2017-03-23 2020-12-23 コニカミノルタ株式会社 放射線画像処理装置及び放射線画像撮影システム
CN109978860B (zh) * 2019-03-27 2021-07-16 北京青燕祥云科技有限公司 胸水判别方法、装置和计算机设备
CN110490925A (zh) * 2019-07-19 2019-11-22 东北大学 一种重构区域划分的方法和装置
JP7468427B2 (ja) * 2021-03-25 2024-04-16 コニカミノルタ株式会社 画像表示装置、制御プログラム、画像表示方法及び画像表示システム

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014192504A1 (ja) * 2013-05-31 2014-12-04 コニカミノルタ株式会社 画像処理装置及びプログラム
JP2015043880A (ja) * 2013-08-28 2015-03-12 コニカミノルタ株式会社 胸部診断支援システム
JP2016214725A (ja) * 2015-05-25 2016-12-22 コニカミノルタ株式会社 動態解析システム
JP2017018681A (ja) * 2016-10-17 2017-01-26 コニカミノルタ株式会社 動態解析システム

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005106474A2 (en) * 2004-04-14 2005-11-10 Edda Technology, Inc. Liver disease diagnosis system, method and graphical user interface
JP5408873B2 (ja) * 2004-05-26 2014-02-05 ベルス・メステヒニーク・ゲーエムベーハー 座標測定装置におけるx線感知装置の校正方法
JP5625799B2 (ja) 2010-11-22 2014-11-19 コニカミノルタ株式会社 動態診断支援情報生成システム
JP2012110400A (ja) * 2010-11-22 2012-06-14 Konica Minolta Medical & Graphic Inc 動態診断支援情報生成システム
JP2012110466A (ja) 2010-11-24 2012-06-14 Canon Inc X線透視撮影装置、方法およびプログラム
WO2014185197A1 (ja) * 2013-05-16 2014-11-20 コニカミノルタ株式会社 画像処理装置及びプログラム
CN105246407A (zh) * 2013-05-28 2016-01-13 柯尼卡美能达株式会社 图像处理装置以及程序
JP6805918B2 (ja) * 2017-03-23 2020-12-23 コニカミノルタ株式会社 放射線画像処理装置及び放射線画像撮影システム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014192504A1 (ja) * 2013-05-31 2014-12-04 コニカミノルタ株式会社 画像処理装置及びプログラム
JP2015043880A (ja) * 2013-08-28 2015-03-12 コニカミノルタ株式会社 胸部診断支援システム
JP2016214725A (ja) * 2015-05-25 2016-12-22 コニカミノルタ株式会社 動態解析システム
JP2017018681A (ja) * 2016-10-17 2017-01-26 コニカミノルタ株式会社 動態解析システム

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020080949A (ja) * 2018-11-16 2020-06-04 コニカミノルタ株式会社 画像処理装置及びプログラム
JP7135759B2 (ja) 2018-11-16 2022-09-13 コニカミノルタ株式会社 画像処理装置及びプログラム
JP2020151232A (ja) * 2019-03-20 2020-09-24 コニカミノルタ株式会社 動態解析装置、動態解析システム、予測率算出方法及びプログラム
JP7074101B2 (ja) 2019-03-20 2022-05-24 コニカミノルタ株式会社 動態解析装置、動態解析システム、予測率算出方法及びプログラム
JP2020168176A (ja) * 2019-04-03 2020-10-15 コニカミノルタ株式会社 画像解析装置、画像解析システム及びプログラム
JP7156142B2 (ja) 2019-04-03 2022-10-19 コニカミノルタ株式会社 画像解析装置、画像解析システム及びプログラム
US11730907B2 (en) 2019-04-03 2023-08-22 Konica Minolta, Inc. Image analysis apparatus, image analysis system, and storage medium
JP2020171427A (ja) * 2019-04-09 2020-10-22 コニカミノルタ株式会社 動態画像解析装置、動態画像解析システム及びプログラム
JP7255329B2 (ja) 2019-04-09 2023-04-11 コニカミノルタ株式会社 動態画像解析装置、動態画像解析システム及びプログラム
JP2021146013A (ja) * 2020-03-23 2021-09-27 コニカミノルタ株式会社 プログラム、動態解析装置及び動態解析システム
JP7487506B2 (ja) 2020-03-23 2024-05-21 コニカミノルタ株式会社 プログラム、動態解析装置、動態解析システム及び情報処理装置

Also Published As

Publication number Publication date
US10878560B2 (en) 2020-12-29
US20180276816A1 (en) 2018-09-27
JP6805918B2 (ja) 2020-12-23

Similar Documents

Publication Publication Date Title
JP6805918B2 (ja) 放射線画像処理装置及び放射線画像撮影システム
EP3264985B1 (en) Tomography imaging apparatus and method of reconstructing tomography image
CN102125439B (zh) 图像处理装置、x射线计算机断层摄像装置以及图像处理方法
JP6597548B2 (ja) 動態解析システム
JP6418091B2 (ja) 胸部画像表示システム及び画像処理装置
JP2009153677A (ja) 動態画像処理システム
US11189025B2 (en) Dynamic image analysis apparatus, dynamic image analysis method, and recording medium
JP2009153678A (ja) 動態画像処理システム
JP6690774B2 (ja) 動態解析システム、プログラム及び動態解析装置
US20230290023A1 (en) Image processing apparatus, image processing system and storage medium
US10918351B2 (en) Radiographic-image processing apparatus
JP2020171476A (ja) 動態解析装置、動態解析システム及びプログラム
US20200379636A1 (en) Image processing apparatus, image processing method, and non-transitory computer-readable storage medium
JP5215050B2 (ja) 画像処理装置
US11484221B2 (en) Dynamic analysis apparatus, dynamic analysis system, expected rate calculation method, and recording medium
JP7143747B2 (ja) 画像表示装置、画像表示方法及び画像表示プログラム
JP6852545B2 (ja) 画像表示システム及び画像処理装置
JP2019180899A (ja) 医用画像処理装置
JP6930515B2 (ja) 画像表示装置、画像表示方法及び画像表示プログラム
US20230263495A1 (en) Pulmonary embolism diagnosis support apparatus, pulmonary embolism diagnosis support method, and storage medium
JP6888721B2 (ja) 動態画像処理装置、動態画像処理プログラム及び動態画像処理方法
US20240120065A1 (en) Medical system for ventilation and imaging
EP4036926A1 (en) Medical system for ventilation and imaging
JP2017217047A (ja) 画像表示システム
JP2023121104A (ja) 肺塞栓症診断支援装置、肺塞栓症診断支援方法及びプログラム

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190620

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200630

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200721

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200903

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20201117

R150 Certificate of patent or registration of utility model

Ref document number: 6805918

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150