CT 画像における分解能測定とその評価 2005.8.27 第 6 回千葉県 CT 研究会 名古屋大学医学部保健学科 市川勝弘.

Slides:



Advertisements
Similar presentations
模型を用いたジェットコターの 力学的原理の検討 06522 住友美香 06534 秦野夏希. 平成22年度 卒業研究発表 山田研究室 研究目的 ジェットコースターのコースは、どのような計算に 基づいて作られているのか、研究を通じて理解し、 計算を用いた模型製作を行う。
Advertisements

Determining Optical Flow. はじめに オプティカルフローとは画像内の明る さのパターンの動きの見かけの速さの 分布 オプティカルフローは物体の動きの よって変化するため、オプティカルフ ローより速度に関する情報を得ること ができる.
X線CTと中性子線CTボリュームデータの
第24回 画像研究会 in 金沢 乳房画像におけるディジタルの 有用性について  名古屋大学医学部保健学科 市川勝弘.
CsIシンチレータと マルチアノードPMTを用いた 硬X線撮像装置の性能測定
HOG特徴に基づく 単眼画像からの人体3次元姿勢推定
第1章 第1節 情報のディジタル化のしくみ 4 音の表現 5 画像の表現
復習.
第1章 第1節 情報のディジタル化のしくみ 4 音の表現 5 画像の表現
アナログとディジタル 五感 視覚、聴覚、味覚、臭覚、触覚 埼玉県立越ヶ谷高等学校・情報科.
画像処理工学 2011年10月27日 担当教員 北川 輝彦.
経営情報 #1 デジタル表現 / 2003 (春) 安田豊 1.
情253 「ディジタルシステム設計 」 (2)modem2
永久磁石を用いた 残留ガスモニターの製作 環境計測 西村荒雄.
9月27日 パラボラミラーによる ミリ波ビーム絞り
CTにおけるMTFとWS,そしてその測定法
デジタル信号処理①
AOによる 重力レンズクェーサー吸収線系の観測 濱野 哲史(東京大学) 共同研究者 小林尚人(東大)、近藤荘平(京産大)、他
首都大学東京 都市教養学部数理科学コース 関谷博之
わかりやすいマルチスライスCTにおける画像再構成
(質問)  体軸分解能を評価するための「SSPの測定」
画像工学 2011年10月6日 担当教員 北川 輝彦.
3.8 m望遠鏡主鏡エッジセンサ 開発進捗 京都大学 理学研究科 M2 河端 洋人.
16 列MDCT による心臓CT の撮影技術 日本放射線技術学会中国・四国部会セミナー 名古屋市立大学病院 中央放射線部
ディジタル回路 1. アナログ と ディジタル 五島 正裕.
1. アナログ と ディジタル 五島 正裕.
フリッカー値とその測定方法について 17070540       松本 春菜.
ー 第1日目 ー 確率過程について 抵抗の熱雑音の測定実験
光の干渉.
臨床画像に密接に関わる画質特性と撮影パラメータ
京大岡山 3.8m 望遠鏡 分割鏡制御に用いる アクチュエータの特性評価
Astro-E2衛星搭載 XISの データ処理方法の最適化
大阪大学 大学院工学研究科 極限光通信工学領域 井上研究室 欅田 直也・橘 遼太郎・隅田 拓也・高 祥史
画像工学 2012年10月3日 担当教員 北川 輝彦.
6/18 Meeting用.
某健保契約の胸部CT検診報告 医療法人新虎の門会 新浦安虎の門クリニック 南 香織 永山 健二 角谷 美佳 大前 利道 斉藤 智 大前 由美
授業展開#3 アナログとデジタル.
ディジタル画像の物理評価 -その実践と活用法- 名古屋市立大学病院 中央放射線部                  市川勝弘.
画像処理 基礎.
X線CTにおけるファントム中の エネルギー変化についての検討
第4回 信号表現とエリアシング.
京大理 高エネルギー研究室 修士課程一年 田口 誠
人工知能特論 9.パーセプトロン 北陸先端科学技術大学院大学 鶴岡 慶雅.
位相カメラの進捗状況 京都大学修士1回 横山 洋海.
細胞の形と変形のための データ駆動型解析手法
東京農業大学 東京情報大学 附属第一高等学校・中等部 附属第二高等学校 附属第三高等学校・中等部
治療用フィルムによる線量分布測定の 基礎的検討Ⅱ
情報機器と情報社会のしくみ Web素材利用
高分子電気絶縁材料の誘電特性計測を用いた劣化診断に関する研究
高分解能ビーム軌道傾きモニターの設計開発
高エネルギー陽子ビームのための高時間分解能 チェレンコフビームカウンターの開発
従来CTとの検査比較( 頭部検査) 1 Slice 16 Slices 1 Slice 16 Slices ルーチン検査 撮影範囲
ディジタル信号処理 Digital Signal Processing
X線CCD検出器 ーCCD‐CREST(deep2)ー の性能評価と性能向上 (京阪修論発表会)
アナログ と ディジタル アナログ,ディジタル: 情報処理の過程: 記録/伝送 と 処理 において, 媒体(メディア)の持つ物理量 と
TNRFにおける中性子ラジオグラフィ システムの高度化 安田 良1、野島健大1、片桐政樹2
名古屋市立大学大学院システム自然科学研究科 MIRU2009: 第12回 画像の認識・理解シンポジウム
ディジタルと情報量を あらわす単位 情報C ③.
はやぶさ試料(RA-QD )の X線CT解析 – X線CT岩石学の適用例 - X線CT解析の結果に基づいて試料を切断し分析
イオン源ビームオンラインモニター・スキャナー
落下水膜の振動特性に関する実験的研究 3m 理工学研究科   中村 亮.
実験結果速報 目的 装置性能の向上 RF入射実験結果 可動リミター挿入 RFパワー依存性 トロイダル磁場依存性 密度依存性
大阪工業大学 情報科学部 情報システム学科 学生番号 B02-014 伊藤 誠
AIRT40+TONIC2 for JARE53/54 Winter-over Observation 新光学系 備忘録
TOBAの現状と今後の計画 坪野研輪講 2012年2月22日 岡田健志.
Simulation study for drift region
23 造波機構における水位計の製作 1 はじめに 4 再現性の低下要因の実験 水位計の再現性の向上を目的としている. 2 実験装置
文字フォントについて.
TES型カロリメータのX線照射実験 宇宙物理実験研究室 新井 秀実.
アナログ と ディジタル アナログ,ディジタル: 情報処理の過程: 記録/伝送 と 処理 において, 媒体(メディア)の持つ物理量 と
Presentation transcript:

CT 画像における分解能測定とその評価 第 6 回千葉県 CT 研究会 名古屋大学医学部保健学科 市川勝弘

報告 X 線コンピュータ断層撮影装置の性能評価に関する基準 ( 第二次勧告 ) 班名 医学放射線学会 CT 性能評価委員会 掲載 日本医師会雑誌, 88(8), 759 ‐ 771(1989) 報告 X 線 CT 装置性能評価に関する基準(案) 班名 日本放射線技術学会 X 線 CT 装置性能評価検討班 掲載 日放技学誌, 47(1) , (1991) 報告 ラセン CT の物理的な画像特性の評価と測定法に関する報告 班名 日本放射線技術学会 ラセン CT 性能評価班 掲載 日放技学誌, 53(11) , (1997) CT の性能(画質)評価に関する勧告,報告 American Association of Physicists in Medicine. Report no. 1 Phantoms for performance evaluation and quality assurance of CT scanners ( 1977 )

分解能とは どこまで小さな物体を分離して識別できるか。 どこまで小さな物体を認識し得るか。 を表す能力。 物体のコントラストなどを規定しない。

ノイズの影響のない状態で どこまで小さな物体を分離して識別できるか。 解像度とは を表す能力 十分なコントラストの物体の分解能とも言える。

CT の場合 ・線量はノイズに大きく影響する ・ディテクタサイズとビュー数は解像度(分解能) に影響する。 を基本とするが.... 実際は... ・再構成関数と後処理も,ノイズと解像度に影響する。 よって,幅広く対応するために 高コントラスト分解能と低コントラスト分解能 が定義された。

