A Wavelet-Based Approach for Estimating Time-Varying Connectivity in Resting-State Functional Magnetic Resonance Imaging.

Introduction: The selection of an appropriate window size, window function, and functional connectivity (FC) metric in the sliding window method is not straightforward due to the absence of ground truth. Methods: A previously proposed wavelet-based method was accordingly adjusted for estimating time-varying FC (TVFC) and was applied to a large high-quality, low-motion dataset of 400 resting-state functional magnetic resonance imaging data. Specifically, the wavelet coherence magnitude and relative phase were averaged across wavelet (frequency) scales to yield TVFC and synchronization patterns. To assess whether the observed fluctuations in TVFC were statistically significant (dynamic FC [dFC]; the distinction between TVFC and dFC is intentional), surrogate data were generated using the multivariate phase randomization (MVPR) and multivariate autoregressive randomization (MVAR) methods to define the null hypothesis of dFC absence. Results: By averaging across all frequencies, core regions of the default mode network (DMN; medial prefrontal and posterior cingulate cortices, inferior parietal lobes, hippocampal formation) were found to exhibit dFC (test-retest reproducibility of 90%) and were also synchronized in activity (-15° ≤ phase ≤15°). When averaging across distinct frequency bands, the same dynamic connections were identified, with the majority of them identified in the frequency range (0.01, 0.198) Hz, though with lower test-retest reproducibility (<66%). Additional analysis suggested that MVPR method better preserved properties (p < 10-10), including time-averaged coherence, of the original data compared with MVAR approach. Conclusions: The wavelet-based approach identified dynamic associations between the core DMN regions with fewer choices in parameters, compared with sliding window method. Impact statement We employed a wavelet-based method, previously used in the literature, and proposed modifications to assess time-varying functional connectivity in resting-state functional magnetic resonance imaging. With this approach, dynamic connections within the default mode network were identified, involving the medial prefrontal and posterior cingulate cortices, inferior parietal lobes, and hippocampal formation, which were also highly consistent in test-retest analysis (test-retest reproducibility of 90%), without the need to select window size, window function, and functional connectivity metric as with the sliding window method, whereby no consensus on the appropriate choices of hyperparameters currently exists in the literature.

Effects of motion related outliers in dynamic functional connectivity using the sliding window method.

BACKGROUND: It has been suggested that the use of window functions, other than the rectangular, in the sliding window method, may be beneficial for reducing the effects of motion-related outliers in the time-series, when assessing dynamic functional connectivity (dFC) in resting-state fMRI (rs-fMRI). METHODOLOGY: Ten window functions for a wide range of window lengths (20-150 s) combined with Pearson and Kendall correlation metrics, were investigated. One hundred high quality rs-fMRI datasets from healthy controls, were used to systematically assess the effect of varying the window function and length on dFC assessment. To this end, two approaches were implemented: a) simulated outliers were added to the experimental data and b) the experimental data were divided into low and high motion subgroups. RESULTS: The presence of experimental motion-noise tended to inflate the number of dynamic connections for longer (≥100 s) wide-shaped windows, while shorter (20-30 s) narrow-shaped windows exhibited increased sensitivity in the presence of simulated outliers. Moreover, window sizes from 60 s to 90 s were mildly affected by motion-related effects. In most cases, the number of dynamic connections increased, and gradually lower frequencies were captured, with an increasing window size. CONCLUSIONS: Subject motion considerably affects the obtained dFC patterns; thus, it is preferable to perform motion artefact removal in the pre-processing stage rather than using alternative window functions to mitigate their effects. Provided that motion-noise is not excessive, the choice of a rectangular window is adequate. Finally, low frequency oscillations in functional connectivity seem to play an important role in the context of dFC assessment.

Assessment of dynamic functional connectivity in resting-state fMRI using the sliding window technique.

INTRODUCTION: Recent studies related to assessing functional connectivity (FC) in resting-state functional magnetic resonance imaging have revealed that the resulting connectivity patterns exhibit considerable fluctuations (dynamic FC [dFC]). A widely applied method for quantifying dFC is the sliding window technique. According to this method, the data are divided into segments with the same length (window size) and a correlation metric is employed to assess the connectivity within these segments, whereby the window size is often empirically chosen. METHODS: In this study, we rigorously investigate the assessment of dFC using the sliding window approach. Specifically, we perform a detailed comparison between different correlation metrics, including Pearson, Spearman and Kendall correlation, Pearson and Spearman partial correlation, Mutual Information (MI), Variation of Information (VI), Kullback-Leibler divergence, Multiplication of Temporal Derivatives and Inverse Covariance. RESULTS: Using test-retest datasets, we show that MI and VI yielded the most consistent results by achieving high reliability with respect to dFC estimates for different window sizes. Subsequent hypothesis testing, based on multivariate phase randomization surrogate data generation, allowed the identification of dynamic connections between the posterior cingulate cortex and regions in the frontal lobe and inferior parietal lobes, which were overall in agreement with previous studies. CONCLUSIONS: In the case of MI and VI, a window size of at least 120 s was found to be necessary for detecting dFC for some of the previously identified dynamically connected regions.

An automatic panoramic image reconstruction scheme from dental computed tomography images.

OBJECTIVES: Panoramic images of the jaws are extensively used for dental examinations and/or surgical planning because they provide a general overview of the patient's maxillary and mandibular regions. Panoramic images are two-dimensional projections of three-dimensional (3D) objects. Therefore, it should be possible to reconstruct them from 3D radiographic representations of the jaws, produced by CBCT scanning, obviating the need for additional exposure to X-rays, should there be a need of panoramic views. The aim of this article is to present an automated method for reconstructing panoramic dental images from CBCT data. METHODS: The proposed methodology consists of a series of sequential processing stages for detecting a fitting dental arch which is used for projecting the 3D information of the CBCT data to the two-dimensional plane of the panoramic image. The detection is based on a template polynomial which is constructed from a training data set. RESULTS: A total of 42 CBCT data sets of real clinical pre-operative and post-operative representations from 21 patients were used. Eight data sets were used for training the system and the rest for testing. CONCLUSIONS: The proposed methodology was successfully applied to CBCT data sets, producing corresponding panoramic images, suitable for examining pre-operatively and post-operatively the patients' maxillary and mandibular regions.

Geometry-based vs. intensity-based medical image registration: A comparative study on 3D CT data.

Spatial alignment of Computed Tomography (CT) data sets is often required in numerous medical applications and it is usually achieved by applying conventional exhaustive registration techniques, which are mainly based on the intensity of the subject data sets. Those techniques consider the full range of data points composing the data, thus negatively affecting the required processing time. Alternatively, alignment can be performed using the correspondence of extracted data points from both sets. Moreover, various geometrical characteristics of those data points can be used, instead of their chromatic properties, for uniquely characterizing each point, by forming a specific geometrical descriptor. This paper presents a comparative study reviewing variations of geometry-based, descriptor-oriented registration techniques, as well as conventional, exhaustive, intensity-based methods for aligning three-dimensional (3D) CT data pairs. In this context, three general image registration frameworks were examined: a geometry-based methodology featuring three distinct geometrical descriptors, an intensity-based methodology using three different similarity metrics, as well as the commonly used Iterative Closest Point algorithm. All techniques were applied on a total of thirty 3D CT data pairs with both known and unknown initial spatial differences. After an extensive qualitative and quantitative assessment, it was concluded that the proposed geometry-based registration framework performed similarly to the examined exhaustive registration techniques. In addition, geometry-based methods dramatically improved processing time over conventional exhaustive registration.