ifanbeam
逆ファンビーム変換
説明
例
サンプル イメージを作成します。関数 phantom はファントム頭部イメージを作成します。
ph = phantom(128);
ファントム頭部イメージのファンビーム変換を作成します。
d = 100; F = fanbeam(ph,d);
ファンビーム表現からファントム頭部イメージを再構成します。元のイメージと再構成したイメージを表示します。
I = ifanbeam(F,d); imshow(ph)
figure imshow(I);
サンプル イメージを作成します。関数 phantom はファントム頭部イメージを作成します。
ph = phantom(128);
イメージのラドン変換を作成します。
P = radon(ph);
変換をパラレル ビーム投影からファンビーム投影に変換します。
[F,obeta,otheta] = para2fan(P,100,... 'FanSensorSpacing',0.5,... 'FanCoverage','minimal',... 'FanRotationIncrement',1);
ファンビーム データからイメージを再構成します。
phReconstructed = ifanbeam(F,100,... 'FanSensorSpacing',0.5,... 'Filter','Shepp-Logan',... 'OutputSize',128,... 'FanCoverage','minimal',... 'FanRotationIncrement',1);
元のイメージと変換後のイメージを表示します。
imshow(ph)
figure imshow(phReconstructed)
入力引数
ファンビーム投影データ。numSensors 行 numAngles 列の数値行列として指定します。numSensors はファンビーム センサーの数で、numAngles はファンビーム回転角の数です。F の各列には 1 つの回転角のファンビーム センサーのサンプルが格納されます。
データ型: double | single
ファンビーム頂点から回転中心までの距離 (ピクセル単位)。正の数値として指定します。ifanbeam は、回転の中心が投影の中心点であると仮定し、ceil(size(F,1)/2) として定義します。
データ型: double | single
名前と値の引数
オプションの引数のペアを Name1=Value1,...,NameN=ValueN として指定します。ここで、Name は引数名で、Value は対応する値です。名前と値の引数は他の引数の後に指定しなければなりませんが、ペアの順序は重要ではありません。
例: I = ifanbeam(F,D,FanRotationIncrement=5) は 5 度のファン回転角増分を指定します。
R2021a より前では、コンマを使用して名前と値をそれぞれ区切り、Name を引用符で囲みます。
例: I = ifanbeam(F,D,"FanRotationIncrement",5) は 5 度のファン回転角増分を指定します。
ファンビーム回転角の増分 (度単位)。正のスカラーとして指定します。
データ型: double
ファンビーム センサーの配置。"arc" または "line" として指定します。
値 | 意味 | 図 |
|---|---|---|
| センサーは、円弧状に等角度で配置されます。弧の中心はファンビーム頂点です。
|
|
| センサーは、x' 軸に平行なラインに沿って等間隔に配置されます。最も近いセンサーは回転中心から距離
|
|
周波数領域でのフィルター処理で使用するフィルター。次の表のいずれかの値として指定します。詳細については、iradon を参照してください。
値 | 説明 |
|---|---|
| Cropped Ram-Lak またはランプ フィルター。このフィルターの周波数応答は | |
| 関数 |
| 関数 |
| ハミング ウィンドウと Ram-Lak フィルターを乗算します。 |
| ハン ウィンドウと Ram-Lak フィルターを乗算します。 |
"None" | フィルター処理を行いません。ifanbeam はフィルター処理されていないデータを返します。 |
データ型: char | string
周波数軸を再スケーリングするスケール係数。(0, 1] の範囲の正の数値を指定します。FrequencyScaling が 1 未満である場合、フィルターは、正規化周波数の周波数範囲 [0, FrequencyScaling] に適合するように圧縮されます。FrequencyScaling を超える周波数については、すべて 0 に設定されます。詳細については、iradon を参照してください。
データ型: double
パラレルビーム データとファンビーム データの間で使用される内挿のタイプ。次のいずれかの値として指定します。
"nearest" - 最近傍点
"linear" — 線形 (既定の設定)
"spline" — 区分的 3 次スプライン
"pchip" — 区分的 3 次エルミート (PCHIP)
データ型: char | string
再構成イメージのサイズ。正の整数として指定します。イメージには同じ数の行と列が含まれます。
OutputSize を指定した場合、ifanbeam は、イメージを大きくするか、または、小さくしますが、データのスケーリングには変化がありません。
メモ
関数 fanbeam を使用して投影を計算した場合、再構成イメージのサイズは元のイメージと異なる場合があります。
OutputSize を指定しない場合、サイズは以下によって自動的に計算されます。
OutputSize = 2*floor(size(R,1)/(2*sqrt(2)))ここで、R は iradon によって使用されるパラレルビーム投影データの長さです。詳細については、アルゴリズムを参照してください。
データ型: double
出力引数
ヒント
アルゴリズム
ifanbeam はファンビーム データをパラレル ビーム投影に変換し、フィルター補正逆投影アルゴリズムを使用して逆ラドン変換を行います。フィルターは周波数領域で直接設計され、投影の FFT が乗算されます。投影は、空間領域でのエイリアシングを避けて FFT の速度を上げるために、フィルター処理を行う前に 2 のべき乗になるようにゼロ パディングされます。
参照
[1]
バージョン履歴
R2006a より前に導入
