首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 93 毫秒
1.
鉴于实际工程应用需求,基于截断投影数据的CT图像局部重建是X射线CT重建的一个重要课题,本文概述了现有的局部重建算法,分析了其优缺点,提出了一种基于改进重建滤波器的局部CT重建算法。通过对滤波函数的改进,使其能量主要集中在中心位置,使旁瓣迅速下降,进而使滤波函数具有局部重建功能。另外改进的滤波函数压制了低频信息,突出了高频信息,具有提高成像空间分辨率的功能,但对图像噪声有所放大,实际实验数据验证了算法的可行性。  相似文献   

2.
本文阐述了用一种新型滤波函数进行CT图象局部重建的原理,并给出了这种重建的二次误差估计.本文选用一种新的滤波函数,它的旁瓣衰减很快,经证明在可能选取的滤波函数中,它的衰减是最快的.新滤波函数由于空域的迅速衰减性,因此具有良好的局部重建性质.同时保持该滤波函数的其它特性,能够有效地抑制噪声,补偿混迭效应,并缓解Gibbs效应.应用新滤波函数,忽略远离计算点的数据,从而直接用CBP实现图像局部重建.在计算机上对模拟和实测数据进行局部重建的实验结果表明本文建议的图像局部重建算法简便快速,重建图象Gibbs效应减小并与用全部数据的CBP重建图像比较仍具有较高的空间及密度分辨率.  相似文献   

3.
CT断层重建中滤波函数设计的新方法   总被引:1,自引:0,他引:1  
滤波函数在CT图像重建过程中起着非常重要的作用,直接关系到重建图像质量。为了提高CT重建图像质量,本文从加权平均的思想出发,根据FBP重建算法的理论基础,提出一种设计滤波函数的新思路,并分析了五点加权平均滤波函数性能比较。最后针对Shepp-Logan模型数据和实际的海螺投影数据,设计出一种新的滤波函数并与S-L和R-L滤波函数的重建效果进行了比较。从比较结果可以看出,新设计的滤波函数重建的图像效果在整体性能上最好,在局部地方,其密度分辨率有所提高。本文为滤波函数的设计提出了一种新的想法和思路。  相似文献   

4.
本文分析并设计了用于CT图像重建的几种预处理滤波器。通过分析不同滤波函数对断层图像重建质量的影响,从几种经典的滤波函数出发,进行改进并与窗函数结合得到了可用于CT图像重建的几种滤波器。有关投影数据的图像重建结果表明,所设计的几种滤波器较好地满足了实际图像重建的要求。  相似文献   

5.
CT图像重建中滤波函数的优化   总被引:1,自引:0,他引:1  
CT图像重建的解析类经典算法为滤波反投影(FBP),其中的滤波环节对于重建质量的影响至关重要,通常情况下滤波函数都是将斜坡函数直接加矩形窗,而理想的矩形窗是造成Gibbs现象的根源所在,Gibbs现象使重建图像出现震荡,不够平滑.根据无穷级数求和理论,将函数展开成级数部分和形式后,可以进一步求解出部分和的算术平均值,该函数即是具有更优收敛性的费耶核(Fejes),本文根据这一理论,提出了一种改进的滤波方法,尝试利用费耶核较好的收敛性来克服Gibbs现象.实验结果表明该方法可以获得更加平滑的结果,对于图像质量的改善取得了明显的效果.  相似文献   

6.
CT图像重建的解析类经典算法为滤波反投影(FBP),其中的滤波环节对于重建质量的影响至关重要,通常情况下滤波函数都是将斜坡函数直接加矩形窗,而理想的矩形窗是造成Gibbs现象的根源所在,Gibbs现象使重建图像出现震荡,不够平滑。根据无穷级数求和理论,将函数展开成级数部分和形式后,可以进一步求解出部分和的算术平均值,该函数即是具有更优收敛性的费耶核(Fejes),本文根据这一理论,提出了一种改进的滤波方法,尝试利用费耶核较好的收敛性来克服Gibbs现象。实验结果表明该方法可以获得更加平滑的结果,对于图像质量的改善取得了明显的效果。  相似文献   

7.
选用四种不同土质试样设计相关试验,利用直接法和间接法结合不同扫描条件下的CT数变化曲线,研究扫描电压、对准模式、滤波函数、试样尺寸等CT扫描参数对岩土试样CT图像的具体影响,研究结果表明:扫描电压和滤波函数对岩土CT图像和CT数影响较大;对准模式和重建矩阵对试验结果影响非常小;滤波函数对岩土CT图像的增强或抑制作用取决于所选择的滤波函数的属性,而与试样的土质无关;最后给出了CT扫描参数选择的原则,借助于合理的CT扫描参数,从而提高CT岩土图像的质量,获得相对准确的岩土CT值。  相似文献   

8.
为了减少X射线对人体的伤害,局部图像重建成为人们研究的重点之一.本文研究了基于一种特殊窗函数的局部图像重建算法.通常,直接用局部区域的投影数据重建局部图像,会使重建图像产生常数偏移,我们将局部区域边缘的投影数据延拓为相应沿径向未知投影数据的方法,改进了局部重建图像,并给出了基于特殊窗函数重建时,该方法的误差分析.另外,本文将该算法应用于扇形束图像重建,并通过数值实验验证了此特殊窗函数应用于扇形束重建时同样有效.  相似文献   

