WO2018083844A1 - X線透視装置およびx線透視方法 - Google Patents

X線透視装置およびx線透視方法 Download PDF

Info

Publication number
WO2018083844A1
WO2018083844A1 PCT/JP2017/026330 JP2017026330W WO2018083844A1 WO 2018083844 A1 WO2018083844 A1 WO 2018083844A1 JP 2017026330 W JP2017026330 W JP 2017026330W WO 2018083844 A1 WO2018083844 A1 WO 2018083844A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
ray
specific part
subject
created
Prior art date
Application number
PCT/JP2017/026330
Other languages
English (en)
French (fr)
Inventor
森 慎一郎
▲高▼橋 渉
Original Assignee
株式会社島津製作所
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 株式会社島津製作所 filed Critical 株式会社島津製作所
Priority to US16/347,048 priority Critical patent/US11045663B2/en
Publication of WO2018083844A1 publication Critical patent/WO2018083844A1/ja

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048Monitoring, verifying, controlling systems and methods
    • A61N5/1049Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/42Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4208Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
    • A61B6/4225Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using image intensifiers
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/44Constructional features of apparatus for radiation diagnosis
    • A61B6/4429Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/485Diagnostic techniques involving fluorescence X-ray imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/01Devices for producing movement of radiation source during therapy
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048Monitoring, verifying, controlling systems and methods
    • A61N5/1064Monitoring, verifying, controlling systems and methods for adjusting radiation treatment in response to monitoring
    • A61N5/1068Gating the beam as a function of a physiological signal
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • G06T7/246Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
    • G06T7/248Analysis of motion using feature-based methods, e.g. the tracking of corners or segments involving reference images or patches
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048Monitoring, verifying, controlling systems and methods
    • A61N5/1049Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
    • A61N2005/1061Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam using an x-ray imaging system having a separate imaging source
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10116X-ray image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10116X-ray image
    • G06T2207/10121Fluoroscopy
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10116X-ray image
    • G06T2207/10124Digitally reconstructed radiograph [DRR]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; Learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30096Tumor; Lesion

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biomedical Technology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Optics & Photonics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Physiology (AREA)
  • Multimedia (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Analysis (AREA)
  • Radiation-Therapy Devices (AREA)

Abstract

制御部30は、被検者の特定部位を含む画像に基づいて作成された複数のテンプレートを記憶する第1画像記憶部32と、被検者の特定部位を含む画像に基づいて作成された複数の正解画像を記憶する第2画像記憶部33とから構成される画像記憶部31と、複数の正解画像に基づいて機械学習により識別器を作成する学習部34と、被検者の特定部位を含む画像を所定のフレームレートで収集した複数の画像を使用して識別器を利用した機械学習により特定部位を含む領域を選択する候補位置選択部35と、候補位置選択部35により選択された特定部位を含む領域に対して複数のテンプレートを利用してテンプレートマッチングを実行することにより特定部位の位置を検出する位置検出部36と、を備える。

Description

X線透視装置およびX線透視方法
 この発明は、X線管から照射され被検者を通過したX線をX線検出器により検出して被検者の特定部位を含む画像を取得し、この被検者の特定部位を含む画像から特定部位の位置を検出し、特定部位の動きを追跡するX線透視装置およびX線透視方法に関する。
 腫瘍などの患部に対してX線や陽子線等の放射線を治療ビームとして照射する放射線治療においては、放射線を患部に正確に照射する必要がある。しかしながら、被検者が体を動かしてしまう場合があるばかりではなく、患部自体に動きが生ずる場合がある。例えば、肺の近くの腫瘍は呼吸に基づき大きく移動する。このため、腫瘍のそばに球形状を有する金製のマーカを留置し、このマーカの位置をX線透視装置により検出して、治療放射線の照射を制御する構成を有する放射線治療装置が提案されている(特許文献1参照)。
 このような放射線治療装置においては、第1X線管と第1X線検出器から成る第1X線透視機構と、第2X線管と第2X線検出器から成る第2X線透視機構とを使用して体内に留置されたマーカを撮影し、第1X線透視機構による二次元の透視画像と第2X線透視機構による二次元の透視画像を利用して三次元の位置情報を得る。そして、連続してX線透視を行い、リアルタイムでマーカの三次元の位置情報を演算することで、移動を伴う部位のマーカを高精度で検出する。そして、検出されたマーカの位置情報に基づいて治療放射線の照射を制御することで、腫瘍の動きに応じた高精度の放射線照射を実行することが可能となる。このマーカの位置情報を得るときには、テンプレートを利用したテンプレートマッチングが実行される。
 ところで、上述したようにマーカを利用して腫瘍の動きを検出するためには、被検者の体内に、予めマーカを留置する必要がある。一方、近年、患者の腫瘍の領域などの特定部位をマーカのかわりに使用することで、マーカの留置を省略するマーカレストラッキングと呼称される方法も提案されている。
特許第3053389号公報
 このX線透視装置は、球状以外の形状を有するマーカや特定部位をも検出し得るものではあるが、特に特定部位を検出するときに、検出精度をさらに向上させることが要請されている。また、多数のテンプレートを利用して類似度の評価を行う場合には、計算コストが大きくなり、リアルタイムでの処理が困難となるという問題が生ずる。
 この発明は上記課題を解決するためになされたものであり、特定部位の検出精度をさらに向上させ、また、被検者の特定部位の画像を使用してテンプレートマッチングを実行するときにも、より容易に計算を実行してリアルタイムで特定部位の位置を検出することが可能なX線透視方法およびX線透視装置を提供することを第1の目的とする。
 また、被検者の特定部位を球状のマーカのかわりに使用するマーカレストラッキングにおいては、被検者の呼吸位相により特定部位の認識される形状が異なることから、特定部位を単一の呼吸位相から撮影した画像を利用してテンプレートを作成しても、これだけでは正確なテンプレートマッチングは実行し得ない。このため、被検者の呼吸位相により特定部位の認識される形状が異なることに対応するため、放射線治療の直前に異なる呼吸位相で撮影された特定部位を含む画像をテンプレートとして複数記憶するとともに、一定時間毎に撮影される画像に対して、予め記憶された複数のテンプレートを利用してテンプレートマッチングを行う必要がある。
 このように複数のテンプレートを使用するときには、オペレータが放射線治療の直前に被検者をX線透視し、X線画像を見ながら呼吸位相毎の腫瘍の位置を指定してテンプレートを作成する必要がある。このため、複数のテンプレートを作成するために時間がかかり、患者にとって苦痛であるばかりではなく、治療のためのスループットの低下をもたらすという問題が生ずる。
 この発明は上記課題を解決するためになされたものであり、放射線治療に先だって、予めCT画像データに基づいて作成されたDRR画像を利用して機械学習を実行することにより、放射線治療の直前にテンプレートの作成に要していた時間をなくして、速やかに動体追跡を実行するとともに、より正確に特定部位の位置を検出することが可能なX線透視装置およびX線透視方法を提供することを第2の目的とする。
 また、画像認識により特定部位の位置を検出するときに、特定部位が骨部付近に配置されているときには、骨部のコントラストが特定部位のコントラストより大きいことから、機械学習を実行するときには、特定部位よりも骨部が学習されてしまい、特定部位の位置を正確に検出できないという問題が生ずる。
 この発明は上記課題を解決するためになされたものであり、特定部位が骨部付近に存在したとしても、正確に特定部位の位置を検出することが可能なX線透視装置およびX線透視方法を提供することを第3の目的とする。
 さらに、放射線治療においては、特定部位の位置が所定の位置にあるだけではなく、被検者の呼吸位相が特定の位相であるときにのみ治療ビームを照射することが好ましい。すなわち、被検者に対する治療計画は、被検者の呼吸位相が特定の位相であることを条件として作成されることから、特定部位の位置が所定の位置にあったとしても、呼吸位相が治療計画時とは異なった場合には、治療計画通りの照射線量とならないためである。ここで、特定部位は呼吸以外の蠕動運動等でも移動することから、特定部位が所定の位置に配置されたとしても、呼吸位相が同一であるという保証はない。
 この発明は上記課題を解決するためになされたものであり、測定部位の位置とともに被検者の呼吸位相をも特定することにより、特定部位が所定の位置にあり、かつ、被検者の呼吸位相が特定の位相となったときに治療ビームを照射することが可能なX線透視装置およびX線透視方法を提供することを第4の目的とする。
 第1の発明は、X線管と、前記X線管から照射され被検者を通過したX線を検出するX線検出器とを備え、前記被検者の特定部位を含む画像を収集することにより、前記特定部位の位置を検出し、前記特定部位の動きを追跡するX線透視装置であって、前記被検者の特定部位を含む画像に基づいて作成されたテンプレートマッチング用のテンプレートと、前記被検者の特定部位を含む画像に基づいて作成された機械学習用の正解画像とを記憶する画像記憶部と、前記画像記憶部に記憶された正解画像に基づいて、前記機械学習により識別器を作成する学習部と、前記被検者の特定部位を含む領域を透視して得た画像に対して、前記テンプレートを利用してテンプレートマッチングを実行するとともに、前記識別器を利用して識別を実行することにより、前記特定部位の位置を検出する位置検出部と、を備えたことを特徴とする。
 第2の発明は、前記位置検出部は、前記テンプレートマッチングと前記識別器を利用した識別とのいずれか一方を実行した後に他方を実行し、あるいは、両方を同時に実行する。
 第3の発明は、X線管と、前記X線管から照射され被検者を通過したX線を検出するX線検出器とを備え、前記被検者の特定部位を含む画像を収集することにより、前記特定部位の位置を検出し、前記特定部位の動きを追跡するX線透視装置であって、前記被検者の特定部位を含む画像に基づいて作成されたテンプレートマッチング用のテンプレートと、機械学習用の正解画像とを記憶する画像記憶部と、前記画像記憶部に記憶された正解画像に基づいて、前記機械学習により識別器を作成する学習部と、前記X線管と前記X線検出器により収集された画像を使用して、前記識別器を利用した識別により特定部位候補位置を選択する候補位置選択部と、前記候補位置選択部により選択された特定部位を含む領域に対して、前記画像記憶部に記憶されたテンプレートを利用してテンプレートマッチングを実行することにより、前記特定部位の位置を検出する位置検出部と、を備えたことを特徴とする。
 第4の発明は、前記機械学習は、SVM、Boostingまたはニューラルネットワークである。
 第5の発明は、前記テンプレートと前記正解画像とは、各々、前記特定部位を含む同一の画像に基づいて作成される。
 第6の発明は、X線管から照射され被検者を通過したX線をX線検出器により検出し、前記被検者の特定部位を含む画像を収集することにより、前記特定部位の位置を検出し、前記特定部位の動きを追跡するX線透視方法であって、前記被検者の特定部位を含む画像に基づいて作成されたテンプレートマッチング用のテンプレートと、前記被検者の特定部位を含む画像に基づいて作成された機械学習用の正解画像とを記憶する記憶工程と、前記記憶工程で記憶された正解画像に基づいて、前記機械学習により識別器を作成する学習工程と、前記X線管から照射され前記被検者を通過したX線を前記X線検出器により検出し、前記被検者の特定部位を含む画像を収集する透視工程と、前記透視工程で得た画像に対して、前記テンプレートを利用したテンプレートマッチングを実行するとともに、前記識別器を利用した識別を実行することにより、前記特定部位の位置を検出する位置検出工程と、を含むことを特徴とする。
 第7の発明は、X線管から照射され被検者を通過したX線をX線検出器により検出し、前記被検者の特定部位を含む画像を収集することにより、前記特定部位の位置を検出し、前記特定部位の動きを追跡するX線透視方法であって、前記被検者の特定部位を含む画像に基づいて作成されたテンプレートマッチング用のテンプレートと、機械学習用の正解画像とを記憶する画像記憶工程と、前記画像記憶工程で記憶された複数の正解画像に基づいて、前記機械学習により識別器を作成する学習工程と、前記X線管から照射され前記被検者を通過したX線を前記X線検出器により検出し、前記被検者の特定部位を含む画像を収集する透視工程と、前記透視工程で得た画像を使用して、前記識別器を利用した識別により前記特定部位の候補位置を選択する候補位置選択工程と、前記候補位置選択工程で選択された特定部位の候補位置に対して、前記画像記憶工程で記憶されたテンプレートを利用してテンプレートマッチングを実行することにより、前記特定部位の位置を検出する位置検出工程と、を含むことを特徴とする。
 第8の発明は、X線管から照射され被検者を通過したX線をX線検出器により検出して被検者の特定部位を含む画像を取得し、前記被検者の特定部位を含む画像から前記特定部位の位置を検出し、前記特定部位の動きを追跡するX線透視装置であって、治療計画記憶時に作成されたCT画像データに対して、前記被検者に対する前記X線管と前記X線検出器との幾何学的透視条件を模擬した仮想的透視投影を行うことにより、前記特定部位を含む複数のDRR画像を作成するDRR画像作成部と、前記DRR画像作成部で作成されたDRR画像に対して機械学習を実行することにより、前記特定部位を認識するための識別器を作成する識別器作成部と、前記X線管から照射され前記被検者を通過したX線を前記X線検出器により検出して得た前記被検者の特定部位を含む画像に対して前記識別器作成部で作成された識別器を使用して識別を行うことにより、前記特定部位の位置を検出する位置検出部と、を備えることを特徴とする。
 第9の発明は、前記DRR画像作成部は、前記幾何学的透視条件の投影座標および角度の少なくとも一方を含むパラメータを変化させ、あるいは、画像の回転、変形及び拡大縮小の少なくとも1つを含む画像処理を施して前記DRR画像を作成する。
 第10の発明は、前記DRR画像作成部は、作成後のDRR画像に対して、コントラスト変化、ノイズ付加及びエッジ強調の少なくとも1つを実行する。
 第11の発明は、前記DRR画像作成部は、治療計画記憶時に作成された連続する複数の呼吸位相における前記特定部位を含む領域の3次元のCT画像データ群からなる4次元CT画像データに基づいて、前記特定部位を含む複数のDRR画像を作成する。
 第12の発明は、前記DRR画像作成部は、治療計画記憶時に作成された連続する複数の呼吸位相における前記特定部位を含む領域の3次元のCT画像データ群からなる4次元CT画像データのうち、前記被検者に治療ビームが照射される呼吸位相を含む複数の呼吸位相のCT画像データに基づいて、前記特定部位を含む複数のDRR画像を作成する。
 第13の発明は、前記DRR画像作成部は、治療計画記憶時に作成されたCT画像データに基づいて、前記特定部位を含み、かつ、前記被検者の骨部が除去されたDRR画像を作成する。
 第14の発明は、前記DRR画像作成部は、治療計画記憶時に作成されたCT画像データにおけるCT値に基づいて前記被検者の骨部が除去されたCT画像データを作成し、この骨部が除去されたCT画像データに基づいて、前記特定部位を含み、かつ、前記被検者の骨部が除去されたDRR画像を作成する。
 第15の発明は、前記DRR画像作成部は、治療計画記憶時に作成されたCT画像データにおけるボクセル値を線積分してDRR画像を作成するときに、このボクセル値に基づいて、前記特定部位を含み、かつ、前記被検者の骨部が除去されたDRR画像を作成する。
 第16の発明は、前記X線管から照射され前記被検者を通過したX線を前記X線検出器により検出し、前記被検者の特定部位を含む画像を収集することにより、前記特定部位を含む複数のX線画像を作成するX線画像作成部と、前記X線画像作成部により作成されたX線画像を利用して、前記識別器作成部において作成した識別器を補正する補正部と、をさらに備える。
 第17の発明は、前記X線管から照射され前記被検者を通過したX線を前記X線検出器により検出し、前記被検者の特定部位を含む画像を収集することにより、前記特定部位を含む複数のX線画像を作成するX線画像作成部をさらに備え、前記識別器作成部は、前記DRR画像作成部で作成されたDRR画像と前記X線画像作成部で作成されたX線画像とに対して機械学習を実行することにより、前記特定部位を認識するための識別器を作成する。
 第18の発明は、前記識別器は、前記特定部位を前記被検者の呼吸位相毎に認識するための多クラス識別器であり、前記位置検出部は、前記X線管から照射され前記被検者を通過したX線を前記X線検出器により検出して得た前記被検者の特定部位を含む画像に対して前記識別器作成部で作成された多クラス識別器を使用して多クラス識別を行うことにより、前記特定部位の位置と前記被検者の呼吸位相とを検出する。
 第19の発明は、X線管から照射され被検者を通過したX線をX線検出器により検出して被検者の特定部位を含む画像を取得し、前記被検者の特定部位を含む画像から前記特定部位の位置を検出し、前記特定部位の動きを追跡するX線透視方法であって、治療計画記憶時に作成されたCT画像データに基づいて、前記特定部位を含む複数のDRR画像を作成するDRR画像作成工程と、前記DRR画像作成工程において作成されたDRR画像に対して機械学習を実行することにより、前記特定部位を認識するための識別器を作成する識別器作成工程と、前記X線管から照射され前記被検者を通過したX線を前記X線検出器により検出し、前記被検者の特定部位を含む画像を得る透視工程と、前記透視工程において得た画像に対して前記識別器作成工程で作成された識別器を使用して識別を行うことにより、前記特定部位の位置を検出する位置検出工程と、を含むことを特徴とする。
 第1から第7の発明によれば、テンプレートマッチングと機械学習とを併用することにより、特定部位の検出精度を向上させることが可能となる。
 第3および第7の発明によれば、被検者の特定部位を含む画像に基づいて作成されたテンプレートを利用したテンプレートマッチングを実行するときにも、より容易に計算を実行してリアルタイムで特定部位の位置を検出することが可能となる。
 第4の発明によれば、機械学習をより迅速に実行することが可能となる。
 第5の発明によれば、複数のテンプレートと正解画像とを、単一の工程により効率的に作成することが可能となる。
 第8から第19の発明によれば、放射線治療に先だって、予めCT画像データに基づいて作成されたDRR画像を利用して機械学習を実行することにより、従来放射線治療の直前にテンプレートの作成に要していた時間をなくすことができ、速やかに動体追跡を実行することが可能となる。
 第9の発明によれば、幾何学的透視条件の投影座標および角度を含むパラメータを変化させ、あるいは、画像の回転、変形、拡大縮小を含む画像処理を施すことから、特定部位が再現性なく移動または変形した場合においても、特定部位の位置を正確に検出することが可能となる。そして、大量のDRR画像を作成することができることから、各患者に対応したオーダーメイドの識別器を学習することができ、さらには、低フレームレートのDRR画像を使用する場合においても、また、特定部位に対して背景が映り込んだ場合においても、特定部位の位置を正確に検出することが可能となる。
 第10の発明によれば、作成後のDRR画像に対して、コントラスト変化、ノイズ付加、エッジ強調を実行することから、DRR画像とX線画像に画質の差異が生じた場合においても、特定部位の位置を正確に検出することが可能となる。
 第11および第12の発明によれば、連続する複数の呼吸位相における前記特定部位を含む領域の3次元のCT画像データ群からなる4次元CT画像データに基づいて特定部位を含む複数のDRR画像を作成することから、被検者の特定部位とともに移動しない部位を除外して機械学習を実行することが可能となり、正確な動体追跡を実行することが可能となる。
 第13から第15の発明によれば、被検者の骨部が除外されたDRR画像を作成することから、被検者の骨部を除外して機械学習を実行することが可能となり、正確な動体追跡を実行することが可能となる。
 第16の発明によれば、X線画像を利用して識別器作成部において作成した識別器を補正することから、多量のX線画像を使用することなく識別器を迅速に補正するとともに、より正確に特定部位の位置を検出することが可能となる。
 第17の発明によれば、CT画像データに基づいて作成されたDRR画像とX線撮影により作成されたX線画像とに対して機械学習を実行することにより、多量のX線画像を使用することなく識別器を迅速に作成するとともに、より正確に特定部位の位置を検出することが可能となる。
 第18の発明によれば、測定部位の位置とともに、そのときの被検者の呼吸位相をも特定することができる。このため、特定部位が所定の位置にあり、かつ、被検者の呼吸位相が特定の位相となったときに、治療ビームを照射することが可能となる。
この発明に係るX線透視装置を、放射線照射装置90とともに示す斜視図である。 この発明に係るX線透視装置の主要な制御系を示すブロック図である。 この発明に係るX線透視装置を使用した動体追跡動作を示すフローチャートである。 テンプレートと正解画像の作成工程を示す説明図である。 被検者Mに対するX線透視領域E1と候補位置E2とを模式的に示す説明図である。 この発明の第2実施形態に係るX線透視装置を使用した動体追跡動作を示すフローチャートである。 この発明の第3実施形態に係るX線透視装置の主要な制御系を示すブロック図である。 この発明の第3実施形態に係るX線透視装置を使用した動体追跡動作を示すフローチャートである。 この発明の第4実施形態に係るX線透視装置によるDRR作成工程を示すフローチャートである。 この発明の第5実施形態に係るX線透視装置によるDRR作成工程を示すフローチャートである。 この発明の第6実施形態に係るX線透視装置の主要な制御系を示すブロック図である。 この発明の第6実施形態に係るX線透視装置を使用した動体追跡動作を示すフローチャートである。 この発明の第7実施形態に係るX線透視装置の主要な制御系を示すブロック図である。 この発明の第7実施形態に係るX線透視装置を使用した動体追跡動作を示すフローチャートである。 この発明の第8実施形態に係るX線透視装置の主要な制御系を示すブロック図である。 この発明の第8実施形態に係るX線透視装置を使用した動体追跡動作を示すフローチャートである。 畳み込みニューラルネットワークの概念を説明するための説明図である。 補間に使用可能な周期関数を示すグラフである。 畳み込みニューラルネットワークを利用して特定部位の位置を認定する動作を示す模式図である。 畳み込みニューラルネットワークを利用して特定部位の位置を認定する動作を示す模式図である。
 以下、この発明の実施の形態を図面に基づいて説明する。図1は、この発明に係るX線透視装置を、放射線照射装置90とともに示す斜視図である。これらのX線透視装置と放射線照射装置90とにより、放射線治療装置が構成される。
 放射線照射装置90は、カウチとも呼称される検診台29上の被検者に対して放射線照射を行うものであり、治療室の床面に設置された基台91に対して揺動可能に設置されたガントリー92と、このガントリー92に配設された治療ビームを出射するヘッド93とを備える。この放射線照射装置90によれば、ガントリー92が基台91に対して揺動することにより、ヘッド93から照射される治療ビームの照射方向を変更することができる。このため、被検者における腫瘍等の患部に対して様々な方向から治療ビームを照射することが可能となる。
 この放射線照射装置90とともに使用されるX線透視装置は、被検者の患部の位置を特定する動体追跡を行うためのX線透視を実行するものである。すなわち、上述した放射線照射装置90を使用した放射線治療時においては、放射線を被検者の体動に伴って移動する患部に正確に照射する必要がある。このため、被検者における腫瘍等の特定の形状を有する部位を特定部位として予め登録し、この特定部位を連続的にX線透視して、特定部位の三次元の位置情報を演算することで、特定部位を高精度で検出する、所謂、動体追跡を行う構成となっている。このように、従来の被検者における患部付近にマーカを留置する代わりに、被検者における腫瘍等の特定部位の画像をマーカとして使用する動体追跡の手法は、マーカレストラッキングと呼称されている。
 このX線透視装置は、第1X線管11a、第2X線管11bと、第1フラットパネルディテクタ21a、第2フラットパネルディテクタ21bとを備える。第1X線管11aから照射されたX線は、検診台29上の被検者を透過した後、第1フラットパネルディテクタ21aにより検出される。第1X線管11aと第1フラットパネルディテクタ21aとは、第1X線撮影系を構成する。第2X線管11bから照射されたX線は、検診台29上の被検者を透過した後、第2フラットパネルディテクタ21bにより検出される。第2X線管11bと第2フラットパネルディテクタ21bとは、第2X線撮影系を構成する。
 図2は、この発明に係るX線透視装置の主要な制御系を示すブロック図である。
 このX線透視装置は、論理演算を実行するプロセッサーとしてのCPU、装置の制御に必要な動作プログラムが格納されたROM、制御時にデータ等が一時的にストアされるRAM等を備え、装置全体を制御する制御部30を備える。この制御部30は、上述した第1X線管11a、第2X線管11bと、第1フラットパネルディテクタ21a、第2フラットパネルディテクタ21bとに接続されている。
 この制御部30は、被検者の特定部位を含む画像に基づいて作成されたテンプレートマッチング用の1または複数のテンプレートを記憶する第1画像記憶部32と、被検者の特定部位を含む画像に基づいて作成された機械学習用の1または複数の正解画像を記憶する第2画像記憶部33とから構成される画像記憶部31を備える。また、この制御部30は、第2画像記憶部33に記憶した正解画像に基づいて機械学習により識別器(後述)を作成する学習部34と、被検者の特定部位を含む画像を所定のフレームレートで収集した複数の画像を使用して識別器を利用した機械学習により特定部位の1または複数の候補位置を選択する候補位置選択部35と、候補位置選択部35により選択された1または複数の候補位置に対して第1画像記憶部32に記憶した複数のテンプレートを利用してテンプレートマッチングを実行することにより特定部位の位置を検出する位置検出部36と、を備える。
 次に、以上のような構成を有するX線透視装置を使用することにより、被検者の体動に伴って移動する特定部位の位置を検出する動体追跡を実行する動作について説明する。図3は、この発明に係るX線透視装置を使用した動体追跡動作を示すフローチャートである。以下の動作は、2つのX線撮影系に対して同様に実行される。
 この発明に係るX線透視装置により動体追跡を実行するときには、最初に、画像記憶工程を実行する(ステップS1)。この画像記憶工程においては、上述した第1X線撮影系、第2X線撮影系を使用して、検診台29上に載置した被検者の特定部位を含む画像を撮影する。そして、この撮影により得たX線画像に基づいて、予め、テンプレートマッチング用の1または複数のテンプレートと、機械学習用の1または複数の正解画像とを作成する。作成されたテンプレートマッチング用の複数のテンプレートは、図2に示す画像記憶部31における第1画像記憶部32に記憶される。また、機械学習用の複数の正解画像は、図2に示す画像記憶部31における第2画像記憶部33に記憶される。また、画像記憶部31における第2画像記憶部33には、機械学習に使用される不正解画像も記憶される。
 図4は、テンプレートと正解画像の作成工程を示す説明図である。
 なお、この図においては、例えば、腫瘍等の被検者の特定部位Cを含む画像100、101、102を3枚連続で撮影した状態を示している。実際には、これらの画像は、所定のフレームレートでより多くの枚数が撮影される。オペレータは、これらの画像100、101、102から特定部位C部分をトリミングする。このようにオペレータがトリミングを行うかわりに、治療計画に用いる4DCT(経時的に撮影された3次元のCT画像データ群からなる4次元CT画像データ)における特定部位の位置から各位相の特定部位Cの投影位置を取得し、画像100、101、102と4DCTの対応関係から、特定部位Cの位置を認識して自動的にトリミングを実行するようにしてもよい。また、4DCTを利用して特定部位Cのおおよその位置を取得し、オペレータがこれを修正するようにしてもよい。
 図4に示すトリミング後の画像105、106、107の各々が、マルチテンプレートマッチングに使用されるテンプレートとなり、第1画像記憶部32に登録される。また、図4に示すトリミング後の画像105、106、107、および/または、これらのトリミング後の画像105、106、107を自動的にわずかに平行移動、回転(例えば、約10度)、変形、拡大縮小した画像が機械学習のための正解画像として生成され、第2画像記憶部33に登録される。トリミングした画像105、106、107を平行移動、回転、変形、拡大縮小した画像も正解画像として使用するのは、腫瘍等の被検者の特定部位Cが被検者の体内で移動または変形した際にも特定部位Cをより確実に追跡できるようにするためである。さらに、画像100のうち、トリミング後の画像105、106、107以外の領域の任意の領域の画像が自動的に抽出されて、不正解画像として生成され、第2画像記憶部33に登録される。なお、不正解画像をその都度作成する代わりに、予め記憶されている不正解画像を使用してもよい。
 上記の通り、複数のテンプレートと複数の正解画像とは、各々、同一の画像に基づいて作成された画像が使用される。すなわち、各テンプレートと各正解画像とは、特定部位Cを含む同一の画像がからトリミングされたものが使用される。これにより、複数のテンプレートと複数の正解画像とを、単一の工程により効率的に作成することが可能となる。
 次に、機械学習における学習工程を実行する(ステップS2)。この学習工程においては、第2画像記憶部33に記憶された複数の正解画像を利用して学習を実行することにより、機械学習のための一つの識別器を作成する。識別器は、計算式とパラメータとから構成されるものであり、機械学習を使用した特定部位Cの位置の追跡に使用される指標である。識別器は、第2画像記憶部33に記憶されている複数の正解画像および複数の不正解画像を使用して機械学習により作成される。計算式は機械学習の種別(後述)によって決まっているが、計算式自体は周知技術であるので、ここでは説明を省略する。この学習工程には、一定の時間が必要である。しかしながら、予め被検者における特定部位Cを含む画像の撮影が完了した後、実際に放射線治療を実行するまでの間の時間に、この学習工程を実行することにより、被検者に負担をかける必要はない。なお、特定部位Cを含む画像の撮影は、後述する放射線治療の直前に実行してもよく、それに先立ち、事前に実行してもよい。
 なお、上述した機械学習の種別としては、例えば、SVM(Support Vector Machine/サポートベクターマシン)を利用することができる。このSVMは、パターン認識を実行するときに、多くの手法の中でも最も迅速性に優れ、かつ、認識性能の高い学習モデルの一つである。また、迅速性に優れた機械学習として、SVMにかえて、Haar‐like特徴量などによるBoosting(ブースティング)や、Deep Learning(深層学習)などのニューラルネットワークを利用してもよい。
 以上の準備工程が完了すれば、被検者を、再度、検診台29上に載置し、この発明に係るX線透視装置によりX線透視による動体追跡を行うとともに、放射線照射装置90から放射線を照射して放射線治療を開始する。
 放射線治療を開始するときには、最初に、X線透視を行う(ステップS3)。このX線透視は、例えば30fps(frame per second)程度の、所定のフレームレートで実行される。そして、所定のフレームレートで取得された複数のX線画像を使用して、予め機械学習により作成された識別器を利用した識別により特定部位Cの1または複数の候補位置を選択する(ステップS4)。
 図5は、被検者Mに対するX線透視領域E1と、特定部位Cの候補位置E2とを模式的に示す説明図である。
 候補位置選択工程(ステップS4)においては、機械学習に基づく識別により、被検者Mに対するX線透視領域E1のうちから、特定部位Cの1または複数の候補位置E2を識別する第1判定を実行する。具体的には、X線透視領域E1の画像の全域を順番に識別器と比較を行い、識別器との類似度に応じて○(正解)か×(不正解)かを判定し、正解の位置を特定部位Cの候補位置と識別する。一般的には、1つのX線透視領域E1の中で特定部位Cの候補位置は2~3個が識別される。なお、識別器との比較において○(正解)か×(不正解)かの2値ではなく、スコア値で結果が出力される場合には、スコア値の閾値を設定しておき、閾値以上のスコア値の位置を特定部位の候補位置とすることができる。
 次に、候補位置選択工程で選択された特定部位Cの1または複数の候補位置E2に対して、画像記憶工程(ステップS1)で記憶された複数のテンプレートを利用してテンプレートマッチングを実行し、特定部位Cについての第2判定を実行することにより、特定部位Cの1または複数の候補位置の中から特定部位Cの位置を一つ検出する(ステップS5)。マルチテンプレートマッチングでは、例えば、X線透視領域E1のうちの特定部位Cの候補位置E2の画像と、第1画像記憶部32に記憶されている複数のテンプレートとを比較して、その類似を求め、類似度が所定の閾値以上の位置が特定部位Cの位置と判定される。この特定部位の位置の検出は、透視のフレームレートに一致するタイミングで実行される。このように、X線透視領域E1の全域ではなく、位置候補選択工程で選択された1または複数の候補位置E2に対してのみ、複数のテンプレートを利用したマルチテンプレートマッチングを実行すればよいことから、従来のように、計算コストが過度に大きくなり、リアルタイムでの処理が困難となるという問題が生ずることを防止することが可能となる。
 そして、特定部位Cの位置が検出できたら、その位置情報に基づいて、放射線照射装置90から被検者Mの患部(腫瘍である特定部位C)に対して放射線を照射する。すなわち、患部の位置が所定の範囲内にあるときに、放射線照射装置90のゲーティングをONとして被検者に放射線を照射する。この特定部位Cの位置の検出と放射線の照射とは、治療が終了するまで繰り返される(ステップS6)。
 以上のように、上述したX線透視装置によれば、機械学習に基づく識別を利用して特定部位Cの候補位置E2を判定し、次に、候補位置E2の中からマルチテンプレートマッチングを利用して特定部位の位置を検出することから、マーカのかわりに被検者Mの特定部位Cの画像を使用して複数のテンプレートを利用したテンプレートマッチングを実行するときにも、より容易に計算を実行してリアルタイムで特定部位Cの位置を検出することが可能となる。
 次に、この発明の他の実施形態について説明する。図6は、この発明の第2実施形態に係るX線透視装置を使用した動体追跡動作を示すフローチャートである。
 上述した実施形態においては、機械学習を利用して特定部位Cの候補位置E2を判定し、次に、候補位置E2の中からマルチテンプレートマッチングを利用して特定部位Cの位置を検出する構成を採用している。これに対して、この第2実施形態においては、テンプレートマッチングと機械学習に基づく識別との両方を利用して、特定部位Cの位置をより正確に検出するようにしている。
 すなわち、テンプレートマッチングは、探索範囲内の最も類似度が高いところを特定部位Cとして検出する構成であることから、特定部位Cが画像の範囲外に配置されていたとしても、類似度が高い部分を特定部位Cとして検出してしまうことがある。これに対して、機械学習では特定部位C以外を検出する可能性は低い。このため、マルチテンプレートマッチングと機械学習との結果を合わせ、両者とも同じ部分を特定部位Cと判定した場合にのみ、それが正しい特定部位Cの位置であると最終的に判定することにより、検出精度を向上させることが可能となる。
 この発明の第2実施形態に係るX線透視装置により動体追跡を実行するときには、最初に、画像記憶工程を実行する(ステップS11)。この画像記憶工程において、テンプレートマッチング用の複数のテンプレートは、図2に示す画像記憶部31における第1画像記憶部32に記憶される。また、機械学習用の複数の正解画像は、図2に示す画像記憶部31における第2画像記憶部33に記憶される。また、画像記憶部31における第2画像記憶部33には、機械学習に使用される不正解画像も記憶される。
 次に、機械学習における学習工程を実行する(ステップS12)。この学習工程においては、第2画像記憶部33に記憶された複数の正解画像を利用して学習を実行することにより、機械学習のための識別器を作成する。
 以上の準備工程が完了すれば、X線透視装置によりX線透視を実行する(ステップS13)ことにより動体追跡を行うとともに、放射線照射装置90から放射線を照射して放射線治療を開始する。
 このときには、画像記憶工程(ステップS11)で記憶された複数のテンプレートを利用してテンプレートマッチングを実行し(ステップS14)、特定部位Cの位置を検出する(ステップS15)。また、これと並行して、予め機械学習によって作成された識別器を利用した識別を実行することにより(ステップS16)により、特定部位Cの位置を検出する(ステップS17)。なお、識別器を利用した識別の手法は先の実施形態と同様である。
 そして、マルチテンプレートマッチングにより検出された特定部位Cの位置と機械学習に基づく識別により検出された特定部位Cの位置とが一致するか否かを判断する(ステップS18)。マルチテンプレートマッチングにより検出された特定部位Cの位置と機械学習に基づく識別により検出された特定部位Cの位置とが一致した場合には、その位置が特定部位Cの位置であると判定してゲーティングをONとする(ステップS19)。一方、マルチテンプレートマッチングにより検出された特定部位Cの位置と機械学習に基づく識別により検出された特定部位Cの位置とが一致しない場合には、特定部位Cの位置の判定は行わずゲーティングをOFFとする(ステップS20)。この特定部位Cの位置の検出と放射線の照射とは、治療が終了するまで繰り返される(ステップS21)。
 なお、図6に示す実施形態においては、テンプレートマッチング(ステップS14、ステップS15)と機械学習による位置検出(ステップS16、ステップS17)とを並行して実行している。しかしながら、テンプレートマッチング(ステップS14、ステップS15)を実行した後に、機械学習による位置検出(ステップS16、ステップS17)を実行するようにしてもよい。また、上述した第1実施形態のように、機械学習による位置検出(ステップS16、ステップS17)を実行した後に、テンプレートマッチング(ステップS14、ステップS15)を実行するようにしてもよい。また、マルチテンプレートマッチングと機械学習との双方で、特定部位Cとの類似度に相当するスコア値を算出し、マルチテンプレートマッチングにより検出されたスコア値と機械学習により検出されたスコア値との和または積が所定の閾値以上のときに特定部位Cの位置と特定するようにしてもよい。上記の実施形態では複数のテンプレートを用いたマルチテンプレートマッチングを行っているが、1つのテンプレートを用いてテンプレートマッチングを行ってもよい。
 次に、この発明のさらに他の実施形態について説明する。図7は、この発明の第3実施形態に係るX線透視装置の主要な制御系を示すブロック図である。
 このX線透視装置は、論理演算を実行するプロセッサーとしてのCPU、装置の制御に必要な動作プログラムが格納されたROM、制御時にデータ等が一時的にストアされるRAM等を備え、装置全体を制御する制御部30を備える。この制御部30は、上述した第1X線管11aおよび第2X線管11bと、第1フラットパネルディテクタ21aおよび第2フラットパネルディテクタ21bとに接続されている。そして、この制御部30は、後述するDRR画像作成部41と、識別器作成部42と、位置検出部43とを備える。
 また、この制御部30は、上述した放射線照射装置90と、治療計画装置99とに接続されている。なお、制御部30と治療計画装置99とは、病院内の被検者管理システムの院内通信である放射線科情報システム(RIS)を介して接続されてもよい。ここで、治療計画装置99は、放射線治療を行うに先だって、治療計画を作成するためのものである。この治療計画装置においては、CT撮影装置により被検者の3次元CT撮影を連続して複数回行うことにより得た、被検者の連続する複数の呼吸位相における特定部位を含む領域の3次元のCT画像データ群からなる4次元CT画像データを記憶している。そして、この4次元CT画像データと被検者のその他のデータとに基づいて、被検者の治療計画が作成される。
 次に、以上のような構成を有するX線透視装置を使用することにより、被検者の体動に伴って移動する特定部位の位置を検出する動体追跡を実行する動作について説明する。図8は、この発明の第3実施形態に係るX線透視装置を使用した動体追跡動作を示すフローチャートである。
 X線透視を実行するときには、最初に、図7に示すDRR画像作成部41により、治療計画記憶時に作成された4次元CT画像データに基づいて、特定部位を含む複数のDRR画像を作成する(ステップS31)。このDRR画像の作成時には、被検者に対するX線管11a、11bとフラットパネルディテクタ21a、21bとの幾何学的透視条件を模擬した仮想的透視投影を行う。ここで、治療計画時に作成される4次元CTデータとは、治療計画記憶時において、連続する複数の呼吸位相において、経時的に連続して撮影される特定部位を含む領域の3次元のCT画像データ群である。
 このDRR画像作成工程においては、治療計画時に作成される4次元CTデータの全てのCT画像データに基づいて、あるいは、治療計画時に作成される4次元CTデータのうち、少なくとも被検者に放射線照射装置90から治療ビームが照射される呼吸位相を含む複数の呼吸位相のCT画像データに基づいて、特定部位を含む複数のDRR画像が作成される。
 そして、DRR画像作成部41により作成された複数のDRR画像を使用して、図7に示す識別器作成部42により、機械学習を利用して特定部位を認識するための識別器を作成する(ステップS32)。このときには、4次元CT画像データにおいて、投影座標や角度などのDRR画像作成のためのパラメータを変化させて大量の正解画像を作成する。このとき、必要に応じ、治療計画で登録された4次元CT画像データにおける特定部位の位置と大きさから、DRR画像上の特定部位の位置と大きさを認識して自動的にトリミングを実行するようにしてもよい。
 この正解画像作成のためのDRR画像の作成時には、4次元CT画像データに対して投影座標及び角度の少なくとも一方を含むDRR画像作成のためのパラメータを変化させてDRR画像を作成する。あるいは、わずかな平行移動、回転、変形及び拡大縮小の少なくとも1つを含む画像処理を実行する。この平行移動、回転、変形、拡大縮小を実行するのは、腫瘍等の被検者の特定部位が被検者の体内で4次元CT画像データに対して再現性のない移動し、または、変形が生じたときにも特定部位をより確実に追跡できるようにするためである。
 ここで、DRR画像のもとになる4次元CT画像データのフレームレートはX線透視画像のフレームレートより小さいが、DRR画像作成のためのパラメータを変化させることにより、DRR画像におけるフレーム間の特定部位を模擬することが可能となる。また、X線透視画像は透過像であるため、特定部位の領域に骨部や横隔膜などの特定部位より高コントラストな背景が重畳して写り込むことになる。これに対して、上述したようにパラメータを変化させると、特定部位に対する背景の重なり方が様々に変化するので、それを学習することで、背景の重なりの影響を防止することが可能となる。
 また、作成されたDRR画像にコントラスト変化、ノイズ付加及びエッジ強調の少なくとも1つを実行する。このコントラスト変化、ノイズ付加、エッジ強調を実行するのは、DRR画像とX線画像の画質の違いを吸収し、特定部位をより確実に追跡できるようにするためである。
 上述した、投影座標や角度などのDRR画像作成のためのパラメータの変化、または、コントラスト変化、ノイズ付加、エッジ強調は、所定の範囲内でランダムに、あるいは、等間隔で様々に変化を与える態様で実施される。これにより、患者一人の4次元CT画像データから、多量のDRR画像を作成することができる。このため、このようにして作成された多量のDRR画像を使用して、各患者に対応したオーダーメイドの識別器を学習することが可能となる。なお、多数の患者のDRR画像を利用して識別器の学習を行うことも可能である。
 正解画像とともに使用される不正解画像は、例えば、以下の方法により作成される。すなわち、不正解画像を作成するときには、DRR画像作成部41により作成された特定部位を含むDRR画像から、特定部位を避けたランダムな位置、すなわち、特定画像の背景となる位置において、複数回のトリミングを行うことで、不正解画像が作成される。また、不正解画像を作成するときには、さらに、特定部位を含まないDRR画像を使用して不正解画像が作成される。
 上述した機械学習作業は、放射線治療に先だって、予め実行しておくことが可能となる。このため、放射線治療の直前にテンプレートを作成する場合のように、被検者を時間的に拘束することなく、また、放射線治療のスループットを向上させることが可能となる。
 この識別器作成工程で利用される機械学習としては、例えば、SVM (Support  Vector Machine/サポートベクターマシン)を利用することができる。このSVMは、パターン認識を実行するときに、多くの手法の中でも最も迅速性に優れ、かつ、認識性能の高い学習モデルの一つである。また、認識性能に優れた機械学習として、SVMにかえて、Haar‐like特徴量などによるBoosting(ブースティング)や、Deep Learning(深層学習)などのニューラルネットワークを利用してもよい。
 以上のような工程により識別器が作成されれば、上述した第1X線撮影系および第2X線撮影系を使用してX線透視を実行する(ステップS33)。そして、図7に示す位置検出部43により、先に作成された識別器を使用して識別を行うことにより、特定部位の位置を検出する(ステップS34)。このX線透視と特定部位の位置の検出を連続して実行することにより、特定部位の位置を追跡する動体追跡が実行される。
 そして、この状態において、放射線照射装置90により治療ビームを照射する(ステップS35)。すなわち、特定部位の位置が特定の位置となったときにのみ、治療ビームを照射して、放射線治療を行う。
 この第3実施形態に係るX線透視装置によれば、機械学習作業を放射線治療に先だって予め実行しておくことができることから、被検者の位置決め後に直ちに放射線治療を実行することが可能となる。
 また、この第3実施形態に係るX線透視装置によれば、X線画像ではなくDRR画像を使用して機械学習を実行することから、機械学習自体を効率的かつ迅速に実行することが可能となる。すなわち、機械学習においては、一般的に、10000枚以上の学習画像がないと実用的な識別精度とならない。しかしながら、X線画像を取得するためには治療直前にX線透視が必要となり時間を要するばかりではなく、正解画像を作成するときにオペレータによる手動トリミングが必要となって多量の正解画像を用意することは困難となる。この第3実施形態によれば、DRR画像を利用することにより、このような問題を解消することが可能となる。
 また、この第3実施形態に係るX線透視装置によれば、治療計画時に作成される4次元CTデータの全てのCT画像データに基づいて、あるいは、治療計画時に作成される4次元CTデータのうち、少なくとも被検者に放射線照射装置90から治療ビームが照射される呼吸位相を含む複数の呼吸位相のCT画像データに基づいて、特定部位を含む複数のDRR画像が作成される。このため、これらのCT画像に共通して存在する特定部位に対して学習ができるとともに、共通して存在しない背景を学習させないようにすることができる。従って、その学習結果に基づいて識別を行うことにより、被検者の特定部位とは異なった動きをする被検者の骨部等が認識されることを防止することが可能となる。また、これらのCT画像データに共通して存在する被検者の特定部位と同じ動きをする背景を学習することができることから、例えば、特定部位のコントラストが極端に小さい場合であっても、特定部位に替えて特定部位と同じ動きをする背景を追跡することが可能となる。このため、この背景を追跡することで、正確な位置に治療ビームを照射して放射線治療を実行することが可能となる。
 次に、この発明のさらに他の実施形態について説明する。図9は、この発明の第4実施形態に係るX線透視装置によるDRR作成工程を示すフローチャートである。
 この第4実施形態は、上述した第3実施形態におけるDRR画像作成工程(ステップS31)のみが異なっている。すなわち、この第4実施形態においては、DRR画像作成部41は、治療計画記憶時に作成されたCT画像データに基づいて、被検者の特定部位を含み、かつ、被検者の骨部が除去されたDRR画像を作成する。
 この第4実施形態におけるDRR画像作成工程(ステップS31)においては、最初に、CT画像データが読み出される(ステップS311)。このCT画像データは、上述した第3実施形態同様、4次元CTデータであってもよく、また、3次元CTデータであってもよい。
 次に、読み出されたCT画像データにおけるCT値が補正される(ステップS312)。すなわち、CT画像においては、被検者の骨部に相当する領域のCT値は1000から1500程度であり、軟部組織に相当する領域のCT値は0から200程度である。このため、CT画像データにおけるCT値が1000以上の領域は骨部領域であると認識し、その領域のCT値に対してその周辺の軟部組織のCT値を外挿することにより、その領域のCT値を軟部組織のCT値に変換する。なお、外挿に替えて、その領域のCT値を代表的な軟部組織のCT値である200に変換してもよい。そして、CT値を補正した後のCT画像データを保存する(ステップS313)。
 しかる後、DRR画像を作成する(ステップS314)。このときには、被検者の骨部領域のCT値がCT画像データから除外されていることから、被検者の骨部の情報は、作成されたDRR画像から除外されていることになる。このため、後段の識別器作成工程(ステップS32)において、被検者の特定部位に着目した機械学習により、より正確な識別器を作成することが可能となる。
 次に、この発明のさらに他の実施形態について説明する。図10は、この発明の第5実施形態に係るX線透視装置によるDRR作成工程を示すフローチャートである。
 この第5実施形態は、上述した第4実施形態と同様、第3実施形態におけるDRR画像作成工程(ステップS31)のみが異なっている。すなわち、この第5実施形態においても、DRR画像作成部41は、治療計画記憶時に作成されたCT画像データに基づいて、被検者の特定部位を含み、かつ、被検者の骨部が除去されたDRR画像を作成する。
 この第5実施形態におけるDRR画像作成工程(ステップS31)においては、最初に、CT画像データが読み出される(ステップS315)。このCT画像データは、上述した第3実施形態同様、4次元CTデータであってもよく、また、3次元CTデータであってもよい。
 次に、読み出されたCT画像データに対して仮想的な投影が実行され、CT画像データにおけるボクセル値が線積分される(ステップS316)。このボクセル値の線積分時に、CT画像データにおけるCT値が1000以上の領域は骨部領域であると認識し、その領域のCT値に対してその周辺の軟部組織のCT値を外挿することにより、その領域のCT値を軟部組織のCT値に変換して線積分する。なお、外挿に替えて、その領域のCT値を代表的な軟部組織のCT値である200に変換して線積分してもよい。
 そして、線積分の結果から、DRR画像を作成する(ステップS317)。このときには、上述した第4実施形態の場合と同様、被検者の骨部の情報は、作成されたDRR画像から除外されていることになる。このため、後段の識別器作成工程(ステップS32)において、被検者の特定部位に着目した機械学習により、より正確な識別器を作成することが可能となる。
 次に、この発明のさらに他の実施形態について説明する。図11は、この発明の第6実施形態に係るX線透視装置の主要な制御系を示すブロック図である。また、図12は、この発明の第6実施形態に係るX線透視装置を使用した動体追跡動作を示すフローチャートである。なお、上述した第3実施形態と同様の部材については、同一の符号を付して詳細な説明を省略する。
 この第6実施形態に係るX線透視装置における制御部30は、図11に示すように、DRR画像作成部41と、X線画像作成部44と、識別器作成部42と、補正部45と、位置検出部43とを備える。そして、この第6実施形態に係るX線透視装置においては、補正部45が、X線画像作成部44により作成されたX線画像を利用して、識別器作成部42において作成した識別器を補正する構成を採用している。
 この第6実施形態に係るX線透視装置においては、上述した第3実施形態の場合と同様に、DRR画像作成部41によりDRR画像を作成した後に(ステップS41)、識別器作成部42により機械学習を利用して識別器を作成する(ステップS42)。この作業は、上述した第3実施形態と同様に、放射線治療に先だって予め実行される。
 次に、図11に示すX線画像作成部44により、被検者の特定部位を含む画像を収集することにより、特定部位を含む複数のX線画像を作成する(ステップS43)。このときには、上述した第1X線撮影系および第2X線撮影系が使用される。X線撮影系により撮影されたX線画像は、必要に応じ、特定部位の領域のみトリミングがなされる。
 そして、DRR画像の場合と同様の手法により、特定部位のX線画像から機械学習に使用される正解画像が作成される。また、DRR画像の場合と同様、特定部位を含むX線画像のうち特定部位を避けたランダムな位置、すなわち、特定画像の背景となる位置において複数回のトリミングを行うことや、特定部位を含まない画像を使用することで不正解画像が作成される。
 しかる後、図11に示す補正部45により、先に作成された識別器が補正される(ステップS44)。これにより、より精度の高い識別器を作成することが可能となる。このときには、機械学習としてSVMを利用した場合には、X線画像で機械学習を行って新たな識別器を作成し、先に作成した識別器と併用する形で新たな識別器を追加する。または、先に作成した識別器と新たな識別器の分離超平面パラメータを合成し、識別器自体を補正する。また、機械学習としてBoostingを利用した場合には、X線画像で機械学習を行って新たな弱識別器を作成し、先に作成した弱識別器の集合である識別器に新たな弱識別器を追加する。また、機械学習としてDeep Learningなどのニューラルネットワークを利用した場合には、先に作成した識別器のパラメータを初期値として、さらにX線画像で機械学習を行い、識別器自体を補正する。この明細書における識別器の補正とは、上述した識別器への新たな識別器の追加と、識別器自体の補正の両方を含む概念である。
 以上のような工程により識別器が補正されれば、X線透視を実行する(ステップS45)。そして、図11に示す位置検出部43により識別を行うことにより、特定部位の位置を検出する(ステップS46)。このX線透視と特定部位の位置の検出を連続して実行することにより、特定部位の位置を追跡する動体追跡が実行される。そして、この状態において、放射線照射装置90により治療ビームを照射する(ステップS47)。
 この第6実施形態に係るX線透視装置によれば、X線画像を利用して識別器作成部42において作成した識別器を補正することから、より正確に特定部位の位置を検出することが可能となる。このとき、予めDRR画像に基づく機械学習を実行していることから、多量のX線画像を使用することなく識別器を補正することが可能となる。このため、上述したような煩雑なオペレータによる手動トリミング作業等を減少させることが可能となり、識別器の補正のための手間と時間を小さなものとすることが可能となる。
 次に、この発明のさらに他の実施形態について説明する。図13は、この発明の第7実施形態に係るX線透視装置の主要な制御系を示すブロック図である。また、図14は、この発明の第7実施形態に係るX線透視装置を使用した動体追跡動作を示すフローチャートである。なお、上述した第3実施形態および第6実施形態と同様の部材については、同一の符号を付して詳細な説明を省略する。
 この第7実施形態に係るX線透視装置における制御部30は、図13に示すように、DRR画像作成部41と、X線画像作成部44と、識別器作成部42と、位置検出部43とを備える。そして、この第6実施形態に係るX線透視装置においては、識別器作成部42が、DRR画像作成部41で作成されたDRR画像とX線画像作成部44で作成されたX線画像とに対して機械学習を実行することにより、特定部位を認識するための識別器を作成する構成を採用している。
 この第7実施形態に係るX線透視装置においては、上述した第3実施形態および第6実施形態の場合と同様に、DRR画像作成部41によりDRR画像を作成する(ステップS51)。また、図13に示すX線画像作成部44により、被検者の特定部位を含む画像を収集することにより、特定部位を含む複数のX線画像を作成する(ステップS52)。
 そして、図13に示す識別器作成部42により機械学習を利用して識別器を作成する(ステップS53)。このときには、第1、第6実施形態の場合と同様に、DRR画像に基づいて作成された正解画像と不正解画像とを利用するとともに、X線画像に基づいて作成された正解画像と不正解画像とを利用して、DRR画像の機械学習とX線画像の機械学習を同時に実行することにより、識別器が作成される。
 このように、この第7実施形態においては、DRR画像作成部41で作成されたDRR画像とX線画像作成部44で作成されたX線画像とに対して機械学習を実行することにより、特定部位を認識するための識別器が作成される。このため、識別器の精度を高いものとすることが可能となる。また、上述したように、機械学習においては多量の学習画像が必要となるが、X線画像を利用した機械学習用のX線画像を取得するためには治療直前にX線透視が必要となり時間を要するばかりではなく、正解画像を作成するときにオペレータによる手動トリミングが必要となって多量の正解画像を用意することは困難となるという問題が生ずる。この第7実施形態によれば、DRR画像を利用することにより、必要となるX線画像を少量とすることができ、識別器の作成をより効率的に実行することが可能となる。 
 なお、上述した説明においては、DRR画像を作成した後にX線画像を作成しているが、この順序は逆であってもよい。また、DRR画像の作成とX線画像の作成とを同時に実行してもよい。
 以上のような工程により識別器が作成されれば、X線透視を実行する(ステップS54)。そして、図13に示す位置検出部43により識別を行うことにより、特定部位の位置を検出する(ステップS55)。このX線透視と特定部位の位置の検出を連続して実行することにより、特定部位の位置を追跡する動体追跡が実行される。そして、この状態において、放射線照射装置90により治療ビームを照射する(ステップS56)。
 この第7実施形態に係るX線透視装置によれば、4次元CT画像データに基づいて作成されたDRR画像とX線撮影により作成されたX線画像とに対して機械学習を実行することにより、より正確に特定部位の位置を検出することが可能となる。このとき、DRR画像に基づく機械学習とX線画像に基づく機械学習とを併用することから、多量のX線画像を使用することなく識別器を迅速に作成することが可能となる。
 次に、この発明のさらに他の実施形態について説明する。図15は、この発明の第8実施形態に係るX線透視装置の主要な制御系を示すブロック図である。なお、上述した第3実施形態と同様の構成については、同一の符号を付して詳細な説明を省略する。
 上述した第3実施形態から第7実施形態においては、機械学習を利用して特定部位を認識するための識別器として2クラス識別器を使用しているのに対し、この第8実施形態においては、特定部位を呼吸位相毎に識別するための多クラス識別器を使用している。すなわち、この第8実施形態においては、上述した第3実施形態における識別器作成部42にとして、特に、多クラス識別器作成部46を採用し、この多クラス識別器作成部46により、特定部位を呼吸位相毎に識別するための多クラスの識別器を作成している。また、この第8実施形態においては、上述した第3実施形態における位置検出部43として位置・位相検出部47を採用し、特定部位の位置と位相とを検出する構成を採用している。
 図16は、この発明の第8実施形態に係るX線透視装置を使用した動体追跡動作を示すフローチャートである。
 X線透視を実行するときには、上述した第3実施形態の場合と同様、最初に、図15に示すDRR画像作成部41により、治療計画記憶時に作成された4次元CT画像データに基づいて、特定部位を含む複数のDRR画像を作成する(ステップS61)。このDRR画像作成工程においては、治療計画時に作成される4次元CTデータの全てのCT画像データに基づいて、あるいは、治療計画時に作成される4次元CTデータのうち、少なくとも被検者に放射線照射装置90から治療ビームが照射される呼吸位相を含む複数の呼吸位相のCT画像データに基づいて、特定部位を含む複数のDRR画像が作成される。
 そして、DRR画像作成部41により作成された複数のDRR画像を使用して、図15に示す多クラス識別器作成部46により、機械学習を利用して、特定部位の位置と呼吸位相とを認識するための多クラス識別器を作成する(ステップS62)。このときには、4次元CT画像データにおいて、投影座標や角度などのDRR画像作成のためのパラメータを変化させて、呼吸位相毎に、大量の正解画像を作成する。
 このときには、上述した第3実施形態の場合と同様に作成した多数の正解画像と不正解画像とを利用して、機械学習により多クラス識別器を作成する。この機械学習作業は、放射線治療に先だって、予め実行しておくことが可能となる。このため、放射線治療の直前にテンプレートを作成する場合のように、被検者を時間的に拘束することなく、また、放射線治療のスループットを向上させることが可能となる。
 この多クラス識別器作成工程で利用される機械学習としては、例えば、畳み込みニューラルネットワーク(CNN/Convolutional Neural Network)に代表されるDeep Learning(深層学習)を利用することができる。この畳み込みニューラルネットワークは、パターン認識を実行するときに、多くの手法の中でも最も認識性能の高い学習モデルの一つである。
 図17は、畳み込みニューラルネットワークの概念を説明するための説明図である。
 畳み込みニューラルネットワークは、X線透視画像から構成される入力層INと,複数のノードから成る出力層OUTとの間に、中間層としての畳み込み層(convolution layer)CLと全結合層(fully connected layer)FLとが配置された構成を有する。この畳み込みニューラルネットワークは、ニューラルネットワークに相当する全結合層FLに対して畳み込み処理が行われる畳み込み層CLを追加した構成を有する。
 出力層OUTは、図17に示すように、被検者の呼吸位相に対応したn個のノードP1~Pnと、背景に対応したノードBkとから構成され、ソフトマックス関数により全ノードの値の合計値が1になるように規格化されている。すなわち、各ノードは入力画像が各クラスに属する確率に相当する値となる。畳み込み層CLと全結合層FLにおいては、学習によってパラメータが決定される。この学習時において、呼吸位相m(mは1以上n以下の整数)の正解画像を入力した時は、出力値の正解を呼吸位相mのノードPmnの値を1、それ以外のノードの値を0とする。また、不正解画像を入力した時は、背景のノードBkの値を1、それ以外のノードの値を0とする。これにより、機械学習が実行される。
 以上のような工程により識別器が作成されれば、上述した第1X線撮影系および第2X線撮影系を使用してX線透視を実行する(ステップS63)。そして、図15に示す位置・位相検出部47により、先に作成された識別器を使用して識別を行うことにより、特定部位の位置とそのときの被検者の呼吸位相とを検出する(ステップS64)。
 この特定部位の検出時には、X線透視画像をサーチした検出窓の画像を入力し、そのときの出力値(確率)の最も大きいノードを検出する。すなわち、背景のノードの値が最大であつた場合は、特定部位を検出していないと判断し、また、呼吸位相nのノードPnの値が最大であった場合は、特定部位を検出し、かつ、そのときの呼吸位相はnであると判断する。
 上述したX線透視と、特定部位の位置および被検者の呼吸位相の検出とを連続して実行することにより、特定部位の位置とそのときの被検者の呼吸位相とを追跡する動体追跡が実行される。そして、この状態において、放射線照射装置90により治療ビームを照射する(ステップS65)。すなわち、特定部位の位置が特定の位置となり、かつ、そのときの被検者の呼吸位相が治療計画時に計画された呼吸位相に一致したときにのみ、治療ビームを照射して、放射線治療を行う。
 なお、被検者の呼吸位相は連続して変化することから、この呼吸位相は、厳密には、出力値が最も大きいノードと2番目に大きいノードとの間の呼吸位相となる。したがって、両者のノード値に基づいて重みづけを行い、呼吸位相を補間してもよい。予めDRR画像上の特定部位の位置を呼吸位相で補間して求めておけば、補間した呼吸位相から予想される特定部位の位置が取得できる。
 図18は、このような補間に使用可能な周期関数を示すグラフである。なお、このグラフは、被検者の呼吸の1サイクルの間に、一定時間毎に10フレームのDRR画像が作成された場合を想定したものである。このグラフにおいては、縦軸は被検者の呼吸位相が最大呼気相となったときの特定部位の位置を基準としたDRR画像上での特定部位のY方向の位置(cm)を示し、横軸は呼吸位相(4次元CT画像データに基づいて作成されたDRR画像のフレーム番号)を示している。ここで、Y方向とは、特定部位の主要な移動方向である。特定部位は、DRR画像上でX、Y方向に移動するが、図18においては、各DRR画像におけるY方向の位置を示している。X方向の位置についても、これと同様である。
 被検者の呼吸位相に伴う特定部位の位置は、図18に示すような被検者の呼吸位相を模擬した周期関数の回帰により、高精度で補間することが可能となる。
 なお、呼吸位相から予想される特定部位の位置と検出された特定部位の位置を比較し、それらの間に大幅に乖離がある場合は誤検出と判断するなど、これを治療ビームの誤照射防止に使用することも可能である。また、両者の特定部位の位置を合成して、最終的な特定部位の位置を算出してもよい。
 また、第一の識別器を、上述した第3実施形態と同様の特定部位と背景の2クラス分類とした機械学習とし、特定部位が検出された場合は、第二の識別として畳み込みニューラルネットワークによる多クラス分類を利用して呼吸位相を検出してもよい。一般に畳み込みニューラルネットワークは計算コストが大きいので、計算コストの小さい機械学習で特定部位の探索を行い、これにより検出された特定部位に対してのみ、またはその周辺のみ、畳み込みニューラルネットワークを行うことで、高速に特定部位の位置と呼吸位相とを検出することが可能となる。また、2クラス分類を複数回繰り返すことにより、多クラス分類を実施する構成としてもよい。
 図19および図20は、畳み込みニューラルネットワークを利用して特定部位の位置を認定する動作を示す模式図である。なお、図19における符号Sは探索窓を示し、符号CはX線透視画像における特定部位の領域を示している。
 図19に示すように、X線透視画像100に対して検出窓Sを走査させることにより、特定部位を探索する。そして、図20に示すように、特定部位である確率、すなわち、全ての呼吸位相のノードの値の和のマップPを作成する。このマップPは、特定部位が存在する画素周辺で確率が大きくなるので、特定部位である確率分布(確率密度関数)に相当する。この確率密度関数が最大値となる位置が、特定部位を検出した位置となる。確率密度関数が最大値となる位置は、ミーンシフト法などのカーネル密度推定によりサブピクセルで高精度に算出することができる。また、確率を重みとしたマップPの重心位置を、特定部位Cを検出した位置としてもよい。
 この第8実施形態に係るX線透視装置によれば、測定部位の位置とともに、そのときの被検者の呼吸位相をも特定することができることから、特定部位が所定の位置にあり、かつ、被検者の呼吸位相が特定の位相となったときに、治療ビームを照射することが可能となる。従って、治療計画に即した放射線照射を実行することが可能となる。
 11a   第1X線管
 11b   第2X線管
 21a   第1フラットパネルディテクタ
 21b   第2フラットパネルディテクタ
 29    検診台
 30    制御部
 31    画像記憶部
 32    第1画像記憶部
 33    第2画像記憶部
 34    学習部
 35    候補位置選択部
 36    位置検出部
 41    DRR画像作成部
 42    識別器作成部
 43    位置検出部
 44    X線画像作成部
 45    補正部
 46    多クラス識別器作成部
 47    位置・位相検出部
 90    放射線照射装置
 99    治療計画装置
 E1    X線透視領域
 E2    選択領域
 C     特定部位
 M     被検者

Claims (19)

  1.  X線管と、前記X線管から照射され被検者を通過したX線を検出するX線検出器とを備え、前記被検者の特定部位を含む画像を収集することにより、前記特定部位の位置を検出し、前記特定部位の動きを追跡するX線透視装置であって、
     前記被検者の特定部位を含む画像に基づいて作成されたテンプレートマッチング用のテンプレートと、前記被検者の特定部位を含む画像に基づいて作成された機械学習用の正解画像とを記憶する画像記憶部と、
     前記画像記憶部に記憶された正解画像に基づいて、前記機械学習により識別器を作成する学習部と、
     前記被検者の特定部位を含む領域を透視して得た画像に対して、前記テンプレートを利用してテンプレートマッチングを実行するとともに、前記識別器を利用して識別を実行することにより、前記特定部位の位置を検出する位置検出部と、
     を備えたことを特徴とするX線透視装置。
  2.  請求項1に記載のX線透視装置において、
     前記位置検出部は、前記テンプレートマッチングと前記識別器を利用した識別とのいずれか一方を実行した後に他方を実行し、あるいは、両方を同時に実行するX線透視装置。
  3.  X線管と、前記X線管から照射され被検者を通過したX線を検出するX線検出器とを備え、前記被検者の特定部位を含む画像を収集することにより、前記特定部位の位置を検出し、前記特定部位の動きを追跡するX線透視装置であって、
     前記被検者の特定部位を含む画像に基づいて作成されたテンプレートマッチング用のテンプレートと、機械学習用の正解画像とを記憶する画像記憶部と、
     前記画像記憶部に記憶された正解画像に基づいて、前記機械学習により識別器を作成する学習部と、
     前記X線管と前記X線検出器により収集された画像を使用して、前記識別器を利用した識別により特定部位候補位置を選択する候補位置選択部と、
     前記候補位置選択部により選択された特定部位を含む領域に対して、前記画像記憶部に記憶されたテンプレートを利用してテンプレートマッチングを実行することにより、前記特定部位の位置を検出する位置検出部と、
     を備えたことを特徴とするX線透視装置。
  4.  請求項1から請求項3のいずれかに記載のX線透視装置において、
     前記機械学習は、SVM、BoostingまたはニューラルネットワークであるX線透視装置。
  5.  請求項1から請求項3のいずれかに記載のX線透視装置において、
     前記テンプレートと前記正解画像とは、各々、前記特定部位を含む同一の画像に基づいて作成されるX線透視装置。
  6.  X線管から照射され被検者を通過したX線をX線検出器により検出し、前記被検者の特定部位を含む画像を収集することにより、前記特定部位の位置を検出し、前記特定部位の動きを追跡するX線透視方法であって、
     前記被検者の特定部位を含む画像に基づいて作成されたテンプレートマッチング用のテンプレートと、前記被検者の特定部位を含む画像に基づいて作成された機械学習用の正解画像とを記憶する記憶工程と、
     前記記憶工程で記憶された正解画像に基づいて、前記機械学習により識別器を作成する学習工程と、
     前記X線管から照射され前記被検者を通過したX線を前記X線検出器により検出し、前記被検者の特定部位を含む画像を収集する透視工程と、
     前記透視工程で得た画像に対して、前記テンプレートを利用したテンプレートマッチングを実行するとともに、前記識別器を利用した識別を実行することにより、前記特定部位の位置を検出する位置検出工程と、
     を含むことを特徴とするX線透視方法。
  7.  X線管から照射され被検者を通過したX線をX線検出器により検出し、前記被検者の特定部位を含む画像を収集することにより、前記特定部位の位置を検出し、前記特定部位の動きを追跡するX線透視方法であって、
     前記被検者の特定部位を含む画像に基づいて作成されたテンプレートマッチング用のテンプレートと、機械学習用の正解画像とを記憶する画像記憶工程と、
     前記画像記憶工程で記憶された複数の正解画像に基づいて、前記機械学習により識別器を作成する学習工程と、
     前記X線管から照射され前記被検者を通過したX線を前記X線検出器により検出し、前記被検者の特定部位を含む画像を収集する透視工程と、
     前記透視工程で得た画像を使用して、前記識別器を利用した識別により前記特定部位の候補位置を選択する候補位置選択工程と、
     前記候補位置選択工程で選択された特定部位の候補位置に対して、前記画像記憶工程で記憶されたテンプレートを利用してテンプレートマッチングを実行することにより、前記特定部位の位置を検出する位置検出工程と、
     を含むことを特徴とするX線透視方法。
  8.  X線管から照射され被検者を通過したX線をX線検出器により検出して被検者の特定部位を含む画像を取得し、前記被検者の特定部位を含む画像から前記特定部位の位置を検出し、前記特定部位の動きを追跡するX線透視装置であって、
     治療計画記憶時に作成されたCT画像データに対して、前記被検者に対する前記X線管と前記X線検出器との幾何学的透視条件を模擬した仮想的透視投影を行うことにより、前記特定部位を含む複数のDRR画像を作成するDRR画像作成部と、
     前記DRR画像作成部で作成されたDRR画像に対して機械学習を実行することにより、前記特定部位を認識するための識別器を作成する識別器作成部と、
     前記X線管から照射され前記被検者を通過したX線を前記X線検出器により検出して得た前記被検者の特定部位を含む画像に対して前記識別器作成部で作成された識別器を使用して識別を行うことにより、前記特定部位の位置を検出する位置検出部と、
     を備えることを特徴とするX線透視装置。
  9.  請求項8に記載のX線透視装置において、
     前記DRR画像作成部は、前記幾何学的透視条件の投影座標および角度の少なくとも一方を含むパラメータを変化させ、あるいは、画像の回転、変形及び拡大縮小の少なくとも1つを含む画像処理を施して前記DRR画像を作成するX線透視装置。
  10.  請求項8に記載のX線透視装置において、
     前記DRR画像作成部は、作成後のDRR画像に対して、コントラスト変化、ノイズ付加及びエッジ強調の少なくとも1つを実行するX線透視装置。
  11.  請求項8に記載のX線透視装置において、
     前記DRR画像作成部は、治療計画記憶時に作成された連続する複数の呼吸位相における前記特定部位を含む領域の3次元のCT画像データ群からなる4次元CT画像データに基づいて、前記特定部位を含む複数のDRR画像を作成するX線透視装置。
  12.  請求項11に記載のX線透視装置において、
     前記DRR画像作成部は、治療計画記憶時に作成された連続する複数の呼吸位相における前記特定部位を含む領域の3次元のCT画像データ群からなる4次元CT画像データのうち、前記被検者に治療ビームが照射される呼吸位相を含む複数の呼吸位相のCT画像データに基づいて、前記特定部位を含む複数のDRR画像を作成するX線透視装置。
  13.  請求項8に記載のX線透視装置において、
     前記DRR画像作成部は、治療計画記憶時に作成されたCT画像データに基づいて、前記特定部位を含み、かつ、前記被検者の骨部が除去されたDRR画像を作成するX線透視装置。
  14.  請求項13に記載のX線透視装置において、
     前記DRR画像作成部は、治療計画記憶時に作成されたCT画像データにおけるCT値に基づいて前記被検者の骨部が除去されたCT画像データを作成し、この骨部が除去されたCT画像データに基づいて、前記特定部位を含み、かつ、前記被検者の骨部が除去されたDRR画像を作成するX線透視装置。
  15.  請求項13に記載のX線透視装置において、
     前記DRR画像作成部は、治療計画記憶時に作成されたCT画像データにおけるボクセル値を線積分してDRR画像を作成するときに、このボクセル値に基づいて、前記特定部位を含み、かつ、前記被検者の骨部が除去されたDRR画像を作成するX線透視装置。
  16.  請求項8に記載のX線透視装置において、
     前記X線管から照射され前記被検者を通過したX線を前記X線検出器により検出し、前記被検者の特定部位を含む画像を収集することにより、前記特定部位を含む複数のX線画像を作成するX線画像作成部と、
     前記X線画像作成部により作成されたX線画像を利用して、前記識別器作成部において作成した識別器を補正する補正部と、
     をさらに備えるX線透視装置。
  17.  請求項8に記載のX線透視装置において、
     前記X線管から照射され前記被検者を通過したX線を前記X線検出器により検出し、前記被検者の特定部位を含む画像を収集することにより、前記特定部位を含む複数のX線画像を作成するX線画像作成部をさらに備え、
     前記識別器作成部は、前記DRR画像作成部で作成されたDRR画像と前記X線画像作成部で作成されたX線画像とに対して機械学習を実行することにより、前記特定部位を認識するための識別器を作成するX線透視装置。
  18.  請求項8に記載のX線透視装置において、
     前記識別器は、前記特定部位を前記被検者の呼吸位相毎に認識するための多クラス識別器であり、
     前記位置検出部は、前記X線管から照射され前記被検者を通過したX線を前記X線検出器により検出して得た前記被検者の特定部位を含む画像に対して前記識別器作成部で作成された多クラス識別器を使用して多クラス識別を行うことにより、前記特定部位の位置と前記被検者の呼吸位相とを検出するX線透視装置。
  19.  X線管から照射され被検者を通過したX線をX線検出器により検出して被検者の特定部位を含む画像を取得し、前記被検者の特定部位を含む画像から前記特定部位の位置を検出し、前記特定部位の動きを追跡するX線透視方法であって、
     治療計画記憶時に作成されたCT画像データに基づいて、前記特定部位を含む複数のDRR画像を作成するDRR画像作成工程と、
     前記DRR画像作成工程において作成されたDRR画像に対して機械学習を実行することにより、前記特定部位を認識するための識別器を作成する識別器作成工程と、
     前記X線管から照射され前記被検者を通過したX線を前記X線検出器により検出し、前記被検者の特定部位を含む画像を得る透視工程と、
     前記透視工程において得た画像に対して前記識別器作成工程で作成された識別器を使用して識別を行うことにより、前記特定部位の位置を検出する位置検出工程と、
     を含むことを特徴とするX線透視方法。
PCT/JP2017/026330 2016-11-02 2017-07-20 X線透視装置およびx線透視方法 WO2018083844A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US16/347,048 US11045663B2 (en) 2016-11-02 2017-07-20 X-ray fluoroscopy device and x-ray fluoroscopy method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2016215438 2016-11-02
JP2016-215438 2016-11-02

Publications (1)

Publication Number Publication Date
WO2018083844A1 true WO2018083844A1 (ja) 2018-05-11

Family

ID=62075942

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2017/026330 WO2018083844A1 (ja) 2016-11-02 2017-07-20 X線透視装置およびx線透視方法

Country Status (3)

Country Link
US (1) US11045663B2 (ja)
JP (1) JP7125700B2 (ja)
WO (1) WO2018083844A1 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020217569A1 (ja) * 2019-04-26 2020-10-29 株式会社日立製作所 患者位置決めシステム、方法およびプログラム
JP2021041090A (ja) * 2019-09-13 2021-03-18 株式会社島津製作所 医用画像処理装置、x線画像処理システム、および、学習モデルの生成方法
WO2023068365A1 (ja) * 2021-10-21 2023-04-27 国立大学法人筑波大学 標的外形推定装置および治療装置

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2536650A (en) 2015-03-24 2016-09-28 Augmedics Ltd Method and system for combining video-based and optic-based augmented reality in a near eye display
WO2019138438A1 (ja) * 2018-01-09 2019-07-18 株式会社島津製作所 画像作成装置
WO2019211741A1 (en) 2018-05-02 2019-11-07 Augmedics Ltd. Registration of a fiducial marker for an augmented reality system
KR102150478B1 (ko) * 2018-06-21 2020-09-01 한국외국어대학교 연구산학협력단 딥러닝 기반의 x-선 영상 뼈 감쇄 장치 및 방법
JP7097794B2 (ja) * 2018-10-18 2022-07-08 富士フイルム医療ソリューションズ株式会社 情報処理システム及び情報処理方法
US11766296B2 (en) 2018-11-26 2023-09-26 Augmedics Ltd. Tracking system for image-guided surgery
JP7311859B2 (ja) * 2019-03-25 2023-07-20 株式会社日立製作所 動体追跡装置とそれを備えた放射線治療システム、プログラム、および動体の追跡方法
JP7330833B2 (ja) * 2019-09-20 2023-08-22 株式会社日立製作所 放射線撮像装置および放射線治療装置
JP7292721B2 (ja) * 2019-10-08 2023-06-19 国立大学法人 筑波大学 標的外形推定装置および治療装置
US11676271B2 (en) 2019-11-01 2023-06-13 Konica Minolta, Inc. Dynamic image analysis apparatus extracting specific frames including a detection target from a dynamic image, dynamic analysis system, and storage medium
US11382712B2 (en) 2019-12-22 2022-07-12 Augmedics Ltd. Mirroring in image guided surgery
US11690581B2 (en) 2020-09-07 2023-07-04 C-Rad Positioning Ab Tumor position determination
US11896445B2 (en) 2021-07-07 2024-02-13 Augmedics Ltd. Iliac pin and adapter

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040116804A1 (en) * 1998-10-23 2004-06-17 Hassan Mostafavi Method and system for radiation application
JP2010187991A (ja) * 2009-02-19 2010-09-02 Hitachi Ltd ベッド位置決めシステム、放射線治療システム及びベッド位置決め方法
JP2016096902A (ja) * 2014-11-19 2016-05-30 株式会社東芝 医用画像処理装置、方法、プログラム及び放射線治療装置
JP2016116659A (ja) * 2014-12-19 2016-06-30 株式会社東芝 医用画像処理装置、治療システム、医用画像処理方法、および医用画像処理プログラム

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3053389B1 (ja) 1998-12-03 2000-06-19 三菱電機株式会社 動体追跡照射装置
JP4329398B2 (ja) 2002-05-10 2009-09-09 ソニー株式会社 顔検出装置及び方法、並びにプログラム及び記録媒体
JP2005078376A (ja) 2003-08-29 2005-03-24 Sony Corp 対象物検出装置、対象物方法、及びロボット装置
JP5159242B2 (ja) * 2007-10-18 2013-03-06 キヤノン株式会社 診断支援装置、診断支援装置の制御方法、およびそのプログラム
JP2010246883A (ja) * 2009-03-27 2010-11-04 Mitsubishi Electric Corp 患者位置決めシステム
JP5791373B2 (ja) * 2011-05-24 2015-10-07 キヤノン株式会社 特徴点位置決定装置、特徴点位置決定方法及びプログラム
JP6004464B2 (ja) * 2012-03-19 2016-10-05 国立大学法人北海道大学 放射線治療制御装置および放射線治療制御プログラム
JP6031258B2 (ja) 2012-05-23 2016-11-24 シャープ株式会社 表示装置
JP6301133B2 (ja) 2014-01-14 2018-03-28 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
WO2016046870A1 (ja) 2014-09-22 2016-03-31 株式会社島津製作所 X線透視装置、放射線治療用動体追跡装置及びx線検出器
US9675317B2 (en) * 2014-12-22 2017-06-13 Toshiba Medical Systems Corporation Interface identification apparatus and method
US11302435B2 (en) * 2016-01-06 2022-04-12 Boston Scientific Scimed, Inc. Systems and methods for planning medical procedures
WO2018076109A1 (en) * 2016-10-24 2018-05-03 Torus Biomedical Solutions Inc. Systems and methods for producing real-time calibrated stereo long radiographic views of a patient on a surgical table

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040116804A1 (en) * 1998-10-23 2004-06-17 Hassan Mostafavi Method and system for radiation application
JP2010187991A (ja) * 2009-02-19 2010-09-02 Hitachi Ltd ベッド位置決めシステム、放射線治療システム及びベッド位置決め方法
JP2016096902A (ja) * 2014-11-19 2016-05-30 株式会社東芝 医用画像処理装置、方法、プログラム及び放射線治療装置
JP2016116659A (ja) * 2014-12-19 2016-06-30 株式会社東芝 医用画像処理装置、治療システム、医用画像処理方法、および医用画像処理プログラム

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020217569A1 (ja) * 2019-04-26 2020-10-29 株式会社日立製作所 患者位置決めシステム、方法およびプログラム
JP2020178948A (ja) * 2019-04-26 2020-11-05 株式会社日立製作所 患者位置決めシステム、方法およびプログラム
JP7287828B2 (ja) 2019-04-26 2023-06-06 株式会社日立製作所 患者位置決めシステム、方法およびプログラム
US11712582B2 (en) 2019-04-26 2023-08-01 Hitachi, Ltd. Patient position determination system, method, and program
JP2021041090A (ja) * 2019-09-13 2021-03-18 株式会社島津製作所 医用画像処理装置、x線画像処理システム、および、学習モデルの生成方法
JP7226207B2 (ja) 2019-09-13 2023-02-21 株式会社島津製作所 医用画像処理装置、x線画像処理システム、および、学習モデルの生成方法
WO2023068365A1 (ja) * 2021-10-21 2023-04-27 国立大学法人筑波大学 標的外形推定装置および治療装置

Also Published As

Publication number Publication date
US20200069967A1 (en) 2020-03-05
US11045663B2 (en) 2021-06-29
JP2018075356A (ja) 2018-05-17
JP7125700B2 (ja) 2022-08-25

Similar Documents

Publication Publication Date Title
JP7125700B2 (ja) X線透視装置およびx線透視方法
JP7099459B2 (ja) 放射線治療用追跡装置、位置検出装置および動体追跡方法
KR102094737B1 (ko) 의료 영상 처리 장치, 의료 영상 처리 방법, 컴퓨터 판독가능한 의료 영상 처리 프로그램, 이동 대상물 추적 장치 및 방사선 치료 시스템
US10143431B2 (en) Medical image processing apparatus and method, and radiotherapeutic apparatus
US9684961B2 (en) Scan region determining apparatus
JP6785485B2 (ja) X線透視装置
CN111918697B (zh) 医用图像处理装置、治疗系统以及存储介质
US9259200B2 (en) Method and system for obtaining a sequence of x-ray images using a reduced dose of ionizing radiation
KR20210020958A (ko) 의용 장치 및 의용 장치의 제어 방법
JP2019037350A (ja) 放射線治療用追跡装置
JP6095112B2 (ja) 放射線治療システム
KR102469141B1 (ko) 의료용 화상 처리 장치, 의료용 화상 처리 방법, 및 프로그램
CN114067994A (zh) 一种目标部位的方位标记方法及系统
JP6792840B2 (ja) X線透視装置及びx線透視方法
JP2021061895A (ja) 医用画像処理装置、医用画像処理プログラム、医用装置、および治療システム
WO2022190846A1 (ja) 医用画像処理装置、医用画像処理方法、医用画像処理プログラム、および放射線治療装置
JP7125703B2 (ja) 医用装置、医用装置の制御方法およびプログラム
CN117726672A (zh) 一种介入设备实时三维成像方法、系统和装置
JP2019170884A (ja) X線透視装置
JP2018079011A (ja) X線透視方法およびx線透視装置
JP2018068584A (ja) X線透視装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17867790

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 17867790

Country of ref document: EP

Kind code of ref document: A1