Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Medical Imaging”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Gamma Ray Detection In Liquid Argon Time Projection Chambers

Gamma rays are photons in excess of 100 keV of energy and can be produced by a variety of mechanisms like radioactive decay and inelastic neutrino-argon interactions. A reconstruc- tion of Compton scattering gamma rays was developed for use in liquid argon time projection chambers. This reconstruction was tested with detectors developed to demonstrate the ca- pabilities of Deep Underground Neutrino Experiment (DUNE) Near Detector Liquid Argon target (ND LAr). In October 2024, a radioactive fluorine-18 source was placed next to the Full Scale Demonstrator. The data obtained was used to test this reconstruction algorithm. It was able to determine the direction of the photons with an error of less than 10 percent. In addition to this work, there has been significant work on designing a liquid argon based Positron Emission Tomography (PET) scanner. PET scanners detect gamma rays emitted by a tracer chemical to help screen for diseases like cancer. In addition to the applications to medical imaging, this thesis briefly explores the applications of this reconstruction to de-excitation photons produced in interactions between neutrinos and argon nuclei. These applications could improve studies focusing on neutrinos originating from accelerators and supernovae.

Murphy, Thomas [Syracuse U.]↗
DOE OSTI2023

An extended focused assessment with sonography in trauma ultrasound tissue-mimicking phantom for developing automated diagnostic technologies

Medical imaging-based triage is critical for ensuring medical treatment is timely and prioritized. However, without proper image collection and interpretation, triage decisions can be hard to make. While automation approaches can enhance these triage applications, tissue phantoms must be developed to train and mature these novel technologies. Here, we have developed a tissue phantom modeling the ultrasound views imaged during the enhanced focused assessment with sonography in trauma exam (eFAST). The tissue phantom utilized synthetic clear ballistic gel with carveouts in the abdomen and rib cage corresponding to the various eFAST scan points. Various approaches were taken to simulate proper physiology without injuries present or to mimic pneumothorax, hemothorax, or abdominal hemorrhage at multiple locations in the torso. Multiple ultrasound imaging systems were used to acquire ultrasound scans with or without injury present and were used to train deep learning image classification predictive models. Performance of the artificial intelligent (AI) models trained in this study achieved over 97% accuracy for each eFAST scan site. We used a previously trained AI model for pneumothorax which achieved 74% accuracy in blind predictions for images collected with the novel eFAST tissue phantom. Grad-CAM heat map overlays for the predictions identified that the AI models were tracking the area of interest for each scan point in the tissue phantom. Overall, the eFAST tissue phantom ultrasound scans resembled human images and were successful in training AI models. Tissue phantoms are critical first steps in troubleshooting and developing medical imaging automation technologies for this application that can accelerate the widespread use of ultrasound imaging for emergency triage.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Neural network for 3D inertial confinement fusion shell reconstruction from single radiographs

In inertial confinement fusion (ICF), x-ray radiography is a critical diagnostic for measuring implosion dynamics, which contain rich three-dimensional (3D) information. Traditional methods for reconstructing 3D volumes from 2D radiographs, such as filtered backprojection, require radiographs from at least two different angles or lines of sight (LOS). In ICF experiments, the space for diagnostics is limited, and cameras that can operate on fast timescales are expensive to implement, limiting the number of projections that can be acquired. To improve the imaging quality as a result of this limitation, convolutional neural networks (CNNs) have recently been shown to be capable of producing 3D models from visible light images or medical x-ray images rendered by volumetric computed tomography. Here, we propose a CNN to reconstruct 3D ICF spherical shells from single radiographs. We also examine the sensitivity of the 3D reconstruction to different illumination models using preprocessing techniques such as pseudo-flatfielding. To resolve the issue of the lack of 3D supervision, we show that training the CNN utilizing synthetic radiographs produced by known simulation methods allows for reconstruction of experimental data as long as the experimental data are similar to the synthetic data. We also show that the CNN allows for 3D reconstruction of shells that possess low mode asymmetries. Further comparisons of the 3D reconstructions with direct multiple LOS measurements are justified.

