医学图像处理的难点与问题.ppt

《医学图像处理的难点与问题.ppt》由会员分享，可在线阅读，更多相关《医学图像处理的难点与问题.ppt（45页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、l遥感图像处理遥感图像处理 应用面窄应用面窄l医学图像处理医学图像处理 应用面宽应用面宽一、视觉信息处理的支撑性一、视觉信息处理的支撑性l高度复杂的成像机制高度复杂的成像机制l综合深厚的数学基础综合深厚的数学基础l精确实用的高新技术精确实用的高新技术l重大需求的应用前景重大需求的应用前景二、医学图像处理的前瞻性二、医学图像处理的前瞻性l生理评判的模糊性生理评判的模糊性l信息生成的随机性信息生成的随机性l物理求解的病态性物理求解的病态性l面对应用的实时性面对应用的实时性三、医学图像处理的特殊性三、医学图像处理的特殊性四、医学图像处理的主要问题四、医学图像处理的主要问题 与我们的工作与我们的工作关

2、关于于成成像像重重建建算算法法现现研研究究热热点点集集中中在在MR与与PET成成像像，而而CT或或螺螺旋旋CT等等的的成成像像研研究究，大大多多是是针针对对具具体体问问题题（如如乳乳房房、心心脏脏）提提出出某某种种特特殊殊技技术术去去解解决，不具普适性。决，不具普适性。1、优质重建、优质重建主要问题主要问题采样密集度的改变常导致图像退化。采样密集度的改变常导致图像退化。相应对策相应对策K-K-空间中的任意重采样技术。空间中的任意重采样技术。工作基础工作基础因缺乏成像前数据，我们未开展此类研究。因缺乏成像前数据，我们未开展此类研究。MRMR图像图像存在问题存在问题图图像像分分辨辨率率低低，噪噪声

3、声大大，缺缺少少结结构构信信息息，无法定位病灶。无法定位病灶。相应对策相应对策 a.a.前前处处理理：基基于于小小波波包包门门限限估估计计的的正正则则化抑噪技术。化抑噪技术。b.b.后后处处理理：提提出出各各种种优优化化方方法法，如如正正则则化方法、最大似然法等滤噪。化方法、最大似然法等滤噪。PETPET图像图像已有工作已有工作 a.a.前前处处理理：由由于于缺缺乏乏成成像像前前数数据据，我我们们未未从从事其研究事其研究 b.b.后后处处理理：我我们们提提出出了了正正则则化化参参数数自自适适应应修修改改模模型型与与 迭迭代代快快速速算算法法（ARCTLSARCTLS），并并从从理理论论上上首首

4、次次证证明明了了迭迭代代解解的的存存在在性性与与稳稳定定性性，重重建建质质量量与与速速度度均均比比国国外外最最近近同同类类算算法（法（RCTLSRCTLS）提高）提高3 3倍以上。倍以上。自自适适应应图图像像重重建建(a)退化图像退化图像(b)正则化总体最小二乘法重建图像正则化总体最小二乘法重建图像(c)局部线化局部线化D.F.P方法重建图像方法重建图像(d)自适应正则化总体最小二乘法重建图像自适应正则化总体最小二乘法重建图像 信噪比与计算耗时的比较信噪比与计算耗时的比较RCTLSRCTLSLL-LL-D.F.PD.F.PARCTLSARCTLSSET1SET1ISNRISNR(dB)(dB)

5、1.171.172.502.504.0124.012C.Time C.Time(s)(s)112.6112.654.354.335.535.5注：实验在注：实验在P-II/166 上实现上实现2、伪影校正、伪影校正MRMR成成成成像像像像中中中中，运运运运动动动动伪伪伪伪影影影影几几几几乎乎乎乎无无无无所所所所不不不不在在在在，损损损损害害害害了了了了图像质量。图像质量。图像质量。图像质量。主要问题主要问题 现现有有各各种种前前处处理理技技术术（快快速速扫扫描描、门门控控扫扫描描、欠欠扫扫描描、过过扫扫描描或或可可矫矫正正重重建建模模式式）的的校校正正效效果不理想果不理想相应对策相应对策必必须

6、须提提出出新新的的方方法法，能能具具普普适适性性，又又有有优优越越的的校正校正效果效果。工作基础工作基础(1)(1)前前处处理理：我我们们最最新新提提出出的的逆逆逆逆向向向向迭迭迭迭代代代代修修修修正正正正算算算算法法法法(IIC)(IIC)(IIC)(IIC)对对对对MRMRMRMR图图图图像像像像中中中中的的的的仿仿仿仿真真真真平平平平移移移移运运运运动动动动伪伪伪伪影影影影处处处处理理理理,比比比比美美美美国国国国最最最最近近近近算算算算法法法法(EF)(EF)(EF)(EF)的的的的信信信信噪噪噪噪比比比比提提提提高高高高约约约约个个个个dB,dB,dB,dB,计计计计算算算算速速速速

