8 research outputs found

    HPC Simulation of Magnetic Resonance Imaging

    Get PDF

    Особливості та умови формування властивостей техногенних ґрунтів

    Get PDF
    This work utilises general numerical magnetic resonance imaging MRI simulations to predict the spatial specificity of the blood oxygenation level-dependent (BOLD) functional MRI (fMRI) signal. A Monte Carlo simulation approach was utilized on a microvascular structure consisting of randomly oriented cylinders representing blood vessels. This framework was employed to numerically investigate the spatial specificity, defined as ratio of pial vessel to microvascular signal, of the spin echo BOLD fMRI signal as a function of field strength, echo time and tissue types [grey matter (GM) and cerebrospinal fluid (CSF), respectively]. Spatial specificity of spin echo BOLD fMRI signal was determined to increase with field strength up to 16 T and with maximal specificity for echo time shorter than tissue T(2). In addition, it was found that, for large pial vessels, the extravascular signal decay could not be described using the oversimplified but nevertheless commonly employed mono-exponential signal decay approximation (MEA). Consequently, a recently proposed model relying on the MEA deviates substantially from our results on the spatial specificity. A refinement of this model is proposed based on an available, more detailed signal description. Finally, the effect of CSF on the spatial specificity was investigated. While a large spatial specificity of the spin echo BOLD fMRI signal is observed if a pial vessel is surrounded by grey matter, this is greatly reduced for a pial vessel situated on a GM/CSF interface, rendering the suppression of pial vessels on the cortex surface unlikely

    Divergent evolution of terrestrial locomotor abilities in extant Crocodylia

    Get PDF
    Extant Crocodylia are exceptional because they employ almost the full range of quadrupedal footfall patterns (“gaits”) used by mammals; including asymmetrical gaits such as galloping and bounding. Perhaps this capacity evolved in stem Crocodylomorpha, during the Triassic when taxa were smaller, terrestrial, and long-legged. However, confusion about which Crocodylia use asymmetrical gaits and why persists, impeding reconstructions of locomotor evolution. Our experimental gait analysis of locomotor kinematics across 42 individuals from 15 species of Crocodylia obtained 184 data points for a wide velocity range (0.15–4.35 ms−1). Our results suggest either that asymmetrical gaits are ancestral for Crocodylia and lost in the alligator lineage, or that asymmetrical gaits evolved within Crocodylia at the base of the crocodile line. Regardless, we recorded usage of asymmetrical gaits in 7 species of Crocodyloidea (crocodiles); including novel documentation of these behaviours in 5 species (3 critically endangered). Larger Crocodylia use relatively less extreme gait kinematics consistent with steeply decreasing athletic ability with size. We found differences between asymmetrical and symmetrical gaits in Crocodylia: asymmetrical gaits involved greater size-normalized stride frequencies and smaller duty factors (relative ground contact times), consistent with increased mechanical demands. Remarkably, these gaits did not differ in maximal velocities obtained: whether in Alligatoroidea or Crocodyloidea, trotting or bounding achieved similar velocities, revealing that the alligator lineage is capable of hitherto unappreciated extreme locomotor performance despite a lack of asymmetrical gait usage. Hence asymmetrical gaits have benefits other than velocity capacity that explain their prevalence in Crocodyloidea and absence in Alligatoroidea—and their broader evolution

    High-performance computing MRI simulations

    No full text
    A new open-source software project is presented, JEMRIS, the Jülich Extensible MRI Simulator, which provides an MRI sequence development and simulation environment for the MRI community. The development was driven by the desire to achieve generality of simulated three-dimensional MRI experiments reflecting modern MRI systems hardware. The accompanying computational burden is overcome by means of parallel computing. Many aspects are covered that have not hitherto been simultaneously investigated in general MRI simulations such as parallel transmit and receive, important off-resonance effects, nonlinear gradients, and arbitrary spatiotemporal parameter variations at different levels. The latter can be used to simulate various types of motion, for instance. The JEMRIS user interface is very simple to use, but nevertheless it presents few limitations. MRI sequences with arbitrary waveforms and complex interdependent modules are modeled in a graphical user interface-based environment requiring no further programming. This manuscript describes the concepts, methods, and performance of the software. Examples of novel simulation results in active fields of MRI research are given

    High performance volume-of-intersection projectors for 3D-PET image reconstruction based on polar symmetries and SIMD vectorisation

    No full text
    For high-resolution, iterative 3D PET image reconstruction the efficient implementation of forward-backward projectors is essential to minimise the calculation time. Mathematically, the projectors are summarised as a system response matrix (SRM) whose elements define the contribution of image voxels to lines-of-response (LORs). In fact, the SRM easily comprises billions of non-zero matrix elements to evaluate the tremendous number of LORs as provided by state-of-the-art PET scanners. Hence, the performance of iterative algorithms, e.g. maximum-likelihood-expectation-maximisation (MLEM), suffers from severe computational problems due to the intensive memory access and huge number of floating point operations.Here, symmetries occupy a key role in terms of efficient implementation. They reduce the amount of independent SRM elements, thus allowing for a significant matrix compression according to the number of exploitable symmetries. With our previous work, the PET REconstruction Software TOolkit (PRESTO), very high compression factors (>300) are demonstrated by using specific non-Cartesian voxel patterns involving discrete polar symmetries. In this way, a pre-calculated memory-resident SRM using complex volume-of-intersection calculations can be achieved. However, our original ray-driven implementation suffers from addressing voxels, projection data and SRM elements in disfavoured memory access patterns. As a consequence, a rather limited numerical throughput is observed due to the massive waste of memory bandwidth and inefficient usage of cache respectively.In this work, an advantageous symmetry-driven evaluation of the forward-backward projectors is proposed to overcome these inefficiencies. The polar symmetries applied in PRESTO suggest a novel organisation of image data and LOR projection data in memory to enable an efficient single instruction multiple data vectorisation, i.e. simultaneous use of any SRM element for symmetric LORs. In addition, the calculation time is further reduced by using simultaneous multi-threading (SMT). A global speedup factor of 11 without SMT and above 100 with SMT has been achieved for the improved CPU-based implementation while obtaining equivalent numerical results

    Spatiotemporal properties of auditory intensity processing in multisensor MEG

    Full text link
    Loudness dependence of auditory evoked potentials (LDAEP) evaluates loudness processing in the human auditory system and is often altered in patients with psychiatric disorders. Previous research has suggested that this measure may be used as an indicator of the central serotonergic system through the highly serotonergic innervation of the auditory cortex. However, differences among the commonly used analysis approaches (such as source analysis and single electrode estimation) may lead to different results. Putatively due to discrepancies of the underlying structures being measured. Therefore, it is important to learn more about how and where in the brain loudness variation is processed. We conducted a detailed investigation of the LDAEP generators and their temporal dynamics by means of multichannel magnetoencephalography (MEG). Evoked responses to brief tones of five different intensities were recorded from 19 healthy participants. We used magnetic field tomography in order to appropriately localize superficial as well as deep source generators of which we conducted a time series analysis. The results showed that apart from the auditory cortex other cortical sources exhibited activation during the N1/P2 time window. Analysis of time courses in the regions of interest revealed a sequential cortical activation from primary sensory areas, particularly the auditory and somatosensory cortex to posterior cingulate cortex (PCC) and to premotor cortex (PMC). The additional activation within the PCC and PMC has implications on the analysis approaches used in LDAEP research
    corecore