3D Reconstruction↗
DOE OSTI2024

Terahertz time-domain spectroscopy imaging of pancreatic ductal adenocarcinoma tissues

Pancreatic ductal adenocarcinoma (PDAC) ranks among the malignancies with the highest fatality and morbidity rates. This is predominantly attributable to an absence of understanding the intricate and diverse microenvironment of the tumor. We use terahertz time-domain spectroscopy (THz-TDS) imaging in transmission geometry to probe ex-vivo the heterogenous microenvironment of the genetically modified murine PDAC tissue that closely resembles the PDAC heterogeneity in human malignancy. We introduced a maximum a-posteriori probability estimation algorithm to objectively the tumor’s heterogenous microenvironment using the average values of refractive index and absorption coefficient within the useable terahertz bandwidth as imaging markers. Furthermore, direct comparison of stained histopathologic images and the refractive index and the absorption coefficient high-resolution, two dimensional maps of the same PDAC samples confirms the high potential of the THz-TDS method for tumor tissue characterization.

absorption↗
DOE OSTI2021

Computed Tomography Scanning and Geophysical Measurements of the Wabash No.1 Core

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia were used to characterize Mt. Simon core from the Wabash CarbonSAFE project, Wabash No. 1 well. The primary impetus of the work in this report is a collaboration between NETL Research and Innovation Center and the Illinois State Geological Survey at the Prairie Research Institute, University of Illinois Urbana-Champaign in Champaign, Illinois. The resultant datasets are presented in this report and can be accessed from NETL's Energy Data eXchange (EDX) online system using the following link: https://edx.netl.doe.gov/dataset/wabash-no-1-well. The Wabash No. 1 well was drilled as part of a U.S. Department of Energy (DOE) funded Carbon Storage Assurance Facility Enterprise field project to assess the feasibility of developing a commercial-scale geological storage complex at the Wabash Valley Resources Integrated Gasification Combined Cycle plant near Terre Haute, Indiana. All equipment and techniques used were non-destructive, enabling future examinations and analyses to be performed on these cores. None of the equipment used was suitable for direct visualization of the pore space in fine-grained structures; fractures, discontinuities, and millimeter scale features were readily detectable with the methods tested. Imaging with the NETL medical CT scanner was performed on the entire core. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), P-wave, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for more detailed analysis. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core; the resulting macro and micro descriptions are relevant to many subsurface energy-related examinations traditionally performed at NETL.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Core Characterization of Bakken Shale from the Bedwell 33-52-1-1H Well, Sheridan County, Montana

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia were used to characterize core of the Bedwell 33-52-1-1H well (API 25-091-21920), a wildcat field in Sheridan County, Montana. The primary impetus of this work is to provide a unique dataset to researchers at universities, national laboratories, geological surveys, and other locations for additional analyses. The resultant datasets are presented in this report and can be accessed from NETL's Energy Data eXchange (EDX) online system using the following link: https://edx.netl.doe.gov/dataset/bedwell-33-52-1-1h-well. The equipment and techniques used to characterize the full core were non-destructive, enabling future examinations and analyses to be performed on these cores. However, none of the equipment used was suitable for direct visualization of the pore space in fine-grained structures such as those found in shales; fractures, discontinuities, and millimeter scale features were readily detectable with the methods tested. Imaging with the NETL medical CT scanner was performed on the entire core. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), P-wave, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for more detailed analysis. A selection of samples (Table 1) were cored, sliced, and powdered for additional characterization analyses. These analyses did involve destructive subsampling of the full core at discrete sections, and the results are presented here in the context of the larger core description. Selected sections of the core were removed from the full core for detailed analyses using the following destructive techniques: bulk inductively coupled plasma-optical emission spectroscopy (ICP-OES), isolation of carbon and sulfur, Fourier transform infrared (FTIR) spectroscopy, scanning electron microscopy (SEM), and higher resolution CT scan with NETL’s micro-CT scanning systems. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core; the resulting macro and micro descriptions are relevant to many subsurface energy related examinations traditionally performed at NETL.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

