Preface

Magnetic resonance imaging (MRI) is one of the most widely employed technologies for the clinical imaging of soft tissues due to its high contrast and spatial resolution, absence of ionizing radiation, and the current breadth of methodological applications. Much of the present diagnostic power of MRI is due to the technology's reliance on the nuclear magnetic moment of hydrogen and its manner of detection. Hydrogen is abundant in soft tissues and possesses a readily detectable magnetic moment in the presence of an externally applied magnetic field. Technical details for generating the field trace their origin to Bloch and Purcell's original method of demonstrating the nuclear magnetic phenomenon. In their technique, the interaction between the nuclear magnetic moment and the magnetic field generated a resonance frequency that was electronically detected by an inducting coil. In its clinical usage, a linear gradient was introduced in the magnetic field, enabling the resonance frequency to be varied as a function of its position along the gradient. This permitted the sampling of biological tissue via a Fourier space along a single axis one point at a time. Measurement of the signal with respect to phase differences along a second axis provided a unique set of values corresponding to each xy set of coordinates on the tissue plane.

This basic technological approach has since provided sufficient latitude to permit the wide-ranging evolution that has led to today's generation of MRI applications. A small sampling includes such well-established procedures as parallel imaging, which has substantially lessened the long acquisition intervals plaguing early clinical use; computational algorithms for preprocessing steps of co-registration and denoising, which are now routine; and functional MRI for brain activity, currently used widely in clinical research.

New advances are benefiting from a synergy of the current procedures together with the parallel evolution now occurring in the computational disciplines. Among the factors propelling the development of functional approaches, for example, are the currently evolved capabilities for rapid data acquisition, enhanced detection, and improved physiological models of function, which together underpin the growing trend toward monitoring dynamic properties in an increasingly data-driven, largescale computational environment. Structural analyses are also undergoing a significant evolution. While MRI has traditionally been used chiefly for structural tissue determinations, the difficulties in managing the burgeoning wealth of imagery data are increasingly requiring computational processing prior to examination by clinicians. Upcoming trends to address this need include new algorithms with enhanced computational power for feature extraction and pattern recognition. Computational developments also include specialized procedures that handle big data analyses, such as MRI fingerprinting and whole brain dynamics.

The current text presents select cases of how these trends are being incorporated into a new generation of MRI applications. The introductory chapter, Chapter 1, covers the current capabilities, providing a foundation for how the next generation of applications will benefit from knowledge previously gained in these standard areas.

The book's focus then shifts to how the current capabilities are incorporated in three domains that are experiencing especially rapid growth: functional and dynamic imaging in cardiac health and disease, advances in computational machine learning and its application to brain networks, and multiparametric, tissue-specific procedures.

determinations provide insight into dynamic brain topology and its relationship to brain states. Because of the large amounts of data generated by these studies, traditional datadriven methods for handling RS-fMRI data, such as independent components analysis and graph theoretic approaches, become unwieldy and lose descriptive power at elevated data levels. The chapter describes, among several other topics, current efforts to address big data handling that can be used to accurately gauge network structure and dynamics.

Among MRI's imaging capabilities is a capacity for distinguishing compositional differences and anomalies in different tissues. This is increasingly used to advantage in the determination of the many diseases affecting bone tissue. Normal bone always contains both fat and red marrow, and these undergo variations in fatty composition/ hematopoietic cellularity as part of a transformation phenomenon associated with these diseases. Marrow fat is distinguished by an elevated signal intensity due to its high proton content, whereas hematopoietic marrow has a much smaller signal. Using magnetic resonance spectroscopy, it is possible to quantify these differences, which enables a precise determination of changes occurring within the bone cavity.

Chapter 6 describes these advanced MRI procedures and their use in a wide range of major bone diseases and other clinical entities. New advances in magnetic field strength, gradient strength, and coils have enabled the development of novel pulse sequences, which have transformed knee MRI from a routine procedure of tissue structure to one of ultrastructure and even postoperative ligament reconstructions, which has previously constituted a challenge. Chapter 7 discusses these advanced procedures for finely delineating knee tissue structure. For example, isotropic threedimensional MRI permits thin-slice MRI with a spatial resolution of less than 1 mm and is seen in oblique planes, which help analyze complex structures like static and dynamic knee stabilizers. Enhanced three-dimensional visualization from various angles is expected to assist preoperative planning and achieve better clinical outcomes

in patients with difficult-to-treat ligament and meniscal injuries.

Diagnostic analyses have traditionally benefitted from studies exploring pulse sequences that optimize image features while minimizing artifact contributions and scanning time. The introduction of non-linear gradients, for example, was undertaken to improve image quality in highly undersampled MRI. Protocols for assessing pulse sequences are of interest for their utility in determining preferred sequences. Chapter 8 presents a pre-diagnostic protocol for assessing sequences that optimize these basic relationships. As a proof-of-principle demonstration, the procedure examines three frequently used pulse sequences: Fast Spin Echo (FSE), Fast Fluid Attenuated Inversion Recovery (FLAIR), and Conventional Spin Echo (CSE).

**Denis Larrivee** 

Pamplona, Spain

Mind and Brain Institute,

University of Navarra Medical School,

**Multiparametric, Tissue-Specific MRI Procedures**

### **Dynamic MRI in Cardiovascular Diagnosis**

