 |
|
Verwenden Sie die folgenden Originalparameter
eines Siemens Computertomographen für Ihre
Programme:
| phi_min |
= |
249° |
| phi_max |
= |
249+289*180/290° |
| P_min |
= |
-258.97245 mm |
| P_max |
= |
258.20055 mm |
Der Bildbereich sei 300 mm x 300 mm groß.
Wählen Sie je nach Leistungsfähigkeit Ihres Rechners
die Parameter n (Feinheit des Gitters), nPhi (Anzahl der
gleichverteilten Drehwinkel zwischen phi_min und phi_max), nP
(Anzahl der gleichverteilten Röntgenquellen zwischen P_min
und P_max) aus der folgenden Tabelle. Versuchen Sie es
zunächst einmal mit den Parametern aus der untersten Zeile.
| n |
nPhi |
nP |
| 256 |
290 |
336 |
| 128 |
290 |
336 |
| 64 |
290 |
336 |
| 64 |
145 |
168 |
| 32 |
145 |
168 |
| 32 |
72 |
84 |
Die Zuordnung der Strahlen zur rechten Seite b sei gemäß
der folgenden Tabelle:
| (phi_min,P_min) |
... |
(phi_min,P_max) |
... |
(phi_max,P_max) |
| b_1 |
... |
b_nP |
... |
b_(nP*nPhi) |
Testen Sie Ihre Programme an eigenen Phantoms sowie an den
folgenden Daten:
|
|
Rekonstruieren Sie damit diese Bilder:
(n,nPhi,nP)=(32,72,84)
|
(n,nPhi,nP)=(64,145,168)
|
(n,nPhi,nP)=(128,290,336)
|
(n,nPhi,nP)=(256,290,336)
|
Verwenden Sie dazu Matlab z.B. die Befehle imagesc und
colormap(gray).
Falls Sie noch kein lauffähiges Programm zur Berechnung der Matrix A
besitzen, benutzen Sie die
function out=erstelleA(n,nPhi,nP).
|
