Minimal Artifacts Research Articles

Particle motion artifacts in equilibrium magnetization measurements of large iron oxide nanoparticles

Iron oxide nanoparticles find many applications due to their response when subjected to externally applied magnetic fields. Equilibrium magnetization measurements (commonly known as magnetization curves) are an essential characterization tool to evaluate if particles display hysteresis, and to obtain magnetic properties such as the saturation magnetization. For superparamagnetic particles, one can obtain a magnetic size distribution by fitting the data to a theoretical model, such as the Langevin function, in what is called magnetogranulometric analysis. If one wishes to use the resulting size estimates as a predictor of particle performance in applications, magnetization data must be obtained under conditions that capture the response of the particles with minimal artifacts. In this paper, we used selected iron oxide nanoparticle batches with physical size ranging from 20 to 45 nm to demonstrate the influence of sample preparation methods on the magnetization data obtained. We show that measurements in powder form and in liquid solvents display varying degrees of particle interaction artifacts at low fields, depending strongly on particle size and on the thickness of the surface coating. In addition, measurements in ‘solid’ waxy hydrocarbon matrices are shown to be susceptible to particle rotation artifacts for large particle sizes. Hard crosslinked polymer matrices are shown to restrict particle motion completely, resulting in magnetization data that follows the Langevin function if the measurement is performed above the blocking temperature of the particles. We end with a discussion of how the presence of matrix-dependent measurement artifacts influence the magnetic diameter fits obtained using magnetogranulometry, and how measuring above and below the blocking temperature can affect fit results.

Direct pixel to pixel principal strain mapping from tagging MRI using end to end deep convolutional neural network (DeepStrain)

Regional soft tissue mechanical strain offers crucial insights into tissue's mechanical function and vital indicators for different related disorders. Tagging magnetic resonance imaging (tMRI) has been the standard method for assessing the mechanical characteristics of organs such as the heart, the liver, and the brain. However, constructing accurate artifact-free pixelwise strain maps at the native resolution of the tagged images has for decades been a challenging unsolved task. In this work, we developed an end-to-end deep-learning framework for pixel-to-pixel mapping of the two-dimensional Eulerian principal strains varvec{{varepsilon }}_{boldsymbol{p1}} and varvec{{varepsilon }}_{boldsymbol{p2}} directly from 1-1 spatial modulation of magnetization (SPAMM) tMRI at native image resolution using convolutional neural network (CNN). Four different deep learning conditional generative adversarial network (cGAN) approaches were examined. Validations were performed using Monte Carlo computational model simulations, and in-vivo datasets, and compared to the harmonic phase (HARP) method, a conventional and validated method for tMRI analysis, with six different filter settings. Principal strain maps of Monte Carlo tMRI simulations with various anatomical, functional, and imaging parameters demonstrate artifact-free solid agreements with the corresponding ground-truth maps. Correlations with the ground-truth strain maps were R = 0.90 and 0.92 for the best-proposed cGAN approach compared to R = 0.12 and 0.73 for the best HARP method for varvec{{varepsilon }}_{boldsymbol{p1}} and varvec{{varepsilon }}_{boldsymbol{p2}}, respectively. The proposed cGAN approach's error was substantially lower than the error in the best HARP method at all strain ranges. In-vivo results are presented for both healthy subjects and patients with cardiac conditions (Pulmonary Hypertension). Strain maps, obtained directly from their corresponding tagged MR images, depict for the first time anatomical, functional, and temporal details at pixelwise native high resolution with unprecedented clarity. This work demonstrates the feasibility of using the deep learning cGAN for direct myocardial and liver Eulerian strain mapping from tMRI at native image resolution with minimal artifacts.

Radiopaque drug-eluting embolisation beads as fiducial markers for stereotactic liver radiotherapy.