Cardiac MRI (CMR) currently offers precise visualization of myocardial tissue structure and function that has proven to be critical for assessing cardiovascular disease. Chapter 2 describes the current broad range of CMR techniques and discusses how a similar range of next-generation applications are building on these abilities. The latter include, for example, updated structural mapping techniques that provide a comprehensive analysis and image of the entire heart with single scans for improved detection of fibrosis and edema that are often missed on older scanning procedures. Novel cardiac methods also include procedures for estimating cardiac tissue strain, obtained using the procedures of displacement encoding with stimulated echocardiography (DENSE) and feature tracking (FT), a post-processing algorithm that calculates myocardial deformation.

Improvements in the management of patients with congenital heart disease have significantly increased the number of patients surviving well into adulthood. Long-term care for such patients requires ongoing follow-up to detect and treat long-term complications. Chapter 3 describes an MRI procedure that measures four-dimensional heart flow in congenital patients. Four-dimensional flow allows a velocity assessment of the whole heart and major vessels with electrocardiogram gating and can be taken in a single acquisition. This is especially relevant for complex heart diagnoses. As one of multiple examples, it can assess the kinetic energy of blood flow, which provides a diagnostic measure for determining ventricular efficiency.

### **Machine Learning and Brain Network Analysis Methods for MRI**

The rapidly increasing volume of MRI data generated by new MRI procedures underscores the imperative for computational resources and serves to drive the development of computational capabilities in virtually all stages of MRI analysis, from acquisition and preprocessing to image reconstruction and feature extraction, and to big data issues. This has resulted in the current proliferation of machine learning methods that have enhanced image acquisition, analysis, and clinical decision-making. Chapter 4 provides a theoretical background for these methods and reviews how various applications are used in specific MRI stages. The chapter discusses, for example, convolutional neural networks, which have proven highly effective for image classification. Advanced computational analysis is of particular benefit for complex functional analyses, such as those employed in analyzing brain networks. Such computational procedures are increasingly applied to functional studies of brain activity.

Functional MRI (fMRI) of the quiescent brain, notably, has revealed that brain activity is highly organized into a group of brain networks (RSNs). Designated resting state, functional, magnetic resonance imaging (RS-fMRI), the method is currently the most powerful tool available for assessing the functional connectivity properties of brain networks. Chapter 5 discusses recent developments in RS-fMRI that seek to build on such functional connectivity determinations by relating causal sources of connectivity changes to brain states and behavior. Coupled with new modeling procedures, these

determinations provide insight into dynamic brain topology and its relationship to brain states. Because of the large amounts of data generated by these studies, traditional datadriven methods for handling RS-fMRI data, such as independent components analysis and graph theoretic approaches, become unwieldy and lose descriptive power at elevated data levels. The chapter describes, among several other topics, current efforts to address big data handling that can be used to accurately gauge network structure and dynamics.

### **Multiparametric, Tissue-Specific MRI Procedures**

Among MRI's imaging capabilities is a capacity for distinguishing compositional differences and anomalies in different tissues. This is increasingly used to advantage in the determination of the many diseases affecting bone tissue. Normal bone always contains both fat and red marrow, and these undergo variations in fatty composition/ hematopoietic cellularity as part of a transformation phenomenon associated with these diseases. Marrow fat is distinguished by an elevated signal intensity due to its high proton content, whereas hematopoietic marrow has a much smaller signal. Using magnetic resonance spectroscopy, it is possible to quantify these differences, which enables a precise determination of changes occurring within the bone cavity.

Chapter 6 describes these advanced MRI procedures and their use in a wide range of major bone diseases and other clinical entities. New advances in magnetic field strength, gradient strength, and coils have enabled the development of novel pulse sequences, which have transformed knee MRI from a routine procedure of tissue structure to one of ultrastructure and even postoperative ligament reconstructions, which has previously constituted a challenge. Chapter 7 discusses these advanced procedures for finely delineating knee tissue structure. For example, isotropic threedimensional MRI permits thin-slice MRI with a spatial resolution of less than 1 mm and is seen in oblique planes, which help analyze complex structures like static and dynamic knee stabilizers. Enhanced three-dimensional visualization from various angles is expected to assist preoperative planning and achieve better clinical outcomes in patients with difficult-to-treat ligament and meniscal injuries.

Diagnostic analyses have traditionally benefitted from studies exploring pulse sequences that optimize image features while minimizing artifact contributions and scanning time. The introduction of non-linear gradients, for example, was undertaken to improve image quality in highly undersampled MRI. Protocols for assessing pulse sequences are of interest for their utility in determining preferred sequences. Chapter 8 presents a pre-diagnostic protocol for assessing sequences that optimize these basic relationships. As a proof-of-principle demonstration, the procedure examines three frequently used pulse sequences: Fast Spin Echo (FSE), Fast Fluid Attenuated Inversion Recovery (FLAIR), and Conventional Spin Echo (CSE).

**Denis Larrivee** 

Mind and Brain Institute, University of Navarra Medical School, Pamplona, Spain

**1**

Section 1

Introduction

Section 1 Introduction

### **Chapter 1**

## Introductory Chapter: New Advances in MRI Clinical Analysis

*Denis Larrivee*

### **1. Introduction**

#### **1.1 Detecting the nuclear magnetic moment**

