JP2019054991A - 解析装置及び解析システム - Google Patents
解析装置及び解析システム Download PDFInfo
- Publication number
- JP2019054991A JP2019054991A JP2017180792A JP2017180792A JP2019054991A JP 2019054991 A JP2019054991 A JP 2019054991A JP 2017180792 A JP2017180792 A JP 2017180792A JP 2017180792 A JP2017180792 A JP 2017180792A JP 2019054991 A JP2019054991 A JP 2019054991A
- Authority
- JP
- Japan
- Prior art keywords
- moving image
- analysis
- average value
- moving
- photographing
- 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.)
- Withdrawn
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 212
- 238000003384 imaging method Methods 0.000 claims abstract description 170
- 238000012937 correction Methods 0.000 claims description 75
- 230000008859 change Effects 0.000 claims description 55
- 230000033001 locomotion Effects 0.000 claims description 30
- 230000005855 radiation Effects 0.000 description 56
- 238000012545 processing Methods 0.000 description 32
- 210000004072 lung Anatomy 0.000 description 18
- 238000003745 diagnosis Methods 0.000 description 12
- 238000004891 communication Methods 0.000 description 11
- 238000000034 method Methods 0.000 description 11
- 238000001514 detection method Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 238000002595 magnetic resonance imaging Methods 0.000 description 6
- 230000017531 blood circulation Effects 0.000 description 5
- 238000010191 image analysis Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 239000004065 semiconductor Substances 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000001678 irradiating effect Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000009423 ventilation Methods 0.000 description 3
- 230000004071 biological effect Effects 0.000 description 2
- 238000013170 computed tomography imaging Methods 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 2
- 238000005206 flow analysis Methods 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 230000002685 pulmonary effect Effects 0.000 description 2
- 230000000241 respiratory effect Effects 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 210000000577 adipose tissue Anatomy 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000037396 body weight Effects 0.000 description 1
- 230000037182 bone density Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000004217 heart function Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012933 kinetic analysis Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004660 morphological change Effects 0.000 description 1
- 230000003387 muscular Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 210000003456 pulmonary alveoli Anatomy 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000012950 reanalysis Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
- YDLQKLWVKKFPII-UHFFFAOYSA-N timiperone Chemical compound C1=CC(F)=CC=C1C(=O)CCCN1CCC(N2C(NC3=CC=CC=C32)=S)CC1 YDLQKLWVKKFPII-UHFFFAOYSA-N 0.000 description 1
- 229950000809 timiperone Drugs 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
【課題】異なる撮影装置で得られた動画像の動態解析結果を比較可能な解析装置を提供する。【解決手段】異なる複数の撮影装置1a〜1cから動画像を受信し、受信した動画像や動画像の解析結果を表示する解析装置3であって、動画像に含まれる構造物Sの移動量を用いて、動画像又は動画像の解析結果を補正する制御部31を備える。【選択図】図1
Description
本発明は、解析装置及び解析システムに関する。
従来のFPD(Flat Panel Detector)等の半導体イメージセンサーを利用した静止画撮影と撮影された静止画による診断に対し、近年ではFPD等を利用した被写体の動画像撮影及びこれにより得られた動画像を解析して診断に提供する解析システムが知られている。
解析システムは、動画像撮影を行う撮影装置と、撮影装置を制御したり撮影された動画像を確認のために表示したりするためのコンソールと、撮影により得られた複数のフレーム画像に対する解析を行う解析装置と、解析結果を表示するための表示装置と、を備えて構成される。
解析システムは、動画像撮影を行う撮影装置と、撮影装置を制御したり撮影された動画像を確認のために表示したりするためのコンソールと、撮影により得られた複数のフレーム画像に対する解析を行う解析装置と、解析結果を表示するための表示装置と、を備えて構成される。
例えば特許文献1には、胸部の動態を撮影して得られた動態画像の複数のフレーム画像間における濃度変化に係る特徴量と、形態変化に係る特徴量とを算出し、濃度変化に係る特徴量と形態変化に係る特徴量の比に基づいて、肺換気機能又は心機能に係る特徴量を算出する解析システムが記載されている。
また、特許文献2には、上記動画撮影を行う撮影装置が、磁気共鳴イメージング装置(以下、MRI(MagneticResonance Imaging)装置という)である場合が例示されている。
また、特許文献2には、上記動画撮影を行う撮影装置が、磁気共鳴イメージング装置(以下、MRI(MagneticResonance Imaging)装置という)である場合が例示されている。
動画像撮影の際に、FPDではマトリクス状に配置された素子に照射される放射線強度に応じた電気信号として各フレーム画像を生成する。システムによっては生成した画像をFPDもしくはコンソール上で対数変換した後、より被写体内の構造物を観察し易くなるよう、LUT等で濃度諧調を変換した動画像をモニター上に表示している。
動態解析では、こうして得られた動画像を構成する複数フレームに含まれる画素値の間で各種演算処理を実施するため、解析元となる動画像の濃度値や動画像の解析結果の信号値が、動態解析の結果に大きく影響することが知られている。
なお、ここで言う「動画像撮影」や「動画像」とは、ビデオカメラ等を用いて行われる通常の動画像撮影や動画像に限定されず、いわば制約付き「動画像撮影」「動画像」と言うべき動態撮影を含む概念である。すなわち、本発明の場合の動画像撮影とは、複数のフレームの放射線画像を時間的に連続して動画状に撮影する動態撮影を含み、通常の動画像撮影におけるフレームレート(例えば1秒間に30フレーム等)よりも小さいフレームレートである場合が少なくなく、また、撮影時間に制約がある場合もある。また、通常の動画像撮影では、通常、撮影した動画像をリアルタイムで表示するように構成されるが、上記の動態撮影のような動画像撮影の場合には、動画像の表示のリアルタイム性は保証されず、撮影と同時に動画像を表示しない(或いは表示できない)ように構成される場合も少なくない。
なお、例えば動態撮影では、放射線照射装置から放射線がパルス状に照射される場合があるが、低線量の放射線が照射され続ける状態で照射される場合もある。そして、放射線画像撮影装置は、通常、前者の場合には、放射線照射装置からの放射線の照射のタイミングと同期させて画像データの読み出し処理を行うように構成され、また、後者の場合には、画像データを読み出すタイミングが適宜調整されて撮影が行われる。
本発明においては、上記何れの場合も「動画像撮影」に含めることができ、時間的に連続して動画状に撮影された複数のフレームの放射線画像を「動画像」と称し、各フレームの放射線画像を「フレーム画像」と称している。
本発明においては、上記何れの場合も「動画像撮影」に含めることができ、時間的に連続して動画状に撮影された複数のフレームの放射線画像を「動画像」と称し、各フレームの放射線画像を「フレーム画像」と称している。
ところで、比較的規模が大きな病院等施設の放射線撮影装置としては、FPD等を利用した透視機能を含む一般撮影装置や血管撮影装置やX線TVの他、CT、MRI等複数の種類の撮影装置が導入されている場合がある。そして、各装置はそれぞれに搭載されたFPD等の放射線検出器や放射線検出器の動作に関わる電気回路とその制御方法等の影響を受けて生成される動画像を、医師等に見やすくすること等を目的として、それぞれ独自に階調処理等の画像処理を施すように構成されており、最終的に各装置で得られる動画像の特性は異なるものとなる。
そのため、例えば、生成される動画像の特性が既知である撮影装置で得られた動画像に対して最適となるパラメーターを用いて動態解析を行う解析装置において、生成される動画像の特性が未知である異なる撮影装置で得られた動画像に対して同じ方法で動態解析を行っても、通常、予期した通りの解析結果は得られない。そのため、各撮影装置で得られた動画像の解析結果同士を比較して診断等に用いようとしても、正しい比較結果が得られない場合が多い。
本発明は、上記の問題点を鑑みてなされたものであり、異なる複数の撮影装置で撮影された動画像の動態解析結果を比較することのできる解析装置を提供することを目的とする。
前記の問題を解決するために、本発明の解析装置は、
異なる複数の撮影装置から動画像を受信し、受信した動画像や動画像の解析結果を表示する解析装置であって、
前記動画像に含まれる構造物の移動量を用いて、前記動画像又は前記動画像の解析結果を補正する補正手段を備えることを特徴とする。
異なる複数の撮影装置から動画像を受信し、受信した動画像や動画像の解析結果を表示する解析装置であって、
前記動画像に含まれる構造物の移動量を用いて、前記動画像又は前記動画像の解析結果を補正する補正手段を備えることを特徴とする。
前記の問題を解決するために、本発明の解析システムは、
動画像を撮影可能である基準となる撮影装置と、
動画像を撮影可能であって前記基準となる撮影装置とは異なる撮影装置と、
請求項1から請求項5のいずれか一項に記載の解析装置と、
を備えることを特徴とする。
動画像を撮影可能である基準となる撮影装置と、
動画像を撮影可能であって前記基準となる撮影装置とは異なる撮影装置と、
請求項1から請求項5のいずれか一項に記載の解析装置と、
を備えることを特徴とする。
本発明によれば、異なる複数の撮影装置で撮影された動画像の動態解析結果を比較可能となる。
以下、本発明に係る解析装置の一実施形態について、図面を参照して説明する。
ただし、以下に述べる実施形態には、本発明を実施するために技術的に好ましい種々の限定が付されているが、本発明の技術的範囲を以下の実施形態および図示例に限定するものではない。
ただし、以下に述べる実施形態には、本発明を実施するために技術的に好ましい種々の限定が付されているが、本発明の技術的範囲を以下の実施形態および図示例に限定するものではない。
[解析装置の位置づけ及び解析システムの全体構成]
本実施形態の解析装置は、解析システムの一部を構成し、複数の異なる撮影装置から動画像を受信し、受信した動画像や動画像の解析結果を表示するものである。
まず、前提として、本実施形態で想定される解析システムの全体構成につき、図1を参照しつつ説明する。
本実施形態の解析装置は、解析システムの一部を構成し、複数の異なる撮影装置から動画像を受信し、受信した動画像や動画像の解析結果を表示するものである。
まず、前提として、本実施形態で想定される解析システムの全体構成につき、図1を参照しつつ説明する。
解析システムは、動画像を撮影可能である異なる複数の撮影装置と、これら異なる複数の撮影装置で撮影された動画像に対して動態解析を行うことが可能な解析装置とを備えている。
本実施形態における解析システム100は、具体的には、図1に示すように、複数の異なる撮影装置1(撮影装置1a,1b,1c)と、解析装置3とがLAN(Local Area Network)等の通信ネットワークNTを介して接続されて構成されている。
複数の撮影装置1は、例えば図示しない一般撮影室に設けられた撮影装置1と、当該撮影装置1を制御するための撮影用コンソール2と、図示しない救急室に設けられた撮影装置1と、当該撮影装置1を制御するための撮影用コンソール2と、可搬型に構成された撮影装置1と、これを制御するための可搬型の撮影用コンソール2等である。
なお、図1では3台の撮影装置1(撮影装置1a,1b,1c)が設けられている例を示しているが、撮影装置1の数は3台に限定されず、2台以上であればよい。さらに、複数の撮影装置1について特に区別しない場合には、単に撮影装置1という。
本実施形態において解析装置3は、診断に供する画像等(動画像や動画像の解析結果)を生成し表示させる診断用コンソールである。
解析システム100を構成する各装置は、DICOM(Digital Image and Communications in Medicine)規格に準じており、各装置間の通信は、DICOMに則って行われる。
なお、解析装置3や撮影装置1は常時通信ネットワークNTに接続されている必要はない。
本実施形態における解析システム100は、具体的には、図1に示すように、複数の異なる撮影装置1(撮影装置1a,1b,1c)と、解析装置3とがLAN(Local Area Network)等の通信ネットワークNTを介して接続されて構成されている。
複数の撮影装置1は、例えば図示しない一般撮影室に設けられた撮影装置1と、当該撮影装置1を制御するための撮影用コンソール2と、図示しない救急室に設けられた撮影装置1と、当該撮影装置1を制御するための撮影用コンソール2と、可搬型に構成された撮影装置1と、これを制御するための可搬型の撮影用コンソール2等である。
なお、図1では3台の撮影装置1(撮影装置1a,1b,1c)が設けられている例を示しているが、撮影装置1の数は3台に限定されず、2台以上であればよい。さらに、複数の撮影装置1について特に区別しない場合には、単に撮影装置1という。
本実施形態において解析装置3は、診断に供する画像等(動画像や動画像の解析結果)を生成し表示させる診断用コンソールである。
解析システム100を構成する各装置は、DICOM(Digital Image and Communications in Medicine)規格に準じており、各装置間の通信は、DICOMに則って行われる。
なお、解析装置3や撮影装置1は常時通信ネットワークNTに接続されている必要はない。
本実施形態において、撮影装置1aは、撮影用コンソール2による撮影制御の下で動作し、放射線を照射する管球等の設定条件や、FPDの画像読み取り条件等が予め解析システム100側において既知であり、当該撮影装置1aによって得られる動画像の特性が解析装置3において既知であるものとする。
このため、本実施形態では、撮影装置1aを「基準となる撮影装置」とし、撮影装置1aによって撮影された動画像を、動画像の特性が既知である「基準動画像D」とする。
また、撮影装置1a以外の撮影装置である撮影装置1b,1cを「基準となる撮影装置とは異なる撮影装置」とし、撮影装置1b,1cによって撮影された動画像を「動画像d」とする。
なお、基準動画像Dや動画像dは、それぞれ各フレーム画像の画像データの集合であり、以下では、基準動画像Dあるいは動画像dの文言を、各フレーム画像の画像データの集合の意味で用いる場合もあるが、動画像を構成する各フレーム画像のうちの1フレーム分のフレーム画像の画像データの意味で用いる場合もある。
このため、本実施形態では、撮影装置1aを「基準となる撮影装置」とし、撮影装置1aによって撮影された動画像を、動画像の特性が既知である「基準動画像D」とする。
また、撮影装置1a以外の撮影装置である撮影装置1b,1cを「基準となる撮影装置とは異なる撮影装置」とし、撮影装置1b,1cによって撮影された動画像を「動画像d」とする。
なお、基準動画像Dや動画像dは、それぞれ各フレーム画像の画像データの集合であり、以下では、基準動画像Dあるいは動画像dの文言を、各フレーム画像の画像データの集合の意味で用いる場合もあるが、動画像を構成する各フレーム画像のうちの1フレーム分のフレーム画像の画像データの意味で用いる場合もある。
[撮影装置1の構成]
撮影装置1は、前述のように、動画像を撮影可能なものである。
本実施形態において「動画像」の撮影とは、被写体に対し、X線等の放射線をパルス状にして所定時間間隔で繰り返し照射するか(パルス照射)、もしくは、低線量率にして途切れなく継続して照射する(連続照射)ことで、被写体の動態を示す複数の画像を取得する動態撮影をいう。すなわち本実施形態における「動画像」とは、このような撮影により得られた、被写体の動態を示す複数の一連の画像を意味する。また、この動画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。
「動画像」の撮影によって撮影可能な被写体(生体)の動態としては、例えば、呼吸運動に伴う肺の膨張及び収縮の形態変化、心臓の拍動等、周期性(サイクル)をもつ人体の動態がある。
以下の実施形態では、パルス照射により胸部(特に肺野)の動態撮影を行う場合を例にとり説明する。
撮影装置1は、前述のように、動画像を撮影可能なものである。
本実施形態において「動画像」の撮影とは、被写体に対し、X線等の放射線をパルス状にして所定時間間隔で繰り返し照射するか(パルス照射)、もしくは、低線量率にして途切れなく継続して照射する(連続照射)ことで、被写体の動態を示す複数の画像を取得する動態撮影をいう。すなわち本実施形態における「動画像」とは、このような撮影により得られた、被写体の動態を示す複数の一連の画像を意味する。また、この動画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。
「動画像」の撮影によって撮影可能な被写体(生体)の動態としては、例えば、呼吸運動に伴う肺の膨張及び収縮の形態変化、心臓の拍動等、周期性(サイクル)をもつ人体の動態がある。
以下の実施形態では、パルス照射により胸部(特に肺野)の動態撮影を行う場合を例にとり説明する。
本実施形態において、例えば撮影装置1aは、図1に示すように、放射線源11、放射線照射制御装置12、放射線検出部13、読取制御装置14等を備える放射線画像撮影装置である。
後述するように、撮影装置1の放射線照射制御装置12及び放射線検出部13は撮影用コンソール2に接続されており、放射線照射制御装置12は撮影用コンソール2から入力された放射線照射条件に基づいて放射線源11を制御して放射線撮影を行う。また、放射線検出部13は、撮影用コンソール2から入力された画像読取条件に基づいて各画素のスイッチング部を制御して、各画素に蓄積された電気信号を読み取ることにより画像データを取得し、取得したフレーム画像の画像データを撮影用コンソール2に出力する。
後述するように、撮影装置1の放射線照射制御装置12及び放射線検出部13は撮影用コンソール2に接続されており、放射線照射制御装置12は撮影用コンソール2から入力された放射線照射条件に基づいて放射線源11を制御して放射線撮影を行う。また、放射線検出部13は、撮影用コンソール2から入力された画像読取条件に基づいて各画素のスイッチング部を制御して、各画素に蓄積された電気信号を読み取ることにより画像データを取得し、取得したフレーム画像の画像データを撮影用コンソール2に出力する。
図1では、撮影装置1aについてのみ詳細な構成を示しているが、撮影装置1b,1cについても同様の構成を備えている。
なお、撮影装置1としては、放射線画像撮影装置の他、CT、MRI等を適用することも可能である。また、X線テレビ等であってもよい。
例えば、撮影装置1a及び撮影装置1bは放射線画像撮影装置であり、撮影装置1cはX線テレビである、というように異なる種類の装置が解析システム内に混在していてもよい。
なお、前述の「異なる複数の撮影装置」とは、複数の撮影装置が異なる個体であることを意味し、撮影装置1a,1b,1cがすべて同じ種類の装置(例えば全て放射線画像撮影装置)である場合を除外するものではない。同じ種類の撮影装置であっても、その設定等によって個体差が生じるため、後述する本解析装置及び解析システムにおける処理を行う意義が存する。
なお、撮影装置1としては、放射線画像撮影装置の他、CT、MRI等を適用することも可能である。また、X線テレビ等であってもよい。
例えば、撮影装置1a及び撮影装置1bは放射線画像撮影装置であり、撮影装置1cはX線テレビである、というように異なる種類の装置が解析システム内に混在していてもよい。
なお、前述の「異なる複数の撮影装置」とは、複数の撮影装置が異なる個体であることを意味し、撮影装置1a,1b,1cがすべて同じ種類の装置(例えば全て放射線画像撮影装置)である場合を除外するものではない。同じ種類の撮影装置であっても、その設定等によって個体差が生じるため、後述する本解析装置及び解析システムにおける処理を行う意義が存する。
放射線源11は、被写体M(被検者)を挟んで放射線検出部13と対向する位置に配置され、放射線照射制御装置12の制御に従って、被写体Mに対し放射線(X線)を照射する。
放射線照射制御装置12は、撮影用コンソール2に接続されており、撮影用コンソール2から入力された放射線照射条件に基づいて放射線源11を制御して放射線撮影を行う。撮影用コンソール2から入力される放射線照射条件は、例えば、パルスレート、パルス幅、パルス間隔、1撮影あたりの撮影フレーム数、X線管電流の値、X線管電圧の値、付加フィルター種等である。パルスレートは、1秒あたりの放射線照射回数であり、後述するフレームレートと一致している。パルス幅は、放射線照射1回当たりの放射線照射時間である。パルス間隔は、1回の放射線照射開始から次の放射線照射開始までの時間であり、後述するフレーム間隔と一致している。
放射線照射制御装置12は、撮影用コンソール2に接続されており、撮影用コンソール2から入力された放射線照射条件に基づいて放射線源11を制御して放射線撮影を行う。撮影用コンソール2から入力される放射線照射条件は、例えば、パルスレート、パルス幅、パルス間隔、1撮影あたりの撮影フレーム数、X線管電流の値、X線管電圧の値、付加フィルター種等である。パルスレートは、1秒あたりの放射線照射回数であり、後述するフレームレートと一致している。パルス幅は、放射線照射1回当たりの放射線照射時間である。パルス間隔は、1回の放射線照射開始から次の放射線照射開始までの時間であり、後述するフレーム間隔と一致している。
放射線検出部13は、FPD(Flat Panel Detector)等の半導体イメージセンサーにより構成される。FPDは、例えば、ガラス基板等を有しており、基板上の所定位置に、放射線源11から照射されて少なくとも被写体Mを透過した放射線をその強度に応じて検出し、検出した放射線を電気信号に変換して蓄積する複数の検出素子(画素)がマトリクス状に配列されている。各画素は、例えばTFT(Thin Film Transistor)等のスイッチング部を備えて構成されている。FPDにはX線をシンチレーターを介して光電変換素子により電気信号に変換する間接変換型、X線を直接的に電気信号に変換する直接変換型があるが、何れを用いてもよい。
放射線検出部13は、被写体Mを挟んで放射線源11と対向するように設けられている。
放射線検出部13は、被写体Mを挟んで放射線源11と対向するように設けられている。
読取制御装置14は、撮影用コンソール2に接続されている。読取制御装置14は、撮影用コンソール2から入力された画像読取条件に基づいて放射線検出部13の各画素のスイッチング部を制御して、当該各画素に蓄積された電気信号の読み取りをスイッチングしていき、放射線検出部13に蓄積された電気信号を読み取ることにより、画像データを取得する。この画像データがフレーム画像である。フレーム画像の画素信号値は濃度値を表す。そして、読取制御装置14は、取得したフレーム画像を撮影用コンソール2に出力する。画像読取条件は、例えば、フレームレート、フレーム間隔、画素サイズ、画像サイズ(マトリックスサイズ)等である。フレームレートは、1秒あたりに取得するフレーム画像数であり、パルスレートと一致している。フレーム間隔は、1回のフレーム画像の取得動作開始から次のフレーム画像の取得動作開始までの時間であり、パルス間隔と一致している。
ここで、放射線照射制御装置12と読取制御装置14は互いに接続され、互いに同期信号をやりとりして放射線照射動作と画像の読み取りの動作を同調させるようになっている。
ここで、放射線照射制御装置12と読取制御装置14は互いに接続され、互いに同期信号をやりとりして放射線照射動作と画像の読み取りの動作を同調させるようになっている。
[撮影用コンソール2の構成]
撮影用コンソール2は、放射線照射条件や画像読取条件を撮影装置1に出力して撮影装置1による放射線撮影及び放射線画像の読み取り動作を制御するとともに、撮影装置1により取得された動態画像を撮影技師等の撮影実施者によるポジショニングの確認や診断に適した画像であるか否かの確認用に表示する。
撮影用コンソール2は、図1に示すように、制御部21、記憶部22、操作部23、表示部24、通信部25を備えて構成され、各部はバス26により接続されている。
図1では、撮影装置1aを制御する撮影用コンソール2についてのみ詳細な構成を示しているが、撮影装置1b,1cを制御する撮影用コンソール2についても同様の構成を備えている。
撮影用コンソール2は、放射線照射条件や画像読取条件を撮影装置1に出力して撮影装置1による放射線撮影及び放射線画像の読み取り動作を制御するとともに、撮影装置1により取得された動態画像を撮影技師等の撮影実施者によるポジショニングの確認や診断に適した画像であるか否かの確認用に表示する。
撮影用コンソール2は、図1に示すように、制御部21、記憶部22、操作部23、表示部24、通信部25を備えて構成され、各部はバス26により接続されている。
図1では、撮影装置1aを制御する撮影用コンソール2についてのみ詳細な構成を示しているが、撮影装置1b,1cを制御する撮影用コンソール2についても同様の構成を備えている。
制御部21は、CPU(Central Processing Unit)、RAM(Random Access Memory)等により構成される。制御部21のCPUは、操作部23の操作に応じて、記憶部22に記憶されているシステムプログラムや各種処理プログラムを読み出してRAM内に展開し、展開されたプログラムに従って後述する撮影制御処理を始めとする各種処理を実行し、撮影用コンソール2各部の動作や、撮影装置1の放射線照射動作及び読み取り動作を集中制御する。
記憶部22は、不揮発性の半導体メモリーやハードディスク等により構成される。記憶部22は、制御部21で実行される各種プログラムやプログラムにより処理の実行に必要なパラメーター、或いは処理結果等のデータを記憶する。例えば、記憶部22は、撮影制御に関する処理を実行するためのプログラムを記憶している。
各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部21は、当該プログラムコードに従った動作を逐次実行する。
また、記憶部22は、撮影の種別や検査対象部位(すなわち撮影部位。ここでは胸部とする)に対応する撮影条件(放射線照射条件及び画像読取条件)を記憶している。
放射線照射条件は、例えば、X線管電流の値、X線管電圧の値、フィルター種、SID(Source to Image−receptor Distance)、動画撮影時のパルスレート、パルス幅、パルス間隔等である。また、画像読取条件は、例えば、画素サイズ、画像サイズ(マトリックスサイズ)、動画撮影時のフレームレート、フレーム間隔等である。フレームレートは、パルスレートと一致している。
さらに、記憶部22には、図示しないRIS(Radiology Information System)等から送信される撮影オーダー情報が記憶されている。撮影オーダー情報には、患者情報、検査情報(検査ID、検査対象部位(ここでは、胸部)、解析の種類(例えば、換気解析、肺血流解析、最大換気量の測定等)、データ属性(緊急、外来一般、病棟経過観察)等)等が含まれる。
各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部21は、当該プログラムコードに従った動作を逐次実行する。
また、記憶部22は、撮影の種別や検査対象部位(すなわち撮影部位。ここでは胸部とする)に対応する撮影条件(放射線照射条件及び画像読取条件)を記憶している。
放射線照射条件は、例えば、X線管電流の値、X線管電圧の値、フィルター種、SID(Source to Image−receptor Distance)、動画撮影時のパルスレート、パルス幅、パルス間隔等である。また、画像読取条件は、例えば、画素サイズ、画像サイズ(マトリックスサイズ)、動画撮影時のフレームレート、フレーム間隔等である。フレームレートは、パルスレートと一致している。
さらに、記憶部22には、図示しないRIS(Radiology Information System)等から送信される撮影オーダー情報が記憶されている。撮影オーダー情報には、患者情報、検査情報(検査ID、検査対象部位(ここでは、胸部)、解析の種類(例えば、換気解析、肺血流解析、最大換気量の測定等)、データ属性(緊急、外来一般、病棟経過観察)等)等が含まれる。
操作部23は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号を制御部21に出力する。また、操作部23は、表示部24の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部21に出力する。
表示部24は、LCD(Liquid Crystal Display)やCRT(Cathode Ray Tube)等のモニターにより構成され、制御部21から入力される表示信号の指示に従って、操作部23からの入力指示やデータ等を表示する。
また、撮影装置1で撮影された画像を表示させてポジショニングの確認等を行うことができるようにしてもよい。
また、撮影装置1で撮影された画像を表示させてポジショニングの確認等を行うことができるようにしてもよい。
通信部25は、LANアダプターやモデムやTA(Terminal Adapter)等を備え、通信ネットワークNTに接続された各装置との間のデータ送受信を制御する。
[解析装置3の構成]
解析装置3は、異なる撮影装置1で撮影された動画像に対して動態解析を行うことが可能なものである。
解析装置3は、図1に示すように、制御部31、記憶部32、操作部33、表示部34、通信部35を備えて構成され、各部はバス36により接続されている。
本実施形態では、解析装置3は、撮影用コンソール2から動画像を受信し、受信した動画像を解析したり、画像処理を行う等して、解析装置3の表示部34(後述)や外部の図示しない表示装置に解析結果等を表示させたり、再解析を行ったりする画像処理装置であり医師による診断を支援する診断用コンソールである。
本実施形態において、解析装置3は、撮影用コンソール2から送信された胸部の動画像に基づいて、胸部の動態解析を行う。
解析装置3は、異なる撮影装置1で撮影された動画像に対して動態解析を行うことが可能なものである。
解析装置3は、図1に示すように、制御部31、記憶部32、操作部33、表示部34、通信部35を備えて構成され、各部はバス36により接続されている。
本実施形態では、解析装置3は、撮影用コンソール2から動画像を受信し、受信した動画像を解析したり、画像処理を行う等して、解析装置3の表示部34(後述)や外部の図示しない表示装置に解析結果等を表示させたり、再解析を行ったりする画像処理装置であり医師による診断を支援する診断用コンソールである。
本実施形態において、解析装置3は、撮影用コンソール2から送信された胸部の動画像に基づいて、胸部の動態解析を行う。
本実施形態の解析システム100では、解析装置3には、3台の撮影装置1a,1b,1cと接続された撮影用コンソール2からそれぞれ撮影装置1a,1b,1cによって撮影された動画像が送られるようになっている。
以下においては、このうち、基準となる撮影装置(すなわち、撮影装置1a)で撮影された動画像を「基準動画像D」とし、基準となる撮影装置(すなわち、撮影装置1a)とは異なる撮影装置(すなわち、撮影装置1b,1c)で撮影された動画像を「動画像d」とするが、特に両者を区別しない場合には単に動画像という。
以下においては、このうち、基準となる撮影装置(すなわち、撮影装置1a)で撮影された動画像を「基準動画像D」とし、基準となる撮影装置(すなわち、撮影装置1a)とは異なる撮影装置(すなわち、撮影装置1b,1c)で撮影された動画像を「動画像d」とするが、特に両者を区別しない場合には単に動画像という。
なお、図1では、診断用コンソールである解析装置3と上述の撮影用コンソール2とが別体の装置である場合を例示しているが、解析システム100の構成はこれに限定されない。
例えば、解析システム100は、撮影用コンソール2及び診断用コンソールとしての機能を備える兼用装置と撮影装置1とが通信ネットワークNTにより接続することで構成されていてもよい。この場合には、当該兼用装置が解析装置として機能する。
例えば、解析システム100は、撮影用コンソール2及び診断用コンソールとしての機能を備える兼用装置と撮影装置1とが通信ネットワークNTにより接続することで構成されていてもよい。この場合には、当該兼用装置が解析装置として機能する。
制御部31は、CPU(Central Processing Unit)、RAM(Random Access Memory)等により構成される。制御部31のCPUは、操作部33の操作に応じて、記憶部32に記憶されているシステムプログラムや、各種処理プログラムを読み出してRAM内に展開し、展開されたプログラムに従って、診断支援処理を始めとする各種処理を実行し、解析装置3の各部の動作を集中制御する。制御部31は、関心領域設定手段、濃度波形生成手段、修正手段、設定条件入力手段として機能する。
具体的には、制御部31は、プログラムとの協働により動画像の受信処理及び解析処理を実行する。さらに、制御部31は、操作部33等の要求に応じて解析結果を表示部34又は図示しない外部の表示装置に表示させる。
特に、本実施形態では、解析装置3の制御部31は、異なる複数の撮影装置1から受信した動画像に含まれる構造物S(図3参照)の移動量を用いて、動画像又は動画像の解析結果を補正する補正手段として機能する。
具体的には、補正手段としての制御部31は、基準となる撮影装置(すなわち、撮影装置1a)で撮影された動画像(すなわち、基準動画像D)に含まれる構造物Sの移動量(例えば横隔膜の移動量(位置の変化量))と、動画像又は動画像の解析結果に含まれる信号値(図3では濃度変化量)との相関関係を用いて、基準となる撮影装置とは異なる撮影装置(すなわち、1b,1c)で撮影された動画像(すなわち、動画像d)又は動画像の解析結果を補正する。
ここで、動画像に含まれる「構造物S」とは、例えば、肺や横隔膜、心臓等、被写体である人体を構成する構造物Sであって、被写体の生体活動に応じてその位置や外縁が変化する(移動する)ものである。
なお、本実施形態の解析装置3の制御部31による補正処理の詳細については後述する。
具体的には、制御部31は、プログラムとの協働により動画像の受信処理及び解析処理を実行する。さらに、制御部31は、操作部33等の要求に応じて解析結果を表示部34又は図示しない外部の表示装置に表示させる。
特に、本実施形態では、解析装置3の制御部31は、異なる複数の撮影装置1から受信した動画像に含まれる構造物S(図3参照)の移動量を用いて、動画像又は動画像の解析結果を補正する補正手段として機能する。
具体的には、補正手段としての制御部31は、基準となる撮影装置(すなわち、撮影装置1a)で撮影された動画像(すなわち、基準動画像D)に含まれる構造物Sの移動量(例えば横隔膜の移動量(位置の変化量))と、動画像又は動画像の解析結果に含まれる信号値(図3では濃度変化量)との相関関係を用いて、基準となる撮影装置とは異なる撮影装置(すなわち、1b,1c)で撮影された動画像(すなわち、動画像d)又は動画像の解析結果を補正する。
ここで、動画像に含まれる「構造物S」とは、例えば、肺や横隔膜、心臓等、被写体である人体を構成する構造物Sであって、被写体の生体活動に応じてその位置や外縁が変化する(移動する)ものである。
なお、本実施形態の解析装置3の制御部31による補正処理の詳細については後述する。
記憶部32は、不揮発性の半導体メモリーやハードディスク等により構成される。記憶部32は、制御部31で診断支援処理を実行するためのプログラムを始めとする各種プログラムやプログラムによる処理の実行に必要なパラメーター、或いは処理結果等のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部31は、当該プログラムコードに従った動作を逐次実行する。
また、記憶部32には、過去に撮影された動画像が識別ID、患者情報(被写体属性情報。例えば、患者ID、患者(被写体)の氏名、身長、体重、年齢、性別等)、検査情報(例えば、検査ID、検査日、検査対象部位(ここでは、胸部)等)等に対応付けて記憶されている。また、記憶部32には、撮影用コンソール2からの受信を開始した各動画像に係る患者情報や検査情報、ステータス(例えば、受信中、解析処理中、解析終了等の進捗状況)情報を含むリスト情報が記憶される。さらに、記憶部32には、動画像に対応付けてその解析結果が記憶される。
操作部33は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、ユーザーによるキーボードに対するキー操作やマウス操作により入力された指示信号を制御部31に出力する。また、操作部33は、表示部34の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部31に出力する。
表示部34は、LCD(Liquid Crystal Display)やCRT(Cathode Ray Tube)等のモニターにより構成され、制御部31から入力される表示信号の指示に従って、各種表示を行う。
なお、表示部34は、解析結果を表示するための表示装置としての機能を有していてもよいし、解析結果を表示するための表示装置は、解析装置3とは別装置として設けられていてもよい。
解析結果を表示するための外部の表示装置としては、例えば、医師の診察室等に設置されている電子カルテや、読影医の読影室等に設置されている読影用の端末装置、手術室等に用意され放射線技師や手術を執刀する外科医等が参照する端末装置(コンソール)等が想定される。このように表示装置が解析装置3とは別体で設けられている場合、解析装置3から解析結果等が表示装置に送信される。この場合、解析結果を受け取った表示装置の操作部からパラメーターの調整指示等が入力されると、入力された指示内容が解析装置3側に送信され、解析装置3は、当該指示に従って調整されたパラメーターによる再解析等を行う。
解析結果を表示するための外部の表示装置としては、例えば、医師の診察室等に設置されている電子カルテや、読影医の読影室等に設置されている読影用の端末装置、手術室等に用意され放射線技師や手術を執刀する外科医等が参照する端末装置(コンソール)等が想定される。このように表示装置が解析装置3とは別体で設けられている場合、解析装置3から解析結果等が表示装置に送信される。この場合、解析結果を受け取った表示装置の操作部からパラメーターの調整指示等が入力されると、入力された指示内容が解析装置3側に送信され、解析装置3は、当該指示に従って調整されたパラメーターによる再解析等を行う。
通信部35は、LANアダプターやモデムやTA等を備え、通信ネットワークNTに接続された各装置との間のデータ送受信を制御する。
[解析装置における補正処理について]
前述のように、本実施形態の解析装置3における補正手段としての制御部31は、動画像に含まれる構造物Sの移動量と、動画像又は動画像の解析結果に含まれる信号値との相関関係を用いて、動画像又は動画像の解析結果を補正する。
本実施形態において、解析装置3には、基準となる撮影装置1(例えば撮影装置1a)の動画像(すなわち、基準動画像D)と、基準となる撮影装置1とは異なる撮影装置1(例えば撮影装置1b,1c)の動画像(すなわち、動画像d)とが送られるようになっており、制御部31は、基準となる撮影装置1(例えば撮影装置1a)で撮影された動画像(すなわち、基準動画像D)に含まれる構造物Sの移動量と、当該動画像又は動画像の解析結果に含まれる信号値との相関関係を用いて、基準となる撮影装置1とは異なる撮影装置1(例えば撮影装置1b,1c)で撮影された動画像(すなわち、動画像d)又は動画像の解析結果を補正する。
前述のように、解析装置3では、基準となる撮影装置1(例えば撮影装置1a)で得られた、構造物Sの移動量と、動画像又は動画像の解析結果に含まれる信号値との相関関係を用いて、基準となる撮影装置1とは異なる撮影装置1(例えば撮影装置1b,1c)の動画像や動画像の解析結果を補正することが可能となる他に、当該相関関係を用いることにより、構造物Sの移動量に対して、基準となる撮影装置1(例えば撮影装置1a)自身の動画像や動画像の解析結果を補正することも可能である。このことは当該相関関係を取得した際の条件(例えばX線の照射条件や環境因子)を正とした場合に、撮影毎に生じる条件のバラつき(例えばX線の照射条件違いやノイズ影響)を、動画像や動画像の解析結果から排除するのに有効である。
前述のように、本実施形態の解析装置3における補正手段としての制御部31は、動画像に含まれる構造物Sの移動量と、動画像又は動画像の解析結果に含まれる信号値との相関関係を用いて、動画像又は動画像の解析結果を補正する。
本実施形態において、解析装置3には、基準となる撮影装置1(例えば撮影装置1a)の動画像(すなわち、基準動画像D)と、基準となる撮影装置1とは異なる撮影装置1(例えば撮影装置1b,1c)の動画像(すなわち、動画像d)とが送られるようになっており、制御部31は、基準となる撮影装置1(例えば撮影装置1a)で撮影された動画像(すなわち、基準動画像D)に含まれる構造物Sの移動量と、当該動画像又は動画像の解析結果に含まれる信号値との相関関係を用いて、基準となる撮影装置1とは異なる撮影装置1(例えば撮影装置1b,1c)で撮影された動画像(すなわち、動画像d)又は動画像の解析結果を補正する。
前述のように、解析装置3では、基準となる撮影装置1(例えば撮影装置1a)で得られた、構造物Sの移動量と、動画像又は動画像の解析結果に含まれる信号値との相関関係を用いて、基準となる撮影装置1とは異なる撮影装置1(例えば撮影装置1b,1c)の動画像や動画像の解析結果を補正することが可能となる他に、当該相関関係を用いることにより、構造物Sの移動量に対して、基準となる撮影装置1(例えば撮影装置1a)自身の動画像や動画像の解析結果を補正することも可能である。このことは当該相関関係を取得した際の条件(例えばX線の照射条件や環境因子)を正とした場合に、撮影毎に生じる条件のバラつき(例えばX線の照射条件違いやノイズ影響)を、動画像や動画像の解析結果から排除するのに有効である。
ここで、動画像そのものに対して補正が行われる場合、「信号値」は動画像における関心領域R(ROI、解析対象領域、図3参照)内の濃度平均値又は濃度平均値の変化量を意味する。
また、動画像の解析結果に対して補正が行われる場合、「信号値」は動画像の解析結果における関心領域R(ROI、解析対象領域)内の信号平均値又は信号平均値の変化量を意味する。
人体内の構造物S(例えば肺や横隔膜、心臓等)が被写体の生体活動によって移動すると、その周辺の臓器等も、構造物Sとある程度の相関関係をもって変化する。
例えば、構造物Sとしての横隔膜が上下動すると、これに同期する呼吸運動に伴い肺も伸縮するため、関心領域Rが肺野である場合、関心領域R内の濃度平均値等が変化する。
前述の「構造物S」の移動量の観測点は、構造物Sを一体的に捉えた平均的移動量の場合もあれば、構造物Sを複数の領域に分割した上での、個々の領域の移動量の場合もある。また、前述の「関心領域R」が、解析対象領域(例えば肺野)全体を示す場合もあれば、解析対象領域を複数の領域に分割した上での、個々の領域を示す場合もある。
動画像そのもの(オリジナルの動画像)において、こうした関心領域R内の濃度平均値が変化すると、当該動画像に対して解析を行った場合に得られる解析結果の信号値もこれに応じて変化する。
制御部31による補正は、動画像そのもの(オリジナルの動画像)に対して行っても、動画像を解析した後の解析結果に対して行っても同様の効果を得ることができ、本実施形態では、上記両方の場合を含むものとする。
また、動画像の解析結果に対して補正が行われる場合、「信号値」は動画像の解析結果における関心領域R(ROI、解析対象領域)内の信号平均値又は信号平均値の変化量を意味する。
人体内の構造物S(例えば肺や横隔膜、心臓等)が被写体の生体活動によって移動すると、その周辺の臓器等も、構造物Sとある程度の相関関係をもって変化する。
例えば、構造物Sとしての横隔膜が上下動すると、これに同期する呼吸運動に伴い肺も伸縮するため、関心領域Rが肺野である場合、関心領域R内の濃度平均値等が変化する。
前述の「構造物S」の移動量の観測点は、構造物Sを一体的に捉えた平均的移動量の場合もあれば、構造物Sを複数の領域に分割した上での、個々の領域の移動量の場合もある。また、前述の「関心領域R」が、解析対象領域(例えば肺野)全体を示す場合もあれば、解析対象領域を複数の領域に分割した上での、個々の領域を示す場合もある。
動画像そのもの(オリジナルの動画像)において、こうした関心領域R内の濃度平均値が変化すると、当該動画像に対して解析を行った場合に得られる解析結果の信号値もこれに応じて変化する。
制御部31による補正は、動画像そのもの(オリジナルの動画像)に対して行っても、動画像を解析した後の解析結果に対して行っても同様の効果を得ることができ、本実施形態では、上記両方の場合を含むものとする。
異なる複数の撮影装置1で動画像が撮影された場合、各撮影装置1で撮影された動画像における関心領域Rの濃度変化は、撮影装置(モダリティ)1ごとの放射線の入力に対する信号出力の特性(すなわち入出力特性)の違いにより差異が生じてしまう。
このため、そのままでは動画像の解析結果を互いに比較することができない。
このような問題は、異なる種類のモダリティ間のみならず、同じ種類のモダリティ間であっても、設置条件の差異等により入出力特性が異なる可能性があるため、同様に生じ得る。
このため、そのままでは動画像の解析結果を互いに比較することができない。
このような問題は、異なる種類のモダリティ間のみならず、同じ種類のモダリティ間であっても、設置条件の差異等により入出力特性が異なる可能性があるため、同様に生じ得る。
一般的に、複数の撮影装置1a,1b,1cに対応する静止画像や動画像を扱う解析装置や解析システムでは、最終的に表示画面上に表示される画像の見た目を調整することにより診断上の比較が可能である。
これに対して、本実施形態における解析装置3により行われる動態解析では、複数の撮影装置1a,1b,1cから取得した動画像における関心領域R(ROI、解析対象領域)内の時間軸変化量を把握し、比較する必要がある。
しかし、例えば、放射線に対する入出力特性(すなわち入力(照射された放射線の線量x)に対する出力(濃度)の特性)の異なる撮影装置1a,1b,1cから取得した動画像の間では、仮に同一の被写体を同一の撮影条件下で撮影して得た動画像であっても、各動画像における各フレーム画像を時系列的に動画化した場合に、関心領域R(ROI、解析対象領域)内の時間軸変化量を正しく把握することができず、これらを比較しても正しい解析結果を得ることができない。
本実施形態では、補正手段としての制御部31は、ある撮影装置1(本実施形態では、基準となる撮影装置1a)で撮影された動画像に含まれる人体の構造物S(例えば、肺や横隔膜や心臓外縁等)の移動量と、動画像の関心領域Rの濃度平均値等や動画像の解析結果に含まれる信号値との相関関係(対応付け)を取得し、これ用いて、この基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)で撮影された動画像又は動画像の解析結果を補正する。
これにより、異種を含む複数の撮影装置(モダリティ)1から取得した動画像の解析結果を比較できるようにする。
これに対して、本実施形態における解析装置3により行われる動態解析では、複数の撮影装置1a,1b,1cから取得した動画像における関心領域R(ROI、解析対象領域)内の時間軸変化量を把握し、比較する必要がある。
しかし、例えば、放射線に対する入出力特性(すなわち入力(照射された放射線の線量x)に対する出力(濃度)の特性)の異なる撮影装置1a,1b,1cから取得した動画像の間では、仮に同一の被写体を同一の撮影条件下で撮影して得た動画像であっても、各動画像における各フレーム画像を時系列的に動画化した場合に、関心領域R(ROI、解析対象領域)内の時間軸変化量を正しく把握することができず、これらを比較しても正しい解析結果を得ることができない。
本実施形態では、補正手段としての制御部31は、ある撮影装置1(本実施形態では、基準となる撮影装置1a)で撮影された動画像に含まれる人体の構造物S(例えば、肺や横隔膜や心臓外縁等)の移動量と、動画像の関心領域Rの濃度平均値等や動画像の解析結果に含まれる信号値との相関関係(対応付け)を取得し、これ用いて、この基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)で撮影された動画像又は動画像の解析結果を補正する。
これにより、異種を含む複数の撮影装置(モダリティ)1から取得した動画像の解析結果を比較できるようにする。
なお、動画像や動画像の解析結果の補正処理を行う対象は、基準となる動画像(基準動画像D)を撮影する撮影装置1aとは異なる撮影装置1b,1cで撮影された動画像dやその解析結果に限定されず、撮影装置1aによって撮影された動画像(基準動画像D)やその解析結果であってもよいし、動画像d及び基準動画像Dの双方や動画像dの解析結果及び動画像Dの解析結果の双方を補正処理の対象としてもよい。
以下、動画像そのものに対して補正が行われる場合、及び動画像の解析結果に対して補正が行われる場合について詳説する。
以下、動画像そのものに対して補正が行われる場合、及び動画像の解析結果に対して補正が行われる場合について詳説する。
《動画像そのものに対して補正が行われる場合》
動画像そのものに対して補正が行われる場合、補正手段である制御部31は、基準となる撮影装置1(本実施形態では撮影装置1a)の動画像に含まれる構造物S(例えば横隔膜)の移動量と、当該基準となる撮影装置1aの動画像における関心領域R内の濃度値の平均値(これを濃度平均値という。)又は濃度平均値の変化量と、の相関関係を用いて、基準となる撮影装置1aとは異なる撮影装置1(本実施形態では撮影装置1b,1c)の動画像に含まれる構造物Sの移動量と、当該異なる撮影装置1b,1cの動画像における関心領域R内の濃度平均値又は濃度平均値の変化量を推定した推定値を求める。
そして、この推定値と、異なる撮影装置1b,1cの動画像から実際に得られた関心領域R内の濃度平均値又は濃度平均値の変化量と、の差分から補正式を求め、当該補正式を適用して異なる撮影装置1b,1cの動画像の全画素データを補正する。
動画像そのものに対して補正が行われる場合、補正手段である制御部31は、基準となる撮影装置1(本実施形態では撮影装置1a)の動画像に含まれる構造物S(例えば横隔膜)の移動量と、当該基準となる撮影装置1aの動画像における関心領域R内の濃度値の平均値(これを濃度平均値という。)又は濃度平均値の変化量と、の相関関係を用いて、基準となる撮影装置1aとは異なる撮影装置1(本実施形態では撮影装置1b,1c)の動画像に含まれる構造物Sの移動量と、当該異なる撮影装置1b,1cの動画像における関心領域R内の濃度平均値又は濃度平均値の変化量を推定した推定値を求める。
そして、この推定値と、異なる撮影装置1b,1cの動画像から実際に得られた関心領域R内の濃度平均値又は濃度平均値の変化量と、の差分から補正式を求め、当該補正式を適用して異なる撮影装置1b,1cの動画像の全画素データを補正する。
図2は、横軸に入力(X線量)、縦軸に出力(濃度)をとった場合の、撮影装置1a,1b,1cの入出力特性を表したものである。図2中、「C1a」は撮影装置1aで撮影された動画像の入出力特性を示し、「C1b」は撮影装置1bで撮影された動画像の入出力特性を示し、「C1c」は撮影装置1cで撮影された動画像の入出力特性を示している。
例えば、基準となる撮影装置1aで撮影された動画像(基準動画像D)の入出力特性が「C1a」として示すようなカーブで表される場合、撮影装置1b,1cにおいて撮影され、特有のLUTによる階調処理が行われる等した場合、動画像dの入出力特性が、撮影装置1aで撮影された動画像(基準動画像D)の入出力特性「C1a」とは異なる「C1b」や「C1c」として示されるようなカーブを描く場合がある。
例えば、基準となる撮影装置1aで撮影された動画像(基準動画像D)の入出力特性が「C1a」として示すようなカーブで表される場合、撮影装置1b,1cにおいて撮影され、特有のLUTによる階調処理が行われる等した場合、動画像dの入出力特性が、撮影装置1aで撮影された動画像(基準動画像D)の入出力特性「C1a」とは異なる「C1b」や「C1c」として示されるようなカーブを描く場合がある。
例えば図2に示す例では、同じX線量(入力)の範囲(図中のA参照)における濃度(出力)をみた場合、撮影装置1aでは、「C1a」で示すように、動画像(基準動画像D)の濃度(出力)が比較的小さめの値となるが変動範囲が比較的大きく現れている(図中のB参照)。これに対し、撮影装置1b,1cでは、「C1b」「C1c」で示すように、動画像dの濃度(出力)が基準動画像Dの濃度(出力)と比べて大きめの値になるが変動範囲は比較的小さく現れている(図中のB参照)。このように、撮影装置1a,1b,1cごとに動画像dの入出力特性が異なる。
図3は、仮に基準となる撮影装置1aとこれとは異なる撮影装置1(例えば撮影装置1b,1c)とで、被写体Mの胸部を同一の撮影条件下で撮影して動画像を得た場合の各動画を構成する各フレームに含まれる構造物Sの移動量と、各動画像における各フレーム内の濃度平均値の時系列変化を示したものである。
図3では、右側中段に、動画像に含まれる構造物Sの移動量の時系列的変化を示すグラフを示している。中段のグラフでは、横軸に時間tを取り、縦軸に構造物Sの位置をとっている。
図3左側に示すように、本実施形態では、撮影装置1によって撮影された胸部の動画像には、関心領域Rとしての肺野と、構造物Sとして横隔膜とが含まれている。図3中段のグラフでは、基準となる撮影装置1aで撮影された基準動画像Dに含まれる横隔膜の位置(移動量)を実線で示し、基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)で撮影された動画像dに含まれる横隔膜の位置(移動量)を破線で示している。
図3では、右側中段に、動画像に含まれる構造物Sの移動量の時系列的変化を示すグラフを示している。中段のグラフでは、横軸に時間tを取り、縦軸に構造物Sの位置をとっている。
図3左側に示すように、本実施形態では、撮影装置1によって撮影された胸部の動画像には、関心領域Rとしての肺野と、構造物Sとして横隔膜とが含まれている。図3中段のグラフでは、基準となる撮影装置1aで撮影された基準動画像Dに含まれる横隔膜の位置(移動量)を実線で示し、基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)で撮影された動画像dに含まれる横隔膜の位置(移動量)を破線で示している。
また、図3では、中段のグラフを挟んで下側に、基準となる撮影装置1aで撮影された動画像(基準動画像D)における各フレームの関心領域R内の濃度平均値(Ave.)を模式的に時系列的に並べ(図3においてDa〜De)、最下段に各画素の濃度平均値の時系列変化を実線のグラフで示している。なお、このグラフでは、横軸に時間tを取り、縦軸に関心領域Rの濃度の平均値(Ave.)をとっている。
また、図3において中段のグラフを挟んで上側には、基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)で撮影された動画像dにおける各フレームの関心領域R内の濃度平均値(Ave.)を模式的に時系列的に並べ(図3においてda〜de)、最上段に各画素の濃度平均値の時系列変化を実線のグラフで示している。なお、このグラフでは、横軸に時間tを取り、縦軸に関心領域Rの濃度の平均値(Ave.)をとっている。
また、図3において中段のグラフを挟んで上側には、基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)で撮影された動画像dにおける各フレームの関心領域R内の濃度平均値(Ave.)を模式的に時系列的に並べ(図3においてda〜de)、最上段に各画素の濃度平均値の時系列変化を実線のグラフで示している。なお、このグラフでは、横軸に時間tを取り、縦軸に関心領域Rの濃度の平均値(Ave.)をとっている。
図3に示すように、基準となる撮影装置1aによって撮影された動画像の各フレーム内の関心領域R内の濃度平均値(Ave.)と基準となる撮影装置1aによって撮影された動画像に含まれる構造物S(ここでは横隔膜)の移動量との間には、構造物S(ここでは横隔膜)の位置が下がると、関心領域R内の濃度平均値が上がる(すなわち画像として濃く黒くなる)との相関関係が認められる。
これは、例えば呼吸をすることによって横隔膜の位置が下がった状態では肺全体が引き伸ばされて放射線を遮る肺胞壁が薄くなり、放射線の透過率が上がるためである。
これは、例えば呼吸をすることによって横隔膜の位置が下がった状態では肺全体が引き伸ばされて放射線を遮る肺胞壁が薄くなり、放射線の透過率が上がるためである。
動画像そのものを補正する場合、補正手段としての制御部31は、基準となる撮影装置1aによって撮影された動画像(オリジナル動画像)における構造物S(例えば横隔膜)の位置や移動量と、関心領域R(例えば肺野)内の濃度平均値や濃度の変化量との相関関係を用いて、基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)によって撮影された動画像(オリジナル動画像)における構造物S(例えば横隔膜)の位置や移動量から関心領域R(例えば肺野)内の濃度平均値や濃度の変化量の予想値を求める。図3右側の一番上において、この予想値の時系列的な変化を表す予想曲線を破線で示す。
さらに、制御部31は、基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)によって撮影された動画像(オリジナル動画像)における関心領域R(例えば肺野)内の実際の濃度平均値や濃度の変化量の実測値(図3右側の一番上のグラフにおいて実線で示す)を取得し、この実測値と予想値との差分を求めて、この差分から補正式を導出する。
そして、制御部31は、この補正式を基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)によって撮影された動画像(オリジナル動画像)dの全画素データに適用することにより、基準となる撮影装置1aによって撮影された動画像(オリジナル動画像)Dと画像の特性が同等となるように補正を行う。
さらに、制御部31は、基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)によって撮影された動画像(オリジナル動画像)における関心領域R(例えば肺野)内の実際の濃度平均値や濃度の変化量の実測値(図3右側の一番上のグラフにおいて実線で示す)を取得し、この実測値と予想値との差分を求めて、この差分から補正式を導出する。
そして、制御部31は、この補正式を基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)によって撮影された動画像(オリジナル動画像)dの全画素データに適用することにより、基準となる撮影装置1aによって撮影された動画像(オリジナル動画像)Dと画像の特性が同等となるように補正を行う。
このように、撮影装置1(例えば撮影装置1b,1c)から得た動画像dについて、基準となる撮影装置1aで撮影された動画像(基準動画像D)の入出力特性と同等となるように動画像の特性を補正することにより、解析装置3において動態解析する対象として、基準となる撮影装置1aで撮影された動画像(基準動画像D)と撮影装置1b,1cから得た動画像dとが混在する場合でも、これらを適切に比較して動態解析を行うことができる。これにより、診断に資する適正な解析結果を得ることが可能となる。
《動画像の解析結果に対して補正が行われる場合》
動画像の解析結果に対して補正が行われる場合、補正手段である制御部31は、基準となる撮影装置1(本実施形態では撮影装置1a)の動画像に含まれる構造物S(例えば横隔膜)の移動量と、当該基準となる撮影装置1aの動画像の解析結果における関心領域R内の信号値の平均値(これを信号平均値という。)又は信号平均値の変化量と、の相関関係を用いて、基準となる撮影装置1aとは異なる撮影装置1(本実施形態では撮影装置1b,1c)の動画像に含まれる構造物Sの移動量と、当該異なる撮影装置1b,1cの動画像の解析結果における関心領域R内の信号平均値又は信号平均値の変化量を推定した推定値を求める。
そして、この推定値と、異なる撮影装置1b,1cの動画像の実際の解析結果における関心領域R内の信号平均値又は信号平均値の変化量と、の差分から補正式を求め、当該補正式を適用して異なる撮影装置1b,1cの動画像の解析結果を補正する。
動画像の解析結果に対して補正が行われる場合、補正手段である制御部31は、基準となる撮影装置1(本実施形態では撮影装置1a)の動画像に含まれる構造物S(例えば横隔膜)の移動量と、当該基準となる撮影装置1aの動画像の解析結果における関心領域R内の信号値の平均値(これを信号平均値という。)又は信号平均値の変化量と、の相関関係を用いて、基準となる撮影装置1aとは異なる撮影装置1(本実施形態では撮影装置1b,1c)の動画像に含まれる構造物Sの移動量と、当該異なる撮影装置1b,1cの動画像の解析結果における関心領域R内の信号平均値又は信号平均値の変化量を推定した推定値を求める。
そして、この推定値と、異なる撮影装置1b,1cの動画像の実際の解析結果における関心領域R内の信号平均値又は信号平均値の変化量と、の差分から補正式を求め、当該補正式を適用して異なる撮影装置1b,1cの動画像の解析結果を補正する。
図4は、仮に基準となる撮影装置1aとこれとは異なる撮影装置1(例えば撮影装置1b,1c)とで、被写体Mの胸部を同一の撮影条件下で撮影して動画像を得た場合の各動画を構成する各フレームに含まれる構造物Sの移動量と、各動画像における各フレーム内の信号値の平均値(信号平均値)の時系列変化を示したものである。
図4では、中段に、動画像に含まれる構造物Sの移動量の時系列的変化を示すグラフを示している。中段のグラフでは、横軸に時間tを取り、縦軸に構造物S(本実施形態では横隔膜)の位置をとっている。このグラフにおいて、撮影装置1aで撮影された動画像に含まれる構造物Sの位置を実線で示し、撮影装置1aとは異なる撮影装置1b,1cで撮影された動画像に含まれる構造物Sの位置を破線で示している。
図4では、中段に、動画像に含まれる構造物Sの移動量の時系列的変化を示すグラフを示している。中段のグラフでは、横軸に時間tを取り、縦軸に構造物S(本実施形態では横隔膜)の位置をとっている。このグラフにおいて、撮影装置1aで撮影された動画像に含まれる構造物Sの位置を実線で示し、撮影装置1aとは異なる撮影装置1b,1cで撮影された動画像に含まれる構造物Sの位置を破線で示している。
図4において、中段のグラフを挟んで下側には、基準となる撮影装置1aで撮影された動画像(基準動画像D)における各フレームの関心領域R内のある画素(図4中の円内の小さい四角参照)の濃度平均値(Ave.)を模式的に時系列的に並べ(図4においてDa〜De)、さらに、その下側に各動画像におけるフレーム間の濃度平均値の差分信号値の時系列変化を示している。なお、図4では、動態の解析処理として、隣接するフレーム間での濃度平均値の差分信号値を求めてその時系列変化をみる例を示しているが、動態の解析処理としては、隣接するフレーム間での差分処理に限定されず、他の解析処理を採用してもよい。
図4の最下段には、差分信号値の時系列変化を実線のグラフで示している。なお、このグラフでは、横軸に時間tを取り、縦軸にフレーム間での関心領域Rの濃度平均値の差分信号値(Ave.)をとっている。
図4の最下段には、差分信号値の時系列変化を実線のグラフで示している。なお、このグラフでは、横軸に時間tを取り、縦軸にフレーム間での関心領域Rの濃度平均値の差分信号値(Ave.)をとっている。
図4に示すように、基準となる撮影装置1aによって撮影された動画像の各フレームの関心領域R内の濃度平均値の差分信号値(Ave.)と基準となる撮影装置1aによって撮影された動画像に含まれる構造物S(ここでは横隔膜)の移動量との間には、構造物S(ここでは横隔膜)の変位量が大きいと、関心領域R内の濃度平均値の差分信号値が大きくなるとの相関関係が認められる。
また、図4において中段のグラフを挟んで上側には、基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)で撮影された動画像dにおける各フレームの関心領域R内のある画素(図4中の円内の小さい四角参照)の濃度平均値(Ave.)を模式的に時系列的に並べ(図4においてda〜de)、さらに、その上側に各動画像におけるフレーム間の濃度平均値の差分信号値の時系列変化を示している。なお、ここでも、動態の解析処理として、隣接するフレーム間での信号値の差分を求めてその時系列変化をみる例を示しているが、動態の解析処理としては、隣接するフレーム間での差分処理に限定されず、他の解析処理を採用してもよい。
図4の最上段には、差分信号値の時系列変化を実線のグラフで示している。なお、このグラフでは、横軸に時間tを取り、縦軸にフレーム間での関心領域Rの濃度平均値の差分信号値(Ave.)をとっている。
図4の最上段には、差分信号値の時系列変化を実線のグラフで示している。なお、このグラフでは、横軸に時間tを取り、縦軸にフレーム間での関心領域Rの濃度平均値の差分信号値(Ave.)をとっている。
動画像の解析結果を補正する場合、補正手段としての制御部31は、基準となる撮影装置1aによって撮影された動画像に含まれる構造物S(例えば横隔膜)の位置や移動量と、動画像を解析して得た関心領域R(例えば肺野)内の信号平均値又は信号平均値の変化量との相関関係を用いて、基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)によって撮影された動画像に含まれる構造物S(例えば横隔膜)の移動量から、動画像を解析して得られるはずの関心領域R(例えば肺野)内の信号平均値又は信号平均値の変化量の予想値を求める。図4右側の一番上において、この予想値の時系列的な変化を表す予想信号曲線を破線で示す。
さらに、制御部31は、基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)によって撮影された動画像の実際の解析結果における関心領域R内の信号平均値又は信号平均値の変化量の実測値(図4右側の一番上のグラフにおいて実線で示す)を取得し、この実測値と予想値との差分を求めて、この差分から補正式を導出する。
そして、制御部31は、この補正式を基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)によって撮影された動画像dの解析結果に適用することにより、基準となる撮影装置1aによって撮影された動画像Dの解析結果と同等となるように補正を行う。
さらに、制御部31は、基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)によって撮影された動画像の実際の解析結果における関心領域R内の信号平均値又は信号平均値の変化量の実測値(図4右側の一番上のグラフにおいて実線で示す)を取得し、この実測値と予想値との差分を求めて、この差分から補正式を導出する。
そして、制御部31は、この補正式を基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)によって撮影された動画像dの解析結果に適用することにより、基準となる撮影装置1aによって撮影された動画像Dの解析結果と同等となるように補正を行う。
このように、撮影装置1(例えば撮影装置1b,1c)から得た動画像dの解析結果について、基準となる撮影装置1aで撮影された動画像(基準動画像D)の解析結果と同等となるように動画像の解析結果を補正することにより、解析装置3において動態解析する対象として、基準となる撮影装置1aで撮影された動画像(基準動画像D)の解析結果と撮影装置1b,1cから得た動画像dの解析結果とが混在する場合でも、これらを適切に比較して動態解析を行うことができる。これにより、診断に資する適正な解析結果を得ることが可能となる。
《被写体の体格情報に応じた補正》
上記のように、基準となる撮影装置1aで撮影された動画像(基準動画像D)と、基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)から得た動画像dとの特性を同様となるように合せる補正を行う場合、補正手段である制御部31は、被写体Mの体格に関する情報(これを体格情報とする。)に基づいて、構造物Sの移動量と、動画像又は動画像の解析結果に含まれる信号値との相関関係における相関係数を補正してもよい。
上記のように、基準となる撮影装置1aで撮影された動画像(基準動画像D)と、基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)から得た動画像dとの特性を同様となるように合せる補正を行う場合、補正手段である制御部31は、被写体Mの体格に関する情報(これを体格情報とする。)に基づいて、構造物Sの移動量と、動画像又は動画像の解析結果に含まれる信号値との相関関係における相関係数を補正してもよい。
構造物Sの移動量と、動画像又は動画像の解析結果に含まれる信号値との相関関係は、被写体Mの体格によって異なる場合がある。このため、体格を考慮して相関係数を補正することで、よりきめ細かい正確な補正処理を行うことが期待できる。
このような体格情報に基づく補正は、上記動画像そのものに対して補正が行われる場合及び動画像の解析結果に対して補正が行われる場合のいずれの場合においても必要に応じて併せて適用することが可能である。
ここで、被写体Mの体格とは、身長、体重、体厚等身体の大きさに関するものや、筋肉質であるか否かや、骨格、骨密度、体脂肪率等の体組成に関するもの等、放射線の透過具合に影響を及ぼす身体的な特徴を広く含む概念とする。
被写体Mの体格情報は、例えば撮影を行う際のオーダー情報に含まれる患者情報等から抽出することができる。
なお、以下では、体格として、体重が重いか否かに応じて補正を行う場合を例示して説明する。
このような体格情報に基づく補正は、上記動画像そのものに対して補正が行われる場合及び動画像の解析結果に対して補正が行われる場合のいずれの場合においても必要に応じて併せて適用することが可能である。
ここで、被写体Mの体格とは、身長、体重、体厚等身体の大きさに関するものや、筋肉質であるか否かや、骨格、骨密度、体脂肪率等の体組成に関するもの等、放射線の透過具合に影響を及ぼす身体的な特徴を広く含む概念とする。
被写体Mの体格情報は、例えば撮影を行う際のオーダー情報に含まれる患者情報等から抽出することができる。
なお、以下では、体格として、体重が重いか否かに応じて補正を行う場合を例示して説明する。
図5は、動画像そのものを補正する場合を例として横軸に時間tを取り、縦軸にフレーム間での関心領域Rの濃度平均値(Ave.)をとったものである。
図5では、各画素の濃度平均値の時系列変化を実線のグラフで示し、構造物Sの移動量と、動画像の濃度値(濃度平均値)との相関関係を用いて求めた濃度平均値や濃度の変化量の予想値の時系列的な変化を表す予想曲線を破線で示している。
さらに、この相関関係における相関係数を被写体Mの体格情報に基づいて、構造物Sの移動量と、動画像又は動画像の解析結果に含まれる信号値との相関関係における相関係数を補正した後の、濃度平均値や濃度の変化量の予想値(体格補正後予想値)の時系列的な変化を表す体格補正済予想曲線を一点鎖線で示している。
なお、体格情報による相関係数の補正は、ある閾値以上であれば補正を行うというように補正を行うか否かの1段階の判断で行ってもよいし、例えば、体重50kg未満であれば補正無し、50kg以上80kg未満であれば第1の補正値で補正、80kg以上であれば第2の補正値で補正するというように、複数段階に閾値と対応付けた補正値を用意しておき、被写体Mの体格情報に応じて適用する補正値を変えるようにしてもよい。
図5では、各画素の濃度平均値の時系列変化を実線のグラフで示し、構造物Sの移動量と、動画像の濃度値(濃度平均値)との相関関係を用いて求めた濃度平均値や濃度の変化量の予想値の時系列的な変化を表す予想曲線を破線で示している。
さらに、この相関関係における相関係数を被写体Mの体格情報に基づいて、構造物Sの移動量と、動画像又は動画像の解析結果に含まれる信号値との相関関係における相関係数を補正した後の、濃度平均値や濃度の変化量の予想値(体格補正後予想値)の時系列的な変化を表す体格補正済予想曲線を一点鎖線で示している。
なお、体格情報による相関係数の補正は、ある閾値以上であれば補正を行うというように補正を行うか否かの1段階の判断で行ってもよいし、例えば、体重50kg未満であれば補正無し、50kg以上80kg未満であれば第1の補正値で補正、80kg以上であれば第2の補正値で補正するというように、複数段階に閾値と対応付けた補正値を用意しておき、被写体Mの体格情報に応じて適用する補正値を変えるようにしてもよい。
補正手段である制御部31は、体格情報に基づいて適宜相関係数を補正し、補正後の相関係数を用いて体格補正後予想値をもとめると、基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)によって撮影された動画像における関心領域R(例えば肺野)内の実際の濃度平均値や濃度の変化量の実測値(図5において実線で示す)を取得し、この実測値と体格補正後予想値との差分を求めて、この差分から補正式を導出する。
そして、制御部31は、この補正式を基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)によって撮影された動画像dの全画素データに適用することにより、基準となる撮影装置1aによって撮影された動画像Dと画像の特性が同等となるように補正を行う。
そして、制御部31は、この補正式を基準となる撮影装置1aとは異なる撮影装置1(例えば撮影装置1b,1c)によって撮影された動画像dの全画素データに適用することにより、基準となる撮影装置1aによって撮影された動画像Dと画像の特性が同等となるように補正を行う。
このように、被写体Mの体格をも考慮した上で、撮影装置1(例えば撮影装置1b,1c)から得た動画像dについて、基準となる撮影装置1aで撮影された動画像(基準動画像D)の入出力特性と同等となるように動画像を補正することにより、解析装置3において動態解析する対象として、基準となる撮影装置1aで撮影された動画像(基準動画像D)と撮影装置1b,1cから得た動画像dとが混在する場合でも、また、同一の撮影装置で異なる被写体Mの動画像が混在する場合でも、更には、これら異なる撮影装置を跨いで撮影された異なる被写体Mの動画像が混在する場合でも、これら動画像やその解析結果の間で、診断に資する適切な比較を行うことが可能となる。
[効果]
以上のように、本実施形態に係る解析装置3及び解析システム100によれば、異なる複数の撮影装置1から動画像を受信し、受信した動画像や動画像の解析結果を表示する場合に、補正手段である制御部31が、動画像に含まれる構造物Sの移動量を用いて、動画像又は動画像の解析結果を補正するようになっている。
以上のように、本実施形態に係る解析装置3及び解析システム100によれば、異なる複数の撮影装置1から動画像を受信し、受信した動画像や動画像の解析結果を表示する場合に、補正手段である制御部31が、動画像に含まれる構造物Sの移動量を用いて、動画像又は動画像の解析結果を補正するようになっている。
このため、複数の撮影装置1a,1b,1cの製造メーカーが異なっている場合や製造メーカーは同じであっても設定等によって撮影装置1a,1b,1cに個体差が生じている場合、また、撮影装置1a,1b,1cとして放射線画像撮影装置、CT、MRI等の複数の種類の装置が混在している場合等、異なる複数の撮影装置1で得られる動画像の特性が互いに異なる(ばらつきがある)場合であっても、これらの撮影装置1a,1b,1cで得られた動画像や動画像の解析結果を適切に補正することにより、各撮影装置1a,1b,1cで得られた動画像の解析結果同士を的確に比較することが可能となる。
本実施形態では、人体に存する構造物Sの移動量と動画像又は前記動画像の解析結果に含まれる信号値との相関関係に着目し、補正手段である制御部31は、この相関関係を用いて動画像や動画像の解析結果を補正する。
このため、簡易かつ適切に動画像の解析結果同士を的確に比較することが可能な状態とすることができる。
このため、簡易かつ適切に動画像の解析結果同士を的確に比較することが可能な状態とすることができる。
本実施形態では、相関関係が、動画像内の構造物Sの移動量と、動画像における関心領域R内の濃度平均値又は濃度平均値の変化量の関係で求められ、補正対象がオリジナルの動画像そのものである場合に補正することができる。
このように、何ら画像処理等が加えられる前の動画像を補正することで、より適切に動画像の解析結果同士を的確に比較することが可能な状態とすることができる。
このように、何ら画像処理等が加えられる前の動画像を補正することで、より適切に動画像の解析結果同士を的確に比較することが可能な状態とすることができる。
また、本実施形態では、相関関係が、動画像内の構造物Sの移動量と、動画像の解析結果における関心領域R内の信号平均値又は信号平均値の変化量の関係で求められ、補正対象が動画像の解析結果である場合でも補正することができる。
このため、撮影装置1等において、既に何らかの画像処理等が加えられた後でも適切に動画像の解析結果同士を的確に比較することが可能な状態とすることができる。
このため、撮影装置1等において、既に何らかの画像処理等が加えられた後でも適切に動画像の解析結果同士を的確に比較することが可能な状態とすることができる。
また、本実施形態では、補正手段である制御部31は、動画像又は動画像の解析結果を補正する際に、動画像の被写体の体格に関する体格情報に基づいて相関関係における相関係数を補正することもできるようになっている。
このため、被写体Mの体格をも考慮した上で、撮影装置1(例えば撮影装置1b,1c)から得た動画像dについて、基準となる撮影装置1aで撮影された動画像(基準動画像D)の入出力特性と同等となるように動画像や動画像の解析結果を補正することにより、より適切に動画像の解析結果同士を的確に比較することが可能な状態とすることができる。
このため、被写体Mの体格をも考慮した上で、撮影装置1(例えば撮影装置1b,1c)から得た動画像dについて、基準となる撮影装置1aで撮影された動画像(基準動画像D)の入出力特性と同等となるように動画像や動画像の解析結果を補正することにより、より適切に動画像の解析結果同士を的確に比較することが可能な状態とすることができる。
[変形例]
なお、本発明が上記の実施形態に限定されず、本発明の趣旨を逸脱しない限り、適宜変更可能であることは言うまでもない。
なお、本発明が上記の実施形態に限定されず、本発明の趣旨を逸脱しない限り、適宜変更可能であることは言うまでもない。
例えば、本実施形態では、解析システム100内の撮影装置1のうち撮影装置1aが基準動画像Dを撮影する基準となる撮影装置である場合を想定しているが、基準動画像Dを撮影する基準となる撮影装置は、必ずしも解析システム100内に存在する実体のある装置である必要はない。
例えば、「基準となる撮影装置」が、予め実験等により入出力特性等が把握された標準的装置であり、これによって撮影される動画像(基準動画像D)内の構造物の移動量と、関心領域内の濃度変化や動態解析結果の信号値の変化との相関関係が工場出荷値(データ)の形で工場出荷時に予め解析装置内(解析装置3の記憶部32内等)に存在していてもよい。
この場合には、解析システム100内に存在する撮影装置1の全てが「基準となる撮影装置とは異なる撮影装置」となり、これら各撮影装置1によって取得される動画像dについて、解析装置3内に予め記憶されている基準動画像D内の構造物の移動量と、関心領域内の濃度変化や動態解析結果の信号値の変化との相関関係と等しくなるように解析装置3の制御部31(すなわち補正手段)において補正処理を行う。
例えば、「基準となる撮影装置」が、予め実験等により入出力特性等が把握された標準的装置であり、これによって撮影される動画像(基準動画像D)内の構造物の移動量と、関心領域内の濃度変化や動態解析結果の信号値の変化との相関関係が工場出荷値(データ)の形で工場出荷時に予め解析装置内(解析装置3の記憶部32内等)に存在していてもよい。
この場合には、解析システム100内に存在する撮影装置1の全てが「基準となる撮影装置とは異なる撮影装置」となり、これら各撮影装置1によって取得される動画像dについて、解析装置3内に予め記憶されている基準動画像D内の構造物の移動量と、関心領域内の濃度変化や動態解析結果の信号値の変化との相関関係と等しくなるように解析装置3の制御部31(すなわち補正手段)において補正処理を行う。
また、本実施形態では、「基準となる撮影装置」により撮影された動画像に合わせて、「基準となる撮影装置とは異なる撮影装置」により撮影された動画像やこの動画像の解析結果を補正する例を主として説明したが、解析装置3(本実施形態では、補正手段として機能する解析装置3の制御部31)により行われる補正処理は、これに限定されない。
例えば補正手段としての解析装置3の制御部31は、基準となる撮影装置(すなわち、撮影装置1a)で撮影された動画像(すなわち、基準動画像D)やその解析結果を、基準となる撮影装置(すなわち、撮影装置1a)とは異なる撮影装置(すなわち、撮影装置1b,1c)で撮影された動画像(すなわち、動画像d)やその解析結果と比較可能となるように、基準となる撮影装置1aで撮影された基準動画像D又はその解析結果に対して補正を施すようにしてもよい。また、基準動画像Dと動画像dと双方に補正を施して両者の中間点に基準動画像Dやその解析結果及び動画像dやその解析結果を近づけるようにしてもよい。
例えば補正手段としての解析装置3の制御部31は、基準となる撮影装置(すなわち、撮影装置1a)で撮影された動画像(すなわち、基準動画像D)やその解析結果を、基準となる撮影装置(すなわち、撮影装置1a)とは異なる撮影装置(すなわち、撮影装置1b,1c)で撮影された動画像(すなわち、動画像d)やその解析結果と比較可能となるように、基準となる撮影装置1aで撮影された基準動画像D又はその解析結果に対して補正を施すようにしてもよい。また、基準動画像Dと動画像dと双方に補正を施して両者の中間点に基準動画像Dやその解析結果及び動画像dやその解析結果を近づけるようにしてもよい。
また、「基準となる撮影装置」により得られた動画像(基準動画像D)の構造物Sの移動量と信号値との相関関係は、同じ被写体Mに適用される場合に限定されない。
例えば、複数の被写体Mについてサンプルを集めて構造物Sの移動量と信号値との相関関係のモデルを作成しておき、これを様々な被写体Mの撮影の際に共通して適用してもよい。多くのサンプルに基づいて相関関係のモデルを作成することでより普遍的なモデルを実現することができる。
例えば、複数の被写体Mについてサンプルを集めて構造物Sの移動量と信号値との相関関係のモデルを作成しておき、これを様々な被写体Mの撮影の際に共通して適用してもよい。多くのサンプルに基づいて相関関係のモデルを作成することでより普遍的なモデルを実現することができる。
また、本実施形態では、関心領域Rが肺野であり、構造物Sが横隔膜である場合を例示したが、関心領域R及び構造物Sは、構造物Sの位置・移動量と関心領域Rの信号値の変化との間に相関関係が認められるようなものであればよく、肺野と横隔膜に限定されない。
例えば、関心領域Rが心臓周辺の血流である血流解析を行う場合には、心臓周辺の血流との間に相関関係が認められる心臓外縁が構造物Sとなり、制御部31は、心臓外縁の位置・移動量と心臓周辺の血流との相関関係を用いて動画像又は動画像の解析結果を補正する。
例えば、関心領域Rが心臓周辺の血流である血流解析を行う場合には、心臓周辺の血流との間に相関関係が認められる心臓外縁が構造物Sとなり、制御部31は、心臓外縁の位置・移動量と心臓周辺の血流との相関関係を用いて動画像又は動画像の解析結果を補正する。
1 撮影装置
1a 基準となる撮影装置
1b,1c 基準となる撮影装置とは異なる撮影装置
3 解析装置
31 制御部
100 解析システム
D 基準動画像
d 動画像
1a 基準となる撮影装置
1b,1c 基準となる撮影装置とは異なる撮影装置
3 解析装置
31 制御部
100 解析システム
D 基準動画像
d 動画像
Claims (6)
- 異なる複数の撮影装置から動画像を受信し、受信した動画像や動画像の解析結果を表示する解析装置であって、
前記動画像に含まれる構造物の移動量を用いて、前記動画像又は前記動画像の解析結果を補正する補正手段を備えることを特徴とする解析装置。 - 前記補正手段は、前記構造物の移動量と、前記動画像又は前記動画像の解析結果に含まれる信号値との相関関係を用いて、前記動画像又は前記動画像の解析結果を補正することを特徴とする請求項1に記載の解析装置。
- 前記補正手段は、前記動画像又は前記動画像の解析結果を補正する際に、前記動画像の被写体の体格に関する情報に基づいて前記相関関係における相関係数を補正することを特徴とする請求項2に記載の解析装置。
- 前記信号値は、前記動画像における関心領域内の濃度平均値又は濃度平均値の変化量であり、
前記補正手段は、基準となる撮影装置の動画像に含まれる前記構造物の移動量と、当該基準となる撮影装置の動画像における関心領域内の濃度平均値又は濃度平均値の変化量との相関関係を用いて、前記基準となる撮影装置とは異なる撮影装置の動画像に含まれる前記構造物の移動量と、当該異なる撮影装置の動画像における関心領域内の濃度平均値又は濃度平均値の変化量を推定した推定値を求め、この推定値と前記異なる撮影装置の動画像から実際に得られた関心領域内の濃度平均値又は濃度平均値の変化量との差分から補正式を求め、当該補正式を適用して前記異なる撮影装置の動画像の全画素データを補正することを特徴とする請求項2又は請求項3に記載の解析装置。 - 前記信号値は、前記動画像の解析結果における関心領域内の信号平均値又は信号平均値の変化量であり、
前記補正手段は、基準となる撮影装置の動画像に含まれる前記構造物の移動量と、当該基準となる撮影装置の動画像の解析結果における関心領域内の信号平均値又は信号平均値の変化量との相関関係を用いて、前記基準となる撮影装置とは異なる撮影装置の動画像に含まれる前記構造物の移動量と、当該異なる撮影装置の動画像の解析結果における関心領域内の信号平均値又は信号平均値の変化量を推定した推定値を求め、この推定値と前記異なる撮影装置の動画像の実際の解析結果における関心領域内の信号平均値又は信号平均値の変化量との差分から補正式を求め、当該補正式を適用して前記異なる撮影装置の動画像の解析結果を補正することを特徴とする請求項2又は請求項3に記載の解析装置。 - 動画像を撮影可能である基準となる撮影装置と、
動画像を撮影可能であって前記基準となる撮影装置とは異なる撮影装置と、
請求項1から請求項5のいずれか一項に記載の解析装置と、
を備えることを特徴とする解析システム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017180792A JP2019054991A (ja) | 2017-09-21 | 2017-09-21 | 解析装置及び解析システム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017180792A JP2019054991A (ja) | 2017-09-21 | 2017-09-21 | 解析装置及び解析システム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2019054991A true JP2019054991A (ja) | 2019-04-11 |
Family
ID=66106237
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017180792A Withdrawn JP2019054991A (ja) | 2017-09-21 | 2017-09-21 | 解析装置及び解析システム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2019054991A (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021199204A1 (ja) * | 2020-03-30 | 2021-10-07 | 株式会社Peco | 方法および電子カルテ管理システム |
-
2017
- 2017-09-21 JP JP2017180792A patent/JP2019054991A/ja not_active Withdrawn
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021199204A1 (ja) * | 2020-03-30 | 2021-10-07 | 株式会社Peco | 方法および電子カルテ管理システム |
JP7005083B1 (ja) * | 2020-03-30 | 2022-02-14 | 株式会社Peco | 方法およびプログラム、電子カルテ管理システム |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11410312B2 (en) | Dynamic analysis system | |
JP2017200565A (ja) | 動態解析装置及び動態解析システム | |
US20180137634A1 (en) | Dynamic image processing system | |
JP6540807B2 (ja) | 撮影用コンソール | |
US20180018772A1 (en) | Dynamic analysis apparatus | |
JP2019058368A (ja) | 動態解析装置及び動態解析システム | |
CN116757991A (zh) | 图像处理装置、图像处理***、图像处理方法及记录介质 | |
JP6848393B2 (ja) | 動態画像処理装置 | |
JP2021194140A (ja) | 画像処理装置及び画像処理方法 | |
US20190180440A1 (en) | Dynamic image processing device | |
JP2019054991A (ja) | 解析装置及び解析システム | |
JP2020000475A (ja) | 動態画像処理装置及びプログラム | |
JP6962025B2 (ja) | 解析装置及び解析システム | |
JP2019005417A (ja) | 動態画像処理装置及び動態画像処理システム | |
JP6885260B2 (ja) | 動態画像解析装置 | |
US10687772B2 (en) | Dynamic analysis apparatus | |
US11049253B2 (en) | Dynamic analysis device and recording medium | |
JP2019010135A (ja) | 動態解析装置及び動態解析システム | |
JP6930638B2 (ja) | 動態解析装置、動態解析プログラム、動態解析方法及び制御装置 | |
JP7424532B1 (ja) | 放射線画像解析装置及びプログラム | |
JP2018064848A (ja) | 動態解析システム | |
US11461900B2 (en) | Dynamic image analysis system and dynamic image processing apparatus | |
US20230059667A1 (en) | Storage medium and case search apparatus | |
JP2020000807A (ja) | 動態画像解析装置及びプログラム | |
JP2021194139A (ja) | 画像表示装置及びプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200318 |
|
A761 | Written withdrawal of application |
Free format text: JAPANESE INTERMEDIATE CODE: A761 Effective date: 20200728 |