临床画像技术学Ⅱ.ppt

課題５　　プログラム Ramp.c を改良して 実空間 SheppLoganフィルタを計算するプログラム SheppLogan.c を作って下さい。 　　　周波数空間SheppLogan の式 周波数空間SheppLogan の式の f0 は高周波遮断周波数。　 f0 をプログラムで設定できるように printf(\n\n Cut off freq. (0～0.5) = ); scanf(%lf,CutOff ); を入れてください。 ここで CutOff の単位は cycles/pixel 。　SheppLoganの式に 使われる f0 の単位は cycles / 画像の１辺長 = cycles / 256 pixel なので　f0 = 256.0 * CutOff ; と変換してから式に代入してください。 サイノグラムの各スライスの１次元配列は、 ３２方向の各々の角度から収集されたデータ。 このデータから、収集された各々の角度に傾いた ２次元透視画像（ Pθ ）を作成する。 //------------ Generate Pth[ ][ ][ ] (Pθ) data --------------- PAI = 3.1415926/180.0; dT = 180.0/(double)MAXFRAME; for(T=0;TMAXFRAME;T++){