While magnetic resonance imaging has been used clinically since the early 1980s [1], the development of the technique traces its origin to Bloch and Purcell's demonstration of the nuclear magnetic phenomenon several decades earlier [2, 3]. Their work showed that in the presence of a constant magnetic field, nuclear magnetic moments generated a resonance frequency that could be electronically detected by means of an inducting coil. The ability to detect an electrical signal induced by the magnetic moment made possible the development of adjunct procedures for spatially profiling analytical specimens. Advances in MRI that have since propelled its clinical application have built on the physical principles underlying the magnetic moment and its detection used in the Bloch and Purcell demonstrations.

A key factor in the clinical use of the nuclear magnetic moment is the quantum mechanical effect present in hydrogen protons, which are abundant in soft tissue and confer the equivalent of an angular "spin" on the nucleus [4]. In the presence of an externally applied magnetic field, the magnetic moment generated by the spin processes about the axis of the applied field at a rate proportional to the external field and at a fixed angle to the field's axis, which defines the Larmor frequency. Within the static field, a second, perpendicular magnetic field applied by radiofrequency can excite the processing nuclear moment, changing the energy level and frequency of the magnetic moment with respect to the static field. In practice, this second field is applied in pulses typically lasting on the order of microseconds. Energy input from the pulses is then re-emitted between pulses as the magnetic moment returns to equilibrium, a process termed free induction decay. Relaxation to thermal equilibrium can occur either longitudinally, with loss of energy to the surrounding environment (a reduction in signal intensity, termed T1), or transversely, with energy exchange (but not loss) with neighboring nuclei (seen in a broadening of signal phase, termed T2). T1 and T2 values are defined as the lengths of time taken for their respective signals to decline by 63% (or to 37% from the equilibrium value), which adopt first order, exponential decay.

#### **1.2 Spatial resolution and localization**

Although the measured signal has its origin in the physical features inherent to the tissue, its detection and subsequent conversion to an anatomical image relies on technical features that amplify, resolve, and reconstruct the unique spatial and physical

properties of the image [5]. Localization of the tissue signal, notably, requires the use of magnetic field gradients, which is achieved by varying the static field strength during signal retrieval. In a typical application, the gradient introduces a linear variation in the static field strength along a tissue axis. Because spin precession frequency is directly proportional to field strength, there is a one-to-one correspondence between frequency and position. Gradients can assume any orthogonal direction and, in principle, can be used to reconstruct 3-dimensional images.

Faster or slower precession, and so higher and lower frequency, is detected as a higher or lower signal. Accordingly, manipulating the magnitude of the field strength range modulates the range of frequencies emitted for signal retrieval along a single axis. Imaging data acquired with small gradients has a small frequency range, limiting the frequency differential that can cover a pixel and so also limiting resolution. Larger gradients expand the frequency range, increasing resolution. When plotted, low frequencies are thus distributed centrally and high frequencies peripherally in a mathematical domain termed K-space [6, 7]; all image points thereby contain information from all frequencies. Planar localization requires encoding information along a second axis, a process also achieved by means of a static field gradient, but in which the phase difference along the second axis is monitored. Because phase differences can only be detected with respect to a reference value, a single scan is required for each phase encoded value. Information from both axes must then be mathematically transformed (by Fourier transform) to yield an x-y, planar image.

Because the typical and early use of the method—termed spin echo—sequentially varied the strength of the static field gradient, an important limitation in the use of MRI imaging has been the scanning time required to produce an adequately resolved image. Scanning time is a function of two factors: the interval between radiofrequency pulses, termed repetition time (TR), and the number of scans required for phase encoding and detection. The product of these yields the total scan time. Employing this protocol, typical scanning times could take on the order of minutes to hours, with higher resolution images needing longer intervals and so corresponding corrections for time related artifacts such as patient movement. Because of such artifacts, scanning time has posed a significant obstacle to image acquisition and quality.

#### **1.3 Evolution in MRI procedures**

Much of the early evolution in MRI methods thus sought to achieve a balance between data acquisition procedures that provided for sufficient spatial and temporal resolution and contrast and the time required to achieve a desired image. An example of this balancing is seen in the line reductions first introduced to the basic spin echo protocol [6]. A drawback to this approach, however, was the attendant loss in either resolution or image size due to the loss of frequency information. In the fast spin, echo planar or multi echo approach, for example, scanning time was reduced by taking multiple lines after the radio frequency pulse. A chief disadvantage of this latter approach, however, was the rapid loss in signal strength due to energy transferal during T2 acquisition, permitting only 3–4 lines per RF pulse, with significant deterioration in image quality.

While the objectives of reducing acquisition time and achieving enhanced spatial and temporal resolution have been of perennial interest in MRI development, new developments for addressing these and other objectives have in turn revealed new needs that invited corresponding solutions. The widespread use of parallel imaging [8], for example, was due to its ability to accelerate data acquisition while maintaining

#### *Introductory Chapter: New Advances in MRI Clinical Analysis DOI: http://dx.doi.org/10.5772/intechopen.113850*

high resolution through reductions in phase encoding lines and the use of algorithms that corrected subsampling aliasing. Such rapid acquisition methods led to the proliferation of multi-parametric approaches, which then facilitated quantification studies for MR fingerprinting [9]. These, and extended MR functional analyses, for instance, whole brain resting state fMRI [10], have engendered the further development of big data methods for the accompanying computational requirements, while data processing needs for processes such as feature extraction [11] propelled the use of machine learning to reduce diagnostic load. Indeed, the range of new techniques extends from advances in image acquisition to data processing and inference to computation and data handling, enhanced spatial resolution, and functional analyses. To illustrate the current proliferation, this introductory chapter to New Advances in Magnetic Resonance Imaging presents a select representation of techniques that trace their evolution from the physical and detection principles first identified by Bloch and Purcell. Chapters that follow each present a single, recently evolved MRI technique in greater detail. It is hoped that the novelty of these advances will be a source of inspiration for the engaged professional and interested scientist alike.

