TheFourierTransformAndItsApplications-Lecture29

Instructor (Brad Osgood):Does somebody know Ė donít you occasionally get bad reactions to flu shots like about a week later? I got a flu shot last week, and I sort of had all these flu symptoms yesterday, like I thought I was gonna die for example. So anyway, I wonít be my usual perky self today, and of course, only half the class is here, so I guess perkiness is in short supply these days. All right. I wish we had one additional day. Well first of all, let me call your attention one more time to the information about the final exam. Itís next Thursday. Thatís a week from tomorrow from 8:30 to 11:30 in Dinkelspiel Auditorium. Iím sure you can find that. It is open books, open notes, and Iíll supply the blue books, and Iíll bring copies of the formula sheets, and so on, and so on. Does anybody have any questions about the exam, about the mechanics of it or anything like that? No? Everybody brimming with confidence, or about shocked and doesnít really wanna think about it right now? Okay. Like I say, I wish we had one additional day in this quarter because I [inaudible] donít quite have enough time to do what I wanna do in full detail, or even if not in full detail, just enough detail so you can see some of the finer points of the arguments. But nevertheless, I still think itís important to try to go over the main highlights of two more topics: 1.) a continuation of what we talked about last time, the general stretch theorem and actually, sort of an application and a physical manifestation of the general stretch; and then the last topic which Iíll begin today and then finish up next time on medical imaging and application of the two-dimensional Fourier transform to medical imaging.

And again, Iím sorry. I really do regret not having one more day somehow that we could do this in a little bit more detail. But to give you a little Ė at least a sense of some of the directions that some of this material goes. So today I wanna talk about Ė and actually, this hearkens back both to yesterday, last time, and stuff we did earlier when we first started talking about the Shah function. I wanna talk about Shahs, lattices, and crystallography. Now Iím doing this because I really think this is a really interesting application. Itís very pretty, and it just shows you somehow how some of these ideas come into play in unexpected ways. This is an application Ė although again unfortunately, we canít see all of the details. I just donít time to give all the details, although itís discussed in the notes. This is an application and say physical manifestation of this idea that in higher dimensions in reciprocal means inverse transpose. Thatís what we talked about when we talked about the general stretch theorem, and Iíll remind you what that says. Itís a physical manifestation Ė physical example. Manifestation is too long a word for a man with a grip Ė physical example of the idea is reciprocal means inverse transpose. We saw that phrase, that aphorism in the context of the generalized stretch theorem, so Iíll remind you what that says. And that comes into the derivation of what Iím going to be doing. So we saw this in the generalized stretch theorem which told us how to take the Fourier transform one and make it change a variable by a matrix, and not just scale the variables independently.

So it said this: it said the Fourier transform of F of AX Ė all right, so you change the variables X by a matrix A, a nonsingular matrix A Ė is one over the determinate of A times the Fourier transform of F evaluated at A inverse transpose at the frequency variable C. Okay? Itís a very interesting formula. We derived it last time, and itís complicated. Itís more complicated than the one-dimensional stretch case, but it includes the one-dimensional stretch case, but what you donít see in one dimensions is this new phenomenon as I say that reciprocal somehow means inverse transpose. Thatís something new. Thatís something different on the scene. All right, so now I wanna show you how this comes in in these contexts. So first of all, I wanna talk about how we might generalize the Shah function to higher dimensions, so let me remind you of the Shah function in one dimension, and some of its main important properties. So I want to generalize the Shah distribution to the sum of deltas. Let me say evenly spaced deltas. So in one dimension we had the model case. That was a Shah function where the deltas are spaced by one, so you put a delta at each integer point. The model case is the basic Shah function where you put a delta function at each integer, zero, one, two, minus one, minus two, and so on. And you define the Shah as just the sum of these delta functions. Shah of X is the sum of K going from minus infinity to infinity at delta X minus K, so thereís a delta function at each integer point.

And the remarkable property Ė and as we said when we were talking about this, this is the deepest property known about the integers. Its equivalent to this fact in number theory called the Poisson summation formula is the Fourier transform of the Shah function is itself. The Fourier transform of Shah is Shah. And that was fundamental in a lot of the work we did. It was fundamental also into the sampling theorem, and also into other applications. We actually looked at this in a slightly more general version, which Iíll remind you of just now, in the context of crystallography, sort of as a way of motivating it, and Iím gonna talk exactly about that in just a second. So we generalize this slightly. We allow for spacing P for other spacing other than the integers, so spacing P. Doesnít matter what I call it. So instead of space to the integers, itís space Ė I put deltas at the points zero, P, 2P, and so on, and so on, minus P, minus 2P. So this time the Shah function is the same sort of definition but the spacing is different. So again itís the sum from minus infinity to infinity of delta X minus K times P. Here the fundamental result is that when you take the Fourier transform, the spacing becomes reciprocal. That is the Fourier transform of Ė let me use the subscript here P to indicate that thatís what the spacing is. Sorry. The Fourier transform of the Shah function with spacing P is actually Ė thereís two places where the reciprocal comes in. It gets scaled out in front by one over P, and then the spacing of the deltas in the frequency domain Ė the spacing of the deltas when you take the Fourier transform is spaced by one over P.

Thatís pretty cool. Itís a very interesting result, and itís a very important example of this sort of reciprocity relationship between the two domains. And this is what I wanna generalize. All of these properties I wanna generalize. Now, Iím gonna stick in two dimensions, but really Ė and thereís a change from one dimension to two dimensions, but thereís very little change from two dimensions to higher dimensions. I can draw the pictures more easily in two dimensions, but really all I have to say will hold in higher dimensions, and in particular three dimensions, which is really the appropriate setting for studying crystals. But as I say just so I can draw the pictures a little more easily, Iím gonna write things and speak about things just in two dimensions. So whatís the generalization of Shah letís say to 2-D. What should it be? Well, what I wanna do is Ė there are two parts to the Shah function. I put a bunch of delta functions, but before I put the delta functions down, I have to have a set of evenly spaced points. That was the whole idea behind the Shah function, so you need Ė thatís the first thing that you want. You want evenly spaced points, and you have to decide what that means, evenly spaced points in the plane, in R2. Now there are a lot of different possibilities here, but again we started with a Shah function with a model case where the deltas or the points were spaced one apart. So for my model of evenly spaced points in the plane, Iím gonna take all the points in two dimensions whose coordinates are integers. So Iíll take a grid in the plane whose points have integer coordinates like so.

So thatís gonna be my model for evenly spaced points in the plane. And of course, in the higher dimensions it would be a similar sort of thing. I would take points with integer coordinates in R3 or in higher dimensions. So a point here, a grid point has coordinates K1 and K2 where K1 and K2 are integers, all these points. Those points are said to form a lattice. Thatís not an unfamiliar term. So this forms whatís called the integer lattice in the plane. Itís denoted by Z2, sort of a bold-faced letter Z. I canít remember if Iíve used that notation before, Z for the integers, to stand for the integers. Z stands for the German word zahlen, which means number. Itís a universally adopted term in Ė pretty much a universally adopted term in certainly mathematics and more and more so in engineering in application just to denote the integers. So Z2 meaning itís in two dimensions. So this is the set of all points, nothing but the set of all points and the lattice points in the plane with integer coordinates. And my Shah function is gonna be defined by putting a delta at each one of those points. Put a two-dimensional Ė put a spike at each one of those points. So define by Shah by putting a delta at each one of those points, each grid Ė Iíll call it lattice points, okay? That is Iíll define the Shah function for this integer lattice as exactly the analog of the one-dimensional case of X. Itís gonna be the sum Ė let me write it like this: K in Z2 of delta X minus K. So again, what Iím trying to do here is Iím trying to make this look as much like the one-dimensional case as possible. So if I write the vector K, what I mean by that is itís a pair, K1, K2 where K1 and K2 are integers. So if I write it like that again it looks like the one-dimensional case. Here again K is just shorthand notation K1, K2, integers.

The basic phenomenon Ė again, generalizing the one-dimensional case is what happens to the Fourier transform of the Shah function for the integer lattice. And this is unfortunately something I donít feel like I really have the time to derive, and it just kills me, but Iíll give you the result. So just as in the 1-D case, you have with the Fourier transform Ė and for the same reasons really, itís a straightforward Ė you set things up so it looks exactly the same as in the one-dimensional case. You have the Fourier transform, the two-dimensional Fourier transform of this lattice Ė the Shah function based on the integers is itself. The proof really again, the derivation of that is the same as in the one-dimensional case. It depends on Ė actually, it depends on higher dimensional Fourier series, but itís really very much the same. You can use the same words and really follow the same argument. The argument is very much the same, very much like the 1-D case. Now, things get more interesting when you ask yourself, ďHow do I change the spacing?Ē In the line, in one dimension, you only have one degree of freedom. The only way to change the spacing is to stretch or shrink in one direction. So how do I change the spacing for a 2-D lattice? What should we mean by that? We still want the points to be somehow evenly spaced, but not the same. Well again, in two dimensions, we have more degrees of freedom than in one dimension, but for our purposes, thereís a very natural way of looking at this. In two dimensions, we obtain a different lattice. Informally, Iíll call it an oblique lattice. Iíll draw a picture in just a second. That is we change the spacing by applying a linear Ė a two by two matrix say A to the integer lattice.

Thatís what it is in words. Let me draw a picture so you see what I mean here because itís not hard to see what I mean. So here let me draw a picture of the integer lattice over here. So this is Z2. And as a matter of fact, let me take Ė letís consider this the origin, and letís take this to be the natural basis of R2, so the two natural basis vectors one zero and zero one. Now I take a linear transformation, a two by two matrix A nonsingular, so nothingís getting collapsed. Then theyíre gonna take those basis vectors E1 and E2 to some other vector, say V1 and V2. And theyíre gonna take Ė now all the points in the integer lattice are just integer combinations of these basis vectors. If I wanted to reach the Point 2 over here, then I just take two times E1. If I wanna reach the Point 2 up here, I take two times E2. If I wanna reach this point, I take E1 plus E2. Okay? One times E1 plus one times E2, and so on and so on. So all the points in the integer lattice are just integer combinations of the basis vectors E1 and E2. Under a matrix or linear transformation, they go over to integer linear combinations of V1 and V2, so what I have over here is a picture that looks like this where the lattice points over here correspond exactly to the lattice points over here, and so on and so on.

So the lattice is sheared if you wanna think of it that way. The square lattice of the integers is made into an oblique lattice by taking a matrix that takes the basis vectors over here, the natural basis vectors and some other pair of vectors. So A of Z2 Ė A applied to Ė I just mean by this applying A to all the different basis Ė all the different points in Z2 results in an oblique lattice L. Okay? And L is just the combination Ė the points in L are combinations of the basis vectors K1 times V1 plus K2 times V2 where K1 and K2 are integers. And here V1 say is A of E1 and V2 is A of E2. Stay with me. Now, one more geometric Ė actually, there are a couple of other things I have to introduce. Itís very common Ė as a matter of fact, crystallographers do this Ė do we have anybody who has actually done any crystallography in here? By the way, Iím just curious. Then this is all new.

But if you do any kind of imaging, youíre gonna see these terms and see these ideas, especially these days in molecular imaging and things like that where these things are coming in in really interesting ways. As a quick geometric description Ė so to jump ahead a little bit, and Iím gonna come back to this in a second, crystals come into this because you consider crystals to be modeled on lattices where thereís an atom at each lattice point. And of course, the idea of considering oblique lattices is because some crystals are modeled on a rectangular lattice or a cubic lattice, but a lot of crystals arenít. Thereís sheer. Theyíre oblique. Thereís a geometric Ė just to get an idea of somehow just to attach a number that describes a lattice in some loose but helpful sense, you often talk about the area of a lattice. So you often speak of the area of really whatís called a fundamental cell, or fundamental parallelogram, fundamental cell in the lattice. And thatís just the area of one of the parallelograms spanned by the basis vectors. So again, if hereís the lattice, and say this is V1, V1 at V2, then the area of the lattice is the area of the Ė weíre gonna come back Ė weíre gonna need this concept in just a second. You talk about the area of the lattice as just the area of the parallelogram determined by the basis vectors, so the area of L is area of the parallelogram.

The reason why crystallographers think about that is because they think of that as the fundamental cell or fundamental building block for the rest of the crystal. That region in space is the fundamental building block for the rest of the crystal, and the question is how big is that. They talk about the volume of the fundamental cell. Everything in crystallography of course is in three dimensions. Iím only doing things in two dimensions, so Iím talking about the area in terms of volume, but thatís what they do. They talk in those terms. So the area of the integer lattice is one because each one of these basis vectors is the length of one, and theyíre perpendicular. The area of that square is just one. Now I wanted to introduce that, but the real thing I want to get to of course is what happens if we consider the Shah function for an oblique lattice, and what is its Fourier transform. So again considering the Shah function for oblique lattice is like considering a model for a crystal where the underlying lattice structure or underlying crystal structure is on an oblique lattice instead of on a square lattice, so we can consider Ė letís call it the Shah function of lattice L. Iím struggling here with my script letters.

So let me just write it like this, itís the sum over all points in the lattice, so those are all the lattice points, the grid points Ė I put a delta at each one of those points, say delta of X minus P. You put a delta function at each lattice point P. Now the real question is what is the Fourier transform. What happens to the Fourier transform? What is the Fourier transform of the Shah function associated with the oblique lattice? A lotís going on here, right? You have to have the notion of a lattice. You have to have the notion of an oblique lattice. You have to have the notion of a Shah function on the lattice. And now youíre asking about the Fourier transform. I mean what levels of complexity have been added to this simple conversation. And the reason why Ė and again, Iíll come back to this in a second, but the reason why you wanna know about the Fourier transform is because in honest to God scientific experiments on crystallography, when youíre shining X-rays through a crystal, you get a bunch of spots, and the pattern of the spots is determined by the Fourier transform of the crystal thatís doing the diffracting.

So if the crystal thatís doing the diffracting is modeled by the Shah function for a lattice, what you know about is you wanna know about the Fourier transform of that Shah function for the lattice because thatís what youíre seeing. Those are the spots that youíre seeing. So now finally, hereís where this reciprocal relationship comes in and itís just great. And again unfortunately, and I really do feel badly about this, I donít have the time to talk too much about the motivation for this, or even derive many of the formulas, so I just wanna show you in some sense what the punch line is without the whole set up, and it just kills me, but Iíve got to. So Iím gonna put two other lattices on that board, but hereís the integer lattice Z2. I have two linear transformations. I have my linear transformation A thatís going to go down to the oblique lattice L. Thatís the sort of the lattice that Iím starting with. Hereís the oblique lattice L. And then the crystallographers Ė and even mathematicians, too, but for completely different reasons as far as I can tell Ė define what they call the dual or reciprocal lattice, L star. So L is A of Z2. L is the transformation A. Any lattice can be realized as a linear transformation applied to the integer lattice. So L is A applied to Z2. Then the crystallographers the reciprocal or dual lattice to be Ė guess what Ė A inverse transpose applied to Z2. So they take A inverse transpose Ė [Crosstalk]

Now thatís also another linear transformation. It goes into another lattice, and I donít have great confidence in actually how I draw it, but something like this say, L star. And they call that the reciprocal lattice. And in fact, as you can easily check, the area of L star Ė that is to say the area of one of these fundamental polygons, polyhedral is Ė so I certainly didnít draw it very well because it doesnít look quite so easy Ė is the reciprocal of the area of the other lattice. Now this is one definition of it, and this is the cleanest definition. Itís completely unmotivated. Iím sorry. There are different ways of motivating it. Thereís one discussion in the notes. That is itís a completely unmotivated definition. Why the hell would you define a dual lattice or reciprocal lattice this way? You just have to take my word for it or look around for other motivations for it, other definitions for it. In fact, people who work in materials and people who work in crystallography actually have a geometric construction of dual lattice. Given a lattice L, they have a way of constructing geometrically the dual lattice, like you know a ruler and compass construction. So itís a time honored important construction in materials to pass from the lattice to the reciprocal lattice.

Now, why? What is the big deal? Because the fundamental on the Shah function is that the Fourier transform for the Shah function for a lattice is the Shah function for the dual lattice except scaled by the area, so thatís a fundamental fact. The Fourier transform Ė Iíve got too many Fs and Ls and everything else in here. The Fourier transform of the Shah function of a lattice L is Ė itís a beautiful formula Ė itís one over the area of L times the Shah function of the dual lattice. Now again in higher dimensions which is where the real applications of crystallography go, youíd have the same result except youíd have volume here instead of area, but itís the same kind of result, and the definition of a dual lattice is the same, except instead of a transformation in R2, instead of everything happening in R2 here, everything would be happening in R3. Youíd have the lattice of integer points in R3. Youíd have a three by three matrix going over to an oblique lattice here in R3. Youíd have the inverse transpose going into another oblique lattice of R3 and so on. And youíd have the same result. This is exactly Ė if you believe it Ė analogous to the one-dimensional formula. This is the analog of the one-dimensional formula. This is the analog to the Fourier transform of the one-dimensional Shah function of with spacing P is one over P times the Shah function with spacing one over P.

There reciprocal just means one over. But see, the thing is you miss something. If you just look at that, you either miss something in the one-dimensional case, or something new and deeper is revealed in the higher dimensional case because in the one-dimensional case, reciprocal means reciprocal. It just means one over. In the higher dimensional case, reciprocal means inverse transpose. Thatís the change in point of view. Thatís sort of the change in intuition. Itís very important. Thatís something you donít see in the one-dimensional case. Thatís a difference between the one-dimensional case and the two-dimensional case in how to understand the notion of reciprocity, to understand what reciprocal means. Itís a mathematical fact, and not only a mathematical fact, a fact of nature that in higher dimensions reciprocal often means Ė can be interpreted as involves, ergo it [inaudible] reciprocal transpose. Now so what about crystals in all of this? I made allusions to this a lot, so let me just say quickly what the situation is for crystals and why this is so important or where this comes up. So again for crystals Ė so again when we did this in the one-dimensional case, the idea is you study a crystal by studying the electron density distribution of a crystal. What is the electron density distribution of a crystal? You see, you conjecture, you measure how the electrons are distributed about an atom, and then you periodize that. Crystals have a periodic structure, so row of X in two dimensions now is the electron density for an atom in a crystal. The electron density distribution for the whole crystal is a periodized version of that. For the crystal as a whole, you take a periodized version. Thatís the whole point about crystals is they have periodic structure. How do you take a periodized version? You convolve with the appropriate Shah function.

So if this is the crystal whose atoms are what you model as a bunch of atoms at lattice points as L, then the density for the lattice, the density for the crystal is the density for a single atom convolved with the Shah function. Thatís sum over all the points in the lattice of row of X minus P, same as in the one-dimensional case. So if thereís a little Ė if the density looks like that there, then that patternís just repeated for the whole lattice, and again I canít draw this too well, but you get the picture. This is now periodic with respect to the lattice. This is also a different phenomenon. It has two periods. Itís periodic in that direction, and itís periodic in that direction. Itís periodic on the whole lattice. The pattern repeats on the whole lattice. Now X-ray crystal experiments, X-ray diffraction experiments measure the Fourier transform. That is an X-ray diffraction experiment produces the magnitude of the Fourier transform, but just basically think of it as the Fourier transform Ė produces the Fourier transform of row of the lattice. Thatís what it does. You see a bunch of spots. What youíre seeing is the Fourier transform of the lattice, and what is that? Well, thatís the Fourier transform of the convolution, so thatís the product of the Fourier transform of row with the Fourier transform of the Shah. This is the Fourier transform of row times the Fourier transform of the Shah function of the lattice.

But we know what that is. Itís the dual lattice. That is this is the Fourier transform of row times one over the area of L times the Shah function of the dual lattice. And if I take the product of the function times a bunch of delta functions, what do I get? I get sum over the points Ė let me call it P star in the dual lattice. Theyíre just delta functions at the lattice points of the dual lattice of the Fourier transform of row at Ė I donít know, S times delta of S minus P star. So the Fourier transform of the density function times the sum of deltas, and remember what happens when you take a function times delta? It just pulls out the value Ė sorry, P star. So once again, you do your experiment, and unless you know the math, youíre not gonna be able to draw the right conclusion. You might think that you do Ė you pass a bunch of X-rays through the crystal, and the spacing of the lattice points that youíre seeing should be proportional to the spacing of the lattice points in the crystal, but itís not. Nature for whatever mysterious reason is taking a Fourier transform. X-ray crystallography, X-ray diffraction takes a Fourier transform for you physically. And what youíre seeing is youíre see the points at Ė youíre seeing the little spots that appear on your X-ray film at points on the dual lattice. And unless you know the math, unless you know that in higher dimensions reciprocal means inverse transpose and so on, you can kiss your Nobel prize goodbye.

But if you sort of follow your pencil through this, and see how the math plays out here, you can draw conclusions about the crystal by knowing this. I think itís really cool. I think itís a really interesting application and physical manifestation of this fact that in higher dimensions inverse transpose means reciprocal. So I donít have time to derive this formula, and Iím sorry for that. The derivation of this formula Ė and that shouldnít surprise you given how all these terms are coming Ė involves exactly the generalized stretch theorem. The reason why inverse transpose comes in here, and the reason why the dual lattice comes in is because in deriving this formula, it uses exactly the generalized stretch theorem. Thatís how it Ė it comes into it in deriving that formula. Isnít that nice? So I say I wish we had one additional day because I would like to go into just a little bit more detail about this, but we donít. We donít because I spent all that time on the damn fast Fourier transform. Everybody wanted to see that damn fast Fourier transform. So full stop, speaking about things I canít spend enough time on, thereís one more topic I wonít be able to spend enough time on. The final topic in the class is application of the two-dimensional Fourier transform to the problems of medical imaging, and particularly to the problems of tomography. So weíll leave this fascinating subject behind, and pass to another fascinating subject. Iíll start off today, and then weíll finish it up next time.

So the final application is Ė the final topic is the application of the two-dimensional Fourier transform to medical imaging, and in particular the problem of tomography. I will tell you what the problem of tomography is in just a second Ė the fundamental problem of tomography. Now let me say at once that what Iím gonna be looking at is really in the realm of CAT scan, of computer tomography, of tomography not MRI. But as it turns out this is actually Ė this is not MRI. As they say in the biz, the modalities of imaging for the two kinds of Ė two approaches are different. However, itís also true for whatever reason Ė well, for a reason that we canít really get into that the final formula that Iím gonna get to solve the problem of tomography is the same as the formula you get when youíre using an MR system to do the imaging. So what weíre doing here, although itís in the realm of tomography, CAT scans and things like that, actually the final formulas that you get and a lot of the considerations that weíre gonna talk about actually apply in the MR case, magnetic resonance case, NMR case, like nuclear magnetic resonance, but again I cannot get into it. Hereís the basic setup. You have a 2-D region, a two-dimensional region, think of it as a slice of your body filled with goop, blood, bones, organs, you know. There you are, and thereís all that goop inside you. And that goop inside you is described by a function of variable density. You can imagine describing the density of the goop. So describe the goop by a function mu of X1 X2. So this is taking place in the X1 X2 plane, and mu gives the density of the goop at each point.

So if you know mu, you know you, so to speak that is to say you wanna get you by taking measurements. You want to recover mu X1 X2 by X-ray measurements or let me just say by whatever. Thatís what you want. You want to reconstruct the function mu. Now hereís how itís done in tomography. The approach to tomography is the following. Tomography means Ė if I get this right. I forgot to look it up to get the precise Ė I think it means section, so the idea is you are taking sections, and in particular actually what Iím gonna do is Iím gonna take a one-dimensional section. That is Iím gonna pass a line through this region. So the approach via tomography Ė a line meaning Iím gonna pass an X-ray. The approach via tomography is you shoot an X-ray through the region along a line and it comes out the other end. It gets diffracted. It get attenuated by all this goop in the middle. So hereís the X-ray, and thereís all this stuff in here, and you measure what happens to it when it comes out. So you know how Ė you know what it was going in, and you measure what it is coming out. Say it goes this way, so you know the intensity of the X-ray going in, and then it gets attenuated, so you measure the intensity Ė letís call this I not, the initial intensity. You measure the intensity going out on the other end. Letís call that I.

