首页/文章/ 详情

让EIT看得更清:物理引导+扩散模型为电阻抗成像带来新突破

14天前浏览305

文献:A Data-Constrained and Physics-Guided Conditional Diffusion Model for Electrical Impedance Tomography Image Reconstruction,Sensors 2026, 26(5), 1728,DOI: 10.3390/s26051728。

电阻抗成像(EIT)是一种“用电流看见内部结构”的成像技术:在物体边界布置电极,注入微小电流,采集电压响应,再反推出内部电导率分布。它无辐射、成本低、速度快,适合肺通气监测、乳腺筛查、工业两相流和柔性触觉传感等场景。但难点也很典型:EIT逆问题高度非线性、病态,测量噪声和电极接触变化都会让图像变糊、边界漂移、形状失真。


这篇论文提出的 MS-CDM(Multi-Source Conditional Diffusion Model)把扩散模型引入EIT重建,并且不是只依赖单一输入,而是同时利用两类信息:一类是边界电压测量,负责把真实传感数据“拉回现场”;另一类是Gauss-Newton粗重建图像,作为物理先验,负责给生成过程提供结构锚点。

方法亮点

• 多源条件融合:把测量域的电压信号和图像域的物理先验统一编码,减少逆问题解空间的歧义。

• 物理增强损失:训练目标同时考虑数据驱动的去噪精度和EIT前向物理一致性,让结果更符合真实测量。

• Hybrid Swin-Mamba Denoising U-Net:浅层与中层用Swin Transformer捕捉局部边界和纹理,深层用双向Mamba建模全局拓扑结构。

• DDIM快速采样:在保持重建质量的同时降低扩散采样步数,为实时传感应用留出空间。


为什么这个组合有效?EIT图像重建最怕两件事:只有物理模型时,图像平滑、边界模糊;只有深度学习时,模型可能在复杂目标或噪声下“猜错形状”。MS-CDM的思路是让两条线互相补位:电压数据保证来源真实,物理先验保证结构不跑偏,扩散模型则通过逐步去噪恢复细节。

实验怎么做

作者基于EIDORS构建了5万组仿真样本,包含单目标、多目标、复杂凹形等不同拓扑;每帧使用16个边界电极,可得到208个独立边界电压测量。为贴近真实设备,还设置了UEF、KTC和自建系统三类激励协议,并加入从无噪声到10 dB的多级高斯噪声。评价指标包括相对误差RE、相关系数CC、结构相似性SSIM和Dice系数。

结果一:仿真重建更准,复杂形状也不容易散

在仿真数据集上,MS-CDM取得最低相对误差0.218±0.094、最高相关系数0.969±0.034、Dice系数0.956±0.041,整体优于传统TR、CNN、RAU-Net、DHU-Net和仅使用电压条件的CDEIT。更重要的是,它在多目标、心形、星形、月牙等复杂结构中仍能保持目标位置、边界连续性和形态一致性。

结果二:噪声来了,也不是“一碰就碎”

在40 dB到无噪声的高信噪比条件下,MS-CDM表现几乎一致;到20 dB时,Dice仍达到0.867±0.120,说明中等噪声下仍能可靠恢复结构。即使在10 dB强噪声条件下,模型虽然边界细节下降,但仍能保留基本定位能力,没有出现完全失效。

结果三:真实水槽和跨系统数据上仍然稳

论文进一步在UEF2017、KTC2023和自建水槽OUR2026三类真实EIT实验数据上验证。三套平台在硬件、激励协议、电极布局和噪声特性上差异明显,因而很适合考验模型的跨系统泛化能力。


在UEF2017上,MS-CDM取得SSIM 0.908±0.027、Dice 0.872±0.062,边界更清晰、伪影更少。

在KTC2023多系统数据上,MS-CDM相对误差降至0.351±0.095,CC达到0.932±0.031,在高噪声和系统差异条件下仍能保持较稳定的形状恢复。

在自建水槽OUR2026数据上,MS-CDM相对误差为0.358±0.122,CC为0.915±0.057,SSIM为0.907±0.023,Dice为0.901±0.055,说明该方法不仅能在公开数据集上表现好,也具备走向实际实验平台的潜力。

从论文到平台:更适合做“可比较、可验证”的EIT算法评估

