Coil-combined split slice-GRAPPA for simultaneous multi-slice diffusion MRI

Objective: To develop a kernel optimization method called coil-combined split slice-GRAPPA (CC-SSG) to improve the accuracy of the reconstructed coil-combined images for simultaneous multi-slice (SMS) diffusion weighted imaging (DWI) data.Methods: The CC-SSG method optimizes the tuning parameters in the k-space SSG kernels to achieve an optimal trade-off between the intra-slice artifact and inter-slice leakage after the root-sum-of-squares (rSOS) coil combining of the de-aliased SMS DWI data. A detailed analysis is conducted to evaluate the contributions of the intra-slice artifact and inter-slice leakage to the total reconstruction error after coil combining.Results: Comparisons of the proposed CC-SSG method with the slice-GRAPPA (SG) and split slice-GRAPPA (SSG) methods are provided using two in-vivo readout-segmented (RS) EPI datasets collected from stroke patients. The CC-SSG method demonstrates improved accuracy of the reconstructed coil-combined images and the estimated diffusion tensor imaging (DTI) maps.Conclusion: CC-SSG strikes a good balance between the intra-slice artifact and inter-slice leakage for rSOS coil combining, and so can yield better reconstruction performance compared to SG and SSG for rSOS reconstruction. The optimal trade-off between the two artifacts is robust to the contrast of SMS data and the choice of the coil combining method.     

Deep learning for radial SMS myocardial perfusion reconstruction using the 3D residual booster U-net

PurposeTo develop an end-to-end deep learning solution for quickly reconstructing radial simultaneous multi-slice (SMS) myocardial perfusion datasets with comparable quality to the pixel tracking spatiotemporal constrained reconstruction (PT-STCR) method.MethodsDynamic contrast enhanced (DCE) radial SMS myocardial perfusion data were obtained from 20 subjects who were scanned at rest and/or stress with or without ECG gating using a saturation recovery radial CAIPI turboFLASH sequence. Input to the networks consisted of complex coil combined images reconstructed using the inverse Fourier transform of undersampled radial SMS k-space data. Ground truth images were reconstructed using the PT-STCR pipeline. The performance of the residual booster 3D U-Net was tested by comparing it to state-of-the-art network architectures including MoDL, CRNN-MRI, and other U-Net variants.ResultsResults demonstrate significant improvements in speed requiring approximately 8 seconds to reconstruct one radial SMS dataset which is approximately 200 times faster than the PT-STCR method. Images reconstructed with the residual booster 3D U-Net retain quality of ground truth PT-STCR images (0.963 SSIM/40.238 PSNR/0.147 NRMSE). The residual booster 3D U-Net has superior performance compared to existing network architectures in terms of image quality, temporal dynamics, and reconstruction time.ConclusionResidual and booster learning combined with the 3D U-Net architecture was shown to be an effective network for reconstructing high-quality images from undersampled radial SMS datasets while bypassing the reconstruction time of the PT-STCR method.     

Data consistency criterion for selecting parameters for k-space-based reconstruction in parallel imaging

k-space-based reconstruction in parallel imaging depends on the reconstruction kernel setting, including its support. An optimal choice of the kernel depends on the calibration data, coil geometry and signal-to-noise ratio, as well as the criterion used. In this work, data consistency, imposed by the shift invariance requirement of the kernel, is introduced as a goodness measure of k-space-based reconstruction in parallel imaging and demonstrated. Data consistency error (DCE) is calculated as the sum of squared difference between the acquired signals and their estimates obtained based on the interpolation of the estimated missing data. A resemblance between DCE and the mean square error in the reconstructed image was found, demonstrating DCE's potential as a metric for comparing or choosing reconstructions. When used for selecting the kernel support for generalized autocalibrating partially parallel acquisition (GRAPPA) reconstruction and the set of frames for calibration as well as the kernel support in temporal GRAPPA reconstruction, DCE led to improved images over existing methods. Data consistency error is efficient to evaluate, robust for selecting reconstruction parameters and suitable for characterizing and optimizing k-space-based reconstruction in parallel imaging.     

A software channel compression technique for faster reconstruction with many channels