Now then it is not hard to show as a reasonable approximation Ė and again, this is derived in the notes. You can show that thereís a relatively simple Ė well, relatively simple. Weíll see. Thereís a relationship between the intensity going out and the intensity coming in. This is I not going in, I coming out, and itís attenuated Ė it drops exponentially, so itís given by intensity coming out is the intensity going in, and then I say it drops exponentially. How? It drops exponentially according to the integral of mu along the line, sort of like the average density. Mu is a variable density, but the total drop in intensity is given by E to the minus integral Ė so integral of L over mu Ė Iím using a shorthand notation here Ė is the line integral of mu along L. Again, you can sort of think of that as the average density or whatever. There are different ways of making that argument, and again I donít wanna take the time now to derive that formula for you, but itís actually not that hard to derive, and itís not even that hard to believe.

Now again, you know I not. You measure I. So all those numbers are known. What you donít know is you donít know mu, and in particular you donít know the integral of mu along a line. What you do is you make lots of measurements. You send X-rays through along all different lines, or along whole families of lines. You send X-rays through along many lines, this way, that way. Usually you think of a parallel line, but then you change the angle. You send them in this way. Donít think of this as a lattice. Iím just changing my X-rays here. You send a bunch of lines in. Then you know, and at every time you know the intensity going in, you measure the intensity coming in, so that means you know the integral along all those different lines of mu. Once again, you write this down. I equals I not E to the minus integral of mu over the line L for all different L. You know this. You measure this. So then you know all the numbers integral over L of mu. When I say all the numbers, you know them along all these different lines. Itís because you made those measurements.

And the fundamental problem of tomography as it is often stated is can you recover mu by knowing all these line integrals. Thatís the problem. So again, I wanna make sure you understand here at least what you know and what youíre trying to figure out. You know the value of these integrals because you know this formula, and you know the intensity going in, and you measure the intensity coming out, so that tells you what this. You just solve for it in terms of log. So the question is Ė the fundamental problem of tomography is can you recover mu by knowing the integral of mu along L for all L, for all lines L. Thatís the question. And the answer is affirmative. The answer is in the affirmative. The answer is yes, and it involves the two-dimensional Fourier transform. Thatís how it comes into the picture. Iím not gonna get to that today, but Iíll tell you how weíre gonna approach it actually. So the answer is yes. Itís not obvious, and itís certainly not obvious that the Fourier transformís gonna come into this in any sort of way, but hereís the approach. This is the approach that weíre gonna take, and hereís how it works. And as I say, itís not obvious. It was a brilliant idea. You consider these numbers as defining a transform of mu. What I mean by that is you start with a line Ė whatís that a function of? It depends on the line. The number depends on the line.

So starting with a line L, you compute the number, the integral of mu along L. L goes to Ė thatís a line Ė goes to the integral over L of mu. Thatís a number. When I say transform is you wanna consider this a transform of mu evaluated on the line L. Now thatís sort of mind expanding for you, but thatís the way you look at it. That is write this as R of mu, the transform of mu at a particular line, so again that is by definition that integral. Thatís the integral along L of mu. Thatís called the Radon transform of mu. Weíre almost out of [inaudible]. Itís called the Radon transform of mu. The fundamental problem of tomography can be stated as saying can you invert the Radon transform. R of mu gives you something [inaudible] L. The question is can you find mu given all the values of R of mu. Can we invert? Is there an inverse Radon transform to find mu? And the answer is yes, and the answer is in terms of the two-dimensional Fourier transform. Itís amazing. Itís just amazing, but thatís whatís gonna happen, so thatís why the two-dimensional Fourier transform, and thatís why Fourier transform techniques come into medical imaging because of exactly this connection.

So we will finish up with that next time, and me, Iím gonna go take a pill and collapse. Thank you all. Iíll see you all on Friday.

[End of Audio]

Duration: 51 minutes