図１本研究の概念図 1.3 評価対象胸部 X 線動画像には肺換気が X 線透過性ピクセル値の変化としてあらわれている 10)~13) これは呼吸により単位容積あたりの肺血管および気管支密度が変化するためである呼吸過程を撮影した胸部 X 線動画像の肺野内で計測したピクセル値を図

Size: px

Start display at page:

Download "図１本研究の概念図 1.3 評価対象胸部 X 線動画像には肺換気が X 線透過性ピクセル値の変化としてあらわれている 10)~13) これは呼吸により単位容積あたりの肺血管および気管支密度が変化するためである呼吸過程を撮影した胸部 X 線動画像の肺野内で計測したピクセル値を図"

わんどなかじゅく
7 years ago
Views:

1 動く軟組織 X 線動画像を対象とした肺換気血流コンプライアンス計測の試み研究責任者金沢大学医薬保健研究域保健学系量子医療技術学講座助教田中利恵共同研究者金沢大学医薬保健研究域保健学系量子医療技術学講座教授真田茂シカゴ大学放射線科医用物理学研究科准教授鈴木賢治金沢大学附属病院放射線部診療放射線技師作田啓太川嶋広貴 1. はじめに 1.1 背景現在肺機能の日常検査はスパイロメーターを用いた肺機能検査によって行われているしかし左右肺の総合的な機能を評価する手法であり局所診断は画像検査に頼らざるを得ない肺機能イメージングとしては肺シンチグラフィ Dynamic CT MRI PET 検査があるこれらの画像検査は治療戦略の決定や病態の質的診断に有用であるものの日常的に繰り返し行えるものではないもし肺機能情報がもっと簡便に取得できるようになれば治療効果判定や術後経過観察に大変有用である申請者らはこれまでに低コスト低被ばく X 線動画像検査法の開発を行ってきた 1)~9) 動画対応 FPD を用いて呼吸過程を撮影し画素値の変化を計測することで肺換気欠損部を X 線透過性変化量の減少領域として検出できることを解明したしかし骨陰影のある通常の胸部 X 線動画像を対象としており解析精度に多くの課題が残っている画像間位置合わせの妨げとなる肋骨や鎖骨を分離できれば肺血管気管支の画像間位置合わせ精度が向上し循環肺機能解析の精度向上が期待される近年 1 枚の胸部 X 線写真から画像処理により骨陰影を除去する技術 (ANN) が共同研究者であるシカゴ大学の鈴木らによって開発されたこの手法を胸部 X 線動画像に適用することで特別な装置や患者被ばくを増やすことなく動く軟組織 X 線画像を作成することができる 1.2 目的本研究では研究代表者らが開発してきた低コスト低被ばく X 線動画像検査法と共同研究者らが開発してきた画像処理による骨陰影除去を融合することで新しい循環呼吸機能イメージングを開発することを目指した ( 図 1) 動く軟組織 X 線動画像を対象とした肺換気血流コンプライアンスの計測を試みたので報告する 105

2 図１本研究の概念図 1.3 評価対象胸部 X 線動画像には肺換気が X 線透過性ピクセル値の変化としてあらわれている 10)~13) これは呼吸により単位容積あたりの肺血管および気管支密度が変化するためである呼吸過程を撮影した胸部 X 線動画像の肺野内で計測したピクセル値を図 2 に示す呼気で X 線透過性が低下ピクセル値は増加し呼気で X 線透過性が上昇ピクセル値は減少しているのが分かる従ってピクセル値の呼吸性変化量を計測することで肺の相対的な含気量を間接的に評価できると考えられるまた胸部 X 線動画像には血流動態も X 線透過性ピクセル値の変化としてあらわれている図 2 において心電図に同調して小刻みに変化する成分が心拍に伴う血流性変化であるこれは心拍出により単位容積あたりの肺血液量が変化する成人男性の平均的な肺内血液量 400~500ml 拍出量に伴う変動量 75ml ためである従ってピクセル値の血流性変化量を計測することで肺の血流動態を間接的に評価できると考えられるただし本法が計測しているのは肺胞レベルで行われているガス交換や肺血流そのものではなくそれらに関連する相対的なパラメータであることを留意しなければならない図２ピクセル値の呼吸性血流性変化発表論文[2] 106

