A Bayesian network reconstruction method based on norm minimization is proposed to address the sparsity and iterative divergence issues in network reconstruction caused by noise and missing values.This method achieves...A Bayesian network reconstruction method based on norm minimization is proposed to address the sparsity and iterative divergence issues in network reconstruction caused by noise and missing values.This method achieves precise adjustment of the network structure by constructing a preliminary random network model and introducing small-world network characteristics and combines L1 norm minimization regularization techniques to control model complexity and optimize the inference process of variable dependencies.In the experiment of game network reconstruction,when the success rate of the L1 norm minimization model’s existence connection reconstruction reaches 100%,the minimum data required is about 40%,while the minimum data required for a sparse Bayesian learning network is about 45%.In terms of operational efficiency,the running time for minimizing the L1 normis basically maintained at 1.0 s,while the success rate of connection reconstruction increases significantly with an increase in data volume,reaching a maximum of 13.2 s.Meanwhile,in the case of a signal-to-noise ratio of 10 dB,the L1 model achieves a 100% success rate in the reconstruction of existing connections,while the sparse Bayesian network had the highest success rate of 90% in the reconstruction of non-existent connections.In the analysis of actual cases,the maximum lift and drop track of the research method is 0.08 m.The mean square error is 5.74 cm^(2).The results indicate that this norm minimization-based method has good performance in data efficiency and model stability,effectively reducing the impact of outliers on the reconstruction results to more accurately reflect the actual situation.展开更多
Data reconstruction is a crucial step in seismic data preprocessing.To improve reconstruction speed and save memory,the commonly used three-dimensional(3D)seismic data reconstruction method divides the missing data in...Data reconstruction is a crucial step in seismic data preprocessing.To improve reconstruction speed and save memory,the commonly used three-dimensional(3D)seismic data reconstruction method divides the missing data into a series of time slices and independently reconstructs each time slice.However,when this strategy is employed,the potential correlations between two adjacent time slices are ignored,which degrades reconstruction performance.Therefore,this study proposes the use of a two-dimensional curvelet transform and the fast iterative shrinkage thresholding algorithm for data reconstruction.Based on the significant overlapping characteristics between the curvelet coefficient support sets of two adjacent time slices,a weighted operator is constructed in the curvelet domain using the prior support set provided by the previous reconstructed time slice to delineate the main energy distribution range,eff ectively providing prior information for reconstructing adjacent slices.Consequently,the resulting weighted fast iterative shrinkage thresholding algorithm can be used to reconstruct 3D seismic data.The processing of synthetic and field data shows that the proposed method has higher reconstruction accuracy and faster computational speed than the conventional fast iterative shrinkage thresholding algorithm for handling missing 3D seismic data.展开更多
Objective To investigate the image quality, radiation dose and diagnostic value of the low-tube-voltage high-pitch dual-source computed tomography(DSCT) with sinogram affirmed iterative reconstruction(SAFIRE) for non-...Objective To investigate the image quality, radiation dose and diagnostic value of the low-tube-voltage high-pitch dual-source computed tomography(DSCT) with sinogram affirmed iterative reconstruction(SAFIRE) for non-enhanced abdominal and pelvic scans. Methods This institutional review board-approved prospective study included 64 patients who gave written informed consent for additional abdominal and pelvic scan with DSCT in the period from November to December 2012. The patients underwent standard non-enhanced CT scans(protocol 1) [tube voltage of 120 k Vp/pitch of 0.9/filtered back-projection(FBP) reconstruction] followed by high-pitch non-enhanced CT scans(protocol 2)(100 k Vp/3.0/SAFIRE). The total scan time, mean CT number, signal-to-noise ratio(SNR), image quality, lesion detectability and radiation dose were compared between the two protocols. Results The total scan time of protocol 2 was significantly shorter than that of protocol 1(1.4±0.1 seconds vs. 7.6±0.6 seconds, P<0.001). There was no significant difference between protocol 1 and protocol 2 in mean CT number of all organs(liver, 55.4±6.3 HU vs. 56.1±6.8 HU, P=0.214; pancreas, 43.6±5.9 HU vs. 43.7±5.8 HU, P=0.785; spleen, 47.9±3.9 HU vs. 49.4±4.3 HU, P=0.128; kidney, 32.2±2.3 HU vs. 33.1±2.3 HU, P=0.367; abdominal aorta, 44.8±5.6 HU vs. 45.0±5.5 HU, P=0.499; psoas muscle, 50.7±4.1 HU vs. 50.3±4.5 HU, P=0.279). SNR on images of protocol 2 was higher than that of protocol 1(liver, 5.0±1.2 vs. 4.5±1.1, P<0.001; pancreas, 4.0±1.0 vs. 3.6±0.8, P<0.001; spleen, 4.7±1.0 vs. 4.1±0.9, P<0.001; kidney, 3.1±0.6 vs. 2.8±0.6, P<0.001; abdominal aorta, 4.1±1.0 vs. 3.8±1.0, P<0.001; psoas muscle, 4.5±1.1 vs. 4.3±1.2, P=0.012). The overall image noise of protocol 2 was lower than that of protocol1(9.8±3.1 HU vs. 11.1±3.0 HU, P<0.001). Image quality of protocol 2 was good but lower than that of protocol 1(4.1±0.7 vs. 4.6±0.5, P<0.001). Protocol 2 perceived 229 of 234 lesions(97.9%) that were detected in protocol 1 in the abdomen and pelvis. Radiation dose of protocol 2 was lower than that of protocol 1(4.4±0.4 m Sv vs. 7.3±2.4 m Sv, P<0.001) and the mean dose reduction was 41.4%. Conclusion The high-pitch DSCT with SAFIRE can shorten scan time and reduce radiation dose while preserving image quality in non-enhanced abdominal and pelvic scans.展开更多
AIM: To evaluate the image quality of hepatic multidetector computed tomography(MDCT) with dynamic contrast enhancement. METHODS: It uses iodixanol 270 mg/m L(Visipaque 270) and 80 kVp acquisitions reconstructed with ...AIM: To evaluate the image quality of hepatic multidetector computed tomography(MDCT) with dynamic contrast enhancement. METHODS: It uses iodixanol 270 mg/m L(Visipaque 270) and 80 kVp acquisitions reconstructed with sinogram affirmed iterative reconstruction(SAFIRE?) in comparison with a standard MDCT protocol. Fiftythree consecutive patients with known or suspected hepatocellular carcinoma underwent 55 CT examinations, with two different four-phase CT protocols. The first group of 30 patients underwent a standard 120 kVp acquisition after injection of Iohexol 350 mg/m L(Accupaque 350~?) and reconstructed with filtered back projection. The second group of 25 patients underwent a dual-energy CT at 80-140 kVp with iodixanol 270. The 80 kVp component of the second group was reconstructed iteratively(SAFIRE?-Siemens). All hyperdense and hypodense hepatic lesions ≥ 5 mm were identified with both protocols. Aorta and portal vessels/liver parenchyma contrast to noise ratio(CNR) in arterial phase, hypervascular lesion/liver parenchyma CNR in arterial phase, hypodense lesion/liver parenchyma CNR in portal and late phase were calculated in both groups.RESULTS: Aorta/liver and focal lesions altogether/liver CNR were higher for the second protocol(P = 0.0078 and 0.0346). Hypervascular lesions/liver CNR was not statistically different(P = 0.86). Hypodense lesion/liver CNR in the portal phase was significantly higher for the second group(P = 0.0107). Hypodense lesion/liver CNR in the late phase was the same for both groups(P = 0.9926).CONCLUSION: MDCT imaging with 80 kVp with iterative reconstruction and iodixanol 270 yields equal or even better image quality.展开更多
AIM:To explore whether computer tomography coronary angiography(CTCA) using iterative reconstruction(IR) leads to significant radiation dose reduction without a significant loss in image interpretability compared to c...AIM:To explore whether computer tomography coronary angiography(CTCA) using iterative reconstruction(IR) leads to significant radiation dose reduction without a significant loss in image interpretability compared to conventional filtered back projection(FBP).METHODS:A consecutive series of 200 patients referred to our institution to undergo CTCA constituted the study population.Patients were sequentially assigned to FBP or IR.All studies were acquired with a 256-slice CT scanner.A coronary segment was considered interpretable if image quality was adequate for evaluation of coronary lesions in all segments ≥ 1.5 mm.RESULTS:The mean age was 56.3±9.6 years and165(83%) were male,with no significant differences between groups.Most scans were acquired using prospective ECG triggering,without differences between groups(FBP 84%vs IR 82%;P=0.71).A total of 3198(94%) coronary segments were deemed of diagnostic quality.The percent assessable coronary segments was similar between groups(FBP 91.7%±4.0% vs IR92.5% ± 2.8%; P=0.12).Radiation dose was significantly lower in the IR group(2.8±1.4 mSvvs 4.6±3.0mSv;P<0.0001).Image noise(37.8±1.4 HUvs 38.2±2.4 HU; P=0.20) and signal density(461.7±51.9HU vs 462.2±51.2 HU; P=0.54) levels did not differ between FBP and IR groups,respectively.The IR group was associated to significant effective dose reductions,irrespective of the acquisition mode.CONCLUSION:Application of IR in CTCA preserves image interpretability despite a significant reduction in radiation dose.展开更多
Sparse-view x-ray computed tomography （CT） imaging is an interesting topic in CT field and can efficiently decrease radiation dose. Compared with spatial reconstruction, a Fourier-based algorithm has advantages in r...Sparse-view x-ray computed tomography （CT） imaging is an interesting topic in CT field and can efficiently decrease radiation dose. Compared with spatial reconstruction, a Fourier-based algorithm has advantages in reconstruction speed and memory usage. A novel Fourier-based iterative reconstruction technique that utilizes non-uniform fast Fourier transform （NUFFF） is presented in this work along with advanced total variation （TV） regularization for a fan sparse-view CT. The proposition of a selective matrix contributes to improve reconstruction quality. The new method employs the NUFFT and its adjoin to iterate back and forth between the Fourier and image space. The performance of the proposed algorithm is demonstrated through a series of digital simulations and experimental phantom studies. Results of the proposed algorithm are compared with those of existing TV-regularized techniques based on compressed sensing method, as well as basic algebraic reconstruction technique. Compared with the existing TV-regularized techniques, the proposed Fourier-based technique significantly improves convergence rate and reduces memory allocation, respectively.展开更多
Purpose: To evaluate the quality of three-dimensional (3D) CT angiography images of the abdominal viscera with small focal spot, low tube voltage, and iterative model reconstruction technique (IMR). Materials and Meth...Purpose: To evaluate the quality of three-dimensional (3D) CT angiography images of the abdominal viscera with small focal spot, low tube voltage, and iterative model reconstruction technique (IMR). Materials and Methods: Seven patients with suspected disease of the pancreatobiliary system had undergone CT with high-quality CTA protocol in the present study. There were 5 men and 2 women, ranging in age from 52 to 80 years (mean: 64 years). Results: Depiction of abdominal small artery, small portal vein was possible in all cases. In two cases that we were able to compare, it was superior to standard CTA in small vascular depiction in CTA made clearly in high quality protocol. Conclusions: Although the use of small focal spot, low tube voltage, and IMR can produce higher-quality images of abdominal vessels than standard CTA, this improvement is not significant at elevated radiation doses.展开更多
Structural shape monitoring plays a vital role in the structural health monitoring systems.The inverse finite element method(iFEM)has been demonstrated to be a practical method of deformation reconstruction owing to i...Structural shape monitoring plays a vital role in the structural health monitoring systems.The inverse finite element method(iFEM)has been demonstrated to be a practical method of deformation reconstruction owing to its unique advantages.Current iFEM formulations have been applied to small deformation of structures based on the small-displacement assumption of linear theory.However,this assumption may be inapplicable to some structures with large displacements in practical applications.Therefore,geometric nonlinearity needs to be considered.In this study,to expand the practical utility of iFEM for large displacement monitoring,we propose a nonlinear iFEM algorithm based on a four-node inverse quadrilateral shell element iQS4.Taking the advantage of an iterative iFEM algorithm,a nonlinear response is linearized to compute the geometrically nonlinear deformation reconstruction,like the basic concept of nonlinear FE analysis.Several examples are solved to verify the proposed approach.It is demonstrated that large displacements can be accurately estimated even if the in-situ sensor data includes different levels of randomly generated noise.It is proven that the nonlinear iFEM algorithm provides a more accurate displacement response as compared to the linear iFEM methodology for structures undergoing large displacement.Hence,the proposed approach can be utilized as a viable tool to effectively characterize geometrically nonlinear deformations of structures in real-time applications.展开更多
Objective To evaluate the feasibility of using a low concentration of contrast medium （Visipaque 270 mgl/mL）, low tube voltage, and an advanced image reconstruction algorithm in head and neck computed tomography ang...Objective To evaluate the feasibility of using a low concentration of contrast medium （Visipaque 270 mgl/mL）, low tube voltage, and an advanced image reconstruction algorithm in head and neck computed tomography angiography （CTA）. Methods Forty patients （22 men and 18 women; average age 48.7 ± 14.25 years; average body mass index 23.9 ± 3.7 kg/m^2） undergoing CTA for suspected vascular diseases were randomly assigned into two groups. Group A （n = 20） was administered 370 mgl/mL contrast medium, and group B （n = 20） was administered 270 mgl/mL contrast medium. Both groups were administered at a rate of 4.8 mL/s and an injection volume of 0.8 mL/kg. Images of group A were obtained with 120 kVp and filtered back projection （FBP） reconstruction, whereas images of group B were obtained with 80 kVp and 80% adaptive iterative statistical reconstruction algorithm （ASiR）. The CT values and standard deviations of intracranial arteries and image noise on the corona radiata were measured to calculate the contrast-to-noise ratio （CNR） and signal-to-noise ratio （SNR）. The beam-hardening artifacts （BHAs） around the skull base were calculated. Two readers evaluated the image quality with volume rendered images using scores from 1 to 5. The values between the two groups were statistically compared. Results The mean CT value of the intracranial arteries in group B was significantly higher than that in group A （P 〈 0.001）. The CNR and SNR values in group B were also statistically higher than those in group A （P 〈 0.001）. Image noise and BHAs were not significantly different between the two groups. The image quality score of VR images of in group B was significantly higher than that in group A （P = 0.001）. However, the quality scores of axial enhancement images in group B became significantly smaller than those in group A （P〈 0.001）. The CT dose index volume and dose-length product were decreased by 63.8% and 64%, respectively, in group B （P 〈 0.001 for both）. Conclusion Visipaque combined with 80 kVp and 80% ASiR provided similar image quality in intracranial CTA with 64% radiation dose reduction compared with the use of lopamidol, 120 kVp, and FBP reconstruc-tion.展开更多
The projection matrix model is used to describe the physical relationship between reconstructed object and projection.Such a model has a strong influence on projection and backprojection,two vital operations in iterat...The projection matrix model is used to describe the physical relationship between reconstructed object and projection.Such a model has a strong influence on projection and backprojection,two vital operations in iterative computed tomographic reconstruction.The distance-driven model（DDM） is a state-of-the-art technology that simulates forward and back projections.This model has a low computational complexity and a relatively high spatial resolution;however,it includes only a few methods in a parallel operation with a matched model scheme.This study introduces a fast and parallelizable algorithm to improve the traditional DDM for computing the parallel projection and backprojection operations.Our proposed model has been implemented on a GPU（graphic processing unit） platform and has achieved satisfactory computational efficiency with no approximation.The runtime for the projection and backprojection operations with our model is approximately 4.5 s and 10.5 s per loop,respectively,with an image size of 256×256×256 and 360 projections with a size of 512×512.We compare several general algorithms that have been proposed for maximizing GPU efficiency by using the unmatched projection/backprojection models in a parallel computation.The imaging resolution is not sacrificed and remains accurate during computed tomographic reconstruction.展开更多
Electrical impedance tomography (EIT) is a technique for determining the electrical conductivity and permittivity distribution inside a medium from measurements made on its surface. The impedance distribution reconstr...Electrical impedance tomography (EIT) is a technique for determining the electrical conductivity and permittivity distribution inside a medium from measurements made on its surface. The impedance distribution reconstruction in EIT is a nonlinear inverse problem that requires the use of a regularization method. The generalized Tikhonov regularization methods are often used in solving inverse problems. However, for EIT image reconstruction, the generalized Tikhonov regularization methods may lose the boundary information due to its smoothing operation. In this paper, we propose an iterative Lavrentiev regularization and L-curve-based algorithm to reconstruct EIT images. The regularization parameter should be carefully chosen, but it is often heuristically selected in the conventional regularization-based reconstruction algorithms. So, an L-curve-based optimization algorithm is used for selecting the Lavrentiev regularization parameter. Numerical analysis and simulation results are performed to illustrate EIT image reconstruction. It is shown that choosing the appropriate regularization parameter plays an important role in reconstructing EIT images.展开更多
Objective: to explore the efficacy of low-dose chest CT scanning combined with KARL3D iterative reconstruction technique. Methods: 100 patients who underwent chest CT examination in our hospital were randomly selected...Objective: to explore the efficacy of low-dose chest CT scanning combined with KARL3D iterative reconstruction technique. Methods: 100 patients who underwent chest CT examination in our hospital were randomly selected as the analysis objects and randomly divided into four groups, each with 25 cases. Group A was reconstructed by FBP algorithm and the pipeline current was 150mA;. The low-dose groups B, C and D were reconstructed by Karl algorithm with tube currents of 80mas, 60mas and 40mas respectively. The radiation dose and subjective and objective scores of the four groups were compared. SPSS210 statistical software was used for data analysis. Results: the radiation dose of group B, C and D using low dose tube current combined with karl3d iterative reconstruction technique was lower than that of group a (p < 0.05);The difference between the objective assessment (SD, SNR) and subjective score of the four groups was all p>0.05. Conclusion: The combination of low-dose iterative reconstruction technique and Karl 3D chest CT scan can obtain good images and reduce radiation dose, which is worthy of promotion in the industry.展开更多
Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique （SIRT） inversion ...Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique （SIRT） inversion is used to deduce the stress redistribution around the longwall mining panel. The mining-induced microseismic events were recorded by mounting an array of receivers on the surface, above the active panel. After processing and filtering the seismic data, the three-dimensional tomography images of the p-wave velocity variations by SIRT passive seismic velocity tomography were provided. To display the velocity changes on coal seam level and subsequently to infer the stress redistribution, these three-dimensional tomograms into the coal seam level were sliced. In addition, the boundary element method （BEM） was used to simulate the stress redistribution. The results show that the inferred stresses from the passive seismic tomograms are conformed to numerical models and theoretical concept of the stress redistribution around the longwall panel. In velocity tomograms, the main zones of the stress redistribution arotmd the panel, including front and side abutment pressures, and gob stress are obvious and also the movement of stress zones along the face advancement is evident. Moreover, the effect of the advance rate of the face on the stress redistribution is demonstrated in tomography images. The research result proves that the SIRT passive seismic velocity tomography has an ultimate potential for monitoring the changes of stress redistribution around the longwall mining panel continuously and subsequently to improve safety of mining operations.展开更多
The simultaneous iterations rithms of the ART family. It is used reconstruction technique （SIRT） widely in tomography because of is one of several reconstruction algoits convenience in dealing with large sparse matr...The simultaneous iterations rithms of the ART family. It is used reconstruction technique （SIRT） widely in tomography because of is one of several reconstruction algoits convenience in dealing with large sparse matrices. Its theoretical background and iteration model are discussed at the beginning of this paper. Then, the implementation of the SIRT to reconstruct the three-dimensional distribution of water vapor by simulation is discussed. The results show that the SIRT can function effectively in water vapor tomography, obtain rapid convergence, and be implemented more easily than inversion.展开更多
In this advanced exploration, we focus on multiple parameters estimation in bistatic Multiple-Input Multiple-Output(MIMO) radar systems, a crucial technique for target localization and imaging. Our research innovative...In this advanced exploration, we focus on multiple parameters estimation in bistatic Multiple-Input Multiple-Output(MIMO) radar systems, a crucial technique for target localization and imaging. Our research innovatively addresses the joint estimation of the Direction of Departure(DOD), Direction of Arrival(DOA), and Doppler frequency for incoherent targets. We propose a novel approach that significantly reduces computational complexity by utilizing the TemporalSpatial Nested Sampling Model(TSNSM). Our methodology begins with a multi-linear mapping mechanism to efficiently eliminate unnecessary virtual Degrees of Freedom(DOFs) and reorganize the remaining ones. We then employ the Toeplitz matrix triple iteration reconstruction method, surpassing the traditional Temporal-Spatial Smoothing Window(TSSW) approach, to mitigate the single snapshot effect and reduce computational demands. We further refine the highdimensional ESPRIT algorithm for joint estimation of DOD, DOA, and Doppler frequency, eliminating the need for additional parameter pairing. Moreover, we meticulously derive the Cramér-Rao Bound(CRB) for the TSNSM. This signal model allows for a second expansion of DOFs in time and space domains, achieving high precision in target angle and Doppler frequency estimation with low computational complexity. Our adaptable algorithm is validated through simulations and is suitable for sparse array MIMO radars with various structures, ensuring higher precision in parameter estimation with less complexity burden.展开更多
The image reconstruction of electrical impedance tomography(EIT)is a nonlinear and ill-posed inverse problem and the imaging results are easily affected by measurement noise,which needs to be solved by using regulariz...The image reconstruction of electrical impedance tomography(EIT)is a nonlinear and ill-posed inverse problem and the imaging results are easily affected by measurement noise,which needs to be solved by using regularization methods.The iterative regularization method has become a focus of the research due to its ease of implementation.To deal with the ill-posed and ill-conditional problems in image reconstruction,the inexact Newton-Landweber iterative method is considered and the Nesterov’s acceleration strategy is introduced.One Nesterov-type accelerated version of the inexact Newton-Landweber iteration is presented to determine the conductivity distributions inside an object from electrical measurements made on the surface.In order to further optimize the acceleration,both the steepest descent step-length and the minimal error step-length are exploited during the iterative image reconstruction process.Landweber iteration and its accelerated version are also implemented for comparison.All algorithms are terminated by the discrepancy principle.Finally,the performance is tested by reporting numerical simulations to verify the remarkable acceleration efficiency of the proposed method.展开更多
Missing data are a problem in geophysical surveys, and interpolation and reconstruction of missing data is part of the data processing and interpretation. Based on the sparseness of the geophysical data or the transfo...Missing data are a problem in geophysical surveys, and interpolation and reconstruction of missing data is part of the data processing and interpretation. Based on the sparseness of the geophysical data or the transform domain, we can improve the accuracy and stability of the reconstruction by transforming it to a sparse optimization problem. In this paper, we propose a mathematical model for the sparse reconstruction of data based on the LO-norm minimization. Furthermore, we discuss two types of the approximation algorithm for the LO- norm minimization according to the size and characteristics of the geophysical data： namely, the iteratively reweighted least-squares algorithm and the fast iterative hard thresholding algorithm. Theoretical and numerical analysis showed that applying the iteratively reweighted least-squares algorithm to the reconstruction of potential field data exploits its fast convergence rate, short calculation time, and high precision, whereas the fast iterative hard thresholding algorithm is more suitable for processing seismic data, moreover, its computational efficiency is better than that of the traditional iterative hard thresholding algorithm.展开更多
In the medical computer tomography （CT） field, total variation （TV）, which is the l1-norm of the discrete gradient transform （DGT）, is widely used as regularization based on the compressive sensing （CS） theory...In the medical computer tomography （CT） field, total variation （TV）, which is the l1-norm of the discrete gradient transform （DGT）, is widely used as regularization based on the compressive sensing （CS） theory. To overcome the TV model＇s disadvantageous tendency of uniformly penalizing the image gradient and over smoothing the low-contrast structures, an iterative algorithm based on the l0-norm optimization of the DGT is proposed. In order to rise to the challenges introduced by the l0-norm DGT, the algorithm uses a pseudo-inverse transform of DGT and adapts an iterative hard thresholding （IHT） algorithm, whose convergence and effective efficiency have been theoretically proven. The simulation demonstrates our conclusions and indicates that the algorithm proposed in this paper can obviously improve the reconstruction quality.展开更多
The order of the projection in the algebraic reconstruction technique(ART)method has great influence on the rate of the convergence.Although many scholars have studied the order of the projection,few theoretical proof...The order of the projection in the algebraic reconstruction technique(ART)method has great influence on the rate of the convergence.Although many scholars have studied the order of the projection,few theoretical proofs are given.Thomas Strohmer and Roman Vershynin introduced a randomized version of the Kaczmarz method for consistent,and over-determined linear systems and proved whose rate does not depend on the number of equations in the systems in 2009.In this paper,we apply this method to computed tomography(CT)image reconstruction and compared images generated by the sequential Kaczmarz method and the randomized Kaczmarz method.Experiments demonstrates the feasibility of the randomized Kaczmarz algorithm in CT image reconstruction and its exponential curve convergence.展开更多
基金supported by the Scientific and Technological Developing Scheme of Jilin Province,China(No.20240101371JC)the National Natural Science Foundation of China(No.62107008).
文摘A Bayesian network reconstruction method based on norm minimization is proposed to address the sparsity and iterative divergence issues in network reconstruction caused by noise and missing values.This method achieves precise adjustment of the network structure by constructing a preliminary random network model and introducing small-world network characteristics and combines L1 norm minimization regularization techniques to control model complexity and optimize the inference process of variable dependencies.In the experiment of game network reconstruction,when the success rate of the L1 norm minimization model’s existence connection reconstruction reaches 100%,the minimum data required is about 40%,while the minimum data required for a sparse Bayesian learning network is about 45%.In terms of operational efficiency,the running time for minimizing the L1 normis basically maintained at 1.0 s,while the success rate of connection reconstruction increases significantly with an increase in data volume,reaching a maximum of 13.2 s.Meanwhile,in the case of a signal-to-noise ratio of 10 dB,the L1 model achieves a 100% success rate in the reconstruction of existing connections,while the sparse Bayesian network had the highest success rate of 90% in the reconstruction of non-existent connections.In the analysis of actual cases,the maximum lift and drop track of the research method is 0.08 m.The mean square error is 5.74 cm^(2).The results indicate that this norm minimization-based method has good performance in data efficiency and model stability,effectively reducing the impact of outliers on the reconstruction results to more accurately reflect the actual situation.
基金National Natural Science Foundation of China under Grant 42304145Jiangxi Provincial Natural Science Foundation under Grant 20242BAB26051,20242BAB25191 and 20232BAB213077+1 种基金Foundation of National Key Laboratory of Uranium Resources Exploration-Mining and Nuclear Remote Sensing under Grant 2024QZ-TD-13Open Fund(FW0399-0002)of SINOPEC Key Laboratory of Geophysics。
文摘Data reconstruction is a crucial step in seismic data preprocessing.To improve reconstruction speed and save memory,the commonly used three-dimensional(3D)seismic data reconstruction method divides the missing data into a series of time slices and independently reconstructs each time slice.However,when this strategy is employed,the potential correlations between two adjacent time slices are ignored,which degrades reconstruction performance.Therefore,this study proposes the use of a two-dimensional curvelet transform and the fast iterative shrinkage thresholding algorithm for data reconstruction.Based on the significant overlapping characteristics between the curvelet coefficient support sets of two adjacent time slices,a weighted operator is constructed in the curvelet domain using the prior support set provided by the previous reconstructed time slice to delineate the main energy distribution range,eff ectively providing prior information for reconstructing adjacent slices.Consequently,the resulting weighted fast iterative shrinkage thresholding algorithm can be used to reconstruct 3D seismic data.The processing of synthetic and field data shows that the proposed method has higher reconstruction accuracy and faster computational speed than the conventional fast iterative shrinkage thresholding algorithm for handling missing 3D seismic data.
文摘Objective To investigate the image quality, radiation dose and diagnostic value of the low-tube-voltage high-pitch dual-source computed tomography(DSCT) with sinogram affirmed iterative reconstruction(SAFIRE) for non-enhanced abdominal and pelvic scans. Methods This institutional review board-approved prospective study included 64 patients who gave written informed consent for additional abdominal and pelvic scan with DSCT in the period from November to December 2012. The patients underwent standard non-enhanced CT scans(protocol 1) [tube voltage of 120 k Vp/pitch of 0.9/filtered back-projection(FBP) reconstruction] followed by high-pitch non-enhanced CT scans(protocol 2)(100 k Vp/3.0/SAFIRE). The total scan time, mean CT number, signal-to-noise ratio(SNR), image quality, lesion detectability and radiation dose were compared between the two protocols. Results The total scan time of protocol 2 was significantly shorter than that of protocol 1(1.4±0.1 seconds vs. 7.6±0.6 seconds, P<0.001). There was no significant difference between protocol 1 and protocol 2 in mean CT number of all organs(liver, 55.4±6.3 HU vs. 56.1±6.8 HU, P=0.214; pancreas, 43.6±5.9 HU vs. 43.7±5.8 HU, P=0.785; spleen, 47.9±3.9 HU vs. 49.4±4.3 HU, P=0.128; kidney, 32.2±2.3 HU vs. 33.1±2.3 HU, P=0.367; abdominal aorta, 44.8±5.6 HU vs. 45.0±5.5 HU, P=0.499; psoas muscle, 50.7±4.1 HU vs. 50.3±4.5 HU, P=0.279). SNR on images of protocol 2 was higher than that of protocol 1(liver, 5.0±1.2 vs. 4.5±1.1, P<0.001; pancreas, 4.0±1.0 vs. 3.6±0.8, P<0.001; spleen, 4.7±1.0 vs. 4.1±0.9, P<0.001; kidney, 3.1±0.6 vs. 2.8±0.6, P<0.001; abdominal aorta, 4.1±1.0 vs. 3.8±1.0, P<0.001; psoas muscle, 4.5±1.1 vs. 4.3±1.2, P=0.012). The overall image noise of protocol 2 was lower than that of protocol1(9.8±3.1 HU vs. 11.1±3.0 HU, P<0.001). Image quality of protocol 2 was good but lower than that of protocol 1(4.1±0.7 vs. 4.6±0.5, P<0.001). Protocol 2 perceived 229 of 234 lesions(97.9%) that were detected in protocol 1 in the abdomen and pelvis. Radiation dose of protocol 2 was lower than that of protocol 1(4.4±0.4 m Sv vs. 7.3±2.4 m Sv, P<0.001) and the mean dose reduction was 41.4%. Conclusion The high-pitch DSCT with SAFIRE can shorten scan time and reduce radiation dose while preserving image quality in non-enhanced abdominal and pelvic scans.
文摘AIM: To evaluate the image quality of hepatic multidetector computed tomography(MDCT) with dynamic contrast enhancement. METHODS: It uses iodixanol 270 mg/m L(Visipaque 270) and 80 kVp acquisitions reconstructed with sinogram affirmed iterative reconstruction(SAFIRE?) in comparison with a standard MDCT protocol. Fiftythree consecutive patients with known or suspected hepatocellular carcinoma underwent 55 CT examinations, with two different four-phase CT protocols. The first group of 30 patients underwent a standard 120 kVp acquisition after injection of Iohexol 350 mg/m L(Accupaque 350~?) and reconstructed with filtered back projection. The second group of 25 patients underwent a dual-energy CT at 80-140 kVp with iodixanol 270. The 80 kVp component of the second group was reconstructed iteratively(SAFIRE?-Siemens). All hyperdense and hypodense hepatic lesions ≥ 5 mm were identified with both protocols. Aorta and portal vessels/liver parenchyma contrast to noise ratio(CNR) in arterial phase, hypervascular lesion/liver parenchyma CNR in arterial phase, hypodense lesion/liver parenchyma CNR in portal and late phase were calculated in both groups.RESULTS: Aorta/liver and focal lesions altogether/liver CNR were higher for the second protocol(P = 0.0078 and 0.0346). Hypervascular lesions/liver CNR was not statistically different(P = 0.86). Hypodense lesion/liver CNR in the portal phase was significantly higher for the second group(P = 0.0107). Hypodense lesion/liver CNR in the late phase was the same for both groups(P = 0.9926).CONCLUSION: MDCT imaging with 80 kVp with iterative reconstruction and iodixanol 270 yields equal or even better image quality.
文摘AIM:To explore whether computer tomography coronary angiography(CTCA) using iterative reconstruction(IR) leads to significant radiation dose reduction without a significant loss in image interpretability compared to conventional filtered back projection(FBP).METHODS:A consecutive series of 200 patients referred to our institution to undergo CTCA constituted the study population.Patients were sequentially assigned to FBP or IR.All studies were acquired with a 256-slice CT scanner.A coronary segment was considered interpretable if image quality was adequate for evaluation of coronary lesions in all segments ≥ 1.5 mm.RESULTS:The mean age was 56.3±9.6 years and165(83%) were male,with no significant differences between groups.Most scans were acquired using prospective ECG triggering,without differences between groups(FBP 84%vs IR 82%;P=0.71).A total of 3198(94%) coronary segments were deemed of diagnostic quality.The percent assessable coronary segments was similar between groups(FBP 91.7%±4.0% vs IR92.5% ± 2.8%; P=0.12).Radiation dose was significantly lower in the IR group(2.8±1.4 mSvvs 4.6±3.0mSv;P<0.0001).Image noise(37.8±1.4 HUvs 38.2±2.4 HU; P=0.20) and signal density(461.7±51.9HU vs 462.2±51.2 HU; P=0.54) levels did not differ between FBP and IR groups,respectively.The IR group was associated to significant effective dose reductions,irrespective of the acquisition mode.CONCLUSION:Application of IR in CTCA preserves image interpretability despite a significant reduction in radiation dose.
基金Projected supported by the National High Technology Research and Development Program of China(Grant No.2012AA011603)the National Natura Science Foundation of China(Grant No.61372172)
文摘Sparse-view x-ray computed tomography （CT） imaging is an interesting topic in CT field and can efficiently decrease radiation dose. Compared with spatial reconstruction, a Fourier-based algorithm has advantages in reconstruction speed and memory usage. A novel Fourier-based iterative reconstruction technique that utilizes non-uniform fast Fourier transform （NUFFF） is presented in this work along with advanced total variation （TV） regularization for a fan sparse-view CT. The proposition of a selective matrix contributes to improve reconstruction quality. The new method employs the NUFFT and its adjoin to iterate back and forth between the Fourier and image space. The performance of the proposed algorithm is demonstrated through a series of digital simulations and experimental phantom studies. Results of the proposed algorithm are compared with those of existing TV-regularized techniques based on compressed sensing method, as well as basic algebraic reconstruction technique. Compared with the existing TV-regularized techniques, the proposed Fourier-based technique significantly improves convergence rate and reduces memory allocation, respectively.
文摘Purpose: To evaluate the quality of three-dimensional (3D) CT angiography images of the abdominal viscera with small focal spot, low tube voltage, and iterative model reconstruction technique (IMR). Materials and Methods: Seven patients with suspected disease of the pancreatobiliary system had undergone CT with high-quality CTA protocol in the present study. There were 5 men and 2 women, ranging in age from 52 to 80 years (mean: 64 years). Results: Depiction of abdominal small artery, small portal vein was possible in all cases. In two cases that we were able to compare, it was superior to standard CTA in small vascular depiction in CTA made clearly in high quality protocol. Conclusions: Although the use of small focal spot, low tube voltage, and IMR can produce higher-quality images of abdominal vessels than standard CTA, this improvement is not significant at elevated radiation doses.
基金supported by the NationalNatural Science Foundation of China(Grant No.11902253)the Fundamental Research Funds for the Central Universities of China.The authors are grateful for this support.
文摘Structural shape monitoring plays a vital role in the structural health monitoring systems.The inverse finite element method(iFEM)has been demonstrated to be a practical method of deformation reconstruction owing to its unique advantages.Current iFEM formulations have been applied to small deformation of structures based on the small-displacement assumption of linear theory.However,this assumption may be inapplicable to some structures with large displacements in practical applications.Therefore,geometric nonlinearity needs to be considered.In this study,to expand the practical utility of iFEM for large displacement monitoring,we propose a nonlinear iFEM algorithm based on a four-node inverse quadrilateral shell element iQS4.Taking the advantage of an iterative iFEM algorithm,a nonlinear response is linearized to compute the geometrically nonlinear deformation reconstruction,like the basic concept of nonlinear FE analysis.Several examples are solved to verify the proposed approach.It is demonstrated that large displacements can be accurately estimated even if the in-situ sensor data includes different levels of randomly generated noise.It is proven that the nonlinear iFEM algorithm provides a more accurate displacement response as compared to the linear iFEM methodology for structures undergoing large displacement.Hence,the proposed approach can be utilized as a viable tool to effectively characterize geometrically nonlinear deformations of structures in real-time applications.
文摘Objective To evaluate the feasibility of using a low concentration of contrast medium （Visipaque 270 mgl/mL）, low tube voltage, and an advanced image reconstruction algorithm in head and neck computed tomography angiography （CTA）. Methods Forty patients （22 men and 18 women; average age 48.7 ± 14.25 years; average body mass index 23.9 ± 3.7 kg/m^2） undergoing CTA for suspected vascular diseases were randomly assigned into two groups. Group A （n = 20） was administered 370 mgl/mL contrast medium, and group B （n = 20） was administered 270 mgl/mL contrast medium. Both groups were administered at a rate of 4.8 mL/s and an injection volume of 0.8 mL/kg. Images of group A were obtained with 120 kVp and filtered back projection （FBP） reconstruction, whereas images of group B were obtained with 80 kVp and 80% adaptive iterative statistical reconstruction algorithm （ASiR）. The CT values and standard deviations of intracranial arteries and image noise on the corona radiata were measured to calculate the contrast-to-noise ratio （CNR） and signal-to-noise ratio （SNR）. The beam-hardening artifacts （BHAs） around the skull base were calculated. Two readers evaluated the image quality with volume rendered images using scores from 1 to 5. The values between the two groups were statistically compared. Results The mean CT value of the intracranial arteries in group B was significantly higher than that in group A （P 〈 0.001）. The CNR and SNR values in group B were also statistically higher than those in group A （P 〈 0.001）. Image noise and BHAs were not significantly different between the two groups. The image quality score of VR images of in group B was significantly higher than that in group A （P = 0.001）. However, the quality scores of axial enhancement images in group B became significantly smaller than those in group A （P〈 0.001）. The CT dose index volume and dose-length product were decreased by 63.8% and 64%, respectively, in group B （P 〈 0.001 for both）. Conclusion Visipaque combined with 80 kVp and 80% ASiR provided similar image quality in intracranial CTA with 64% radiation dose reduction compared with the use of lopamidol, 120 kVp, and FBP reconstruc-tion.
基金supported by the National High Technology Research and Development Program of China(Grant No.2012AA011603)the National Natural Science Foundation of China(Grant No.61372172)
文摘The projection matrix model is used to describe the physical relationship between reconstructed object and projection.Such a model has a strong influence on projection and backprojection,two vital operations in iterative computed tomographic reconstruction.The distance-driven model（DDM） is a state-of-the-art technology that simulates forward and back projections.This model has a low computational complexity and a relatively high spatial resolution;however,it includes only a few methods in a parallel operation with a matched model scheme.This study introduces a fast and parallelizable algorithm to improve the traditional DDM for computing the parallel projection and backprojection operations.Our proposed model has been implemented on a GPU（graphic processing unit） platform and has achieved satisfactory computational efficiency with no approximation.The runtime for the projection and backprojection operations with our model is approximately 4.5 s and 10.5 s per loop,respectively,with an image size of 256×256×256 and 360 projections with a size of 512×512.We compare several general algorithms that have been proposed for maximizing GPU efficiency by using the unmatched projection/backprojection models in a parallel computation.The imaging resolution is not sacrificed and remains accurate during computed tomographic reconstruction.
文摘Electrical impedance tomography (EIT) is a technique for determining the electrical conductivity and permittivity distribution inside a medium from measurements made on its surface. The impedance distribution reconstruction in EIT is a nonlinear inverse problem that requires the use of a regularization method. The generalized Tikhonov regularization methods are often used in solving inverse problems. However, for EIT image reconstruction, the generalized Tikhonov regularization methods may lose the boundary information due to its smoothing operation. In this paper, we propose an iterative Lavrentiev regularization and L-curve-based algorithm to reconstruct EIT images. The regularization parameter should be carefully chosen, but it is often heuristically selected in the conventional regularization-based reconstruction algorithms. So, an L-curve-based optimization algorithm is used for selecting the Lavrentiev regularization parameter. Numerical analysis and simulation results are performed to illustrate EIT image reconstruction. It is shown that choosing the appropriate regularization parameter plays an important role in reconstructing EIT images.
文摘Objective: to explore the efficacy of low-dose chest CT scanning combined with KARL3D iterative reconstruction technique. Methods: 100 patients who underwent chest CT examination in our hospital were randomly selected as the analysis objects and randomly divided into four groups, each with 25 cases. Group A was reconstructed by FBP algorithm and the pipeline current was 150mA;. The low-dose groups B, C and D were reconstructed by Karl algorithm with tube currents of 80mas, 60mas and 40mas respectively. The radiation dose and subjective and objective scores of the four groups were compared. SPSS210 statistical software was used for data analysis. Results: the radiation dose of group B, C and D using low dose tube current combined with karl3d iterative reconstruction technique was lower than that of group a (p < 0.05);The difference between the objective assessment (SD, SNR) and subjective score of the four groups was all p>0.05. Conclusion: The combination of low-dose iterative reconstruction technique and Karl 3D chest CT scan can obtain good images and reduce radiation dose, which is worthy of promotion in the industry.
文摘Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique （SIRT） inversion is used to deduce the stress redistribution around the longwall mining panel. The mining-induced microseismic events were recorded by mounting an array of receivers on the surface, above the active panel. After processing and filtering the seismic data, the three-dimensional tomography images of the p-wave velocity variations by SIRT passive seismic velocity tomography were provided. To display the velocity changes on coal seam level and subsequently to infer the stress redistribution, these three-dimensional tomograms into the coal seam level were sliced. In addition, the boundary element method （BEM） was used to simulate the stress redistribution. The results show that the inferred stresses from the passive seismic tomograms are conformed to numerical models and theoretical concept of the stress redistribution around the longwall panel. In velocity tomograms, the main zones of the stress redistribution arotmd the panel, including front and side abutment pressures, and gob stress are obvious and also the movement of stress zones along the face advancement is evident. Moreover, the effect of the advance rate of the face on the stress redistribution is demonstrated in tomography images. The research result proves that the SIRT passive seismic velocity tomography has an ultimate potential for monitoring the changes of stress redistribution around the longwall mining panel continuously and subsequently to improve safety of mining operations.
基金supported by the National Natural Science Foundation of China(40974018)Nationa l863 Plan Projects(2009AA12Z307)
文摘The simultaneous iterations rithms of the ART family. It is used reconstruction technique （SIRT） widely in tomography because of is one of several reconstruction algoits convenience in dealing with large sparse matrices. Its theoretical background and iteration model are discussed at the beginning of this paper. Then, the implementation of the SIRT to reconstruct the three-dimensional distribution of water vapor by simulation is discussed. The results show that the SIRT can function effectively in water vapor tomography, obtain rapid convergence, and be implemented more easily than inversion.
基金supported in part by the National Natural Science Foundation of China(No.62071476)in part by China Postdoctoral Science Foundation(No.2022M723879)in part by the Science and Technology Innovation Program of Hunan Province,China(No.2021RC3080)。
文摘In this advanced exploration, we focus on multiple parameters estimation in bistatic Multiple-Input Multiple-Output(MIMO) radar systems, a crucial technique for target localization and imaging. Our research innovatively addresses the joint estimation of the Direction of Departure(DOD), Direction of Arrival(DOA), and Doppler frequency for incoherent targets. We propose a novel approach that significantly reduces computational complexity by utilizing the TemporalSpatial Nested Sampling Model(TSNSM). Our methodology begins with a multi-linear mapping mechanism to efficiently eliminate unnecessary virtual Degrees of Freedom(DOFs) and reorganize the remaining ones. We then employ the Toeplitz matrix triple iteration reconstruction method, surpassing the traditional Temporal-Spatial Smoothing Window(TSSW) approach, to mitigate the single snapshot effect and reduce computational demands. We further refine the highdimensional ESPRIT algorithm for joint estimation of DOD, DOA, and Doppler frequency, eliminating the need for additional parameter pairing. Moreover, we meticulously derive the Cramér-Rao Bound(CRB) for the TSNSM. This signal model allows for a second expansion of DOFs in time and space domains, achieving high precision in target angle and Doppler frequency estimation with low computational complexity. Our adaptable algorithm is validated through simulations and is suitable for sparse array MIMO radars with various structures, ensuring higher precision in parameter estimation with less complexity burden.
基金National Natural Science Foundation of China(12101204,12261021)Heilongjiang Provincial Natural Science Foundation of China(LH2023A018)Modern Numerical Method Course for Research Program on Teaching Reform of Degree and Postgraduate Education of Heilongjiang University(2024)。
文摘The image reconstruction of electrical impedance tomography(EIT)is a nonlinear and ill-posed inverse problem and the imaging results are easily affected by measurement noise,which needs to be solved by using regularization methods.The iterative regularization method has become a focus of the research due to its ease of implementation.To deal with the ill-posed and ill-conditional problems in image reconstruction,the inexact Newton-Landweber iterative method is considered and the Nesterov’s acceleration strategy is introduced.One Nesterov-type accelerated version of the inexact Newton-Landweber iteration is presented to determine the conductivity distributions inside an object from electrical measurements made on the surface.In order to further optimize the acceleration,both the steepest descent step-length and the minimal error step-length are exploited during the iterative image reconstruction process.Landweber iteration and its accelerated version are also implemented for comparison.All algorithms are terminated by the discrepancy principle.Finally,the performance is tested by reporting numerical simulations to verify the remarkable acceleration efficiency of the proposed method.
基金supported by the National Natural Science Foundation of China (Grant No.41074133)
文摘Missing data are a problem in geophysical surveys, and interpolation and reconstruction of missing data is part of the data processing and interpretation. Based on the sparseness of the geophysical data or the transform domain, we can improve the accuracy and stability of the reconstruction by transforming it to a sparse optimization problem. In this paper, we propose a mathematical model for the sparse reconstruction of data based on the LO-norm minimization. Furthermore, we discuss two types of the approximation algorithm for the LO- norm minimization according to the size and characteristics of the geophysical data： namely, the iteratively reweighted least-squares algorithm and the fast iterative hard thresholding algorithm. Theoretical and numerical analysis showed that applying the iteratively reweighted least-squares algorithm to the reconstruction of potential field data exploits its fast convergence rate, short calculation time, and high precision, whereas the fast iterative hard thresholding algorithm is more suitable for processing seismic data, moreover, its computational efficiency is better than that of the traditional iterative hard thresholding algorithm.
文摘In the medical computer tomography （CT） field, total variation （TV）, which is the l1-norm of the discrete gradient transform （DGT）, is widely used as regularization based on the compressive sensing （CS） theory. To overcome the TV model＇s disadvantageous tendency of uniformly penalizing the image gradient and over smoothing the low-contrast structures, an iterative algorithm based on the l0-norm optimization of the DGT is proposed. In order to rise to the challenges introduced by the l0-norm DGT, the algorithm uses a pseudo-inverse transform of DGT and adapts an iterative hard thresholding （IHT） algorithm, whose convergence and effective efficiency have been theoretically proven. The simulation demonstrates our conclusions and indicates that the algorithm proposed in this paper can obviously improve the reconstruction quality.
基金National Natural Science Foundation of China(No.61171179,No.61171178)Natural Science Foundation of Shanxi Province(No.2010011002-1,No.2010011002-2and No.2012021011-2)
文摘The order of the projection in the algebraic reconstruction technique(ART)method has great influence on the rate of the convergence.Although many scholars have studied the order of the projection,few theoretical proofs are given.Thomas Strohmer and Roman Vershynin introduced a randomized version of the Kaczmarz method for consistent,and over-determined linear systems and proved whose rate does not depend on the number of equations in the systems in 2009.In this paper,we apply this method to computed tomography(CT)image reconstruction and compared images generated by the sequential Kaczmarz method and the randomized Kaczmarz method.Experiments demonstrates the feasibility of the randomized Kaczmarz algorithm in CT image reconstruction and its exponential curve convergence.
