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:
-
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.
-
From your chosen sample image, create the projections according
to your system matrix. COmpute the "exact" values, as well
as some noisy data.
-
Implement the straightforward ART algorithm, and
reconstruct the image.
Keep the program for the system matrix. You will be using it
in your "project".