3 2. 内容 2.1 画像の取得動画対応フラットパネルディテクタ (FPD) (CXDI-50RF) を搭載したポータブル X 線撮影装置 ( 試作機キヤノン ) を用いて最大努力呼吸の過程を立位正面背腹方向にて撮影した撮影条件は 120 kv 0.1 mas/pulse SID 1.0~1.2 m 5.0 fps とし 10 秒間 ( 吸気 5 秒 + 呼気 5 秒 ) に 50 フレームの胸部 X 線動画像を取得したこれらの撮影条件は被検者への総被ばく線量が通常の胸部単純 X 線撮影の2 方向 ( 正面 + 側面 ) の合計線量以下となるように設定した取得画像のマトリックスサイズは pixels ピクセルサイズは μm 撮像視野は cm 階調数は 16bits グレースケールである本研究は本学医学部の倫理委員会の承認を得て行なわれ被検者には撮影に関する十分な説明を行い同意を得た今回は合計 23 症例 (31-91 歳中央値 64 M:F=13:10) の胸部 X 線動画像を取得した対象は肺癌間質性肺炎肺挫傷皮下気腫気腫性嚢胞などの基礎疾患のある症例である 2.2 軟組織 X 線動画像および骨 X 線動画像の作成肋骨陰影低減 (BS) 処理は胸部 X 線写真における肺結節の検出精度向上を目的に開発された画像処理技術である 14) 米国では既に臨床実用され肺結節の検出率が 16.8% 向上したとの報告もある 15) 16) 本研究ではこの BS 処理を動画適応し世界で初めて肺血管および気管支などの軟部組織と肋骨や鎖骨などの骨陰影をそれぞれ分離した X 線動画像の作成を行った図 3 に BS 処理で作成した軟組織 X 線動画像と骨 X 線動画像を示す BS 処理ではその副産物として取り除いた骨陰影の画像も生成される本研究の主たる目的は BS 処理によって作成した軟組織 X 線動画像を対象とした肺機能解析だが, 研究過程で骨 X 線動画像得から肺機能評価をサポートする新しい診断情報の取得に成功したので併せて報告する図 3 胸部 X 線動画像を対象とした BS 処理 2.3 肺機能解析肺血流および肺換気によるピクセル値の変化は微小であり肉眼での評価は極めて困難であるそこでフレーム間差分および差分値の可視化が有用である図 4 は1 心拍のピクセル値の血流性変化量を可視化したマッピング画像である心室収縮期には血液が心室から送り出される様子を心室拡張期には血液が心室に流入する様子をとらえることができている図 5 に呼吸器疾患症例 ( 嚢胞性肺気腫 31M) の吸気過程の 2 フレーム間差分で作成したマッピング画像を示す肺シンチグラフィ上で確認される肺換気障害部はピクセル値変化量の減少領域として描出されている図 6 にマッピング画像作成アルゴリズムを示すまず肺野領域を認識後横隔膜動態から呼吸位相を推定し呼吸位相と息止め位相の画像に分ける 17) 18) 呼吸位相の画像を対象に肺機能を息止め位相の画像を対象に心機能をそれぞれ評価したフレーム間差分を行い算出した差分値をその値の大きさに応じて表示することで図 4-5 に示すような肺換気 / 肺血流マッピング画像を作成した健常者の肺換気および肺血流分布は左右対称であり立位では肺基底部ほど大きくなる傾向がある 8) 19) したがって正常分布からの逸脱同一被検者における左右肺での比較経時変化の有無などにより異常は検出されるこれま 107

