Search

We study the problem of dynamically maintaining the connected components of an undirected graph subject to edge insertions and deletions. We give the first parallel algorithm for the problem which is work-efficient, supports batches of updates, runs in polylogarithmic depth, and uses only linear total space. The existing algorithms for the problem either use super-linear space, do not come with strong theoretical bounds, or are not parallel. On the empirical side, we provide the first implementation of the cluster forest algorithm, the first linear-space and poly-logarithmic update time algorithm for dynamic connectivity. Experimentally, we find that our algorithm uses up to 19.7x less space and is up to 6.2x faster than the level-set algorithm of HDT, arguably the most widely-implemented dynamic connectivity algorithm with strong theoretical guarantees., Comment: This is a preliminary version of a paper that will appear at VLDB'25

CT is a main modality for imaging liver diseases, valuable in detecting and localizing liver tumors. Traditional anomaly detection methods analyze reconstructed images to identify pathological structures. However, these methods may produce suboptimal results, overlooking subtle differences among various tissue types. To address this challenge, here we employ generative diffusion prior to inpaint the liver as the reference facilitating anomaly detection. Specifically, we use an adaptive threshold to extract a mask of abnormal regions, which are then inpainted using a diffusion prior to calculating an anomaly score based on the discrepancy between the original CT image and the inpainted counterpart. Our methodology has been tested on two liver CT datasets, demonstrating a significant improvement in detection accuracy, with a 7.9% boost in the area under the curve (AUC) compared to the state-of-the-art. This performance gain underscores the potential of our approach to refine the radiological assessment of liver diseases., Comment: 10 pages, 11 figures

Ultra-high resolution images are desirable in photon counting CT (PCCT), but resolution is physically limited by interactions such as charge sharing. Deep learning is a possible method for super-resolution (SR), but sourcing paired training data that adequately models the target task is difficult. Additionally, SR algorithms can distort noise texture, which is an important in many clinical diagnostic scenarios. Here, we train conditional denoising diffusion probabilistic models (DDPMs) for PCCT super-resolution, with the objective to retain textural characteristics of local noise. PCCT simulation methods are used to synthesize realistic resolution degradation. To preserve noise texture, we explore decoupling the noise and signal image inputs and outputs via deep denoisers, explicitly mapping to each during the SR process. Our experimental results indicate that our DDPM trained on simulated data can improve sharpness in real PCCT images. Additionally, the disentanglement of noise from the original image allows our model more faithfully preserve noise texture., Comment: 5 pages, 4 figures

While drastic improvements in CT technology have occurred in the past 25 years, spatial resolution is one area where progress has been limited until recently. New photon counting CT systems, are capable of much better spatial resolution than their (energy integrating) predecessors. These improvements have the potential to improve the evaluation obstructive coronary artery disease by enabling more accurate delineation between calcified plaque and coronary vessel lumen. A new set of vessel phantoms has been designed and manufactured for quantifying this improvement. Comparisons are made between an existing clinical CT system, a prototype photon counting system, with images from a micro CT system being used as the gold standard. Scans were made of the same objects on all three systems. The resulting images were registered and the luminal cross section areas were compared. Luminal cross-sections near calcified plaques were reduced due to blooming, but this effect was much less pronounced in images from the prototype photon counting system as compared to the images from the clinical CT system.

Recent development of photon-counting CT (PCCT) brings great opportunities for plaque characterization with much-improved spatial resolution and spectral imaging capability. While existing coronary plaque PCCT imaging results are based on detectors made of CZT or CdTe materials, deep-silicon photon-counting detectors have unique performance characteristics and promise distinct imaging capabilities. In this work, we report a systematic simulation study of a deep-silicon PCCT scanner with a new clinically-relevant digital plaque phantom with realistic geometrical parameters and chemical compositions. This work investigates the effects of spatial resolution, noise, motion artifacts, radiation dose, and spectral characterization. Our simulation results suggest that the deep-silicon PCCT design provides adequate spatial resolution for visualizing a necrotic core and quantitation of key plaque features. Advanced denoising techniques and aggressive bowtie filter designs can keep image noise to acceptable levels at this resolution while keeping radiation dose comparable to that of a conventional CT scan. The ultrahigh resolution of PCCT also means an elevated sensitivity to motion artifacts. It is found that a tolerance of less than 0.4 mm residual movement range requires the application of accurate motion correction methods for best plaque imaging quality with PCCT., Comment: 13 figures, 5 tables

Radiomics uses quantitative medical imaging features to predict clinical outcomes. Currently, in a new clinical application, finding the optimal radiomics method out of the wide range of available options has to be done manually through a heuristic trial-and-error process. In this study we propose a framework for automatically optimizing the construction of radiomics workflows per application. To this end, we formulate radiomics as a modular workflow and include a large collection of common algorithms for each component. To optimize the workflow per application, we employ automated machine learning using a random search and ensembling. We evaluate our method in twelve different clinical applications, resulting in the following area under the curves: 1) liposarcoma (0.83); 2) desmoid-type fibromatosis (0.82); 3) primary liver tumors (0.80); 4) gastrointestinal stromal tumors (0.77); 5) colorectal liver metastases (0.61); 6) melanoma metastases (0.45); 7) hepatocellular carcinoma (0.75); 8) mesenteric fibrosis (0.80); 9) prostate cancer (0.72); 10) glioma (0.71); 11) Alzheimer's disease (0.87); and 12) head and neck cancer (0.84). We show that our framework has a competitive performance compared human experts, outperforms a radiomics baseline, and performs similar or superior to Bayesian optimization and more advanced ensemble approaches. Concluding, our method fully automatically optimizes the construction of radiomics workflows, thereby streamlining the search for radiomics biomarkers in new applications. To facilitate reproducibility and future research, we publicly release six datasets, the software implementation of our framework, and the code to reproduce this study., Comment: 33 pages, 4 figures, 4 tables, 2 supplementary figures, 3 supplementary table, submitted to Medical Image Analysis; revision