### **2. Acquisition procedures: advances in parallel imaging**

Parallel imaging is now used in nearly every clinical MRI scan for rapid data acquisition, for numerous reasons [8]. Many abdominal and cardiac scans, for example, are taken with patients holding their breath, with its obvious need for short scanning times. In other cases, for instance, multiline sequences following excitation pulses (e.g., in turbo spin echo), blurring artifacts are introduced due to the substantial T2 decay occurring during line retrieval. In still other applications, rapid data acquisition is essential because of the need to acquire large data sets [12, 13].

In parallel imaging, scanning time is reduced because the phased array coils yield unique views of the tissue objective, eliminating the scanning time for a significant portion of the region subject to gradient encoding [8, 14]. Due to the rapid decline in sensitivity of each coil element with increasing distance from the coil, this limits data acquisition to a clearly delineated, tissue profile. Individual images are combined to yield a comprehensive image. As a matter of principle, the maximum acceleration factor is related to the number of coils. Since most parallel imaging often employs 4–8 coil arrays and arrays containing 32, or even 128, channels are known (e.g., cardiac imaging), the reduction in scan time can be substantial.

Protocols for parallel imaging are currently classed according to whether aliased pixels are segregated in the imaging domain (SENSitivity Encoding or SENSE) or in K-space (Generalized Autocalibrating Partial Parallel Acquisition or GRAPPA) [14]. Techniques that act on the image domain first reconstruct and then correct, while those that act on the frequency domain first correct and then reconstruct. For those segregating in the imaging domain, prior knowledge of the coil sensitivity profiles enables separation of folded pixels from the undersampled image to recover the full image. Recent improvements have included phase-constrained SENSE. Whereas in conventional SENSE unknown variables are complex value, in phase-constrained SENSE, values are real, reducing the variable number by half. In contrast to the SENSE techniques, the GRAPPA algorithm is a K-space technique, which operates on acquired frequency information that is embedded in K-space. The technique is based on the principle that K-space information is shared between points in K-space due to the variation in multiple static field gradients; hence, missing information can be

computationally reconstructed from acquired K-space data. In order to reconstruct information from the missing points, acquired data points must be adjusted by weighting using an autocalibration signal obtained from another region of K-space.

Two recent modifications of the GRAPPA technique have included the 2D CAIPIRINHA and Wave-CAIPI techniques [15]. In 2D CAIPIRINHA, acceleration in the Ky and Kz directions is accompanied by a phase offset along Kz. This generates unique frequency patterns that require less computational resources for resolving aliasing. The Wave-CAIPI technique builds on 2D CAIPIRINHA by adding sinusoidal gradients along the Ky and Kz axes with a 180° shift. These additions amplify the acceleration by a factor of 9 over the 2D CAIPIRINHA technique. Unlike other very fast acquisitions, Wave-CAIPI is not subject to blurring from data gridding or artifacts from distortion due to uneven static magnetic fields.

### **3. Reconstruction procedures for image analysis**

The initial stages of image acquisition, preprocessing, and segmentation prepare the data for extraction of meaningful information. By definition, these steps involve removal of nonmeaningful or noise-based signals. One common source of noise, for example, is due to patient movement. Motion artifacts are corrected by registering sequential images, which can be carried out using available algorithms suited for medical imaging. A current standard for MRI image registration is the InsightToolKit (ITK), which contains a suite of algorithms including such processes as transformations, similarity metrics, and contrast normalization [11].

Recent trends in preprocessing and segmentation (e.g., denoising) have employed machine-based learning applications. Due to the labor involved, another major trend in machine-based learning procedures has been that of feature recognition and classification [16]. The high quantity of imaging data acquired from existing MRI scanning, especially, has made clinical diagnoses based on MRI images increasingly laborious, driving efforts for automated data abstraction and analysis. Computationally, machine learning relies on algorithms derived from neural network structures, which are composed of nodes joined by weighted edges. Inputs to nodes are weighted by a set of parameters and multiplied by transfer functions, for example, sigmoid and hyperbolic tangent functions, which transform weighted inputs. Among the most widely used of these deep neural networks (DNN) are the convolutional neural network (CNN), ResNet, the generative adversarial neural networks (GANs), and the U-nets.

Deep network models used for processing whole images are highly complicated, which significantly amplifies the processing time. For example, a CNN training model can often create millions of parameters in searches for feature classification. This drawback has been addressed by the use of image representations of smaller size that increases the processing efficiency. Most current feature extraction applications utilize wavelet transformation techniques in conjunction with neural network processing for MRI images [17]. The wavelet transform is used to remove extraneous detail and make the image more efficient for network processing; hence, approximated images have denser information content than original images. Convolutional layers then apply an initial filtering to produce an initial feature map, that is, what the network deems as unique features, which is refined by further network processing.

### **4. Multi-information sourcing**