7、度度度度提提提提高高高高一一一一倍倍倍倍以以以以上上上上,且且且且稳稳稳稳健性好健性好健性好健性好(2)(2)后后处处理理：目目前前我我们们正正想想从从统统计计学学入入手手来来展展开此类研究。开此类研究。a.模板模板原始原始图图像像 b.平移平移伪伪影影图图像像 c.能量聚焦法修正能量聚焦法修正结结果果 d.逆向迭代修正逆向迭代修正结结果果MR平移运动伪影消除平移运动伪影消除a、e为原始图像；为原始图像；b、f为模拟伪影图像；为模拟伪影图像；c、g为能量聚焦法修正的结果；为能量聚焦法修正的结果；d、h为自动逆向迭代法修正的结果为自动逆向迭代法修正的结果EF和和IIC在信噪比和运行时间上的对比在

8、信噪比和运行时间上的对比 方法方法方法方法模板模板图像图像头颅头颅图像图像腹部腹部图像图像信噪比信噪比（dB）EF2.0361.9982.141IIC15.7168.6519.121时间时间（min）EF12.3013.6515.71IIC5.256.868.43注：实验在注：实验在P-III/1G 上实现上实现3、图像配准、图像配准多多多多模模模模态态态态医医医医学学学学图图图图像像像像信信信信息息息息融融融融合合合合是是是是当当当当前前前前研研研研究究究究热热热热点点点点，而前提是配准技术。而前提是配准技术。而前提是配准技术。而前提是配准技术。主要问题主要问题当当前前国国际际主主流流算算法

9、法的的误误配配率率较较大大，难难以以实实用。用。相应对策相应对策提提出出模模糊糊或或广广义义模模糊糊相相似似准准则则来来提提高高配配准准精度。精度。工作基础工作基础(1)(1)模模糊糊弹弹性性配配准准相相似似性性测测度度方方法法正正在在研研究中究中 (2)(2)我我们们提提出出的的医医医医学学学学图图图图像像像像刚刚刚刚性性性性配配配配准准准准的的的的模模模模糊糊糊糊相相相相似似似似性性性性测测测测度度度度函函函函数数数数（FSFSFSFS），稳稳稳稳健健健健性性性性与与与与计计计计算算算算精精精精度度度度明明明明显显显显优优优优于于于于当当当当前前前前同同同同类类类类研研研研究究究究（MIM

10、IMIMI），配准精度提高配准精度提高配准精度提高配准精度提高2 2 2 2倍倍倍倍。实实实实 验验验验 采采采采 用用用用 美美美美 国国国国 VanderbiltVanderbilt大大大大 学学学学Retrospective Retrospective Registration Registration Evaluation Evaluation ProjectProject（RREPRREP）项项项项目目目目组组组组提提提提供供供供的的的的国国国国际际际际通通通通用用用用刚刚刚刚性性性性配配配配准准准准图图图图像像像像数数数数据据据据，评评评评估估估估方方方方式式式式采采采采用用用用配

11、配配配准准准准获获获获得得得得的的的的结结结结果果果果与与与与项项项项目目目目组组组组已已已已有有有有的的的的金金金金标标标标准准准准进进进进行比较。行比较。行比较。行比较。基于模糊边缘场的基于模糊边缘场的CT与与MR刚性配准刚性配准FS与与MI配准方法的误配率比较配准方法的误配率比较方法方法CT-MRPET-MRMI20%44%FS5%22%4、目标分割、目标分割医医医医学学学学图图图图像像像像目目目目标标标标分分分分割割割割也也也也是是是是当当当当前前前前国国国国际际际际研研研研究究究究热热热热点点点点，是是是是精确量化诊断的重要依据。精确量化诊断的重要依据。精确量化诊断的重要依据。精确量

12、化诊断的重要依据。主要问题主要问题分分割割精精度度不不高高，灰灰度度归归类类算算法法没没有有考考虑虑空空间间特性。特性。相应对策相应对策要要提提出出模模糊糊随随机机分分割割模模型型，以以提提高高分分割割的的效效果。果。工作基础工作基础我我我我们们们们提提提提出出出出的的的的广广广广义义义义模模模模糊糊糊糊GibbsGibbsGibbsGibbs分分分分割割割割算算算算法法法法（FGSFGSFGSFGS）可可可可自自自自动动动动获获获获取取取取最最最最优优优优分分分分割割割割效效效效果果果果，且且且且算算算算法法法法稳稳稳稳健健健健性性性性优优优优于国外同类算法（于国外同类算法（于国外同类算法（

