2

我想知道是否有人遇到過這個問題。matlab中濾波反投影/ iradon的定性和定量分析

我可以從matlab中重建類似於原始圖像的圖像,但實際值總是不同的。

例如,原始圖像的值在0到1的矩陣範圍內,而我的重建圖像範圍從-0.2到0.4。

儘管重建圖像看起來與原始圖像相似,但是圖像中的數據具有不同的比例。

這是我的意思的示例代碼。

p=phantom(64); 
theta=0:1:179; 
r=radon(p,theta); 
ir=iradon(r,theta); 

figure 
subplot(1,2,1);imagesc(p) 
subplot(1,2,2);imagesc(ir) 

回答

0

這些結果並不完全是我發現的。

>> min(min(ir)) 

-0.0583 

>> max(max(ir)) 

0.9658 

請記住,逆氡變換隻能近似重建原始圖像。只有180個觀點,肯定會有一些分歧。

由於像素必須投影到新的座標系並重新投影 - 在投影和背投影過程中,Radon變換本質上會導致某些信息丟失。這會導致重建圖像稍微降級。 Radon變換不像傅里葉變換那樣具有相同的可逆性。

爲了獲得更好的效果,請嘗試使用更大的圖像尺寸和更多的視角。

p=phantom(256); 
theta=0:0.01:179; 

,並嘗試使用不同的過濾器(在F中F.B.P.)如Shepp-洛根,這減少了高頻率和減輕過沖。

ir=iradon(r,theta,'linear','Shepp-Logan');