Compared to methods for single parameter data analysis, models using multiparametric techniques offer the chief advantage of assessing correlations between multiple quantitative parameters of interest [9], with the potential for significantly greater accuracy. Various relevant MRI quantities include, but are not limited to, longitudinal (T1) and transverse (T2) relaxation times as well as retrospective synthesis of conventional MR contrasts, which are monitored in conjunction with novel techniques for rapid data acquisition and computational analysis.

Multi-parametric, analytical techniques currently employ similar strategies [18] with simultaneous sampling of the parameter and K-spaces, in which transient-state data are obtained by varying the acquisition parameters and undersampled K-space snapshots are taken after each excitation. Parametric maps are then computed using a physical model *via* the Bloch equations. Several procedures that have evolved from this approach include magnetic resonance fingerprinting (MRF) [18] and quantitative transient-state imaging (QTI).

In MRF, modulating the MRI sequence parameters across the time domain yields a time series of weighted MRI images, with each tissue having a unique MRI signal fingerprint. Such fingerprints can be computationally simulated and a dictionary of tissue specific fingerprints built from the simulations. During image reconstruction, putative fingerprints are matched to this dictionary. The fingerprint with the greatest correlation in the dictionary provides the MRI parameters for a given voxel. After all voxels are analyzed, parametric maps are then constructed. Due to its ability to recognize very specific structural elements, MRF has the potential for diagnosis of a wide variety of clinical conditions [19–22]. Novel techniques like quantitative sequencing enable clinical approaches capable of accelerated quantitative mapping of dynamic physiological processes. For example, such techniques have been used to examine blood flow, with computation made of velocity scalars or vectors that could be employed in cardiac assessments. In one study, computations of scalar velocities were computed in a direction perpendicular to a vessel slice based on multi-parametric T1, T2, and proton density recordings [23]. A drawback to such methods is that of their reliance on physical models of the physiological events the quantitative mapping is intended to simulate, with the potential for loss of valuable data. Moreover, the complexity of such models can demand increasing computational resources, significantly extending the data acquisition period.

### **5. Functional approaches using MRI**

Blood oxygenation-dependent contrast, termed BOLD or fMRI imaging [24], was developed to indirectly assess neural activity in the brain by monitoring activityinduced changes in blood oxygenation. BOLD exploits the neurovascular mechanism of hyperemia, whereby localized brain activity recruits increased blood oxygenation within the region of activation. In typically used protocols, the image sequence relies on T2\* weighting with scan times under 5 seconds to measure the hemodynamic response function. Since its discovery, use of this procedure has undergone extraordinary growth [13, 24].

The dependence of the BOLD signal on neurovascular mechanisms, however, has meant that fMRI is also constrained by limitations inherent in the hemodynamic response. Chief among these is the much slower response time than the underlying neural processes being measured. Temporal information of spiking events is therefore heavily blurred, requiring the use of mathematical processing, like that of the general linear model or experimental block protocols, to infer event related, signal activity [24]. With processing, temporal resolution in the 100 ms range can be achieved, which is roughly tenfold slower than the neural events being monitored.

#### **5.1 High-strength imaging**

Another challenge to fMRI is that of the low signal to noise ratio (SNR), a consequence of data acquisition, and constraints limiting the extent of preprocessing that can be performed. One of the technical improvements now being attempted to improve the SNR is the use of high-strength magnetic fields, which enhances the anatomical specificity of imaging. Whereas most scanning is done using 3 T fields, equipment using 7 T fields is becoming increasingly prevalent. At these higher field strengths, it has been shown that less spatial smoothing is required and neural activity in cases of resting state networks displays higher correlation coefficients, indicating greater spatial resolution [25–27]. On the other hand, use of higher field strengths has several drawbacks, including longer sampling intervals, inhomogeneous magnetic field properties, and the logarithmic growth in specific absorption rate (SAR) with increasing field strength [28, 29].

#### **5.2 Multimodal studies (w EEG/MEG)**

Among the techniques used to overcome the temporal limitations of fMRI have been multimodal approaches, which combine fMRI with such methods as the EEG or MEG. Both EEG and MEG display rapid temporal responses, with the capability for resolving neural events at millisecond scales. Use of these approaches in conjunction with fMRI is thus premised on the greatly improved temporal resolution offered by these procedures. Advanced technologies have now been developed to simultaneously record EEG and fMRI signal, which help to understand the relationship between the spatial and temporal characteristics of physiological signals [30]. Nonetheless, in comparison to fMRI alone, combined approaches have seen limited use. In the case of EEG, spatial resolution is greatly inferior to that of fMRI and MEG approaches that suffer from source localization issues. This means that the experimental design or clinical assessment must clearly attribute the signal sources prior to drawing experimental and/or clinical conclusions in such combined approaches.

#### **5.3 DIANA fMRI**

The persistence of interpretive difficulties with multimodal approaches has generated a long-standing interest in the development of alternative methods capable of both high spatial and temporal resolution. One recently developed method has merged the detection of ultra-weak magnetic fields generated by neural electrical activity with the fMRI detection of the hemodynamic response [31]. This approach, termed Direct Imaging of Neuronal Activity for functional MRI (DIANA-fMRI), interleaves K-space lines used for imaging the hemodynamic response with a K-space line that directly measures the ultra-weak magnetic field. Millisecond resolution is achieved using fast, low-angle shot (FLASH); gradient-echo imaging; and short repetition intervals (e.g., 5 ms). When carried out at high field strengths (9.4 Tesla),

