JP2023034671A - X線診断装置及び医用情報処理装置 - Google Patents
X線診断装置及び医用情報処理装置 Download PDFInfo
- Publication number
- JP2023034671A JP2023034671A JP2021141010A JP2021141010A JP2023034671A JP 2023034671 A JP2023034671 A JP 2023034671A JP 2021141010 A JP2021141010 A JP 2021141010A JP 2021141010 A JP2021141010 A JP 2021141010A JP 2023034671 A JP2023034671 A JP 2023034671A
- Authority
- JP
- Japan
- Prior art keywords
- error
- ray
- scattered radiation
- index value
- bone
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000010365 information processing Effects 0.000 title claims description 15
- 238000004364 calculation method Methods 0.000 claims abstract description 118
- 210000000988 bone and bone Anatomy 0.000 claims abstract description 85
- 238000003384 imaging method Methods 0.000 claims description 106
- 230000005855 radiation Effects 0.000 claims description 103
- 238000005259 measurement Methods 0.000 claims description 27
- 238000012937 correction Methods 0.000 claims description 23
- 239000012491 analyte Substances 0.000 abstract 3
- 230000006870 function Effects 0.000 description 192
- 238000012545 processing Methods 0.000 description 78
- 230000037182 bone density Effects 0.000 description 64
- 229910052500 inorganic mineral Inorganic materials 0.000 description 31
- 239000011707 mineral Substances 0.000 description 31
- 230000015654 memory Effects 0.000 description 26
- 238000010586 diagram Methods 0.000 description 20
- 238000000034 method Methods 0.000 description 19
- 230000009977 dual effect Effects 0.000 description 15
- 210000004872 soft tissue Anatomy 0.000 description 11
- 238000001739 density measurement Methods 0.000 description 10
- 238000009547 dual-energy X-ray absorptiometry Methods 0.000 description 8
- 230000007246 mechanism Effects 0.000 description 8
- 238000001514 detection method Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 239000000463 material Substances 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000013473 artificial intelligence Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 239000002355 dual-layer Substances 0.000 description 2
- 238000002594 fluoroscopy Methods 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 230000007721 medicinal effect Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 210000000689 upper leg Anatomy 0.000 description 2
- 208000001132 Osteoporosis Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 210000004705 lumbosacral region Anatomy 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/46—Arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
- A61B6/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus 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/505—Apparatus 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 bone
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5294—Devices using data or image processing specially adapted for radiation diagnosis involving using additional data, e.g. patient information, image labeling, acquisition parameters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4241—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5282—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Human Computer Interaction (AREA)
- Physiology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
【課題】被検体の骨の状態について、時系列的な変化を適切に評価できるようにすること。【解決手段】実施形態に係るX線診断装置は、算出部と、表示制御部とを備える。算出部は、異なる2種のエネルギーのX線に対応する被検体の撮影画像及び撮影条件のうち少なくとも一方に基づいて、被検体の骨の状態を評価する指標値の誤差を算出する。表示制御部は、算出された指標値の誤差に基づく情報を表示させる。【選択図】図1
Description
本明細書及び図面に開示の実施形態は、X線診断装置及び医用情報処理装置に関する。
従来、骨密度(BMD:Bone Mineral Density)や骨塩量(BMC:Bone Mineral Content)等、骨の状態を表す指標を測定する技術として、DXA(Dual-energy X-ray Absorptiometry)法が知られている。DXA法では、異なる2種のエネルギーのX線に対応する被検体の撮影画像データを用いて骨と軟組織を弁別し、それによって骨密度や骨塩量を算出する。
骨の状態を表す指標を測定する装置としては、専用のX線診断装置が知られている(以下、専用機と記す)。専用機は、細い短冊状のX線照射領域を、X線照射領域の短辺方向に順次移動させて撮影することで、散乱線の影響を低減している。また、近年、専用機ではないX線診断装置(以下、汎用機と記す)においても、DXA法による骨密度・骨塩量の測定が行われている。
なお、上述した専用機や汎用機は、骨密度・骨塩量を表示するが、かかる表示では、骨密度・骨塩量の時系列的な変化を適切に評価できない。
「大腿骨近位部BMD測定マニュアル」、[令和3年3月1日検索]、インターネット<http://www.josteo.com/ja/guideline/doc/4_1.pdf>
本明細書及び図面に開示の実施形態が解決しようとする課題の一つは、被検体の骨の状態について、時系列的な変化を適切に評価できるようにすることである。ただし、本明細書及び図面に開示の実施形態により解決しようとする課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を他の課題として位置づけることもできる。
実施形態に係るX線診断装置は、算出部と、表示制御部とを備える。算出部は、異なる2種のエネルギーのX線に対応する被検体の撮影画像及び撮影条件のうち少なくとも一方に基づいて、前記被検体の骨の状態を評価する指標値の誤差を算出する。表示制御部は、算出された前記指標値の誤差に基づく情報を表示させる。
以下、図面を参照して、X線診断装置及び医用情報処理装置の実施形態について詳細に説明する。なお、本願に係るX線診断装置及び医用情報処理装置は、以下に示す実施形態によって限定されるものではない。また、以下の説明において、同様の構成要素には共通の符号を付与するとともに、重複する説明を省略する。
(第1の実施形態)
第1の実施形態に係るX線診断装置における構成について説明する。本実施形態に係るX線診断装置は、骨の状態を評価する指標をDXA法によって測定するものである。骨の状態を評価する指標としては、BMD(骨密度)やBMC(骨塩量)等がある。X線診断装置は、骨の状態を測定する専用機のほか、X線TV装置やX線一般撮影装置などの汎用機がある。なお、第1の実施形態では、汎用機であるCアームタイプのX線TV装置を一例として説明する。
第1の実施形態に係るX線診断装置における構成について説明する。本実施形態に係るX線診断装置は、骨の状態を評価する指標をDXA法によって測定するものである。骨の状態を評価する指標としては、BMD(骨密度)やBMC(骨塩量)等がある。X線診断装置は、骨の状態を測定する専用機のほか、X線TV装置やX線一般撮影装置などの汎用機がある。なお、第1の実施形態では、汎用機であるCアームタイプのX線TV装置を一例として説明する。
図1は、第1の実施形態に係るX線診断装置1の構成の一例を示すブロック図である。図1に示すように、X線診断装置1は、X線高電圧装置11と、X線管12と、X線絞り器13と、天板14と、Cアーム15と、X線検出器16と、メモリ17と、ディスプレイ18と、入力インターフェース19と、処理回路20とを備える。
X線高電圧装置11は、処理回路20による制御の下、X線管12に高電圧を印加する。例えば、X線高電圧装置11は、変圧器(トランス)、整流器等の電気回路を有し、X線管12に印加する高電圧を発生する高電圧発生装置と、X線管12が照射するX線に応じた出力電圧の制御を行うX線制御装置とを有する。なお、高電圧発生装置は、変圧器方式であってもよいし、インバータ方式であってもよい。
X線管12は、熱電子を発生する陰極(フィラメント)と、熱電子の衝突を受けてX線を発生する陽極(ターゲット)とを有する真空管である。X線管12は、X線高電圧装置11から印加される高電圧を用いて、陰極から陽極に向けて熱電子を照射することにより、X線を発生する。
X線絞り器13は、X線管12により発生されたX線の照射範囲を絞り込むX線絞りと、X線管12から曝射されたX線を調節するフィルタとを有する。
X線絞り器13におけるX線絞りは、例えば、スライド可能な4枚の絞り羽根を有する。X線絞りは、絞り羽根をスライドさせることで、X線管12が発生したX線を絞り込んで被検体Pに照射させる。ここで、絞り羽根は、鉛などで構成された板状部材であり、X線の照射範囲を調整するためにX線管12のX線照射口付近に設けられる。また、絞り羽根は、対向する羽根が非対称に移動可能に形成されてもよく、或いは、対向する羽根が対称移動のみ可能に形成されてもよい。
X線絞り器13におけるフィルタは、被検体Pに対する被曝線量の低減とX線画像データの画質向上を目的として、その材質や厚みによって透過するX線の線質を変化させ、被検体Pに吸収されやすい軟線成分を低減したり、X線画像データのコントラスト低下を招く高エネルギー成分を低減したりする。また、フィルタは、その材質や厚み、位置などによってX線の線量及び照射範囲を変化させ、X線管12から被検体Pへ照射されるX線が予め定められた分布になるようにX線を減衰させる。
例えば、X線絞り器13は、モータ、アクチュエータ等の駆動機構を有し、処理回路20による制御の下、駆動機構を動作させることによりX線の照射を制御する。例えば、X線絞り器13は、処理回路20から受け付けた制御信号に応じて駆動電圧を駆動機構に印加することにより、X線絞りの絞り羽根の開度を調整して、被検体Pに対して照射されるX線の照射範囲を制御する。また、例えば、X線絞り器13は、処理回路20から受け付けた制御信号に応じて駆動電圧を駆動機構に印加することにより、フィルタの位置を調整することで、被検体Pに対して照射されるX線の線量の分布を制御する。
天板14は、被検体Pを載せるベッドであり、図示しない寝台の上に配置される。例えば、寝台は、モータ、アクチュエータ等の駆動機構を有し、処理回路20による制御の下、駆動機構を動作させることにより、天板14を移動・傾斜させる。
Cアーム15は、X線管12及びX線絞り器13と、X線検出器16とを、被検体Pを挟んで対向するように保持する。例えば、Cアーム15は、モータ、アクチュエータ等の駆動機構を有し、処理回路20から受け付けた制御信号に応じて駆動電圧を駆動機構に印加することにより、X線管12及びX線絞り器13と、X線検出器16とを被検体Pに対して回転・移動させ、X線の照射位置や照射角度を制御する。
X線検出器16は、例えば、マトリクス状に配列された検出素子を有するX線平面検出器(Flat Panel Detector:FPD)である。X線検出器16は、X線管12から照射されて被検体Pを透過したX線を検出して、検出したX線量に対応した検出信号を処理回路20へと出力する。なお、X線検出器16は、グリッド、シンチレータアレイ及び光センサアレイを有する間接変換型の検出器であってもよいし、入射したX線を電気信号に変換する半導体素子を有する直接変換型の検出器であっても構わない。
メモリ17は、例えば、RAM(Random Access Memory)等の半導体メモリ素子により実現される。メモリ17は、処理回路20による処理結果を一時的に記憶する。例えば、メモリ17は、処理回路20によって収集されたX線画像データなどの各種データを受け付けて一時記憶する。ここで、本願におけるX線画像データは、X線検出器16によって検出された検出信号、検出信号に基づいて生成された投影データ、投影データに基づいて生成されたX線画像を含む。また、メモリ17は、処理回路20によって読み出されて実行される各種機能に対応するプログラムを記憶する。
ディスプレイ18は、各種の情報を表示する。例えば、ディスプレイ18は、処理回路20による制御の下、操作者の指示を受け付けるためのGUIや、各種のX線画像を表示する。また、ディスプレイ18は、処理回路20による処理結果を表示する。例えば、ディスプレイ18は、被検体Pの骨の状態を評価する指標(骨密度や骨塩量)の測定値及び当該測定値の誤差等を表示する。
入力インターフェース19は、操作者からの各種の入力操作を受け付け、受け付けた入力操作を電気信号に変換して処理回路20に出力する。例えば、入力インターフェース19は、マウスやキーボード、トラックボール、スイッチ、ボタン、ジョイスティック、操作面へ触れることで入力操作を行うタッチパッド、表示画面とタッチパッドとが一体化されたタッチスクリーン、光学センサを用いた非接触入力回路、音声入力回路等により実現される。なお、入力インターフェース19は、装置本体と無線通信可能なタブレット端末等で構成されることにしても構わない。また、入力インターフェース19は、マウスやキーボード等の物理的な操作部品を備えるものだけに限られない。例えば、装置とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、この電気信号を処理回路20へ出力する電気信号の処理回路も入力インターフェース19の例に含まれる。
処理回路20は、X線診断装置1全体の動作を制御する。また、処理回路20は、メモリ17に記憶されたプログラムを読み出して実行することにより、収集機能201、制御機能202及び算出機能203として機能する。制御機能202は、表示制御部の一例である。また、算出機能203は、算出部の一例である。
ここで、被検体の骨の状態について時系列的な変化を適切に評価することが重要である。しかしながら、過去の検査における測定値と、今回の検査における測定値とを比較する際に、測定された指標値の表示のみでは、その時系列的な変化が有意差と言えるのか否かを一見して把握できず、被検体の骨の状態について時系列的な変化を適切に評価できない。
以上のように構成されたX線診断装置1は、被検体の骨の状態を評価する指標値の他に当該指標値の誤差を表示する。これにより、被検体の骨の状態について時系列的な変化を適切に評価できるようにする。
以下、図2を参照してX線診断装置1が実行する処理を説明する。図2は、X線診断装置1による処理手順の概要を示すフローチャートである。
収集機能201は、デュアルエナジーによる撮影のX線条件を決定する。例えば、収集機能201は、入力インターフェース19から入力されたX線条件、メモリ17に記憶されたX線条件(例えば、前回検査時に設定された)を撮影のX線条件として決定する。収集機能201は、入力インターフェース19を介して操作者によって撮影の実行指示がなされると、決定したX線条件に応じてデュアルエナジーによる撮影を実行する。
なお、収集機能201は、デュアルエナジーによる撮影のX線条件を透視による位置決めに基づいて決定することもできる。かかる場合には、収集機能201は、まず、操作者による透視の実行操作に応じて位置決めを行うための透視画像の収集を順次実行する。ここで、収集機能201は、収集される透視画像について、例えば、画素値の平均値を閾値と比較し、比較結果を次の透視画像のX線条件にフィードバックさせるABC(Automatic Brightness Control)制御を実行する。
そして、収集機能201は、ABC制御によってX線条件が安定した状態に至っている場合に、当該X線条件に基づく体厚の情報を被検体の「体厚に関する情報」としてメモリ17に保持させるとともに、「体厚に関する情報」に基づいてデュアルエナジーによる撮影(異なる2種類の管電圧による撮影)のX線条件を決定する。なお、ABC制御によって安定した状態に至ったX線条件に基づいて、デュアルエナジーによる撮影のX線条件を決定する場合でもよい。収集機能201は、操作者によって撮影の実行指示がなされると、決定したX線条件に応じた撮影を実行する(ステップS101)。
例えば、収集機能201は、腰椎や大腿骨近位部といったROI(Region of Interest:関心領域)を含む領域を対象として第1の管電圧(高電圧)による撮影を行い、第1の管電圧に対応するX線画像データを収集する。また、収集機能201は、上記ROIを含む同一の領域を対象として第2の管電圧(低電圧)による撮影を行い、第2の管電圧に対応するX線画像データを収集する。
なお、異なる2種のエネルギーのX線に対応する被検体の撮影画像を収集するための撮影は、上記したデュアルエナジーによる撮影に限られない。例えば、連続X線エネルギーを分光して低エネルギーのX線と、高エネルギーのX線とを検出する2層検出器を用い、連続X線エネルギーのX線を1回照射する撮影の場合でもよい。
収集機能201が撮影を実行すると、算出機能203は、デュアルエナジーによる撮影で収集された2種類のX線画像データから被検体の骨の状態を評価する指標値を算出する。制御機能202は、算出された指標値(例えば、骨密度や骨塩量等)を表示する(ステップS102)。例えば、算出機能203は、2種類のX線画像データを基に骨画像を生成し、生成した骨画像の関心領域における骨密度や骨塩量等を計測する。制御機能202は、計測された骨密度や骨塩量等をディスプレイ18に表示させる。
また、ステップS102においては、算出機能203は、算出された指標値(骨密度や骨塩量等)の誤差を算出する。制御機能202は、算出された指標値(骨密度や骨塩量等)の誤差を表示する。
以下、ステップS102における処理の一例を、図3を用いて説明する。ここでは、被検体の骨の状態を評価する指標についてBMD(骨密度)を一例として説明する。
算出機能203は、デュアルエナジーによる撮影で収集された2種類のX線画像データから骨画像を生成する。また、算出機能203は、生成した骨画像に対してROIを設定する(ステップS201)。
算出機能203は、操作者によって指定されたROIを設定してもよい。また、算出機能203は、既存のセグメンテーション処理の結果として骨画像から抽出されたROIを設定してもよい。なお、算出機能203によってROIが設定されるX線画像データは、第1の管電圧(高電圧)による撮影で収集された高kV画像、第2の管電圧(低電圧)による撮影で収集された低kV画像、或いは、後述する骨強調画像のいずれの画像でもよい。骨の領域がより明瞭に示されている高kV画像及び骨強調画像を用いることで、より正確なROIの設定を行うことができる。また、骨強調画像が用いられる場合、ROIを設定する前に、後述する物質弁別処理が実行される。
続いて、算出機能203は、第1の管電圧(高電圧)による撮影で収集された投影データ(高kV画像)、及び、第2の管電圧(低電圧)による撮影で収集された投影データ(低kV画像)を基に、骨密度を測定する(ステップS202)。
具体的には、算出機能203は、2種類の投影データそれぞれについて線減弱係数の分布を求め、線減弱係数の分布の各位置(各画素)について、2つの物質(骨と骨以外(軟組織))の線減弱係数と混合量とからなる連立方程式を解くことで、各位置における2つの物質の混合量や混合割合を算出する。また、算出機能203は、各位置における骨と軟組織との混合量に基づいて、骨と軟組織それぞれに対応する2種類の投影データを生成する。なお、収集機能201によって生成された投影データは、メモリ17によって記憶される。
また、算出機能203は、各位置(各画素)について、2つの物質(骨と骨以外(軟組織))の質量減弱係数と密度とからなる連立方程式を解くことで、各位置における2つの物質の密度をそれぞれ算出する。例えば、算出機能203は、ROIにおける骨密度(BMDm)を算出する。
算出機能203は、異なる2種のエネルギーのX線に対応する被検体の撮影画像及び撮影条件の少なくとも一方に基づいて、指標値(骨密度)の誤差を算出する。ここで、算出機能203は、単一のDXA撮影、或いは、2層検出器を用いた単一の撮影において得られた骨密度の測定値の誤差を算出する。なお、第1の実施形態に係る骨密度の誤差は、撮影画像に含まれる散乱線に起因する誤差と、量子ノイズ等に起因する誤差とを有する。
以下、撮影画像及び撮影条件を基に指標値(骨密度)の誤差を算出する例について説明する。なお、撮影条件は、単一のDXA撮影(或いは、2層検出器を用いた単一の撮影)におけるX線条件(管電圧、管電流、パルス幅などを含む)と、幾何学的撮影条件と、被検体の体厚に関する情報とを含む。
(散乱線に起因する誤差について)
算出機能203は、撮影画像に含まれる散乱線に起因する誤差を求める(ステップS203)。撮影画像に含まれる散乱線に起因する誤差は、図4に示すように、散乱線が低減されていない画像(元の画像)を基に測定された骨密度と、散乱線が低減された画像を基に測定された骨密度との差分によって求められる。
算出機能203は、撮影画像に含まれる散乱線に起因する誤差を求める(ステップS203)。撮影画像に含まれる散乱線に起因する誤差は、図4に示すように、散乱線が低減されていない画像(元の画像)を基に測定された骨密度と、散乱線が低減された画像を基に測定された骨密度との差分によって求められる。
より具体的には、算出機能203は、第1の管電圧(高電圧)に対応するX線画像データ(高kV画像)についての散乱線画像(高kV散乱線画像)、第2の管電圧(低電圧)に対応するX線画像データ(低kV画像)についての散乱線画像(低kV散乱線画像)を取得する。散乱線画像は、X線条件(kV及びmAs)と、幾何学的撮影条件と、被検体の体厚に関する情報とに基づいて、各画素における散乱線量(散乱線画像)が取得される。
幾何学的撮影条件は、絞り羽根の位置情報、SID(Source to Image receptor Distance)、SOD(Source-to-Object Distance)(又は検出器・天板間距離)等が挙げられる。絞り羽根の位置情報は、例えば、非対称に移動可能な絞り羽根である場合、4つのパラメータにより表される。また、絞り羽根の位置情報は、対称移動のみ可能な絞り羽根の場合、2つのパラメータにより表される。撮影における幾何的撮影条件は、撮影された高kV画像及び/又は低kV画像のDICOMタグに記録される。
また、被検体の体厚に関する情報は、ABC制御によって取得された「体厚に関する情報」を用いてもよく、或いは、デュアルエナジー撮影によって収集されたX線画像データを用いて算出する場合でもよい。
メモリ17には、例えば、X線条件、幾何学的撮影条件及び体厚と、散乱線との関係を示す散乱線データが予め記憶される。算出機能203は、デュアルエナジー撮影のそれぞれについて、散乱線データをもとに散乱線量を推定して、「高kV散乱線画像」及び「低kV散乱線画像」をそれぞれ取得する。算出機能203は、高電圧での撮影におけるX線条件と、幾何学的撮影条件と、被検体の体厚に関する情報とに基づいて、図4に示す「高kV散乱線画像」を推定する。同様に、算出機能203は、低電圧での撮影におけるX線条件と、幾何学的撮影条件と、被検体の体厚に関する情報とに基づいて、図4に示す「低kV散乱線画像」を推定する。
そして、算出機能203は、高kV画像から高kV散乱線画像を除外し、低kV画像から低kV散乱線画像を除外することにより、各差分画像(差分画像高kV、差分画像低kV)を生成する。算出機能203は、各差分画像から生成された骨画像(散乱線が低減された画像)に対して、ステップS201で設定されたROIの位置・大きさに基づくROIを設定し、設定したROIにおける骨密度を算出する。算出機能203は、散乱線が低減されていない画像のROIにおける骨密度と、散乱線が低減されていない画像のROIにおける骨密度との差分を、撮影画像に含まれる散乱線に起因する誤差として求める。
(量子ノイズ等に起因する誤差について)
図3に戻って、上述したように散乱線に起因する誤差を算出すると、算出機能203は、量子ノイズ等に起因する測定値の誤差を求める(ステップS204)。量子ノイズ等による誤差としては、量子ノイズや回路ノイズによる誤差がある。量子ノイズによる誤差は、X線の照射経路で確率的に生じるX線のばらつきによる誤差であり、回路ノイズによる誤差は、X線診断装置1に含まれる回路で確率的に生じる動作のばらつきによる誤差である。
図3に戻って、上述したように散乱線に起因する誤差を算出すると、算出機能203は、量子ノイズ等に起因する測定値の誤差を求める(ステップS204)。量子ノイズ等による誤差としては、量子ノイズや回路ノイズによる誤差がある。量子ノイズによる誤差は、X線の照射経路で確率的に生じるX線のばらつきによる誤差であり、回路ノイズによる誤差は、X線診断装置1に含まれる回路で確率的に生じる動作のばらつきによる誤差である。
算出機能203は、撮影画像に基づく値のばらつきを評価する値を、測定値の誤差として算出する。例えば、算出機能203は、撮影画像に複数設定されたROI(例えば、椎体)について、ROI毎に測定された骨密度の標準偏差や分散等の統計情報を、量子ノイズ等に起因する誤差として算出してもよい。
また、算出機能203は、量子ノイズ量と回路ノイズ量との和を量子ノイズ等に起因する誤差として算出してもよい。例えば、算出機能203は、X線画像(例えば、デュアルエナジー撮影によって収集された高kV画像或いは低kV画像など)の画素値と回路ノイズ量とに基づいてX線画像における光子数を求め、光子数に基づいて量子ノイズ量を推定する。回路ノイズ量は、X線を照射しない状態で取得された信号に基づいて算出される。
また、算出機能203は、物質弁別処理に基づいて生成された骨強調画像及び軟組織強調画像に基づく量子ノイズ量を、量子ノイズ等に起因する誤差として算出してもよい。例えば、算出機能203は、骨強調画像及び軟組織強調画像を用いたシミュレーションによって各画素に入射した光子数を求め、光子数に基づいて量子ノイズ量を推定する。
このように、量子ノイズ等に起因する誤差が算出されたことに応じて、制御機能202は、算出された前記指標値の誤差に基づく情報を表示させる。例えば、制御機能202は、骨密度の測定値を誤差とともにディスプレイ18に表示させる(ステップS205)。誤差は、撮影画像に含まれる散乱線に起因する誤差と、量子ノイズ等に起因する誤差とを有する。
図5は、ディスプレイ18に表示される情報の一例を示す図である。縦軸は骨密度(BMD)を示し、横軸は検査日(Date)を示す。
例えば、制御機能202は、「Date:a」及び「Date:b」について、骨密度の測定値を表示するとともに、「散乱線に起因する誤差」及び「量子ノイズ等に起因する誤差」の少なくとも一方を加味したエラーバーを表示させる。エラーバーの測定値に対して高い側は、散乱線が含まれることで骨密度が低く見積もられていることから、「散乱線に起因する誤差」を含む。また、エラーバーの測定値に対して高い側及び低い側は、「量子ノイズ等に起因する誤差」を含む。なお、制御機能202によって表示される指標値の誤差に基づく情報は、上記した誤差の値を示す情報に限られず、誤差に基づいて導出される情報の場合でもよい。例えば、制御機能202は、算出された誤差と、誤差の基準値とを比較して、比較結果(例えば、誤差が基準値を超えている等)を表示する場合でもよい。
なお、測定値の誤差を推定する条件を変更した場合の誤差を表示してもよい。以下、撮影条件を基に指標値(骨密度)の誤差を算出する例について説明する。かかる場合には、算出機能203が、幾何学的撮影条件を変更した場合における測定値の誤差を求める。具体的には、算出機能203は、撮影条件における幾何学的撮影条件に基づいて、骨密度の測定値の誤差を求める。例えば、算出機能203は、幾何学的撮影条件に基づいて、ROI内に入射する散乱線の量を推定し、推定した散乱線の量に基づいて、測定値の誤差を推定する。一例を挙げると、算出機能203は、絞り羽根の位置情報、SID、SOD等の条件が変更された場合の「散乱線に起因する誤差」を算出する。制御機能202は、幾何学的撮影条件を変更した場合の測定値の誤差を表示させる。
図6は、幾何学的撮影条件を変更した場合にディスプレイ18に表示される情報の一例を示す図であり、図5における「Date:b」における幾何学的撮影条件を変更した場合の測定値の誤差を表示させる場合について示す。
例えば、制御機能202は、「Date:b」の測定値及び誤差の表示と並べて、幾何学的撮影条件を変更した場合の測定値の誤差(点線で示したエラーバー)を表示させる。なお、幾何学的撮影条件を変更した場合をシミュレーションした測定値の誤差を表示させる態様は、図6に示す例に限られず、その他、種々の態様で表示させることができる。例えば、制御機能202は、実際の条件に基づく誤差と識別可能となるように異なる色や形状で示したエラーバーを測定値に重ねて表示させることができる。また、制御機能202は、幾何学的撮影条件を変更するためのGUIとともに、別ウィンドウで表示させることもできる。
上述した実施形態では、被検体の骨の状態を評価する指標値として、骨密度を用いる場合について説明したが、実施形態はこれに限定されるものではなく、被検体の骨の状態を評価する指標値として、骨塩量を用いる場合でもよい。
上述したように、第1の実施形態によれば、算出機能203は、異なる2種のエネルギーのX線に対応する被検体の撮影画像及び撮影条件のうち少なくとも一方に基づいて、被検体の骨の状態を評価する指標値の誤差を算出する。制御機能202と、算出された指標値の誤差に基づく情報を表示させる。したがって、第1の実施形態に係るX線診断装置1は、骨密度・骨塩量の検査において、1回の検査ごとの測定値の誤差を表示させることができ、被検体の骨の状態について、時系列的な変化を適切に評価できるようにすることを可能にする。
また、第1の実施形態によれば、算出機能203は、撮影条件における幾何学的撮影条件に基づいて、被検体の骨の状態を評価する指標値の誤差を算出する。したがって、第1の実施形態に係るX線診断装置1は、幾何学的な撮影条件を種々変更可能な汎用機によって骨密度及び骨塩量の測定値を算出する場合でも、その誤差を精度よく推定することを可能にする。
例えば、汎用機により骨密度・骨塩量検査を行う場合、その幾何的な撮影条件(絞り羽根の開度、SID、SOD等)によってROIに入射する散乱線の量が変化するため、これに起因して、測定値おける誤差量が変化することとなる。第1の実施形態に係るX線診断装置1は、幾何学的撮影条件に基づいて誤差を推定することができるため、このような汎用機による測定においてもその誤差を精度よく算出することができる。
また、第1の実施形態によれば、算出機能203は、撮影画像における統計情報を、被検体の骨の状態を評価する指標値の誤差として算出する。したがって、第1の実施形態に係るX線診断装置1は、撮影画像に基づいて測定値の誤差を算出することを可能にする。
また、第1の実施形態によれば、算出機能203は、幾何学的撮影条件に基づいて、ROI内に入射する散乱線の量を推定し、推定した散乱線の量に基づいて、指標値の誤差を推定する。したがって、第1の実施形態に係るX線診断装置1は、散乱線による誤差を考慮した指標値の誤差を提示することを可能にする。
また、第1の実施形態によれば、算出機能203は、量子ノイズに基づく誤差を含む誤差を推定する。したがって、第1の実施形態に係るX線診断装置1は、量子ノイズによる誤差を考慮した指標値の誤差を提示することを可能にする。
また、第1の実施形態によれば、算出機能203は、幾何学的撮影条件を変更した場合における指標値の誤差を算出する。制御機能202は、幾何学的撮影条件を変更した場合の指標値の誤差を表示させる。したがって、第1の実施形態に係るX線診断装置1は、指標値の誤差について、幾何学的な撮影条件に関する種々のシミュレーションを実行することを可能にする。
(変形例1)
上述した実施形態では、散乱線補正前の撮影画像(図4における高kV画像及び低kV画像)に基づいて算出した指標値(骨密度)を測定値とし、散乱線補正後の撮影画像(図4における差分画像高kV及び差分画像低kV)に基づいて算出した指標値(骨密度)と測定値との差分を当該測定値の誤差として求める例を説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、散乱線補正後の撮影画像(図4における差分画像高kV及び差分画像低kV)に基づいて算出した指標値(骨密度や骨塩量等)を測定値とし、散乱線補正前の撮影画像に基づいて算出した指標値と測定値との差分を当該測定値の誤差として求めてもよい。かかる場合には、散乱線に起因する誤差が、エラーバーの測定値に対して低い側に含まれることとなる。なお、散乱線補正後の撮影画像に基づいて算出した指標値を表示する場合、散乱線に起因する誤差を含めず、量子ノイズ等に起因する誤差のみを表示する場合でもよい。
上述した実施形態では、散乱線補正前の撮影画像(図4における高kV画像及び低kV画像)に基づいて算出した指標値(骨密度)を測定値とし、散乱線補正後の撮影画像(図4における差分画像高kV及び差分画像低kV)に基づいて算出した指標値(骨密度)と測定値との差分を当該測定値の誤差として求める例を説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、散乱線補正後の撮影画像(図4における差分画像高kV及び差分画像低kV)に基づいて算出した指標値(骨密度や骨塩量等)を測定値とし、散乱線補正前の撮影画像に基づいて算出した指標値と測定値との差分を当該測定値の誤差として求めてもよい。かかる場合には、散乱線に起因する誤差が、エラーバーの測定値に対して低い側に含まれることとなる。なお、散乱線補正後の撮影画像に基づいて算出した指標値を表示する場合、散乱線に起因する誤差を含めず、量子ノイズ等に起因する誤差のみを表示する場合でもよい。
(変形例2)
上述した実施形態では、X線条件、幾何学的撮影条件及び体厚と、散乱線との関係を示す散乱線データに基づいて散乱線量を推定して、「高kV散乱線画像」及び「低kV散乱線画像」をそれぞれ取得する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、AI(Artificial Intelligence)により、「高kV散乱線画像」及び「低kV散乱線画像」をそれぞれ取得する場合でもよい。かかる場合には、例えば、撮影画像と、当該撮影画像に基づく散乱線画像とを学習データとして学習済みモデルが予め生成され、メモリ17に記憶される。算出機能203は、当該学習済みモデルに対して撮影画像(高kV画像及び低kV画像)を入力することで、「高kV散乱線画像」及び「低kV散乱線画像」をそれぞれ取得する。これにより、撮影条件を用いることなく、撮影画像のみを用いて、指標値の誤差を表示させることができる。
上述した実施形態では、X線条件、幾何学的撮影条件及び体厚と、散乱線との関係を示す散乱線データに基づいて散乱線量を推定して、「高kV散乱線画像」及び「低kV散乱線画像」をそれぞれ取得する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、AI(Artificial Intelligence)により、「高kV散乱線画像」及び「低kV散乱線画像」をそれぞれ取得する場合でもよい。かかる場合には、例えば、撮影画像と、当該撮影画像に基づく散乱線画像とを学習データとして学習済みモデルが予め生成され、メモリ17に記憶される。算出機能203は、当該学習済みモデルに対して撮影画像(高kV画像及び低kV画像)を入力することで、「高kV散乱線画像」及び「低kV散乱線画像」をそれぞれ取得する。これにより、撮影条件を用いることなく、撮影画像のみを用いて、指標値の誤差を表示させることができる。
(第2の実施形態)
上述した第1の実施形態では、指標値の誤差が、散乱線に起因する誤差と、量子ノイズ等に起因する誤差とを含む場合について説明した。第2の実施形態では、指標値の誤差が、散乱線量の推定精度上の誤差と、量子ノイズ等に起因する誤差とを含む場合について説明する。なお、第2の実施形態に係るX線診断装置1は、第1の実施形態に係るX線診断装置1と比較して、算出機能203による処理内容が異なる。以下、これを中心に説明する。
上述した第1の実施形態では、指標値の誤差が、散乱線に起因する誤差と、量子ノイズ等に起因する誤差とを含む場合について説明した。第2の実施形態では、指標値の誤差が、散乱線量の推定精度上の誤差と、量子ノイズ等に起因する誤差とを含む場合について説明する。なお、第2の実施形態に係るX線診断装置1は、第1の実施形態に係るX線診断装置1と比較して、算出機能203による処理内容が異なる。以下、これを中心に説明する。
以下、図7を用いて、第2の実施形態に係る処理手順を説明する。なお、図7は、図2におけるステップS102における処理の詳細を示す。ここでは、被検体の骨の状態を評価する指標についてBMD(骨密度)を一例として説明する。
例えば、第2の実施形態に係るX線診断装置1においては、図7に示すように、まず、算出機能203が、第1の実施形態と同様に、収集されたX線画像に対してROIを設定する(ステップS301)。
続いて、算出機能203は、デュアルエナジー撮影で収集された2種類のX線画像データそれぞれについて、各画素での散乱線量と散乱線の推定精度上の誤差をそれぞれ算出する(ステップS302)。
(散乱線の推定精度上の誤差について)
算出機能203は、撮影条件を基に散乱線の推定精度上の誤差を算出する。具体的には、算出機能203は、X線条件や体厚に関する情報において、散乱線量の推定に用いた値と実際の値との乖離などに起因する「散乱線の推定精度上の誤差」を算出する。例えば、条件ごとに、散乱線量として推定されうる範囲を示す推定精度データが予め生成され、メモリ17に記憶される。算出機能203は、デュアルエナジー撮影のそれぞれについて、推定精度データをもとに散乱線量としてとりうる範囲を算出する。第2の実施形態に係る算出機能203は、上記した散乱線量としてとりうる範囲に基づいて、測定値の誤差を算出する。
算出機能203は、撮影条件を基に散乱線の推定精度上の誤差を算出する。具体的には、算出機能203は、X線条件や体厚に関する情報において、散乱線量の推定に用いた値と実際の値との乖離などに起因する「散乱線の推定精度上の誤差」を算出する。例えば、条件ごとに、散乱線量として推定されうる範囲を示す推定精度データが予め生成され、メモリ17に記憶される。算出機能203は、デュアルエナジー撮影のそれぞれについて、推定精度データをもとに散乱線量としてとりうる範囲を算出する。第2の実施形態に係る算出機能203は、上記した散乱線量としてとりうる範囲に基づいて、測定値の誤差を算出する。
例えば、算出機能203は、「高kV画像」の撮影条件に基づいて算出した散乱線量として推定されうる範囲と、「低kV画像」の撮影条件に基づいて算出した散乱線量として推定されうる範囲とを用いた統計処理により、測定値の誤差を推定する。一例を挙げると、算出機能203は、「高kV画像」の撮影条件に基づく範囲と「低kV画像」の撮影条件に基づく範囲とにおいて、最大となる範囲、或いは、最小となる範囲を、測定値の誤差として推定する。
続いて、算出機能203は、散乱線によるX線量のズレを補正して、骨強調画像及び軟組織強調画像を算出する(ステップS303)。具体的には、算出機能203は、高kV画像及び低kV画像からそれぞれ対応する散乱線画像を減算することにより、各差分画像(差分画像高kV、差分画像低kV)を生成する。そして、各差分画像を用いた物質弁別処理に基づいて、骨強調画像と軟組織強調画像とが生成される。なお、高kV画像及び低kV画像からそれぞれ対応する散乱線画像は、第1の実施形態に記載のいずれかの手法を用いて実行される。
そして、算出機能203は、「量子ノイズ等に起因する誤差」を算出する(ステップS304)。なお、このステップは、図3のステップS204と同様に、第1の実施形態に記載のいずれかの手法を用いて実行される。
その後、算出機能203は、「散乱線量の推定精度上の誤差」と「量子ノイズ等に起因する誤差」とに基づいて、骨強調画像及び軟組織強調画像の誤差を算出する(ステップS305)。具体的には、算出機能203は、「散乱線量の推定精度上の誤差」及び「量子ノイズ等に起因する誤差」を用いた、数値的手法又は解析的手法によって、骨強調画像及び軟組織強調画像の誤差を算出する。
そして、算出機能203は、骨密度の測定値及び誤差を算出する(ステップS306)。具体的には、算出機能203は、ステップS303において実行した物質弁別処理に基づいて、骨密度の測定値を算出する。また、算出機能203は、ステップS305において算出した強調画像の誤差に基づいて、骨密度の測定値の誤差を算出する。
上述したように、骨密度の測定値及び誤差が算出されると、制御機能202は、骨密度の測定値及び誤差をディスプレイ18に表示させる(ステップS307)。図8は、ディスプレイ18に表示される情報の一例を示す図である。縦軸は骨密度(BMD)を示し、横軸は検査日(Date)を示す。
例えば、制御機能202は、「Date:a」及び「Date:b」について、骨密度の測定値を表示するとともに、「散乱線量の推定精度上の誤差」と「量子ノイズ等に起因する誤差」とに基づいて算出された骨密度の測定値の誤差を示すエラーバーを表示させる。
なお、制御機能202は、第1の実施形態と同様に、幾何学的撮影条件を変更した場合における測定値の誤差を表示することができる。すなわち、算出機能203は、絞り羽根の位置情報、SID、SOD等の条件が変更された場合の「散乱線量の推定精度上の誤差」を算出する。そして、算出機能203は、算出した「散乱線量の推定精度上の誤差」を用いて骨密度の測定値の誤差を算出する。
制御機能202は、幾何学的撮影条件を変更した場合の測定値の誤差を表示させる。図9は、幾何学的撮影条件を変更した場合にディスプレイ18に表示される情報の一例を示す図であり、図8における「Date:b」における幾何学的撮影条件を変更した場合の測定値の誤差を表示させる場合について示す。例えば、制御機能202は、「Date:b」の測定値及び誤差の表示と並べて、幾何学的撮影条件を変更した場合の測定値の誤差(点線で示したエラーバー)を表示させる。なお、第2の実施形態においても、第1の実施形態と同様に、幾何学的撮影条件を変更した場合をシミュレーションした測定値の誤差を表示させる態様は、種々の態様で表示させることができる。
上述した実施形態では、被検体の骨の状態を評価する指標値として、骨密度を用いる場合について説明したが、実施形態はこれに限定されるものではなく、被検体の骨の状態を評価する指標値として、骨塩量を用いる場合でもよい。
上述したように、第2の実施形態によれば、算出機能203は、散乱線量の推定精度上の誤差を用いて指標値の誤差を推定する。したがって、第2の実施形態に係るX線診断装置1は、散乱線量の推定誤差に応じた誤差を算出することを可能にする。
(第3の実施形態)
第3の実施形態では、絞り羽根によってX線が遮蔽された領域の画素値を用いて散乱線を推定する場合について説明する。すなわち、第3の実施形態では、画像を基に散乱線に起因する誤差を算出する。図10は、第3の実施形態に係るX線診断装置1aの構成の一例を示すブロック図である。なお、第3の実施形態に係るX線診断装置1aは、第1の実施形態に係るX線診断装置1と比較して、処理回路20aが補正機能204を新たに実行する点と、算出機能203による処理内容が異なる。以下、これらを中心に説明する。なお、補正機能204は、補正部の一例である。
第3の実施形態では、絞り羽根によってX線が遮蔽された領域の画素値を用いて散乱線を推定する場合について説明する。すなわち、第3の実施形態では、画像を基に散乱線に起因する誤差を算出する。図10は、第3の実施形態に係るX線診断装置1aの構成の一例を示すブロック図である。なお、第3の実施形態に係るX線診断装置1aは、第1の実施形態に係るX線診断装置1と比較して、処理回路20aが補正機能204を新たに実行する点と、算出機能203による処理内容が異なる。以下、これらを中心に説明する。なお、補正機能204は、補正部の一例である。
補正機能204は、X線絞りによって定められたX線の照射領域以外の領域における画素値に基づいて、ROI内の散乱線の量を補正する。具体的には、補正機能204は、まず、絞り羽根によってX線が遮蔽された領域において検出された検出信号を散乱線に起因する信号として、X線が遮蔽された領域における散乱線量を算出する。
そして、補正機能204は、算出機能203によって推定された散乱線画像に基づいて、X線の照射領域(絞り羽根によってX線が遮蔽されていない領域)における座標依存性の散乱線関数を生成し、生成した散乱線関数を、X線が遮蔽された領域における散乱線量を用いて補正する。例えば、補正機能204は、絞り羽根によるX線照射の境界(X線が照射される位置と遮蔽される位置との境界)において、上記した座標依存性の散乱線関数に基づく散乱線量が連続するように、散乱線関数に定数を乗じる。
そして、補正機能204は、補正された散乱線関数に基づいて、X線の照射領域の各位置における散乱線量を算出する。算出機能203は、補正機能204によって算出された散乱線量を用いて、測定値の誤差を算出する。なお、算出機能203によって推定される散乱線画像は、撮影条件を基に推定されたものでもよく、或いは、AIを用いて推定されたものでもよい。
上述したように、第3の実施形態に係るX線診断装置1aは、絞り羽根によってX線が遮蔽された領域の画素値を用いて散乱線を推定することができる。さらに、第3の実施形態に係るX線診断装置1aでは、ROIの位置に基づいて、散乱線量を推定してもよい。
具体的には、算出機能203は、絞り羽根の開度又はROIから絞り羽根までの距離に応じた散乱線の推定を実行する。例えば、絞り羽根の開度又はROIから絞り羽根までの距離が閾値よりも大きい場合、算出機能203は、第1の実施形態において説明した手法により散乱線量を推定する。ここで、算出機能203は、上記した補正機能204によって算出された補正済みの散乱線量を用いる場合でもよい。
一方、絞り羽根の開度又はROIから絞り羽根までの距離が閾値よりも小さい場合、算出機能203は、X線絞りによって定められたX線の照射領域以外の領域における画素値に基づいて、ROI内の散乱線の量を算出する。すなわち、算出機能203は、絞り羽根によってX線が遮蔽された領域において検出された検出信号に基づいて、ROI内の散乱線量を推定する。例えば、算出機能203は、X線が遮蔽された領域において検出された検出信号に基づいて当該領域における散乱線量を算出して、算出した散乱線量をROI内の散乱線量として用いる。
上述したように、第3の実施形態によれば、補正機能204は、X線絞りによって定められたX線の照射領域以外の領域における画素値に基づいて、ROI内の散乱線の量を補正する。したがって、第3の実施形態に係るX線診断装置1aは、実際に検出される散乱線量に基づいて、推定した散乱線量を補正することができ、より正確な誤差を提示することを可能にする。
また、第3の実施形態によれば、算出機能203は、ROIの位置に基づいて、ROI内の散乱線の量を推定する。したがって、第3の実施形態に係るX線診断装置1aは、ROIの位置に応じた散乱線量の推定を行うことを可能にする。
また、第3の実施形態によれば、算出機能203は、X線絞りによって定められたX線の照射領域以外の領域における画素値に基づいて、ROI内の散乱線の量を算出する。したがって、第3の実施形態に係るX線診断装置1aは、実際に検出される散乱線量を利用した散乱線量の推定を行うことを可能にする。
(その他の実施形態)
これまで第1~第3の実施形態について説明したが、上述した第1~第3の実施形態以外にも、種々の異なる形態にて実施されてよいものである。
これまで第1~第3の実施形態について説明したが、上述した第1~第3の実施形態以外にも、種々の異なる形態にて実施されてよいものである。
上述した実施形態では、デュアルエナジー撮影後に、指標値(骨密度・骨塩量等)の測定値と誤差を表示させる場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、撮影前に撮影条件に基づく誤差を算出して、撮影条件が適切か否かを判定して報知する場合でもよい。
かかる場合には、算出機能203は、単一のDXA撮影に用いられる撮影条件(或いは、2層検出器による撮影に用いられる撮影条件)に基づいて、指標値(骨密度・骨塩量等)の誤差を算出し、算出した誤差と当該誤差に対応する指標値(骨密度・骨塩量等)の過去の測定値とに基づいて、撮影条件に関する判定を実行する。制御機能202は、判定結果を表示させる。
ここで、算出機能203は、上記した判定として、絶対値との比較による撮影条件の判定と、過去の指標値(骨密度・骨塩量等)からの変化の評価に関する条件の判定とを実行することができる。以下、これらを順に説明する。
まず、絶対値との比較による撮影条件の判定を行う場合、算出機能203は、幾何学的撮影条件とX線条件とに基づいて、指標値(骨密度・骨塩量等)の測定値における代表的な誤差量を算出する。例えば、算出機能203は、検査対象の被検体の過去の体厚に関する情報を取得し、取得した過去の体厚に関する情報と、今回の撮影における幾何学的撮影条件とX線条件とに基づいて、指標値(骨密度・骨塩量等)の測定値の誤差を算出する。なお、被検体の体厚の関する情報は、透視による位置決めに基づいて推定されたものが用いられる場合でもよい。また、撮影条件に基づく代表的な誤差量の算出は、撮影条件ごとに代表的な誤差量が対応付けられたデータが予めメモリ17に記憶され、当該データに基づいて実行される場合でもよい。
そして、算出機能203は、算出した代表的な誤差量を過去の指標値(骨密度・骨塩量等)の測定値に加算した場合に、絶対値との比較において問題となり得るか否かを判定する。ここで、絶対値とは、若年成人の骨密度を100%とした場合にどの程度低下しているかを示すYAM(Young Adult Mean)値や、被検体の年齢に対応する年齢層の平均値、骨粗しょう症診断閾値などである。すなわち、算出機能203は、このような絶対値と測定値とを比較する際に、誤差によって比較に問題が生じるか否かを判定する。
また、算出機能203は、撮影条件の変更に応じて代表的な誤差量を再算出し、絶対値との比較における問題が解消するか否かを判定することもできる。なお、上記した代表的な誤差量は、今回測定される指標値(骨密度・骨塩量等)の測定値の過去からの変化量を見込んで、大き目に設定される場合でもよい。
図11は、その他の実施形態に係る算出機能203及び制御機能202による処理の一例を説明するための図である。縦軸は骨密度(BMD)を示し、横軸は検査日(Date)を示す。また、図11は、「Date:b」の検査前に撮影条件に関する判定を実行する場合の例を示す。
例えば、算出機能203は、「Date:b」の検査で用いる撮影条件「パラメータA:1111、パラメータB:2222」に基づいて、代表的な誤差量を算出する。そして、算出機能203は、「Date:a」の検査におけるBMDの測定値に対して代表的な誤差量を付加した際に、誤差の範囲内に閾値「e」が含まれるか否かを判定する。すなわち、算出機能203は、閾値「e」を超えるか否かが誤差によって変化するか否かを判定する。
例えば、算出機能203は、図11の上段の図に示すように、誤差の範囲内に閾値「e」が含まれる場合、閾値との比較において問題となると判定する。制御機能202は、判定結果として、例えば、図11の上段に示すグラフをディスプレイ18に表示させる。ここで、制御機能202は、閾値との比較において問題となると判定された場合に、図11の上段で示すグラフの表示の他、撮影条件が適切ではないことを示す情報(例えば、アラートなど)をさらに表示させることができる。
そして、算出機能203は、図11の下段の図に示すように、撮影条件のうち「パラメータA:1111」を「パラメータA:3333」に変更する処理に応じて、代表的な誤差量を再算出する。さらに、算出機能203は、「Date:a」の検査におけるBMDの測定値に対して代表的な誤差量を付加した際に、誤差の範囲内に閾値「e」が含まれるか否かを判定する。ここで、図11の下段の図に示すように、誤差の範囲内に閾値「e」が含まれていない場合に、算出機能203は、撮影条件が適切であると判定する。制御機能202は、撮影条件の変更に応じて、図11の下段に示すグラフをディスプレイ18に表示させる。
なお、撮影条件の変更は、操作者による入力操作に基づいて実行されてもよく、或いは、算出機能203が、予め設定された情報に基づいて自動で実行する場合でもよい。
次に、過去の指標値(骨密度・骨塩量等)の測定値からの変化の評価に関する条件の判定を実行する場合、算出機能203は、過去の指標値(骨密度・骨塩量等)の測定値と測定の目的とに基づいて、推定される誤差に関する閾値を決定する。すなわち、算出機能203は、過去の測定値と、目的と、推定される誤差に基づいて、撮影条件が適切であるか否かを判定するための閾値(図11における閾値eなど)を決定する。
例えば、算出機能203は、指標値(骨密度・骨塩量等)の測定値が小さいほど、閾値を小さく設定する。また、算出機能203は、測定の目的が健診である場合には閾値を大きく設定し、目的が薬効効果判定である場合には閾値を小さく設定し、薬効効果判定を含まない経過観察の場合には閾値をこれらの中間に設定する。
なお、上述した、絶対値との比較による撮影条件の判定と、過去の指標値(骨密度・骨塩量等)の測定値からの変化の評価に関する条件の判定とは、両方実行される場合でもよく、どちらか一方のみが実行される場合でもよい。
また、上述した実施形態では、指標値(骨密度・骨塩量等)の誤差として、散乱線に起因する誤差と量子ノイズ等に起因する誤差とを含む場合、及び、散乱線量の推定精度上の誤差と量子ノイズ等に起因する誤差とを含む場合について、説明した。しかしながら、実施形態はこれに限定されるものではなく、散乱線に起因する誤差、量子ノイズ等に起因する誤差、或いは、散乱線量の推定精度上の誤差のいずれかを指標値(骨密度・骨塩量等)の誤差とする場合でもよい。
上述した実施形態では、X線診断装置によって各種処理が実行される場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、上述した各処理が医用情報処理装置によって実行される場合でもよい。
図12は、その他の実施形態に係る医用情報処理装置3の構成の一例を示すブロック図である。図12に示すように、医用情報処理装置3は、ネットワーク2を介してX線診断装置1に接続される。そして、医用情報処理装置3は、通信インターフェース31と、メモリ32と、入力インターフェース33と、ディスプレイ34と、処理回路35とを有する。なお、医用情報処理装置3は、例えば、タブレット端末や、ワークステーションなどの情報処理装置である。
通信インターフェース31は、処理回路35に接続され、ネットワークを介して接続されたX線診断装置1等との間で行われる各種データの伝送及び通信を制御する。例えば、通信インターフェース31は、ネットワークカードやネットワークアダプタ、NIC(Network Interface Controller)等によって実現される。
メモリ32は、処理回路35に接続され、各種データを記憶する。例えば、メモリ32は、RAM(Random Access Memory)、フラッシュメモリ等の半導体メモリ素子や、ハードディスク、光ディスク等によって実現される。本実施形態では、メモリ32は、X線診断装置1から受信したX線画像データ(位置決めのために収集された透視画像や、デュアルエネジー撮影によって収集されたX線画像)などを記憶する。また、メモリ32は、処理回路35の処理に用いられる種々の情報や、処理回路35による処理結果等を記憶する。
入力インターフェース33は、種々の設定などを行うためのトラックボール、スイッチボタン、マウス、キーボード、操作面へ触れることで入力操作を行うタッチパッド、表示画面とタッチパッドとが一体化されたタッチモニタ、光学センサを用いた非接触入力回路、及び音声入力回路等によって実現される。入力インターフェース33は、処理回路35に接続されており、操作者から受け付けた入力操作を電気信号に変換して処理回路35に出力する。なお、本明細書において入力インターフェース33は、マウス、キーボードなどの物理的な操作部品を備えるものだけに限られない。例えば、装置とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、この電気信号を制御回路へ出力する処理回路も入力インターフェースの例に含まれる。
ディスプレイ34は、処理回路35に接続され、処理回路35から出力される各種情報及び各種画像を表示する。例えば、ディスプレイ34は、液晶モニタやCRT(Cathode Ray Tube)モニタ、タッチモニタ等によって実現される。例えば、ディスプレイ34は、操作者の指示を受け付けるためのUI(User Interface)や、種々の画像、処理回路35による種々の処理結果を表示する。
処理回路35は、入力インターフェース33を介して操作者から受け付けた入力操作に応じて、医用情報処理装置3が有する各構成要素を制御する。処理回路35は、図12に示すように、例えば、制御機能351と、算出機能352と、補正機能353とを実行する。ここで、例えば、図12に示す処理回路35の構成要素である制御機能351と、算出機能352と、補正機能353が実行する各処理機能は、コンピュータによって実行可能なプログラムの形態でメモリ32内に記録されている。処理回路35は、例えば、プロセッサであり、メモリ32から各プログラムを読み出し、実行することで読み出した各プログラムに対応する機能を実現する。換言すると、各プログラムを読み出した状態の処理回路35は、図12の処理回路35内に示された各機能を有することとなる。
制御機能351は、医用情報処理装置3の全体を制御する。また、制御機能351は、X線診断装置1とのデータの送受信や、上述した制御機能202と同様の処理を実行する。算出機能352は、上述した算出機能203と同様の処理を実行する。補正機能353は、上述した補正機能204と同様の処理を実行する。
各実施形態において説明したX線診断装置においては、各処理機能がコンピュータによって実行可能なプログラムの形態でメモリ17へ記憶されている。処理回路20は、メモリ17からプログラムを読み出して実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読み出した状態の処理回路20は、読み出したプログラムに対応する機能を有することとなる。なお、上記した各実施形態においては、各処理機能が単一の処理回路20によって実現される場合を示したが、実施形態はこれに限られるものではない。例えば、処理回路20は、複数の独立したプロセッサを組み合わせて構成され、各プロセッサが各プログラムを実行することにより各処理機能を実現するものとしても構わない。また、処理回路20が有する各処理機能は、単一又は複数の処理回路に適宜に分散又は統合されて実現されてもよい。
上記説明において用いた「プロセッサ」という文言は、例えば、CPU(Central Processing Unit)、GPU(Graphics Processing Unit)、あるいは、特定用途向け集積回路(Application Specific Integrated Circuit:ASIC)、プログラマブル論理デバイス(例えば、単純プログラマブル論理デバイス(Simple Programmable Logic Device:SPLD)、複合プログラマブル論理デバイス(Complex Programmable Logic Device:CPLD)、及びフィールドプログラマブルゲートアレイ(Field Programmable Gate Array:FPGA))等の回路を意味する。プロセッサはストレージ111に保存されたプログラムを読み出し実行することで機能を実現する。
なお、上述した各実施形態においては、メモリ17が各処理機能に対応するプログラムを記憶するものとして説明した。しかしながら、複数のメモリ17を分散して配置し、処理回路20は、個別のメモリ17から対応するプログラムを読み出す構成としても構わない。また、メモリ17にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。
上述した実施形態に係る各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。更に、各装置にて行われる各処理機能は、その全部又は任意の一部が、CPU及び当該CPUにて解析実行されるプログラムにて実現され、あるいは、ワイヤードロジックによるハードウェアとして実現されうる。
また、上述した実施形態で説明した制御方法は、予め用意された制御プログラムをパーソナルコンピュータやワークステーション等のコンピュータで実行することによって実現することができる。この制御プログラムは、インターネット等のネットワークを介して配布することができる。また、この制御プログラムは、ハードディスク、フレキシブルディスク(FD)、CD-ROM、MO、DVD等のコンピュータで読み取り可能な非一過性の記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。
以上説明した少なくとも一つの実施形態によれば、被検体の骨の状態について、時系列的な変化を適切に評価できるようにすることができる。
いくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、実施形態同士の組み合わせを行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
1、1a X線診断装置
3 医用情報処理装置
20、20a、35 処理回路
202、351 制御機能
203、352 算出機能
204、353 補正機能
3 医用情報処理装置
20、20a、35 処理回路
202、351 制御機能
203、352 算出機能
204、353 補正機能
Claims (11)
- 異なる2種のエネルギーのX線に対応する被検体の撮影画像及び撮影条件のうち少なくとも一方に基づいて、前記被検体の骨の状態を評価する指標値の誤差を算出する算出部と、
算出された前記指標値の誤差に基づく情報を表示させる表示制御部と、
を備える、X線診断装置。 - 前記算出部は、前記撮影条件における幾何学的撮影条件に基づいて、前記被検体の骨の状態を評価する指標値の誤差を算出する、請求項1に記載のX線診断装置。
- 前記算出部は、前記撮影画像における統計情報を、前記被検体の骨の状態を評価する指標値の誤差として算出する、請求項1又は2に記載のX線診断装置。
- 前記算出部は、前記幾何学的撮影条件に基づいて、関心領域内に入射する散乱線の量を推定し、推定した散乱線の量に基づいて、前記指標値の誤差を算出する、請求項2に記載のX線診断装置。
- 前記算出部は、量子ノイズに基づく誤差を含む前記誤差を推定する、請求項2~4のいずれか1つに記載のX線診断装置。
- 前記算出部は、関心領域の位置に基づいて、関心領域内の散乱線の量を推定し、推定した散乱線の量に基づいて、前記指標値の誤差を算出する、請求項2に記載のX線診断装置。
- 前記算出部は、前記撮影条件における幾何学的撮影条件を変更した場合における前記指標値の誤差を算出し、
前記表示制御部は、前記幾何学的撮影条件を変更した場合の前記指標値の誤差を表示させる、請求項2、4~6のいずれか1つに記載のX線診断装置。 - X線絞りによって定められたX線の照射領域以外の領域における画素値に基づいて、関心領域内の散乱線の量を補正する補正部をさらに備え、
前記算出部は、補正後の散乱線の量に基づいて、前記指標値の誤差を算出する、請求項4~6のいずれか1つに記載のX線診断装置。 - 前記算出部は、X線絞りによって定められたX線の照射領域以外の領域における画素値に基づいて、関心領域内の散乱線の量を算出する、請求項1に記載のX線診断装置。
- 前記算出部は、前記撮影画像の撮影条件に基づいて、前記指標値の誤差を算出し、算出した誤差と当該誤差に対応する過去の測定値とに基づいて、前記撮影画像の撮影に用いられる条件に関する判定を実行し、
前記表示制御部は、判定結果を表示させる、請求項1~9のいずれか1つに記載のX線診断装置。 - 異なる2種のエネルギーのX線に対応する被検体の撮影画像及び撮影条件のうち少なくとも一方に基づいて、前記被検体の骨の状態を評価する指標値の誤差を算出する算出部と、
算出された前記指標値の誤差に基づく情報を表示させる表示制御部と、
を備える、医用情報処理装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2021141010A JP2023034671A (ja) | 2021-08-31 | 2021-08-31 | X線診断装置及び医用情報処理装置 |
US17/822,093 US20230069062A1 (en) | 2021-08-31 | 2022-08-24 | X-ray diagnosis apparatus, medical information processing apparatus and method, and storage medium |
CN202211031851.3A CN115721327A (zh) | 2021-08-31 | 2022-08-26 | X射线诊断装置、医用信息处理装置和方法以及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2021141010A JP2023034671A (ja) | 2021-08-31 | 2021-08-31 | X線診断装置及び医用情報処理装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2023034671A true JP2023034671A (ja) | 2023-03-13 |
Family
ID=85286221
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2021141010A Pending JP2023034671A (ja) | 2021-08-31 | 2021-08-31 | X線診断装置及び医用情報処理装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20230069062A1 (ja) |
JP (1) | JP2023034671A (ja) |
CN (1) | CN115721327A (ja) |
-
2021
- 2021-08-31 JP JP2021141010A patent/JP2023034671A/ja active Pending
-
2022
- 2022-08-24 US US17/822,093 patent/US20230069062A1/en active Pending
- 2022-08-26 CN CN202211031851.3A patent/CN115721327A/zh active Pending
Also Published As
Publication number | Publication date |
---|---|
CN115721327A (zh) | 2023-03-03 |
US20230069062A1 (en) | 2023-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10235766B2 (en) | Radiographic image analysis device and method, and storage medium having stored therein program | |
US8295570B2 (en) | Methods and apparatus for measuring body circumference | |
CN107007294A (zh) | X射线图像诊断装置以及骨密度测量方法 | |
US20220287665A1 (en) | Estimation device, estimation method, and estimation program | |
US20110060247A1 (en) | Methods and apparatus for measuring bone lengths | |
RU2711250C1 (ru) | Система кт-визуализации и способ для системы кт-визуализации | |
JP2014239840A (ja) | 放射線画像診断装置 | |
JP5388476B2 (ja) | 画像表示装置、乳房x線撮影装置および画像表示プログラム | |
US20210279918A1 (en) | Medical information processing apparatus, medical image diagnosis apparatus, and medical information processing method | |
US20030048873A1 (en) | Method for improving a radiological examination and device therefor | |
US20220358652A1 (en) | Image processing apparatus, radiation imaging apparatus, image processing method, and storage medium | |
JP2023034671A (ja) | X線診断装置及び医用情報処理装置 | |
US20220079543A1 (en) | Radiation image processing device and radiation image processing method | |
US20230136752A1 (en) | System and method for imaging a subject | |
US11596374B2 (en) | X-ray diagnostic apparatus and medical-information processing apparatus | |
CN110881993B (zh) | X射线ct装置、医用图像处理装置以及x射线ct系统 | |
EP3673807A1 (en) | Medical processing apparatus, x-ray diagnosis apparatus, and medical processing method | |
US20230104524A1 (en) | X-ray image processing apparatus, x-ray diagnosis apparatus, method, and storage medium | |
US20230309942A1 (en) | Radiation imaging apparatus, radiation imaging method, and non-transitory computer-readable storage medium | |
US11191512B2 (en) | Radiographic image processing device and radiographic image processing method | |
US20240020849A1 (en) | Radiation image processing apparatus, operation method of radiation image processing apparatus, and non-transitory computer readable medium | |
JP2023056167A (ja) | 骨情報測定装置、通知装置、x線診断装置、骨情報測定方法及びプログラム | |
US11937970B2 (en) | System and method for calibrating a camera feature detection system of an x-ray system | |
JP7321715B2 (ja) | トモシンセシス用ファントム及び測定方法 | |
JP7504739B2 (ja) | 医用画像表示装置及び医用画像表示方法 |