9.
在传统的CT系统中,系统的硬件成本和计算量都是非常巨大的.本文深入分析了一种利用偏移放置的面阵探测器的锥束CT系统,在这种系统中X射线束仅仅覆盖被扫描物体的一半体积,投影数据在探测器方向上是截断的,探测器尺寸和投影数据量都减少为传统CT系统的一半.在这种扫描方式下,现有的重建方法是首先利用重排算法获得180度范围内的平行束投影数据,然后再利用滤波反投影(FBP)算法重建出物体的三维图像.但是重排算法不可避免地会引入误差,降低重建图像的空间分辨率.本文提出了一种反投影滤波(BPF)形式的直接反投影重建方法,该方法不需要对投影数据重排,直接反投影滤波重建出最终的图像.因此,该算法在数学上更简洁,计算速度更快,能够更多地保留重建图像的高频信息.最后,数值模拟实验结果验证了该系统和重建算法能够获得高质量的CT图像.  相似文献   

10.
本文基于锥束CT滤波反投影重建的FDK算法,通过两种算法改进并结合基于共享内存的OpenMP并行技术和代码优化,实现了锥束CT图像的快速重建.基于锥束CT实际投影数据的重建结果表明,图像重建速度得到了较大的提高,断层图像重建质量与FDK原型算法相当.  相似文献   

11.
In the linear filter method of interpreting resistivity sounding data, as developed by Ghosh (1971), it appears that the filter function in the x-domain approaches an oscillating function for both large positive and large negative abscissa values. In the present note the reason for this oscillating behaviour is derived, and a possible practical application is indicated.  相似文献   

12.
三维电磁偏移数值滤波器实现及参数分析   总被引:2,自引:1,他引:1       下载免费PDF全文
基于满足扩散方程的电磁偏移场解析式,推导了能够应用于实际数据处理的电磁偏移数值滤波器.为了合理选取滤波器参数,本文基于偏移微分方程,构建了满足偏移微分方程的解析测试场,通过测试场与偏移数值滤波器处理结果比较,重点研究了滤波器窗口大小、采样距离、偏移深度等滤波器参数对数值偏移滤波效果的影响,确定了数值滤波器参数选择原则.考虑到实际应用,在计算偏移电导率时,本文利用背景场定义了反射函数.建立了三维典型地电模型做数值计算,利用电磁偏移数值滤波器处理模拟数据,得到了合理的目标体位置和形态.测试场比较和数值试验证明,电磁偏移成像是一种可靠稳定的电磁资料解释技术.  相似文献   

13.
根据实际工作需要,为确保数字地震资料震相分析的准确性,基于MATLAB GUI设计了可视化滤波器。实现相应的数字低通、高通、带通、带阻滤波器,界面友好,具有可实用性。半年来,该滤波器在加格达奇地震台利用率高,发挥了较大作用。  相似文献   

14.
数字地震观测中的线性和最小相位滤波   总被引:4,自引:0,他引:4  
线性相位滤波和最小相位滤波是数字地震观测中经常都要遇到的两类滤波器。现代地震数据采集器中就预置有这两种滤波系数供用户选用。本文以EDAS—C24地震数据采集器为例,通过比较两者滤波特性的差异后得出:为什么在以震相判读分析为主要用途的数字地震观测中必须选用最小相位滤波器,最好不要使用线性相位滤波器。同时还指出,最小相位滤波对优势频谱较高的地方微震波形影响较大,且振荡尾巴较长,分析研究时应注意其影响。  相似文献   

15.
采用代数重建技术重建三维物体,提出了一种基于splatting原理的ART权值计算方法,以增强重建图像的空间分辨率.该算法基于三维锥束扫描方式,采用双三次样条局部基函数作为插值核,并利用其足印函数对物体进行插值来计算像素的权值.与常量插值相比,这种方法对物体的灰度值有更精确的近似,在相同的实验条件下可重建出分辨率高的物体.给出了适用于三维锥束重建的splatting算法的具体实现,并对实验数据进行了实物重建及误差分析.  相似文献   

16.
Koefoed has given practical procedures of obtaining the layer parameters directly from the apparent resistivity sounding measurements by using the raised kernel function H(λ) as the intermediate step. However, it is felt that the first step of his method—namely the derivation of the H curve from the apparent resistivity curve—is relatively lengthy. In this paper a method is proposed of determining the resistivity transform T(λ), a function directly related to H(λ), from the resistivity field curve. It is shown that the apparent resistivity and the resistivity transform functions are linearily related to each other such that the principle of linear electric filter theory could be applied to obtain the latter from the former. Separate sets of filter coefficients have been worked out for the Schlumberger and the Wenner form of field procedures. The practical process of deriving the T curve simply amounts to running a weighted average of the sampled apparent resistivity field data with the pre-determined coefficients. The whole process could be graphically performed within an quarter of an hour with an accuracy of about 2%.  相似文献   

17.
Kalman filtering for stochastic dynamic tidal models, is a hyperbolic filtering problem. The questions of observability and stability of the filter as well as the effects of the finite difference approximation on the filter performance are studied. The degradation of the performance of the filter, in case an erroneous filter model is used, is investigated. In this paper we discuss these various practical aspects of the application of Kalman filtering for tidal flow identification problems. Filters are derived on the basis of the linear shallow water equations. Analytical methods are used to study the performance of the filters under a variety of circumstances.  相似文献   

18.
In the field of seismic exploration, ground roll seriously affects the deep effective reflections from subsurface deep structures. Traditional curvelet transform cannot provide an adaptive basis function to achieve a suboptimal denoised result. In this paper, we propose a method based on empirical curvelet transform (ECT) for ground roll attenuation. Unlike the traditional curvelet transform, this method not only decomposes seismic data into multiscale and multi-directional components, but also provides an adaptive filter bank according to frequency content of seismic data itself. So, ground roll can be separated by using this method. However, as the frequency of reflection and ground roll components are close, we apply singular value decomposition (SVD) in the curvelet domain to differentiate the ground roll and reflection better. Examples of synthetic and field seismic data reveal that the proposed method based ECT performs better than the traditional curvelet method in terms of the suppression of ground roll.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号