方法学与计算模型基础
关于MioAtlas 3D心脏模拟器的数学基础、坐标变换、电生理参数及科学文献来源的技术文档。
1. 本文档的目的与范围
本文档阐述了MioAtlas的数学与生物物理学公式推导。在交互式医学与心电模拟中，算法透明度至关重要：每个定量参数、导联向量和传导常数都必须基于经过同行评审的直接文献，或明确注明为几何简化。
公开本规范使学生、教育工作者和研究人员能够准确验证12导联心电图波形是如何合成的、电激动波前是如何插值的，以及冠状动脉闭塞是如何与体表心电图模式相对应的。
开发编年史：模拟器的首个可用版本——3D心脏、冠状动脉、缺血场景与12导联心电图——在一天之内构建完成。传导系统、43种心律失常目录、床旁颈静脉搏动与心音图面板以及增强现实模式是在2026年九月9日至九月12日之间的随后几天内开发的。
2. 冠状动脉解剖重建
冠状动脉血管树在心外膜表面建模了11支动脉：左主干（LM）、左前降支（LAD）、回旋支（LCx）、右冠状动脉（RCA）、第一和第二对角支（D1、D2）、第一和第二钝缘支（OM1、OM2）、锐缘支（AM）、后降支（PDA）以及后侧支（PLB）。
为了在解剖沟中保持连续的曲率切线，结合了两种样条函数公式：
-
三次Bezier样条（LAD和LM）：使用位置控制点（
co）、入切线手柄（handleLeft）及出切线手柄（handleRight）进行重建，通过THREE.CubicBezierCurve3计算并聚合为平滑的THREE.CurvePath。 -
向心Catmull-Rom样条（LCx、RCA及侧支）：从测量的顶点序列中插值生成，并通过
THREE.CatmullRomCurve3(points, false, 'centripetal')进行退化点滤波，防止高曲率线段中的环状伪影。
规范坐标变换
原始曲线在Blender中构建（+Z向上，右手定则）。对于glTF标准下的WebGL渲染（+Y向上，右手定则），在运行时应用规范线性变换：
coordenadas_glTF = [x_blender, z_blender, -y_blender]
该转换在模拟器初始化时通过自动化断言（assertCoordinateConversion）进行验证，确保跨解剖轴的空间对齐一致性。
心肌缺血区域的解剖学验证
6个ST段抬高型心肌梗死（STEMI）区域的划分基于AHA/ACC 17节段模型（Cerqueira et al., 2002）。
缺血体积的校准在标准化条件下逐像素进行了审核：心脏暂停于舒张末期（0.0 s）、100%心肌不透明度以及充分演变的STEMI分期（强度0.98）。利用基线与病理状态之间的对比差分图，优化了LAD近段、LAD中段、LCx近段、RCA近段、RCA中段及PDA近段的椭球中心和半径，排除了坏死向邻近组织的渗漏。
3. 心脏传导系统建模
心脏电激动传导系统在展示全部四个心腔及室间隔的冠状切面上呈现。它包含：
- 17个具解剖学命名的主干结构：窦房结（SAN，界沟）、前结间束、中结间束（Wenckebach束）、后结间束（Thorel束）、通向左心房的Bachmann束、房室结（AVN，Koch三角）、慢径路与快径路、His束、右束支（RBB）、节制带、左束支主干（LBB）、左前分支（LAF）、左后分支（LPF）、左间隔分支（LSF）以及Kent旁路。
- 84支Purkinje纤维分支：将电冲动分布至双侧心室心肌的心内膜下传导网。
组织传导速度
沿样条曲线的电冲动传播遵循基准电生理传导速度：
- 窦房结（传出）：0.05 m/s。
- 心房心肌与结间束：1.0 m/s。
- 房室结（递减延搁）：0.05 m/s。
- His束与His-Purkinje系统：2.0 至 4.0 m/s。
- 心室心肌：0.4 至 0.5 m/s。
传导时相学（窦性心律）
从窦房结始发（t = 0 ms）开始，模型应用生理性电激动时间窗（conductionTimingMs）：
| 节段 / 结构 | 属性 | 时程 | 累积时间窗 |
|---|---|---|---|
| 窦房结传出 | saExit |
5 ms | 0 – 5 ms |
| 心房扩散 | atrialSpread |
60 ms | 5 – 65 ms |
| 房室结延搁 | avDelay |
80 ms | 60 – 140 ms |
| His束 | his |
15 ms | 140 – 155 ms |
| 主束支（RBB/LBB） | branches |
20 ms | 155 – 175 ms |
| 分支（LAF/LPF/LSF） | fascicles |
15 ms | 175 – 190 ms |
| Purkinje纤维网 | purkinje |
15 ms | 190 – 205 ms |
| 心室心肌 | ventricularMyocardium |
85 ms | 205 – 290 ms |
4. 心电图的数学合成
心电监护仪实时合成12个标准体表导联（I、II、III、aVR、aVL、aVF、V1–V6）及5个扩展导联（V7–V9、V3R 与 V4R）。
确定性 120 Hz 累加器
引擎与浏览器刷新率解耦，运行一个固定的 120 Hz 时间累加器（Δt = 8.333 ms），以确保严格的临床校准精度：
- 走纸速度：25 mm/s（1 mm = 0.04 s；5 mm = 0.20 s）。
- 电压：10 mm/mV（1 mm = 0.1 mV；10 mm = 1.0 mV）。
- 12导联面板：2.5 秒（120 Hz 下采样 300 点，72 bpm 下对应 3.0 个心动周期）。
- 连续节律记录带（II）：10.0 秒（120 Hz 下 1,200 个样本）。
模型参考区间
在 72 bpm 的正常窦性心律下，该目录应用经典基准数值（referenceIntervalsMs）：
- P波（
p）：90 ms（参考范围：80–100 ms）。 - PR间期（
pr）：160 ms（参考范围：120–200 ms）。 - QRS波群（
qrs）：90 ms（参考范围：80–100 ms）。 - 基线QT间期（
qt）：在 72 bpm 时为 390 ms（参考范围：350–430 ms）。
QT心率适应：Fridericia公式与Bazett公式对比
为了使QT间期适应变化的心率，模拟器采用Fridericia公式（RR的立方根）而非Bazett公式（平方根）：
QT(RR) = QT_base · (RR / RR_72)^(1/3)
原理依据：Bazett公式在心率加快时过度校正（在心动过速时产生虚假的QTc延长），而在心动过缓时低估。Fridericia公式在全心率谱范围内具有更优越的生理稳定性，将动态QT限制在 200 至 600 ms 之间。
3D胸腔向量与心脏电偶极子
躯干电极放置在真实的3D坐标处：RA [-0.080, +0.070, -0.025]、LA [+0.080, +0.065, +0.025]、LL [0.000, -0.090, +0.055] 以及胸前导联 V1–V6。瞬时导联电压由动态心室电偶极子与单位导联向量之间的标量点积得出（Frank理论，1956）。
5. 模拟器未能涵盖的内容：生理学简化
为了在标准网络浏览器中确保流畅的 60 FPS 性能，MioAtlas相对于真实的生物学复杂性采用了明确的技术简化：
1. 未求解有限元离子微分方程
引擎建模运动学波前传播与连续几何积分，而非通过双域模型（Hodgkin-Huxley 或 Luo-Rudy）计算跨膜离子通道动力学。
2. 优势型建模为三种离散模式，不包含开口异常
模拟器建模了三种经典的冠状动脉优势型配置（右优势型 ~85%、左优势型 ~8%、均衡型 ~7%），改变从属分支起点（PDA与PLB）及受累缺血区域。然而，它不包含开口起源异常、心肌桥或复杂的先天性解剖变异。
3. 缺血通过几何体积建模，而非微血管流体动力学
心肌梗死区域和ST段偏移受校准椭球体和矩阵表格控制，而非体循环血压流体动力学或侧支循环血流。
4. 机械耦联通过变形目标（morph targets）实现，而非主动弹性
心室收缩（收缩期、心尖扭转、心房收缩辅助）由解析形态键插值驱动，而不计算Frank-Starling曲线环或后负荷依赖性心室弹性。
5. 向量投影未建立非均匀躯干容积导电体模型
导联波形通过将心脏电偶极子投影到导联向量上合成，而未在不同电导率介质（肺、肋骨、骨骼肌）之间求解3D拉普拉斯方程。
6. 体格检查通过时相事件表生成
颈静脉搏动与心音图使用时钟同步的时相波形函数生成，而非来自体内心腔压力与声学换能器的实际记录。
6. 科学来源与参考书目文献
-
Kligfield P, Gettes LS, Bailey JJ, et al. (2007). "Recommendations for the standardization and interpretation of the electrocardiogram: Part I: The electrocardiogram and its technology". Circulation. 2007;115(10):1306-1324.
-
Surawicz B, Childers R, Deal BJ, et al. (2009). "Recommendations for the standardization and interpretation of the electrocardiogram: Part III: Intraventricular conduction disturbances". Circulation. 2009;119(10):e235-e240.
-
Rautaharju PM, Surawicz B, Gettes LS, et al. (2009). "Recommendations for the standardization and interpretation of the electrocardiogram: Part IV: The ST segment, T and U waves, and the QT interval". Circulation. 2009;119(10):e241-e250.
-
Wagner GS, Macfarlane P, Wellens H, et al. (2009). "AHA/ACCF/HRS Recommendations for the Standardization and Interpretation of the Electrocardiogram: Part VI: Acute Ischemia/Infarction". Circulation. 2009;119(10):e262-e270.
-
Cerqueira MD, Weissman NJ, Dilsizian V, et al. (2002). "Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart". Circulation. 2002;105(4):539-542.
-
Anderson RH, Yanni J, Boyett MR, Chandler NJ, Dobrzynski H. (2009). "The anatomy of the cardiac conduction system". Clin Anat. 2009;22(1):99-113.
-
Frank E. (1956). "An accurate, clinically practical system for spatial vectorcardiography". Circulation. 1956;13(5):737-749.
-
Fridericia LS. (1920). "Die Systolendauer im Elektrokardiogramm bei normalen Menschen und bei Herzkranken". Acta Med Scand. 1920;53:469-486.
-
Bazett HC. (1920). "An analysis of the time-relations of electrocardiograms". Heart. 1920;7:353-370.
-
Sovilj S, Magjarević R, Lovell NH, Dokos S. (2013). "A simplified 3D model of whole heart electrical activity and 12-lead ECG generation". Comput Math Methods Med. 2013;2013:134208.
-
Thygesen K, Alpert JS, Jaffe AS, et al. (2019). "Fourth Universal Definition of Myocardial Infarction (2018)". Eur Heart J. 2019;40(3):237-269.
7. 审校状态与技术合作
未经心血管专科医师审校。
MioAtlas完全由医学学士、人工智能顾问 Guillermo Nimubona Castrillo 独立设计并实现。本项目尚未经临床伦理委员会或医学学会审核。
如果您是心血管内科、心脏电生理或心血管解剖学领域的专科医师，并有意愿审核这些模型，您的审校贡献将署名并注明您的专业及所属医疗机构。联系方式：contacto@mioatlas.com。
本内容仅供医学教学参考。并非临床诊断工具，不能替代专业医务人员的临床评估。 医疗声明.