高コントラスト分解能と低コントラスト分解能

高コントラスト分解能

高コントラストな物体に対する分解能 ノイズの影響がない場合の分解能=解像度 よって 用途 ・装置の限界の解像度を見る。 ・再構成関数の特性を見る。

ノイズの影響がない場合の分解能=解像度 そんな状態があるの?実現できるの? 高コントラストならそういう状態なの?

シミュレーションプログラムで確認します。 ノイズが増えても高コントラストファントムの画像とプロファイルは ほとんど影響受けません。低コンは,認識不可能となっています。 ・高コン,低コン共に,同量のノイズを付加。

ノイズ付加による画像の変化 ノイズは、低コントラストな物体にだけに影響

アクリル円筒 エアーホール 0.3mm ~ 2mm 100mm スキャン画像 高コントラスト分解能測定ファントム

MTFと高コンファントム画像の関係 mm 1.0cycles/mm は, 0.5mm の解像を示す。 MTF に応じて,画像も 忠実に変化している。

MTF による画像の連続的変化を シミュレーションプログラムで見てみます。 (割愛)

水 0.1mmφ 銅ワイヤ 直径 50mm 程度の 樹脂製シリンジ (例:オムニパーク ml ) 活栓 MTF測定ファントム 樹脂製 150ml シリンジは,太さも手頃。 0.1mm の銅ワイヤー はビニール線の中から 1 本抜いて使用。

プロファイル 仮想スリット PSF Z 方向 金属ワイヤ 約 0.1mmφ ワイヤーファントム MTF の計算過程 これをフーリエ変換

MTF 測定への影響因子 ・ノイズはあまり関係しない。 (ファントムサイズに見合った線量(電流)) ・電圧は,ほぼ関係しない。 ・ビュー数も,ほとんど影響しない。 ・再構成関数は大きく影響する。 ・焦点サイズも影響する。 ・周辺部では,回転速度が影響する。 ・ QQ オフセットやフライングフォーカスも影響する。

関数による MTF の変化

AB20 AB30 AB35 AB40 SOMATOM VZ の測定結果 番号順になっているとは限らなかった。 よってMTF測定が有用であった。

AB30 AB35 Q:どちらが高解像度でしょうか?

A . 解像度は,どこまで細かく解像するか(解像する 最高周波数)を示すので,答えは AB35

Center90mm off-center

Spatial frequency (cycles/mm) sec 0.5sec Modulation transfer factor 回転速度による MTF の変化 周辺部の MTF

低コントラスト分解能

低コントラストな物体に対する分解能 ノイズの影響を見る=ノイズを測る よって 用途 ・装置のノイズ特性(感度)を見る。 ・再構成関数の影響を見る。 低コントラスト分解能

アクリル円筒 周囲 CT 値 -5 直径 20,10,5,4,3,2mm 140mm ポリウレタン等 ( CT 値 50 付近) スキャン画像 低コントラスト分解能測定ファントム

SD=5 電流なら 200mA スライス厚なら 10mm 厚 SD=7 100mA 5mm 厚 SD=13 約 50mA 2mm 厚 CT 画像におけるノイズ

標準関数シャープな関数ソフトな関数 SD は関数の違い(解像度の違い)を考慮すべき。 同じ画像を関数を変えて表示 SD=5SD=2SD=8

ノイズの影響を見る=ノイズを測る ・ほんとにノイズしか見てないの。 ・物体に対する分解能を見るので,解像度 は影響しないの?

ソフトな関数 =解像度悪い シャープな関数 =解像度良い 10mm 5mm 2mm 解像度に関わらず、認識できる=解像度の影響をあまり うけない。 低コンファントム、ノイズがない状態では? 周囲=0物体= -20 周囲=0物体= -20

通常関数高解像度関数 ノイズが加わると 高解像度でも認識能力は変わらない⇒ノイズ特性だけが見れる