Computed Tomography Scanning and Geophysical Measurements of UW Enterprises LP 1-250512-129 Well in Southwestern Indiana

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia, were used to characterize Illinois Basin core of the Upper Devonian-Early Mississippian New Albany Shale Formation from Posey County, Indiana. The primary impetus of this work is a collaboration between Indiana Geological and Water Survey at Indiana University Bloomington, NETL, and Woolsey Operating Company LLC to characterize and make publicly available core information from the New Albany Shale of the Illinois Basin. Core characterization of this unconventional oil/gas well will aid in understanding the lithology changes and the fracture complexity of the New Albany Shale. There is a potential for this formation to be developed in the future for critical mineral and rare earth element (CM/REE) extraction. The resultant datasets are presented in this report and can be accessed from NETL's Energy Data eXchange (EDX) online system using the following link: https://edx.netl.doe.gov/dataset/uw-enterprises. All equipment and techniques used were non-destructive, enabling future examinations and analyses to be performed on this core. None of the equipment used was suitable for direct visualization of the pore space in the fine-grained structures studied; however, fractures, discontinuities, and millimeter-scale features were readily detectable with the methods tested. Imaging with the NETL medical CT scanner was performed on the entire core. Targeted higher resolution CT scanning of select sections was performed with NETL’s industrial and micro-CT scanner. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), P-wave, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for more detailed analysis. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core; the resulting macro and micro descriptions are relevant to many subsurface energy related examinations traditionally performed at NETL.

58 GEOSCIENCES↗
DOE OSTI2023

Computed Tomography Scanning and Geophysical Measurements of Appalachian Basin Core from the Jones and Laughlin #1 Well, Beaver County, PA

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia, were used to characterize Appalachian Basin core from Beaver County, Pennsylvania. The primary impetus of this work is a collaboration between the U.S. Department of Energy (DOE) and the Pennsylvania Geological Survey to characterize and make publicly available core information from the Onondaga-Huntersville formations of the Appalachian Basin. This stratigraphic well and the core data produced in this report will aid in understanding the structural complexities of the Onondaga-Huntersville formations. The resultant datasets are presented in this report and can be accessed from NETL's Energy Data eXchange (EDX) online system using the following link: https://edx.netl.doe.gov/dataset/jonesandlaughlin1well. All equipment and techniques used were non-destructive, enabling future examinations and analyses to be performed on these cores. Fractures, discontinuities, and millimeter-scale features were readily detectable with imaging performed with the NETL medical CT scanner over the entire core. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for further study. Targeted higher resolution CT scanning of select sections was performed with NETL’s micro-CT scanner. The combination of methods used provides a multiscale analysis of the core; the resulting macro and micro descriptions are relevant to many subsurface energy related examinations traditionally performed at NETL.

58 GEOSCIENCES↗
DOE OSTI2023

Health Hazards of Exposures to Radioiodine