13、于国外同类算法（MLMLMLML与与与与CGSCGSCGSCGS）。）。）。）。实实验验数数据据采采用用哈哈佛佛大大学学附附属属医医院院提提供供在在Internet 网网上上的的20幅幅正正常常脑脑部部MR图图像像，分分割割结结果果与与其其提提供供的的手手动动分分割割结结果果（即即金金标标准准）进行对比。进行对比。(a)原始图像原始图像,(b)ML 分割结果分割结果,(c)CGS 分割结果分割结果,(d)FGS分割结果分割结果 （a）（b）（c）（d）脑部冠状面磁共振图像脑部脑部MR图像的分割算法比较图像的分割算法比较 (a)为白质的重叠率为白质的重叠率 (b)为灰质的重叠率为灰质的重叠率5、

14、运动估计、运动估计心心心心脏脏脏脏序序序序列列列列图图图图像像像像分分分分析析析析中中中中，左左左左心心心心室室室室内内内内壁壁壁壁的的的的廓廓廓廓线线线线运动跟踪是心功能量化评价的重要依据。运动跟踪是心功能量化评价的重要依据。运动跟踪是心功能量化评价的重要依据。运动跟踪是心功能量化评价的重要依据。主要问题主要问题轮廓跟踪精度低，误差大。轮廓跟踪精度低，误差大。相应对策相应对策要提出新跟踪算法，解决精度与速度问题。要提出新跟踪算法，解决精度与速度问题。工作基础工作基础最最新新提提出出的的广广义义模模糊糊粒粒子子滤滤波波（GFPF）算算法法，跟跟踪踪效效果果明明显显优优于于当当前前国国际际上上最

15、最好好的的Unscented PF（UPF），跟跟踪踪精精度度提提高高约约一一倍多。倍多。试试验验数数据据由由美美国国NIH心心肺肺血血液液研研究究所所提提供供，对对心心脏脏图图像像的的处处理理结结果果与与其其公公布布的的手手工工勾勾画画左心房内壁（即金标准）轮廓进行比较。左心房内壁（即金标准）轮廓进行比较。心脏序列图像的左心房边缘（心脏序列图像的左心房边缘（NIHNIH金标准）金标准）上行上行:GFPF的跟踪结果的跟踪结果;下行下行:UPF的跟踪结果的跟踪结果跟踪结果误差对比跟踪结果误差对比 均方均方误误差差ErrorError（1 1）ErrorError（3 3）ErrorError（5

16、 5）ErrorError（AverAver）GFPF0.29740.29740.30790.30790.34230.34230.31140.3114UPF0.54010.54010.66440.66440.77080.77080.68320.6832主要问题主要问题医医学学序序列列图图像像的的帧帧内内冗冗余余与与帧帧间间冗冗余余度度很很大大，失失去去一一些些无无关关细细节节不不影影响响诊诊断断，但但迄迄今今医医学学图图像像的的有有损损压缩没有一个公认的标准。压缩没有一个公认的标准。相应对策相应对策要要提提出出新新的的有有损损标标准准及及算算法法，实实现现高高压压缩缩比比下下的的高高信噪比优质

17、解码图像。信噪比优质解码图像。6、图像压缩、图像压缩工作基础工作基础提提提提出出出出了了了了模模模模糊糊糊糊聚聚聚聚类类类类的的的的软软软软分分分分类类类类问问问问题题题题的的的的最最最最优优优优分分分分类类类类测测测测度度度度与与与与算算算算法法法法（OFCOFC），自自自自动动动动实实实实现现现现优优优优化化化化分分分分类类类类，且且且且运运运运算算算算速速速速度度度度比比比比国国国国际通用算法（际通用算法（际通用算法（际通用算法（LBGLBG）提高）提高）提高）提高5656倍。倍。倍。倍。(a)OFC (b)LBGCT_head 码码本大小本大小为为980 (c)OFC (d)LBGMR

18、_head 码码本大小本大小为为970基基于于最最优优模模糊糊矢矢量量量量化化的的图图像像压压缩缩实验结果与分析实验结果与分析图图图图像像像像 码码码码本本本本大小大小大小大小OFCOFC（PSNRPSNR）LBG LBG（PSNRPSNR）图图图图像像像像 码码码码本本本本大小大小大小大小OFC OFC（PSNRPSNR）LBG LBG（PSNRPSNR）MRMR97097041.3041.3032.8232.82CTCT1549154941.7841.7831.6431.6452252235.8035.8031.4031.4098098033.4933.4927.9827.98300300

