JP5948950B2 - 動態解析システム - Google Patents
動態解析システム Download PDFInfo
- Publication number
- JP5948950B2 JP5948950B2 JP2012038085A JP2012038085A JP5948950B2 JP 5948950 B2 JP5948950 B2 JP 5948950B2 JP 2012038085 A JP2012038085 A JP 2012038085A JP 2012038085 A JP2012038085 A JP 2012038085A JP 5948950 B2 JP5948950 B2 JP 5948950B2
- Authority
- JP
- Japan
- Prior art keywords
- analysis
- image data
- image
- data
- imaging
- 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.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims description 458
- 238000003384 imaging method Methods 0.000 claims description 162
- 238000003860 storage Methods 0.000 claims description 137
- 238000000034 method Methods 0.000 claims description 133
- 230000008569 process Effects 0.000 claims description 113
- 230000005855 radiation Effects 0.000 claims description 89
- 238000001514 detection method Methods 0.000 claims description 42
- 230000001678 irradiating effect Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 description 133
- 210000004072 lung Anatomy 0.000 description 60
- 238000003745 diagnosis Methods 0.000 description 40
- 238000004891 communication Methods 0.000 description 39
- 230000008859 change Effects 0.000 description 33
- 238000009423 ventilation Methods 0.000 description 27
- 238000012937 correction Methods 0.000 description 26
- 230000033001 locomotion Effects 0.000 description 23
- 230000017531 blood circulation Effects 0.000 description 20
- 238000004364 calculation method Methods 0.000 description 19
- 210000000038 chest Anatomy 0.000 description 18
- 238000007726 management method Methods 0.000 description 16
- 230000006870 function Effects 0.000 description 13
- 230000029058 respiratory gaseous exchange Effects 0.000 description 13
- 238000011156 evaluation Methods 0.000 description 12
- 238000013500 data storage Methods 0.000 description 11
- 230000000284 resting effect Effects 0.000 description 10
- 230000007423 decrease Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000003434 inspiratory effect Effects 0.000 description 7
- 238000004321 preservation Methods 0.000 description 7
- 238000012546 transfer Methods 0.000 description 6
- 230000002159 abnormal effect Effects 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 5
- 239000008280 blood Substances 0.000 description 5
- 210000004369 blood Anatomy 0.000 description 5
- 238000010191 image analysis Methods 0.000 description 5
- 230000010349 pulsation Effects 0.000 description 5
- 230000000241 respiratory effect Effects 0.000 description 5
- 230000002861 ventricular Effects 0.000 description 5
- 230000000747 cardiac effect Effects 0.000 description 4
- 238000013399 early diagnosis Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 230000036387 respiratory rate Effects 0.000 description 4
- 239000004065 semiconductor Substances 0.000 description 4
- 238000004904 shortening Methods 0.000 description 4
- 206010047571 Visual impairment Diseases 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 3
- 238000012790 confirmation Methods 0.000 description 3
- 230000002950 deficient Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 238000002438 flame photometric detection Methods 0.000 description 3
- 238000002594 fluoroscopy Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 210000000481 breast Anatomy 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000007599 discharging Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 239000010408 film Substances 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 210000000115 thoracic cavity Anatomy 0.000 description 2
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 description 1
- 208000001871 Tachycardia Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- OJIJEKBXJYRIBZ-UHFFFAOYSA-N cadmium nickel Chemical compound [Ni].[Cd] OJIJEKBXJYRIBZ-UHFFFAOYSA-N 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 238000011976 chest X-ray Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 229910001416 lithium ion Inorganic materials 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 229910052987 metal hydride Inorganic materials 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 229910052759 nickel Inorganic materials 0.000 description 1
- PXHVJJICTQNCMI-UHFFFAOYSA-N nickel Substances [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 description 1
- -1 nickel metal hydride Chemical class 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000002250 progressing effect Effects 0.000 description 1
- 230000002685 pulmonary effect Effects 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000007634 remodeling Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 125000003003 spiro group Chemical group 0.000 description 1
- 230000006794 tachycardia Effects 0.000 description 1
- 208000008203 tachypnea Diseases 0.000 description 1
- 206010043089 tachypnoea Diseases 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 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
- 230000007704 transition Effects 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Processing Or Creating Images (AREA)
Description
特許文献1に記載の技術では、一回の撮影により生成されるRAWデータが15フレームであるとすると、解析処理により更に15フレームの表示用データが追加生成される。即ち、従来の1枚の静止画像による診断と比べて約30倍の量のデータが関与することになり、既存の医療施設において装備されているデータ保存容量では対処に不安が残る。
更に、必要以上の高フレームレートで撮影を行った場合、動態撮影を行って取得したフレーム画像群のうち、例えば、奇数フレームのみ診断に使用すれば充分(従って偶数フレームは未使用)というケースも想定される。
被写体に連続的に放射線を照射可能な放射線源と、2次元状に配置された複数の検出素子を有し前記複数の検出素子のそれぞれにおいて前記放射線源により連続的に照射され前記被写体を透過した放射線を順次検出する放射線検出器と、を有し、前記放射線源と前記放射線検出器を用いて前記被写体の動態を撮影し複数のフレーム画像データを生成する撮影手段と、
前記撮影手段により生成された複数のフレーム画像データに基づいて前記被写体の動態を解析する解析手段と、
を備える動態解析システムであって、
前記複数のフレーム画像データを、前記解析手段による解析に使用する複数の画像データと解析に使用しない複数の画像データに弁別する弁別手段と、
前記解析手段により前記解析に使用する複数の画像データを用いて解析が行われた後、前記解析に使用された複数の画像データを保存する第1の保存手段と、
前記解析手段による解析に使用されなかった未使用の複数の画像データを保存する第2の保存手段と、
を備え、
前記弁別手段は、前記撮影手段により生成された複数のフレーム画像データの中から一定のフレーム間隔で抽出した複数のフレーム画像データを前記解析手段による解析に使用する複数の画像データとして弁別する。
前記第1の保存手段は、更に、前記解析手段により生成された解析結果データと前記解析手段において解析の過程で生成された中間生成データとを前記解析に使用された画像データに対応付けて保存する。
検索条件となる特徴量を指定するための入力手段と、
前記入力手段により指定された検索条件を満たす解析結果データ及びこれに対応する使用画像データを前記第1の保存手段から検索する検索手段と、
を備える。
請求項4に記載の発明は、請求項1〜3の何れか一項に記載の発明において、
前記弁別手段は、前記複数のフレーム画像データのフレームレート又は解像度に基づいて分別を行う。
請求項5に記載の発明は、請求項1〜4の何れか一項に記載の発明において、
前記撮影手段は、前記被写体の動態を1周期以上にわたり撮影する。
被写体に連続的に放射線を照射可能な放射線源と、2次元状に配置された複数の検出素子を有し前記複数の検出素子のそれぞれにおいて前記放射線源により連続的に照射され前記被写体を透過した放射線を順次検出する放射線検出器と、を有し、前記放射線源と前記放射線検出器を用いて前記被写体の動態を撮影し複数のフレーム画像データを生成する撮影手段と、
前記撮影手段により生成された複数のフレーム画像データに基づいて前記被写体の動態を解析する解析手段と、
を備える動態解析システムであって、
前記複数のフレーム画像データを、前記解析手段による解析に使用する複数の画像データと解析に使用しない複数の画像データに弁別する弁別手段と、
前記解析手段により前記解析に使用する複数の画像データを用いて解析が行われた後、前記解析に使用された複数の画像データを保存する第1の保存手段と、
前記解析手段による解析に使用されなかった未使用の複数の画像データを保存する第2の保存手段と、
を備え、
前記弁別手段は、前記複数のフレーム画像データのフレームレート又は解像度に基づいて分別を行う。
前記第1の保存手段は、更に、前記解析手段により生成された解析結果データと前記解析手段において解析の過程で生成された中間生成データとを前記解析に使用された画像データに対応付けて保存する。
請求項8に記載の発明は、請求項7に記載の発明において、
検索条件となる特徴量を指定するための入力手段と、
前記入力手段により指定された検索条件を満たす解析結果データ及びこれに対応する使用画像データを前記第1の保存手段から検索する検索手段と、
を備える。
請求項9に記載の発明は、請求項6〜8の何れか一項に記載の発明において、
前記撮影手段は、前記被写体の動態を1周期以上にわたり撮影する。
動態撮影とは、被写体に対し、X線等の放射線をパルス的に連続照射して複数の画像を取得(即ち、連続撮影)することをいう。動態撮影では、例えば、呼吸運動に伴う肺の膨張及び収縮の形態変化、心臓の拍動等の、周期性(サイクル)を持つ被写体の動態を撮影する。この連続撮影により得られた一連の画像を動態画像と呼ぶ。また、動態画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。
静止画撮影とは、従来のフィルム方式やCR方式と同様に撮影部位の濃度分解能に基づく診断に使用されるもので、被写体に対し、X線等の放射線を1回照射して一枚の静止画像を取得することをいう。
撮影室R1には、例えば、立位撮影用のブッキー装置1と、臥位撮影用のブッキー装置2と、放射線源3aと、クレードル4と、コンソール5と、操作卓6と、アクセスポイントAPと、が備えられている。
撮影室R2には、例えば、立位撮影用のブッキー装置1と、臥位撮影用のブッキー装置2と、放射線源3b、3cと、クレードル4と、コンソール5と、操作卓6と、アクセスポイントAPと、が設けられている。
撮影室R3には、例えば、立位撮影用のブッキー装置1と、臥位撮影用のブッキー装置2と、放射線源3bと、クレードル4と、コンソール5と、操作卓6と、アクセスポイントAPと、が設けられている。
図2に、ブッキー装置1の機能構成例を示す。図2に示すように、ブッキー装置1は、制御部11と、検出器装着部12と、通信I/F13と、駆動部14とを備えて構成されている。
例えば、制御部11は、検出器装着部12にFPD9a又は9bが装着されると、装着されたFPDにコネクター12bを介してFPDID(FPDの識別情報)の送信要求を行い、FPDIDが受信されると、自己の識別番号であるブッキーIDをFPDIDに対応付けて通信I/F13を介してコンソール5に送信する。また、受信されたFPDIDを一時的にRAMに記憶する。
また、例えば、制御部11は、検出器装着部12からFPDが抜き取られると、通信I/F13を介してコンソール5にFPDIDを送信し、当該FPDIDの(撮影管理テーブル521からの)消去要求を行う。
ブッキー装置2は、制御部21と、検出器装着部22と、通信I/F23と、駆動部24とを備えて構成されている。制御部21、検出器装着部22、通信I/F23、駆動部24は、それぞれ上述の制御部11、検出器装着部12、通信I/F13、駆動部14と同様の構成であるので説明を援用する。更に、ブッキー装置2は、被写体を載置するための被写体台26を備える。
コンソール5は、当該FPDIDを取得すると、以後、当該FPDを自己の制御下におき、起動やスリープ遷移等を制御する。
なお、本実施の形態においては、FPDを撮影室内に持ち込んだ際及び持ち出す際にクレードル4に装着することで、FPDの撮影室への侵入及び持ち出しをクレードル4を介してコンソール5で検知できるようになっている。なお、FPDの撮影室への侵入及び持ち出しは、上記のクレードル方式以外も採用可能で、例えば、WO2008/111355号公報に開示されているRFID方式等を用いることができる。
例えば、制御部51は、通信I/F55を介してFPDID及びブッキーIDが受信されると、記憶部52の撮影管理テーブル521(図4参照)の受信されたブッキーIDに対応する領域にFPDIDを書き込む。また、制御部51は、通信I/F55を介してクレードル4からFPDIDが受信されると、記憶部52の撮影管理テーブル521のブッキーIDが対応付けられていない領域にFPDIDを書き込む。また、制御部51は、ブッキー装置1又は2を介して画像データが受信されると、画像受信時刻を撮影管理テーブル521の送信元のブッキー装置のブッキーIDに対応する領域に記憶させる。
また、例えば、制御部51は、後述する撮影・解析処理を実行し、HIS/RIS7から取得した撮影オーダー情報に基づいて、このコンソール5が設置されている撮影室でそのオーダーの撮影が可能であるか否かを判断し、判断結果を表示する。そして、撮影が可能である場合、コンソール5は撮影に用いられる放射線源及び撮影に用いられるFPDを制御して撮影を行わせる。
例えば、記憶部52には、各撮影室における撮影を管理するための撮影管理テーブル521が記憶されている。
図4に、撮影管理テーブル521のデータ格納例を示す。図4に示すように、撮影管理テーブル521には、「ブッキーID」、「管球タイプ」、「FPDID」、「画像受信時刻」等の項目が設けられている。「ブッキーID」、「管球タイプ」の領域には、当該コンソール5が設置されている撮影室に設けられているブッキー装置及び放射線源のタイプの情報が予め設定されている。ブッキーIDに対応付けられている「FPDID」の領域は、そのブッキーIDに装着されているFPDを管理するための領域であり、ブッキー装置からFPDID及びブッキーIDが受信された際に、受信されたFPDIDがブッキーIDに対応付けて格納される。また、ブッキーIDが対応付けられていない「FPDID」の領域は、撮影室に存在するFPDを管理するための領域であり、クレードル4からFPDIDが受信された際に、受信されたFPDIDが格納される。なお、ブッキー装置からFPDが抜き取られてFPDIDの消去要求が受信された場合には、制御部51により消去要求されたFPDIDのうち当該ブッキーIDに対応して格納されているFPDIDは消去される。また、ブッキーIDが対応付けられていない「FPDID」の領域に格納されているFPDIDがクレードル4から受信されると、制御部51によりそのFPDIDのFPDは撮影室から持ち出された(即ち、撮影室に存在しなくなった)と判断され、そのFPDIDは撮影管理テーブル521から消去される。また、「画像受信時刻」には、通信I/F55から画像データが受信された際に、その時刻が格納される。
尚、上記の画素サイズ情報は後述するブロック化(ビニング)処理時に使用される。
なお、表示部54の画面上に、透明電極を格子状に配置した感圧式(抵抗膜圧式)のタッチパネル(図示せず)を形成し、表示部54と入力部53とが一体に構成されるタッチスクリーンとしてもよい。この場合、タッチパネルは、手指やタッチペン等で押下された力点のXY座標を電圧値で検出し、検出された位置信号が操作信号として制御部51に出力されるように構成される。なお、表示部54は、一般的なPC(Personal Computer)に用いられるモニタよりも高精細のものであってもよい。
図6に、FPD9aの機能構成例を示す。図6に示すように、FPD9aは、制御部91、検出部92、記憶部93、コネクター94、バッテリー95、無線通信部96等を備えて構成され、各部はバス97により接続されている。
例えば、制御部91は、コネクター94を介して接続されたブッキー装置1又は2からの要求に応じてFPD9aの識別情報であるFPDIDを記憶部93から読み出して要求元のブッキー装置に送信する。
また、例えば、制御部91は、コンソール5から入力された画像読取条件に基づいて検出部92のスイッチング部を制御して、各放射線検出素子(以下、検出素子)に蓄積された電気信号の読み取りをスイッチングしていき、検出部92に蓄積された電気信号を読み取ることにより、画像データ(静止画像又はフレーム画像)を生成する。そして、制御部91は、生成した画像データを、順次コネクター94及びブッキー装置1又は2を介してコンソール5に出力する。尚、撮影により取得された各フレーム画像は、一旦FPD9aの記憶部93に記憶され、全撮影終了後に、纏めてFPD9aからコンソール5に出力することとしても良い。
なお、FPD9aは、ブッキーに装填されない単体使用時には、バッテリー駆動及び無線通信する構成であるが、動態撮影の場合には、登録第4,561,730号公報に開示されているように、ブッキーを介した外部電力供給及び有線通信する構成に変更することが好ましい。これは、静止画撮影に比べ、データ転送容量や転送時間が圧倒的に増えるので、一のフレーム画像の転送中に他のフレーム画像の撮影(読取り)に対しノイズを与えないため、且つ、転送時間自体を短くするためであり、さらに、一連の撮影途中のバッテリー切れを防止するためである。
なお、生成された静止画像又はフレーム画像を構成する各画素は、検出部92の各検出素子のそれぞれから出力された信号値(ここでは、濃度値と呼称する)を示す。
また、単体での使用は勿論のこと、ブッキー装置に装填しても使用可能で、ブッキー装置に装填時には、コネクター接続により、バッテリー/無線方式から、有線/電力供給方式に切り替えることができる。従い、複数の患者を連続的に静止画撮影する場合に於いても、バッテリー切れを気にする必要がなくなる。
PACS10のサーバー装置は、第1の保存手段として、解析に使用された使用画像データや解析結果データを撮影オーダー情報に対応付けて記憶する。
具体的に、PACS10のサーバー装置は、解析結果データ(画像として生成されたもの)、解析に使用された使用画像データ及び解析の過程で生成された中間生成データを格納する領域、及び格納されたデータを管理するための情報管理テーブル、解析項目毎の解析結果テーブル等を有している。
情報管理テーブルは、サーバー装置に保存されている画像データを管理するためのテーブルであり、撮影オーダー情報に含まれる「撮影ID」、「撮影日時」、「患者情報」、「撮影部位」、「撮影方向」、「***(立位、臥位)」、「撮影方法」、「解析項目」のほか、「画像データの保存先パス」(使用画像データの保存先のパス、中間生成データの保存先パス、解析結果データ(画像として生成されたもの)の保存先パス)、「識別ID」、「解析ID」等の項目の格納領域を有する。解析IDは、保存時に付与されるユニークなIDである。
解析結果テーブルは、動態画像の解析結果(算出された特徴量等)の値を解析項目毎に格納するテーブルであり、例えば、「解析ID」、「解析日時」、「結果値1」〜「結果値n」の項目の格納領域を有している。「結果値1」〜「結果値n」は、各解析項目で算出される特徴量等を表しており、解析項目毎に異なる。例えば、解析項目が後述の「(3)換気−気流速度」である場合、解析結果テーブルの結果値に係る項目は、「平均値」、「標準偏差」、「下位5%」、「上位5%」となる。これにより、解析結果(算出された特徴量)の結果値をキーとして、解析IDを介してRAWデータ、中間生成データ、解析結果の画像データを検索することができるので、同じような結果値の症例に係るデータを読み出して医師が診断に使用することが可能となっている。例えば、PACS10の読影端末(入力手段)から特徴量の種類及び値の範囲が検索条件として指定されると、PACS10のサーバー装置は、指定された検索条件で検索を行い、検索条件を満たす解析結果データを読影端末に出力する検索手段として機能する。
例えば、被検者の再診までのサイクルが2ケ月であるとすると、2ヶ月間に発生するデータを保存可能な容量をもち、当該データ満杯になると、古いデータから順次削除して、新規データを上書き保存する。
図7Aに、動態解析システム100において実行される撮影・解析処理の流れを示す。図7Aのコンソール5側の処理は、コンソール5の制御部51と記憶部52に記憶されているプログラムとの協働により実行される。解析サーバー8側の処理は、解析サーバー8の制御部と記憶部に記憶されている解析プログラムとの協働により実行される。
なお、この場合、一旦、撮影室内でFPDの入替操作を行って、再度コンソールまで戻り、再度ステップS1の処理から行わせ、確実性を高めるフローとしても良い。
図7Bに、ステップS4において実行される撮影制御処理のフローチャートを示す。図7Bに示す撮影制御処理は、制御部51と記憶部52に記憶されているプログラムとの協働により実行される。
まず、指定された撮影オーダー情報の撮影が可能な放射線源及びFPDが起動され、使用されるブッキー装置に応じて放射線源の向き及び位置が調整される。撮影技師により被写体に応じてFPDやブッキー装置の位置等が調整されると、それに応じて放射線源の向き及び位置が調整される(ステップS101)。また、記憶部52から撮影する部位や動態撮影か静止画撮影かの別に応じた放射線照射条件及び画像読取条件が読み出され、放射線源に放射線照射条件が設定されるとともに、ブッキー装置を介してFPDに画像読取条件が設定される(ステップS102)。動態撮影の結果を用いて解析を行う場合は、診断に使用し得る解析精度を確保するため、フレームレートが3.75枚/秒以上に設定される。ここで、操作者は、肺野の動態撮影の場合、安静呼吸の動態を撮影するために被検者に楽にするように指示し、安静呼吸を続けるよう促す。撮影準備が整った時点で、技師は前室へ移動し、操作卓6を操作して放射線照射指示を入力する。
動態撮影の場合は、ステップS102で設定されたパルス間隔で放射線源3aにより放射線が照射され、ステップS102で設定されたフレームレートでFPD9aによりフレーム画像が取得される。予め定められたフレーム画像数の撮影が終了すると、制御部51により放射線源3a及びFPD9aに撮影終了の指示が出力され、撮影動作が停止される。撮影されるフレーム画像数は、少なくとも動態1サイクルが撮影できる枚数である。撮影により取得されたフレーム画像はFPD9aからブッキー装置を介して順次コンソール5に入力される。順次送信される各フレーム画像を、補正処理(オフセット・ゲイン・欠陥補正等)や画像処理(諧調処理)を施すことなく、コンソール5の表示部54に表示可能に単純間引き処理を行い順次切替表示するだけでも、技師等の操作者は被検者の姿勢がくずれた(体動が発生した)か否かや撮影対象となる動態が1周期以上含まれているか等をほぼリアルタイムで視認可能となり好ましい。この視認結果に基づき、直ちに再撮影にとりかかることができる。
なお、撮影オーダー情報により解析が指示されていない場合には、オフセット補正用のダーク画像を読取り、コンソール5に入力することとしても良い。
静止画撮影の場合は、ステップS102で設定された条件で被写体の1枚の静止画像及びオフセット補正用の1または数枚のダーク画像が撮影される。撮影により取得された静止画像及びダーク画像は、FPDからブッキー装置を介してコンソール5に入力される。
本件発明者等は、動態解析においては、静止画のような個々の画素の絶対的出力値はあまり重要ではなく、個々の画素におけるフレーム間の相対的出力値(変動成分)に基づく特徴量の算出が基本となっているため、上記の補正処理の一部或いは全部を省略しても、補正処理を行った場合と略同等の解析結果を得ることが可能である知見を得た。従い、解析結果を得るまでの時間を短縮するため、一部或いは全部の補正処理を省略することもできる。
入力部53の所定の操作により撮影OKを示す判断結果が入力されると(ステップS110;YES)、処理は図7AのステップS5に移行する。
使用画像データと未使用画像データの弁別は、例えば、受信されたフレーム画像群のフレームレート、解像度(画素数)等に基づいて行われる。例えば、一連のフレーム画像群が15fpsで動態1周期分(ここでは動態周期3秒間とする)の撮影により得られた画像データである場合、15×3=45フレームとなる。しかし、呼吸の換気機能についての解析には、1周期につき均等な時間間隔で撮影された20フレーム程度で充分であるので、例えば、偶数(奇数)フレームのみ使用し、奇数(偶数)フレームは未使用に弁別される。また、23fpsで撮影された場合、23×3=69フレームとなるが、この場合には、1番目、4番目、7番目、・・・のように3フレーム毎に使用画像データが抽出(弁別)され、残りのフレーム画像は未使用画像データに弁別される。
解析処理の内容は、撮影対象部位毎に異なる。ここでは、肺野の解析を行う場合を例にとり説明する。
フレーム間差分画像は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のローパスフィルタ処理→フレーム間差分処理→ノイズ除去処理
ビニング処理は、上述のように、各フレーム画像において、画像領域を所定サイズの画素ブロック単位の小領域に分割し、小領域毎に領域内画素の信号値の代表値、例えば、平均値を算出する(平均化する)処理である。なお、代表値としては、平均値に限らず、中央値、平均値、最頻値としてもよい。画素ブロックのサイズは、解析対象となる部位、及び/又は解析により算出される特徴量に応じたものとすることが解析精度向上の点から好ましい。
時間軸方向のローパスフィルタ処理は、換気による信号値の時間変化を抽出するための処理であり、例えば、カットオフ周波数0.5Hzでフィルタリングする。
フレーム間差分処理は、一連のフレーム画像の同じ画素位置の小領域(FPDの同じ位置の検出素子から出力された領域)を互いに対応付け、各小領域毎に、隣接するフレーム画像間で信号値の差分値を算出し、フレーム間差分画像を作成する処理である。
なお、フレーム間差分画像の静止画像を作成する場合は、肺野全体の濃度変化又は横隔膜の位置の変化を解析することにより一連のフレーム画像における吸気期間と呼気期間を算出し、小領域毎に吸気期間については正のフレーム間差分値の絶対値を積算し、呼気期間については負のフレーム間差分値の絶対値を積算した画像を作成する。
換気−波形描画は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のローパスフィルタ処理→波形描画
ビニング処理及び時間軸方向のローパスフィルタ処理は上述のとおりである(以下同じ)。
波形描画処理は、一連のフレーム画像の同じ画素位置の小領域(FPDの同一位置の検出素子から出力された画素ブロックの領域)を互いに対応付け、各小領域毎に、横軸を撮影開始からの経過時間、縦軸を画素の平均信号値とした座標平面を作成して、各フレーム画像の撮影開始からの経過時間とその小領域について算出された平均信号値が交わる点をプロットすることより、換気量を示す信号値の時間変化を示す波形を描画する処理である。
気流速度は、各小領域の肺のやわらかさ(肺コンプライアンス)を示す特徴量である。気流速度は、一連のフレーム画像に、以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のローパスフィルタ処理→フレーム間差分処理→フレーム間差分値の代表値(最大値もしくは平均値)を算出
フレーム間差分値)の最大値をそれぞれ算出し、その比(最大流速比)の分布を示すヒス
トグラムを作成するとともに、何れか1つのフレーム画像上の各小領域を最大流速比に応じた輝度もしくは色で示す画像を作成し、両者を並べて解析結果として提供してもよい。この解析を最大流速比のヒストグラム解析と呼ぶ。
最大流速比のヒストグラム解析では、図11A〜図11C等に示すように、各小領域の吸気気流速度の最大値(絶対値)と呼気気流速度の最大値(絶対値)との比の値がヒストグラム表示されるとともに、COPDであるか否かの指標となる肺野全体での平均値や標準偏差が表示された画像が生成される。また、併せて静止画像上の各小領域を比の値に応じた輝度もしくは色で示すことで、異常個所の分布を医師が容易に把握できるような診断情報提供する。
換気量の振幅は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のローパスフィルタ処理→一連のフレーム画像の同じ画素位置の小領域(FPDの同一位置の検出素子から出力された画素ブロックの領域)を互いに対応付け、各小領域毎に、呼吸1サイクル中の最大信号値(極大値)−最小信号値(極小値)の算出
吸気遅延時間は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のローパスフィルタ処理→一連のフレーム画像の同じ画素位置の小領域(FPDの同一位置の検出素子から出力された画素ブロックの領域)を互いに対応付け、肺野全体の濃度変化又は横隔膜位置の変化を解析することにより安静呼気位のフレーム画像を抽出し、各小領域毎に、安静呼気位のフレーム画像から、吸気時において安静呼気位の信号値との差が所定の閾値以上となるまでの時間の算出
吸気時間、呼気時間は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のローパスフィルタ処理→一連のフレーム画像の同じ画素位置の小領域(FPDの同一位置の検出素子から出力された画素ブロックの領域)を互いに対応付け、各小領域毎に、呼吸1サイクル中の最大信号値(極大値)、最小信号値(極小値)の算出→最大信号値から最小信号値までの時間を呼気時間、最小信号値から最大信号値までの時間を呼気時間として算出
フレーム間差分画像は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のハイパスフィルタ処理→フレーム間差分処理→ノイズ除去
時間軸方向のハイパスフィルタ処理は、血流による信号値の時間変化を抽出するための処理であり、例えば、カットオフ周波数0.7Hzでフィルタリングする。その他は、上述の(1)換気−フレーム間差分画像で説明した処理と同様である。
血流−波形描画は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のハイパスフィルタ処理→各小領域毎の波形描画
ビニング処理及び時間軸方向のハイパスフィルタ処理は上述のとおりである(以下同じ)。
波形描画処理は、上述の(2)換気−波形描画で説明した処理と同様の処理である。
血流量の振幅は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のハイパスフィルタ処理→一連のフレーム画像の同じ画素位置の小領域(FPDの同一位置の検出素子から出力された画素ブロックの領域)を互いに対応付け、各小領域毎に、心拍1サイクル中の最大信号値(極大値)−最小信号値(極小値)の算出
吸気遅延時間は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のハイパスフィルタ処理→一連のフレーム画像の同じ画素位置の小領域(FPDの同一位置の検出素子から出力された画素ブロックの領域)を互いに対応付け、心室領域の濃度変化又は心壁位置の変化を解析することにより心室拡張期の終わりに相当するフレーム画像を抽出し、各小領域毎に、心室拡張期の終わりに相当するフレーム画像から、心室収縮期において心室拡張期の終わりの信号値との差が所定の閾値以上となるまでの時間の算出
横隔膜移動量は、撮影された各フレーム画像に、以下の処理を施すことによって算出される。
画像解析により各フレーム画像から横隔膜の位置を抽出→各フレーム画像の横隔膜の位置を追跡し、移動量を算出
(12)胸郭移動量解析
胸郭移動量は、撮影された各フレーム画像に、以下の処理を施すことによって算出される。
画像解析により各フレーム画像から上部胸郭(上肋骨(第2〜第6肋骨))、下部胸郭(下肋骨(第7〜第10肋骨))の位置を抽出→各フレーム画像の上胸郭、下胸郭の位置を追跡し、移動量を算出する。
(13)呼吸数、呼吸周期
呼吸数、呼吸周期は、撮影された各フレーム画像に、以下の処理を施すことによって算出される。
各フレーム画像から画像解析により求めた横隔膜位置(肺尖から横隔膜までの距離)の変化若しくはローパスフィルタ処理後の肺野全体の信号変化(信号値(平均信号値)の極大値から次の極小値までの時間間隔)より呼吸周期を求め、呼吸周期の逆数から単位時間当たりの呼吸数を算出する。
(14)心拍数、心周期
心拍数、心周期は、撮影された各フレーム画像に、以下の処理を施すことによって算出される。
各フレーム画像から画像解析により求めた心壁位置の変化若しくはハイパスフィルタ処理後の肺野全体の信号変化(信号値(平均信号値)の極大値から次の極小値までの時間間隔)より心周期を求め、心周期の逆数から単位時間当たりの呼吸数を算出する。
(15)スパイロ検査相当の値の算出
横隔膜位置の時間変化波形を算出し、FEV1.0%(一秒率)相当の値を算出する。
胸郭と横隔膜位置の変化から肺野面積の変化量を算出し、別途測定した胸厚の変化量を乗ずることで、VC(肺活量)相当の値を算出する。
コンソール5においては、ネットワーク通信部56により解析結果データ、使用画像データ、及び中間生成データが受信されると、受信された解析結果データ、使用画像データ及び中間生成データが撮影オーダー情報と対応付けられてPACS10に送信される(ステップS11)。
そして、撮影・解析処理は終了する。なお、特徴量の解析に使用した各フレーム画像データに関しては、従来の静止画のような濃度諧調ベースの病変部読影には使用できないものであるので、保存データ容量削減の観点から、保存せず、算出された特徴量データのみを保存することが好ましい。また、内部保存する使用画像データに於いても、任意の一フレーム画像を単発で従来の静止画同様の観察に提供することは逆に誤診を招く恐れがあり、誤診防止の観点からは、例えば、最大肺野面積フレーム画像(安静吸気位画像)と最小肺野面積フレーム画像(安静呼気位画像)の比較等の複数のフレーム画像の相対比較に用いること以外には、単発フレームの使用を禁止することが好ましい。
また、解析サーバー8にて、解析結果データ、使用画像データ及び中間生成データと、コンソール5から受信した撮影オーダー情報とを対応付け、これらのデータを解析サーバー8からPACS10に送信する構成としても良い。
このワーピング処理を行うには、一のフレーム画像を複数の小領域に分割し、当該一のフレーム画像の各小領域に描画された構造物の部分と同一の部分を描画した小領域を各フレーム画像毎に抽出しなければならない。ワーピング処理では、一般的には肺野内構造物による空間的な濃度変化をもとに位置合わせを行うため、各フレーム画像にわたり当該構造物の濃度が忠実に(均一に)再現される必要があり、従って、検出器の各画素の出力バラツキ等は極力抑える必要があり(従い、オフセット補正処理、ゲイン補正処理、欠陥画素補正処理、ラグ補正処理等の各種補正処理を行って、バラツキを補正する必要がある)、補正処理に時間を要し、更に、高精度のワーピング処理を行うには、それに従い分解能の細かい画像を必要とするので画素サイズの小さな検出器が必要となり、各フレーム画像毎のデータ容量が増え、処理対象となる全体データ容量が大幅増加となる。そのため、大容量メモリや高速処理CPU等のハードウエアが必要であり、かつ処理時間も必要となっていた。
まず、体厚方向(z方向;側面)における信号値の変化ついて説明する。
図8は、安静呼気位である時刻T1における肺野の体厚方向(z方向)を模式的に示す図、時刻T1から吸気していくことにより安静吸気位となった時刻T2における肺野の体厚方向を模式的に示す図、及びFPDの検出素子位置(体軸方向(y方向))を模式的に示す図である。図8においては、吸気することにより肺胞a及び肺胞bのy方向の位置は下方へ移動し、時刻T1における肺胞bのy方向の位置と時刻T2における肺胞aのy方向の位置が一致した例を示している。
図9に示すように、通常、吸気においては、肺胞は左肺野の場合左下方向に、右肺野の場合右下方向に移動する。これを鉛直方向(y方向)への移動と水平方向(x方向)への移動に分解する。y方向への肺胞の移動に対するワーピング処理は前述のとおりである。以下に、安静換気時の、肺胞のx方向への移動について説明する。
また、従来のシステムにおいて、原画像から間引き画像が作成されていたが、これは被写体の関心領域が解析可能な範囲内に収まっているかの確認や、検出アルゴリズムに適合する濃度範囲に検出対象画像の濃度範囲を収めるための階調処理条件の算出等、検出処理の前準備段階で使用されるだけであり、CADによる検出処理においては使用されることなく、廃棄されるものであった。
ここでは、換気量のフレーム間差分画像、最大流速比のヒストグラム解析、血流量のフレーム間差分画像を対象項目とした場合を例として示している。なお、各フレーム画像は、照射線量及びフレームレート等の放射線照射条件及び画像読取条件を一定として撮影された画像を用いている。図10では、フレームレート7.5枚/秒、10秒間での入射表面線量は0.2mGy、FPD9aの画素サイズは200μmである。なお、図10(及び後述する図11A及び図11B)において、評価結果△以上は診断に使用できる精度で特徴量の算出が可能であることを示し、△→○→◎の順に解析精度が低い→高いとなっていることを示す。また、処理時間×は処理時間が長く、あまり実用化には向かないことを示し、△以上は処理時間が実用上許容できる範囲であることを示す。なお、この処理時間の場合、Xレベルであっても、全画素データを使う場合にくらべると、かなりの処理時間短縮は可能にはなっているので、使用不可というわけでは無い。
図11A〜図11Cに、ブロックサイズをそれぞれ2mm角、5mm角、10mm角としたときの最大流速比のヒストグラム解析の解析結果の一例を示す。
なお、ブロックサイズを0.5mm角、1mm角とした場合は、図11Aに示すヒストグラムとほとんど変化はみられなかった。
図11Aに示すブロックサイズ2mm角の解析結果と図11Bに示すブロックサイズ5mm角の解析結果を比較すると、若干ヒストグラムの形状に変化が見られる。図11Aに示すブロックサイズ2mm角の解析結果と図11Cに示すブロックサイズ10mm角の解析結果を比較すると、ヒストグラムの形状に大幅な変化が見られ、解析精度が大幅に低下することがわかる。よって、解析精度の点からいえば、ブロックサイズ5mm角以下、特に、2mm角以下であることが必要である。一方、処理時間の点からすると、上述の図10に示すように、ブロックサイズが小さいほど処理時間は長くなる。
よって、解析精度の維持及び処理時間短縮化の観点から、肺野の換気の解析ではブロックサイズは2mm角から5mm角程度とすることが好ましい。
以上の観点からも、肺野の換気に関する解析処理を行う場合のブロックサイズは2mm角〜5mm角程度が好ましい。
尚、従来から行われている間接撮影(縮小画像に基づく診断)による1次スクリーニング、異常疑い者に対しては直接撮影(ライフサイズ画像に基づく診断)による2次精査というフローをとる所謂集団検診に胸部動態解析を流用する場合、左右の肺野領域を上・中・下にそれぞれ3分割するような比較的大きなブロックサイズとすることで、各ブロックサイズ毎の解析結果を1次スクリーニングの代替とし、疑い被検者に対しては、解析用のブロックサイズを再編して小ブロック化し、当該小ブロック毎に局所的な、言い換えると精密な解析を行い2次精査の代替とすることも可能である。このフローの特徴は、従来のフローと異なり再撮影が不要で被検者の負担低減がはかれるとともに早期診断に繋げることができる。
図12A〜図12Bは、FPD9aのフレームレートを2枚/秒〜20枚/秒と変化させたときの上記解析処理の代表的な項目の解析結果、処理時間、間引きによる影響、患者被曝の程度、の評価結果を示す図である。図12Aは、各フレーム画像における照射線量を一定にした場合の評価結果を示す。図12Bは、一回の動態撮影のトータルの照射線量を一定(撮影10秒間での入射照射線量0.2mGy相当)にした場合の評価結果を示す。図12A、図12Bとも画素サイズは200μm、ブロックサイズは2mm角である。
図13A〜図13Eに、動態撮影においてトータルの照射線量を一定とし、フレームレートを2枚/秒〜30枚/秒と変化させたときの最大流速比のヒストグラム解析の解析結果の一例を示す。
図13Aに示すように、フレームレート2枚/秒とした場合、フレームレート3.75枚/秒とした場合に比べ、平均値は1.15→1.32(15%増)、分散値は0.38→0.52(37%増)と解析結果に差異が生じている。また、異常値となる小領域が増えており(図13A〜図13Eの肺野領域の中で静止画上に最大流速比に応じた輝度が重畳されていない領域)、解析精度は低下している。ここで、吸気、呼気のそれぞれにおいてフレーム間差分値の最大値を算出する際、ノイズの影響を抑制する目的で、フレーム間差分値の最大値が所定閾値より小さい場合は異常値としてその小領域は解析対象から除外している。より具体的には、吸気時は信号値が増加するため、フレーム間差分値の最大値が正の所定閾値以上である場合を正常値、呼気時は信号値が減少するため、フレーム間差分値の最大値が負の所定閾値以下である場合を正常値、それ以外を異常値と判定している。
例えば、SID=200cm(被写体厚を20cmと想定)、フレームレート7.5枚/秒における胸部動態撮影の撮影条件の一例としては、管電圧100kV、管電流50mA、パルス幅2ms、付加フィルタAl0.5mm+Cu0.1mmが挙げられる。この
とき、1フレーム画像あたりのX線パルスの照射線量は0.1mAs(=50mA×0.002s)であり、トータル線量が一定となるようにフレームレートに応じて1フレーム画像あたりの照射線量を変えた場合、フレームレート15枚/秒における1フレーム画像あたりのX線パルスの照射線量は0.05mAs、フレームレート30枚/秒における1フレーム画像あたりのX線パルスの照射線量は0.025mAsとなる。しかし、X線パルス幅を1msに制御することは難しく、また、管電流を下げようとすると、X線管のアノード・カソード間に存在するコンデンサに蓄えられた電荷を放電する時定数が長くなり、その結果、管電圧の下降が穏やかになるため、上記電荷の放電を急峻に行うための回路が必要となり、装置コストが増大する。従って、X線パルスの制御の観点では、低フレームレートの方が有利となる。
また、一般的に、FPDにおける残像(ラグ)は、時間(撮影間隔)に対して指数関数的に減少するため、トータル線量が一定となるようにフレームレートに応じて1フレーム画像当たりの照射線量を変えた場合でも、即ち、フレームレートに反比例させて線量を減少させた場合でも、残像の視認性の観点では、フレームレートが小さいほうが有利となる。
従って、換気の解析では、フレームレートを3.75枚/秒〜7.5枚/秒とすることが好ましい。
フレームレート3.75枚/秒とした場合、フレーム画像間隔が離れすぎているため、場合によっては心臓からの拍出タイミングを捉えられない場合がある。フレームレート7.5枚/秒、15枚/秒とした場合、図14B、図14Cに示すように、平均的な脈拍数(安静時50回〜100回/分。60回/分とすると心拍周期は1.0秒。)の成人の画像であれば血液の拍出タイミングを捉えることができる。ただし、レアケースな頻脈患者(単位時間当たりの脈拍数が多い患者。100回〜120回/分。120回/分とすると心拍周期は0.5秒。)の場合、15枚/秒以上が好ましい。フレームレート30枚/秒とした場合は、1フレーム画像当たりの線量が少ないためにノイズが多く、フレーム間差分画像の画質が低下している。
従って、血流の解析では、フレームレート7.5枚/秒〜15枚/秒とすることが好ましい。
また、ポジショニング確認用にビニング処理及び単純間引き処理された画像を解析用に再利用することができるので、処理時間を大幅に短縮化することができる。
これは、既存の動画撮影対応、例えば透視用のFPD等に於いては、FPD自体が、例えば、2×2画素等のビニング処理を行ったデータを表示装置に出力してリアルタイム表示を行い、手術具等の位置確認等に使用するものがあるが、このような型式の出力信号を、そのまま特徴量の解析に使用することが可能となり、既存の撮影装置の出力信号を使用しても動態に係る特徴量の解析が行えることを示している。
例えば、一回の動態撮影で15のフレーム画像が生成され、これらに基づき上述のフレーム間差分画像の静止画像を作成した場合、撮影により15フレーム生成され、解析により中間生成データが14フレーム、最終的成果物としての静止画像(解析結果出力表示用画像)が1つ生成される。即ち、静止画像に比べて約30倍もの画像データが生成される。これら全てを施設内部で全て保存するには、膨大な保存設備が必要となる。一方、診断に使用された画像については医師法により明確に保存義務があるが、新しい手法である動態解析に伴い生成される画像群の保存義務対象は未定で、今後の法規制は流動的である。そこで、上記実施の形態においては、将来の保存義務化の予想される診断のエビデンスとなる解析結果データ、及び保存義務化の可能性のある解析結果データの生成に係る使用画像データ、及び中間生成データを施設内部にあるPACS10で保存し、保存義務化の可能性が低い解析結果の生成に使用されていない未使用画像データを外部保存サーバー60に保存する構成としている。
まず、大規模施設の例について説明する。ここでは、解析結果データを用いて診断を行う(診断に際して一連の動態画像の表示(各フレーム画像のパラパラ捲り表示)はしない)大規模施設を想定している。
図1Aに示す動態解析システム100のように、動態解析を内部で実施する大規模施設Aでは、上記実施の形態における保存形態の他、以下の(a−1)〜(a−7)の保存形態が挙げられる。なお、使用画像データと未使用画像データの弁別は解析サーバー8で実施される。
(a−1)解析結果データ、使用画像データ、中間生成データを施設内のPACS10で内部保存する。未使用画像データは外部のデータ保存サービスを利用して外部保存サーバー60に保存する。
(a−2)解析結果データ、使用画像データ、中間生成データを施設内のPACS10で内部保存する。未使用画像データは削除する。これは従来の静止画撮影系に於ける再撮影時と同様の運用方法である。
特に、胸部の動態撮影においては、以下の(a−3)、(a−4)の保存形態も考えられる。これらは、解析結果データとともに、使用画像データのうち肺野の動態に係る最も重要な画像データである、最大肺野面積フレーム画像及び最小肺野面積フレーム画像のみを内部保存する保存形態であり、内部のPACS10に保存するデータ量を更に抑えることができる。この保存方法は、既に施設内で運用されているPACS10に係る画像保存方法と最も親和性が高い。
(a−3)解析結果データ、及び使用画像データのうち最大肺野面積フレーム画像及び最小肺野面積フレーム画像を施設内のPACS10で内部保存する。それ以外の使用画像データ、未使用画像データ及び中間生成データは外部保存サーバー60に保存する。
この保存方法は、ブロックサイズを変更して再度解析処理等を行いたい場合に有効となる。
(a−4)解析結果データ、及び使用画像データのうち最大肺野面積フレーム画像及び最小肺野面積フレーム画像を施設内のPACS10で内部保存する。それ以外の使用画像データ、未使用画像データ及び中間生成データは削除する。
更に、内部のPACS10に保存するデータ量を更に抑えるには、以下の(a−5)、(a−6)の保存形態が考えられる。
(a−5)解析結果データのみを施設内のPACS10で内部保存する。使用画像データ、中間生成データは外部保存サーバー60に保存する。未使用画像データは削除する。
(a−6)解析結果データのみを施設内のPACS10で内部保存する。使用画像データ及び中間生成データは解像度を下げて(単純間引きやビニング処理のいずれでも良い)外部保存サーバー60に保存する。未使用データは削除する。
(a−7)解析結果データのみを施設内のPACS10で内部保存する。使用画像データ及び中間生成データは1画素あたりのビット数を減らして(例えば、各画素に非線形なLUTを適用して、1画素12ビットの画素値を8ビットに変換する等)、外部保存サーバー60に保存する。未使用データは削除する。
(a−8)解析結果データのみを施設内のPACS10で内部保存する。それ以外は全てそのまま、或いは解像度を下げて(単純間引きやビニング処理のいずれでも良い)、或いは、1画素あたりのビット数を減らして、或いは、解像度及び1画素あたりのビット数の双方を減らして、外部保存サーバー60に保存する。
ここで、従来の濃度コントラストに基づく静止画ベースの診断に於いては、FPDの各素子出力は12bit分解能(0〜4095諧調)を有することが一般的であった。
一方、本実施の形態における動態の解析処理に於いては、隣接フレーム間の差分値に基づく解析であるので、フレーム間差分値に於ける各素子の出力範囲は12bitよりもかなり小さな範囲に収まる。従って、解析に使用した一連の画像データの保存に際しては、1枚の画像のみRAWデータとし、他は当該画像を基準とするフレーム間差分値に置換して保存することで、データ保存容量を圧縮することができる。また、胸部の場合、最大肺野面積フレーム画像、及び/又は最小肺野面積フレーム画像をRAWデータで保存し、他フレーム画像は前記画像群を基準としたフレーム間差分値として保存することでも良い。いずれの方式でも、解析に使用した一連の画像データの再現が可能である。
図16に、施設内に動態撮影装置(例えば、既存の透視装置)は有するものの、解析サーバー8を有さない大規模施設Bを対象とした動態解析システム101の全体構成を示す。図16に示すように、動態解析システム101は、動態撮影システム40、PACS10を有する施設内システム31、外部保存サーバー61、解析センターの解析サーバー8、外部保存サーバー60を備えて構成されている。施設内システム31と解析センターの解析サーバー8はインターネットPNを介してデータ送受信可能である。外部保存サーバー61は、大規模施設Bと外部保存サービスの契約を締結している会社の外部保存サーバーであり、インターネットPNを介して大規模施設Bの施設内システムとデータ送受可能である。外部保存サーバー62は、解析センターと外部保存サービスの契約を締結している会社の外部保存サーバーであり、インターネットPNを介して解析センターの解析サーバー8とデータ送受可能である。なお、図16において、図1A〜図1Bに示す動態解析システム100を構成する装置、システムと同様の構成及び機能を有する装置、システムは、同名同符号で示している(図17A、図17B、図18についても同様である)。
(b−1)動態撮影後、施設内システム31から解析サーバー8へは全RAWデータを送信し(送信済みRAWデータはコピーを一時保存しておく)、解析サーバー8で使用画像データと未使用画像データを弁別し、解析結果データとともに使用画像データ情報(例えば、使用したフレーム画像の撮影順を示す番号等)を返信する。施設内システム31では、解析結果データ、及び一連のフレーム画像のうち使用されたフレーム画像をPACS10に保存する。その他、中間生成データ及び未使用画像データは解析サーバー8により外部保存サーバー62に送信し保存する。
(b−2)施設内システム31から解析サーバー8へは全RAWデータを送信し、解析サーバー8で使用画像データと未使用画像データを弁別し、解析結果データとともに使用画像データ及び中間生成データを返信する。施設内システム31では、解析結果データ、使用画像データ及び中間生成データをPACS10に保存する。未使用画像データは解析サーバー8により外部保存サーバー62に保存する。
(b−3)施設内システム31から解析サーバー8に解析に使用する1周期当たりのフレーム数を問い合わせ、解析サーバー8の回答に応じて解析に使用するフレーム画像のみを抽出して解析サーバー8に送信する。使用しなかったフレーム画像は施設内システム31から外部保存サーバー61に保存する。解析サーバー8では解析結果データのみを施設システム31へ返信し、施設内システム31では解析結果データをPACS10に保存する。使用画像データ及び中間生成データは解析サーバー8により外部保存サーバー62に保存する。
(b−4)施設内システム31から解析サーバー8に解析に使用する1周期当たりのフレーム数を問い合わせ、解析サーバー8の回答に応じて解析に使用するフレーム画像のみを抽出して解析サーバー8に送信する。使用しなかった未使用画像データは施設内システム31から外部保存サーバー61に保存する。解析サーバー8では解析結果データ、使用画像データ及び中間生成データを施設システム31へ返信し、施設内システム31では解析結果データ、使用画像データ及び中間生成データをPACS10に保存する。
なお、外部保存サーバー62は、上述のPACS10のサーバー装置のように、情報管理テーブルや解析結果テーブル等を有し、特徴量の値や患者情報等をキーとして解析サーバー8から送信された使用画像データや解析結果データを検索可能に保存しており、解析サーバー8から要求されたデータを検索して読み出し、解析サーバー8に送信することが可能である。解析サーバー8は、施設内システム31からのリクエストに応じて、外部保存サーバー62に保存されている使用画像データ等を読み出して追加の特徴量を算出して施設内システム31に送信することが可能である。過去の動態解析に係るデータ群(解析結果含む)を参照画像として利用する際に有効となる。
小規模施設は、開業医やクリニック等の小規模な医療施設である。施設内に装備されているデータ保存容量は一般的には小さく、勿論、解析サーバーは備えていない。従って、基本的には解析結果データのみしか施設内で保存できない可能性が大きい。このような状況に適合するデータフローを提供する必要がある。
図17Aに、施設内に動態撮影システム40を有する小規模施設Cを対象とした動態解析システム200の全体構成を示す。図17Aに示すように、動態解析システム200は、動態撮影システム40、HDD70を有する施設内システム32、解析センターの解析サーバー8、外部保存サーバー62を備えて構成されている。施設内システム32と解析センターの解析サーバー8はインターネットPNを介してデータ送受信可能である。解析サーバー8と外部保存サーバー62とは、インターネットPNを介してデータ送受可能である。
図17Aに示す動態解析システム200においては、施設内システム32で動態撮影が終了すると、解析サーバー8に一連のRAWデータを送信する。解析サーバー8は解析結果データのみを施設内システム32に返信し、施設内システム32では解析結果データのみをHDD70に保存する。なお、胸部の動態の解析の場合は、解析サーバー8は、解析結果データとともに、使用画像データのうち肺野の動態に係る最も重要な画像データである最大肺野面積フレーム画像及び最小肺野面積フレーム画像を施設内システム32に返信し、施設内システム32のHDD70で解析結果データとともに最大肺野面積フレーム画像及び最小肺野面積フレーム画像を保存することとしてもよい。
解析サーバー8は、下記の何れかにより、使用画像データ、中間生成データ及び未使用画像データを処理する。
・使用画像データ及び中間生成データを外部保存サーバー62に保存し、未使用画像データは削除する。
・使用画像データ及び中間生成データを解析時より解像度を下げて、或いは、1画素あたりのビット数を減らして、或いは、解像度及び1画素あたりのビット数の双方を減らして、或いは、使用画像データのうちの所定のフレーム画像(最大肺野面積フレーム画像及び/又は最小肺野面積フレーム画像)とフレーム間差分値データとの組合せで、外部保存サーバー62に保存し、未使用画像データは削除する。ここで、解析に使用した各フレーム画像の各画素が12ビットの諧調範囲に相当するデータを有しているものとすれば、これらの隣接フレーム画像間の差分値が示す値は、12ビット諧調範囲に比べ、狭い(少ない)諧調範囲内に収まるので、保存の為の容量を削減することができる。
・全てを外部保存サーバー62に保存する。
・胸部動態の解析の場合は、使用画像データのうち、最大肺野面積フレーム画像及び最小肺野面積フレーム画像を外部保存サーバー62に保存し、他は削除する。
その他、解析サーバー8では、施設からの要求に応じて、患者説明用の動態画像(RAWデータ、又は見やすいように加工した画像データ)、及び/又は教示用の動態画像を併せて施設に送信することとしてもよい。教示用の動態画像とは、動態画像による読影に習熟していない医師が動態画像の読影力を挙げるための動態画像で、例えば、典型的な症例の動態画像(RAWデータ、又は見やすいように加工した画像データ)等である。
図17Bに、施設内に動態撮影システムを有さない小規模施設Dを対象とした動態解析システム201の全体構成を示す。図17Bに示すように、動態解析システム201は、施設内システム33、イメージングセンターに設置され動態撮影を行うイメージングシステム41、解析センターの解析サーバー8、外部保存サーバー62を備えて構成されている。施設内システム33と解析センターの解析サーバー8とはインターネットPNを介してデータ送受信可能である。イメージングシステム41と解析サーバー8とはインターネットPNを介してデータ送受信可能である。解析サーバー8と外部保存サーバー62とは、インターネットPNを介してデータ送受可能である。なお、イメージングシステム41は、上述の動態撮影システム40と略同様の構成を有する。
図17Bに示す動態解析システム101においては、医師による被検者の問診後、医師より指示された被検者が訪れた外部のイメージングセンターのイメージングシステム41で動態撮影が行われ、一連のRAWデータが解析センターの解析サーバー8に送信され、解析が行われる。解析結果データは施設内システム33に送信される。また、外部保存サーバー62へ一連のRAWデータ、解析結果データ、中間生成データが保存される。
解析センターでは、外部保存サーバー62と提携し、施設側からの様々なリクエストに対応できるようになっている。例えば、患者説明用の動態画像や教示用の動態画像が施設内システム33から要求された場合、解析サーバー8は、外部保存サーバー62に保存されている、該当するRAWデータを読み出して、RAWデータ又は見やすいように加工した画像データを施設内システム33に送信する。また、追加の特徴量解析が依頼された場合、解析サーバー8は、外部保存サーバー62に保存されている該当患者のRAWデータを読み出して追加特徴量の解析を行い、解析結果データを施設内システム33へ送信する。また、例えば、施設から紹介状のある別の施設(大病院等)への解析結果データ及びRAWデータの送信が依頼された場合、解析サーバー8は、外部保存サーバー62から解析結果データ及びRAWデータを読み出し、指定されたアドレスに解析結果データ及びRAWデータを送信する。なお、紹介先で解析が行えるようであれば、RAWデータのみを送信することとしてもよい。
外部保存サーバー62における上記データの保存期間は、施設側から指定することとしてもよいし、一定期間リクエストがない場合は削除することとしてもよい。保存の基点は解析結果データを施設に送信したときである。
例えば、上記実施の形態においては、コンソール5において間引き処理を行うこととして説明したが、例えば、特許第4,546,174号公報のように、動態画像を撮影する撮影装置側(本実施の形態ではFPD9a)でビニング処理や単純間引き処理を行い、処理済みのフレーム画像をコンソール5に送信することとしてもよい。このようにすれば、FPDとコンソールとの間での画像データの転送時間も短縮することができるので、更に好ましい。
を適用することが可能である。また、本発明に係るプログラムのデータを通信回線を介して提供する媒体として、キャリアウエーブ(搬送波)も適用される。
1 ブッキー装置
11 制御部
12 検出器装着部
13 通信I/F
14 駆動部
15 バス
2 ブッキー装置
21 制御部
22 検出器装着部
23 通信I/F
24 駆動部
25 バス
3a 放射線源
3b 放射線源
3c 放射線源
4 クレードル
5 コンソール
51 制御部
52 記憶部
521 撮影管理テーブル
53 入力部
54 表示部
55 通信I/F
56 ネットワーク通信部
57 バス
8 解析サーバー
81 制御部
82 記憶部
83 入力部
84 表示部
85 通信部
86 バス
9a FPD
9b FPD
91 制御部
92 検出部
93 記憶部
94 コネクター
95 バッテリー
96 無線通信部
97 バス
6 操作卓
7 HIS/RIS
10 PACS
70 HDD
101、200、201、300 動態解析システム
60、61、62、63 外部保存サーバー
Claims (9)
- 被写体に連続的に放射線を照射可能な放射線源と、2次元状に配置された複数の検出素子を有し前記複数の検出素子のそれぞれにおいて前記放射線源により連続的に照射され前記被写体を透過した放射線を順次検出する放射線検出器と、を有し、前記放射線源と前記放射線検出器を用いて前記被写体の動態を撮影し複数のフレーム画像データを生成する撮影手段と、
前記撮影手段により生成された複数のフレーム画像データに基づいて前記被写体の動態を解析する解析手段と、
を備える動態解析システムであって、
前記複数のフレーム画像データを、前記解析手段による解析に使用する複数の画像データと解析に使用しない複数の画像データに弁別する弁別手段と、
前記解析手段により前記解析に使用する複数の画像データを用いて解析が行われた後、前記解析に使用された複数の画像データを保存する第1の保存手段と、
前記解析手段による解析に使用されなかった未使用の複数の画像データを保存する第2の保存手段と、
を備え、
前記弁別手段は、前記撮影手段により生成された複数のフレーム画像データの中から一定のフレーム間隔で抽出した複数のフレーム画像データを前記解析手段による解析に使用する複数の画像データとして弁別する動態解析システム。 - 前記第1の保存手段は、更に、前記解析手段により生成された解析結果データと前記解析手段において解析の過程で生成された中間生成データとを前記解析に使用された画像データに対応付けて保存する請求項1に記載の動態解析システム。
- 検索条件となる特徴量を指定するための入力手段と、
前記入力手段により指定された検索条件を満たす解析結果データ及びこれに対応する使用画像データを前記第1の保存手段から検索する検索手段と、
を備える請求項2に記載の動態解析システム。 - 前記弁別手段は、前記複数のフレーム画像データのフレームレート又は解像度に基づいて分別を行う請求項1〜3の何れか一項に記載の動態解析システム。
- 前記撮影手段は、前記被写体の動態を1周期以上にわたり撮影する請求項1〜4の何れか一項に記載の動態解析システム。
- 被写体に連続的に放射線を照射可能な放射線源と、2次元状に配置された複数の検出素子を有し前記複数の検出素子のそれぞれにおいて前記放射線源により連続的に照射され前記被写体を透過した放射線を順次検出する放射線検出器と、を有し、前記放射線源と前記放射線検出器を用いて前記被写体の動態を撮影し複数のフレーム画像データを生成する撮影手段と、
前記撮影手段により生成された複数のフレーム画像データに基づいて前記被写体の動態を解析する解析手段と、
を備える動態解析システムであって、
前記複数のフレーム画像データを、前記解析手段による解析に使用する複数の画像データと解析に使用しない複数の画像データに弁別する弁別手段と、
前記解析手段により前記解析に使用する複数の画像データを用いて解析が行われた後、前記解析に使用された複数の画像データを保存する第1の保存手段と、
前記解析手段による解析に使用されなかった未使用の複数の画像データを保存する第2の保存手段と、
を備え、
前記弁別手段は、前記複数のフレーム画像データのフレームレート又は解像度に基づいて分別を行う動態解析システム。 - 前記第1の保存手段は、更に、前記解析手段により生成された解析結果データと前記解析手段において解析の過程で生成された中間生成データとを前記解析に使用された画像データに対応付けて保存する請求項6に記載の動態解析システム。
- 検索条件となる特徴量を指定するための入力手段と、
前記入力手段により指定された検索条件を満たす解析結果データ及びこれに対応する使用画像データを前記第1の保存手段から検索する検索手段と、
を備える請求項7に記載の動態解析システム。 - 前記撮影手段は、前記被写体の動態を1周期以上にわたり撮影する請求項6〜8の何れか一項に記載の動態解析システム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012038085A JP5948950B2 (ja) | 2012-02-24 | 2012-02-24 | 動態解析システム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012038085A JP5948950B2 (ja) | 2012-02-24 | 2012-02-24 | 動態解析システム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2013172782A JP2013172782A (ja) | 2013-09-05 |
JP5948950B2 true JP5948950B2 (ja) | 2016-07-06 |
Family
ID=49266319
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012038085A Active JP5948950B2 (ja) | 2012-02-24 | 2012-02-24 | 動態解析システム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5948950B2 (ja) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6407573B2 (ja) * | 2014-06-03 | 2018-10-17 | キヤノンメディカルシステムズ株式会社 | 医用画像処理装置 |
JP6430865B2 (ja) * | 2015-03-17 | 2018-11-28 | 株式会社日立製作所 | X線透視撮影装置、透過x線画像変換方法およびプログラム |
JP2018110637A (ja) * | 2017-01-10 | 2018-07-19 | コニカミノルタ株式会社 | 動態画像処理装置 |
JP2018153297A (ja) * | 2017-03-16 | 2018-10-04 | コニカミノルタ株式会社 | X線動画像処理装置 |
JP7066362B2 (ja) * | 2017-09-25 | 2022-05-13 | キヤノンメディカルシステムズ株式会社 | 医用情報処理装置、システム及び方法 |
JP7077899B2 (ja) | 2018-09-27 | 2022-05-31 | コニカミノルタ株式会社 | 撮影制御装置、放射線画像撮影装置、放射線画像撮影システム及びプログラム |
CN117174260B (zh) * | 2023-11-02 | 2024-01-30 | 四川省肿瘤医院 | 医疗影像数据管理***及数据分析方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004164320A (ja) * | 2002-11-13 | 2004-06-10 | Canon Inc | データ保存装置、データ保存方法、プログラム、および記憶媒体 |
JP4728626B2 (ja) * | 2004-11-05 | 2011-07-20 | 株式会社東芝 | X線画像診断装置 |
JP5365578B2 (ja) * | 2010-05-14 | 2013-12-11 | 株式会社島津製作所 | 画像処理方法およびそれを用いた放射線撮影装置 |
JP2011255032A (ja) * | 2010-06-10 | 2011-12-22 | Shimadzu Corp | 画像処理方法およびそれを用いた放射線撮影装置 |
-
2012
- 2012-02-24 JP JP2012038085A patent/JP5948950B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
JP2013172782A (ja) | 2013-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10278661B2 (en) | Dynamic diagnosis support information generation system | |
JP5919717B2 (ja) | 動態医用画像生成システム | |
JP5948950B2 (ja) | 動態解析システム | |
JP5765343B2 (ja) | 胸部診断支援システム及びプログラム | |
JP6413927B2 (ja) | 動態解析装置及び動態解析システム | |
JP5672147B2 (ja) | 胸部診断支援情報生成システム | |
JP5874636B2 (ja) | 診断支援システム及びプログラム | |
US8948543B2 (en) | Exposure condition decision support system and method | |
JP5617559B2 (ja) | 放射線画像撮影システム | |
JP2009153677A (ja) | 動態画像処理システム | |
WO2017017985A1 (ja) | 撮影用コンソール及び放射線画像撮影システム | |
JP5625800B2 (ja) | 動態診断支援情報生成システム | |
Spin-Neto et al. | Radiographic observers' ability to recognize patient movement during cone beam CT | |
JP2019051322A (ja) | 動態解析システム | |
JP2016087325A (ja) | 胸部動態撮影支援システム | |
JP5625799B2 (ja) | 動態診断支援情報生成システム | |
JP5617577B2 (ja) | 胸部診断支援情報生成方法 | |
JP5772653B2 (ja) | 動態診断支援情報生成システム、動態診断支援情報生成方法及び動態解析装置 | |
JP2017169830A (ja) | 動態解析装置 | |
JP2014147844A (ja) | 動態診断支援情報生成システム | |
JP2012110398A (ja) | 診断支援情報生成システム | |
JP2016209267A (ja) | 医用画像処理装置及びプログラム | |
JP2020000475A (ja) | 動態画像処理装置及びプログラム | |
JP7143747B2 (ja) | 画像表示装置、画像表示方法及び画像表示プログラム | |
JP2013102850A (ja) | 医用画像撮影システム、医用画像処理装置及びプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20140807 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20150430 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150602 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20150728 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20151104 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160203 |
|
A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20160218 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20160510 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20160523 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5948950 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |