Homework Spring 2007: Set 2

In this homework set we are creating a system matrix for 2-D imaging. You can write your own set-up routine, or somehow make use of the MATLAB Image toolbox function "radon".

To test your programs apply them to the Shepp-Logan phantom (MATLAB image toolbox function "Phantom"), or use any one of the following (grey-scale) sample images:

  1. Find a way to compute the system matrix (directly related to the Radon transform), and a way to efficiently store this matrix or easily recompute its entries.
  2. From your chosen sample image, create the projections according to your system matrix. COmpute the "exact" values, as well as some noisy data.
  3. Implement the straightforward ART algorithm, and reconstruct the image.
Keep the program for the system matrix. You will be using it in your "project".