Objective:To determine the feasibility of using radiopaque (RO) beads as direct tumour surrogates for image-guided radiotherapy (IGRT) in patients with liver tumours after transarterial chemoembolisation (TACE).Methods:A novel vandetanib-eluting RO bead was delivered via TACE as part of a first-in-human clinical trial in patients with either hepatocellular carcinoma or liver metastases from colorectal cancer. Following TACE, patients underwent simulated radiotherapy imaging with four-dimensional computed tomography (4D-CT) and cone-beam CT (CBCT) imaging. RO beads were contoured using automated thresholding, and feasibility of matching between the simulated radiotherapy planning dataset (AVE-IP image from 4D data) and CBCT scans assessed. Additional kV, MV, helical CT and CBCT images of RO beads were obtained using an in-house phantom. Stability of RO bead position was assessed by comparing 4D-CT imaging to CT scans taken 6–20 days following TACE.Results:Eight patients were treated and 4D-CT and CBCT images acquired. RO beads were visible on 4D-CT and CBCT images in all cases and matching successfully performed. Differences in centre of mass of RO beads between CBCT and simulated radiotherapy planning scans (AVE-IP dataset) were 2.0 mm mediolaterally, 1.7 mm anteroposteriorally and 3.5 mm craniocaudally. RO beads in the phantom were visible on all imaging modalities assessed. RO bead position remained stable up to 29 days post TACE.Conclusion:RO beads are visible on IGRT imaging modalities, showing minimal artefact. They can be used for on-set matching with CBCT and remain stable over time.Advances in knowledge:The role of RO beads as fiducial markers for stereotactic liver radiotherapy is feasible and warrants further exploration as a combination therapy approach.

Oesophageal IGRT considerations for SBRT of LA-NSCLC: barium-enhanced CBCT and interfraction motion

BackgroundTo determine the optimal volume of barium for oesophageal localisation on cone-beam CT (CBCT) for locally-advanced non-small cell lung cancers (NSCLC) and quantify the interfraction oesophageal movement relative to tumour.MethodsTwenty NSCLC patients with mediastinal and/or hilar disease receiving radical radiotherapy were recruited. The first five patients received 25 ml of barium prior to their planning CT and alternate CBCTs during treatment. Subsequent five patient cohorts, received 15 ml, 10 ml and 5 ml. Six observers contoured the oesophagus on each of the 107 datasets and consensus contours were created. Overall 642 observer contours were generated and interobserver contouring reproducibility was assessed. The kappa statistic, dice coefficient and Hausdorff Distance (HD) were used to compare barium-enhanced CBCTs and non-enhanced CBCTs. Oesophageal displacement was assessed using the HD between consensus contours of barium-enhanced CBCTs and planning CTs.ResultsInterobserver contouring reproducibility was significantly improved in barium-enhanced CBCTs compared to non-contrast CBCTs with minimal difference between barium dose levels. Only 10 mL produced a significantly higher kappa (0.814, p = 0.008) and dice (0.895, p = 0.001). The poorer the reproducibility without barium, the greater the improvement barium provided. The median interfraction HD between consensus contours was 4 mm, with 95% of the oesophageal displacement within 15 mm.Conclusions10 mL of barium significantly improves oesophageal localisation on CBCT with minimal image artifact. The oesophagus moves substantially and unpredictably over a course of treatment, requiring close daily monitoring in the context of hypofractionation.

Design and Validation of a Multi-Point Injection Technology for MR-Guided Convection Enhanced Delivery in the Brain.

Convection enhanced delivery (CED) allows direct intracranial administration of neuro-therapeutics. Success of CED relies on specific targeting and broad volume distributions (VD). However, to prevent off-target delivery and tissue damage, CED is typically conducted with small cannulas and at low flow rates, which critically limit the maximum achievable VD. Furthermore, in applications such as gene therapy requiring injections of large fluid volumes into broad subcortical regions, low flow rates translate into long infusion times and multiple surgical trajectories. The cannula design is a major limiting factor in achieving broad VD, while minimizing infusion time and backflow. Here we present and validate a novel multi-point cannula specifically designed to optimize distribution and delivery time in MR-guided intracranial CED of gene-based therapeutics. First, we evaluated the compatibility of our cannula with MRI and common viral vectors for gene therapy. Then, we conducted CED tests in agarose brain phantoms and benchmarked the results against single-needle delivery. 3T MRI in brain phantoms revealed minimal susceptibility-induced artifacts, comparable to the device dimensions. Benchtop CED of adeno-associated virus demonstrated no viral loss or inactivation. CED in agarose brain phantoms at 3, 6, and 9 μL/min showed >3x increase in volume distribution and 60% time reduction compared to single-needle delivery. This study confirms the validity of a multi-point delivery approach for improving infusate distribution at clinically-compatible timescales and supports the feasibility of our novel cannula design for advancing safety and efficacy of MR-guided CED to the central nervous system.

