Rapid analysis of non-uniformly sampled pulsed field gradient data for velocity estimation

Bretthorst's recent generalization of the Lomb-Scargle periodogram shows th at a sufficient statistic for frequency estimation from non-uniformly, but simultaneously sampled quadrature data is equivalent to the FFT of those da ta with the missing samples replaced by zeros. We have applied this concept to the rapid analysis of pulsed field gradient MRI data which have been no n-uniformly sampled in the velocity encoding wave vector q. For a small num ber of q samples, it is more computationally efficient to calculate the per iodogram directly rather than using the FFT algorithm with a large number o f zeros. The algorithm we have implemented for finding the peak of the gene ralized periodogram is simple and robust; it involves repeated apodization and grid searching of the periodogram until the desired velocity resolution is achieved. The final estimate is refined by quadratic interpolation. We have tested the method for fully developed Poiseuille flow of a Newtonian f luid and have demonstrated substantial improvement in the precision of velo city measurement achievable in a fixed acquisition time with non-uniform sa mpling. The method is readily extendible to multidimensional data. Analysis of a 256 by 256 pixel image with 8 q samples and an effective velocity res olution of better than 1/680 of the Nyquist range requires approximately I minute computation time on a 400 MHz SUN Ultrasparc II processor. (C) 2001 Elsevier Science Inc. All rights reserved.