Quantifying the Statistical Impact of GRAPPA in fcMRI Data with a Real-Valued Isomorphism

The interpolation of missing spatial frequencies through the generalized auto-calibrating partially parallel acquisitions (GRAPPA) parallel magnetic resonance imaging (MRI) model implies a correlation is induced between the acquired and reconstructed frequency measurements. As the parallel image reconstruction algorithms in many medical MRI scanners are based on the GRAPPA model, this study aims to quantify the statistical implications that the GRAPPA model has in functional connectivity studies. The linear mathematical framework derived in the work of Rowe , 2007, is adapted to represent the complex-valued GRAPPA image reconstruction operation in terms of a real-valued isomorphism, and a statistical analysis is performed on the effects that the GRAPPA operation has on reconstructed voxel means and correlations. The interpolation of missing spatial frequencies with the GRAPPA model is shown to result in an artificial correlation induced between voxels in the reconstructed images, and these artificial correlations are shown to reside in the low temporal frequency spectrum commonly associated with functional connectivity. Through a real-valued isomorphism, such as the one outlined in this manuscript, the exact artificial correlations induced by the GRAPPA model are not simply estimated, as they would be with simulations, but are precisely quantified. If these correlations are unaccounted for, they can incur an increase in false positives in functional connectivity studies

Incorporating Relaxivities to More Accurately Reconstruct MR Images

Purpose To develop a mathematical model that incorporates the magnetic resonance relaxivities into the image reconstruction process in a single step. Materials and methods In magnetic resonance imaging, the complex-valued measurements of the acquired signal at each point in frequency space are expressed as a Fourier transformation of the proton spin density weighted by Fourier encoding anomalies: T2⁎, T1, and a phase determined by magnetic field inhomogeneity (∆B) according to the MR signal equation. Such anomalies alter the expected symmetry and the signal strength of the k-space observations, resulting in images distorted by image warping, blurring, and loss in image intensity. Although T1 on tissue relaxation time provides valuable quantitative information on tissue characteristics, the T1 recovery term is typically neglected by assuming a long repetition time. In this study, the linear framework presented in the work of Rowe et al., 2007, and of Nencka et al., 2009 is extended to develop a Fourier reconstruction operation in terms of a real-valued isomorphism that incorporates the effects of T2⁎, ∆B, and T1. This framework provides a way to precisely quantify the statistical properties of the corrected image-space data by offering a linear relationship between the observed frequency space measurements and reconstructed corrected image-space measurements. The model is illustrated both on theoretical data generated by considering T2⁎, T1, and/or ∆B effects, and on experimentally acquired fMRI data by focusing on the incorporation of T1. A comparison is also made between the activation statistics computed from the reconstructed data with and without the incorporation of T1 effects. Result Accounting for T1 effects in image reconstruction is shown to recover image contrast that exists prior to T1 equilibrium. The incorporation of T1 is also shown to induce negligible correlation in reconstructed images and preserve functional activations. Conclusion With the use of the proposed method, the effects of T2⁎ and ∆B can be corrected, and T1 can be incorporated into the time series image-space data during image reconstruction in a single step. Incorporation of T1 provides improved tissue segmentation over the course of time series and therefore can improve the precision of motion correction and image registration

The SENSE-Isomorphism Theoretical Image Voxel Estimation (SENSE-ITIVE) Model for Reconstruction and Observing Statistical Properties of Reconstruction Operators

The acquisition of sub-sampled data from an array of receiver coils has become a common means of reducing data acquisition time in MRI. Of the various techniques used in parallel MRI, SENSitivity Encoding (SENSE) is one of the most common, making use of a complex-valued weighted least squares estimation to unfold the aliased images. It was recently shown in Bruce et al. [Magn. Reson. Imag. 29(2011):1267-1287] that when the SENSE model is represented in terms of a real-valued isomorphism,it assumes a skew-symmetric covariance between receiver coils, as well as an identity covariance structure between voxels. In this manuscript, we show that not only is the skew-symmetric coil covariance unlike that of real data, but the estimated covariance structure between voxels over a time series of experimental data is not an identity matrix. As such, a new model, entitled SENSE-ITIVE, is described with both revised coil and voxel covariance structures. Both the SENSE and SENSE-ITIVE models are represented in terms of real-valued isomorphisms, allowing for a statistical analysis of reconstructed voxel means, variances, and correlations resulting from the use of different coil and voxel covariance structures used in the reconstruction processes to be conducted. It is shown through both theoretical and experimental illustrations that the miss-specification of the coil and voxel covariance structures in the SENSE model results in a lower standard deviation in each voxel of the reconstructed images, and thus an artificial increase in SNR, compared to the standard deviation and SNR of the SENSE-ITIVE model where both the coil and voxel covariances are appropriately accounted for. It is also shown that there are differences in the correlations induced by the reconstruction operations of both models, and consequently there are differences in the correlations estimated throughout the course of reconstructed time series. These differences in correlations could result in meaningful differences in interpretation of results