如果把这类算法落到工程平台中,关键不是只展示一张好看的重建图,而是要能统一训练、统一推理、统一评价不同网络。多网络架构EIT对比评估平台可以承担这个角色:在训练端配置网络、批大小和学习率;在测试端加载权重,比较原始电压、GN重建、真实标签和深度学习预测结果。


一句话总结

MS-CDM的价值在于:它没有把EIT重建简单交给“黑箱网络”,而是把真实边界电压、物理粗重建和生成式扩散过程放在同一个框架内协同工作。对EIT这种典型病态逆问题来说,这种“数据约束+物理先验+生成式逐步恢复”的路线,既提升了图像清晰度,也增强了噪声鲁棒性和跨系统适应能力。未来如果进一步加速采样、轻量化部署,并扩展到三维和动态EIT,它有望服务于床旁肺监测、工业过程检测和可穿戴传感等更多实时场景。

技术服务支持

围绕EIT图像重建与智能传感应用,我们也可以提供完整技术服务支持:包括MATLAB仿真建模、EIDORS有限元正问题与逆问题建模、不同电极/激励协议下的数据集构建,以及基于Python和PyTorch的网络训练、模型复现、指标评估与可视化对比。无论是论文方法复现、实验平台验证,还是面向具体应用场景的算法定制,都可以进一步开展合作。

来源:320科技工作室
非线性MATLABpythonElectric
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-08-11
最近编辑:14天前
320科技工作室
硕士 | 结构工程师 微❤️ CAE320
获赞 240粉丝 445文章 458课程 0
点赞
收藏
作者推荐

第十四届“宇瞳杯”光学设计大赛赛题2“AR 眼镜拍摄镜头设计”

