Reduction of aliasing in modulation transfer function measurements

The images formed by many radiological systems are difficult to sample at spatial intervals small enough to avoid aliasing in the calculation of the system's modulation transfer function. However, if a system's response can be assumed to be symmetrical, this assumption can be used to effectively double the sampling density and to double the frequency limit before aliasing occurs. To accomplish this, a more complex algorithm is required. In this work, the formula for the calculation of the modulation transfer function from a symmetrical, one-dimensional line spread function is derived and a similar result for a symmetrical, two-dimensional point spread function is presented. The effect of noisy data and errors in the estimation of the offset of the center of the line spread function from a sampling point are investigated by simulation studies. For low noise (relative standard deviation of 1%) and an offset error of no more than 2% or 3% of a sampling interval, reasonable precision is obtained. These conditions appear to be achievable, especially when the noise is Poisson distributed.