Dual-phase Au-Pt alloys free from magnetic susceptibility artifacts in magnetic resonance imaging

PurposeMagnetic resonance imaging (MRI) devices are frequently used in image-based diagnosis. In the case of large artifacts, which are generated in magnetic resonance (MR) images when magnetic materials, such as metals, are present in the body, these devices are less useful. This study aimed to develop a dual-phase Au-Pt alloy that does not generate artifacts in MR images and has high workability to prepare medical devices. Materials and methodsA processing method to produce a dual-phase Au-Pt alloy was established, and the magnetic susceptibility and artifacts of different alloy compositions were determined using a SQUID (superconducting quantum interference device) flux meter and a 1.5 T-MRI system. The crystallographic phases of the prepared alloy samples were identified using X-ray diffraction. Sample cross-sections were observed using a metallurgical microscope. Furthermore, a thinning test was conducted to examine alloy workability. ResultsDual-phase Au-Pt alloys Au70Pt30 and Au67Pt33—the former heat-treated at 800 and 850 °C and the latter heat-treated at 900 °C—generated minimal artifacts when imaged in a 1.5 T-MRI system. Their volume magnetic susceptibility increased as the heat-treatment temperature decreased. The alloy surfaces were observed to be uniform. Moreover, the workability of the dual-phase alloy was considerably better than that of the single-phase alloy. ConclusionVolume magnetic susceptibility could be controlled by changing the composition and processing temperature of the Au-Pt alloys. Dual-phase Au-Pt alloys those do not generate magnetic susceptibility artifacts in MRI images and have good workability could be prepared. The alloys are expected to be used in the preparation of various implantable medical devices.

Hepatic Iron Quantification Using a Free-Breathing 3D Radial Gradient Echo Technique and Validation With a 2D Biopsy-Calibrated R2* Relaxometry Method.

Hepatic iron content (HIC) is an important parameter for the management of iron overload. Non-invasive HIC assessment is often performed using biopsy-calibrated two-dimensional breath-hold Cartesian gradient echo (2D BH GRE) R2* -MRI. However, breath-holding is not possible in most pediatric patients or those with respiratory problems, and three-dimensional free-breathing radial GRE (3D FB rGRE) has emerged as a viable alternative. To evaluate the performance of a 3D FB rGRE and validate its R2* and fat fraction (FF) quantification with 3D breath-hold Cartesian GRE (3D BH cGRE) and biopsy-calibrated 2D BH GRE across a wide range of HICs. Retrospective. Twenty-nine patients with hepatic iron overload (22 females, median age: 15 [5-25] years). Three-dimensional radial and 2D and 3D Cartesian multi-echo GRE at 1.5 T. R2* and FF maps were computed for 3D GREs using a multi-spectral fat model and 2D GRE R2* maps were calculated using a mono-exponential model. Mean R2* and FF values were calculated via whole-liver contouring and T2* -thresholding by three operators. Inter- and intra-observer reproducibility was assessed using Bland-Altman and intraclass correlation coefficient (ICC). Linear regression and Bland-Altman analysis were performed to compare R2* and FF values among the three acquisitions. One-way repeated-measures ANOVA and Wilcoxon signed-rank tests, respectively, were used to test for significant differences between R2* and FF values obtained with different acquisitions. Statistical significance was assumed at P < 0.05. The mean biases and ICC for inter- and intra-observer reproducibility were close to 0% and >0.99, respectively for both R2* and FF. The 3D FB rGRE R2* and FF values were not significantly different (P > 0.44) and highly correlated (R2 ≥ 0.98) with breath-hold Cartesian GREs, with mean biases ≤ ±2.5% and slopes 0.90-1.12. In non-breath-holding patients, Cartesian GREs showed motion artifacts, whereas 3D FB rGRE exhibited only minimal streaking artifacts. Free-breathing 3D radial GRE is a viable alternative in non-breath-hold patients for accurate HIC estimation. 3 TECHNICAL EFFICACY: Stage 2.

High fidelity fiber orientation density functions from fiber ball imaging.

