主要内容
このページの翻訳は最新ではありません。ここをクリックして,英語の最新版を参照してください。
iradon
逆ラドン変換
説明
例
フィルター補正ありとフィルター補正なしの逆投影法の比較
ファントムのイメージを作成します。イメージを表示します。
P =幻影(128);imshow (P)标题(原始图像的）
イメージのラドン変換を実行します。
R =氡(P, 0:179);
フィルター補正ありの逆投影法を実行します。
I1 = iradon (R, 0:179);
フィルター補正なしの逆投影法を実行します。
I2 = iradon (R, 0:179,“线性”，“没有”)；
再構成後のイメージを表示します。
图subplot(1,2,1) imshow(I1，[]) title(过滤后的反向投影的) subplot(1,2,2) imshow(I2，[]) title(未经过滤的反向投影的）
単一角度での逆投影法の検証
ファントムのイメージを作成します。
P =幻影(128);
イメージのラドン変換を実行し,45度の投影に対応する投影ベクトルを求めます。
R =氡(P, 0:179);r45 = R(:, 46岁);
この単一の投影ベクトルに対する逆ラドン変換を実行します。iradonの構文は,θがスカラーの場合,インクリメントとして扱われるため,これを直接計算することができません。投影ベクトルをコピーし,その結果を2で割ることで,タスクを実行できます。
I = iradon([r45 r45]， [45 45])/2;
結果を表示します。
imshow(我,[])标题(“45度反投影”）
入力引数
出力引数
アルゴリズム
iradonは,回転の中心が投影の中心点装天花板(大小(R, 1) / 2)であると仮定しています。
iradonは,フィルター補正逆投影法アルゴリズムを使って,逆ラドン変換を実行します。フィルターは周波数領域で直接設計され,投影のFFTが乗算されます。投影は,空間領域でのエイリアシングを避けてFFTの速度を上げるために,フィルター処理を行う前に2のべき乗になるようにゼロパディングされます。
参照
[1] Kak, a.c.， M. Slaney，计算机层析成像原理，纽约，纽约，IEEE出版社，1988。
拡張機能
R2006aより前に導入
你也可以从以下列表中选择一个网站:
