2016-07-29 6 views
2

私は、MATLABのImage Processing Toolboxに組み込まれているラドン関数を使用しています。今日まで、私は期待した結果を得たカスタム関数を使用していました。具体的には、Point Spread Function(PSF)の投影をいくつかの方向(基線は0/45/90/135度)で取得する数学的モデルを開発しています。あなたが例を実行する場合イメージ対角線を横切ってラドン変換した後のアーティファクト?

I = zeros(1000,1000); 
I(250:750, 250:750) = 1; 
theta = [0 45 90 135]; 
[R,xp] = radon(I,theta); 
figure;plot(R);legend('0°','45°','90°','135°'); 

は、あなたが135分の45°(対角線)のためのプロットは、アーティファクトを示していることがわかります。私は試していた問題が表示されます、本当に簡単な例を用意しました

カーブに沿った鋸歯状の形状をしています。最初は、私が使用しているサンプリンググリッド(偶数点)と関係があると思いました。ただし、奇数点のグリッドを使用する場合、問題はそこに残ります。ラドン変換はいくつかの方向に渡って累積的な積分でしかないので、私はこの結果をあまり理解していません。したがって、私はこの「鋸歯」パターンを得るべきではありません。

私は結果について本当に混乱しています。誰かが同じ問題を実験しましたか?

ありがとうございます。

答えて

0

randon()機能で実装されていると思う単純なフォワードプロジェクタを使用すると、エイリアシングアーティファクトが発生します。このアーティファクトを削除するには、サンプリング数を増やす必要があります(randon()はおそらく同じ数のファントムのサンプリングを使用しますが、その数をファントムサンプリングの数の2倍以上に増やすことをお勧めします) GEのCT画像再構成ソフトウェアで使用される駆動型プロジェクタなど、

関連する問題