4 では通常の胸部 X 線動画像を解析対象としてきたが本研究では肋骨陰影を除去した軟組織 X 線動画像を対象からマッピング画像を作成し肺換気肺血流肺コンプライアンスなどの肺機能に関連するパラメータの取得を試みた 2.4 肋骨動態解析図 7 に骨 X 線動画像を対象とした局所移動ベクトル計測の解析プロセスを示す. 格子状に区切られた領域ごとに隣り合うフレーム間で移動ベクトルを計測しその結果をオリジナルの画像上に重ね合わせて表示した. 得られたベクトル画像からは肋骨運動の向きと大きさを直感的に評価することができる図 4 1 心拍のピクセル値の血流性変化量を可視化した肺血流マッピング画像 (22F, 正常 ) 発表論文 [2] 図 7 骨 X 線動画像を対象とした局所移動ベクトル計測図 5 吸気過程の 2 フレーム間差分で作成した肺換気マップ ( 左 : 従来画像, 右 : 肋骨除去画像 ) と肺換気シンチグラフィ (31M, 嚢胞性肺気腫 ) 図 6 画像解析アルゴリズム発表論文 [2] 3. 成果 3.1 軟組織 X 線動画像を対象とした肺機能評価図 5 に軟組織 X 線動画像から作成した肺換気マッピング画像を示す従来画像では肋骨による動きアーチファクトが評価の妨げとなっていたが軟組織 X 線動画像から作成した肺換気マップでは肋骨による動きアーチファクトが軽減していることが分かる図 8 に肺癌症例 (66 歳男性 ) の肺血流マップおよび肺血流シンチグラフィを示す肺血流シンチグラフィでは左肺全体の血流低下を示す所見がみられたフレーム間差分値を可視化した肺血流マップでは右肺に比べ左肺でフレーム間変化の減少を示す分布を示した従来画像では鎖骨による動きアーチファクトが評価の妨げとなっていたが肋骨除去画像では鎖骨に 108

5 よる動きアーチファクトが軽減していることが分かる以上より肋骨除去処理が解析精度の向上に有用であることが明らかとなった肺換気マップ上で変化の大きな領域は呼吸により肺血管気管支に密度が大きく変化した領域であるすなわち肺換気や肺コンプライアンスが十分大きい領域であると推察される一方, 肺血流マップ上で変化の大きな領域は肺血液量が大きく変化した領域である今後更に症例を重ねてこれらの関連性の解明を行いたい図 8 吸気過程の 2 フレーム間差分で作成した肺換気マップ ( 左 : 従来画像, 右 : 肋骨除去画像 ) と肺換気シンチグラフィ (66M, 肺癌 ) 3.2 骨 X 線動画像を対象とした肋骨動態解析図 9 に骨 X 線動画像を対象とした局所移動ベクトル計測の結果を示す副産物として生成した骨 X 線動画像だがこの骨動画像を対象に肋骨動態を解析したところ肺機能評価をサポートする新しい診断情報としての有用であることが明らかとなった ( 下図 : 側弯症症例の肋骨運動低下 ) 側弯症( 先天的な背骨の湾曲 ) や外傷による肋骨損傷症例では肋骨運動が制約され呼吸機能障害をきたすことが知られているすなわち生命維持に欠かせない呼吸機能評価を可能にする肋骨動態情報は大変有用であるしかし通常の胸部 X 線動画像では骨気管支血管の複雑な動きを分離できず肋骨動態の単独評価は不可能であった技術的困難を理由に実用化を断念したが骨動画像を対象とすることでこれまで不可能だった肋骨単体の動態解析が可能になった図 9 肋骨動態解析の結果 (68F, 左肺癌側弯症 ) (a) 通常の胸部 X 線動画像,(b) 骨 X 線動画像 3.3 軟組織 X 線動画像を対象とした肺癌の動態追跡放射線治療分野での応用を期待できる成果も得られた図 10 に示すのはオリジナル動画像および軟組織 X 線動画像を対象とした標的追跡の結果であるオリジナル画像では肋骨陰影の影響で追跡エラーが発生しているが軟組織 X 線動画像では標的を正確に追跡することができたこのように BS 処理の動画応用の有用性が示された現時点では 1 枚あたりの処理に 15 秒かかっているため動画応用に向けた処理速度の高速化が課題である図 10 オリジナル動画像および軟組織動画像を対象とした標的追跡の結果 (84F 右肺癌 ) 3.4 既存技術と比較した利点近年の技術進歩で CT や MRI を用いて詳細な 3 次元形態情報が得られさらに動画像に近い疑似動態画像の取得も可能になったしかし CT はレントゲン検査の約 100 倍の被ばくがあり MRI は検査に約 1 時間かかり循環器呼吸器の 109

