# XCT 检测：原理、重建与三维数据分析 版本：2026-10-28 三维图像的每个体素都经过采集与计算。沿着投影、重建、分割、统计这条顺序，学会判断图像是否可靠，再把它变成有单位、可复查的结果。 这是原创课程阅读稿。自编案例和数值为教学设定；公式的假设写在正文中。网页交互不包含在 Markdown 文件内。部分原文需机构权限，标准以项目采用版次为准。 ## 目录 - 01 · XCT 得到了什么：投影、切片和体数据 - 02 · X 射线衰减：从 Beer–Lambert 定律到投影值 - 03 · 几何放大、体素与分辨率 - 04 · 扫描设计：曝光、投影数、转动与校准 - 05 · 重建原理：Radon 变换、FBP、FDK 与迭代法 - 06 · 伪影排查：束硬化、环、条纹、运动和截断 - 07 · 数据导入：位深、字节序、切片顺序和空间单位 - 08 · 孔隙分割：材料 ROI、阈值与 Otsu - 09 · 连通域、去噪和形态学：孔的数量为何会变 - 10 · 缺陷表怎么读：体积、等效径、形状与位置 - 11 · 孔隙率与分布：统计总体、代表体积和不确定度 - 12 · 裂纹、壁厚与尺寸：为什么不能套用孔分析全部步骤 - 13 · 完整手算：从 100 万体素得到孔隙报告 - 14 · 软件实操路线：导入、分割、对象表与导出 - 15 · 原位 XCT 与 DVC：从多次扫描研究损伤演化 - 16 · 文献导读与最终报告：重建、分割、检出各自验证 ## 01 · XCT 得到了什么：投影、切片和体数据 先区分 X 射线投影与三维重建，明确检测对象、数据层级和能够回答的问题。 ### 一张透射图不是一个截面 X 射线穿过物体后，各路径上的衰减叠加在探测器上。普通投影图把厚度方向结构重叠起来；CT 从多个角度测量，再通过重建估计三维衰减分布。看到二维探测器图，并不意味着直接看到了其中某一层。 XCT、micro-CT 和工业 CT 的叫法反映应用或尺度，并不自动保证某种分辨率。应说明设备几何、能谱、焦点、探测器和实际成像能力。同步辐射、实验室锥束和高能工业系统的条件可以差别很大。 ### 保留三个数据层 第一层是投影及暗场、平场、几何和采集记录；第二层是带空间标定的重建灰度体；第三层是分割标签、表面网格、缺陷表和统计报告。第三层是经算法解释后的结果，不能替代前两层。颜色鲜艳的三维渲染只是显示方式。 灰度值通常与特定系统和条件下的有效衰减有关。未经物理标定，不应把工业 CT 灰度直接标成密度、元素浓度或医学 HU。多种材料可能灰度重叠，同一种材料也会因伪影产生灰度梯度。 ### 从问题倒推扫描 若关心孔隙率，需要有代表性的体积与明确外包络；若关心小裂纹，则需足够的裂纹张开衬度和合适取向；若关心尺寸，要做几何与尺度校准。一个设置不能同时无限提高视场、空间分辨率和扫描速度。 研究前列出目标量、希望识别的最小特征、材料与最大穿透厚度、允许扫描时间和独立验证方法。先做合理的试扫检查最厚方向是否仍有有效信号，再决定正式参数。 ### 初学者阅读路线 Kak 与 Slaney 的教材解释投影为何可以重建；Withers 等的综述连接实验、分析和局限；ASTM E1441 提供 CT 的术语与原则入口。先理解数据形成，再打开分割软件，较容易发现导入与阈值错误。 本专题最后的孔隙算例和下载练习使用合成体素。它们有已知结构，适合学习算法和单位；不代表真实设备的噪声、检出能力或某一材料的质量水平。 ### 本章来源 - [Kak 与 Slaney：Principles of Computerized Tomographic Imaging](https://www.slaney.org/pct/pct-toc.html)：作者提供的断层成像教材阅读入口。 - [Withers 等（2021）：X-ray computed tomography](https://doi.org/10.1038/s43586-021-00015-4)：XCT 实验、重建、分析与可重复性的系统综述。 - [ASTM E1441](https://store.astm.org/standards/e1441)：CT 原理及术语指南 ## 02 · X 射线衰减：从 Beer–Lambert 定律到投影值 逐步解释透过强度、线积分、暗平场校正和多能谱偏离，避免把黑白图直接当作材料标签。 ### 单能窄束的出发点 理想单能窄束中，均匀材料满足 I=I₀exp(−μL)，μ 为线性衰减系数，L 为路径长度。非均匀路径上，p=−ln(I/I₀)=∫μ(s)ds。μ 的单位是长度的倒数，p 无量纲。射线经过更长路径或更大 μ 时，透过强度通常更低。 衰减涉及吸收和散射等相互作用，依赖光子能量和材料组成。质量衰减系数 μ/ρ 与线性系数不同，需乘密度 ρ 才能换算。不能仅按密度给所有材料排序，也不能把一个能量点的系数当作整个 X 射线管谱。 ### 原始探测器信号怎么变成 p 实际记录有暗信号 D 与空场/平场 F。简化归一化为 T=(I−D)/(F−D)，随后 p=−ln T。暗平场应对应合适的曝光和设备状态，并处理坏点与漂移。若 F−D≤0 或 I−D≤0，说明需要检查饱和、缺光、坏点或校正，不能直接取对数。 教学设定 D=100、F=1100、I=600，则 T=0.5，p≈0.6931。若已知均匀单能路径 L=5 mm，则 μ≈0.1386 mm⁻¹。这里忽略散射、能谱和探测器非线性，实际系统不能只凭这一式完成全部标定。 ### 多能谱为何更复杂 X 射线管产生宽能谱，检测强度可写成不同能量贡献的积分：I=∫S(E)R(E)exp[−∫μ(E,s)ds]dE。低能光子往往更容易被吸收，穿透后谱变硬，因此 −ln(I/I₀) 不再严格对应一个固定 μ 的线积分。 滤片可减少低能成分，但会降低总通量。提高管电压、管电流或曝光时间的作用不同，还受焦点尺寸、功率与探测器限制。参数选择应依据实际材料和厚度验证，不应提供一个所有金属通用的 kV 值。 ### X 射线管与光子相互作用 实验室射线管中，电子经电压加速后撞击靶材，减速产生连续制动辐射，内层电子跃迁产生特征谱线。以 keV 表示的最大光子能量数值由以 kV 表示的管电压上限决定，但大多数光子不在该最大能量处。管电流主要影响光子通量，功率和靶热负荷又会限制焦点及工作条件。 在常见工业 CT 能量与材料范围内，光电吸收和 Compton 散射等过程影响衰减与对比。光电响应对元素和能量较敏感，散射还受电子密度等因素影响；具体比例应查能量与材料对应的数据。NIST XCOM 可以查询各项系数，使用时核对质量系数/线性系数和单位。 ### 探测器把光子变成灰度 常见平板探测器可通过闪烁体先将 X 射线转换为可见光，再由光电传感器读出；其他系统也可能采用直接转换或计数方式。闪烁光扩散、像素孔径、电子噪声、响应不均与饱和都会影响图像。探测器像素很小并不保证系统点扩散同样小。 量子噪声在理想独立计数模型中标准差约为 √N，故相对噪声约为 1/√N；实际读出与散射等噪声还需另计。把曝光加四倍可在该理想模型下把相对量子噪声减半，但也增加时间和热漂移影响，不能保证最终计量误差减半。 ### 吸收对比与相位对比 吸收对比主要利用强度衰减差。X 射线相位也会因折射率变化而改变，在适当相干性与传播距离下，波传播可把相位变化转成边缘强度变化，形成相位对比。低吸收差材料可能因此更清楚，但边缘增强不能直接按纯吸收灰度解释。 相位恢复需要相应模型与假设，可能影响边界锐度和定量。实验室及同步辐射系统是否采用相位成像，应在采集和重建记录中写明。不能把一个有明暗边缘的孔仅凭纯衰减阈值处理后，宣称已得到无偏的开口宽度。 ### 从公式反查图像 若同一种材料中心与边缘灰度不同，先检查束硬化、散射和重建，而不是立刻宣布存在成分梯度。若厚方向投影接近暗场，重建条纹可能来自光子不足；后期滤波不能恢复已经缺失的信息。 练习：透过率从 0.5 降至 0.25，p 从 ln2 变为 ln4，增加 ln2。在上述单能均匀模型中，这对应路径长度加倍；在多能谱系统中，不能保证仍严格成比例。 ### 本章来源 - [Kak 与 Slaney：Principles of Computerized Tomographic Imaging](https://www.slaney.org/pct/pct-toc.html)：作者提供的断层成像教材阅读入口。 - [Withers 等（2021）：X-ray computed tomography](https://doi.org/10.1038/s43586-021-00015-4)：XCT 实验、重建、分析与可重复性的系统综述。 - [NIST XCOM：Photon Cross Sections Database](https://physics.nist.gov/PhysRefData/Xcom/Text/chap4.html)：元素、化合物和混合物的光子相互作用与衰减数据。 ## 03 · 几何放大、体素与分辨率 用源—物—探测器关系理解体素采样，讲清焦点模糊、运动、MTF 与检出能力。 ### 几何放大如何计算 点源近似下，放大倍数 M=SDD/SOD，其中 SDD 是源到探测器距离，SOD 是源到物体参考平面的距离。探测器像素间距 pdet 对应物方采样约 pdet/M。锥束物体不同深度的实际投影放大略有差别，重建依赖完整几何校准。 教学设定 SDD=600 mm、SOD=100 mm、pdet=120 μm，则 M=6，物方采样约 20 μm。把物体移近源可以增加放大，但可能缩小可覆盖视场、增加碰撞风险，并使有限焦点模糊更重要。 ### 体素尺寸为什么不是实际分辨率 有限焦点尺寸 f 在物方造成的几何模糊量级约为 f(M−1)/M；探测器模糊、机械漂移、运动、重建滤波和信噪比还会共同作用。重建软件设置 5 μm 体素，不会自动让系统分辨两个相隔 5 μm 的结构。 空间分辨率可通过适用试样的边缘、线对或 MTF 等方法评价，检出还同时取决于目标衬度、形状和噪声。ASTM E1695 的性能测量将分辨能力与对比敏感性联系起来，不能只报告体素间距。 ### 部分体积效应 一个体素可能同时包含金属和空气，重建值会体现模糊与混合。细裂纹即使开口小于一个体素，也可能造成弱灰度变化，但这不意味着可以准确测量它的宽度。相反，一个几体素大小的低衬度孔也可能不可靠。 不存在对所有设备与材料都成立的“三个体素必检出”规则。应区分可见、可检出、可分割和可准确测量，并用已知特征或独立方法评估目标任务。 ### 如何设计代表性试验 尺寸大时，全件扫描提供覆盖，小 ROI 高倍率扫描提供细节，两者可互补。ROI 的取样位置应基于研究问题，不能只选最致密或最疏松区域代表整件。若截取小试样，需考虑切割是否释放残余应力或改变裂纹开口。 记录体素间距、重建矩阵、视场、系统性能评价与最小可靠计量尺度，并说明依据。后面的统计应对低于可靠计量尺度的对象保持谨慎，而不是把所有一体素暗点都算作真实孔。 ### 本章来源 - [Withers 等（2021）：X-ray computed tomography](https://doi.org/10.1038/s43586-021-00015-4)：XCT 实验、重建、分析与可重复性的系统综述。 - [ASTM E1695](https://store.astm.org/standards/e1695)：CT 系统性能测量 ## 04 · 扫描设计：曝光、投影数、转动与校准 从研究目标制定采集方案，说明信噪比、角采样、时间与运动伪影的取舍。 ### 先核对穿透与动态范围 最厚方向应仍有可用透射信号，最薄或空场区域也不能饱和。曝光增加通常可降低相对量子噪声，但热漂移、样品运动和总时长可能变得更重要。平均多帧时应检查帧间稳定性，不能无条件假设噪声都是独立随机量。 试扫可比较不同电压、滤片和曝光下的对比与伪影，同时记录焦点/功率限制。选择依据应是目标缺陷的可见性与测量稳定性，不是让三维渲染“更漂亮”。 ### 角采样不能无限稀疏 二维理想平行束中，若物体横向包含约 N 个有效分辨单元，常用数量级估算是 180° 范围需约 πN/2 个投影。它建立在理想采样假设上，不是实际锥束扫描的通用强制数目。锥束通常采用完整旋转或满足相应条件的采集轨迹。 角度过少会使投影信息欠采样，产生条纹并影响小结构。增加插值后的切片数不能补充缺失角度。应根据设备、物体和任务验证投影数，而不是仅按重建矩阵像素数机械套公式。 ### 固定、几何和参考数据 使用适当的低衰减支撑并保证整个旋转过程稳定。记录旋转轴、源探测距离、探测器偏置和几何校准；中心偏差会造成双边缘等问题。暗场、平场和参考体的采集时间也要保存，便于发现漂移。 原位试验中，夹具本身会衰减并限制角度，加载状态在扫描期间变化会违反静态重建假设。扫描时长与损伤时间尺度应匹配，否则一个体数据可能混合多个状态。 ### 采集记录的最小字段 保存样品编号、坐标、尺寸、材料、目标缺陷、kV/μA、滤片、曝光、帧平均、角度序列、几何、探测器模式、校准文件和异常事件。操作者不应只留下一个重建结果文件。 练习：两次扫描体素同为 20 μm，但一次厚方向几乎无透射，另一次有足够信号。两者检测能力不相等；体素大小相同不能作为可比性的全部依据。 ### 本章来源 - [Kak 与 Slaney：Principles of Computerized Tomographic Imaging](https://www.slaney.org/pct/pct-toc.html)：作者提供的断层成像教材阅读入口。 - [Withers 等（2021）：X-ray computed tomography](https://doi.org/10.1038/s43586-021-00015-4)：XCT 实验、重建、分析与可重复性的系统综述。 - [ASTM E1441](https://store.astm.org/standards/e1441)：CT 原理及术语指南 ## 05 · 重建原理：Radon 变换、FBP、FDK 与迭代法 从线积分的反问题讲起，理解反投影为什么模糊、滤波的作用和锥束算法的条件。 ### 从多个线积分恢复空间分布 二维平行束投影可写 pθ(s)=∫μ(s cosθ−t sinθ, s sinθ+t cosθ)dt。θ 是投影方向，s 是探测器坐标。正向投影把二维图像积成一维曲线，多个角度曲线组成正弦图。一个离中心的点在正弦图中沿正弦轨迹移动。 傅里叶切片定理说明：一个方向投影的一维傅里叶变换，对应原二维图像频域中的一条径向切片。角度不足会留下未充分采样的频域信息，解释了某些条纹与方向性模糊。 ### 为何简单反投影不够 直接把每条投影值沿射线路径涂回去，会让一点的贡献在许多方向扩散，得到模糊图。滤波反投影 FBP 先对投影滤波，再回投影，补偿简单反投影的频率权重。常用 ramp 滤波也会放大高频噪声，因此窗函数与带宽会改变锐度和噪声。 本页二维投影动画使用均匀圆盘与孔的解析弦长，只演示正向投影和角度关系，不声称实现工业三维重建。改变角度可看到一个孔如何在探测器坐标上移动。 ### FDK 的位置 Feldkamp、Davis 与 Kress 1984 年提出的实用锥束算法通常包含几何加权、滤波与反投影，广泛用于圆轨迹锥束 CT。它是实用近似；圆轨迹对一般三维对象并非完备采样，离中心平面较远或锥角较大时可能出现误差。 因此“算法运行成功”不等于所有体素都精确重建。软件的束硬化、环伪影与锐化设置也会影响结果，需保存具体版本与参数以便重复。 ### 迭代方法与正则化 离散形式可写 p≈Aμ，A 描述各射线穿过各体素的路径权重。迭代法最小化数据不一致并加入约束或正则化，例如非负性与平滑性。它可以改善某些噪声或欠采样情况，但先验可能抹平微裂纹或改变边界。 比较算法时用相同的目标量和参考真值检查偏差。若目标是孔隙率，不能只用视觉平滑度选参数；若目标是细裂纹，强平滑后“噪点消失”可能也意味着裂纹消失。 ### 本章来源 - [Kak 与 Slaney：Principles of Computerized Tomographic Imaging](https://www.slaney.org/pct/pct-toc.html)：作者提供的断层成像教材阅读入口。 - [Feldkamp、Davis 与 Kress（1984）：Practical cone-beam algorithm](https://doi.org/10.1364/JOSAA.1.000612)：圆轨迹锥束 CT 常用近似重建算法。 - [Withers 等（2021）：X-ray computed tomography](https://doi.org/10.1038/s43586-021-00015-4)：XCT 实验、重建、分析与可重复性的系统综述。 ## 06 · 伪影排查：束硬化、环、条纹、运动和截断 通过原始投影与三正交切片区分伪影和真实结构，避免在错误重建上继续精细统计。 ### 束硬化与散射 同质圆柱中心灰度低于边缘的杯状效应可能与束硬化有关，高衰减结构之间也可出现条纹。散射会给探测信号叠加不满足理想线积分模型的成分。可考虑适当能谱过滤、参数优化和经验证的校正，但校正本身也需验证。 若中心“低密度区”随扫描参数明显变化，应首先怀疑成像效应。不能把所有暗带都分割成真实孔隙，也不能仅凭某个校正滑块消除了暗带就宣布结果准确。 ### 环伪影与探测器问题 探测器某些像素响应异常或平场不匹配，可能在重建中围绕旋转轴形成环。观察不同高度切片、正弦图和空场可帮助识别。环与真实同心结构可能外观相近，需要结合几何和采集记录判断。 过强的环校正可能删除真实的圆柱形或同心结构。比较校正前后同一坐标的图像，并保存参数。不要只展示校正后的最终截图。 ### 运动、截断与金属条纹 运动会造成边缘重复、模糊或不一致，夹具漂移与加载过程都可能参与。物体超出探测器视场导致投影截断，可能影响内部灰度和边界。高衰减路径光子不足时，对数处理会严重放大误差。 这些问题有的需要重新设计扫描，而不是后期换阈值。必须区分几何/采集错误与分割算法选择。让阈值“吸收”采集缺陷会使结果失去物理解释。 ### 一份排查顺序 先看全角度投影有没有饱和、缺光或运动，再查暗平场和几何；随后看三正交切片，比较结构是否沿真实材料边界连续；最后测试合理的重建参数变化。每次只改变可解释的一组因素。 教学练习：一圈暗环在每张切片上都以旋转轴为中心，而样品偏心放置。它更值得作为环伪影候选检查，不能仅因位于材料内部就计入孔隙。 ### 本章来源 - [Withers 等（2021）：X-ray computed tomography](https://doi.org/10.1038/s43586-021-00015-4)：XCT 实验、重建、分析与可重复性的系统综述。 - [Kak 与 Slaney：Principles of Computerized Tomographic Imaging](https://www.slaney.org/pct/pct-toc.html)：作者提供的断层成像教材阅读入口。 - [ASTM E1695](https://store.astm.org/standards/e1695)：CT 系统性能测量 ## 07 · 数据导入：位深、字节序、切片顺序和空间单位 从文件检查开始建立可重复的数据流程，防止软件导入成功但几何和灰度已经错误。 ### 先识别文件与元数据 重建体可能是 TIFF 序列、RAW 加头文件、DICOM、NRRD 或其他格式。必须知道 Nx、Ny、Nz、数据类型、是否有符号、字节序及头长度。RAW 本身通常不完整保存这些信息，缺少元数据时不能靠“看起来像”确认正确。 切片文件名应按实际 z 顺序排序，避免 1、10、100、11 的字典序错误。重复切片、缺失切片和反向排列会造成形状错误。记录导入前文件数量与预期矩阵，并核查首、中、末切片。 ### 空间标定优先于统计 设置 Δx、Δy、Δz 和单位，必要时记录方向矩阵与原点。体素可以各向异性，体积为 ΔxΔyΔz。显示软件为方便观察可能拉伸图像，统计应使用物理间距，而不是截图像素。 教学体素为 20×20×40 μm³，一个体素体积为 16,000 μm³=1.6×10⁻⁵ mm³。若误设成 20 μm 立方体，所有体积会少一半，z 方向长度也会失真。 ### 保留原灰度与原坐标 不要为方便截图先转成 8 位再做定量。位深转换、窗口范围截断和逐片自动对比度会改变灰度分布，使统一阈值失去意义。显示窗宽窗位可以调整，但应与真实数值变换区分。 裁剪后保存原体中的偏移与坐标变换；后续缺陷表要能定位回原始体。重采样会改变灰度和边界，标签图通常需与连续灰度体采用不同插值策略，且必须核对配准。 ### 软件操作的通用检查 在三正交切片中检查已知几何是否正确，例如圆柱是否被误拉成椭圆、标称长度量级是否合理。再看灰度直方图和若干固定坐标值，确认未发生位深、字节序或符号错误。导入阶段通过后才进入滤波。 保存项目或参数文本：源文件、导入规则、单位、裁剪区域、方向及版本。软件菜单可能变化，最重要的是每项操作都对应一个可检查的数据属性。 ### 本章来源 - [Withers 等（2021）：X-ray computed tomography](https://doi.org/10.1038/s43586-021-00015-4)：XCT 实验、重建、分析与可重复性的系统综述。 - [3D Slicer：Segment Editor 官方文档](https://slicer.readthedocs.io/en/latest/user_guide/modules/segmenteditor.html)：Threshold、Islands 与分割几何的操作依据。 - [ImageJ User Guide：Analyze Menu](https://imagej.net/ij/docs/guide/146-30.html)：标定、测量与颗粒分析的经典文档。 ## 08 · 孔隙分割：材料 ROI、阈值与 Otsu 把外部空气、材料外包络和内部孔隙分开，解释自动阈值的目标与局限。 ### 先定义统计分母 要计算孔隙率，先定义被研究材料的外包络 ROI，包括其中材料和要统计的内部孔隙，排除外部空气、支架以及按研究定义排除的设计孔道。孔隙率 φ=Vvoid/VROI。分母不是整幅矩形图像，也不是仅剩的实体材料体积。 复杂表面、开口孔或设计流道会使外包络定义变得困难。应在分析前规定是否计入开口孔、表面粗糙谷和设计空腔，并保存 ROI 标签。不同分母的结果不能直接比较。 ### 灰度阈值是分类规则 在简单高衬度材料中，可在固定 ROI 内把 g