首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 78 毫秒
1.
CT图象重建的算法优化和代码优化   总被引:4,自引:1,他引:3  
卷积反投影是CT图象重建的最重要算法作为一种线性处理,其基本算法并不复杂,但由于CT的数据量庞大,使得图象重建的计算十分耗时.国外CT大都采用基于64位CPU的工作站并配以专门为CT建象设计的阵列处理板,这种昂贵的造价数十万元的计算机系统对一般医院用户来说是很大的经济负担.相反我们一直采用通用PC工作站,不仅降低了成本,而且计算机可以得到迅速的升级.这样,为满足医院临床诊断和治疗对建象速度的需要,我们就必须对该算法进行优化,优化的效果也成为CT软件能否适应市场需要,能否最终产品化的关键.本文总结了针对CT图象重建分别在算法层次和代码层次上所做的优化工作.第一,在算法实现的层次上,(1)在反投影算法中,进行内外循环交换,将最内层循环的计算量最大限度地减少,(2)利用反投影算法的对称性,只计算出部分图象数据并对称得到其它图象数据,(3)根据不同的视野范围自动调整反投影的计算量,(4)充分利用FFT的特点,提高实数线性卷积的效率.(5)数据的行列交换,用于提高CPU对数据的存取速度.第二,在代码编程的层次上,充分利用新一代微机CPU--Pentium III提供的新资源、新指令和相应的代码优化工具软件,重新编写算法关键部分的程序代码,以提高程序的实际执行效率.在东大阿尔派全身CT扫描机CT-C2000上的实际测试表明,通过上述算法优化和代码优化,CT图象重建软件的运行速度提高了8倍以上,己达到了医院实际使用的要求.  相似文献   

2.
唐巍  王彦飞 《地球物理学报》2018,61(11):4598-4612
利用X射线计算机断层成像(CT)方法对页岩的扫描成像是一种无损的,对研究页岩微纳孔隙结构有重要意义的方法.传统的CT重构通常使用的是显式的滤波反投影(Filtered Back Projection,FBP)方法,该算法具有较快的成像速度,但常伴随有伪影或不稳定等情况.对于纳米CT而言,可以通过迭代优化的方法对投影数据进行成像,传统的迭代成像有收敛速度慢导致的计算时间长等缺点.Kaczmarz算法作为一种重要的代数重建技术(ART),由于其几何意义明显,操作容易等优点,在CT重构中起着重要的作用,我们可以通过块状迭代或随机迭代的方式对其收敛速度进行改进.对于所求解问题的不适定性,代数重建过程中需要引入正则化的技巧来改善解的稳定性.本文根据实际问题的需要,使用页岩数值模型,验证了正则化Kaczmarz方法的有效性,并对重庆漆辽龙马溪组页岩样品的实际数据进行了处理,得到了较好的效果.  相似文献   

3.
本文基于重力梯度张量密度反演基本理论,建立了模型约束正则化密度反演矩阵方程.分析了代数重构算法(ART)中迭代初始值、松弛因子和终止条件三个关键参数的影响;与最小二乘求逆法对应比较分析了算法的时间和精度.结果表明:在地震、地质等地球物理手段提供初值、边界等约束较多的情况下,ART可以克服方程的不适定进行直接求解,并且合理的松弛因子和终止条件可有效提高反演效率.当初始信息不足时,添加光滑假设、深度加权等模型约束,正则化方程可以提高反演结果的可靠性.ART的行迭代可有效避免观测误差的积累和矩阵求逆的计算,从而使计算精度和速度提高数倍.最后基于GOCE地球重力场模型所得重力梯度,以地震层析成像所得速度模型为约束,对华北克拉通密度结构进行了反演,并与该区已有密度研究结果进行了对比.结果表明:利用GOCE重力场系数计算重力梯度扰动,以速度模型为约束,基于代数重构算法进行重力梯度反演所得密度模型与重力-地震联合反演所得密度模型具有很好的对应性.ART算法为重力梯度张量反演中大规模复杂问题的快速计算提供了又一种有效手段.  相似文献   

