共查询到18条相似文献，搜索用时 173 毫秒
1.
在图像重建中,Landweber迭代算法是图像重建算法中的重要方法.本文将针对Landweber分块迭代算法中松弛参数的选取进行研究.在重建过程中采取对投影矩阵按投影角度分块的方法,选取特定的松弛参数.通过数值实验得出结论:对于按角度分块的块迭代算法,松弛参数选取为λ乘以块矩阵与其共轭转置矩阵乘积的最大特征值分之一,当采集完全投影数据,且λ接近(1/6)～(1/7)时效果最好.另外,本文按角度分块的做法和松弛参数的选取方法对于有限角度图像重建问题也是可行的. 相似文献
2.
由投影重建图像的对称网格迭代算法 总被引：1，自引：0，他引：1
本文对于工业CT检测中常用的代数迭代重建算法提出了改进,利用投影射线之间存在的几何对称结构,提出了图像重建的对称网格迭代算法(简写为SM-IRT).该算法简化了投影系数矩阵的计算,调整了迭代算法逐线校正的迭代顺序.对模拟数据和工业CT实测数据进行了重建图像的数值实验,结果表明:与常规算法比较,本文提出的新算法重建速度快,成像精度高. 相似文献
3.
处理由不完全数据重建图象的方法之一是代数重建法。但是,代数重建法的一个不足是关于重建图象的收敛速度很慢;本文利用代数重建法有限次迭代的重建图象反映原图象的分布趋势和向原图象变化的这一特点,提出一种优化的图象重建外插算法。数值模拟结果表明,新算法在重建图象的收敛速度和质量方面均有改进。 相似文献
4.
在用计算机断层成像方法由EUV观测图像重建等离子体层全球密度分布时,地球的遮挡和有限角度都会导致投影数据不完备,从而无法精确重建出等离子体层的密度分布.本文针对该问题,提出一种基于图像总变差极小化的代数迭代算法.通过重建等离子体层投影数据缺失最为严重的中心子午面,证明该算法能够显著提高重建图像的质量. 并且在IMAGE卫星仅能达到90°的有限投影角度下,此算法重建图像的相关系数可达0.760,而代数迭代算法的相关系数仅为0.696. 相似文献
5.
本文中介绍了一种新型的正电子型核医学仪器——正电子心脑功能仪,并从该仪器的系统要求和几何特点出发,提出了一种特殊的不完备投影的图象重建方法,称为变权反投影迭代重建方法.该方法基于代数重建法的基本理论,利用在反投影过程中的对校正次数不均匀性的补偿,使算法的收敛速度加快,采用该方法进行图象重建,与传统的加法代数重建法相比,收敛速度提高约2倍. 相似文献
6.
序列子集联合代数重建技术 总被引：2，自引：0，他引：2
图像重建迭代算法的主要缺点是计算量大,重建速度慢。为减少计算时间,Hudson等提出了有序子集算法,由于该算法在每次迭代时使用固定的子集个数,重建图像的质量主要依赖于迭代步中的子集数。本文提出序列子集联合代数重建技术,在每次迭代后减少使用的子集个数,这样在加速图像收敛的同时恢复重建图像的各种频率元素。实验结果表明序列子集联合代数重建技术可在少数次迭代后提供较高质量的重建图像,且对噪声数据不敏感。 相似文献
7.
代数重建算法中的投影访问顺序对重建图像的收敛速度和精度有很大的影响。针对三维锥束代数重建算法,比较了几种投影序列排序方案的重建图像收敛速度:顺序访问方案,多水平级排序方案,素数分解排序方案,加权距离排序方案。通过仿真实验得到的实验结果表明:在三维锥束ART算法中按照后三种投影序列排序方案排序后的图像重建比顺序访问方案的收敛速度快,并且WDS排序方案的收敛速度最好。 相似文献
8.
扇束图像重建中OSEM算法及子集划分的研究 总被引：5，自引：2，他引：3
目的：将一种快速迭代图像重建方法OSEM用于仿真Phantom模型的扇束投影数据的重建中,验证其重建参数的有效取值范围,分析比较不同子集选取对扇束OSEM重建的图像质量以及收敛速度的影响。方法：本研究采用事先计算特定几何扇束扫描的概率矩阵和最大间距划分投影子集的方法,选取图像大小为128*128。通过对模拟的Phantom扇束投影数据的重建,结果：OSEM重建效果好于传统滤波反投影(FBP)重建的结果,尤其是在有噪声情况下,在一定范围内重建质量与迭代次数和子集划分个数成正比,而计算时间仅仅与迭代次数正比。结论：恰当地选取子集个数,扇束OSEM重建可用较少迭代次数和较短的计算时间获取好的重建图像。 相似文献
9.
基于最小平方的Fourier地震数据重建方法最终转化为求解一个线性方程组, 其系数矩阵是Toeplitz矩阵,可以用共轭梯度法求解该线性方程组.共轭梯度法的迭代次数受系数矩阵病态程度的影响,地震数据的非规则采样程度越高,所形成的系数矩阵病态程度越高,就越难收敛和得到合理的计算结果.本文研究了基于Toeplitz矩阵的不同预条件的构造方法,以及对共轭梯度法收敛性的影响.通过预条件的使用,加快了共轭梯度法的迭代速度, 改进了共轭梯度算法的收敛性,提高了计算的效率.数值算例和实际地震数据重建试验证明了预条件共轭梯度法对计算效率有很大的提高. 相似文献
10.
一种改进投影系数计算的快速ART算法 总被引：1，自引：0，他引：1
代数重建(ART)是图像重建领域中的重要方法。为了提高ART算法的重建速度,本文提出了一个新的通过判别射线与网格的相交状况快速计算投影系数的算法。该算法主要使用简单的诸如加减法以及比较运算来计算投影系数。避免了传统算法中的求交排序计算,因而重建速度大大提升。仿真实验使用Shepp-Logan模型,实验结果显示该算法重建效率优于传统算法。 相似文献
11.
数字合成X射线体层成像技术能利用有限角度下的投影数据重建物体任意断层的图像.在数宁合成X射线体层成像重建算法研究中,模拟投影数据是重要的步骤,本文提出了一种基于GPU光线投射算法的数字合成X射线体层成像投影数据模拟方法.比较传统CPU模拟手段,GPU模拟方法计算速率快,且基于硬件支持的三线性插值能够得到更加接近实际的投... 相似文献
12.
一种基于数据外插改进的ART迭代算法 总被引：1，自引：0，他引：1
本文针对有限角度的投影数据的CT图像重建问题,提出了一种基于数据外插改进的ART算法.该算法的基本思想足运用已知角度的投影数据来补全未知角度的投影数据,再用ART算法进行图像重建.最后用模拟的投影数据进行了重建图像的数值实验.实验结果表明该算法不但提高了重建图像质量,同时也提高了图像达代的收敛速度. 相似文献
13.
本研究主要目的在于探讨与评估在同步辐射光源显微成像时使用不同的图像重建算法时的图像质量与重建速度.此次实验我们采用鸟类的羽毛做为受测物体,并且使用位于新竹科学工业园区内的同步辐射研究中心的同步辐射光源及高分辨X光摄影机来取得X光的影像.本研究使用两种不同的算法来重建图像,分别为滤波反投影法与梯度最大似然迭代重建法,作为这次评估的算法.由重建图像的结果显示,我们可以发现梯度最大概似迭代式重建法在边缘的地方有模糊化的现象,另外在图像质量方面,滤波反投影算法的影像上的环状伪影较梯度最大似然迭代算法明显.因此本研究就针对目前同步辐射的图像重建法做一探讨,并对两种不同的重建算法做一适当的评估,以期望提供不同的图像重建法给予不同需求的使用者. 相似文献
14.
15.
Hua Zhang Su Diao Wei Chen Guangnan Huang Hongxing Li Min Bai 《Geophysical Prospecting》2019,67(5):1201-1218
Seismic data reconstruction, as a preconditioning process, is critical to the performance of subsequent data and imaging processing tasks. Often, seismic data are sparsely and non-uniformly sampled due to limitations of economic costs and field conditions. However, most reconstruction processing algorithms are designed for the ideal case of uniformly sampled data. In this paper, we propose the non-equispaced fast discrete curvelet transform-based three-dimensional reconstruction method that can handle and interpolate non-uniformly sampled data effectively along two spatial coordinates. In the procedure, the three-dimensional seismic data sets are organized in a sequence of two-dimensional time slices along the source–receiver domain. By introducing the two-dimensional non-equispaced fast Fourier transform in the conventional fast discrete curvelet transform, we formulate an L1 sparsity regularized problem to invert for the uniformly sampled curvelet coefficients from the non-uniformly sampled data. In order to improve the inversion algorithm efficiency, we employ the linearized Bregman method to solve the L1-norm minimization problem. Once the uniform curvelet coefficients are obtained, uniformly sampled three-dimensional seismic data can be reconstructed via the conventional inverse curvelet transform. The reconstructed results using both synthetic and real data demonstrate that the proposed method can reconstruct not only non-uniformly sampled and aliased data with missing traces, but also the subset of observed data on a non-uniform grid to a specified uniform grid along two spatial coordinates. Also, the results show that the simple linearized Bregman method is superior to the complex spectral projected gradient for L1 norm method in terms of reconstruction accuracy. 相似文献
16.
To carry out a 3D prestack migration of the Kirchhoff type is still a task of enormous computational effort. Its efficiency can be significantly enhanced by employing a fast traveltime interpolation algorithm. High accuracy can be achieved if secondorder spatial derivatives of traveltimes are included in order to account for the curvature of the wavefront. We suggest a hyperbolic traveltime interpolation scheme that permits the determination of the hyperbolic coefficients directly from traveltimes sampled on a coarse grid, thus reducing the requirements in data storage. This approach is closely related to the paraxial ray approximation and corresponds to an extension of the wellknown method to arbitrary heterogeneous and complex media in 3D. Application to various velocity models, including a 3D version of the Marmousi model, confirms the superiority of our method over the popular trilinear interpolation. This is especially true for regions with strong curvature of the local wavefront. In contrast to trilinear interpolation, our method also provides the possibility of interpolating source positions, and it is 56 times faster than the calculation of traveltime tables using a fast finitedifference eikonal solver. 相似文献
17.
本文基于锥束CT滤波反投影重建的FDK算法,通过两种算法改进并结合基于共享内存的OpenMP并行技术和代码优化,实现了锥束CT图像的快速重建.基于锥束CT实际投影数据的重建结果表明,图像重建速度得到了较大的提高,断层图像重建质量与FDK原型算法相当. 相似文献
18.
岩样中包体构造的CT实验探测 总被引：1，自引：1，他引：1
CT技术图象重建方法可分为变换法和代数法两大类。本文讨论了几种代数重建法的特点,开展了地球物理CT技术在岩石实验中应用的探讨性工作。研究表明,实验数据的测量误差直接影响着反演结果,当测量误差比较大时,即使数据量很大,也难以提高反演结果的分辨。岩样的吸收、换能器与样品的耦合、岩样均匀性等因素对实验数据的测量精度都有着重要的影响。文中还针对实验条件进行了数值模拟,以检验本方案的可行性,并对当前流行的BPT、ART和SIRT算法及其组合算法的成象能力结合本实验的条件进行了性能的比较。结果表明:BPT方法计算简单,但所给出的图象偏粗糙,特别是在井间观测时水平方向的分辨率较差;ART方法收敛速度快,图象的反差大,反演结果的高频成分比较丰富,但边界异常大;SIRT方法收敛速度慢,图象的低频比较好,结构平滑,边界异常小,反演受测量误差的影响比较小。综合了ART和SIRT两种反演方法特点的ART-SIRT联合反演方法比较理想,发挥出了各自的长处。 相似文献