医学三维重建中的表面阶梯效应与几何保真：单例工程验证及自动化质量控制

摘要

目的：解释医学三维模型表面阶梯效应的形成环节，考察以受控几何变化改善表面表达的可行性，并明确自动化质量控制的证据边界。方法：使用一例既有去标识化跟骨训练数据，在修正目标标签遗漏后固定同一基线网格，设置 4、10、20 次滤波遍历的 Taubin 候选及 40 次无体积约束 Laplacian 负对照。计算对应顶点位移、总体和分连通分量体积变化、体积质心偏移、离散粗糙度代理指标以及网格完整性，按预设工程条件选择遍历次数最少的合格候选。结果：三档候选均通过本轮工程筛选；4 次候选的粗糙度代理指标降低 42.64%，最大顶点位移为 0.2493 mm，总体积变化为 −0.00259%，小连通分量体积变化为 −0.23247%。负对照虽然使粗糙度代理指标降低 93.10%，却导致小分量体积减少 58.94%，被规则拒绝。结论：本例支持将平滑作为可追溯、受约束的候选生成步骤，而非逐例凭外观反复修改。结果仅证明相对基线的几何变化可被量化，不证明解剖精度提升或临床可用。

关键词：医学三维重建；阶梯效应；Taubin 平滑；几何保真；连通分量；质量控制

一、研究问题：表面更光滑，是否意味着模型更可靠？

数字模型从 CT 体素到分割标签、等值面网格，再进入网页预览和制造文件，经历多个转换环节。用户看到的“粗糙”可能来自采样间距、分割边界、网格三角片，也可能受到材质、透明度与法线表达影响。若把这些现象统一交给“平滑”处理，容易遗漏真正的结构缺失，并掩盖局部尺寸变化。

医学建模质量需要贯穿影像、分割、模型处理和输出核验。RSNA 3D 打印工作组的原始共识建议关注 STL 生成过程中的平滑与简化，并把导出轮廓重新叠加到源影像检查。^[1] 本文沿用这一质量控制思路，聚焦一个具体工程问题：能否在保留原网格和独立小骨块的条件下，生成有明确偏差记录的表面候选？

本文不提出新的平滑算法。贡献在于将一次实际异常拆分为结构完整性、几何变更和显示效果三个可核对问题，报告同一模型上的参数对照及失败负对照，并给出可自动执行的候选筛选框架。

二、现象定位与形成机制

现场先发现目标标签遗漏，造成主体结构未完整导出。该问题属于导出语义错误：缺少目标结构时，任何平滑都无法恢复缺失的解剖信息。因此，本轮只对标签范围已修正的独立候选进行表面试验，保留早期失败文件，未将错误导出作为有效基线。

修正后的模型仍呈现规律性层纹。源体数据的体素间距约为 0.5391 × 0.5391 × 1.0000 mm；分割标签属于离散体素，直接提取其等值面会继承采样网格的空间特征。等值面提取时应正确使用与数组轴顺序一致的间距，否则还会叠加尺度错误。^[4] 本例的体素各向异性与二值边界离散化，是解释层纹的合理机制；本轮没有做控制变量试验来量化二者各自的贡献。

前期界面核对中，将透明度从约 25% 调至 100% 后仍能看到层纹，提示其并非完全由透明混合造成。法线平均及着色策略能够改变视觉连续性，但不会因此补足缺失标签或证明几何正确。^[5] 本轮保持几何对照条件一致；对比图可见局部棱角趋于缓和，仍存在残余纹理，并未把“消除所有可见层纹”设为目标。

三、材料与试验设计

试验单位为一例既有去标识化训练 CT 对应的数字网格，体数据尺寸为 512 × 512 × 274。本文不披露患者身份、医疗机构、影像截图、原始分割或可下载的病例网格，也不以训练数据测试冒充前瞻性临床研究。侧别及拟使用范围仍需工程师与具备相应资质的临床人员核对。

在同一份分割的目标标签集合修正后，获得包含 64,284 个三角面的基线 STL。网格有两个按面邻接识别的连通分量，体积分别约为 53,413.8355 mm³ 和 211.6299 mm³，总体积为 53,625.4655 mm³。基线通过封闭性和面绕序一致性检查。本轮未进行补洞、减面、骨块复位、重新分割或模型再训练。

Taubin 方法将表面平顺化表述为离散表面信号滤波。^[2] 本轮实际使用 Trimesh 4.12.2 的 filter_taubin，收缩系数 λ＝0.50，膨胀参数 ν＝0.51，比较 4、10、20 次遍历。此处“次”对应软件的单次交替遍历，分别是 2、5、10 组收缩—膨胀配对；不应误读为完整配对的次数。^[3] 三档均从同一原始基线独立生成，没有串行累积处理。