4.
大地电磁测深二维反演是现有MT数据处理解释的主流,而反演的速度依赖于正演的计算精度和速度.本文在简要介绍CUDA并行计算模型的基础上,对MT二维有限单元数值模拟进行了分析,采用基于CUDA并行技术和有限元数值模拟技术,实现了基于CUDA并行计算的大地电磁二维有限元数值模拟方法.通过一个三层模型的CDUA并行算法计算结果与解析解对比,及3个复杂模型的CDUA并行算法与CPU串行算法计算结果对比分析,结果表明,这种新的数值模拟方法在保证计算精度的同时能显著提高计算速度,加速比可达20多倍.  相似文献   

5.
Feldkamp-type-VOI锥束CT重建算法加速的研究   总被引:1,自引:0,他引:1  
本文具体介绍了一种Feldkamp-type感兴趣区域锥束CT重建方法.首先,我们分析了算法的特点:继而利用锥束扫描模式存在的对称性的反投影优化性质进行算法优化;最后,通过Shepp-Logan数值模拟实验显示优化方法可有效提高图像重建的速度且几乎不影响图像质量.  相似文献   

6.
多层螺旋CT中,采用大螺距可以大大提高扫描速度,在工业和医学中颇受青睐。但在给定探测器尺寸条件下,螺距过大会出现明显的伪影。Noo等1999年提出的单层重排算法(SSRB)能够适应较大的螺距,且有较高的重建图像质量。本文采用新的重排方法对SSRB算法进行了改进,理论分析及仿真实验表明,相同探测器大小时,改进的SSRB算法(ISSRB)可适用的螺距是SSRB的(1+tan^2δ)倍(δ射线源半张角),且有相近或更高的图像质量。  相似文献   

7.
在实际的地球物理数据采集工作中,会因为多方面的客观原因导致数据缺失,对缺失数据进行插值重构是地球物理数据处理和解释的基础问题。基于地球物理数据自身或在变换域内的稀疏性,将地球物理数据的重构转化为稀疏优化问题可提高数据重构的精确度与稳定性。本文建立了LO范数最小化的地球物理数据稀疏重构模型,针对不同规模、不同特征的地球物理数据引入了两种不同类型的LO范数最小优化问题的近似求解算法,即基于LO范数最小化的迭代再加权最小二乘算法与具有快速收敛性的快速迭代硬阈值法。理论分析与数值试验表明,将迭代再加权最小二乘算法应用到位场数据重构中可发挥其收敛速度快,计算时间短,精度高的优势,而快速迭代硬阈值法更适合处理地震数据,相对于传统的迭代硬阈值法计算效率有了很大的提高。  相似文献   

8.
基于Abaqus软件的并行计算集群平台构建与优化方法   总被引:1,自引:0,他引:1  
根据有限元显式算法和隐式算法的特点,研究了岩土工程动力分析并行计算集群系统的硬件要求、集群系统的构建方法,构建了基于EM64T硬件构架、双路Intel Xeon处理器、Linux操作系统和64位Abaqus软件的32CPU并行计算集群平台,测试了存储子系统对集群性能的影响,比较了两种千兆以太网络作为、集群子网络的性能优化方法。以地下结构的地震反应分析为例,测试了优化前后该集群系统的计算速度,发现两种以太网络性能优化方法都可以有效提高集群计算速度。列举了该集群系统在深水桥梁基础流固耦合动力分析、地下结构地震反应分析和快速轨道交通环境振动分析中的应用,显示了该集群在显式算法、隐式算法及小规模、大规模数值计算问题中的并行计算效率,证明所构建的Abaqus数值模拟并行计算集群平台能够满足计算规模、计算精度和时效性的要求。  相似文献   

9.
基于工业CT测量数据的逆向建模技术广泛应用于军事、汽车制造、航空、航天等领域。目前,在基于工业CT点云数据的逆向重构精度验证和标定上,国内外缺少相关研究。本文提出了一种用于标定工业CT设备点云提取精度的配准及评价算法。基于CT扫描和三维光学扫描点云数据,在研究多种配准算法的基础上,总结手动选点粗配准和自动全局精配准相结合的配准算法:手动选择三个匹配点,利用坐标变换完成粗配准;用ICP算法完成自动精配准,并利用修正的Hausdorff距离为目标函数来评价配准结果。对基于七参数法和ICP算法的计算步骤做了详细说明,并赋以实例。最后,以自行设计的标准件为例,获得CT测量设备点云提取精度标定值,验证配准流程的适用性。   相似文献   