Iodine is a chemical element with atomic number 53. Iodine-127 is stable (non-radioactive) and commonly found in nature. Elemental iodine is a purple-colored solid at room temperature and pressure, but spontaneously sublimates (turns into vapor). Iodine is an essential element for life, and is required for proper functioning of the thyroid. Iodine is present in many foods, and is readily absorbed by the body and concentrated in the thyroid gland. A fraction of iodine ingested or inhaled is rapidly removed by the kidneys. The rest of the inhaled or ingested iodine is absorbed the by thyroid and retained for many months. Iodine has a biological half-life of approximately 120 days in health individuals. The biological half-life can be shorter in individuals with hyperthyroidism, and longer in individuals with hypothyroidism. Iodine has a number of radioactive isotopes, most of which have relatively short half-lives (days or weeks). Short half-life iodine isotopes are useful for a variety of medical applications, including imaging and cancer therapy. For example, Iodine-123 (half-life 13 hours) is commonly used for medical imaging of the thyroid, while iodine-131 (half-life 8 days) is used for suppressing thyroid function in individuals with hyperthyroidism or ablating (killing) thyroid cells to treat thyroid cancer. Iodine-125 (half-life 59 days) is produced in nuclear reactors, and has medical uses. Although iodine-125 can be used for thyroid imaging, Iodine-123 is more commonly used for that purpose because of its shorter half-life and higher-energy emissions. Iodine-125 is more commonly used for cancer treatment, and can be processed into small metal pellets (seeds) inserted directly into a tumor. Iodine-125 emits low-energy x-rays which can kill tumor cells and generally cannot escape the tumor, sparing other tissues. Medical iodine for imaging or treatment is typically administered orally in the form a pill or liquid solution. A typical adult thyroid scan using iodine-123 involves having the patient swallow between one and four 0.1 millicuries pills, with the exact dose dependent on the patient’s weight. This results in a whole-body committed effective dose of 80 – 320 mrem, and a thyroid equivalent dose of 1443 – 5772 mrem. Note that the whole-body effective dose relates to the overall cancer risk, while the larger equivalent dose to the thyroid only indicates that most of this risk is the result of exposure to the thyroid. These doses are considered safe, although the procedure is not recommended for pregnant or breastfeeding women. In contrast, the quantities of iodine-131 used for treatment of hyperthyroidism and thyroid cancers are much higher. For treatment of hyperthyroidism, 4 – 10 millicuries are administered, while for thyroid cancer the administration can range from 50 – 150 millicuries of I-131. In addition to medical exposures, large populations were exposed to radioiodine as result of the atomic bombings of Hiroshima and Nagasaki in Japan, and the Chernobyl nuclear accident. These populations have been carefully followed for many years to assess the effect of their radiation exposures on cancer risk. As a result, a great deal is known about the cancer risks associated with radioiodine exposure. Because iodine is concentrated in the thyroid, the principal risk of exposure to radioiodine is thyroid cancer. Children have the highest risk of thyroid cancer after exposure to radioactive iodine. According to a large study of Japanese atomic bomb survivors, an effective dose of one Sievert (100,000 mrem) has been observed to increase the risk of thyroid cancer by a factor of 9.5 in children aged zero to nine years old, by a factor of 3 in children aged 10 to 19 years old, and by barely detectable amounts in adults. Another way of quantifying the risk from radioiodine exposure is from risk coefficients, which provide the risk per unit intake of radionuclides in terms of both morbidity (any cancer) and mortality (death). Both morbidity (risk of cancer) and mortality (death) risks are shown in the table below. Note that because thyroid cancer is almost never fatal, the morbidity coefficients are much larger than the mortality coefficients.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2025

Modeling inter‐reader variability in clinical target volume delineation for soft tissue sarcomas using diffusion model

Abstract Background Accurate delineation of the clinical target volume (CTV) is essential in the radiotherapy treatment of soft tissue sarcomas. However, this process is subject to inter‐reader variability due to the need for clinical assessment of risk and extent of potential microscopic spread. This can lead to inconsistencies in treatment planning, potentially impacting treatment outcomes. Most existing automatic CTV delineation methods do not account for this variability and can only generate a single CTV for each case. Purpose This study aims to develop a deep learning‐based technique to generate multiple CTV contours for each case, simulating the inter‐reader variability in the clinical practice. Methods We employed a publicly available dataset consisting of fluorodeoxyglucose positron emission tomography (FDG‐PET), x‐ray computed tomography (CT), and pre‐contrast T1‐weighted magnetic resonance imaging (MRI) scans from 51 patients with soft tissue sarcoma, along with an independent validation set containing five additional patients. An experienced reader drew a contour of the gross tumor volume (GTV) for each patient based on multi‐modality images. Subsequently, two additional readers, together with the first one, were responsible for contouring three CTVs in total based on the GTV. We developed a diffusion model‐based deep learning method that is capable of generating arbitrary number of different and plausible CTVs to mimic the inter‐reader variability in CTV delineation. The proposed model incorporates a separate encoder to extract features from the GTV masks, leveraging the critical role of GTV information in accurate CTV delineation. Results The proposed diffusion model demonstrated superior performance with the highest Dice Index (0.902 compared to values below 0.881 for state‐of‐the‐art models) and the best generalized energy distance (GED) (0.209 compared to values exceeding 0.221 for state‐of‐the‐art models). It also achieved the second‐highest recall and precision metrics among the compared ambiguous image segmentation models. Results from both datasets exhibited consistent trends, reinforcing the reliability of our findings. Additionally, ablation studies exploring different model structures and input configurations highlighted the significance of incorporating prior GTV information for accurate CTV delineation. Conclusions The proposed diffusion model successfully generates multiple plausible CTV contours for soft tissue sarcomas, effectively capturing inter‐reader variability in CTV delineation.

