Attenuation Data Sets for Radiation Treatment Planning Simulation

You are free to use materials from all parts of this site for the educational purposes for which it was created. You are not free to take anything. It is your responsibility to give credit to this author when you use it. If you find this material useful, it would be courteous to let me know. If you find errors or omissions, it would be responsible to let me know. (Statement taken and modified from http://carbon.cudenver.edu/~hgreenbe/morality.html)

K. A. Bachman
University of Colorado at Denver and Health Sciences Center
Mathematics Department
kbachman@math.cudenver.edu
http://www-math.cudenver.edu/~kbachman/

MAIN PAGE

The links, Head Series and Prostate Series, are to sets of attenuation matrices available for download in zipped format. Each set consists of four resolutions (64x64, 128x128, 256x256, 512x512) with 10, 20, and 40 angles for each resolution. The relevant data in each file is stored in the variable A_set. A_set is a cell array with the following construction using the head examples.

A_set{1}{1}=body
A_set{1}{2}=brain
A_set{1}{3}=eyes
A_set{1}{4}=tumor
A_set{1}{5}=attenuation

The '1' in A_set{1} refers to the slide, or image, number. In all cases, only one slide is used. The second index of the first four cells refers to the n²x1 binary matrix corresponding to the indicated structure, i.e., body, brain, etc. n indicates the length of one side. The fifth cell contains an n²xm attenutation matrix, where m refers to the number of beamlets. In the case of the 64x64 10-angle head example, 9 beamlets per angle were used to completely cover the tumor. Therefore, m=90, so that is the number of columns. There are 4096, n², rows since it is a 64x64 image. For the prostate, the structure is the same, except that there are three critical structures instead of two. Therefore, the sixth cell contains the attenuation matrix.

A_set{1}{1}=body
A_set{1}{2}=rectum
A_set{1}{3}=prostate
A_set{1}{4}=bones
A_set{1}{5}=tumor
A_set{1}{6}=attenuation

In order to assign the matrix of a particular cell, say attenuation matrix for head, to a new matrix M, use the following syntax (MATLAB): M=A_set{1}{5}.

Matrix Nomenclature
Note that in the following links, the matrices are named as follows: Consider A064010. The 'A' indicates that it is an attenuation matrix. The three digits following the 'A' indicate the matrix size. In this case, it is 64x64. The remaining digits indicate the number of angles. In this case, 10 angles were used. In the case of matrices from prostate images, a 'p' is appended to the name, e.g. A064010p.

Compressed Head and Prostate Series
These sets are identical to the above but are stripped of beamlets that miss the tumor and pixels that are not intersected by any beamlet.

Zipped compressed Head and Prostate Series created in MATLAB (Release 13).

The head files contain the following matrices:
A (A_set{1}{5}) - attenuated beams
A_new (compressed A_set{1}{5})
I_comp (List (single-column matrix) of indices of the rows from the original (uncompressed) complete data set that are retained in the compressed set.
colindx (List (single-row matrix) of indices of the columns from the original (uncompressed) attenuation data set that are retained in the compressed set.
b (A_set{1}{1}) - body
b_new (compressed A_set{1}{1})
br (A_set{1}{2}) - brain
br_new (compressed A_set{1}{2})
e (A_set{1}{3}) - eyes
e_new (compressed A_set{1}{3})
t (A_set{1}{4}) - tumor
t_new (compressed A_set{1}{4})

The prostate files contain the following matrices:
A (A_set{1}{^})
A_new (compressed A_set{1}{6}) - attenuated beams
I_comp (List (single-column matrix) of indices of the rows from the original (uncompressed) complete data set that are retained in the compressed set.
colindx (List (single-row matrix) of indices of the columns from the original (uncompressed) attenuation data set that are retained in the compressed set.
b (A_set{1}{1}) - body
b_new (compressed A_set{1}{1})
a (A_set{1}{2}) - rectum
a_new (compressed A_set{1}{2})
p (A_set{1}{3}) - prostate
p_new (compressed A_set{1}{3})
f (A_set{1}{4}) - bones
f_new (compressed A_set{1}{4})
t (A_set{1}{5}) - tumor
t_new (compressed A_set{1}{5})