18 - Diagnostic Medical Image Processing (DMIP) 2010/11 [ID:1267]
50 von 979 angezeigt

Okay, so please apologize.

So we start.

Good morning for the Tuesday session.

Today, 60 minutes on image reconstruction.

And yesterday, I was at the beginning a little bit annoyed about the fact that the slides

were corrupted and I couldn't read it.

But I think we got a very good intuition what is actually going on in CT reconstruction

or to say it in other words, what are the basic problems that are required to be solved

to do a CT reconstruction.

So we have to know about the line integral, the projection ray.

We know how to characterize the 3D projection ray given a point source.

Projection matrices is basically what we need for that.

We have discussed these things in detail.

We have seen that we basically have to solve a system of integral equations or in other

words, we measure a set of line integrals of a function and the problem is to compute

the original function.

That's CT reconstruction.

And then we discussed the Fourier slice theorem.

We did not prove it.

We discussed it.

And we have referenced, let me say we have referenced a very nice result saying that

the 2D Fourier transform of the slice we want to reconstruct has a relationship to the projections

and their Fourier transform.

And basically the picture we have to keep in mind is this one here.

We have our slice.

We want to reconstruct.

We have our, let me draw it this way, our detector and then we have our projection rays.

And if I consider the Fourier transform of this two dimensional slice I want to reconstruct,

I know by the Fourier slice theorem that basically the Fourier transform of this 1D signal, what

we measure is identical to the Fourier transform, the 2D Fourier transform sampled along the

line that is parallel to the detector.

And what we started out to discuss yesterday at the end of the lecture, we looked at the

very specific case that the y coordinate is zero of the slice we want to reconstruct or

let me say it this way, the rotation angle of our detector is zero.

And we were able to see that the Fourier transform of this 1D signal is identical to this one.

The final argument is still missing.

We can have a look at it closely in a few minutes.

And the second stage of the proof of the Fourier slice theorem is basically the fact that it

makes no difference whether I rotate the function and compute the Fourier transform of it or

whether I compute the Fourier transform of a function and rotate it afterwards.

That means if I have a 2D function and I rotate the function and then I compute the Fourier

transform, the result is the same as if I compute the Fourier transform of the original

function and then rotate the Fourier transform.

So this commutes.

And what we basically can do then is if we know this relationship here between this theta

angle is zero and this line here, we basically rotate back the Fourier transform and consider

always this straight line here and successively sample the Fourier space by rotating the Fourier

transform or rotating the detector around it.

That's the argument that we apply.

So this is the picture.

Zugänglich über

Offener Zugang

Dauer

00:56:39 Min

Aufnahmedatum

2010-12-21

Hochgeladen am

2011-04-11 13:53:29

Sprache

de-DE

Einbetten
Wordpress FAU Plugin
iFrame
Teilen