Warum gibt iradon negative Pixelwerte zurück?

Question

Ich nahm 200 Projektionen bei einem Schrittwinkel von 1,8 Grad mit LabVIEW-Software. Die Größe des Bildes ist 2748 x 2748 Pixel, Uint16. Dann lade ich mit Matlab die Projektionsbilder ein, mache die Flachfeldkorrektur, verändere das Bild um 1/3 und speichere die Bilder als .mat-Datei. Dann führe ich den folgenden Code für die gefilterte Rückprojektion aus.

interp='linear'; %set interpolation: nearest, linear, spline, pchip, v5cubic 
filter='Hann'; %set filter: Ram-Lak, Shepp-Logan, Cosine, Hamming, Hann, None 

for s=1:916 

    for i=1:200 
    a(i,:)=proj065(:,s,i); 
    end 
    a=a'; 
    %figure(3), imagesc(a) 
    b=iradon(a,1.8,interp,filter); 
    imagesc(b); 
    recon(:,:,s)=b; 
    s 
    clear a 
end 

Wenn ich einen Filter in diesem Code verwende, werde ich negative Pixelwerte erhalten.

Aber wenn ich den Code ohne den Filter ausführen, werde ich positive Pixelwerte erhalten.

Gibt es eine Idee, warum iradon negative Pixelwerte in gefilterter Rückprojektion liefert?

Vielen Dank.

Nurul

Answer 1

Ja, die FBP (gefilterten Rückprojektion) Algorithmus wird das tun. Falsche Rekonstruktion von Voxeln mit negativen Werten aufgrund von Rauschen und Diskretisierung der Daten. Nichts, was Sie dagegen tun können, als nur diese Werte zuzuschneiden.

Da meine Doktorarbeit über Tomographie-Rekonstruktionsalgorithmen ist, fühle ich mich vertraglich verpflichtet (scherzhaft), die Verwendung von iterativen Algorithmen zu empfehlen, um bessere Bilder zu erhalten (nie schlechter, oft wesentlich besser). Überprüfen Sie SART/SIRT oder CGLS für dieses Problem.

Sie rufen jedoch Ihre Funktion falsch an! In der Tomographie reicht die Schrittweite nicht aus, um ein Bild zu rekonstruieren, im Allgemeinen brauchen Sie die exakten Winkel, also akzeptiert iradon keine Schrittweite als Eingabe, es akzeptiert eine Anordnung von Winkeln.

in Ihrem Fall sollte thetatheta=linspace(0,360-200/360,200) sein, und Sie sollten für Ihren Vorschlag iradon(a,theta,...)