signal to noise ratios are reported to be in the range of 20 to 1. To date, the technique has been used only on animal models.

While promising, the technique suffers from its own unique set of interpretive uncertainties. For example, electromagnetic effects based on neuronal current models appear to be ruled out as these oppose the direction of the observed DIANA response. Moreover, in humans, increasing the stimulus duration does not lead to correlative signal changes, suggesting that the signal may be confounded by other interactions such as inflow effects and subject motion. These uncertainties suggest that further development of the approach will require improved understanding of the biophysical factors contributing to the DIANA signal.

#### **5.4 Resting state fMRI**

One of the most significant fMRI developments is the use of fMRI during rest, coined resting state-fMRI (RS-fMRI). RS-fMRI focuses on spontaneous lowfrequency fluctuations (< 0.1 Hz) in the BOLD signal that occur in the absence of task-related activities. The functional significance of these fluctuations was first recognized by Biswal et al. [32] in a study in which subjects were told not to perform any cognitive, language, or motor tasks. After determining the correlation between the BOLD time course of a seed region identified by bilateral finger tapping and that of all other areas in the brain, the authors found that fluctuations in the left somatosensory cortex were highly correlated with homologous areas in the contralateral hemisphere. This observed correlation led to their conclusion that such "resting networks" manifested the functional connectivity of the brain. The observation of spontaneous, synchronous fluctuations occurring between brain regions has since led to studies that have identified as many as 7 to 17 other stable networks [33], with 7 consistently agreed upon.

Because characterization of resting state networks (RSNs) in the human brain relies on the analysis of temporal fluctuations in the blood oxygenation leveldependent (BOLD) signal, the delineation of RSNs has been directly dependent on the ability of fMRI to detect neural activity [6]. The dependence on the BOLD signal means that RS-fMRI shares advantages that accrue to fMRI—the ability to monitor neural activity, albeit indirectly—as well as disadvantages that characterize its use. Chief among these limitations is fMRI's temporal resolution, which is dependent on the hemodynamic response time [34]. Accordingly, a key factor in the use of RS-fMRI is the measurement of neural activity fluctuations rather than spiking events per se.

Early studies of RSN functional connectivities, like that of Biswal et al. [32, 35], relied on the selection of regions of interest based on investigator preferences. However, while the simplicity and interpretability of the ROI technique make it procedurally facile and a frequently adopted approach, the method relies entirely on user-defined ROIs and so is limited for network discovery by its a priori, selected criteria. Due to this caveat and coupled with the evolution of mathematical models and improved computational capabilities, there has been a paradigm shift from that of imposing initial conditions on the data to that of extracting patterns of brain activity directly from the raw time series. The main example of this approach is independent component analysis [36]. In this approach, the time series signal is assumed to be due to multiple spatiotemporal processes that are statistically independent of each other. By extracting the independent signals, various time courses of specific brain regions can be constructed and grouped into maps representative of their spatial distribution. Another approach to the interpretation of RS-fMRI datasets employs graph theory,

where activity sources comprise nodes and connectivity defines the edges that link these nodes [10, 36]. Unlike ICA, which focuses chiefly on the strength of correlation between different domains, graph theory characterizes the features of network topology. The graph theory approach describes the interaction between nodes by means of such graph parameters as average path length, clustering coefficients, node degree, centrality measures, and level of modularity. Graph theory is thus a promising technique for exploring the integration and segregation of networks in the brain. Among the topological features studied, modularity, the assessment of the presence of functionally independent units or modules that compose resting state networks, has increasingly been used to characterize functional adjustments occurring during behavior, network perturbations, or pathologies that affect network function, revealing significant alterations in such pathologies as stroke [37] and psychiatric disease [38, 39].

In principle, inferences of causality from directed functional connectivity determinations can be extended to brain-wide neuronal dynamics. Empirical studies from RS-fMRI, for example, show that RSNs are differentiated on the basis of their metastability and synchrony [40]. These and similar observations have led to models of brain function and behavior that predict that the human brain at rest operates at maximum metastability, that is, in a state of maximal network switching. The demonstration of RSN properties like metastability thus suggests that directed connectivity changes may be used to assess the construction of brain states. The methodological question that arises is that of generating a descriptive approach relating functional neuroimaging data to whole brain dynamics. Recent attempts to address this question have adopted two approaches. The first employs a BOLD, data-driven, computational method that leverages the method of *recurrence structure analysis* (RSA), a mathematical procedure derived from Poincare's recurrence theorem [41]. This "recurrent" behavior can be described by a *recurrence plot method* (RP), which allows a matrixbased visualization of recurrent states. The second approach posits the governance of RSN dynamics by a ground state global attractor. This global ground state is mathematically described as a stable stationary solution representing a point of maximal stability in a landscape of stationary points (nodes) that information flows toward or away from [42]. This theoretical framework has been shown to successfully account for the highly structured dynamics arising from spontaneous brain activity in RSNs [43].

### **6. Conclusion and future directions**

Building on Bloch and Purcell's method for detecting the nuclear magnetic resonance phenomenon, clinical applications have evolved at virtually all stages along the MRI procedural pipeline. Long intervals for data acquisition that plagued early clinical use have been substantially lessened; algorithms for preprocessing with coregistration and denoising have become routine; machine learning for tissue-specific fingerprinting and feature extraction has been introduced, and functional methods for cardiac and brain diagnosis are reliably and regularly used. Future directions will see the roles of adjunct computational resources continue to grow. MRI has traditionally been a procedure for structural tissue determinations, and the assessing of structural differences between normal and trauma or disease states has typically been performed by clinicians. Aiding them will be new algorithms and enhanced computational power for improved feature extraction and pattern recognition capabilities.