A scale-model room as a practical teaching experiment

A practical experiment is described which was used to help university students increase their understanding of the effect of construction methods and window design on passive solar heating and electrical heating. A number of one tenth scale model rooms were constructed by students and sited out-of-doors in the late autumn. The models were fabricated to mimic available commercial construction techniques with careful consideration being given to window size and placement for solar access. Each model had a thermostatically controlled electric heating element. The temperatures and electricity use of the models were recorded using data-loggers over a two week period. The performances of the models based on energy consumption and internal temperature were compared with each other and with predictions based upon thermal mass and R-values. Examples of questions used by students to facilitate this process are included. The effect of scaling on thermal properties was analysed using Buckingham&rsquo;s p-theorem.<br /

Weak Gravitational Flexion

Flexion is the significant third-order weak gravitational lensing effect responsible for the weakly skewed and arc-like appearance of lensed galaxies. Here we demonstrate how flexion measurements can be used to measure galaxy halo density profiles and large-scale structure on non-linear scales, via galaxy-galaxy lensing, dark matter mapping and cosmic flexion correlation functions. We describe the origin of gravitational flexion, and discuss its four components, two of which are first described here. We also introduce an efficient complex formalism for all orders of lensing distortion. We proceed to examine the flexion predictions for galaxy-galaxy lensing, examining isothermal sphere and Navarro, Frenk & White (NFW) profiles and both circularly symmetric and elliptical cases. We show that in combination with shear we can precisely measure galaxy masses and NFW halo concentrations. We also show how flexion measurements can be used to reconstruct mass maps in 2-D projection on the sky, and in 3-D in combination with redshift data. Finally, we examine the predictions for cosmic flexion, including convergence-flexion cross-correlations, and find that the signal is an effective probe of structure on non-linear scales.Comment: 17 pages, including 12 figures, submitted to MNRA

CFHTLenS: the Canada–France–Hawaii Telescope Lensing Survey – imaging data and catalogue products

We present data products from the Canada–France–Hawaii Telescope Lensing Survey (CFHTLenS). CFHTLenS is based on the Wide component of the Canada–France–Hawaii Telescope Legacy Survey (CFHTLS). It encompasses 154 deg^2 of deep, optical, high-quality, sub-arcsecond imaging data in the five optical filters u*g′r′i′z′. The scientific aims of the CFHTLenS team are weak gravitational lensing studies supported by photometric redshift estimates for the galaxies. This paper presents our data processing of the complete CFHTLenS data set. We were able to obtain a data set with very good image quality and high-quality astrometric and photometric calibration. Our external astrometric accuracy is between 60 and 70 mas with respect to Sloan Digital Sky Survey (SDSS) data, and the internal alignment in all filters is around 30 mas. Our average photometric calibration shows a dispersion of the order of 0.01–0.03 mag for g′r′i′z′ and about 0.04 mag for u* with respect to SDSS sources down to i_(SDSS) ≤ 21. We demonstrate in accompanying papers that our data meet necessary requirements to fully exploit the survey for weak gravitational lensing analyses in connection with photometric redshift studies. In the spirit of the CFHTLS, all our data products are released to the astronomical community via the Canadian Astronomy Data Centre at http://www.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/CFHTLens/query.html. We give a description and how-to manuals of the public products which include image pixel data, source catalogues with photometric redshift estimates and all relevant quantities to perform weak lensing studies

Prelaunch testing of the GEOS-3 laser reflector array

The prelaunch testing performed on the Geos-3 laser reflector array before launch was used to determine the lidar cross section of the array and the distance of the center of gravity of the satellite from the center of gravity of reflected laser pulses as a function of incidence angle. Experimental data are compared to computed results

Vector coherent state representations, induced representations, and geometric quantization: II. Vector coherent state representations

It is shown here and in the preceeding paper (quant-ph/0201129) that vector coherent state theory, the theory of induced representations, and geometric quantization provide alternative but equivalent quantizations of an algebraic model. The relationships are useful because some constructions are simpler and more natural from one perspective than another. More importantly, each approach suggests ways of generalizing its counterparts. In this paper, we focus on the construction of quantum models for algebraic systems with intrinsic degrees of freedom. Semi-classical partial quantizations, for which only the intrinsic degrees of freedom are quantized, arise naturally out of this construction. The quantization of the SU(3) and rigid rotor models are considered as examples.Comment: 31 pages, part 2 of two papers, published versio

An exactly solvable model of a superconducting to rotational phase transition

We consider a many-fermion model which exhibits a transition from a superconducting to a rotational phase with variation of a parameter in its Hamiltonian. The model has analytical solutions in its two limits due to the presence of dynamical symmetries. However, the symmetries are basically incompatible with one another; no simple solution exists in intermediate situations. Exact (numerical) solutions are possible and enable one to study the behavior of competing but incompatible symmetries and the phase transitions that result in a semirealistic situation. The results are remarkably simple and shed light on the nature of phase transitions.Comment: 11 pages including 1 figur