The fiber orientation density function (fODF) in white matter is a primary physical quantity that can be estimated with diffusion MRI. It has often been employed for fiber tracking and microstructural modeling. Requirements for the construction of high fidelity fODFs, in the sense of having good angular resolution, adequate data to avoid sampling errors, and minimal noise artifacts, are described for fODFs calculated with fiber ball imaging. A criterion is formulated for the number of diffusion encoding directions needed to achieve a given angular resolution. The advantages of using large b-values (≥6000 s/mm2 ) are also discussed. For the direct comparison of different fODFs, a method is developed for defining a local frame of reference tied to each voxel's individual axonal structure. The Matusita anisotropy axonal is proposed as a scalar fODF measure for quantifying angular variability. Experimental results, obtained at 3 T from human volunteers, are used as illustrations.

Feasibility of Abdominal Magnetic Resonance Imaging in Patients With Residual Concentrated Enteric Contrast After Fluoroscopic Abdominal Examination.

This study aimed to evaluate the image quality, image artifacts, radiologist confidence, and ability to provide definitive diagnosis for all patients with magnetic resonance imaging (MRI) performed after an abdominal fluoroscopic examination and to determine the utility of MRI in this setting. Thirty-one MRI examinations performed a median of 2 days after fluoroscopic bowel evaluation (barium, n = 13; iodine, n = 18), 20 within 3 days of MRI, were retrospectively reviewed. The image quality, artifact emanating from bowel, inhomogeneity artifact, radiologist confidence, ability to render a definitive diagnosis, and identification of emergent or important findings for all MRI examinations were assessed. These same features were evaluated on 5 computed tomographies performed after fluoroscopy (before the MRI) in the same cohort. All 31 MRI examinations performed after fluoroscopic studies with concentrated barium or iodine solutions were diagnostic for answering the clinical question according to radiologist and report review, regardless of magnet strength and type of fluoroscopic contrast ingested. Magnetic resonance imaging after fluoroscopy had excellent overall image quality (mean score, 4.74/5), minimal to no artifact emanating from bowel (mean, 4.63/5), minimal inhomogeneity artifact (mean, 4.38/5), and excellent diagnostic confidence (mean, 4.98/5). No additional imaging was necessary for diagnosis after MRI. Computed tomography after fluoroscopy had lower overall image quality, more image artifacts, and lower diagnostic confidence (P < 0.05). Magnetic resonance imaging is a useful tool for evaluating patients with retained concentrated enteric contrast from recent fluoroscopic examinations. In the absence of contraindication, MRI should be considered in the evaluation of urgent clinical problems in patients who recently underwent a fluoroscopic bowel evaluation.

Minimal artifact actively shimmed metallic needles in MRI.

Signal voids caused by metallic needles pose visualization and monitoring challenges in many MRI applications. In this work, we explore a solution to this problem in the form of an active shim insert that fits inside a needle and corrects the field disturbance (ΔB0 ) caused by the needle outside of it. The ΔB0 induced by a 4 mm outside-diameter titanium needle at 3T is modeled and a two-coil orthogonal shim set is designed and fabricated to shim the ΔB0 . Signal recovery around the needle is assessed in multiple orientations in a water phantom with four different pulse sequences. Phase stability around the needle is assessed in an ex-vivo porcine tissue dynamic gradient echo experiment with and without shimming. Additionally, heating of the shim insert is assessed under 8 min of continuous operation with 1A current and concurrent imaging. An average recovery of ~63% of lost signal around the needle across orientations is shown with active shimming with a maximum current of 1.172 A. Signal recovery and correction of the underlying ΔB0 is shown to be independent of imaging sequence. Needle-induced phase gradients outside the perceptible signal void are also minimized with active shimming. Temperature rise of up to 0.9° Celsius is noted over 8 min of continuous 1A active shimming operation. A sequence independent method for minimization of metallic needle induced signal loss using an active shim insert is presented. The method has potential benefits in a range of qualitative and quantitative interventional MRI applications.

A Hybrid Synthetic Aperture and Time-Reversal MUSIC Algorithm for Subwavelength Radar Imaging