In magnetic resonance imaging, highly parallel imaging using coil arrays with a large number of elements is an area of growing interest. With increasing channel numbers for parallel acquisition, the increased reconstruction time and extensive computer memory requirements have become significant concerns. In this work, principal component analysis (PCA) is used to develop a channel compression technique. This technique efficiently reduces the size of parallel imaging data acquired from a multichannel coil array, thereby significantly reducing the reconstruction time and computer memory requirement without undermining the benefits of multichannel coil arrays. Clinical data collected with a 32-channel cardiac coil are used in all of the experiments. The performance of the proposed method on parallel, partially acquired data, as well as fully acquired data, was evaluated. Experimental results show that the proposed method dramatically reduces the processing time without considerable degradation in the quality of reconstructed images. It is also demonstrated that this PCA technique can be used to perform intensity correction in parallel imaging applications.     

Parallel magnetic resonance imaging using coils with localized sensitivities

The purpose of this study was to present clinical examples and illustrate the inefficiencies of a conventional reconstruction using a commercially available phased array coil with localized sensitivities. Five patients were imaged at 1.5 T using a cardiac-synchronized gadolinium-enhanced acquisition and a commercially available four-element phased array coil. Four unique sets of images were reconstructed from the acquired k-space data: (a) sum-of-squares image using four elements of the coil; localized sum-of-squares images from the (b) anterior coils and (c) posterior coils and a (c) local reconstruction. Images were analyzed for artifacts and usable field-of-view. Conventional image reconstruction produced images with fold-over artifacts in all cases spanning a portion of the image (mean 90 mm; range 36-126 mm). The local reconstruction removed fold-over artifacts and resulted in an effective increase in the field-of-view (mean 50%; range 20-70%). Commercially available phased array coils do not always have overlapping sensitivities. Fold-over artifacts can be removed using an alternate reconstruction method. When assessing the advantages of parallel imaging techniques, gains achieved using techniques such as SENSE and SMASH should be gauged against the acquisition time of the localized method rather than the conventional sum-of-squares method.     

APIR4EMC: Autocalibrated parallel imaging reconstruction for extended multi-contrast imaging

PurposeTo improve image quality of multi-contrast imaging with the proposed Autocalibrated Parallel Imaging Reconstruction for Extended Multi-Contrast Imaging (APIR4EMC).MethodsAPIR4EMC reconstructs multi-contrast images in an autocalibrated parallel imaging reconstruction framework by adding contrasts as virtual coils. Compensation of signal evolution along the echo train of different contrasts is performed to improve signal prediction for missing samples. As a proof of concept, we performed prospectively accelerated phantom and in-vivo brain acquisitions with T1, T1-fat saturated (Fatsat), T2, PD, and FLAIR contrasts. The k-space sampling patterns of these acquisitions were jointly optimized. Images were jointly reconstructed with the proposed APIR4EMC method as well as individually with GRAPPA. Root mean square error (RMSE) to fully sampled reference images and g-factor maps were computed for both methods in the phantom experiment. Visual evaluation was performed in the in-vivo experiment.ResultsCompared to GRAPPA, APIR4EMC reduced artifacts and improved SNR of the reconstructed images in the phantom acquisitions. Quantitatively, APIR4EMC substantially reduced noise amplification (g-factor) as well as RMSE compared to GRAPPA. Signal evolution compensation reduced artifacts. In the in-vivo experiments, 1 mm³ isotropic 3D images with contrasts of T1, T1-Fatsat, T2, PD, and FLAIR were acquired in as little as 7.5 min with the acceleration factor of 9. Reconstruction quality was consistent with the phantom results.ConclusionCompared to single contrast reconstruction with GRAPPA, APIR4EMC reduces artifacts and noise amplification in accelerated multi-contrast imaging.     

A 32-channel coil system for MR vessel wall imaging of intracranial and extracranial arteries at 3T

PurposeTo develop a RF coil system for joint imaging of intracranial and extracranial arterial vessel wall at 3T.Materials and methodThe coil system consists of a 24-channel head coil combined with an 8-channel carotid coil. It is compared with a standard coil configuration (12-channel head coil + 4-channel neck coil + 8-channel carotid coil) for SNR and g-factors in phantoms and healthy volunteers. The clinical relevance of the proposed coil system is also evaluated in patients.ResultsIn phantom experiments, the SNR of the proposed coil system is 53% higher than the maximum SNR of the standard coil configuration at the center of the phantom which usually corresponds to the intracranial region of the head. The g-factors of the proposed coil system in the sagittal plane are lower than the standard coil configuration (by 10.8% and 26.6% for R = 2 and 4 respectively) in the same experiment. In healthy volunteer experiments, 55% of the pixels have SNR above 100 for the proposed coil system, which is 33% more than that of the standard coil configuration. The maximum g-factors in the standard configuration are higher than those from the new coil design by 12% at R = 2 and up to 36% at R = 4 in the sagittal plane. In patients, in-vivo intracranial and extracranial arterial wall images at an isotropic spatial resolution of 0.6 mm can be acquired using the proposed coil system. Plaques are well depicted from the images.ConclusionsThe performance of the proposed coil set is superior to the standard coil configuration, providing high SNR, low g-factor and good spatial coverage needed for simultaneous high resolution imaging of intracranial and extracranial arterial walls. Images acquired in 7.6 min using the proposed coil system can achieve an isotropic spatial resolution of 0.6 mm and can be used to depict plaques on the intracranial and extracranial arterial walls in patients.     