Dong, Yafei [Yale Biomedical Imaging Institute Yal↗
DOE OSTI2022

NaroNet: Discovery of tumor microenvironment elements from highly multiplexed images

Understanding the spatial interactions between the elements of the tumor microenvironment -i.e. tumor cells. fibroblasts, immune cells- and how these interactions relate to the diagnosis or prognosis of a tumor is one of the goals of computational pathology. We present NaroNet, a deep learning framework that models the multi-scale tumor microenvironment from multiplex-stained cancer tissue images and provides patient-level interpretable predictions using a seamless end-to-end learning pipeline. Trained only with multiplex-stained tissue images and their corresponding patient-level clinical labels, NaroNet unsupervisedly learns which cell phenotypes, cell neighborhoods, and neighborhood interactions have the highest influence to predict the correct label. To this end, NaroNet incorporates several novel and state-of-the-art deep learning techniques, such as patch-level contrastive learning, multi-level graph embeddings, a novel max-sum pooling operation, or a metric that quantifies the relevance that each microenvironment element has in the individual predictions. We validate NaroNet using synthetic data simulating multiplex-immunostained images where a patient label is artificially associated to the -adjustable- probabilistic incidence of different microenvironment elements. We then apply our model to two sets of images of human cancer tissues: 336 seven-color multiplex-immunostained images from 12 high-grade endometrial cancer patients; and 382 35-plex mass cytometry images from 215 breast cancer patients. In both synthetic and real datasets, NaroNet provides outstanding predictions of relevant clinical information while associating those predictions to the presence of specific microenvironment elements.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Computed Tomography Scanning and Petrophysical Measurements of Eastern Williston Basin Twin Buttes and Hagel Formations

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the U.S. Department of Energy’s (DOE) National Energy Technology Laboratory (NETL) in Morgantown, West Virginia, were used to characterize core from two wells that represent coal resources across North Dakota. These include the MC23080C Well in Mercer County and the 23-B001 Well in Oliver County. The primary impetus of this work was to capture a detailed digital representation of the core from the MC23080C and 23-B001 Wells. The collaboration between the NETL and the Energy and Environment Research Center (EERC) enables other research entities to access information about this potential carbon ore, rare earth, and critical mineral resource plays in the Williston Basin. All equipment and techniques used were non-destructive, enabling future examinations and analyses to be performed on these cores. Fractures, discontinuities, and millimeter-scale features were readily detectable with the medical CT scanner acquired images. Imaging with the NETL medical CT scanner was performed on entire cores. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), gamma density, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for potential future analysis. Higher-resolution industrial and micro-CT images were acquired from selected zones along the depth of the core to visualize the structure in higher detail. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core, with the resulting macro- and micro-descriptions relevant to many subsurface energy-related examinations traditionally performed at NETL.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

High-Density Polyethylene Custom Focusing Lenses for High-Resolution Transient Terahertz Biomedical Imaging Sensors