In clinical CT system, the x-ray tube emits polychromatic x-rays, and the x-ray detectors operate in the current-integrating mode. This physical process is accurately described by an energy-dependent non-linear integral equation. However, the non-linear model is not invertible with a computationally efficient solution, and is often approximated as a linear integral model in the form of the Radon transform. Such approximation basically ignores energy-dependent information and would generate beam hardening artifacts. Dual-energy CT (DECT) scans one object using two different x-ray energy spectra for the acquisition of two spectrally distinct projection datasets to improve imaging performance. Thus, DECT can reconstruct energy and material-selective images, realizing monochromatic imaging and material decomposition. Nevertheless, DECT would increase radiation dose, system complexity, and equipment cost relative to single-spectrum CT. In this paper, a machine-learning-based CT reconstruction method is proposed to perform monochromatic image reconstruction using a single-spectrum CT scanner. Specifically, a residual neural network (ResNet) model is adapted to map a CT image to a monochromatic counterpart at a pre-specified energy level. This ResNet is trained on clinical dual-energy data, showing an excellent convergence to a minimal loss. The trained network produces high-quality monochromatic images on testing data, with a relative error of less than 0.2%. This work has great potential in clinical DECT applications such as tissue characterization, beam hardening correction and proton therapy planning.

Deep learning (DL) has shown unprecedented performance for many image analysis and image enhancement tasks. Yet, solving large-scale inverse problems like tomographic reconstruction remains challenging for DL. These problems involve non-local and space-variant integral transforms between the input and output domains, for which no efficient neural network models have been found. A prior attempt to solve such problems with supervised learning relied on a brute-force fully connected network and applied it to reconstruction for a $128^4$ system matrix size. This cannot practically scale to realistic data sizes such as $512^4$ and $512^6$ for three-dimensional data sets. Here we present a novel framework to solve such problems with deep learning by casting the original problem as a continuum of intermediate representations between the input and output data. The original problem is broken down into a sequence of simpler transformations that can be well mapped onto an efficient hierarchical network architecture, with exponentially fewer parameters than a generic network would need. We applied the approach to computed tomography (CT) image reconstruction for a $512^4$ system matrix size. To our knowledge, this enabled the first data-driven DL solver for full-size CT reconstruction without relying on the structure of direct (analytical) or iterative (numerical) inversion techniques. The proposed approach is applicable to other imaging problems such as emission and magnetic resonance reconstruction. More broadly, hierarchical DL opens the door to a new class of solvers for general inverse problems, which could potentially lead to improved signal-to-noise ratio, spatial resolution and computational efficiency in various areas.

Known for its high morbidity and mortality rates, lung cancer poses a significant threat to human health and well-being. However, the same population is also at high risk for other deadly diseases, such as cardiovascular disease. Since Low-Dose CT (LDCT) has been shown to significantly improve the lung cancer diagnosis accuracy, it will be very useful for clinical practice to predict the all-cause mortality for lung cancer patients to take corresponding actions. In this paper, we propose a deep learning based method, which takes both chest LDCT image patches and coronary artery calcification risk scores as input, for direct prediction of mortality risk of lung cancer subjects. The proposed method is called Hybrid Risk Network (HyRiskNet) for mortality risk prediction, which is an end-to-end framework utilizing hybrid imaging features, instead of completely relying on automatic feature extraction. Our work demonstrates the feasibility of using deep learning techniques for all-cause lung cancer mortality prediction from chest LDCT images. The experimental results show that the proposed HyRiskNet can achieve superior performance compared with the neural networks with only image input and with other traditional semi-automatic scoring methods. The study also indicates that radiologist defined features can well complement convolutional neural networks for more comprehensive feature extraction., Comment: IEEE conference format

Hybrid quantum systems have been developed with various mechanical, optical and microwave harmonic oscillators. The coupling produces a rich library of interactions including two mode squeezing, swapping interactions, back-action evasion and thermal control. In a multimode mechanical system, coupling resonators of different scales (both in frequency and mass) leverages the advantages of each resonance. For example: a high frequency, easily manipulated resonator could be entangled with a low frequency massive object for tests of gravitational decoherence. Here we demonstrate coherent optomechanical state swapping between two spatially and frequency separated resonators with a mass ratio of 4. We find that, by using two laser beams far detuned from an optical cavity resonance, efficient state transfer is possible through a process very similar to STIRAP (Stimulated Raman Adiabatic Passage) in atomic physics. Although the demonstration is classical, the same technique can be used to generate entanglement between oscillators in the quantum regime.