Retrospective correction of bias in diffusion tensor imaging arising from coil combination mode

PurposeTo quantify and retrospectively correct for systematic differences in diffusion tensor imaging (DTI) measurements due to differences in coil combination mode.BackgroundMulti-channel coils are now standard among MRI systems. There are several options for combining signal from multiple coils during image reconstruction, including sum-of-squares (SOS) and adaptive combine (AC). This contribution examines the bias between SOS- and AC-derived measures of tissue microstructure and a strategy for limiting that bias.MethodsFive healthy subjects were scanned under an institutional review board-approved protocol. Each set of raw image data was reconstructed twice—once with SOS and once with AC. The diffusion tensor was calculated from SOS- and AC-derived data by two algorithms—standard log-linear least squares and an approach that accounts for the impact of coil combination on signal statistics. Systematic differences between SOS and AC in terms of tissue microstructure (axial diffusivity, radial diffusivity, mean diffusivity and fractional anisotropy) were evaluated on a voxel-by-voxel basis.ResultsSOS-based tissue microstructure values are systematically lower than AC-based measures throughout the brain in each subject when using the standard tensor calculation method. The difference between SOS and AC can be virtually eliminated by taking into account the signal statistics associated with coil combination.ConclusionsThe impact of coil combination mode on diffusion tensor-based measures of tissue microstructure is statistically significant but can be corrected retrospectively. The ability to do so is expected to facilitate pooling of data among imaging protocols.     

Improving dynamic contrast-enhanced MRI of the lung using motion-weighted sparse reconstruction: Initial experiences in patients

PurposeThe purpose of this study was to evaluate the performance of motion-weighted Golden-angle RAdial Sparse Parallel MRI (motion-weighted GRASP) for free-breathing dynamic contrast-enhanced MRI (DCE-MRI) of the lung.MethodsMotion-weighted GRASP incorporates a soft-gating motion compensation algorithm into standard GRASP reconstruction, so that motion-corrupted motion k-space (e.g., k-space acquired in inspiratory phases) contributes less to the final reconstructed images. Lung MR data from 20 patients (mean age = 57.9 ± 13.5) with known pulmonary lesions were retrospectively collected for this study. Each subject underwent a free-breathing DCE-MR scan using a fat-statured T₁-weighted stack-of-stars golden-angle radial sequence and a post-contrast breath-hold MR scan using a Cartesian volumetric-interpolated imaging sequence (BH-VIBE). Each radial dataset was reconstructed using GRASP without motion compensation and motion-weighted GRASP. All MR images were visually evaluated by two experienced radiologists blinded to reconstruction and acquisition schemes independently. In addition, the influence of motion-weighted reconstruction on dynamic contrast-enhancement patterns was also investigated.ResultsFor image quality assessment, motion-weighted GRASP received significantly higher visual scores than GRASP (P < 0.05) for overall image quality (3.68 vs. 3.39), lesion conspicuity (3.54 vs. 3.18) and overall artifact level (3.53 vs. 3.15). There was no significant difference (P > 0.05) between the breath-hold BH-VIBE and motion-weighted GRASP images. For assessment of temporal fidelity, motion-weighted GRASP maintained a good agreement with respect to GRASP.ConclusionMotion-weighted GRASP achieved better reconstruction performance in free-breathing DCE-MRI of the lung compared to standard GRASP, and it may enable improved assessment of pulmonary lesions.     

Optimizing Image Reconstruction in SENSE Using GPU