通常関数高解像度関数 ノイズが加わると ある程度コントラストがある場合は,ソフトな関数は見やすい

シミュレーションプログラムで 確認してみましょう。 (割愛)

関数(解像度)が違いSDが同じ(= 5.3 )の場合 軟部用標準高解像度 Q:さて,どれが一番低線量でしょうか?

A. 軟部関数 軟部関数は,もともと SD が小さい。その小さい SD を持ち上げるのは,線量を下げることになる。 よって,低線量。 高解像度関数は, SD が高い。その SD を下げるには 線量増加となり,高線量。

AB20 AB30 AB35 AB40 SOMATOM VZ の測定結果 AB30 AB35 AB40 低周波を強調する関数は,目障りなノイズ が増える。 視覚的に番号順になっているとは限らなか った。よってMTF測定が有用であった。

低コントラスト分解能(ノイズ)への影響因子 ・高コン(解像度)ほとんど関係しない。 ・よって再構成関数もほとんど関係しない。 ・すなわち,ノイズ特性を測ることができる。 ・ただし,線量との関係を見る場合には,解像度 を統一するべきである。 Q:それはなぜでしょうか?

A. 関数によって, SD が変わる。すなわち解像度に よって SD が変わる。機種が違うものを比較する 場合に,標準関数だからといっても解像度が違えば SD による比較は不可。

ディジタル画像としての性質 今は亡き,ハム太

CT 画像データと解像との関係 CT 画像の解像度= CT の性能の解像度 と ピクセルの解像度 マトリクス数 512 X 512 階調数 12 ビット~ 16 ビット 4096 階調以上 ピクセルサイズ= FOV に依存 FOV=350mm -> ピクセルサイズ =0.68mm FOV=200mm -> ピクセルサイズ =0.4mm FOV=120mm -> ピクセルサイズ =0.23mm

Raw データの情報量を満たすピクセルサイズ 0.5mm の解像が必要な場合、 0.5mm は 1.0cycles/mm と考え、それの忠実再現には 2cycles/mm のサンプリングが必要。 (サンプリングの定理) 最大の FOV は 500mm 500mm / 0.5mm = 1000 ローデータの情報量をそのままに表す画像は 0.5mm ピクセルサイズで 1000X1000 ピクセルあればよい. *究極は 1024X1024 再構成か?! 0.5mm

ズーミング FOV200mm 単純拡大 FOV200mm 0.5mm 必要 -> 余裕を見て 200mm のズーミング ....しかし、 0.5mm の解像はあるか? 0.5mm 1.00 標準関数なら単純拡大で十分である。( 0.7mm x 512 = 350mm で十分)

FOV200mmFOV350mm 単純拡大 ピクセルサイズ= 0.4 ピクセルサイズ= mm 1.00 胸部 High-resolution CT では有効 ....

ディジタル画像としての CT 画像(体軸方向) スライス厚 → サンプリングアパーチャ スライス間隔 → サンプリング間隔

Z 方向 スライス厚のスムーズフィルタ こんな被写体をとってるのと同じ スライス厚 Z 軸のプロファイル

TOPO 画像を使ってもっと具体的に 10mm5mm 3mm1mm 1mm のスライス厚では,構造をやや正確に表せることがわかる。

スライス厚 1 mmで、再構成間隔1mmは一般的...でも 0.5mm の解像には、 0.5mm 以下の間隔が被写体情報を 満足した間隔だといえる。 Z Z スライス厚 1mm の分解能は? ー>少なくとも 0.5mm はある。 情報量で考えるなら、再構成間隔はせまく 再構成間隔

MTF スライス厚 0.5mm なら,情報量を生かすことができる。 スライス 1.0mm では, 0.5cycles/mm ( 1mm )までしか再現できない。

スライス間隔 =1mm スライス間隔 =0.5mm

ディジタル画像としてのCT画像 ・ピクセルサイズ=解像したいサイズ ・スライス厚 /2 =情報量を生かす再構成間隔

ご清聴ありがとうございました。