A promising chipless RFID approach uses millimeter-wave synthetic aperture radar (SAR) to image metal ink-printed ID tags from a meter or more away. Due to printing cost, it is desirable to minimize the size and spacing of metal patches within a tag, preferably into the subwavelength regime. Although circular SAR (CSAR) has a sharply peaked point response in 2D, its side lobes of closely-spaced targets interfere strongly with each other to distort the image. An alternative 2D subwavelength imaging approach with minimal side lobes is Time-Reversal MUSIC (TR-MUSIC). Traditional TR-MUSIC, however, requires a large number of transmitters and receivers. We propose a hybrid synthetic aperture TR-MUSIC algorithm (SATR-MUSIC) that combines the benefits of both approaches. Using relatively few transceivers, SATR-MUSIC is able to resolve objects separated by approximately in 2D with minimal background artifacts. It does so by averaging TR-MUSIC’s imaging kernel incoherently over the synthetic aperture.

Effects of opioids on respiration assessed by a contact-free unconstraint respiratory monitor with load cells under the bed in patients with advanced cancer.

Nocturnal periodic breathing of chronic opioid users has been predominantly documented by the use of polysomnography. No previous studies have assessed the opioid effects of respiratory rhythms throughout the day without the use of physical restraint. We recently developed a contact-free unconstraint vital sign monitoring system with four load cells placed under the bed legs, which allows continuous measurements of respiratory change at the center of gravity on the bed. We aimed to reveal details of the patient's 24-h respiratory status under a monitoring system and to test the hypothesis that respiratory rhythm abnormalities are opioid dose-dependent and worsen during the night time. Continuous 48-h respiratory measurements were successfully performed in 51 patients with advanced cancer (12 opioid-free patients and 39 opioid-receiving patients). Medians of respiratory variables with minimal body movement artifacts were calculated for each 8-h split time period. Compared with opioid-free patients, opioid-receiving patients had slower respiratory rate with higher respiratory rate irregularity without changing tidal centroid shift regardless of the time period. Irregular ataxic breathing was only identified in opioid-receiving patients (33%, P = 0.023) whereas incidence rate of periodic breathing did not differ between the groups. Multivariate regression analyses revealed that opioid dose was an independent risk factor for occurrence of irregular breathing [odds ratio 1.81 (95% CI: 1.39-2.36), P < 0.001], and ataxic breathing [odds ratio 2.08 (95% CI: 1.60-2.71), P < 0.001]. Females developed the ataxic breathing at lower opioid dose compared with males. We conclude that respiratory rhythm irregularity is a predominant specific feature of opioid dose-dependent respiratory depression particularly in female patients with advanced cancer.NEW & NOTEWORTHY Through usage of a novel contact-free unconstraint vital sign monitoring system with four load cells placed under the bed legs allowing continuous measurements of respiratory changes of center of gravity on the bed, this study is the first to assess detailed respiratory characteristics throughout day and night periods without interference of daily activities in patients with advanced cancer receiving opioids. Respiratory rhythm irregularity is a predominant specific feature of opioid dose-dependent respiratory depression particularly in female patients with advanced cancer.

T2*-weighted MRI produces viable fetal "Black-Bone" contrast with significant benefits when compared to current sequences.

Fetal "black bone" MRI could be useful in the diagnosis of various skeletal conditions during pregnancy without exposure to ionizing radiation. Previously suggested susceptibility-weighted imaging (SWI) is not available in the suggested form on all scanners leading to long imaging times that are susceptible to motion artefacts. We aimed to assess if an optimized T2*-weighted GRE sequence can provide viable "black bone" contrast and compared it to other sequences in the literature. A retrospective study was conducted on 17 patients who underwent fetal MRI. Patients were imaged with an optimized T2*-weighted GRE sequence, as well as at least one other "black-bone" sequence. Image quality was scored by four blinded observers on a five-point scale. The T2*-weighted GRE sequence offered adequate to excellent image quality in 63% of cases and scored consistently higher than the three other comparison sequences when comparing images from the same patient. Image quality was found to be dependent on gestational age with good image quality achieved on almost all patients after 26 weeks. T2*-weighted GRE imaging can provide adequate fetal "black bone" contrast and performs at least as well as other sequences in the literature due to good bone to soft tissue contrast and minimal motion artefacts. T2*-weighted fetal "black-bone" imaging can provide excellent bone to soft tissue contrast without using ionizing radiation. It is as good as other "black bone" sequences and may be simpler and more widely implemented, with less motion artefacts.