Clinical diagnoses will increasingly rely on dynamic and physiological parameters associated with functional performance.

In the nervous system, the drive to consolidate improved temporal resolution with the current abilities for high spatial resolution, for instance, will seek to enlist methods successfully used for high spatial resolution in novel ways that match neuronal activity timescales, such as the interleaving of K-space data coupled with theoretical models that enable direct neural activity monitoring. Noninvasive techniques will also introduce perturbation studies to monitor how these influence brain dynamics. Together, these and other innovations will amplify the already leading clinical role of MRI-based applications.

### **Author details**

Denis Larrivee1,2

1 Mind and Brain Institute, University of Navarra Medical School, Spain

2 Loyola University Chicago, USA

\*Address all correspondence to: sallar1@aol.com

© 2023 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## **References**

[1] Smith FW, Hutchison JM, Mallard JR, et al. Oesophageal carcinoma demonstrated by whole-body nuclear magnetic resonance imaging. British Medical Journal (Clinical Research Edition). 1981;**282**(6263):510-512

[2] Bloch F, Hansen WW, Packard ME. Nuclear induction. Physics Review. 1946;**69**:127

[3] Purcell EM, Torrey HC, Pound RV. Resonance absorption by nuclear magnetic moments in a solid. Physics Review. 1946;**69**:37-38

[4] Grover VPB, Tognarelli JM, Crossey MME, et al. Magnetic resonance imaging: Principles and techniques: Lessons for clinicians. Journal of Clinical and Experimental Hepatology. 2015;**5**(3):246-255

[5] Moratal D, Vallés-Luch A, Martí-Bonmatí L, Brummer M. k-space tutorial: An MRI educational tool for a better understanding of k-space. Biomedical Imaging Intervention Journal. 2008;**4**(1):e15. DOI: 10.2349/biij.4.1.e15

[6] Mezrich R. A perspective on K space. Radiology. 1995;**195**(2):297-315

[7] Sýkora S. K space formulation of MRI. In: Sykora S, editor. Stan's Library. Vol. Volume I. Castano Primo, Italy: Permalink; 2005

[8] Deshmane A, Eng M, Gulani V, et al. Parallel MR imaging. Journal of Magnetic Resonance Imaging. 2012;**36**(1):55-72. DOI: 10.1002/jmri.23639

[9] Gomez A, Molina-Romera M, Buonincontri G, et al. Designing contrasts for rapid, simultaneous parameter quantification and flow visualization with quantitative transient-state imaging. Journal of Clinical and Experimental Hepatology. 2015;**5**(3):246-255

[10] Yang J, Gohel S, Vachha B. Current methods and new directions in resting state fMRI. Clinical Imaging. 2020;**65**:47- 53. DOI: 10.1016/j.clinimag.2020.04.004

[11] Vadmal V, Junno G, Badye C, et al. MRI image analysis methods and applications. Neuro-Oncology Advances. 2020;**2**(1):1-13

[12] Gordon Y, Partovi S, Müller-Eschner M, et al. Dynamic contrast-enhanced magnetic resonance imaging: Fundamentals and application to the evaluation of the peripheral perfusion. Cardiovascular Diagnostic Therapy. 2014;**4**(2):147-164. DOI: 10.3978/j. issn.2223-3652.2014.03.01

[13] Bandettini P. The spatial, temporal, and interpretive limits of functional MRI. In: Davis K, Charney D, Coyle JT, Nemeroff C, editors. Neuropsychopharmacology: The Fifth Generation of Progress. Philadelphia: Lippincott, Williams, and Wilkins; 2002

[14] Hamilton J, Franson D, Seiberlich N. Recent advances in parallel imaging for MRI. Progress in Nuclear Magnetic Resonance Spectroscopy. 2017;**101**:71-95. DOI: 10.1016/j.pnmrs.2017.04.002

[15] Calogero C. Recent advances in parallel imaging for MRI: WAVE-CAIPI technique. Journal of Advanced Health Care. 2022;**4**(1):23-26

[16] Khalila M, Ayada H, Adiba A. Performance evaluation of feature extraction techniques in MR-brain image *Introductory Chapter: New Advances in MRI Clinical Analysis DOI: http://dx.doi.org/10.5772/intechopen.113850*

classification system. Procedia Computer Science. 2018;**127**:218-225

[17] Fayaz M, Torokeldiev N, Turdumamatov S, et al. An efficient methodology for brain MRI classification based on DWT and convolutional neural network. Sensors. 2021;**21**:7480. DOI: 10.3390/ s21227480

[18] Ma D, Gulani V, Seiberlich N, et al. Magnetic resonance fingerprinting. Nature. 2013;**495**:187-192

[19] Zhao B, Haldar JP, Setsompop K, Wald LL. Optimal experiment design for magnetic resonance fingerprinting. In: 2016 38th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC). Lake Buena Vista, FL. pp. 453-456. DOI: 10.1109

[20] Liao C, Bilzic B, Manhard MK, et al. 3D MR fingerprinting with accelerated stack-of-spirals and hybrid slidingwindow and GRAPPA reconstruction. NeuroImage. 2017;**162**:13-22. DOI: 10.1016/j.neuroimage.2017.08.030

