卷积反投影重建(二维)
反投影重建算法
反投影重建算法
反投影重建算法（FBP）是一种计算机断层扫描成像（CT）重建图像
的方法。
该算法基于通过旋转体与X射线束的物理原理，将多个方向
的X射线透射数据进行积分，并使用反投影算法将数据重构成一张图像。
FBP算法分为两个基本部分：投影操作和反投影操作。
投影操作是一
种从图像中提取片段的技术，而反投影操作则是将这些片段重构成图像。
FBP重建算法的本质是一种频域过滤操作，其通过滤波技术提取
图像中的高频信息，并使用反投影技术将其还原为一张二维图像。
反投影重建算法的主要优点是其速度和适应性。
这种算法能够轻松地
生成高质量的图像，并且对于许多不同的应用程序都可以使用不同的
滤波模式。
目前，FBP算法被广泛应用于医学成像、工业检测和材料
科学等领域。
需要注意的是，FBP算法并不是完美的。
由于其基于体积的积分，因
此它可能受到一个“锐角偏差”问题的影响。
锐角偏差问题是指，当
图像中存在锐利的边缘或角落时，算法可能会出现伪影或失真的问题。
为了应对这个问题，一些改进算法被提出，例如金刚簇算法（来自中
国科技大学），基于块的迭代顺序最小化算法和模糊模式匹配算法等。
总之，反投影重建算法是一种实用的成像算法，对于许多不同的领域都具有广泛的适用性。
虽然这种算法具有其局限性，但是通过改进算法可以进一步提高它的可靠性和精度。
卷积和反卷积的计算公式
卷积和反卷积的计算公式一、卷积计算公式。
（一）离散卷积（一维情况）设离散序列x[n]和h[n]，它们的卷积y[n]定义为：y[n]=∑_m =-∞^∞x[m]h[n - m]（二）离散卷积（二维情况）对于二维离散信号x[m,n]和h[m,n]，其卷积y[m,n]为：y[m,n]=∑_k =-∞^∞∑_l=-∞^∞x[k,l]h[m - k,n - l]（三）连续卷积（一维情况）对于连续函数x(t)和h(t)，它们的卷积y(t)定义为：y(t)=∫_-∞^∞x(τ)h(t-τ)dτ二、反卷积计算公式。
反卷积（也称为去卷积）是卷积的逆运算。
在离散情况下，如果已知y[n]（卷积结果）和h[n]，求x[n]，可以通过求解以下方程（在某些条件下）：y[n]=∑_m =-∞^∞x[m]h[n - m]1. 频域方法（离散情况）- 对y[n]、h[n]分别进行离散傅里叶变换（DFT），得到Y[k]和H[k]。
- 根据卷积定理Y[k]=X[k]H[k]，则X[k]=(Y[k])/(H[k])（假设H[k]≠0）。
- 再对X[k]进行逆离散傅里叶变换（IDFT）得到x[n]。
2. 迭代算法（离散情况）- 一种简单的迭代算法是假设初始的x^0[n]=y[n]/h[0]（当h[0]≠0时）。
- 然后通过迭代公式x^i + 1[n]=x^i[n]+frac{y[n]-∑_m =-∞^∞x^i[m]h[n - m]}{∑_m =-∞^∞h[m]h[n - m]}逐步逼近真实的x[n]，其中i表示迭代次数。
在连续情况下，反卷积的求解更加复杂，通常也可以利用频域方法，通过傅里叶变换将问题转换到频域，利用Y(ω)=X(ω)H(ω)，得到X(ω)=(Y(ω))/(H(ω))（假设H(ω)≠0），再通过逆傅里叶变换得到x(t)，但在实际应用中要考虑到函数的性质、收敛性等诸多问题。
计算机断层扫描成像(CT)
1引言自七十年代初第一台电子计算机断层扫描装置问世以来，成像技术发展异常迅速，设备不断更新。
以医学成像为例，已实现了三大飞跃，即脏器清晰图像的获得，把生化病理研究推向分子结构的水平和直接提供有关成像组织的化学成分的信息，步入了断层显像的新时代。
计算机断层扫描和图像重建技术，是在不破坏物体情况下，将物体每一个断层面上的结构和组份的分布情况显示出来的一种实验方法，都是利用计算机图像重建的方法来得到物体内部的信息。
人们对射线成像的最早认识是从x 光机开始的。
医用x 光机成像技术的发展和应用已有近百年的历史，它是利用x 射线的物理性能和生物效应，来对人体器官组织进行检查。
由于普通x 光机只能把人体内部形态投影在二维平面上，因此会引起成像器官和骨骼等的前后重叠，造成影像模糊。
为了克服这一缺点，英国ENI 公司的工程师豪恩斯菲尔德(G.N.Hounsfield)运用了美国物理学家科马克(Cormack)于1963年发表的图像重建数学模型，推出了第一台x 射线计算机断层图像重建技术（X-CT ）装置，并1977年9月在英国Ackinson Morleg 医院投入运行。
1979年该技术的发明者Hounsfield 和Cormack 为此获得了诺贝尔医学奖。
X-CT 的出现是X 射线成像技术的一个重大突破。
经过多代的发展，X-CT 已获得广泛的应用。
在医学上，目前已可用来诊断脊柱和头部损伤，颅内肿病，脑中血凝块，及肌体软组织损伤，胃肠疾病，腰部和骨盆恶性病变等等。
目前X-CT 除了广泛应用于临床诊断、生命科学和材料科学以外，还在工业和交通等方面也有重要的应用，例如，在线实时无损检测工业CT 等。
2CT 成像实验原理2.1概述数学上可以证明，通过对物体进行多次投影就可得到该物体的几何形状。
CT 的基本思想是：让一束γ射线投射在物体上，通过物体对γ射线的吸收（多次投影）便可获得物体内部的物质分布信息。
当强度为0I 的一个窄束γ射线穿过吸收系数为μ的物体时，其强度满足指数衰减关系0ut I I e -=(1)式中t 为射线所穿过物质层厚度。
反投影重建算法
反投影重建算法反投影重建算法是一种常用于医学影像领域的图像重建方法，其原理基于X射线的穿透性质和成像原理。
该算法可以从一系列投影图像中恢复出被测物体的内部结构，为医学诊断提供重要的信息。
在了解反投影重建算法之前，我们先来了解一下X射线成像的基本原理。
X射线成像是一种无创检查方法，通过将X射线束通过被测物体后，测量X射线在不同方向上的衰减情况，从而获取被测物体的内部结构信息。
在传统的X射线成像中，使用的是传统的透射成像方法，即将X射线通过被测物体后，直接在胶片或探测器上记录吸收信息。
然而，传统透射成像存在一些局限性，如对于密度相近的组织无法进行有效区分。
反投影重建算法的出现就是为了解决传统透射成像的局限性。
该算法利用了被测物体对X射线的吸收特性，通过在不同方向上测量X 射线的衰减情况，然后将这些投影数据进行反投影操作，即根据投影数据在不同方向上的衰减信息，将其反向投射回原来的空间位置。
这样，就可以通过反投影操作恢复出被测物体的内部结构。
具体来说，反投影重建算法的步骤如下：1. 收集投影数据：通过将X射线束从不同方向照射到被测物体上，并在胶片或探测器上记录下不同方向上的吸收信息，形成一系列投影数据。
2. 反投影操作：将投影数据进行反投影操作，即将各个投影数据反向投射回原来的空间位置。
这个过程可以看作是将每个投影数据沿着其对应的投影方向，从胶片或探测器上反向传播到被测物体内部的过程。
3. 累加重建：将所有反投影后的数据进行累加，得到最终的重建图像。
在反投影操作中，相交的区域会累加，而非相交的区域则不会累加，从而在最终的重建图像中得到被测物体的内部结构信息。
反投影重建算法是一种简单而有效的图像重建方法，在医学影像领域得到了广泛的应用。
与传统的透射成像相比，反投影重建算法可以提供更多的信息，能够更好地区分密度相近的组织，从而提高了诊断的准确性。
此外，反投影重建算法还可以应用于其他领域，如安全检查、材料检测等。
儿童CT检查操作的特殊要求
儿童CT检查操作的特殊要求1 正当性判断——儿童CT检查的正当性判断应比成人更加谨慎，临床适应症更严格，逐例评估收益-风险比。
——在严重外伤等紧急情况下，对儿童实施CT检查具有正当性。
在条件许可的情况下，儿童的影像检查应优先考虑使用无电离辐射的超声或磁共振成像等方法。
——儿童心脏大血管CT成像较心导管造影更具正当性。
——不使用CT对儿童的疾病进展进行连续监测或健康筛查。
2 扫描前准备——详细了解患儿病史，获取既往检查的影像资料。
——婴儿检查前采取喂食与襁褓束缚的方式来减少运动。
提供儿童主题装饰、音乐、玩具、数数以及与父母保持语音通话的方式，分散幼儿注意力，安抚恐慌焦虑情绪，取得患儿的配合。
——必要时予以镇静或麻醉。
设备条件允许时采用双源CT大螺距扫描或超宽探测器采集，以及尽可能进行不注射对比剂的平扫，可减少对镇静或麻醉的需求。
——远离扫描部位的辐射敏感器官采用铅制围裙、眼镜、围脖等防护用具遮盖，尽可能予以包裹。
扫描区域内不使用铋制防护用品。
——备齐儿童不同年龄段适用的复苏囊、喉镜等抢救用物型号。
3 扫描方案优化——为不同分组的儿童定制各检查部位的儿童专用扫描协议，最合理的分组依据为受检部位体厚，其次为患儿体重，前两者信息不易获得时根据患儿年龄分组。
——根据临床指征或检查目的进行个性化的扫描参数优化，在满足诊断需求前提下尽可能降低患儿辐射剂量。
——为获得一致的图像质量，不同厂家甚至同一厂家不同型号的CT应定制不同的扫描协议。
3.1 扫描体位——使用安全绑带将儿童固定于检查床。
儿童胸腹部检查两臂无法上举抱头时，伸直上肢置于身体两侧。
——常规仰卧位，根据检查及静脉穿刺部位选择头先进或足先进。
不能配合呼吸的患儿进行肺部双期相扫描时，取侧卧位，贴紧检查床的一侧表现为呼气相，另一侧为吸气相。
——检查部位置于机架旋转中心。
3.2 定位扫描——80kV～100kV，正位定位扫描时，球管置于检查床下。
定位扫描范围限制目标检查部位。
MRI设备详细介绍
mri设备详细介绍mri, 设备MRI设备MRI设备是利用生物体的磁性核（主要是氢核）在磁场中所表现出的MR特性来进行成像的设备。
随着超导技术、磁体技术、电子技术、计算机技术和材料科学的进步，MRI设备得到飞速的发展。
MRI设备已成为最先进、最昂贵的现代化诊断设备之一。
MRI设备既是评价医院综合能力的一项重要指标，又是医院现代化程度和诊断水平的标志。
我国现有600多台MRI设备正在运行，并以每年几十台的速度增长（含临床应用型和临床研究型）。
本章将以临床应用型永磁开放式MRI设备为例，系统地介绍MRI设备的构成和工作原理。
第一节概述一、发展简史MR现象于1946年第一次由布洛赫（F.Bloch）领导的斯坦福大学研究小组和伯塞尔（E.Purcell）领导的哈佛大学研究小组分别在水与石蜡中独立地观察到。
因此，布洛赫和伯塞尔共同获得了1952年的诺贝尔物理学奖。
随后，人们利用MRI技术进行了多领域的应用。
MRI设备早期集中在物理和化学方面，用来确定化学成分、分子结构和反应过程。
1967年，第一次用MRI设备测试人体活体。
1971年，达马丁（Damadian）发现了MRI的一个重要参数—T1。
肿瘤组织的T1值远大于相应正常组织的T1值。
此结果预示着MRI设备在医学诊断中的广阔应用前景。
1973年，受CT图像重建的启示，纽约州立大学的劳特布尔（Lauterbur）在《Nature》杂志上发表了MRI设备空间定位方法（均匀静磁场上迭加梯度磁场）。
利用MRI模型（两个并排在一起的充水试管）的四个一维投影，成功的获得了第一幅MRI模型的二维图像。
1974年，曼斯菲尔德（Mansfield）研究出脉冲梯度法选择成像断层的方法。
1975年，恩斯特（Ernst）研究出相位编码的成像方法。
1977年，爱特斯坦（Edelstein）、赫切逊（Hutchison）等研究出自旋扭曲（Spin Warp）成像法。
1977年，达马丁完成了首例动物活体肿瘤检测成像，并获得首张人体活体MRI设备图像。
CT系统参数标定模板的设计
CT系统参数标定模板的设计饶勤;张明辉【摘要】本文针对二维CT系统的参数标定问题,在医学数据的基础上进行探究,首先采用普通模板对未知CT系统进行标定,确定了该CT系统旋转中心、探测器单元之间的距离以及X射线的方向;然后基于反投影滤波模型,在MATLAB软件上编写算法实现了CT成像,并运用了更加高效方便的卷积函数代替了滤波过程,以及采用汉宁窗算法减少了傅里叶变换的数据损失;最后分析整个标定与重建模型的精度和稳定性,得出误差主要源于数据的离散性的结论,并进一步总结出标定模板的设计要点,提出了一个高效、易用的标定模板.该模板经过验证,其吸收率函数具有易于拟合的特征.新的模版与标定算法有效提升了标定精度,而且对噪声有着更高的容忍性.【期刊名称】《电声技术》【年(卷),期】2019(043)006【总页数】4页(P18-21)【关键词】CT系统参数标定;反投影滤波模型;汉宁窗算法;MATLAB【作者】饶勤;张明辉【作者单位】东南大学自动化学院,江苏南京211189;东南大学自动化学院,江苏南京211189【正文语种】中文【中图分类】TH7891 引言CT(电子计算机断层扫描术)是一种利用样品对射线能量的吸收特性对生物组织和工程材料的样品进行断层成像，由此获取样品内部的结构信息的扫描技术。
因其检测技术对样品无损伤，现已广泛应用于医疗诊断、工业检测、安检等诸多领域当中[1-2]。
在CT系统中，采用以平行X射线垂直入射探测器平面的方式的二维CT系统最为典型，其把等间距排列的探测器单元近似看作为一系列的接收点，探测器与X射线发射器的相对位置保持不变的CT系统。
在检测过程中，其绕旋转中心逆时针旋转180次，通过512个探测器单元来测量样品吸收衰减后的X射线能量，然后通过增益等操作得到180组接收信息。
而在CT系统的安装过程中往往存在着一定的误差，从而使系统成像质量较差。
由此，对安装完成的CT系统进行参数定标就显得十分重要。
反投影重建算法
反投影重建算法介绍反投影重建算法是一种常用于医学图像重建的技术，它通过将多个投影图像反向投影到空间中，来重建出原始物体的三维形态。
该算法基于X射线摄影和计算机图形学的理论，被广泛应用于医学成像领域，尤其是计算机断层扫描(CT)。
基本原理反投影重建算法的基本原理是通过测量某物体在不同角度下的投影数据，然后将这些数据反向投影到空间中，得到物体的三维重建图像。
具体步骤如下：1.采集投影数据：使用X射线或其他成像设备，获得物体在不同角度下的投影数据。
这些投影数据可以通过测量X射线的衰减程度来获取。
2.投影数据的逆变换：将投影数据的方向进行反转，得到反投影数据。
这一步是算法中最重要的一步，它将每个测量点的投影数据映射回三维空间。
3.反投影重建：将反投影数据从二维空间转换为三维空间。
这可以通过将反投影数据相加来实现。
在这个过程中，需要注意对数据进行插值，以确保重建图像的精确度。
算法细节反投影重建算法的具体细节有多种实现方式，下面列举几种常见的算法：1. 直接反投影算法直接反投影算法是一种最简单的算法，它将每个投影像素直接投影回三维空间。
具体步骤如下：1.对每个投影像素点，计算它在三维空间中的坐标。
这可以通过测量点在投影平面上的位置和投影矩阵来实现。
2.将每个投影像素点的坐标值累加到三维空间的相应位置上。
这一步是反投影的核心过程，通过将每个像素点的坐标值相加，最终可以得到物体的三维重建图像。
3.可选的后处理步骤：根据需要，可以对重建图像进行滤波、增强或其他处理，以提高图像的质量。
2. 迭代反投影算法迭代反投影算法是一种更复杂但更精确的算法，它通过多次迭代求解，逐渐优化重建图像的质量。
具体步骤如下：1.初始化重建图像：将重建图像的像素初始化为某一固定值。
2.迭代求解：重复以下步骤，直到算法收敛或达到最大迭代次数：–计算投影数据的理论值：根据当前重建图像，计算预测的投影数据。
–计算误差：比较预测的投影数据与实际测量数据之间的差异，得到误差。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
卷积反投影图像重建1 反投影重建基本介绍[1]设待重建图像为),(y x a ，它的二维傅氏变换为^12(,)(,)A A ωωρφ=。
根据中心切片定理，^(,)A ρφ可通过),(y x a 在不同视角φ下的投影()rp x φ的一维傅氏变换求得。
即：待建图像：12^1212()12122^2cos()02cos()2cos()(,)(,)(,)1(,)4(,)(,)(,)i x y i r i r i r a r a x y F A A e d d A e d d P e d d d P e d ωωππρθφππρθφππρθφθωωωωωωπρφρρφρφρρφφρρφρ-∞∞+-∞-∞∞--∞∞--∞∞--∞======⎰⎰⎰⎰⎰⎰⎰⎰[](1.1)因为cos()r x r θφ=-，所以有：122(cos sin )22cos()r x y x y x r ωωπρφφπρπρθφ+=+==- 同时：12d d J d d ωωρφ=11222//2cos 2sin 4//2sin 2cos J ωρωφπφπρφπρωρωφπφπρφ∂∂∂∂-===∂∂∂∂ (1.2)先来看该式的第二个积分：[]22cos()cos()cos()cos()(,)(,)|()(,)|(,)|cos(),r r r r i x i r x r r r x r r x r P ed Pe d h x p x g x g r πρπρθφθφθφθφρρφρρρφρφφθφφ∞∞-=--∞-∞=-=-==*==-⎰⎰ (1.3)式中：(,)()(,)r r r g x h x p x φφ=*(1.4)^121(,)(,)()()(,)rA A F p x PP φφωωρφρρφ⎡⎤====⎣⎦式(3.10)的物理意义是投影(,)r p x φ经过传递函数为1[()]rF h x ρ=的滤波器后得到的修正后的投影(,)r g x φ在满足cos()r x r θφ=-时的值。
将(3.11)代入(3.8)，得到：^(,)[cos(),]a r g r d πθθφφφ=-⎰ (1.5)称为滤波反投影方程，其物理意义是经过给定点(,)r θ的所有滤波后的投影在0φ=~π范围内的累加—反投影重建，得出(,)r θ点的像素值。
可见，滤波（卷积）反投影算法的具体包含三大步：(1) 把在固定视角i φ下测得的投影(,)r p x φ经过滤波，得到滤波后的投影(,)r g x φ；(2) 对每一个i φ，把(,)r i g x φ反投影于满足cos()r i x r θφ=-的射线上的所有各点(,)r θ；(3) 将步骤(2)中的反投影值对所有0<φπ≤进行累加（积分），得到重建后的图像。
2 重建流程2.1 首先我们利用phantom （）函数产生一个头部幻影图像，用以检测二维重建算法，代码如下：I=phantom(256); subplot(2,2,1) imshow(I,[]);title('256*256原始图像'); 效果图如图1所示，图1为一个大椭圆和几个小椭圆。
2.2 初始参数设置重建采用的是平移加旋转的扫描方式，射线源在某一角度下水平移动，将物体全部照射后旋转一角度，如此重复，在这个过程中探测器相应地运动以接收X 射线。
根据此原理,将重建程序的初始参数设置如下:[N,N]=size(I);z=2*ceil(norm(size(I)-floor((size(I)-1)/2)-1))+3;% radon变换默认平移点数/角度Nt=360; % 角度采样点数Nd=N; % 平移数x=pi/180; % 角度增量d=N/Nd; % 平移步长theta = 1:Nt;a=zeros(N);2.3 产生无噪声投影数据[R,xp] = radon(I,theta);e=floor((z-Nd)/2)+2;R=R(e:(Nd+e-1),:);R1=reshape(R,256,360);radon(I,theta)产生I投影,默认z点/角度,即使指定N点也是z 点.所以为避免重建图像放大或缩小,下面计算取投影时需补偿,补偿量e,如对256的图像,补偿为55,即pm的第55个点作为计算用的第一个投影。
2.4 添加噪声并将所有噪声平行投影进行显示[mm,nn]=size(R1);di=lognrnd(0,0.15,mm,nn);R1= 10*(R1-min(R1(:)))/( max(R1(:))-min(R1(:)));I0 = 1.5e5; % incident photons; decrease this for simulating "low dose" scansrand('state', 0), randn('state', 0);yi= poissrnd(I0 * di.*exp(-R1))+3*randn(size(R1));if any(yi(:) == 0)warn('%d of %d values are 0 in sinogram!', ...sum(yi(:)==0), length(yi(:)));endR1 = log(I0 ./ max(yi,0.01)); % noisy sinogramR1=max(R1,0);% 显示ff=2;uu=22000;v=ff*exp(R1/uu);subplot(2,2,2)imagesc(R1);title('256*360有噪声平行投影');colormap(gray)colorbarQ=reshape(R1,256,360);效果图如图-2: 图-2 2.5 滤波器的选择与设计最基本的从投影重建图像的滤波器：1971年提出的R-L重建滤波器(下图中实线表示)。
图-3 R-L 滤波器示意图空域表达式为:(2.1)其中B为截至频率。
若，d为空间采样间隔，可解出离散的滤波器空间脉冲响应：(2.2)MATLAB 代码如下:g=-(Nd/2-1):(Nd/2);for i=1:256if g(i)==0hl(i)=1/(4*d^2);else if mod(g(i),2)==0hl(i)=0;elsehl(i)=(-1)/(pi^2*d^2*(g(i)^2));endendendk=Nd/2:(3*Nd/2-1);2.6 重建图像代码设计:for m=1:Ntpm=Q(:,m);u=conv(hl,pm);pm=u(k);Cm=((N-1)/2)*(1-cos((m-1)*x)-sin((m-1)*x));for i=1:Nfor j=1:NXrm=Cm+(j-1)*cos((m-1)*x)+(i-1)*sin((m-1)*x);if Xrm<1n=1;t=abs(Xrm)-floor(abs(Xrm));elsen=floor(Xrm);t=Xrm-floor(Xrm);endif n>(Nd-1)n=Nd-1;endp=(1-t)*pm(n)+t*pm(n+1);a(N+1-i,j)=a(N+1-i,j)+p;endendend重建后的图像如图-3 所示: 图-3 程序中还包含对重建结果的评价/归一化均方距离判据/归一化平均绝对距离判据以及程序的运行时间等.不再进行详细介绍。
程序的最终运行效果如下：图-4 程序最终运行效果图3、分析反投影重建方法包括卷积反投影重建的缺点是会产生星状伪迹，原因分析如下：断层平面中某一点的密度值可以看作是这一平面内所有经过该点的射线的投影值之和（的均值）。
整幅重建图像可以看作是所有方向下的投影累加而成。
射线标号示于图5中，像素值（代表密度）分别1x ,2x ,3x ,4x ， 赋值如下:15x =,20x =,32x =,418x =图5 图像重建像素值根据投影的定义(某条射线投影值为该条射线穿过的所有的像素值之和)，每条射线的投影i p (1,2i =)为：1215p x x =+=, 23420p x x =+=,3137p x x =+= 42418p x x =+=, 532p x ==, 61423p x x =+= 720p x ==根据反投影重建算法的物理意义，重建图像中各像素，得到： 113635x p p p =++=，214723x p p p =++=， 323529x p p p =++=424661x p p p =++=，(a) 原图像像素值 (b)反投影重建后图像 (c)求平均后图像图 6 反投影示例重建后的图像如图6(b)所示，可以看出原图像中像素值不为零的点反投影重建后仍较突出，但原图中像素值为零的点，经反投影重建后不再为零，即有伪迹。
有时为了使重建后图像的像素值更接近于原图的像素值，在求反投影时，把数据除以投影的数目（即射线数），如图6(c)所示。
因此有：,11pn k k ii px pn ==∑(3.1)计算机图像处理技术结课（报告）共9页第11页该式可作为反投影重建算法的计算式。
其中k x 表示像素k 的值，i k p ,表示经过像素k 的第i 条射线投影，pn 表示图像内的射线条数。
图7(a)表示空间中一个孤立点源A ，密度为1。
经过A 点的三条射线也示于图中。
射线束理论上可以很多，取三条示意。
不经过A 点的射线投影为零，经过A 点的射线投影值均为1， 1231p p p ===。
(a) 孤立点源A 及三条射线 (b)相应的反投影重建图，有星状伪迹图 7 孤立点源的反投影重建及星状伪迹经反投影重建后，得到A 点的像素值为123()/31A f p p p =++=。
A 点以外的像素值原来为零，经反投影重建后不等于零，而是等于1/3。
所以，经反投影重建后的图像除保留A 点的像外，还有像素值为1/n 的灰雾背景，后者称为星状伪迹。
产生星状伪迹的原因在于：反投影的本质是把取自有限物体空间的投影均匀地回抹（反投影）到射线所及的无限空间的各点上，包括原先像素值为零的点。
参考文献[1] 滤波反投影,Zhengdong Zhou.2012(1)(2)(3)。