Development of Feed Forward and Conditional Models Based on Rodent Imaging Data to Predict Effects of Nitrogen Mustard on Pulmonary Function

Nitrogen mustard (NM; (bis(2-chloroethyl)methylamine) is a cytotoxic vesicant known to cause pulmonary damage. In these studies, we developed mathematical models that allow prediction of NM-induced functional deficits in the lung based on structural alterations identified by magnetic resonance imaging (MRI) and computed tomography (CT). Male Wistar rats were exposed to PBS (CTL) or NM (0.125 mg/kg) via i.t. instillation. CT and MRI scans were performed prior to exposure and 28 d post exposure. Pulmonary function was assessed using a SCIREQ flexiVent at a positive end expiratory pressure (PEEP) of 3 cm H2O. Because of accuracy and minimal artifacts, MRI data was used to calculate total lung volume (W). As CT is capable of differentiating between regions of air and water within the lung, it was used as a direct measure of fluid density within the tissue. Relative lung volumes calculated from CT data were used to summarize hyperinflation (I), normal lung, and consolidated lung tissue (B) using voxel density. Using these elements, a feed forward model was constructed to predict respiratory impedance (Z), a measure of airflow, with the equation: z = (A(η-j))/ωα where A = B/(I+W). Calculated values of A from imaging data confirmed alterations in lung structure following NM exposure (1.38 ± 0.08, n=4, CTL vs. 0.99 ± 0.16, n=3, NM) as a result of hyperinflation (I). The predicted Z spectra derived from imaging data was compared to the measured Z spectra obtained from the flexiVent data. A strong correlation was observed between predicted and measured Z spectra at 28 d post NM (R2=0.997 ± 0.01, n=7). Predicted and measured Z spectra were consistent with previous findings that NM caused a loss of ventilatory function at high frequencies. To further evaluate lung function, components of the Z spectra were used to conditionally model resistance, RL = (a+bf)/(c+f) and elastance, EL = E0 + ΔE(1 – е-βf). This model allowed us to differentiate functional contributions of tissue and airway components at low and high frequencies in imaging and flexiVent data sets. Together, our results demonstrate that structural imaging data can be used to assess lung function as a toxicological endpoint following exposure to mustard vesicants. The use of this model will aid in the rapid evaluation of potential countermeasures to mitigate lung toxicity.

The Surgical Treatment of Pyogenic Spondylodiscitis using Carbon-Fiber–Reinforced Polyether Ether Ketone Implants: Personal Experience of a Series of 81 Consecutive Patients

Pyogenic spondylodiscitis (PSD) is a complex disorder that often required postoperative imaging. Carbon-fiber-reinforced polyether ether ketone (CFRP) is radiolucent and offers an optimal assessability of anatomic structures. A retrospective file review of patients who were operated on for PSD using CFRP implants was performed to assess the clinical outcome, implant-associated complications, and revision surgery. A minimum follow-up of 3 months was required for evaluation of clinical and radiographic data, which included computed tomography and magnetic resonance imaging (MRI) assessment, to determine implant stability and assessability of soft tissue and nerve structures using a grading system. Eighty-one consecutive patients with a mean of 69.5 years were identified. Debridement and stabilization were performed in 8 cervical, 17 thoracic, and 57 lumbar procedures; 72 interbody fusion procedures using cages were performed. Intraoperatively, no implant-associated complication was noted. The mean follow-up was 7 months, at which 52 patients attended. Improved mobility and reduced pain levels were reported by 87%, and MRI assessability was graded ideal. Residual sign of infection was seen in 5 cases, which influenced antibiotic therapy. Asymptomatic radiolucent zones were identified in 13 patients (16%) and screw loosening in 2 (2.4%). In 1 patient, the pedicle screw tip broke and remained within the vertebral body. A repeated procedure because of progressive vertebral body destruction, implant loosening, or subsidence was performed in 5 patients (6.1%). The surgical treatment of PSD using CFRP is safe. The repeat procedure rate as a result of implant loosening is 6.1%. Minimal artifacts offer ideal assessability of soft tissue structures on an MRI.

Enhanced UnaG With Minimal Labeling Artifact for Single-Molecule Localization Microscopy.

We introduced enhanced UnaG (eUnaG), a ligand-activatable fluorescent protein, for conventional and super-resolution imaging of subcellular structures in the mammalian cells. eUnaG is a V2L mutant of UnaG with twice brighter bulk fluorescence. We previously discovered the reversible fluorescence switching behavior of UnaG and demonstrated the high photon outputs and high localization numbers in single-molecule localization microscopy (SMLM). In this study, we showed that the fluorescence of eUnaG can be switched off under blue-light illumination, while a high concentration of fluorogenic ligands in the buffer can efficiently restore the fluorescence, as in UnaG. We demonstrated the capacity of eUnaG as an efficient protein label in mammalian cells, as well as for SMLM by utilizing its photoswitchable nature. While cytosolic UnaG proteins showed aggregated patches and fluorescence reduction at high expression levels, eUnaG-labeled protein targets successfully formed their proper structures in mammalian cells without notable distortion from the endogenous structure in the majority of transiently expressing cells. In particular, eUnaG preserved the vimentin filament structures much better than the UnaG. eUnaG provided similarly high single-molecule photon count distribution to UnaG, thus also similarly high resolution in the super-resolution images of various subcellular structures. The sampling coverage analysis of vimentin filaments in SMLM images showed the improvement of labeling efficiency of eUnaG. eUnaG is a high-performance fluorescent protein for fluorescence and single-molecule localization imaging in green emission with minimal labeling artifact.

Fast Single-Image Dehazing Algorithm Based on Piecewise Transformation in Optical Model

Single-image dehazing techniques are extensively used in outdoor optical image acquisition equipment. Most existing methods pay attention to use various priors to estimate scene transmission. In this paper, a fast single-image dehazing algorithm is proposed based on a piecewise transformation model between the minimum channels of the hazy image and the haze-free image in optical model. The minimum channel of the haze-free image is obtained by the piecewise transformation, which is a quadratic function model that we establish for the dark region, and a linear transformation model is established for the bright region. Using the minimum channels of the hazy image and the haze-free image, a transmission estimation model is established based on the haze optical model with adjustment variables. To obtain an accurate estimation of atmospheric light, we estimate the atmospheric light twice. Finally, the haze-free image is restored. Experimental results show that the proposed algorithm has minimal halo artifacts and color distortion in various depths of field, flat areas and sky areas. From the subjective evaluation, objective evaluation and running time analysis, it can be seen that the algorithm in this paper is superior to most existing technologies.

A task-relevant experimental pain model to target motor adaptation.

Motor adaptation is thought to be a strategy to avoid pain. Current experimental pain models do not allow for consistent modulation of pain perception depending on movement. We showed that low-frequency sinusoidal stimuli delivered at painful intensity result in minimal habituation of pain perception (over 60s) and minimal stimulation artefacts on electromyographic signals. When the amplitude of the low-frequency sinusoidal stimuli was modulated based on the vertical force participants applied to the ground with their right leg while standing upright, we demonstrated a strong association between perceived pain and motor adaptation. By enabling task-relevant modulation of perceived pain intensity and the recording electromyographic signals during electrical painful stimulation, our novel pain model will permit direct experimental testing of the relationship between pain and motor adaptation. Contemporary pain adaptation theories predict that motor adaptation occurs to limit pain. Current experimental pain models, however, do not allow for pain intensity modulation according to one's posture or movements. We developed a task-relevant experimental pain model using low-frequency sinusoidal electrical stimuli applied over the infrapatellar fat pad. In fourteen participants, we compared perceived pain habituation and stimulation-induced artefacts in vastus medialis electromyographic recordings elicited by sinusoidal (4, 10, 20 and 50Hz) and square electrical waveforms delivered at constant peak stimulation amplitude. Next, we simulated a clinical condition where perceived knee pain intensity is proportional to the load applied on the leg by controlling sinusoidal current amplitude (4Hz) according to the vertical force the participants applied with their right leg to the ground while standing upright. Pain ratings habituated over a 60s period for 50Hz sinusoidal and square waveforms but not for low-frequency sinusoidal stimuli (P<0.001). EMG filters removed most stimulation artefacts for low-frequency sinusoidal stimuli (4Hz). While balancing upright, participants' pain ratings were correlated with the force applied by the right leg (R2 =0.65), demonstrating task-relevant changes in perceived pain intensity. Low-frequency sinusoidal stimuli can induce knee pain of constant intensity for 60s with minimal EMG artefacts while enabling task-relevant pain modulation when controlling current amplitude. By enabling task-dependent modulation of perceived pain intensity, our novel experimental model replicates key temporal aspects of clinical musculoskeletal pain while allowing quantification of neuromuscular activation during painful electrical stimulation. This approach will enable researchers to test the predicted relationship between movement strategies and pain.

MRI Compatible, Customizable, and 3D-Printable Microdrive for Neuroscience Research.

The effective connectivity of brain networks can be assessed using functional magnetic resonance imaging (fMRI) to quantify the effects of local electrical microstimulation (EM) on distributed neuronal activity. The delivery of EM to specific brain regions, particularly with layer specificity, requires MRI compatible equipment that provides fine control of a stimulating electrode’s position within the brain while minimizing imaging artifacts. To this end, we developed a microdrive made entirely of MRI compatible materials. The microdrive uses an integrated penetration grid to guide electrodes and relies on a microdrilling technique to eliminate the need for large craniotomies, further reducing implant maintenance and image distortions. The penetration grid additionally serves as a built-in MRI marker, providing a visible fiducial reference for estimating probe trajectories. Following the initial implant procedure, these features allow for multiple electrodes to be inserted, removed, and repositioned with minimal effort, using a screw-type actuator. To validate the design of the microdrive, we conducted an EM-coupled fMRI study with a male macaque monkey. The results verified that the microdrive can be used to deliver EM during MRI procedures with minimal imaging artifacts, even within a 7 Tesla (7T) environment. Future applications of the microdrive include neuronal recordings and targeted drug delivery. We provide computer aided design (CAD) templates and a parts list for modifying and fabricating the microdrive for specific research needs. These designs provide a convenient, cost-effective approach to fabricating MRI compatible microdrives for neuroscience research.

Azimuthal averaging–reconstruction filtering techniques for finite-difference general circulation models in spherical geometry

Abstract. When solving hydrodynamic equations in spherical or cylindrical geometry using explicit finite-difference schemes, a major difficulty is that the time step is greatly restricted by the clustering of azimuthal cells near the pole due to the Courant–Friedrichs–Lewy condition. This paper adapts the azimuthal averaging–reconstruction (ring average) technique to finite-difference schemes in order to mitigate the time step constraint in spherical and cylindrical coordinates. The finite-difference ring average technique averages physical quantities based on an effective grid and then reconstructs the solution back to the original grid in a piecewise, monotonic way. The algorithm is implemented in a community upper-atmospheric model, the Thermosphere–Ionosphere Electrodynamics General Circulation Model (TIEGCM), with a horizontal resolution up to 0.625∘×0.625∘ in geographic longitude–latitude coordinates, which enables the capability of resolving critical mesoscale structures within the TIEGCM. Numerical experiments have shown that the ring average technique introduces minimal artifacts in the polar region of general circulation model (GCM) solutions, which is a significant improvement compared to commonly used low-pass filtering techniques such as the fast Fourier transform method. Since the finite-difference adaption of the ring average technique is a post-solver type of algorithm, which requires no changes to the original computational grid and numerical algorithms, it has also been implemented in much more complicated models with extended physical–chemical modules such as the Coupled Magnetosphere–Ionosphere–Thermosphere (CMIT) model and the Whole Atmosphere Community Climate Model with thermosphere and ionosphere eXtension (WACCM-X). The implementation of ring average techniques in both models enables CMIT and WACCM-X to perform global simulations with a much higher resolution than that used in the community versions. The new technique is not only a significant improvement in space weather modeling capability, but it can also be adapted to more general finite-difference solvers for hyperbolic equations in spherical and polar geometries. Highlights. The ring average technique is adapted to solve the issue of clustered grid cells in polar and spherical coordinates with a finite-difference method. The ring average technique is applied to develop a 0.625∘×0.625∘ high-resolution TIEGCM and more complicated geoscientific models with polar and spherical coordinates as well as finite-difference numerical schemes. The high-resolution TIEGCM shows good capability in resolving mesoscale structures in the ionosphere–thermosphere (I–T) system.