基于约束共轭梯度的高能闪光照相3维重建算法
利用空间正交约束的相机自标定和三维重建
是二次曲线 C 在成像平面上的投影 , v1 、v2 、v3 是图
像上互相垂直的 3 个方向的消失点 , 则它们形成了
一个自极三角形 ,如图 2 所示 ,且满足下列条件
v1T wv2 = v1T wv3 = v2T wv3 = 0
(6)
式 (6) 为关于 w 的 3 个线性方程组 , 因此可由空间 3
直线将在像平面相交于一点 ,称为消失点. 移动相机
使视角发生变化 ,消失点位置也随之发生变化 ,因此
可由消失点的位置变化计算出相机的参数.
11 2 消失点位置测算
正确测算消失点的位置是相机标定的前提[5] ,
常用的 Ho ugh 变换直线提取算法适合轮廓分明的
几何结构. 本文采用 Burns[6] 相位编组的改进算法 ,
在不同方向采用不同阈值的线段长度进行过滤 , 去
除那些由大纹理造成的短直线.
在提取的直线束中 ,以直线的倾角θ为 x 轴 ,相
同倾角的直线条数为 y 轴作直方图 , 找到 3 个局部
最大角度θi ,去除其他方向的直线段 ,得到对应 3 个
空间正交方向的 3 组直线束 , 这些提取出的直线束
就是空间景物中需要找到的平行直线.
相互正交的单位矢量 , 在成像平面上检测出的 3 个 消失点分别为 ( xi , yi ) i = 1 ,2 ,3 , 由于 vi 相互正交 , 则满 足 vi vj = 0 ( i ≠j) ,即
Abstract : A novel met hod for t he t hree dimensio nal (3D) reco nst ructio n f ro m two uncalibrated images was described. The camera int rinsic parameters can be estimated linearly by using t hree pair s of correspo nding vanishing point s f ro m t he t hree mut ual ort hogo nal space directio ns , t hen t he camera motio n parameter s be2 t ween t wo views can be estimated f ro m t hree gro up s of lines t hat are mut ually o rt hogo nal in 3D space. Af2 ter calculating t he camera p rojectio n mat rices , t he coordinates in 3D space can be calculated using t riangu2 lar measurement . Co mparing wit h t he active camera self2calibratio n met hod using t hree mut ual o rt hogo nal t ranslatio ns , t he camera ort hogo nal motio n co nst raint s are t ransformed into t he ort hogo nal co nst raint s of t he spatial st ruct ure in o ur app roach , so t hat it is effective to be realized wit h higher adaptabilit y and t he 3D model s of t he scene can be recovered f ro m 2D images taken by an uncalibrated handhold digital camera. This app roach has been applied to t he real images of architect ural scenes , and a 3D model of t he building was reco nst ructed wit h bet ter performance. The new images generated f ro m t he reco nst ructed 3D mo del for t he new viewpoint s are co nsistent wit h t he perceptio n of t he real scene , and t he measure error of t he plane angle bet ween t he reco nst ructed 3D model and t he real scene is wit hin 1. 5 % - 2. 6 %. Keywords : t h ree di mensional reconst ruction ; camera cal ibration ; com p uter v ision
单幅图像三维表面重建中的共轭梯度算法
单幅图像三维表面重建中的共轭梯度算法
俞鸿波;赵荣椿;王兵;江泽涛
【期刊名称】《计算机工程与应用》
【年(卷),期】2003(039)035
【摘要】从单幅图像获得物体的表面高度是计算机视觉中的一个重要研究领域,迭代算法的计算精确度高,但收敛速度较慢.该文对于几种常用的共轭梯度优化算法进行了详细分析,提出了在三维表面重建过程中实现共轭梯度算法的具体步骤和计算方法,并评价了算法的性能和优缺点.对合成图像进行仿真,并将表面恢复结果和算法收敛速度与传统的变分迭代方法比较,验证了算法的可行性和实时性.
【总页数】4页(P24-26,56)
【作者】俞鸿波;赵荣椿;王兵;江泽涛
【作者单位】南昌航空工业学院测控系,南昌,330034;西北工业大学计算机系,西安,710072;西北工业大学计算机系,西安,710072;西北工业大学计算机系,西
安,710072;南昌航空工业学院测控系,南昌,330034
【正文语种】中文
【中图分类】TP391
【相关文献】
1.OpenInventor在基于单幅图像三维重构中的应用 [J], 张晓峰;侯国柱;吴凤和
2.不标定三维欧氏重构中单幅彩色伪随机编码图像特征点提取方法研究 [J], 廖素引;吴先良;卫敏;马修水;李桂华
3.单幅图像形貌恢复在三维数字加工中的应用 [J], 宋丽梅;曲兴华;杨燕罡;叶声华
4.基于单幅图像的海南生态资源三维图像重构与建模方法 [J], 吴丽华;冯建平;何书前
5.单幅图像三维表面重建的算法研究与实现 [J], 孙林丽;李言;郑建明
因版权原因，仅展示原文概要，查看原文内容请购买。
基于影像几何约束匹配策略的三维重建方法研究与实现
摘要在计算机视觉领域，基于图像的三维重建是近年来的研究热点。
其核心是利用影像信息还原被摄场景几何结构，通过图像投影获得真实纹理，构建三维模型。
图像间匹配点对(同名点)的质量最终决定三维模型的完整性和精度。
本文在现有研究基础上，研究无定姿影像提高匹配点对数量与正确性的方法，利用图像几何约束恢复被摄物体空间坐标，达到更好地三维场景重建的目标。
主要工作内容如下：(1)高分辨率影像快速特征提取与初匹配针对高分辨率影像在SIFT特征提取过程中，计算规模过大、效率低下的问题，通过改进并行分块策略缩小数据规模，实现单张影像104~105数量级特征点在102秒内的提取。
针对随之大大增加的特征匹配计算量，分析特征点尺度分布规律，提出大量特征点的尺度筛选方法，在10秒内快速建立影像之间可靠的初步匹配。
(2)基于特征点尺度分布、对极约束与视差分布的密集匹配方法大量特征点增大了相似特征出现的概率，导致仅通过特征匹配建立的匹配点对数量十分有限。
根据初步匹配得到的信息，建立两视图对极约束、特征尺度比例约束与已有匹配点附近的视差约束，缩小特征匹配的搜索范围。
采用并行方法加速匹配过程，实现数百秒内104数量级的特征点密集匹配。
(3)改进以三视图为单元的图像三维重建方法两视图对极约束将匹配点搜索范围限制在对极线中，仍可能造成误匹配。
本文进一步将三视图对极约束约束用于扩展当前匹配点集合。
根据影像属性数据建立欧式空间中的三视图投影矩阵，并反向投影生成被摄对象三维点云。
对大范围地形数据，提出平面投影空间三角网生成方法，简化建模过程。
形成针对无定姿影像的以三视图为单元，快速构建结构与纹理信息丰富的三维模型的方法。
(4)镜面成像等特殊光环境的处理方法分析镜面反射、光场强度不均、折射与遮挡等特殊情况对建模过程的干扰，研究满足现有三维重建计算假设的镜面反射，推导相关数学模型，减少其对三维重建结果的负面影响，并利用其提供的新约束优化模型。
(5)无定姿影像三维重建软件原型系统研发研发无定姿影像三维重建软件原型系统，并进行实验。
基于结构光的三维图像重建方法、装置及存储介质[发明专利]
![基于结构光的三维图像重建方法、装置及存储介质[发明专利]](https://img.taocdn.com/s3/m/177717f45ebfc77da26925c52cc58bd63186939c.png)
(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号 (43)申请公布日 (21)申请号 201711344608.6(22)申请日 2017.12.15(71)申请人 中国科学院深圳先进技术研究院地址 518000 广东省深圳市南山区深圳大学城学苑大道1068号(72)发明人 宋展 曾海 (74)专利代理机构 深圳青年人专利商标代理有限公司 44350代理人 傅俏梅(51)Int.Cl.G06T 7/73(2017.01)G06T 15/00(2011.01)G06T 3/00(2006.01)(54)发明名称基于结构光的三维图像重建方法、装置及存储介质(57)摘要本发明适用图像处理领域，提供了基于结构光的三维图像重建方法、装置及存储介质，该方法包括：提取输入的物体编码图像的主编码特征点；构建主编码特征点的拓扑网络，根据拓扑网络提取物体编码图像中包括的所有编码元素图像；对每个编码元素图像中的图形特征点的初始位置进行定位，计算编码元素图像的辅助编码特征点；根据主编码特征点和辅助编码特征点，使用深度学习网络对所有编码元素图像进行识别；根据预设的极线编码策略，将识别到得的每个编码元素图像对应的编码信息与预先存储的编码图案对应的编码信息进行匹配，以实现解码；根据预先获取的三维图像重建系统标定参数和匹配得到的解码信息对物体进行三维图像重建，以得到物体的三维图像。
权利要求书2页 说明书7页 附图4页CN 108122254 A 2018.06.05C N 108122254A1.一种基于结构光的三维图像重建方法，其特征在于，所述方法包括下述步骤：当接收到用户输入的物体三维图像重建请求时，提取输入的物体编码图像的主编码特征点，所述物体编码图像中的编码元素图形具有旋转对称性且包含预设数量的辅助编码特征点；根据所述物体编码图像的主编码特征点，构建所述主编码特征点的拓扑网络，根据所述拓扑网络提取所述物体编码图像中包括的所有编码元素图像；使用预设的角点检测算法对每个所述编码元素图像中的图形特征点的初始位置进行定位，根据所述初始位置以及所述编码元素图像的灰度值计算所述编码元素图像的辅助编码特征点；根据所述主编码特征点和所有所述编码元素图像的辅助编码特征点，使用预先训练好的深度学习网络对所述所有编码元素图像进行识别；根据预设的极线编码策略，将识别到得的每个所述编码元素图像对应的编码信息与预先存储的编码图案对应的编码信息进行匹配，以实现对应的主编码特征点及辅助编码特征点的解码；根据预先获取的三维图像重建系统标定参数和所述匹配得到的解码信息对所述物体进行三维图像重建，以得到所述物体的三维图像。
由共轭梯度得到三维磁像
由共轭梯度得到三维磁像
Pilki.,M;刘心铸
【期刊名称】《物探化探译丛》
【年(卷),期】1997(000)A02
【总页数】5页(P84-88)
【作者】Pilki.,M;刘心铸
【作者单位】不详;不详
【正文语种】中文
【中图分类】P631.2
【相关文献】
1.基于磁梯度张量的共轭梯度3D约束反演 [J], 李金朋;张英堂;范红波;李志宁;尹刚
2.基于三维等效磁网络法计算混合励磁爪极发电机负载特性 [J], 王群京;倪有源;张学;李争
3.三维地质模型中地震波共轭梯度非线性走时反演 [J], 李飞;张雪梅;姬运达;何少林
4.非线性共轭梯度三维反演在山西河津铁矿音频大地电磁数据勘察中的应用 [J], 邹宗霖; 翁爱华; 周子琨; 连鑫葆; 郭俊豪
5.一种改进的三维子空间极小化共轭梯度法 [J], 刁新柳; 刘红卫; 赵婷
因版权原因，仅展示原文概要，查看原文内容请购买。
基于共轭梯度算法的重力梯度数据三维聚焦反演研究
基于共轭梯度算法的重力梯度数据三维聚焦反演研究高秀鹤;黄大年【期刊名称】《地球物理学报》【年(卷),期】2017(060)004【摘要】重力梯度数据相对于传统重力数据,能够更细致、准确地描述地球浅部构造和研究矿产资源分布等信息.本文采用共轭梯度算法,在加权密度域求解重力梯度数据三维聚焦反演最优化问题,以恢复地下三维密度分布,目标函数包括数据不拟合函数和最小支撑稳定函数.首先,在推导目标函数对加权密度的一阶导数时,为了得到更合理的计算公式,我们考虑变加权函数中含有密度变量;此外,本文通过密度上下限约束,改善了传统聚焦反演中聚焦因子选取困难的问题.新算法获得的反演结果,对聚焦因子的选择约束较少,相比传统聚焦算法,能够更容易的获得理想结果.将方法应用于理论模型验证其有效性和正确性,并应用本文方法处理文顿盐丘地区的航空全张量重力梯度数据,得到了与已知地质信息匹配的密度分布,表明本文方法具有处理实际数据的能力.%Gravity gradient data can be used to resolve the earth shallow structure in more detail and study the distribution of mineral resources more accurately compared with traditional gravity data.In order to recover the subsurface 3D density distribution we solve the optimization problem about the focusing inversion of gravity gradient tensor data in a space of the weighted density by the conjugate gradient algorithm.The objective function includes a misfit functional between the observed data and the predicted data and a minimum support stabilizing functional.At first,when deriving the gradients of the objective function,we consider thedensity variable in the variable weighting matrix to obtain a more reasonable formula.At the same time,by performing penalization of the density we successfully overcome the problem that the inversion result is so sensitive to the focusing factor pared with the conventional focusing algorithm we can obtain the focusing results more easily.To investigate the validity and practicability,we apply the method to synthetic gravity gradient anomalies on typical models.Then we apply the inversion method to gravity gradient data from the Vinton saltdome,Louisiana,USA.The results show a good match to the known geologic information.【总页数】13页(P1571-1583)【作者】高秀鹤;黄大年【作者单位】吉林大学地球探测科学与技术学院,长春 130000;吉林大学地球探测科学与技术学院,长春 130000;国土资源都应用地球物理重点实验室,长春 130000【正文语种】中文【中图分类】P631【相关文献】1.重力和重力梯度数据联合聚焦反演方法 [J], 秦朋波;黄大年2.基于深度加权的重力梯度张量数据的3D 聚焦反演 [J], 杨娇娇;刘展;陈晓红;徐凯军3.基于地质体空间位置优化约束的航空重力梯度数据三维物性反演 [J], 张楠;吴燕冈;周帅;孙鹏飞4.基于数据空间和稀疏约束的三维重力和重力梯度数据联合反演 [J], 张镕哲;李桐林;刘财;李福元;邓馨卉;石会彦5.基于多级混合并行算法的重力梯度数据三维密度反演研究（英文） [J], 侯振隆;黄大年;王恩德;程浩因版权原因，仅展示原文概要，查看原文内容请购买。
基于共轭梯度的工业 CT 图像重建
基于共轭梯度的工业 CT 图像重建
孟颖
【期刊名称】《科学技术与工程》
【年(卷),期】2008(008)008
【摘要】为了保证不完全投影数据的重建图像质量,通过共轭梯度法对工业CT图像进行重建.通过对模拟数据和实际工件断层扫描数据进行图像重建,估计了算法的有效性.结果表明,与最速下降法相比,此算法适用于不完全投影数据的图像重建,在保证重建图像拟合度的同时,大大提高了重建速度.
【总页数】4页(P2085-2088)
【作者】孟颖
【作者单位】北京理工大学珠海学院信息科学与技术学院,珠海,519085
【正文语种】中文
【中图分类】TP391.41
【相关文献】
1.基于加权奇异值分解截断共轭梯度的电容层析图像重建 [J], 陈宇;高宝庆;张立新;陈德运;于晓洋
2.一种基于约束共轭梯度的闪光照相图像重建算法 [J], 景越峰;刘瑞根;董维申
3.基于信赖域技巧的共轭梯度方法的图像重建算法研究 [J], 李晓杰
4.基于改进共轭梯度法的ERT图像重建 [J], 李秀艳;韩倩;汪剑鸣;王琦;段晓杰
5.一种基于共轭梯度的正则化OT图像重建算法 [J], 金永明;王加俊
因版权原因，仅展示原文概要，查看原文内容请购买。
改进的约束共轭梯度闪光照相图像重建算法
高能 闪光 照相 是一 项利 用 ｘ射 线 的强穿 透 能力 以及 它与 物质 相互 作用 的性 质 ， 据成 像平 面上 Ｘ射 线能 根 量 分 布对 客体 及其 内部 结构 进行 定量 测 量与诊 断 的测试 技 术 。高能 闪光 照相 对研 究客 体在 爆轰 过程 中 的行为 具有 重要 意义 。美 国 、 国 、 罗斯等 建 立 了一 系列 的闪光 ｘ射 线 照相 实 验平 台 ， 法 俄 获得 了较高 质 量 的客 体 密度 测量结 果 ｌ 。为 了获得更 高 精度 的 闪光照 相密 度测 量数 据 ， 国斥 巨资建 立质 子 闪光照 相 系统 ， 期待 获得 更 】 】 美 并
算 法 的 计 算 量
１ 闪 光 照 相 图像 重 建 模 型
闪光照 相成 像 系统 如 图 １所 示 ， 像 过 程 ｘｒ 成 ＇ｖ ａ 。１ｍ ｔ 。ｉａ ｒ １ ｏ ｏ ｅ ｄｔ ｔ ｂ ｊ ｔ ｅｃ ｅｏ
Ｘ达探面像器录Ｅ 三 二ｌ 一 体的强布图收记Ｏ ＝ ＝二 １ 射空度曼辫 艰 二：… ｌ 后图测，接件 … ＝ 三 … 线间分由 到像平 Ｕ三 二 Ｇ 二 ｅ 蛙 （ ４）
Ａ 一｛ Ａ （ ２ｔ
【△ ２ ｒ
式 中 ：， ∈｛ ， ， ， ， ｝ 』 为单一 断层 的采 样 点数 。 ｉＪ ０ １ … 』 一１ ， ， ＼ ＼
２ 约 束 共 轭梯 度 重 建 算 法
该 算 法 在 约 束 共 轭 梯 度 迭 代 重 建 的 基 础 上 ， 出 了 新 的 预 优 矩 阵 选 取 方 案 ， 小 了重 建 图 像 的轴 线 噪 声 ， 用 提 减 利 松 弛 迭 代 步长 代 替 共 轭 梯 度 法 中 的最 优 迭 代 步 长 ， 采 用 了 新 的收 敛 准 则 ， 保 证 算 法 收 敛 的 同 时 ， 少 了 重 并 在 减 建 算 法 的计 算 量 。数 值 试 验 表 明 ， 与传 统 约束 共 轭 梯 度 重 建 算 法 相 比 ， 进 算 法 稳 定 收 敛 ， 改 迭代 速度 更 快 ， 能 并 有效提高重建质量 。
改进的预处理共轭梯度快速算法在三维重力梯度数据反演中的应用（英文）
改进的预处理共轭梯度快速算法在三维重力梯度数据反演中的应用（英文）王泰涵;黄大年;马国庆;孟兆海;李野【期刊名称】《应用地球物理：英文版》【年(卷),期】2017(14)2【摘要】随着全张量重力梯度(FTG)测量技术的不断发展,重力梯度数据的三维反演技术在油气和矿产勘探中日益受到重视与关注。
为了快速处理和解释大规模的高精度数据,图形处理器GPU(Graphics Processing Unit)和预处理分解技术(Preconditioning methods)在地球物理反演中的使用变得十分重要。
本文结合对称逐次超松弛(SSOR)技术与不完全乔列斯基分解共轭梯度算法(ICCG)提出改进的预处理共轭梯度法,并考虑到方法预处理分解占用额外的时间,开发该算法的GPU 并行算法来提高加速效果。
然后通过含噪的模型数据反演来证明改进的并行预处理方法在三维全张量重力梯度数据反演中的适应性。
由此,基于NVIDIATesla C2050 GPU的并行SSOR-ICCG算法和在2.0GHz CPU上的串行程序比较,达到了大约25倍的加速比。
最后,我们将该算法应用于美国路易斯安那州南方Vinton盐丘的实测航空重力梯度数据反演中,反演出良好的反演结果,验证了该方法在三维全张量重力梯度数据快速反演中的优势和可行性。
【总页数】14页(P301-313)【关键词】全张量重力梯度(FTG);ICCG方法;共轭梯度法;重力梯度反演;CPU和GPU【作者】王泰涵;黄大年;马国庆;孟兆海;李野【作者单位】吉林大学地球探测科学与技术学院;天津航海仪器研究所;吉林省电力勘测设计院【正文语种】中文【中图分类】P【相关文献】1.基于共轭梯度算法的重力梯度数据三维聚焦反演研究 [J], 高秀鹤;黄大年2.预处理共轭梯度法在VVP三维风场反演中的应用 [J], 张少波;胡明宝;张鹏3.基于改进的预处理共轭梯度法和三维插值法精确和快速解算GRACE地球重力场[J], 郑伟;许厚泽;钟敏;员美娟4.基于多级混合并行算法的重力梯度数据三维密度反演研究（英文） [J], 侯振隆;黄大年;王恩德;程浩5.利用Extrapolation Tikhonov正则化算法进行重力梯度三维密度反演（英文）[J], 刘金钊;柳林涛;梁星辉;叶周润因版权原因，仅展示原文概要，查看原文内容请购买。
高能X光照相中密度重建的共轭梯度算法
高能X光照相中密度重建的共轭梯度算法
许海波;魏素花
【期刊名称】《计算物理》
【年(卷),期】2006(23)2
【摘要】通过对高能X光照相过程的物理分析,得到了辐射照相的点扩展函数,给出了估价函数的表达式.并以典型的法国试验客体模型为例,将共轭梯度算法应用到高能X光照相的密度重建中,得到了较满意的结果.这种方法从模拟辐射照相过程出发,将辐射照相的模拟结果与实际测量结果做比较并找出最大可能解.可在一定程度上克服其它密度重建算法通过反卷积过程消除模糊效应所引起的不确定性.
【总页数】7页(P144-150)
【关键词】X光照相;点扩展函数;估价函数;共轭梯度算法;密度重建
【作者】许海波;魏素花
【作者单位】北京应用物理与计算数学研究所
【正文语种】中文
【中图分类】O443
【相关文献】
1.共轭梯度算法在高能X光照相密度重建中的应用 [J], 许海波;魏素花
2.改进的约束共轭梯度闪光照相图像重建算法 [J], 景越峰;刘军;管永红
3.基于约束共轭梯度的高能闪光照相图像复原算法 [J], 景越峰;张小琳;管永红;刘瑞根
4.基于约束共轭梯度的高能闪光照相3维重建算法 [J], 景越峰;刘瑞根;管永红;张小琳;危才华
5.一种基于约束共轭梯度的闪光照相图像重建算法 [J], 景越峰;刘瑞根;董维申因版权原因，仅展示原文概要，查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第 ｌ 卷 第 ５期 ９ ２０ ０ ７年 ５月
强 激 光 与 粒 子 束
ＨＩ ＧＨ Ｐ ＯＷ ＥＲ ＬＡＳ ＥＲ ＡＮＤ ＡＲＴＩ Ｐ ＣＬＥ ＢＥ Ｍ Ｓ Ａ
ＶＯＩ １ Ｎｏ ５ ． ９， ．
Ｍａ ２ ７ ｙ， ００
【 ａ１ ａｎｃ ａ、 ｎ Ｓｌ ｅ ｏｌｒ ｐｒ ． ｉ ｏｌ Ｌ ｏｉ ｅｃ ｉ
，ｃｒ ｓ ｇｍ ｏｓ
一
ｎｃ
＼
／
’…
－＿ １
Ｉ
～． 一 ） ｌ一曩 ＜ 二 ～一ｒ 一 ～’ 一一 ～！ 一 Ｌ ～ ｜
・
。 ｄ ｅ ｃ ｃ ｃ ’
截面 ， 然后 通过 累加各 个 断层 截面 而得 客体 ３ 数据 。而 ３维 算法 具 有将 客 体 空 间数 据整 体 性 一 次 重建 出来 维
的优 点 ， ３维重 建方 程 规模很 大 ， 但 因而计算 量很 大 。 目前 存 在 的 高能 闪 光 照相 图像 重 建 算 法都 是 ２维 的 ， 本 文提 出了一种 基 于 ３ 理论 的 ３维 约束共 轭 梯度 算法 （ Ｄ Ｃ Ｇ） 维 ３ － Ｃ 。根 据单 轴 闪光 照相 的特 点 ，Ｄ Ｃ Ｇ采 用 同 ３ －Ｃ 心 圆环 网格有 效解 决 了 ３维算法 计 算量 太 大的缺 点 ， 引入 了相邻 断层 重建 数据 间 的光滑 性约 束 ， 并 有效 克服 了
ｍ
＼ ＼
／ ／
Ｆｇ １ Ｒａ ｉｇ ａ ｈ ｙ ｔ ｍ ａ ｄ ｔ ｅｓ ｃｉｎ ｏ ｈ ｂ ｅｔ ｉ． ｄｏ ｒ ｐ ｙｓ ｓｅ ｎ ｈ ｅ ｔ ｆ ｅｏ ｊ ｃ ｏ ｔ
图 １ 闪光 照 相 系统 及 客 体 截 面 网格 简 图
＊ 收 稿 日期 ：０ ６０ —９ ２０ —９２ ； 修 订 日期 ：０ ７０ —１ ２ ０—２０ 基 金 项 目 ： 国 工 程 物 理 研 究 院 双 百 人 才基 金 资 助 课 题 中 作者简介 ： 景越 峰 （ ９Ｏ ） 男 ， 龙 江 肇州 人 ， 士 ， 事 图 像处 理 研究 工 作 ； ｙ＠ ｅｏ ．ｏ 。 １ ８一 ， 黑 硕 从 ｊｆ ｙ ｕ ｃｒ ｎ
重建 图像 的轴 线噪声 和 伪影 。数 值试 验表 明 ，Ｄ－Ｃ ３ Ｃ Ｇ具 有 良好 的 收敛性 、 定性 以及 较高 的重 建精 度 。 稳
１ 闪光 照 相 ３维 重 建 模 型
考虑 闪光 照相 系统 为单轴 照 相 ， 因而假 定被 照 客体 是轴 对 称 的 。图 １ 闪光照 相 系统 及 客 体 同 心 圆环 网 是 格截 面简 图 。对 于轴 对 称客体 ， 们定 义 客体存 在 Ｍ 个 重建 截 面 ， 每个 截面 上 ， 我 在 重建 模 型采用 等 间距 的 同 心 圆环 网格 ， 每一个 截 面分 为与 采 样 点 数 相对 应 的 Ｎ 层 ， 每层 的 间距 为 ２ ｒ Ａ 。这样 闪光 照 相 系统 在 单 能 、 模 无 糊 、 散射 的情况 下 ， 图像 ３维重 建 的代数 方程 组 为 无 其
维普资讯
强 激 光 与 粒 子 束
ｚ— Ａ
第 １ ９卷
（） １
式 中 ： 为 测量得 到 的 ＭＮ 维 量化 逐层 排 序 的投影 图像 矢量 ； ｚ 为所 求 的 ＭＮ 维 逐层 排序 的客 体空 间线 吸 收系 数 分 布矢 量 ； 平行光 束 投影下 ， Ｘ射 线穿 过客 体形 成 投影 图 像 的正 向 ３维 投 影 矩 阵 （ 在 Ａ为 ＭＮ×ＭＮ 的块 对
有 滤波 反投影 重建 （ Ｂ ） ］正 则 化 预优 约 束 共 轭梯 度 重 建 （ Ｐ Ｃ ［ 等 。Ｆ Ｐ算 法重 建 速 度 快 ， 重 建 图 Ｆ Ｐ 口、 Ｒ Ｃ Ｇ）。 Ｂ 但
像 分辨 率较低 ， 噪声抵 抗 能 力 较 差 。Ｒ Ｃ Ｇ 算 法 可 以较 好 地 抑 制 图像 噪声 ， 对 ＰＣ 但重 建 图像 轴 线 仍 有 较 大 噪 声 ， 重建 中引入 的正则 化准 则使 得重 建 图像存 在一 定 程度 的模 糊 。 且 ３ 重建 是 目前 图像重 建 中十分 活跃 的一个 领域 。一 般 ２维算 法 的原理 是 首先 重 建 出客 体各 个 ２维 断层 维
摘
要 ： 针 对 高 能 闪 光 照 相 系 统 成 像 信 噪 比低 的 特 点 ， 出 了 ３维 约 束 共 轭 梯 度 算 法 ， 算 法 基 于 ３维 提 该
重 建 理 论 来 重 建 方 程 ， 用 约 束 共 轭 梯 度 法 迭 代 求 图像 重 建 的 最 优解 。 数 值 模 拟 结 果 表 明 ： 用 ３ 约 束 共 轭 利 采 维
梯度算 法能有效 改善高能闪光照相重建 图像轴 线噪声 大的 问题 ． 一种具有很 高抗噪 能力 的闪光 照相 图像 重 是
建算 法。
关键词 ： 约束共 轭梯度 ； ３ 图像 重建； 闪光照相 ； 数值模拟 ； 信噪比 维
中 图分 类 号 ： Ｔ ３ １ Ｐ ９ 文献 标 识 码 ： Ａ
闪光 照相是 一项 用途 广 泛的无 损 检测 技 术 ， 利用 Ｘ射 线 的穿 透 能力 及 其 与 物质 相 互作 用 的性 质 ， 对 它 来
客体 内部结 构进 行诊 断 。高能 闪光 照 相的重 要 目标 是重 建客 体 的空 间密 度分 布 。由于受 到模 糊 、 散射 、含有 很大 的 误差 ， 严 重影 响 了图 像重建 质 量 。 目前 针对 高能 闪光照 相 的主要 重 建方 法 这
文 章 编 号 ： １０ —３ ２ ２ ０ ） ５０ ６ —５ ０ １４ ２ （０ ７ ０ —８ ３０
基 于 约束 共轭 梯 度 的高 能 闪光 照 相 ３维 重 建 算 法
景越峰， 刘瑞根 ， 管永红， 张小琳 ， 危才华
（ 国工 程 物理 研 究 院 流 体 物 理 研 究 所 ，四 川 绵 阳 ６ １ ０ ） 中 ２ ９ ０