核医学图像重建快速迭代算法OSEM
OSEM算法在编码板成像中的应用
２ ０ １ ６年
第 ９期
９月
Байду номын сангаас
核 电子学 与探 测技 术
Ｎｕ ｃ ｌ ｅ ａ ｒ Ｅｌ ｅ ｃ ｔ ｒ ｏ ｎｉ ｃ ｓ＆ Ｄｅ ｔ ｅ ｃ ｔ ｉ ｏ ｎ Ｔｅ ｃ ｈ ｎ ｏ ｌ ｏ ｇ ｙ
Ｖ ｏ １ ． ３ ６ Ｎ ｏ ． ９
Ｓ ｅ ｐ ｔ ． ２ ０ １ ６
“ ’ （ ， Ｙ ）＝ ／ ” （ ， Ｙ ） ￥
／ ” （ ） ｈ ， ， ， ）～ Ｙ （ ｘ ＇ ｙ ）
，
川Ｊ
］（ １ ）
９ ０ ３
式中 为第 ｋ次图像估 计值 ； ｐ （ ｚ 。 Ｙ ） 为 实测投 影值 ； ＾ （ ， Ｙ ） 为 编 码 函数 ； “ ” 表 示 卷
１ 算法简介
１ ． １ ＭＬ Ｅ Ｍ 与 ＯＳ Ｅ Ｍ 迭 代算 法
算量大 、 迭代收敛速度慢 的问题 ， Ｈ ． Ｍ ． Ｈ ｕ ｄ ｓ ｏ ｎ
等人 提 出对 Ｍ Ｌ Ｅ Ｍ 算法改进 的迭代加 速算 法， 即有序子集 一 期望最 大化（ Ｏ Ｓ Ｅ Ｍ） 算法 － ２ ］ 。
ＯＳ Ｅ Ｍ 算 法在 编 码 板成 像 中 的应 用
李汉平 ， 王锋 ， 艾宪芸
（ 防化研究院 国民核生 化灾害防护 国家重点实验室 ， 北京 １ ０ ２ ２ ０ ５ ）
摘要 ： 为提 高 辐 射编码成 像中 ＭＬ Ｅ Ｍ算法在 图像重建 时的速度 及质 量 ， 论 文基 于核 医学 Ｅ Ｃ Ｔ图 像重建 中的 Ｏ Ｓ Ｅ Ｍ算 法思想在图像重建 时对 编码成像所得 投影数 据进 行分组 ， 并 采用精 细采样 平衡相 关法的重建结果作 为初 值进行迭代 。通 过实验 对该算 法进行验 证。实 验表 明 ： Ｏ Ｓ Ｅ Ｍ 算法 可 以有 效提 高收敛速度 ， 改 善重建图像 质量 。在 辐射编码 成像中采用 Ｏ Ｓ Ｅ Ｍ算法是有效 且可行的 。 关键词 ： 编码板 ； 辐射成像 ； 图像重建 ； 最大似然期 望最 大化 算法 ； 有 序子集期 望最 大化算 法 ； Ｍ Ｃ模
迭代重建算法
迭代重建算法1. 算法概述迭代重建算法是一种用于图像处理和计算机视觉领域的重要算法。
它的主要目标是通过多次迭代来逐步改进图像的质量或者提取出感兴趣的特征。
迭代重建算法在许多应用中都有广泛的应用，如图像增强、图像去噪、图像分割等。
2. 算法原理迭代重建算法通常基于一个初始估计值，并通过多次迭代来不断改进这个估计值。
每次迭代过程中，根据一定的准则函数计算出一个更新值，并将这个更新值与当前估计值进行合并，得到新的估计值。
通过不断迭代，最终得到一个收敛的估计值。
具体来说，迭代重建算法可以分为以下几个步骤：2.1 初始化在开始进行迭代之前，需要对初始估计值进行初始化。
初始估计值可以根据具体问题进行选择，也可以随机生成。
2.2 迭代更新在每一轮迭代中，根据问题需求选择合适的准则函数来评估当前估计值和真实值之间的差距。
根据准则函数的结果，计算出一个更新值，并将更新值与当前估计值进行合并，得到新的估计值。
2.3 收敛判断在每一轮迭代后，需要判断算法是否已经收敛。
可以通过比较当前估计值和上一轮迭代的估计值之间的差距来进行判断。
如果差距小于某个阈值，则认为算法已经收敛，可以停止迭代。
2.4 结果输出当算法收敛后，最终得到的估计值就是我们想要的结果。
根据具体问题的需求，可以将结果输出为图像、特征向量等形式。
3. 算法优缺点3.1 优点•迭代重建算法能够逐步改进估计值，从而提高图像质量或者提取感兴趣的特征。
•算法具有较好的灵活性，可以根据具体问题选择合适的准则函数和更新策略。
•算法通常能够在较少的迭代次数内收敛，并且具有较高的精度。
3.2 缺点•迭代重建算法通常需要进行大量的计算，在处理大规模数据时可能会面临计算时间过长的问题。
•算法的收敛性和稳定性可能受到初始估计值的选择和准则函数的设计等因素的影响。
•算法对噪声和异常值比较敏感，可能会导致结果不准确。
4. 应用案例迭代重建算法在图像处理和计算机视觉领域有着广泛的应用。
磁共振成像技术中的图像重建算法
磁共振成像技术中的图像重建算法磁共振成像技术是一种用于观察人体内部结构的非侵入性医学成像技术。
它通过对人体内部的磁场进行扫描，可以得到高分辨率的图像信息，从而帮助医生进行诊断。
在磁共振成像技术中，图像重建算法是非常重要的一环。
它负责从扫描得到的原始数据中重建出人体内部的结构信息，并生成可视化的图像用于医学诊断。
目前，磁共振成像技术的图像重建算法主要分为两类：频域算法和空域算法。
下面将分别对这两种算法进行介绍。
一、频域算法频域算法将磁共振信号转换到频域进行处理，然后再将处理后的数据转换回时域，得到最终的图像。
其中，最常用的频域算法是快速傅里叶变换（FFT）。
它可以将磁共振信号快速地转换到频域进行处理，然后再进行反变换，得到重建后的图像。
虽然快速傅里叶变换的速度很快，但是这种算法存在一定的局限性。
例如，磁共振信号中存在很多不同频率的信号，而快速傅里叶变换对信号的不同频率处理效果不能很好地区分，从而影响图像的质量。
二、空域算法空域算法是通过对原始数据进行处理，直接得到重建后的图像。
其中，最常用的空域算法是反向投影算法。
这种算法可以将不同方向的扫描数据按照一定的规则投影到图像平面上，然后将所有的投影结果叠加起来，得到最终的重建图像。
反向投影算法的优点是可以处理不同方向的扫描数据，其中还可以添加一些先验信息，从而提高图像质量。
然而，这种算法也存在一些问题，比如有时会出现伪影情况。
此外，还有一些其他的空域算法，比如基于大脑并行矩形图像重建的算法（BART）和基于稀疏表示的重建算法（CS-MRI）。
这些算法可以在一定程度上提高图片的质量，并降低成像时间。
总结起来，磁共振成像技术的图像重建算法是非常复杂的，需要结合理论和实践进行优化。
随着计算机技术和算法的不断发展，未来有望实现更快速、更准确、更高质量的图像重建算法，从而实现更好的医学诊断效果。
核医学图像重建快速迭代算法OSEM
核医学图像重建快速迭代算法OSEM一、引言核医学影像设备如单光子断层扫描仪(SinglePositronEmissionComputeTomography，SPECT)、正电子发射断层扫描仪(PositronEmissionTomo-graphy，PET)融合了当今最高层次的核医学技术，是目前医学界公认的极为先进的大型医疗诊断成像设备，在肿瘤学、心血管疾病学和神经系统疾病学研究中，以及新医药学开发研究等领域中已经显示出它卓越的性能。
随着核医学断层影像设备的广泛应用和计算机技术的迅速发展，图像重建方法作为该类设备中的一个关键技术，其研究工作越来越受到人们的重视。
本文概述了传统的图像重建方法，并详细介绍了一种具有较高图像质量和较短计算时间的重建算法—有序子集最大期望值方法(Ord-eredSubsetsExpectationMaximization，OSEM)在核医学影像设备中的应用。
二、传统的图像重建方法在核医学影像设备中，需要根据物体某一层面在不同探测器上检测到的投影值来重建该断层图像层面，即二维图像重建。
传统的图像重建方法主要分为解析法和迭代法。
解析法是以中心切片定理(CentralSliceTheorem)为理论基础的求逆过程。
常用的一种解析法称为滤波反投影法(FilteredBack-Projection，FBP)。
FBP法首先在频率空间对投影数据进行滤波，再将滤波后的投影数据反投影得到重建断层图像。
滤波器选为斜坡函数和某一窗函数的乘积，窗函数用于控制噪声，其形状权衡着统计噪声和空间分辨。
常用的窗函数有Hanning窗，Hamming 窗，Butterworth窗以及Shepp-Logan窗。
解析法的优点是速度快，可用于临床实时断层重建。
但当测量噪声较大或采样不充分时，这类算法的成像效果不甚理想，尤其是在核医学断层图像重建中对小尺寸源的成像效果差(即所谓偏体积效应)。
在滤波中如果对高频信号不做抑制，截止频率高，此时空间分辨最好，但所重建的图像不平滑，易产生振荡和高频伪影;反之，采用较低截止频率，过多压抑高频成分的低通窗函数会造成重建图像的模糊，故在变换法中低噪声和高分辨对滤波器的要求是矛盾的，需折衷选择。
核医学显像数据重建和分析的计算方法和策略
核医学显像数据重建和分析的计算方法和策略核医学显像是一种通过将放射性同位素或放射性药物注入体内，利用放射性同位素的辐射来检测和诊断疾病的影像学方法。
核医学显像技术是一种非侵入性的检测方法，广泛应用于心血管疾病、癌症、神经系统疾病和骨科领域等的临床诊断中。
核医学显像数据重建和分析的计算方法和策略是指利用计算机技术对核医学显像所获得的原始数据进行重建和分析的技术方法和策略。
这些方法和策略的目的是提取和分析有用的信息，以帮助医生进行疾病的诊断和评估治疗效果。
核医学显像数据重建是指将从患者身上获得的放射性同位素或放射性药物的信号转化为图像的过程。
在重建过程中，需要考虑到放射性同位素或放射性药物的衰变特性、探测器的响应特性以及数据采集的方式等。
目前常用的核医学显像数据重建方法有滤波反投影算法、最大似然算法、迭代算法等。
这些算法通过对原始数据进行处理，重建出代表患者内部生物分布的图像。
在重建出图像之后，还需要对图像进行进一步的分析。
核医学显像数据分析包括定量分析和定性分析两个方面。
定量分析主要是通过对图像中不同区域的计数进行测量，来获得有关组织功能和代谢状态的信息。
定性分析则是通过对图像的形态、分布和连续性等特征进行观察和分析，来获得有关病变的信息。
在核医学显像数据分析中，还有一些常用的计算方法和策略。
例如，感兴趣区域（ROI）分析是一种常用的方法，通过在图像中选择特定的区域来提取该区域的定量参数。
另外，还有基于统计学的方法，如标准化摄取值（SUV）等，可以用来评估肿瘤的代谢活性。
此外，还有一些定量分析方法，如动力学分析和时间-活动曲线分析，可以用来研究药物在人体内的代谢和转化过程。
值得注意的是，在进行核医学显像数据重建和分析时，数据质量的保证非常重要。
因为用于重建和分析的数据具有辐射成像的特点，必须确保数据的准确性和可靠性。
同时，隐私保护也是一个重要的问题，必须确保患者的个人隐私不被泄露。
总之，核医学显像数据重建和分析的计算方法和策略是一项关键的技术，对于提高核医学诊断的准确性和可靠性具有重要意义。
临床PET图像重建与处理简介
对于一幅断层图像，将各个角度上的投影数据，反投影 叠加。核医学图像的有用信息多集中在低频区域，为消 除投影数据中的高频噪声，在反投影之前，先将投影数 据经过一个低通滤波器，保留投影数据中有用的低频信 号，去除高频噪声。
FBP算法流程
OS-EM法
迭代法求解过程: a. 假定一初始图像； b. 计算该图像投影； c. 同测量投影值对比； d. 计算校正系数并更新初始图像值； e. 满足停步规则时，迭代中止 ，否则 以新的重建图像作为初始图像从b 步开始。
PET图像重建方法与途径
软件中提供的重建方法有滤波反投影法（Filter Back-Projection，FBP）有序子集最大期望值法 （Ordered Subset Expectation Maximization， OS-EM）。 PET软件中图像重建途径有离线重建和在线重建。 两种途径的区别是：在线重建能够在多床位或动 态扫描时，每一帧投影数据采集完后即重建该帧 图像，而不需要等所有投影数据都采集完后再进 行重建。
重建模块
重建模块参数说明
1. 反投影运算中心相对于检查床的便移量和XY轴旋转角 度 2. 正弦图预处理方法：无，单层重组，傅立叶重组，组 3. 数据校正：散射校正、衰变校正、归一化、几何校正、 角度压缩 4. 重建方法：FBP、OS-EM 5. 图像信息：图像大小、层面选取、放大因子 6. OS-EM选项：迭代次数、子集个数 7. 滤波器选项：R（径向：减小反投影的1/n模糊）；Z （轴向平滑） ；对估测图像做平滑（三维重建时使用）
协议创建窗口
离线重建
离线重建：重建硬盘上已经存在的投影数据。 首先从Operator Prompt模块，选择重建需要的投 影数据和衰减文件。然后添加Reconstruction模块， 设置合适的参数。
扇束重建中OSEM迭代算法与角度采样的模拟研究 OSEM
扇束重建中OSEM迭代算法与角度采样的模拟研究刘力 印胤 吴晓锋（中国科学院高能物理研究所，北京 100039）摘要：本文通过模拟生成扇束断层扫描成像中不同角度采样数目和不同角度采样范围的扇束投影数据，对比研究了不同角度采集条件下图像重建OSEM（有序子集最大似然法）算法的成像效果，验证了在扇束采集与重建中的规律，显示了OSEM迭代重建算法相对于FBP（滤波反投影法）的优越性，并探讨了如何控制角度采样来获取较好的重建效果，同时尽量减少采集与处理的时间。
关键词：CT 图像重建OSEM FBP中图分类号：在CT (Computerized Tomography) 断层扫描与核医学ECT（Emission Computerrized Tomography）断层扫描中，扇束扫描是一种常用采集方式，其重建方法就有解析算法和迭代算法。
有序子集最大期望值方法（Ordered Subsets Expectation Maximization, OSEM）是一种迭代重建算法，它具有较好的重建图像质量和较短的计算时间，已应用在商业ECT的重建和CT图像重建的研究中。
本文通过研究仿真Phantom模型数据分析比较了在不同角度采集条件下的扇束OSEM重建结果，并分别在无噪声和有噪声的情况下同FBP重建图像进行对比，显示了OSEM迭代重建算法的特点和相对于FBP（滤波反投影法）的优越性。
1 CT模型与OSEM重建算法本文模拟了一个扇束CT原型，见图1：图1 模拟CT原型示意图其中被测物体范围（FOV ）是一个以O 为圆心，半径为R 的圆（单层平面），S 为点X 光源，D 为圆弧状等角间距排列的探测器（bin ）阵列，以S 为其圆心；dso 为S-O 距离，max γ为扇束半张角，一个扇束覆盖整个FOV 。
假设整个系统保持物体不动，光源－探测器系统绕O 点进行旋转采样。
o 360 我们选取图像大小为128×128，以相邻两象素间的距离为长度单位，dso=247.3，max γ为，D 中bin 的个数为128。
锥束CT的图像分块OSEM重建算法
锥束CT的图像分块OSEM重建算法
曾理;冀东江;邹晓兵
【期刊名称】《计算机工程与应用》
【年(卷),期】2007(043)032
【摘要】在计算机断层成像(CT)中,有序子集最大期望值算法(OSEM)能够在较短的时间内重建出高质量的图像.对含有噪声的投影数据,投影旋转分度子集划分的不同会影响到图像的重建质量和收敛速度.针对三维锥束CT情况,研究了一种基于图像分块的变子集OSEM重建方法,该方法将图像空间分割成等大小的图像块,然后在迭代过程中,对于不同的图像块用变化的子集进行图像重建.计算机仿真实验表明:该方法在锥束CT图像重建中,能够在抑制噪声的同时提高重建图像的收敛速度.【总页数】4页(P201-203,248)
【作者】曾理;冀东江;邹晓兵
【作者单位】重庆大学,光电技术及系统教育部重点实验室ICT研究中心,重
庆,400044;重庆大学,数理学院,重庆,400044;重庆大学,数理学院,重庆,400044;重庆大学,数理学院,重庆,400044
【正文语种】中文
【中图分类】TP391
【相关文献】
1.基于优化子集顺序的三维OSEM图像重建算法 [J], 孔慧华;潘晋孝;吴琨
2.基于OSEM算法的X射线荧光CT重建算法研究 [J], 刘亚楠
3.OSEM重建算法中几个关键问题的研究 [J], 凌松云;曹文田;包尚联
4.OSEM重建算法及其改进算法的研究和比较 [J], 杨娟;王明泉;石浪;侯慧玲
5.图像重建算法FBP与OSEM在工业CT应用中的对比研究 [J], 贠明凯; 刘力因版权原因，仅展示原文概要，查看原文内容请购买。
医学图像重建算法概述-PPT
反投影 20
20 20
20
1.25 1.00 1.50 1.25
有序子集迭代方法示例
这个图像不知道
5 10 15 20
15 35 知道这些
25 20 30
25
猜测图像
第一子集
1次迭代
10 10
20 0.75
7.5
7.5
0.75
10 10
20 1.75
17.5 17.5
1.75
有序子集迭代方法示例
典型迭代重建示例
这个图像不知道
5 10 15 20
15 35 知道这些
猜测图像
10 10 10 10
25 20 30
25
20 0.75 20 1.75 1次 迭代
9.375 14.0625 21.875 32.8125
20 0.75 20 1.75
20 20 20
20
1.25 1.00 1.50 1.25
比尔定律：描述X射线穿过物体时发生衰减的规律，如下：
I0
均匀物质，衰减系数μ
Id
初始强度
长度L
穿透强度
Id I 0L
决定了是什Байду номын сангаас物质
Id I 0L
实际情况：物质并非均匀，组成复杂，因此，我们将物质分 成许多小份，即有了体素（像素）的概念，如下d：
I0
初始强度
u1 u2 u u4 3
ddd
un
In
45 1 456 2 101 23 4
5
22 67 2 334 1 43 111 12 45 45 88 1 445 3.3 34 134
23 4 23 1
2 345 2.3 111
(医学课件)医学图像重建算法概述
深度神经网络法
利用深度神经网络对医学图像进行自动重建，这种方法需要大量的训练数据和计算资源。
自编码器法
利用自编码器将输入图像进行编码，再解码得到重建的图像，这种方法较为新颖，具有较大的发展潜力。
基于深度学习的图像重建算法
医学图像重建算法的关键技术
03
1
图像数据的预处理技术
2
3
去除无效和冗余的图像数据，提高算法准确性。
医学图像重建算法的实际应用与挑战
04
辅助医生诊断
通过医学图像重建算法，医生可以获得更加清晰、准确的医学图像，提高诊断的准确性和效率。
医学图像重建算法在临床诊断中的实际应用
定制治疗方案
医学图像重建算法可以帮助医生根据患者的个体差异，制定更加精准的治疗方案，提高治疗效果。
评估治疗效果
医学图像重建算法可以对治疗后的效果进行评估，帮助医生判断治疗效果，及时调整治疗方案。
xx年xx月xx日
医学图像重建算法概述
医学图像重建算法简介医学图像重建算法的分类与特点医学图像重建算法的关键技术医学图像重建算法的实际应用与挑战医学图像重建算法的实验分析与比较
contents
目录
医学图像重建算法简介
01
医学图像重建是一种将图像从投影数据或测量数据中恢复或重构的技术，通常在医学诊断和治疗中需要高分辨率和直观的医学图像。
02
模型优化
针对特定问题对模型进行调参，提高模型收敛速度和重建精度。
利用GPU加速技术，提高算法的计算效率。
高性能计算的实现方法
并行计算
将计算任务分配给多个计算节点，加速大规模数据处理。
分布式计算
采用高性能存储设备，提高算法的存储效率和读写速度。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、引言核医学影像设备如单光子断层扫描仪(SinglePositronEmissionComputeTomography，SPECT)、正电子发射断层扫描仪(PositronEmissionTomo-graphy，PET)融合了当今最高层次的核医学技术，是目前医学界公认的极为先进的大型医疗诊断成像设备，在肿瘤学、心血管疾病学和神经系统疾病学研究中，以及新医药学开发研究等领域中已经显示出它卓越的性能。
随着核医学断层影像设备的广泛应用和计算机技术的迅速发展，图像重建方法作为该类设备中的一个关键技术，其研究工作越来越受到人们的重视。
本文概述了传统的图像重建方法，并详细介绍了一种具有较高图像质量和较短计算时间的重建算法—有序子集最大期望值方法(Ord-eredSubsetsExpectationMaximization，OSEM)在核医学影像设备中的应用。
二、传统的图像重建方法在核医学影像设备中，需要根据物体某一层面在不同探测器上检测到的投影值来重建该断层图像层面，即二维图像重建。
传统的图像重建方法主要分为解析法和迭代法。
解析法是以中心切片定理(CentralSliceTheorem)为理论基础的求逆过程。
常用的一种解析法称为滤波反投影法(FilteredBack-Projection，FBP)。
FBP法首先在频率空间对投影数据进行滤波，再将滤波后的投影数据反投影得到重建断层图像。
滤波器选为斜坡函数和某一窗函数的乘积，窗函数用于控制噪声，其形状权衡着统计噪声和空间分辨。
常用的窗函数有Hanning窗，Hamming窗，Butterworth窗以及Shepp-Logan窗。
解析法的优点是速度快，可用于临床实时断层重建。
但当测量噪声较大或采样不充分时，这类算法的成像效果不甚理想，尤其是在核医学断层图像重建中对小尺寸源的成像效果差(即所谓偏体积效应)。
在滤波中如果对高频信号不做抑制，截止频率高，此时空间分辨最好，但所重建的图像不平滑，易产生振荡和高频伪影;反之，采用较低截止频率，过多压抑高频成分的低通窗函数会造成重建图像的模糊，故在变换法中低噪声和高分辨对滤波器的要求是矛盾的，需折衷选择。
且难以在重建中引入各种校正和约束，如衰减校正等。
迭代法是从一个假设的初始图像出发，采用迭代的方法，将理论投影值同实测投影值进行比较，在某种最优化准则指导下寻找最优解。
迭代求解方法的基本过程是: (1)假定一初始图像f(0); (2)计算该图像投影d;
(3)同测量投影值d对比; (4)计算校正系数并更新f值; (5)满足停步规则时，迭代中止;
(6)由新的f作为f(0)从(2)重新开始。
该方法最大优点之一是可以根据具体成像条件引入与空间几何有关的或与测量值大小有关的约束和条件因子，如可进行对空间分辨不均匀性的校正、散射衰减校正、物体几何形状约束、平滑性约束等控制迭代的操作。
其中实现对比的方法有多种，施加校正系数的方法也有多种。
在某些场合下，比如在相对欠采样、低计数的核医学成像中可发挥其高分辨的优势。
但是迭代法收敛速度慢，运算时间长，运算量大，而且重建图像会随着迭代次数的增加而趋于“老化”甚至发散，出现高频伪影，这些缺点极大地限制了它在临床中的应用。
[!--empirenews.page--]三、OSEM迭代算法为了加快收敛速度，减少运算时间，提高图像质量，人们提出了很多快速算法，其中有序子集最大期望值法是很有应用前景的一种快速迭代重建算法，它是在最大似然期望法(MaximumLike-lihoodExpectationmaximization，MLEM)的基础上发展起来的。
MLEM方法旨在寻找与测量的投影数据具有最大似然性(ML)的估计解，其迭代过程是由最大期望值算法(EM)来实现的。
由于是以统计规律为基础，MLEM重建法具有很好的抗噪声能力，是目前公认为最优秀的迭代重建算法之一，尤其是在处理统计性差的数据时，更能显示出它相对于解析法的优越性，但是这种方法仍然存在迭代法的运算量大、运算时间长等缺点。
MLEM方法在每一次迭代过程中，使用所有的投影数据对重建图像每一个象素点的值进行校正，重建图像只被替换一次。
OSEM方法在每一次迭代过程中将投影数据分成N个子集，每一个子集对重建图像各象素点值校正以后，重建图像便被更新一次，所有的子集运算一遍，称为一次迭代过程，它所需要的运算时间与FBP重建的时间基本相等。
在ML-EM方法一次迭代过程中，重建图像被更新一次，而在OSEM方法中重建图像被更新N次，所以OSEM方法具有加快收敛的作用。
OSEM 算法中子集的选取和划分有很多种，在SPECT中投影数据可以根据每个采样角度实时地进行划分和重建，在PET中由于各个探测器上测得的投影数据是在符合判选之后同时获得的，因此可以在全部投影数据采集完成之后划分子集。
不同子集的重建顺序也可以有选择的进行，如
可将两个位于相对垂直的角度上的子集按相邻顺序进行重建，以加快收敛速度。
四、数据模拟与临床实验结果分别采用FBP法、MLEM法和OSEM法对仿真模型和临床数据进行图像重建。
仿真模型类似Jaszczak模型，在64×64的Phantom切片中间的圆形区域上分布着大小不等、呈指数衰减的点状源。
选取观测角度个数为32，探测器单元(Bin)的个数为64，模拟实际投影矩阵，投影数据符合泊松随机分布。
临床PET的Transmission投影数据由美国密西根大学J.Fessler教授提供，观测角度为192个，探测器Bin个数为160，PET为CTIECATEXACT。
图1为采用不同方法对临床(人体模型)投影数据的重建结果，其中FBP法选用的滤波器为But-terworth滤波，陡度因子N=2，截止频率为0.2，OSEM法为N=16一次迭代重建结果;图2为不同子集划分情况下一次迭代重建结果;图3为不同子集划分情况下经过适当迭代次数的重建结果。
由上述Phantom模拟数据与临床投影数据的重建结果可以看出，FBP法重建图像分辨率较低，MLEM法和OSEM法经过多次迭代运算以后可以得到较高分辨率的图像。
比较MLEM法和OSEM法，两者一次迭代运算的时间相同，且子集划分数目N与迭代次数乘积相同时，两者重建图像质量相近，相当于运算速度提高了N倍。
对不同子集划分的研究结果表明，由于OSEM的一次迭代定义为所有的子集通过一次投影校正，所以无论将投影数据划分为多少个子集，其一次迭代运算时间是相同的，所不同的是迭代过程中的收敛速度。
子集划分个数增加，迭代收敛速度增加，迭代次数减少，但子集划分个数并不是越多越好，它和图像重建质量之间存在均衡关系，在临床应用中必须选取合适的子集个数，这样才能在提高运算速度的同时，确保重建图像质量。
[!--empirenews.page--]目前在西门子的HR+等PET的升级软件中都含有OSEM重建算法，由于OSEM算法在保证成像质量的同时可以大大减少运算时间，在核医学影像设备的图像重建中具有广阔的应用前景。
恰当地选取子集个数，可用较少的OSEM迭代运算获得图像质量好于FBP的重建结果，所花费的运算时间并不比FBP法多出许多。
此外，由于迭代法可以合并一些校正计算，使得选用迭代法重建更有效、更合算。
因此，OSEM迭代将代替FBP法而成为核医学断层影像设备的主流重建算法。