[21] Marty B, Carlier PG. Physiological and pathological skeletal muscle T1 changes quantified using a fast inversion-recovery radial NMR imaging sequence. Science Reports. 2019;**9**:6852. DOI: 10.1038/

s41598-019-43398-x

[22] Buonincontri G, Sawiak S. Three-dimensional MR fingerprinting with simultaneous B1 estimation. Magnetic Resonance Medicine. 2015;**2015**:1-9

[23] Cao X et al. A model-based velocity mapping of blood flow using MR fingerprinting. Proceedings of the International Society for Magnetic Resonance in Medicine. 2017;**25**:0941

[24] Loued-Khenissi L, Doll O, Preuschoff K. An overview of functional magnetic resonance imaging techniques for organizational research. Organizational Research Methods. 2019;**22**(1):17-45

[25] Yacoub E, Van De Moortele PF, Shmuel A, et al. Signal and noise characteristics of Hahn SE and GE BOLD fMRI at 7 T in humans. NeuroImage. 2005;**2005**(24):738-750. DOI: 10.1016/j. neuroimage

[26] Pohmann R, Speck O, Scheffler K. 2016. Signal-to-noise ratio and MR tissue parameters in human brain imaging at 3, 7, and 9.4 Tesla using current receive coil arrays. Magnetic Resonance Medicine. 2016;**75**:801-809. DOI: 10.1002/mrm.25677

[27] Branco P, Seixas D, Castro SL. Temporal reliability of ultra-high field resting- state MRI for single-subject sensorimotor and language mapping. NeuroImage. 2018;**168**:499-508. DOI: 10.1016/j.neuroimage.2016.11.029

[28] Raimondo L, Icaro AF, Jurjen HO, et al. Advances in resting state fMRI acquisitions for functional connectomics. NeuroImage. 2021;**243**:118503

[29] Van de Moortele PF, Auerbach EJ, Olman C, et al. T1 weighted brain images at 7 Tesla unbiased for proton density, T2<sup>∗</sup> contrast and RF coil receive B1 sensitivity with simultaneous vessel visualization. NeuroImage. 2009;**46**:432-446. DOI: 10.1016/j. neuroimage.2009.02.009

[30] Fleury M, Figueiredo P, Vourvopoulos A, and Lécuyer A. Two is better? Combining EEG and fMRI for BCI and Neurofeedback: A systematic review. 2023. hal-04038069

[31] Toi PT, Jang JH, Kyeong SM, et al. In vivo direct imaging of neuronal activity at high temporospatial resolution. Science. 2022;**378**(6616): 160-168. DOI: 10.1126/science.abh43

[32] Biswal B, Yetkin FZ, Haughton VM, Hyde JS. Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magnetic Resonance Medicine. 1995;**34**(4): 537-541. DOI: 10.1002/mrm.1910340409

[33] Damoiseaux SA, Rombouts RB, Barkhof F, Beckman CF. Consistent resting-state networks across healthy subjects. Proceedings of the National Academy of Sciences. 2006;**103**(37):13848-13853. DOI: 10.1073

[34] Ogawa S, Tank DW, Menon R, et al. Intrinsic signal changes accompanying sensory stimulation: Functional brain mapping with magnetic resonance imaging. Proceedings of the National Academy of Science USA. 1992;**89**:5951-5955

[35] Yeo BTT, Krienen FM, Sepulcre J, et al. The organization of the human cerebral cortex estimated by intrinsic functional connectivity. Journal of Neurophysiology. 2011;**106**:1125-1165. DOI: 10.1152/jn.00338.2011

[36] Smitha KA, Akhil RK, Arun KM, et al. Resting state fMRI: A review on methods in resting state connectivity analysis and resting state networks. The Neuroradiology Journal. 2017;**30**(4):305-317. DOI: 10.1177/1971400917697342

[37] Corbetta M, Siegel JS, Schulman GL. On the low dimensionality of behavioral deficits and alterations of brain network connectivity after focal injury. Cortex. 2018;**107**:229-237

[38] Lerman-Sinkoff DB, Barch DM. Network community structure alterations in adult schizophrenia: Identification and localization of alterations. Neuroimage Clinic. 2016;**10**:96-106. DOI: 10.1016/j. nicl.2015.11.011

[39] Bullmore ET, Sporns O. Complex brain networks: Graph theoretical analysis of structural and functional systems. Nature Review Neuroscience. 2009;**10**:186-198. DOI: 10.1038/nrn2575

[40] Lee WH, Frangou S. Linking functional connectivity and dynamic properties of resting-state networks. Scientific Reports. 2017;**7**:16610. DOI: 10.1038/s41598-017-16789-1

[41] Beim Graben P, Jimenez-Marin A, Diez I, Cortes JM, et al. Metastable resting state brain dynamics. Frontiers in Computational Neuroscience. 2019;**13**:62. DOI: 10.3389/fncom.2019.00062

[42] Carvalho A, Langa J, Robinson J. Attractors for Infinite-Dimensional Nonautonomous Dynamical Systems. New York: Springer; 2012

[43] Gonzalez L, Panda R, Ponce-Alvarez A, et al. Loss of consciousness reduces the stability of brain hubs and the heterogeneity of brain dynamics. Communications Biology. 2021;**4**(1):34489535

## Section 2