设置无体积约束的 Laplacian 平滑作为收缩负对照：λ＝0.50，40 次遍历。它用于检查筛选规则能否拒绝“更平滑但变形明显”的输出，并非与 Taubin 进行等强度算法优劣评比。候选遍历次数、筛选阈值及最小改动选择规则，在本轮候选计算前写入本地试验计划；这属于内部预设，并非公开预注册。

四、测量定义与筛选条件

对应顶点位移定义为 dᵢ＝‖vᵢ′−vᵢ‖₂，报告最大值、均值和第 95 百分位数。保持相同面连接关系，可避免将最近邻顶点距离误当作一一对应的变形量。本文的 dᵢ 只描述相对于基线的处理位移，不能解释为相对于真实骨表面的误差。

总体及每个连通分量分别计算 ΔV＝100% × (V′−V)/V，并记录体积质心的欧氏偏移。连通分量通过基线面索引保持对应。虽然滤波不主动连接独立分量，分量形状和质心仍可能改变，因此“保留骨块”不能仅凭数量相同来判断。

粗糙度代理指标定义为 R＝meanᵢ‖(AV)ᵢ−vᵢ‖₂，其中 A 为基线网格等权邻接平均算子；相对降幅为 100% × (R₀−R₁)/R₀。该指标反映邻域坐标偏离，不是经临床验证的表面质量评分，也不是曲率、真实噪声或分割准确率。它依赖网格密度及连接关系，本研究只在相同连接关系下比较。

候选还须保持三角面连接关系、封闭性及面绕序一致性，不出现零面积面；另外统计候选面法线与对应基线面法线点积不大于零的面。该法线检查只标记相对于原方向的显著改变，不能替代严格的单元翻转或自相交检测。二进制 STL 写出后重新读取，另行记录数值舍入造成的微小坐标变化。

表 1 的阈值是本轮探索性筛选条件，由本地工程试验设定，没有被解释为临床容差、设备精度要求或适用于所有骨骼的统一参数。选择规则为：满足全部条件时，选取遍历次数最少的候选；若无候选满足，保留原网格并转入复核，不通过临时放宽阈值强行产生结果。

五、结果：轻度候选与负对照

三档 Taubin 候选均通过本轮工程条件，故按预设规则选择 4 次候选。其平均顶点位移为 0.08375 mm，第 95 百分位为 0.16108 mm，最大值为 0.24930 mm。粗糙度代理指标从 0.18524 mm 降至 0.10626 mm，降幅为 42.64%。本文不对单例参数对照进行显著性检验，也不把数万个顶点当作独立病例扩大样本量。

4 次候选的总体积为 53,624.0746 mm³，相对变化为 −0.00259%；大分量变化为 −0.00168%，小分量变化为 −0.23247%。小分量体积质心偏移为 0.00414 mm，包围盒三轴尺寸变化的最大绝对值为 0.11994 mm。两个连通分量及 64,284 个三角面保留，未检出零面积面或上述相对法线方向异常。表中主指标在写出前的浮点网格上测得。

负对照的粗糙度代理降幅达到 93.10%，同时最大顶点位移达到 2.15043 mm，总体积减少 2.70531%，小分量体积减少 58.93791%，并有 15 个面触发相对法线方向检查。因此负对照被拒绝。这表明仅优化外观或全局平滑指标，可能对小结构造成远大于总体积指标所显示的影响。

针对选定候选生成了轴位、冠状位、矢状位各三个位置的内部 CT 轮廓叠加图。该九截面抽查用于查看原网格与候选的空间关系，不是全切片复核、专家盲评或独立金标准验证。原始 CT、分割和基线文件的核验记录保留；输入哈希保护、拒绝覆盖、拒绝把输入当输出、参数遍历约束、负对照拒绝及源文件保持六项程序契约检查通过。

六、讨论：为什么不应每例手工反复平滑

本例说明，可以把表面处理变成一次自动候选生成与质量检查，而不必把反复手工调参设为默认流程。但“每例自动检查”与“每例必须平滑”是两项不同要求。原网格已满足任务需要、关键解剖边缘不能被扰动、输入质量不足或候选超出偏差条件时，都可能应当保留原网格并人工复核。

