TY - JOUR
T1 - Farey Sequences and Discrete Radon Transform Projection Angles
AU - Svalbe, Imants
AU - Kingston, Andrew
PY - 2003/3
Y1 - 2003/3
N2 - This paper examines how the minimal set of digital projection angles for the Discrete Radon Transform (DRT) is selected from the known sequence of Farey fractions. A square array of prime size p defines a unique direction for each digital projection, m, through an integer ratio xm/ym. Here xm and ym define the nearest neighbour distance, dm, between projection samples under a modulus p sampling rule. We show the maximum gap length, dmax, on square and hexagonal lattices is < √(2p/√3) and ≤ √p respectively. The DRT angles are shown to replicate the entire Farey fraction sequence over the interval 1 ≤ dm ≤ √p for a square lattice. For the interval √p < dm < dmax, the DRT set skips some Farey angles. The skipping of Farey angles is the result of a geometric restriction on the distance minimisation process. The DRT angle set varies significantly with p for those projections with dm near dmax. This complicates the comparison of DRT projections over similar sized arrays. The detailed angle distributions of the DRT and the Farey sequences reflect, in different ways, the variable gaps between prime numbers. Thanks are due to Graham Farr, Computer Science, Monash University, for pointing us to the literature on finding minimal vectors in lattices. IS acknowledges support from the Centre for X-Ray Physics and Imaging for this project. AK is a Monash University postgraduate student in receipt of an Australian Postgraduate Award scholarship, provided through the Australian Government.
AB - This paper examines how the minimal set of digital projection angles for the Discrete Radon Transform (DRT) is selected from the known sequence of Farey fractions. A square array of prime size p defines a unique direction for each digital projection, m, through an integer ratio xm/ym. Here xm and ym define the nearest neighbour distance, dm, between projection samples under a modulus p sampling rule. We show the maximum gap length, dmax, on square and hexagonal lattices is < √(2p/√3) and ≤ √p respectively. The DRT angles are shown to replicate the entire Farey fraction sequence over the interval 1 ≤ dm ≤ √p for a square lattice. For the interval √p < dm < dmax, the DRT set skips some Farey angles. The skipping of Farey angles is the result of a geometric restriction on the distance minimisation process. The DRT angle set varies significantly with p for those projections with dm near dmax. This complicates the comparison of DRT projections over similar sized arrays. The detailed angle distributions of the DRT and the Farey sequences reflect, in different ways, the variable gaps between prime numbers. Thanks are due to Graham Farr, Computer Science, Monash University, for pointing us to the literature on finding minimal vectors in lattices. IS acknowledges support from the Centre for X-Ray Physics and Imaging for this project. AK is a Monash University postgraduate student in receipt of an Australian Postgraduate Award scholarship, provided through the Australian Government.
KW - Farey sequences
KW - digital projection
KW - discrete Radon transform
UR - http://www.scopus.com/inward/record.url?scp=34247276118&partnerID=8YFLogxK
U2 - 10.1016/S1571-0653(04)00482-2
DO - 10.1016/S1571-0653(04)00482-2
M3 - Article
SN - 1571-0653
VL - 12
SP - 154
EP - 165
JO - Electronic Notes in Discrete Mathematics
JF - Electronic Notes in Discrete Mathematics
ER -