6 機能検査として日常的に行えるものではない申請者らがこれまで開発してきた FPD による低コスト低被ばく X 線動画像検査法によれば従来のレントゲン検査と同等のコスト被ばく時間で造影剤を使用することなく循環器呼吸器の機能評価が可能である回診車への搭載によりベッドサイド手術室集中治療室屋外での利用が可能になる救急医療や災害時救急における治療方針決定を簡便かつ迅速に行う一手法となることが予想される 4. まとめ動画対応 FPD を用いた機能イメージングによれば従来の単純 X 線検査時に付加的に機能情報を取得できる横隔膜動態や心機能など形態変化となって画像上に投影される機能情報は数値による定量化が有用であるまたピクセル値の変化となって画像上にあらわれる肺換気や肺血流情報はフレーム間差分とマッピング技術による可視化が有用である正常値からの逸脱や左右肺の比較により機能異常は検出可能である 2011 年には動画と静止画に対応する可搬型 FPD が開発されポータブルでの動画撮影も可能となった本研究の最終目標は動画対応 FPD を用いたポータブル X 線肺機能イメージング ( 診る聴診器 ) の開発である ( 図 11) 診断基準の確立臨床評価異常検出のアルゴリズム開発などが今後の課題である図 11 低コスト低被ばくポータブル X 線肺機能イメージング ( 診る聴診器 ) の概念図謝辞画像データの取得にあたりご協力いたキヤノン ( 株 ) および金沢大学附属病院放射線部のスタッフの皆様肋骨陰影低減処理にご尽力いただいた ( 株 ) 東陽テクニカの岸谷康氏に心から感謝申し上げますまた本研究助成の成果は平成 25 年度第 3 回技術交流助成金を支給いただき国際光工学学会第 26 回国際シンポジウム医用画像学会 2014 年 (SPIE2014 )( 2014/2/15 ~ 2014/2/20, サンディエゴ, 米国 ) にて発表したこのような機会を与えていただいた輕部征夫理事長をはじめ中谷財団関係の皆様に厚く御礼申し上げます参考文献 1) Rie Tanaka, Shigeru Sanada, Masaki Fujimura, et al. Ventilatory impairment detection based on distribution of respiratory-induced changes in pixel values in dynamic chest radiography: a feasibility study. IJCARS, 6(1), , ) Rie Tanaka, Shigeru Sanada, Masaki Fujimura a, et al. Development of pulmonary blood flow evaluation method with a dynamic flat-panel detector (FPD): quantitative correlation analysis with findings on perfusion scan. Radiological physics and technology, 3(1), 40-45, ) Rie Tanaka, Shigeru Sanada, Masaki Fujimura a, et al. Pulmonary blood flow evaluation using a dynamic flat-panel detector: Feasibility study with pulmonary diseases. IJCARS, 4(5); , ) 田中利恵, 真田茂, 藤村政樹, 他. 動画対応フラットパネルディテクタによる肺機能画像診断法 : 肺シンチグラフィ所見との比較. 日本放射線技術学会雑誌.65(6); ,

7 5) 田中利恵, 藤村政樹, 安井正英, 他. 胸部 X 線動態撮影法による換気 - 血流 mismatch 症例を対象とした V/Q study. 医用画像情報学会雑誌.26(3);68-72, ) Rie Tanaka, Shigeru Sanada, Masaki Fujimura, a, et al. Development of functional chest imaging with a dynamic flat-panel detector (FPD). Radiological physics and technology. 1(2), , ) Rie Tanaka, Shigeru Sanada, Katsumi Tsujioka a, et al. Development of a cardiac evaluation method using a dynamic flat-panel detector (FPD) system: a feasibility study using a cardiac motion phantom. Radiological physics and technology. 1(1), 27-32, ) Rie Tanaka, Shigeru Sanada, Nobuo Okazaki a, et al.detectability of Regional Lung Ventilation with Flat-panel Detector-based Dynamic Radiography. J.Digit.Imag. 21(1), , ) Rie Tanaka, Shigeru Sanada, Nobuo Okazaki a, et al. Evaluation of pulmonary function using breathing chest radiography with a dynamic flat-panel detector (FPD): Primary results in pulmonary diseases. Invest Radiol, 41(10), , ) Goodman LR: Felson's Principles of Chest Roentgenology. A programmed text, 3rd ed. Philadelphia, London, Toronto: W B Saunders Co., ) Squire LF, Novelline RA: Fundamentals of Radiology, 4th ed. Cambridge, Massachusetts, and London: Harvard University Press, ) Toffolo RR, Beerel FR: The autofluoroscope and 133-Xe in dynamic studies of pulmonary perfusion and ventilation. Radiology, 94: , ) George RB, Weill H, Tahir AH: Fluorodensimetric evaluation of regional ventilation in chronic obstructive pulmonary disease. South Med J. 64, , ) Suzuki K, Abe H, MacMahon H and Doi K: Image-processing technique for suppressing ribs in chest radiographs by means of massive training artificial neural network (MTANN), IEEE Trans. Med. Imaging. 25, , ) Freedman MT, Lo SC, Seibel JC and Bromley CM: Lung nodules: improved detection with software that suppresses the rib and clavicle on chest radiographs, Radiology 260, , ) Li F, Hara T, Shiraishi J, Engelmann R, MacMahon H and Doi K: Improved detection of subtle lung nodules by use of chest radiographs with bone suppression imaging: receiver operating characteristic analysis with and without localization, Am. J Roentgenol. 196, , ) Tanaka R, Sanada S, Kobayashi T, Suzuki M, Matsui T, Matsui O: Computerized methods for determining respiratory phase on dynamic chest radiographs obtained by a dynamic flat-panel detector (FPD) system. Journal of Digital Imaging, 19, 41-51, ) Xu XW, Doi K: Image Feature Analysis for Computer-Aided Diagnosis: Accurate Determination of Ribcage Boundary in Chest Radiographs. Medical Physics. 22, , ) West JB: Ventilation how gas gets to the alveoli. In:West JB Ed. Respiratory Physiology The Essentials, 6thedn. Philadelphia: Lippincott Williams & 111