参数更强时，粗糙度代理指标持续下降，处理位移以及小分量体积变化却随之增大。总体积变化很小，也不能推出每处结构都保持不变。本文因此采用分量级检查与最小改动原则，而不选择视觉上最圆滑或代理指标最低的结果。即使如此，当前算法没有识别关节面、骨折线、薄骨皮质或具有临床意义的尖锐特征；几何门槛通过不等于这些特征已获保护。

显示效果与制造几何应分别管理。改善顶点法线可以使网页呈现更连续，但导出文件的顶点未必改变；改变网格顶点则可能影响切片结果。一个合理的平台应明确标注原始网格、表面候选及已核验制造版本，避免把显示参数变化解释成模型质量已被修复。

结构遗漏与表面层纹也需要分别治理。标签集合、左右侧、空间坐标和文件版本首先决定“导出的是什么”；平滑只影响既有表面的坐标。若目标主体已缺失，后续几何处理没有足够信息补回它。相应地，界面默认选项也不能作为病例真实部位或侧别的证据，应依据同一病例、同一影像的任务记录回显，缺少记录时要求明确选择。

七、面向平台的可追溯自动化流程

建议在导出链中设置独立候选步骤：核对目标标签与坐标；锁定输入指纹；生成有限组表面候选；输出分量级偏差报告；按事先确定的规则选择候选或停止。候选应拥有独立版本及算法参数记录，并与原始分割、原始网格建立可追溯关系。自动化负责重复计算和异常提示，工程及临床复核负责确认具体任务的解剖范围与可接受偏差。

本轮已实现独立处理器及本例验证，尚未证明跨病例通用效果，也未将该处理器部署为正式平台的默认导出策略。通用导出标签策略、前端任务状态和表面候选管理仍需分别完成集成验收。该边界应随发布记录保留，不能从一例技术通过推断整个平台已经获得临床有效性证据。

后续集成应允许不平滑、自动生成候选和人工退回三种结果，并在同一病例中并列查看版本。失败时保留原始文件、错误原因和参数；用户重新扫描或刷新页面不应无提示地重复建模。打印前仍需按目标设备、材料和实机参数完成相应核验。

八、局限性与后续验证

第一，本研究仅一例训练数据，只有一种采样间距组合和一种网格构建来源；参数选择在同一例上完成，没有独立验证集。不得据此宣称对所有跟骨、所有扫描协议或所有部位有效。

第二，基线是经过目标标签范围修正的算法产物，不是专家共识分割、实体计量或真实解剖金标准。0.2493 mm 的最大位移说明“改动多少”，不说明“距离正确答案多远”。源分割中的误差可能被保留，也可能被平滑扩大。

第三，未进行自相交、分量间最小间隙和关键解剖区域的专项检查，也没有全面专家阅片、打印件扫描复测、材料收缩评估或临床使用反馈。封闭性、面绕序及法线检查的通过只能覆盖其各自定义的网格属性。

第四，粗糙度代理指标随网格离散方式变化，不宜跨模型直接比较；九截面抽查存在漏检风险。后续应在预先规定的病例范围中加入独立人工参考、关键区域分层测量、不同采样间距和退化输入测试，再确定按用途配置的参数范围。研究报告应同时披露失败比例及不适用情形。

本文没有据以声明伦理批准或豁免的文件，不作相关认定。若开展正式临床研究、提交学术期刊或公开病例影像，需要另行落实数据使用权限、适用的伦理审查和发表要求。官网本次仅公开方法及去标识化数值汇总，不开放原始医学数据。

九、结论

本例的表面阶梯效应提示，医学三维重建需要同时检查目标结构是否完整、表面修改是否受控以及显示是否清楚。轻度 Taubin 候选以约 0.2493 mm 的最大对应顶点位移，将既定粗糙度代理指标降低 42.64%；负对照则说明，更低的粗糙度并不保证小结构得到保留。

可推进的工程方向，是把有记录、有偏差约束、有失败出口的表面候选生成纳入流程。每例都需要质量检查，但是否进行几何平滑应由输入、任务和验证结果决定。当前成果是一例可复核的技术候选，解剖复核、实机验证及跨病例验证仍是后续工作。

研究透明性声明

数据与代码：本轮本地记录保留输入输出指纹、固定参数、三个候选、负对照及检查结果。由于包含医学数据关联信息，原始影像、分割、病例网格和内部路径不随官网文章发布。公开表格足以核对本文所述数值关系，但不能替代在获得授权的数据上重新运行。

利益关系：文章由子殷科技官网发布，讨论本团队正在研发的流程，存在机构利益相关性。本文不声称获得第三方认证或外部独立验证。写作与程序执行使用 AI 辅助，数值来自实际程序输出；未编造专家评审、伦理批准、基金资助或临床效果。