<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE document PUBLIC "-//CNX//DTD CNXML 0.5//EN" "http://cnx.rice.edu/technology/cnxml/schema/dtd/0.5/cnxml_plain.dtd">
<document xmlns="http://cnx.rice.edu/cnxml" xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:bib="http://bibtexml.sf.net/" id="id28760614">
<name>Reconstruction</name>
<metadata>
  <md:version>1.1</md:version>
  <md:created>2005/12/12 00:05:35.932 US/Central</md:created>
  <md:revised>2005/12/12 20:17:09.585 US/Central</md:revised>
  <md:authorlist>
      <md:author id="dmiller7">
      <md:firstname>Deborah</md:firstname>
      <md:othername>Marie</md:othername>
      <md:surname>Miller</md:surname>
      <md:email>dmiller7@rice.edu</md:email>
    </md:author>
      <md:author id="wscott">
      <md:firstname>Warren</md:firstname>
      <md:othername>R</md:othername>
      <md:surname>Scott</md:surname>
      <md:email>wscott@rice.edu</md:email>
    </md:author>
  </md:authorlist>

  <md:maintainerlist>
    <md:maintainer id="dmiller7">
      <md:firstname>Deborah</md:firstname>
      <md:othername>Marie</md:othername>
      <md:surname>Miller</md:surname>
      <md:email>dmiller7@rice.edu</md:email>
    </md:maintainer>
    <md:maintainer id="wscott">
      <md:firstname>Warren</md:firstname>
      <md:othername>R</md:othername>
      <md:surname>Scott</md:surname>
      <md:email>wscott@rice.edu</md:email>
    </md:maintainer>
  </md:maintainerlist>
  
  

  <md:abstract/>
</metadata>
<content>

<section id="fbp">
<name>Theory of Filtered Backprojection Algorithm (FBP)</name>
<para id="id27778402">The FBP algorithm allows us to take the
projections, PӨ(t), developed in the previous sections and
reconstruct the original image, f(x,y).</para>
<para id="id28778068">A key idea is the Fourier Slice Theorem. It
says that the Fourier Transform of a projection at an angle theta
is equivalent to the values of the 2-dimensional Fourier Transform
of the image evaluated along a radial line of the same angle.
Knowing this fact, we are able to obtain the Fourier Transform of
the image, F(u,v), from projections taken at multiple
angles.</para>
<para id="id3027203">We start with 2-dimensional Inverse Fourier
Transform:</para>
<para id="id2664234"><figure id="eq123">
<media type="image/jpeg" src="rec_eq1.gif"/>
</figure></para>
<para id="id27954060">Since we have projections for given angles, a
change to polar coordinates is useful.</para>
<para id="id28770907"><figure id="eq456">
<media type="image/jpeg" src="rec_eq2.gif"/>
</figure></para>
<para id="id3643267">Using symmetry, this simplifies to:</para>
<para id="id2"><figure id="eq7">
<media type="image/jpeg" src="rec_eq31.gif"/>
</figure></para>
<para id="id28316353">Using the Fourier Slice Theorem, we
substitute in the Fourier Transform of the projection,
SӨ(ω).</para>
<para id="id14741629"><figure id="eq4441">
<media type="image/jpeg" src="rec_eq4.gif"/>
</figure></para>
<para id="id28773823">With this formula, we are now able to
reconstruct the original image. We now see that that the FPB
algorithm has certain benefits. We can begin reconstructing the
image after the first projection has been calculated, since the
image is built up by summing over all the angles. This could
increase speed and practicality for real time applications.</para>
</section>
</content>
</document>