Parallel magnetic resonance imaging (MRI) (pMRI) uses multiple receiver coils to reduce the MRI scan time. To accelerate the data acquisition process in MRI, less amount of data is acquired from the scanner which leads to artifacts in the reconstructed images. SENSitivity Encoding (SENSE) is a reconstruction algorithm in pMRI to remove aliasing artifacts from the undersampled multi coil data and recovers fully sampled images. The main limitation of SENSE is computing inverse of the encoding matrix. This work proposes the inversion of encoding matrix using Jacobi singular value decomposition (SVD) algorithm for image reconstruction on GPUs to accelerate the reconstruction process. The performance of Jacobi SVD is compared with Gauss–Jordan algorithm. The simulations are performed on two datasets (brain and cardiac) with acceleration factors 2, 4, 6 and 8. The results show that the graphics processing unit (GPU) provides a speed up to 21.6 times as compared to CPU reconstruction. Jacobi SVD algorithm performs better in terms of acceleration in reconstructions on GPUs as compared to Gauss–Jordan method. The proposed algorithm is suitable for any number of coils and acceleration factors for SENSE reconstruction on real time processing systems.     

FPGA implementation of real-time SENSE reconstruction using pre-scan and Emaps sensitivities

Sensitivity Encoding (SENSE) is a widely used technique in Parallel Magnetic Resonance Imaging (MRI) to reduce scan time. Reconfigurable hardware based architecture for SENSE can potentially provide image reconstruction with much less computation time. Application specific hardware platform for SENSE may dramatically increase the power efficiency of the system and can decrease the execution time to obtain MR images. A new implementation of SENSE on Field Programmable Gate Array (FPGA) is presented in this study, which provides real-time SENSE reconstruction right on the receiver coil data acquisition system with no need to transfer the raw data to the MRI server, thereby minimizing the transmission noise and memory usage. The proposed SENSE architecture can reconstruct MR images using receiver coil sensitivity maps obtained using pre-scan and eigenvector (E-maps) methods. The results show that the proposed system consumes remarkably less computation time for SENSE reconstruction, i.e., 0.164 ms @ 200 MHz, while maintaining the quality of the reconstructed images with good mean SNR (29 + dB), less RMSE (< 5 × 10^− 2) and comparable artefact power (< 9 × 10^− 4) to conventional SENSE reconstruction. A comparison of the center line profiles of the reconstructed and reference images also indicates a good quality of the reconstructed images. Furthermore, the results indicate that the proposed architectural design can prove to be a significant tool for SENSE reconstruction in modern MRI scanners and its low power consumption feature can be remarkable for portable MRI scanners.     

Rapid MR imaging by sensitivity profile indexing and deconvolution reconstruction (SPID)

A new parallel MR imaging technique, which uses localized information from the elements of a multi-coil array to accelerate imaging, is described. The technique offers an alternative reconstruction approach to currently available techniques (e.g., SMASH and SENSE). Following a partial k-space data acquisition, image reconstruction in this approach proceeds in two steps: first, fitting the measured coil sensitivities to a set of partially localized target functions, a blurred intermediate image of the studied object is produced. Blurring is obtained in a systematic manner, forming images of the studied object convolved with a known convolution kernel. Full spatial resolution is then recovered by deconvolution of the blurred images with the known kernel function. The technique offers flexibility in the arrangement of the acquired signal data k-lines, and a mechanism for controlling reconstruction quality through the convolution the deconvolution procedure. The technique was validated in phantom and in vivo imaging experiments demonstrating high time reduction factors.     

GPU based parallel framework for receiver coil sensitivity estimation in SENSE reconstruction

Magnetic Resonance Imaging (MRI) uses non-ionizing radiations and is safer as compared to CT and X-ray imaging. MRI is broadly used around the globe for medical diagnostics. One main limitation of MRI is its long data acquisition time. Parallel MRI (pMRI) was introduced in late 1990's to reduce the MRI data acquisition time. In pMRI, data is acquired by under-sampling the Phase Encoding (PE) steps which introduces aliasing artefacts in the MR images. SENSitivity Encoding (SENSE) is a pMRI based method that reconstructs fully sampled MR image from the acquired under-sampled data using the sensitivity information of receiver coils. In SENSE, precise estimation of the receiver coil sensitivity maps is vital to obtain good quality images. Eigen-value method (a recently proposed method in literature for the estimation of receiver coil sensitivity information) does not require a pre-scan image unlike other conventional methods of sensitivity estimation. However, Eigen-value method is computationally intensive and takes a significant amount of time to estimate the receiver coil sensitivity maps. This work proposes a parallel framework for Eigen-value method of receiver coil sensitivity estimation that exploits its inherent parallelism using Graphics Processing Units (GPUs). We evaluated the performance of the proposed algorithm on in-vivo and simulated MRI datasets (i.e. human head and simulated phantom datasets) with Peak Signal-to-Noise Ratio (PSNR) and Artefact Power (AP) as evaluation metrics. The results show that the proposed GPU implementation reduces the execution time of Eigen-value method of receiver coil sensitivity estimation (providing up to 30 times speed up in our experiments) without degrading the quality of the reconstructed image.     