10.
波动方程数值模拟在油气勘探中发挥着重要作用,如何提高它的计算效率一直是人们研究的课题.目前多核计算机已非常普及,而现有程序基本都是采用MPI实现并行计算的.这种进程级粒度的并行计算方式在PC-Cluster之类的分布式计算机上效果很好,可是在单个节点上却受内存限制,往往只能使用少数几个甚至单个计算核心,多核处理器的效能难以得到有效发挥.本文基于MATLAB科学计算平台构建波动方程数值模拟算法,将占主要计算量的波场外推式计算分解为矩阵-向量乘法、向量对应元素相乘或相减运算,通过MEX文件机制,在MATLAB中引入了OpenMP多线程并行计算,解决了MPI进程级并行在内存使用受限的情况下,多核利用效率低的问题.在四核计算机上测试表明,右端项计算平均加速3.37倍,解对角线方程平均加速1.66倍,波场外推式的计算平均加速3.11倍,使正演计算的总体计算速度提高了近3倍,有效提高了计算效率.  相似文献   

11.
12.
正SCIENCE CHINA Earth Sciences,an academic journal cosponsored by the Chinese Academy of Sciences and the National Natural Science Foundation of China,and published by Science China Press and Springer,is committed to publishing high-quality,original results in both basic and applied research.  相似文献   

13.
正SCIENCE CHINA Earth Sciences,an academic journal cosponsored by the Chinese Academy of Sciences and the National Natural Science Foundation of China,and published by Science China Press and Springer,is committed to publishing high-quality,original results in  相似文献   

14.
In natural waters arsenic concentrations up to a few milligrams per litre were measured. The natural content of arsenic found in soils varies between 0.01 mg/kg and a few hundred milligrams per kilogram. Anthropogenic sources of arsenic in the environment are the smelting of ores, the burning of coal, and the use of arsenic compounds in many products and production processes in the past. A lot of arsenic compounds are toxic and cause acute and chronic poisoning. In aqueous environment the inorganic arsenic species arsenite (As(III)) and arsenate (As(V)) are the most abundant species. The mobility of these species is influenced by the pH value, the redox potential, and the presence of adsorbents such as oxides and hydroxides of Fe(III), Al(III), Mn(III/IV), humic substances, and clay minerals.  相似文献   

15.
16.
17.
18.
正SCIENCE CHINA Earth Sciences,an academic journal cosponsored by the Chinese Academy of Sciences and the National Natural Science Foundation of China,and published by Science China Press and Springer,is committed to publishing high-quality,original results in both basic and applied research.  相似文献   

19.
正Director:Shangfu Kuang,China Vice-directors:Chunhong Hu,China Duihu Ning,China Guangquan Liu,China The International Research and Training Center on Erosion and Sedimentation(IRTCES)was jointly set up by the Government of China and UNESCO on July 21,1984.It aims at the promotion of international exchange of knowledge and cooperation in the studies of erosion and  相似文献   

20.
Feedback mechanisms, which operate upstream through drawdown and backwater effects and downstream through sediment discharge are responsible for channel evolution. By combining these mechanisms with channel processes it euables a dynamic process-response model to be developed to simulate the initial evolution of straight gravel-bed channels. When erosion commences on a land surface, sediment entrained in the headwater reach by hydraulic action is selectively transported, deposited and reworked. This produces a damped oscillation between degradation and aggradation as the channel and valley respond to spatial and temporal variations in sediment calibre and hydraulic conditions. The initial cut and fill phases are responsible for valley incision and floodplain development while secondary and subsequent activity can produce river terraces. Eventually sediment entrainment in the headwaters declines as slopes are reduced. Subsequent channel evolution is relatively insignificant because it is dependent on local weathering activity producing material that can be transported on declining slopes. Therefore landforms produced during the initial phase of development, when local weathering was non-limiting, dominate the landscape.  相似文献   

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

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