8 Wilkins, 2000,pp 11Y19 発表論文 [1] Rie Tanaka, Shigeru Sanada, Keita Sakuta, Hiroki Kawashima. Quantitative analysis of rib kinematics based on dynamic chest bone images: preliminary results. Journal of Medical Imaging,2015 (In press) [2] 田中利恵, 真田茂. 低コスト低被ばくポータブル X 線肺機能イメージング ( 診る聴診器 ) の開発. 医用画像情報学会誌, 2014(In press) [3] Rie Tanaka, Shigeru Sanada, Makoto Oda, Keita Sakuta, Hiroki Kawashima. Functional chest radiography based on dynamic analysis of soft-tissues in bone suppression images. Computer Assisted Radiology and Surgery 28 th International Congress and Exhibition, Journal of computer assisted radiology and surgery,9(1):s22, [4] Rie Tanaka, Shigeru Sanada, Makoto Oda, Mitsutaka Suzuki, Keita Sakuta, Hiroki Kawashima. Application of bone suppression technique to real-time tracking radiotherapy. RSNA th99 Scientific Assembly and Annual Meeting pdf 2013 [5] Rie Tanaka, Kenji Suzuki, Shigeru Sanada, Makoto Oda, Mitsutaka Suzuki, Keita Sakuta, Hiroki Kawashima. Quantitative analysis of rib movement based on dynamic chest bone images: Preliminary results. The international society for optical engineering. Medical imaging Proc. of SPIE, Proc. of SPIE Vol. 9034, , 1-6. [6] Rie Tanaka, Shigeru Sanada, Makoto Oda, Mitsutaka Suzuki, Keita Sakuta, Hiroki Kawashima and Hiroji Iida. Improved accuracy of IGRT with bone suppression image-processing. Journal of computer assisted radiology and surgery, 8(1):S41-42,2013 [7] Rie Tanaka, Shigeru Sanada, Makoto Oda, Mitsutaka Suzuki, Keita Sakuta, Hiroki Kawashima and Hiroji Iida. Functional chest radiography with a portable dynamic FPD: Detection of blood flow impairment based on Circulation map during cardiac pumping. Journal of computer assisted radiology and surgery, 8(1):S ,2013 [8] Rie Tanaka, Shigeru Sanada, Makoto Oda, Mitsutaka Suzuki, Keita Sakuta, Hiroki Kawashima. Improved accuracy of image guided radiation therapy (IMRT) based on bone suppression technique, IEEE 2013 Nuclear Science Symposium and Medical Imaging Conference Record (NSS/MIC), CFP13NSS-DVD, ISBN:

動くレントゲン

動くレントゲン FPD 1985 1998 2000 2002 FPD X X X FPD X CAD X FPD X CAD X X 1970 Autofluoroscope Toffolo RR, et al. The autofluoroscope and 133-Xe in dynamic studies of pulmonary perfusion and ventilation. Radiology,