Dual-domain cascade of U-nets for multi-channel magnetic resonance image reconstruction

The U-net is a deep-learning network model that has been used to solve a number of inverse problems. In this work, the concatenation of two-element U-nets, termed the W-net, operating in k-space (K) and image (I) domains, were evaluated for multi-channel magnetic resonance (MR) image reconstruction. The two-element network combinations were evaluated for the four possible image-k-space domain configurations: a) W-net II, b) W-net KK, c) W-net IK, and d) W-net KI. Selected four element (WW-nets) and six element (WWW-nets) networks were also examined. Two configurations of each network were compared: 1) each coil channel was processed independently, and 2) all channels were processed simultaneously. One hundred and eleven volumetric, T1-weighted, 12-channel coil k-space datasets were used in the experiments. Normalized root mean squared error, peak signal-to-noise ratio and visual information fidelity were used to assess the reconstructed images against the fully sampled reference images. Our results indicated that networks that operate solely in the image domain were better when independently processing individual channels of multi-channel data. Dual-domain methods were better when simultaneously reconstructing all channels of multi-channel data. In addition, the best cascade of U-nets performed better (p < 0.01) than the previously published, state-of-the-art Deep Cascade and Hybrid Cascade models in three out of four experiments.     

Accelerated electron paramagnetic resonance imaging using partial Fourier compressed sensing reconstruction

PurposeElectron paramagnetic resonance (EPR) imaging has evolved as a promising tool to provide non-invasive assessment of tissue oxygenation levels. Due to the extremely short T₂ relaxation time of electrons, single point imaging (SPI) is used in EPRI, limiting achievable spatial and temporal resolution. This presents a problem when attempting to measure changes in hypoxic state. In order to capture oxygen variation in hypoxic tissues and localize cycling hypoxia regions, an accelerated EPRI imaging method with minimal loss of information is needed.MethodsWe present an image acceleration technique, partial Fourier compressed sensing (PFCS), that combines compressed sensing (CS) and partial Fourier reconstruction. PFCS augments the original CS equation using conjugate symmetry information for missing measurements. To further improve image quality in order to reconstruct low-resolution EPRI images, a projection onto convex sets (POCS)-based phase map and a spherical-sampling mask are used in the reconstruction process. The PFCS technique was used in phantoms and in vivo SCC7 tumor mice to evaluate image quality and accuracy in estimating O₂ concentration.ResultsIn both phantom and in vivo experiments, PFCS demonstrated the ability to reconstruct images more accurately with at least a 4-fold acceleration compared to traditional CS. Meanwhile, PFCS is able to better preserve the distinct spatial pattern in a phantom with a spatial resolution of 0.6 mm. On phantoms containing Oxo63 solution with different oxygen concentrations, PFCS reconstructed linewidth maps that were discriminative of different O₂ concentrations. Moreover, PFCS reconstruction of partially sampled data provided a better discrimination of hypoxic and oxygenated regions in a leg tumor compared to traditional CS reconstructed images.ConclusionsEPR images with an acceleration factor of four are feasible using PFCS with reasonable assessment of tissue oxygenation. The technique can greatly enhance EPR applications and improve our understanding cycling hypoxia. Moreover this technique can be easily extended to various MRI applications.     

Channel-combination method for phase-based |B1+| mapping techniques

PurposeThe aim of this study was to propose a channel combination method for |B₁⁺| mapping methods using phase difference to reconstruct |B₁⁺| map.Theory and methodsPhase-based |B₁⁺| mapping methods commonly consider the phase difference of two scans to measure |B₁⁺|. Multiple receiver coils acquire a number of images and the phase difference at each channel is theoretically the same in the absence of noise. Affected by noise, phase difference is approximately governed by Gaussian distribution. Considering data from all channels as samples, estimation can be achieved by maximum likelihood method. With this method, all phase differences at each channel are combined into one. In this study, the proposed method is applied with Bloch-Siegert shift |B₁⁺| mapping method. Simulations are performed to illustrate the phase difference distribution and demonstrate the feasibility and facility of the proposed method. Phantom and vivo experiments are carried out at 1.5 T scanner equipped with 8-channel receiver coil. In all experiments, the proposed method is compared with weighted averaging (WA) method.ResultsSimulations revealed appropriateness of approximating the distribution of phase difference to Gaussian distribution. Compared with WA method, the proposed method reduces errors of |B₁⁺| calculation. Phantom and vivo experiments provide further validation.ConclusionConsidering phase noise distribution, the proposed method achieves channel combination by finding the estimation from data acquired by multiple receivers coil. The proposed method reduces |B₁⁺| reconstruction errors caused by noise.     