1.设计指标要求表 1.1设计指标要求 2.系统设计重难点分析AR眼镜拍摄镜头的设计受到多重严苛约束,是微型光学系统设计中的重点与难点,核心矛盾聚集在体积限制、性能要求与量产可行性之间,具体难点如下。首先是空间限制带来的结构性挑战。受AR眼镜外观轻薄与佩戴舒适性约束,镜头系统总长通常被限制在15mm以内,同时需要实现大视场、短焦距且预留足够后截距放置传感器、滤光片及装配间隙,这使得像差校正难度大幅增加,球差、彗差等问题更为突出,加之镜片数量受限,难以实现性能平衡。其次是高精度成像与广角像差的冲突。AR镜头需极高几何精度,畸变需控制在极低水平,否则会导致定位漂移、虚实叠加错位,但大视场广角设计天生易产生桶形畸变、场曲和色差。在极小体积内校正多种像差,需依赖高折射率玻璃、非球面甚至自由曲面镜片,大幅提升了设计与加工成本及难度。再者是与传感器的严格匹配要求。为保证边缘亮度均匀,镜头主光线入射角需严格匹配传感器微透镜,角度偏差会导致边缘照度骤降、偏色,影响算法稳定性。该AR眼镜拍摄镜头设计的核心在于在空间与重量限制下,实现高性能光学成像。设计的核心矛盾在于“长焦”与“轻小”的冲突——既要在有限的15mm总长内达成长焦特性,满足AR眼镜的使用距离;又要通过轻量化设计控制总重,同时维持F/1.8的大光圈进光。这直接导致了设计的核心难点:需在紧凑结构中精准校正大光圈引入的高级像差,兼顾0.8μm极小像素对应的超高解像力与多波长色差校正。因此在对系统进行设计时,在校正系统像差的同时需要考虑长度、重量等对系统的限制。3.设计思路通过查找相关论文、公开专利、镜头手册等资料,寻找与设计参数较为相似的初始结构。该AR眼镜拍摄物镜共有六片光学透镜,第2、4、5为非球面镜片,第1、3片为球面镜片,孔径光阑位于第一面镜子前表面。系统焦距为14.0mm;F数为1.4;波长范围为405~656nm;全像高为9.8mm,系统总长为22.06mm。表3.1系统初始结构参数面序号曲率半径厚度材料折射率材料阿贝数1(光阑)无限-1.97 28.031.991.59368.3423101.534(偶次非球面)14.983.471.71353.9395(偶次非球面)-27.940.776-12.990.671.80822.7647423.240.778(偶次非球面)-5.461.541.71729.5019(偶次非球面)-10.670.3210(偶次非球面)4.853.081.83542.72111(偶次非球面)7.837.12保护玻璃无限0.71.51764.199无限0.114(像面)无限0 图 3.1初始结构光学系统图在完成初始结构选型之后,首先对系统内部的基本参数进行约束,并且对系统进行简单优化。优化后结果如下图所示: 图 3.2初步优化后的初始结构在完成上述优化操作之后选择波前+RMS、波前+PTV等默认优化操作数与像差操作数进行配合,对光学系统像质进行多次迭代优化,逐步校正系统内部像差,并且提高系统的传递函数MTF。因为系统采用了大量的塑料偶次非球面,因此在边缘位置产生了反曲使得系统难以校正。因此在系统所有指标均满足设计要求的情况下,对系统面型进行矫正,使其满足加工工艺要求。4.设计结果及分析系统2D图图 3.13系统结构指标完成情况图 4.2指标操作数 系统点列图 图 4.3系统点列图传递函数MTF(600lp/mm)图 4.4 600lp/mmMTF(全视场)图 4.5 600lp/mmMTF(中心视场)传递函数MTF(500lp/mm)图4.7 500lp/mmMTF(φ5.4视场内)FFT MTFvs.Field图 4.8 FFT MTFvs.Field场曲/畸变图 4.9场曲/畸变相对照度图 4.10相对照度 鬼像分析鬼像是指非预期反射光线在像面形成干扰像的现象,是影响光学系统成像质量的重要因素之一。其主要成因是光线在光学工作表面上的偶次反射,因此在对系统进行鬼像分析时,要重点观察光线在光学表面的反射情况。使用Tracepro进行鬼像分析,系统要求在对系统进行鬼像分析时对90°范围内灯光进行分析,并且鬼像相对能量小于1*10e-6,因此首先设置不同角度的灯光,并且设置能量阈值,分别如下图所示:图 4.11 设置光源光学系统中,大部分杂散光从系统侧面出射,在实际成像时并不会干扰成像的,即不会形成鬼像。所以只需要对这部分杂散光进行吸收或者遮挡,因此对透镜侧面镀吸收膜层,对侧面杂散光进行处理然后设计镜筒及其材料,机械外壳如下图4.14所示:图 4.14机械外壳在完成上述设置之后对光线进行追迹,得到结构以及辐照度/光照图如下图4.15、图4.16所示:图 4.15 2D图图 4.16辐照度/光照度图通过多次对光学系统内部杂散光进行分析,并且对发生反射较多的光学表面度增透膜消除反射后,可以发现系统鬼像减少,且能量较低,集中在1*e-5~1*e-7左右,对系统进行分析并且对第五片透镜前后表面镀0.995膜层。对光线追迹后系统系统2D图、辐照度/光照度图如下图4.25、图4.26所示,可以看出系统中杂散光得到有效抑制,并且鬼像相对能量均低于1*e-6,满足指标要求。图 4.17 系统2D图图 4.18 辐照度/光照度图 5.可行性分析公差分析光学系统最终的公差分配如下表5.2所示:表 5.2公差分配表参数表面公差等级曲率半径(光圈)1-4,6-10Q35Q1厚度(mm)1-2,5-100.002mm3-40.0015mm表面偏心X/Y(mm)1-10Q1表面倾斜X/Y(度)1-10Q3不规则度1-10Q3元件偏心X/Y(mm)1-5Q1元件倾斜X/Y(度)1-5Q3折射率1-10Q3阿贝误差(%)1,3,5,7,9Q3首先对600lp/mm处中心视场MTF进行分析。在Zemax中进行公差分析设置,进行蒙特卡洛分析,运行两百次。如下图5.2所示:图 5.2公差分析在进行公差分析后,得到的公差分析结果里,良品率如下表5.3所示表 5.3公差分析结果概率衍射MTF平均90%0.1917857480%0.2106714650%0.2427752220%0.2734314910%0.28401757由上述分析结果可以看出,AR眼镜拍摄物镜在该公差范围内变化时,衍射MTF平均有90%概率在0.191以上,符合指标要求。然后对600lp/mm处中心视场MTF进行分析,重复上述操作:在进行公差分析后,得到的公差分析结果中,良品率如下表5.4所示表 5.4公差分析结果概率衍射MTF平均90%0.1568401280%0.1747467350%0.2072839220%0.2329036110%0.24329197由上述分析结果可以看出,AR眼镜拍摄物镜在该公差范围内变化时,衍射MTF平均有90%概率在0.156以上,符合指标要求。光学系统图光学系统图如下图5.5所示:图 5.3光学系统图来源:320科技工作室

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