Transient terahertz time-domain spectroscopy (THz-TDS) imaging has emerged as a novel non-ionizing and noninvasive biomedical imaging modality, designed for the detection and characterization of a variety of tissue malignancies due to their high signal-to-noise ratio and submillimeter resolution. We report our design of a pair of aspheric focusing lenses using a commercially available lens-design software that resulted in about 200 × 200-μm2 focal spot size corresponding to the 1-THz frequency. The lenses are made of high-density polyethylene (HDPE) obtained using a lathe fabrication and are integrated into a THz-TDS system that includes low-temperature GaAs photoconductive antennae as both a THz emitter and detector. The system is used to generate high-resolution, two-dimensional (2D) images of formalin-fixed, paraffin-embedded murine pancreas tissue blocks. The performance of these focusing lenses is compared to the older system based on a pair of short-focal-length, hemispherical polytetrafluoroethylene (TeflonTM) lenses and is characterized using THz-domain measurements, resulting in 2D maps of the tissue refractive index and absorption coefficient as imaging markers. For a quantitative evaluation of the lens effect on the image resolution, we formulated a lateral resolution parameter, R2080, defined as the distance required for a 20–80% transition of the imaging marker from the bare paraffin region to the tissue region in the same image frame. The R2080 parameter clearly demonstrates the advantage of the HDPE lenses over TeflonTM lenses. The lens-design approach presented here can be successfully implemented in other THz-TDS setups with known THz emitter and detector specifications.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Computed Tomography Scanning and Petrophysical Measurements of the Lively Grove #1 Well Core

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia were used to characterize core from the Lively Grove #1 Well (API 121892494700), drilled near Marissa, Washington County, Illinois. Core from the well was obtained as part of the Illinois Storage Corridor’s Carbon Storage Assurance Facility Enterprise (CarbonSAFE) project (DE FE0031892). The primary impetus of this work was to capture a detailed digital representation of the core from the Lively Grove #1 Well. The collaboration between the U.S. Department of Energy’s (DOE) NETL and the Illinois State Geological Survey (ISGS) enables other research entities to access information about this potential carbon storage location and Cambro-Ordovician Storage Complex formations of the Illinois Basin. The resultant datasets are presented in this report and can be accessed from NETL’s Energy Data eXchange (EDX) online system. All equipment and techniques used were non-destructive, enabling future examinations and analyses to be performed on these cores. Fractures, discontinuities, and millimeter-scale features were readily detectable with the medical CT scanner-acquired images. Imaging with the NETL medical CT scanner was performed on entire cores. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), gamma density, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for potential future analysis. Higher-resolution industrial and micro-CT images were acquired from selected zones along the depth of the core to visualize the structure in higher detail. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core; with the resulting macro- and micro-descriptions relevant to many subsurface energy related examinations traditionally performed at NETL.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

High-dimensional multivariate autoregressive model estimation of human electrophysiological data using fMRI priors

Multivariate autoregressive (MVAR) model estimation enables assessment of causal interactions in brain networks. However, accurately estimating MVAR models for high-dimensional electrophysiological recordings is challenging due to the extensive data requirements. Hence, the applicability of MVAR models for study of brain behavior over hundreds of recording sites has been very limited. Prior work has focused on different strategies for selecting a subset of important MVAR coefficients in the model to reduce the data requirements of conventional least-squares estimation algorithms. Here we propose incorporating prior information, such as resting state functional connectivity derived from functional magnetic resonance imaging, into MVAR model estimation using a weighted group least absolute shrinkage and selection operator (LASSO) regularization strategy. The proposed approach is shown to reduce data requirements by a factor of two relative to the recently proposed group LASSO method of Endemann et al (Neuroimage 254:119057, 2022) while resulting in models that are both more parsimonious and more accurate. The effectiveness of the method is demonstrated using simulation studies of physiologically realistic MVAR models derived from intracranial electroencephalography (iEEG) data. The robustness of the approach to deviations between the conditions under which the prior information and iEEG data is obtained is illustrated using models from data collected in different sleep stages. This approach allows accurate effective connectivity analyses over short time scales, facilitating investigations of causal interactions in the brain underlying perception and cognition during rapid transitions in behavioral state.

62 RADIOLOGY AND NUCLEAR MEDICINE↗