19、32.8832.8829.9629.9668368329.9729.9726.4126.41（在有损压缩比相同条件下）（在有损压缩比相同条件下）基于模糊聚类优化的序列图像分形压缩基于模糊聚类优化的序列图像分形压缩 测测试试图图像像为为标标准准QCIFQCIF格格式式的的运运动动图图像像序序列列，包包括括由由150150帧帧图图像像(176144)(176144)组组成成的的标标准准序序列列和和2222帧帧医医学学图图像像序序列列（128128128128）。只只对对亮亮度度信信号号进进行行编编码码，色度信号的压缩原理完全相同。色度信号的压缩原理完全相同。表表 基于基于OFCOFC算法不同运动图

20、像序列实验结果比较算法不同运动图像序列实验结果比较测试图像测试图像ForemanForemanForemanForemanNewsNewsNewsNewsContainerContainerContainerContainerQCIF-CardioQCIF-CardioQCIF-CardioQCIF-Cardio算法算法经典经典经典经典算法算法算法算法OFCOFCOFCOFC经典经典经典经典算法算法算法算法OFCOFCOFCOFC经典经典经典经典算法算法算法算法OFCOFCOFCOFC经典经典经典经典算法算法算法算法OFCOFCOFCOFC平均平均PSNRPSNR（dBdB）37.1937.1

21、937.1937.1937.1537.1537.1537.1534.5634.5634.5634.5633.2533.2533.2533.2536.8936.8936.8936.8936.7136.7136.7136.7136.0136.0136.0136.0133.4933.4933.4933.49平平 均均 编编 码码 时时 间间（s s）48.248.248.248.24.24.24.24.215.615.615.615.62.32.32.32.39.339.339.339.332.62.62.62.629.829.829.829.87.67.67.67.6压缩比（倍）压缩比（倍）149

22、.94149.94149.94149.94142.53142.53142.53142.5398.1398.1398.1398.1389.5689.5689.5689.5672.2772.2772.2772.2786.7586.7586.7586.7563.6463.6463.6463.6462.8962.8962.8962.89 两种算法下重建运动图像序列两种算法下重建运动图像序列ForemanForeman的每帧的每帧PSNRPSNR比较比较 1st frame 11th frame 7、边缘检测、边缘检测提出了广义模糊边缘检测算子，是唯一满提出了广义模糊边缘检测算子，是唯一满足国际边缘检测

23、三准则的快速算法足国际边缘检测三准则的快速算法广广义义模模糊糊边边沿沿检检测测对对比比图图（采采用用国国际际标标准准图图像像测测试试）原图象原图象(带加性噪声带加性噪声)GFO算法边缘检测结果算法边缘检测结果Roberts算法边缘检测结果算法边缘检测结果LOG算法边缘检测结果算法边缘检测结果主要问题主要问题医医学学MR图图像像由由于于图图像像采采集集过过程程中中射射频频磁磁场场的的非非均均匀匀性性，造造成成图图像像中中灰灰度度的的非非均均匀匀分分布布，表表现现为为图图像像中中局部过亮与过暗，掩盖了有效的诊断信息。局部过亮与过暗，掩盖了有效的诊断信息。相应对策相应对策要提出非均匀性校正后处理方法

24、。要提出非均匀性校正后处理方法。8、MR非均匀性校正非均匀性校正工作基础工作基础提提出出了了基基于于Gibbs随随机机场场理理论论与与EM算算法法的的MR图像非均匀性校正算法。图像非均匀性校正算法。MR非均匀性的矫正非均匀性的矫正9、血管血管重建重建由由由由于于于于DSADSA只只只只能能能能提提提提供供供供2D2D数数数数据据据据，故故故故无无无无法法法法对对对对病病病病灶灶灶灶与与与与介介介介入入入入路路路路径径径径进进进进行行行行空空空空间间间间定定定定位位位位，为为为为患患患患者者者者带带带带来来来来很很很很大大大大痛痛痛痛苦。苦。苦。苦。主要问题主要问题仅利用仅利用2D2D信息不能导出信息不能导出3D3D血管空间结构。血管空间结构。相应对策相应对策国国际际上上缺缺乏乏此此类类研研究究，我我们们提提出出从从人人体体血血管管标本标本 模模型型出出发发，建建立立3D3D血血管管网网络络数数字字化化模模型型作作为为知知识识工工程程；然然后后以以患患者者2D2D（两两个个以以上上不不同同方方位位）DSADSA图图像像为为引引导导，采采用用智智能能分分析析与与优优化化路径搜寻策略来重建患者真路径搜寻策略来重建患者真3D3D环境环境。工作基础工作基础 由于工程较大，我们现正开始研究由于工程较大，我们现正开始研究精品课件精品课件！精品课件精品课件！谢谢 谢谢