Joint Reconstruction of Multi-contrast Images and Multi-channel Coil Sensitivities

Magnetic resonance imaging (MRI) has an important feature that it provides multiple images with different contrasts for complementary diagnostic information. However, a large amount of data is needed for multi-contrast images depiction, and thus, the scan is time-consuming. Many methods based on parallel magnetic resonance imaging (pMRI) and compressed sensing (CS) are applied to accelerate multi-contrast MR imaging. Nevertheless, the image reconstructed by sophisticated pMRI methods contains residual aliasing artifact that degrades the quality of the image when the acceleration factor is high. Other methods based on CS always suffer the regularization parameter-selecting problem. To address these issues, a new method is presented for joint multi-contrast image reconstruction and coil sensitivity estimation. The coil sensitivities can be shared during the reconstruction due to the identity of coil sensitivity profiles of different contrast images for imaging stationary tissues. The proposed method uses the coil sensitivities as sharable information during the reconstruction to improve the reconstruction quality. As a result, the residual aliasing artifact can be effectively removed in the reconstructed multi-contrast images even if the acceleration factor is high. Besides, as there is no regularization term in the proposed method, the troublesome regularization parameter selection in the CS can also be avoided. Results from multi-contrast in vivo experiments demonstrated that multi-contrast images can be jointly reconstructed by the proposed method with effective removal of the residual aliasing artifact at a high acceleration factor.     

Optimization of an 8-Channel Loop-Array Coil for a 7 T MRI System with the Guidance of a Co-Simulation Approach

Minimizing coupling between coil elements is technically challenging in designing large-sized, volume-type phased-array coils for human head imaging at ultrahigh fields, e.g., 7 T. As a widely used decoupling method, the capacitive decoupling method has shown excellent performance for loop array. However, building a multi-channel loop array with capacitive decoupling method is laborious that tuning frequency and matching of one coil element will affect adjacent elements and even next adjacent elements. In this study, we made an 8-channel loop-array transmit/receive radio-frequency coil on a 7 T magnetic resonance imaging system with the guidance of frequency domain three-dimensional electromagnetic and radio-frequency circuit co-simulation. The position of decoupling capacitors was investigated and values of all capacitors were predicted from co-simulation. The co-simulation approach cost about 2 days and the error of the predicted and practical capacitance was <5 %. To demonstrate the accuracy of simulation, we evaluated the simulated and measured S-parameter matrixes and B ₁ ⁺ profiles in a birdcage-like excitation mode on a cylindrical water phantom. In addition, B ₁ ⁺ maps and images of human head were shown with the fabricated coil. To demonstrate the parallel imaging performance of this coil array, GRE images using GRAPPA acceleration with the reduction factor R of 1, 2, 3, and 4 were acquired.     

光子计数集成成像的压缩与重构

集成成像技术作为一种重要的裸眼三维显示技术,在完整记录三维场景信息的同时,庞大的数据量给传输和存储带来了压力。为了实现图像的有效压缩和重构,根据光子计数集成成像的特点,基于分布式压缩感知理论,提出用于图像压缩与重构的方案。该方案将图像分为参考图像和非参考图像两类,对其设置不同的测量率并分别进行重构。为保证非参考图像的重构质量,提出一种联合重构算法。该算法首先对非参考图像进行分块测量,依据与参考图像之间的相关性进行图像块分类,然后结合参考图像测量值信息构建新的测量矢量,利用新的测量矢量完成初次图像重构。为了进一步提升图像重构质量,对初次重构结果进行二次残差补偿重构,获得最终重构结果。最后通过设置不同的测量率进行了大量实验,实验结果表明,所提算法在测量率为0.25时,图像重构质量可以达到30 dB,测量率为0.4时,图像质量可以达到35 dB,算法性能具有一定的优越性。