《医学图像的配准与融合》由会员分享，可在线阅读，更多相关《医学图像的配准与融合（120页珍藏版）》请在金锄头文库上搜索。
1、第八章第八章 医学图像的配准与融合医学图像的配准与融合 8.1 配准与融合的应用背景介绍配准与融合的应用背景介绍 8.2 医学图像配准概述医学图像配准概述 8.3 图像配准的理论基础图像配准的理论基础 8.4 常用的医学图像配准方法常用的医学图像配准方法 8.5 图像配准的评估图像配准的评估步喧枕潭党肄宫躁钉校坐澈套悔桥抒腋蓉芯轨胞干息筷即沼宝钵冈毫莫涉医学图像的配准与融合医学图像的配准与融合第八章第八章 医学图像的配准与融合医学图像的配准与融合 8.6 图像融合概述图像融合概述 8.7 常用的图像融合方法常用的图像融合方法 8.8 图像融合效果的评价图像融合效果的评价压钉壕前灌蹬评瑰训冉骆恃
2、吗描穗姥候恒赃服篮脚稿君守靶降供禽巧扬盟医学图像的配准与融合医学图像的配准与融合 随着计算机技术的高速发展，医学成像技术日随着计算机技术的高速发展，医学成像技术日新月异，新月异，为临床医学提供了各种形态和功能的影像为临床医学提供了各种形态和功能的影像信息。信息。 但是各种成像技术和检查方法都有它的优势与但是各种成像技术和检查方法都有它的优势与不足，并非一种成像技术可以适用于人体所有器官不足，并非一种成像技术可以适用于人体所有器官的检查和疾病诊断，也不是一种成像技术能取代另的检查和疾病诊断，也不是一种成像技术能取代另一种成像技术，而是相辅相成、相互补充。一种成像技术，而是相辅相成、相互补充。 8
3、.1 8.1 应用背景介绍应用背景介绍倔匠召喀耳喘财蛀资栖夏襟晨背剔戚晤侣掏肠胰婚惟极薪咕裕庚盘促撇叫医学图像的配准与融合医学图像的配准与融合根据医学图像所提供的信息内涵，分为两大类：根据医学图像所提供的信息内涵，分为两大类：l 解剖结构图像（解剖结构图像（CTCT、MRIMRI、B B超等）超等）l 功能图像（功能图像（SPECTSPECT、PETPET等）等）解剖图像以较高的分辨率提供了脏器的解剖形态信息（功能解剖图像以较高的分辨率提供了脏器的解剖形态信息（功能图像无法提供脏器或病灶的解剖细节），但无法反映脏器的图像无法提供脏器或病灶的解剖细节），但无法反映脏器的功能情况。功能情况。功能图
4、像分辨率较差，但它提供的脏器功能代谢信息是解剖功能图像分辨率较差，但它提供的脏器功能代谢信息是解剖图像所不能替代的；图像所不能替代的；8.1 8.1 应用背景介绍应用背景介绍牢靖在盒肯文持逮傍咳栏玩任蟹疚倪畴炯抬烦辨涨癸渣粳芹续疮蛙墒希鄂医学图像的配准与融合医学图像的配准与融合n 目前这两类成像设备的研究都已取得了很大的进步，图像的目前这两类成像设备的研究都已取得了很大的进步，图像的空间分辨率和图像质量有很大的提高，但由于成像原理不同所空间分辨率和图像质量有很大的提高，但由于成像原理不同所造成的图像信息局限性造成的图像信息局限性, ,使得单独使用某一类图像的效果并不理使得单独使用某一类图像的效
5、果并不理想。想。n 因此，为了提高诊断正确率，需要综合利用患者的各种图像因此，为了提高诊断正确率，需要综合利用患者的各种图像信息。信息。 8.1 8.1 应用背景介绍应用背景介绍转技讽基筒窟伎藐邱佣芋耶说农羚愧扦喇刨卧艳盂剖焉雕淫调艾汞三翰宁医学图像的配准与融合医学图像的配准与融合n 最有效的解决方法：以医学图像配准技术为基础，利用信息最有效的解决方法：以医学图像配准技术为基础，利用信息融合技术，将这两种图像结合起来，利用各自的信息优势，在融合技术，将这两种图像结合起来，利用各自的信息优势，在一幅图像上同时表达来自人体的多方面信息。一幅图像上同时表达来自人体的多方面信息。n 更加直观地提供了人
6、体解剖、生理及病理等信息。其中配准更加直观地提供了人体解剖、生理及病理等信息。其中配准技术是图像融合的先决条件技术是图像融合的先决条件, ,必须先进行配准必须先进行配准, ,才能实现准确地才能实现准确地融合。融合。8.1 8.1 应用背景介绍应用背景介绍厦郴粤糊晨糠椎歉揩兆犹肄奔幽诅晦掉巡菊碾球翌砰储泼券郭尘纳黑麦谆医学图像的配准与融合医学图像的配准与融合二、医学图像配准在临床上的应用二、医学图像配准在临床上的应用 医医学学图图像像配配准准具具有有很很重重要要的的临临床床应应用用价价值值。对对使使用用各各种种不不同同或或相相同同的的成成像像手手段段所所获获得得的的医医学学图图像像进进行行配配准
7、准不不仅仅可可以以用用于于医医疗疗诊诊断断，还还可可用用于于手手术术计计划划的的制制定定、放放射射治治疗疗计计划划的的制制定定、病病理理变变化化的的跟跟踪踪和和治治疗疗效效果果的的评评价价等等各各个个方方面。面。彼描慈乌跌巩夸滩范菊屠藏友拘伶副董责寂词洋锥案耕镇租揭笺眶屹补肚医学图像的配准与融合医学图像的配准与融合临床应用举例临床应用举例 1 1、计计算算机机辅辅助助手手术术中中，外外科科医医生生根根据据配配准准的的CT/MR/DSACT/MR/DSA图图像像精精确确定定位位病病灶灶及及周周围围相相关关的的解解剖剖结结构构信信息息，设设计计出出缜缜密密的的手手术术计计划划。在在手手术术过过程程
8、中中，利利用用三三维维空空间间定定位位系系统统使使术术前前计计划划的的虚虚拟拟病病人人、手手术术台台上上的的真真实实病病人人和和手手术术器器械械三三者者精确联系起来进行手术跟踪。精确联系起来进行手术跟踪。 二、医学图像配准在临床上的应用二、医学图像配准在临床上的应用庭赘序屿调婉筋戚磺浦瓮眼写虑男空讨狞蝴锭妆热咸菊级驾速使该霖劳慧医学图像的配准与融合医学图像的配准与融合临床应用举例临床应用举例 2 2、在在癫癫痫痫病病的的治治疗疗中中，一一方方面面需需要要通通过过CTCT，MRIMRI等等图图像像获获得得病病人人的的解解剖剖信信息息，另另一一方方面面又又需需要要通通过过SPECTSPECT或或P
9、ETPET等等得得到到病病人人的的功功能能信信息息，这这两两方方面面的的结结合合将将有有助助于于对对病病人人的的精确治疗。精确治疗。二、医学图像配准在临床上的应用二、医学图像配准在临床上的应用虐堕假坷瓜亲胰新阴臣风要怨昔闹氏狱榆甲逛腰砍皑快农托魁陵互谣严朽医学图像的配准与融合医学图像的配准与融合临床应用举例临床应用举例 3 3、放放射射治治疗疗中中，应应用用CTCT和和MRMR图图像像的的配配准准和和融融合合来来制制定定放放疗疗计计划划和和进进行行评评估估，用用CTCT图图像像精精确确计计算算放放射射剂剂量量，用用MRMR图图像像描描述述肿肿瘤瘤的的结结构构，用用PETPET和和SPECTSP
10、ECT图图像像对对肿肿瘤瘤的的代代谢谢、免免疫疫及及其其他他生生理理方方面面进进行行识识别别和和特特性性化化处处理理，整整合合的的图图像像可可用于改进放射治疗计划或立体定向活检或手术。用于改进放射治疗计划或立体定向活检或手术。 二、医学图像配准在临床上的应用二、医学图像配准在临床上的应用血唆挪焙拭恋具筹寅缴寻傣弄焦坛衣朵决老郝碧遥糯淋症孵蹬堰卒缴谊潘医学图像的配准与融合医学图像的配准与融合8.2 8.2 医学图像配准概述医学图像配准概述一、医学图像配准的概念一、医学图像配准的概念 二、医学图像配准方法的分类二、医学图像配准方法的分类三、医学图像配准的基本过程三、医学图像配准的基本过程 卞传壳禹
11、钳祝趣垂扼雄番莱嫌锁货戒又戊省僧秒屁枝烫件霖裁带减拣望拟医学图像的配准与融合医学图像的配准与融合一、医学图像配准的概念一、医学图像配准的概念 医医学学图图像像配配准准是是指指对对于于一一幅幅医医学学图图像像寻寻求求一一种种（或或一一系系列列）空空间间变变换换，使使它它与与另另一一幅幅医医学学图图像像上上的的对对应应点点达达到到空空间间上上的的一一致致。这这种种一一致致是是指指人人体体上上的的同同一一解解剖剖点点在在两两张张匹匹配配图图像像上上有有相相同同的的空空间间位位置置（位位置置一一致致，角角度度一一致致、大大小一致）。小一致）。 配配准准的的结结果果应应使使两两幅幅图图像像上上所所有有的
12、的解解剖剖点点，或或至至少少是是所有具有诊断意义的点及手术感兴趣的点都达到匹配。所有具有诊断意义的点及手术感兴趣的点都达到匹配。 很瘁秤邪避硬亢诲滦唱霄顶砂春宠吨萧地青迎物落饯丝今挚湍肆篓摈涝澈医学图像的配准与融合医学图像的配准与融合医学图像配准示意图医学图像配准示意图一、医学图像配准的概念一、医学图像配准的概念装疟彬率坏贵栗乒邻羞渭傈起嘱颗铁尹准绕包剩遏脏没盒柴沉酷握瘸哉馋医学图像的配准与融合医学图像的配准与融合二、医学图像配准方法的分类二、医学图像配准方法的分类 到到目目前前为为止止，图图像像配配准准方方法法的的分分类类始始终终没没有有一一个个统统一一的的说说法法。目目前前比比较较流流行行
13、的的是是19931993年年Van Van den den ElsenElsen等等人人对医学图像配准进行的分类，归纳了七种分类标准。对医学图像配准进行的分类，归纳了七种分类标准。锚曹徊聊她磕绩苛抒笺包搁逾暂讹墓云攀奠介庆腐资觉俊萎灭裂蛹钨弘葛医学图像的配准与融合医学图像的配准与融合（一）按图像维数分类（一）按图像维数分类 按图像维数分为按图像维数分为2D/2D2D/2D，2D/3D2D/3D，以及，以及3D/3D3D/3D配准。配准。2D/2D2D/2D配准通常指两个断层面间的配准；配准通常指两个断层面间的配准；2D/3D2D/3D配准通常配准通常指空间图像和投影图像（或者是单独的一个层面）
14、间的指空间图像和投影图像（或者是单独的一个层面）间的直接配准；直接配准；3D/3D3D/3D配准指配准指2 2幅三维空间图像间的配准。幅三维空间图像间的配准。二、医学图像配准方法的分类二、医学图像配准方法的分类铆脉迄否坑创违锨舰寒裁惜症篷卞碾锹黍笆臼婿矽潭勿牧悍浓犀酒钟苛休医学图像的配准与融合医学图像的配准与融合（二）根据医学图像的模态分类（二）根据医学图像的模态分类 单模态医学图像配准：是指待配准的两幅图像是用同一单模态医学图像配准：是指待配准的两幅图像是用同一种成像设备获取的。一般用在生长监控、减影成像等。种成像设备获取的。一般用在生长监控、减影成像等。 多模态图像配准：是指待配准的两幅图
15、像来源于不同的多模态图像配准：是指待配准的两幅图像来源于不同的成像设备，主要应用于神经外科的诊断、手术定位及放成像设备，主要应用于神经外科的诊断、手术定位及放疗计划设计等。疗计划设计等。 二、医学图像配准方法的分类二、医学图像配准方法的分类崩善决栗薯蒋秋藉衡解霉盒窃内庙蔡倡耘嚷沦樟卿哇难倍烩撰堕超杨侩礼医学图像的配准与融合医学图像的配准与融合（三）根据变换性质分类（三）根据变换性质分类 刚性变换：只包括平移和旋转。刚性变换：只包括平移和旋转。 仿射变换：将平行线变换为平行线。仿射变换：将平行线变换为平行线。 投影变换：将直线映射为直线。投影变换：将直线映射为直线。 曲线变换：将直线映射为曲线曲
16、线变换：将直线映射为曲线。二、医学图像配准方法的分类二、医学图像配准方法的分类误参整傣告铜擞晰炳聘怕慢窄嘉呸赴珍竹皮迁茨鹿兢秤倔掀赴角括澎瘦恐医学图像的配准与融合医学图像的配准与融合（四）根据用户交互性的多少分类（四）根据用户交互性的多少分类 自动配准：用户只需提供相应的算法和图像数据。自动配准：用户只需提供相应的算法和图像数据。 半自动配准：用户需初始化算法或指导算法半自动配准：用户需初始化算法或指导算法( (如拒绝或接如拒绝或接受配准假设受配准假设) )； 交互配准：用户在软件的帮助下进行配准交互配准：用户在软件的帮助下进行配准 二、医学图像配准方法的分类二、医学图像配准方法的分类咆符泅傀
17、牢弹哟撒遍缔异臻吭廖秒镣遵完循质秒租氨拇萝误腻体缮皮脉泌医学图像的配准与融合医学图像的配准与融合（五）根据配准所基于的图像特征分类（五）根据配准所基于的图像特征分类 基于外部特征的图像配准：是指在研究对象上设置一些基于外部特征的图像配准：是指在研究对象上设置一些 标志点，使这些标记点能在不同的影像模式中标志点，使这些标记点能在不同的影像模式中 显示，然后再用自动、半自动或交互式的方法显示，然后再用自动、半自动或交互式的方法 用标记将图像配准。用标记将图像配准。 基于内部特征的配准：主要包括三个方面：基于标记的基于内部特征的配准：主要包括三个方面：基于标记的 配准方法、基于分割的配准方法、基于像
18、素特配准方法、基于分割的配准方法、基于像素特 性的配准性的配准。 二、医学图像配准方法的分类二、医学图像配准方法的分类侧洱妙馋锈锭谐偷胺膀霸止贰等器肿寂寺坝洼岂岔栖运肆胺盂切挑贤宿锦医学图像的配准与融合医学图像的配准与融合（六）根据配准过程中变换参数确定的方式分类（六）根据配准过程中变换参数确定的方式分类 1 1、通过直接计算公式得到变换参数的配准：限制在基于特征信息、通过直接计算公式得到变换参数的配准：限制在基于特征信息( (例如小数目的特征点集、二维曲线、三维表面例如小数目的特征点集、二维曲线、三维表面) )的配准应用中。的配准应用中。 2 2、通过在参数空间中寻求某个函数的最优解得到变换
19、参数的配准：、通过在参数空间中寻求某个函数的最优解得到变换参数的配准： 所有的配准都变成一个能量函数的极值求解问题。所有的配准都变成一个能量函数的极值求解问题。 二、医学图像配准方法的分类二、医学图像配准方法的分类拓稳钠肝拿话卯陶赤痛藻肃媳涩喂征陌除蚁栈兑鸵缴敝范竣原座萤壤围头医学图像的配准与融合医学图像的配准与融合（七）根据主体分类（七）根据主体分类 1 1、同一患者、同一患者(Intrasubject)(Intrasubject)的的 配准：指将来自同一个病配准：指将来自同一个病人的待配准图像，用于任何种类的诊断中；人的待配准图像，用于任何种类的诊断中； 2 2、不同患者、不同患者(Int
20、ersubject) (Intersubject) 的配准：指待配准图像来自的配准：指待配准图像来自不同病人，主要用在三维头部图像不同病人，主要用在三维头部图像(MR(MR、CT)CT)的配准中的配准中 3 3、 患者与图谱的患者与图谱的(Atlas)(Atlas)图像配准。是指待配准图像一幅图像配准。是指待配准图像一幅来自病人，一幅来自图谱。来自病人，一幅来自图谱。二、医学图像配准方法的分类二、医学图像配准方法的分类滥删蹭素退材庄侥胆找峪澳订酬于浓闲殆箱耐烹惶弛版透粪荡鞠亢迭甚态医学图像的配准与融合医学图像的配准与融合三、医学图像配准的基本过程三、医学图像配准的基本过程 1 1、根据待配准图
21、像（浮动图像）、根据待配准图像（浮动图像）I I2 2与参考图像（基准图与参考图像（基准图像）像）I I1 1，提取出图像的特征信息组成特征空间；，提取出图像的特征信息组成特征空间； 2 2、根据提取出的特征空间确定出一种空间变换，使待配、根据提取出的特征空间确定出一种空间变换，使待配准图像准图像I I2 2经过该变换后与参考图像经过该变换后与参考图像I I1 1能够达到所定义的能够达到所定义的相似性测度；相似性测度； 3 3、在确定变换的过程中，还需采取一定的搜索策略也就、在确定变换的过程中，还需采取一定的搜索策略也就是优化措施以使相似性测度更快更好地达到最优值。是优化措施以使相似性测度更快
22、更好地达到最优值。幕笆院圣秋颁狄掺踞审津乙蔷扭物博锥扣圣尽选悍监哦缅训瞻戚牢邑疼荤医学图像的配准与融合医学图像的配准与融合NY待配准图像I2几何变换T变换后的图像I*=T(I2)相似性测度评价T最优？最优T更新T初始变换T参考图像I1图像配准的流程图图像配准的流程图三、医学图像配准的基本过程三、医学图像配准的基本过程哮王凳咕价崩擒个尉扭纂俄淑赴逢欢脱缕蛰偿囊越捣标畸亚撤硒硒古功涤医学图像的配准与融合医学图像的配准与融合8.3 图像配准的理论基础图像配准的理论基础一、图像配准原理一、图像配准原理 二、空间变换二、空间变换 三、参数的优化搜索三、参数的优化搜索 四、插值方法四、插值方法 五、相似性
23、测度五、相似性测度 驴啃碘惦塞靳了藤玛缘酗讶鸽囱牲潭统埂食份童剃辨寓舅萝润萝名涣勤盾医学图像的配准与融合医学图像的配准与融合一、图像配准原理一、图像配准原理 对于在不同时间或对于在不同时间或/ /和不同条件下获取的两幅图像和不同条件下获取的两幅图像A(x)A(x)和和B(x)B(x)的配准，就是要定义一个相似性测度并寻找的配准，就是要定义一个相似性测度并寻找一个空间变换关系，使得经过该空间变换后两幅图像间一个空间变换关系，使得经过该空间变换后两幅图像间的相似性达到最大（或者差异性最小）。即使图像的相似性达到最大（或者差异性最小）。即使图像A A上的上的每一个点在图像每一个点在图像B B上都有唯
24、一的点与之对应，并且这两点上都有唯一的点与之对应，并且这两点应对应同一解剖位置。应对应同一解剖位置。 届权鹅撰赶磐振噬松拥煮断我绳棱惺磕掷牲任奏纷思遵弘遂拂翘诗架伏搜医学图像的配准与融合医学图像的配准与融合S S是相似性测度，配准的过程归结为寻求最佳空间变换的过程。是相似性测度，配准的过程归结为寻求最佳空间变换的过程。 公式表示公式表示: : 由于空间变换包含多个参数，是一个多参数最优化问题，一由于空间变换包含多个参数，是一个多参数最优化问题，一般由迭代过程实现：般由迭代过程实现：一、图像配准原理一、图像配准原理车涪什憨田趣拔嗡藤戈茵销斗县速烘影送佃拉娶推镣闪揖甩步指凉焉肋阔医学图像的配准与融
25、合医学图像的配准与融合二、空间变换二、空间变换 图像图像A A和和B B的配准就是寻找一种映射关系的配准就是寻找一种映射关系T T：X XAAX XB B，使得使得X XA A上的每一点在上的每一点在X XB B上都有唯一的点与之对应。上都有唯一的点与之对应。 这种映射关系表现为一组连续的空间变换，如整幅图像这种映射关系表现为一组连续的空间变换，如整幅图像应用相同的空间变换，则称之为全局变换（应用相同的空间变换，则称之为全局变换（global global transformationtransformation），否则，称之为局部变换（），否则，称之为局部变换（local local tra
26、nsformationtransformation）。）。 囚雕霞见法佛茨尼巧吐粘素肥彬渣贴熔姓疥耪奏尹传掸桶妖榨放氛拭禾捕医学图像的配准与融合医学图像的配准与融合 图像配准的基本变换图像配准的基本变换二、空间变换二、空间变换痉瓜批眷颖敬戒篓见寸垫罢潦捐床加受若臆碾悉寸它餐匙不浇驰羔搏琴酸医学图像的配准与融合医学图像的配准与融合二、空间变换二、空间变换刚体变换刚体变换 刚体：刚体：是指物体内部任意两点间的距离保持不变。是指物体内部任意两点间的距离保持不变。 刚体变换：刚体变换：使得一幅图像中任意两点间的距离在变换前使得一幅图像中任意两点间的距离在变换前 后保持不变。后保持不变。 例如：人体的头
27、部由坚硬的颅骨支撑，在处理时通常忽例如：人体的头部由坚硬的颅骨支撑，在处理时通常忽 略头部皮肤的微小变形，将整个人脑看作是一个略头部皮肤的微小变形，将整个人脑看作是一个 刚体。刚体。绎赔醋高量蜘晒宜眶试担访毗什魄烦角玉谅狙九癣健做汹豌债起刹东裕长医学图像的配准与融合医学图像的配准与融合两幅图像之间的刚体变换可由一个刚体模型描述：两幅图像之间的刚体变换可由一个刚体模型描述：s s是比例变换因子是比例变换因子。 是图像之间沿是图像之间沿x x，y y，z z方向上的平移量。方向上的平移量。R R是是3 33 3的旋转矩阵，满足约束条件：的旋转矩阵，满足约束条件： 二、空间变换二、空间变换刚体变换刚
28、体变换佃投竣埃寞锋繁耙般钧宋捍绚阜饭橇混幼鸿请氢逢诸仕遁吹十据挠玖详迅医学图像的配准与融合医学图像的配准与融合相对笛卡尔坐标系的三个坐标轴，相对笛卡尔坐标系的三个坐标轴，R R有三种不同的形式：有三种不同的形式： 分别表示围绕分别表示围绕坐标轴的旋转角度坐标轴的旋转角度 二、空间变换二、空间变换刚体变换刚体变换骑屡簧焊底君晃梗捐吃睹水粮冲慷阻寅脯竞跺筐厂闻劲涣证碳吵到芝疽阴医学图像的配准与融合医学图像的配准与融合二、空间变换二、空间变换- -仿射变换仿射变换仿射变换：仿射变换：将直线映射为直线，并保持平行性。将直线映射为直线，并保持平行性。不满足：不满足：志建当象脊九喷联此瞪蜘委企外蒲报厨宽溜
29、挎椿钻赂层德掂铝朋票鸥俩茧医学图像的配准与融合医学图像的配准与融合在笛卡儿坐标系下，二维仿射变换的旋转矩阵在笛卡儿坐标系下，二维仿射变换的旋转矩阵R R表示为：表示为：三维：三维：二、空间变换二、空间变换- -仿射变换仿射变换富悔沮屎财监带私肖瞥棵辱啃汗集饥殖瘸痔窖笋颧膀诬枚燎沫礼室蛤厩瞩医学图像的配准与融合医学图像的配准与融合二、空间变换二、空间变换投影变换投影变换投影变换：投影变换：将直线映射为直线，但不保持平行性质。将直线映射为直线，但不保持平行性质。投影变换主要用于二维投影图像与三维体积图像的配准。投影变换主要用于二维投影图像与三维体积图像的配准。 二维投影变换按照下式将图像二维投影变
30、换按照下式将图像映射至图像映射至图像：是依赖于图像本身的常数。是依赖于图像本身的常数。 辰迢嘉搅衰靡驮鸵蛙唯条妙乡攀雇痘愉奏恶饥萤掏沈传镶典术凡恤粗拄车医学图像的配准与融合医学图像的配准与融合二、空间变换二、空间变换非线性变换非线性变换非线性变换是把直线变换为曲线。它反映的是图像中组织非线性变换是把直线变换为曲线。它反映的是图像中组织或器官的严重变形或位移。典型的非线性变换是多项式函或器官的严重变形或位移。典型的非线性变换是多项式函数，如二次、三次函数及薄板样条函数。有时也使用指数数，如二次、三次函数及薄板样条函数。有时也使用指数函数。非线性变换多用于使解剖图谱变形来拟合图像数据函数。非线性变
31、换多用于使解剖图谱变形来拟合图像数据或对有全局性形变的胸、腹部脏器图像的配准。或对有全局性形变的胸、腹部脏器图像的配准。拥桑脾掘肩粉罚拂扇烧绞升烹庸玩病滋侨俯砰胶丽陶诅顺找褂邮寝瞄椽地医学图像的配准与融合医学图像的配准与融合二、空间变换二、空间变换非线性变换非线性变换1. 1. 二阶多项式变换二阶多项式变换 吮赦灵擞域功卞狐溪朴弦翁犬阎乳实淀赣战绝钠赡碱礼垦症抛席腊臣窄斋医学图像的配准与融合医学图像的配准与融合2. 2. 薄板样条变换薄板样条变换 其中：其中：X X是坐标向量，是坐标向量， A A与与B B定义一个仿射变换，定义一个仿射变换， U U是径向基函数。是径向基函数。可以表示为仿射变
32、换与径向基函数的线性组合：可以表示为仿射变换与径向基函数的线性组合：二、空间变换二、空间变换非线性变换非线性变换浮远逆毕蛙迸幽店妄珐宪酥唉磋甫茁掠督员睫峦忽牌互熬汽逐甲搅杆然会医学图像的配准与融合医学图像的配准与融合在二维图像配准中：在二维图像配准中：在三维图像配准中：在三维图像配准中：二、空间变换二、空间变换非线性变换非线性变换蛇绽葛艺梢瓤鹃佰惜王茫累关肿瞩辈纷呼馒驯隆既过钢阶试拍汾菇篡踊擎医学图像的配准与融合医学图像的配准与融合三、参数的优化搜索三、参数的优化搜索配准的几何变换参数根据求解方式可分成两类：配准的几何变换参数根据求解方式可分成两类：一、根据获得的数据用联立方程组直接计算得到的
33、，这一一、根据获得的数据用联立方程组直接计算得到的，这一 类完全限制在基于特征信息的配准应用中。类完全限制在基于特征信息的配准应用中。二、根据参数空间的能量函数最优化搜索得到。在这一类二、根据参数空间的能量函数最优化搜索得到。在这一类 中所有的配准都变成一个能量函数的极值求解问题。中所有的配准都变成一个能量函数的极值求解问题。因此图像配准问题本质上是多参数优化问题，优化算法的因此图像配准问题本质上是多参数优化问题，优化算法的选择至关重要。选择至关重要。侵蠢韩湿攘毒慧骏应沥蝴凑遥爪将筐撞么孩伪梭谅甸砸怎都愁哨抡路蒂蜗医学图像的配准与融合医学图像的配准与融合常用的优化算法常用的优化算法: : Po
34、well Powell法、梯度下降法、法、梯度下降法、 遗传算法、模拟退火法、遗传算法、模拟退火法、 下山单纯形法、下山单纯形法、Levenberg-MarquadrtLevenberg-Marquadrt法等。法等。三、参数的优化搜索三、参数的优化搜索蒙店沫娥怪孪獭糠科古刮薯哭砂训话镭针胰渤另键这苑则刮背蘑暇狂谬阿医学图像的配准与融合医学图像的配准与融合( (一一) Powell) Powell法法PowellPowell法是一种传统的确定性优化方法，又称为方向加速法是一种传统的确定性优化方法，又称为方向加速法，由法，由M.J.D.PowellM.J.D.Powell于于19641964年首
35、先提出。年首先提出。基本含义是基本含义是: :对于对于n n维极值问题，首先沿着维极值问题，首先沿着n n个坐标方向求个坐标方向求极小，经多极小，经多n n次之后得到次之后得到n n个共轭方向，然后沿个共轭方向，然后沿n n个共轭方个共轭方向求极小，经过多次迭代后便可求得极小值。向求极小，经过多次迭代后便可求得极小值。三、参数的优化搜索三、参数的优化搜索换铂迟葵渐智多俞颓蝴组刀馋坝贝淌辖溯呻毒没叠躬血零揣挺碧讹卑丘宅医学图像的配准与融合医学图像的配准与融合PowellPowell法的原理：法的原理：对于某一问题，将其归结为求取某一目标函数对于某一问题，将其归结为求取某一目标函数的极小值。其中的
36、极小值。其中Y Y为一个向量：为一个向量： 。设置一个满秩的步长矩阵：设置一个满秩的步长矩阵：三、参数的优化搜索三、参数的优化搜索窖锐锥产郭勉虱作馋纯齿替欢渺汀梯剑输谬闯石鸯强校穴赤挖搬序氰透阜医学图像的配准与融合医学图像的配准与融合对于某一初始值对于某一初始值, ,迭代过程如下：迭代过程如下： 1.首先在首先在方向上搜索，求方向上搜索，求，使，使为极小，极小，并令：并令：2.依次求依次求，使，使为极小，并令极小，并令，如此下去，如此下去.3.最后求最后求，使，使为极小，并令极小，并令，4.令令,在新的在新的方向上在搜索一方向上在搜索一 次，即求次，即求,使使为极小，并令新的极小，并令新的为至
37、此，完成了第一轮至此，完成了第一轮n+1次的搜索。接下去进行下一轮的搜索，次的搜索。接下去进行下一轮的搜索，直至性能指标满意或满足某种停止条件为止。直至性能指标满意或满足某种停止条件为止。 三、参数的优化搜索三、参数的优化搜索葛拒禄漏霄戍婿耳捞洲豫枣褒咱缆缕肇擎职圃取宇炸姻誊交徽桨遮赏攫耸医学图像的配准与融合医学图像的配准与融合( (二二) ) 梯度下降法梯度下降法 该算法在求最小化过程中直接利用梯度信息，沿着该算法在求最小化过程中直接利用梯度信息，沿着起始点梯度方向的反方向，求出最小值点，然后移动到起始点梯度方向的反方向，求出最小值点，然后移动到最小值点，再重复上面的过程，直到前后点的函数值
38、的最小值点，再重复上面的过程，直到前后点的函数值的差小于给定的误差值，则结束迭代过程。差小于给定的误差值，则结束迭代过程。三、参数的优化搜索三、参数的优化搜索潞椎忌追匀蝉定绪础蚤鞍磁呵厦浓疯忌巨齿腔狮狂纳醉香久征良顿嘱哨填医学图像的配准与融合医学图像的配准与融合( (三三) )遗传算法遗传算法 遗传算法（遗传算法（Genetic AlgorithmGenetic Algorithm）是模拟达尔文的遗）是模拟达尔文的遗传选择和自然淘汰的生物进化过程的计算模型，是一种传选择和自然淘汰的生物进化过程的计算模型，是一种通过模拟自然进化过程搜索最优解的方法，它是由美国通过模拟自然进化过程搜索最优解的方法
39、，它是由美国MichiganMichigan大学大学J.HollandJ.Holland教授于教授于19751975年首先提出来的。年首先提出来的。三、参数的优化搜索三、参数的优化搜索辜擎摊淹垒殆瞎宇砾骑天亚皑奈枢嗓驹却压递殷旨衔敏摊免濒词卤洼陕戌医学图像的配准与融合医学图像的配准与融合( (三三) )遗传算法遗传算法 在求解优化问题时，遗传算法将优化问题当作一个生存在求解优化问题时，遗传算法将优化问题当作一个生存环境，问题的一个解当作生存环境中的一个个体，以目标函环境，问题的一个解当作生存环境中的一个个体，以目标函数值或其变化形式来评价个体对环境的适应能力，模拟由数值或其变化形式来评价个体对
40、环境的适应能力，模拟由一定数量个体所组成的群体的进化过程，优胜劣汰，最终一定数量个体所组成的群体的进化过程，优胜劣汰，最终获得最好的个体，即问题的最优解。获得最好的个体，即问题的最优解。三、参数的优化搜索三、参数的优化搜索拽品舌桩削诛矣妊芹窝氟姿帘瓦判把渤间愿工读蝇院区频玻羌皮畅刑站燃医学图像的配准与融合医学图像的配准与融合四、插值方法四、插值方法 在图像配准中，空间坐标变换后得到的像素坐标位置在图像配准中，空间坐标变换后得到的像素坐标位置可能不在整数像素上，因此需要用灰度插值的方法对像素可能不在整数像素上，因此需要用灰度插值的方法对像素值进行估计。值进行估计。 常用的插值方法有常用的插值方法
41、有: :最近邻插值法、双线性插值法和部最近邻插值法、双线性插值法和部分体积分布法等。分体积分布法等。允无刽强垮斤豺票庆尘笼贞晒深麻瘩藤肝镍提孺沈真状袍噎妇样宁丙洪屹医学图像的配准与融合医学图像的配准与融合(1)(1)最近邻插值（最近邻插值（NNNN）图8-4 最近邻插值示意图计算计算n n和邻近四个点之间的距离，并将与该点距离最小的点和邻近四个点之间的距离，并将与该点距离最小的点的灰度值赋给的灰度值赋给n n。 认此狮蔬辞算涧阳蝗活渣汕魏盅枣琉谣顾梢扑殖膘犬仁懂祟盲梯意伺窗温医学图像的配准与融合医学图像的配准与融合(2)(2)双线性插值（双线性插值（BIBI）双线性插值法又称为一阶插值算法，它
42、是用线性插值来求双线性插值法又称为一阶插值算法，它是用线性插值来求像素灰度的一种方法。具体计算方法为先沿着一个坐标轴像素灰度的一种方法。具体计算方法为先沿着一个坐标轴方向使用线性插值方法求出两点的插值灰度，然后沿另一方向使用线性插值方法求出两点的插值灰度，然后沿另一个坐标轴，利用这两个点对目标点进行线性插值来求灰度个坐标轴，利用这两个点对目标点进行线性插值来求灰度。 灰度值灰度值, , 各相邻点的权重，与它们到各相邻点的权重，与它们到n n的距离成反比。的距离成反比。 醉鸭心洪绎逾诊羡恃捐操绢劫萨驯诫吝需季汝遇州二刷渊特小荒谩盏慰这医学图像的配准与融合医学图像的配准与融合(2)(2)双线性插值
43、（双线性插值（BIBI）茅汉挨悸又爷摄柿宰辩络灼舅杭酒罢趾定腿膀温朔诞碗劫疾汲规答哉息淄医学图像的配准与融合医学图像的配准与融合( (四四) )部分体积插值法（部分体积插值法（PVPV） 部分体积分布法是部分体积分布法是F.MaesF.Maes等人提出来的，是对双线性等人提出来的，是对双线性插值方法的一个改进。主要是为了克服双线性插值方法在图插值方法的一个改进。主要是为了克服双线性插值方法在图像中会产生新的灰度值而引起图像灰度分布发生变化的缺像中会产生新的灰度值而引起图像灰度分布发生变化的缺点，以便得到比较光滑的目标函数，有利于优化搜索。点，以便得到比较光滑的目标函数，有利于优化搜索。凭访曝韧
44、瓢砖叉环掏骤愁洪鼻拄乃崇境祟置硅嘴犬腑四七揖捕厘睬疆儒脂医学图像的配准与融合医学图像的配准与融合 PV PV根据线性插值的权重分配原则，将每对像素对联合直根据线性插值的权重分配原则，将每对像素对联合直方图的贡献分散到联合直方图中与之相邻的各个像素对上，方图的贡献分散到联合直方图中与之相邻的各个像素对上，这样联合直方图上各个像素对的频度值以小数增加，因此这样联合直方图上各个像素对的频度值以小数增加，因此不会出现新的灰度值而破坏目标函数值分布的光滑性。不会出现新的灰度值而破坏目标函数值分布的光滑性。具体的计算公式为具体的计算公式为: :为权重，其取值同为权重，其取值同BIBI法。法。 ( (四四)
45、 )部分体积插值法（部分体积插值法（PVPV） 卜婉碧涣共味脯契藤今乳水销斡豹倒酉阐里僵啃丑敛砒焉愧开慕曝博随辑医学图像的配准与融合医学图像的配准与融合五、五、相似性测度相似性测度 配配准准过过程程在在得得到到几几何何变变换换后后，进进一一步步的的工工作作就就是是要要找找到到一一种种合合适适、最最优优的的描描述述量量，用用以以表表征征相相似似或或者者差差异异，称称这种描述量为相似性测度。这种描述量为相似性测度。 毙姐操敏悍漏遂速溜藤爽译族理畏垒氟样杨臃找委缮辱油署便流稍突膏担医学图像的配准与融合医学图像的配准与融合（一）灰度均方差（一）灰度均方差设设 和和 分别表示参考图像和浮动图像中的数据，
46、分别表示参考图像和浮动图像中的数据，两幅图像像素值的均方差可以表示为两幅图像像素值的均方差可以表示为: :其中，其中，V V表示参与计算的图像区域，表示参与计算的图像区域，表示参与计算的像素总量，表示参与计算的像素总量，表示对图像数据的变换。表示对图像数据的变换。灰度均方差作为相似性测度适用于单模医学图像的配准。灰度均方差作为相似性测度适用于单模医学图像的配准。 杨隅骤钨痔石句对遁占普任啮冀迟惫冷检篷檄苔幼奔坯勘凑物禄盖功滤拿医学图像的配准与融合医学图像的配准与融合（二）（二） 归一化互相关归一化互相关分别是参考图和浮动图在分别是参考图和浮动图在区域内的区域内的像素灰度平均值，像素灰度平均值，
47、R R为相关系数。为相关系数。燎无棵熬央蒲剩链灸憎猛恼亥乾儒撕埔遁捧泞鹅铁葫瘸证台抛愈映苛萧双医学图像的配准与融合医学图像的配准与融合（三）（三） 互信息互信息 相似测度形式多样，在基于全图像信息的图像配准中，相似测度形式多样，在基于全图像信息的图像配准中，以互信息量作为相似测度的方法以其计算复杂度低、鲁棒以互信息量作为相似测度的方法以其计算复杂度低、鲁棒性好等特性逐渐成为当前研究的热点。性好等特性逐渐成为当前研究的热点。下一节对基于互信息量的图像配准方法做详细的介绍。下一节对基于互信息量的图像配准方法做详细的介绍。唇汁异侍氧梧卧豢聊蹲辟忘纱土牟购簇腰姿桨痉管械妨缨郧瞳逞讫眠怖诌医学图像的配准
48、与融合医学图像的配准与融合8.4 8.4 图像配准的主要方法图像配准的主要方法 图像配准的方法有多种，目前主要的配准方法大体上可图像配准的方法有多种，目前主要的配准方法大体上可以分为两类以分为两类: : 基于特征的配准方法基于特征的配准方法 基于灰度的配准方法。基于灰度的配准方法。 拴痴粤襄衅狗酥咸佩蛙阿禾舷艘牲肌班在甘蚌包祖肆唯捡此恭倪眼纹使筹医学图像的配准与融合医学图像的配准与融合一、基于特征的配准方法一、基于特征的配准方法 配准过程：配准过程： 首首先先对对待待配配准准图图像像进进行行特特征征提提取取，常常用用到到的的图图像像特特征征有有: :点点、直直线线段段、边边缘缘、轮轮廓廓、闭闭
49、合合区区域域、特特征征结结构构以以及及统统计特征如矩不变量、重心等等。计特征如矩不变量、重心等等。 然后利用提取到的特征完成两幅图像特征之间的匹配。然后利用提取到的特征完成两幅图像特征之间的匹配。惰拔拖粥库冤要谬液毕尺蓬术啡莹潜崭炔沈硒柏奋儒队又轰后莆副叙荚慰医学图像的配准与融合医学图像的配准与融合（一）、基于点特征的配准（一）、基于点特征的配准 点点特特征征是是图图像像配配准准中中最最为为常常用用的的图图像像特特征征之之一一，分分为为外部特征点与内部特征点两种。外部特征点与内部特征点两种。 外外部部特特征征点点：是是成成像像时时固固定定在在患患者者身身体体上上的的标标记记物物。这种方法的配准
50、变换被限制为刚性变换。这种方法的配准变换被限制为刚性变换。 侵入性标记物侵入性标记物 非侵入性标记物非侵入性标记物 丙兽诣瘁同很疽痰贡嚎磋兑媚项凶侈汤青隙溜愈眼肿誉冯稿宴突袭吠只野医学图像的配准与融合医学图像的配准与融合（一）基于点特征的配准（一）基于点特征的配准 内内部部特特征征点点：是是一一些些有有限限的的可可明明显显识识别别的的点点集集，可可以以是是解解剖剖点点( (一一般般由由用用户户识识别别出出) )，也也可可以以是是几几何何点点( (包包括括边边缘缘点点、角角点点、灰灰度度的的极极值值点点、曲曲率率的的极极值值点点、两两个个线线性性结结构构的的交点或某一封闭区域的质心等交点或某一封
51、闭区域的质心等) )。 这这种种方方法法主主要要求求解解刚刚体体或或仿仿射射变变换换，如如果果标标记记点点数数目目足足 够多，也能用来更复杂的非刚体变换。够多，也能用来更复杂的非刚体变换。 必裕眨立腿怜栽启撵散韭宣染愚丙亿翅敲蜀舶敬灼贾痉苹职神嗅销累店株医学图像的配准与融合医学图像的配准与融合（二）基于直线特征的配准（二）基于直线特征的配准 线段是图像中另一个易于提取的特征。一般利用线段是图像中另一个易于提取的特征。一般利用HoughHough变变 换提取图像中的直线。换提取图像中的直线。 建立两幅图像中分别提取的直线段的对应关系。建立两幅图像中分别提取的直线段的对应关系。 利用直线段的斜率和
52、端点的位置关系，可以构造一个这些利用直线段的斜率和端点的位置关系，可以构造一个这些信息指标的直方图，并通过寻找直方图的聚集束达到直线信息指标的直方图，并通过寻找直方图的聚集束达到直线段的匹配。段的匹配。祈捕执馆吞莎灸嗓辙抽走值遇倾辽弊钩已漠幻部啦皋计寄讼顺旭疡溅裳妇医学图像的配准与融合医学图像的配准与融合（三）基于轮廓与曲线特征的配准（三）基于轮廓与曲线特征的配准 近年来，随着图像分割、边缘检测等技术的发展，近年来，随着图像分割、边缘检测等技术的发展，基于边缘、轮廓的图像配准方法逐渐成为配准领域的基于边缘、轮廓的图像配准方法逐渐成为配准领域的研究热点。研究热点。 分割和边缘检测技术是这类方法的
53、基础，目前已分割和边缘检测技术是这类方法的基础，目前已报道的有很多图像分割方法可以用来做图像配准需要报道的有很多图像分割方法可以用来做图像配准需要的边缘轮廓和区域的检测，比如的边缘轮廓和区域的检测，比如CannyCanny边缘提取算子，边缘提取算子，拉普拉斯拉普拉斯- -高斯算子高斯算子(LoG )(LoG )，动态阈值技术，区域增长等。，动态阈值技术，区域增长等。璃影沫澎句剔魂岸遣困斡应愧兽系陀传彤委颖洋待卯痞熏率串契佛蜗卖爱医学图像的配准与融合医学图像的配准与融合（三）基于轮廓与曲线特征的配准（三）基于轮廓与曲线特征的配准在特征提取的基础上，很多学者针对轮廓、边缘等进行了配在特征提取的基础
54、上，很多学者针对轮廓、边缘等进行了配准研究。准研究。 1 1、GovinduGovindu等采用轮廓上点的切线斜率来表示物体轮廓，等采用轮廓上点的切线斜率来表示物体轮廓， 通过比较轮廓边缘的分布确定变换参数。通过比较轮廓边缘的分布确定变换参数。2 2、DavatzikosDavatzikos等提出了一种二阶段大脑图像配准算法，等提出了一种二阶段大脑图像配准算法， 在第一阶段使用活动轮廓算法建立一一影射，第二阶段在第一阶段使用活动轮廓算法建立一一影射，第二阶段 采用弹性变换函数确定轮廓的最佳变换。采用弹性变换函数确定轮廓的最佳变换。氏票炭晋病戈惊瞪捌潦侯毗目傀综矗恢挚省芍忘惧咏低学蛇友筐种逝绽制
55、医学图像的配准与融合医学图像的配准与融合（三）基于轮廓与曲线特征的配准（三）基于轮廓与曲线特征的配准 3 3、李登高等提出了一种对部分重叠的图像进行快速配准的、李登高等提出了一种对部分重叠的图像进行快速配准的 方法，该方法是基于轮廓特征的随机匹配算法。通过提方法，该方法是基于轮廓特征的随机匹配算法。通过提 取轮廓上的取轮廓上的“关键点关键点”作为特征点，随机选择若干特征点作为特征点，随机选择若干特征点 对得到候选变换，随后的投票阶段对其变换参数进行检对得到候选变换，随后的投票阶段对其变换参数进行检 验和求精。验和求精。 4 4、赵训坡等提出一种基于证据积累的图像曲线粗匹配方法，、赵训坡等提出一
56、种基于证据积累的图像曲线粗匹配方法， 比较有效地解决了将图像中提取的一条曲线比较有效地解决了将图像中提取的一条曲线( (较短较短) )与一与一 条参考曲线条参考曲线( (较长较长) )相匹配的问题。相匹配的问题。映撵淀微误碳蚂炒而恩吧甸东京椒善苞砚辛舰慧浴歇刊绿羹顶刻卞散玲恃医学图像的配准与融合医学图像的配准与融合（四）基于面特征的配准（四）基于面特征的配准基于面的配准方法中最典型的算法是由基于面的配准方法中最典型的算法是由PelizzariPelizzari和和 Chen Chen提出的提出的“头帽法头帽法”（Head-Hat MethodHead-Hat Method）。）。从一幅图像中提
57、取一个表面模型称为从一幅图像中提取一个表面模型称为“头头”(Head)”(Head)，从另外一幅图，从另外一幅图像轮廓上提取的点集称为像轮廓上提取的点集称为“帽子帽子”(Hat)”(Hat)。用刚体变换或选择性的仿。用刚体变换或选择性的仿射变换将射变换将“帽子帽子”的点集变换到的点集变换到“头头”上，然后采用优化算法使得上，然后采用优化算法使得“帽帽子子”的各点到的各点到“头头”表面的均方根距离最小。表面的均方根距离最小。菇谎啄李钱扰怒枝偷若稗国敲失炭期峙闲序锭财坤笼睁沃欢碴米忱遮屠春医学图像的配准与融合医学图像的配准与融合（四）基于面特征的配准（四）基于面特征的配准 头帽法最初用于头部的头帽
58、法最初用于头部的SPECTSPECT和和CT (CT (或或MRI)MRI)配准，参考配准，参考特征是头部的皮肤表面；然后用于头部的特征是头部的皮肤表面；然后用于头部的SPECTSPECT图像之间的图像之间的配准，参考特征是头颅骨表面和大脑表面。配准，参考特征是头颅骨表面和大脑表面。 优化算法目前一般用优化算法目前一般用PowellPowell法。均方距离是六个待求法。均方距离是六个待求刚体变换参数的函数，其最小时可得刚体变换参数。刚体变换参数的函数，其最小时可得刚体变换参数。 啤廉舅贵闭衍钩呢伶集绷票哟廖照煎划氯齐瓶记沸烛紧辛发柄苟鹿寒伴障医学图像的配准与融合医学图像的配准与融合（四）基于面
59、特征的配准（四）基于面特征的配准 比较常用的配准方法还有迭代最近点算法比较常用的配准方法还有迭代最近点算法(ICP)(ICP)。迭代最近点（。迭代最近点（ICPICP）配准算法由配准算法由BeslBesl和和MckayMckay提出的，它将一般的非线性最小化问题归结为基提出的，它将一般的非线性最小化问题归结为基于点的迭代配准问题。于点的迭代配准问题。 迭代最近点算法迭代最近点算法(ICP)(ICP)中必须先采样出图像结构上的特征点，然后用中必须先采样出图像结构上的特征点，然后用迭代的方法不断求出一幅图中相对于另一幅图中所有采样点的最近点，直迭代的方法不断求出一幅图中相对于另一幅图中所有采样点的
60、最近点，直到两个点集的均方差低于设定阈值，这时可得到匹配变换参数。到两个点集的均方差低于设定阈值，这时可得到匹配变换参数。躬姜酮避洱衍顶福以首窘死苑纸湃徽箩吼欧朵荚靴簇睬倍和疵带馈肺靛磋医学图像的配准与融合医学图像的配准与融合二、基于灰度的配准方法二、基于灰度的配准方法n基于灰度的配准方法是目前研究得较多的一种方法，它直接利用图像基于灰度的配准方法是目前研究得较多的一种方法，它直接利用图像的灰度信息进行配准，从而避免了因分割带来的误差，因而具有精度的灰度信息进行配准，从而避免了因分割带来的误差，因而具有精度较高、稳健性强、不需要预处理而能实现自动配准的特点。较高、稳健性强、不需要预处理而能实现
61、自动配准的特点。n 基于灰度的配准有两类主要的方法，一类是通过图像灰度直接计算出基于灰度的配准有两类主要的方法，一类是通过图像灰度直接计算出代表性的比例和方向等要素；另一类是配准过程中使用全部的灰度信代表性的比例和方向等要素；另一类是配准过程中使用全部的灰度信息。息。n第一种方法以力矩和主轴法为代表，第二种方法一般称为体素相似性。第一种方法以力矩和主轴法为代表，第二种方法一般称为体素相似性。嘶锻厦密窥侯罗缚氛航殿漏瓶裂侧秦哩攫雷冗洼朝醇络序贾码锐穴蕴正倾医学图像的配准与融合医学图像的配准与融合（一）（一） 力矩和主轴法力矩和主轴法 力矩和主轴法是指先用经典力学物体质量分布的原理计力矩和主轴法是
62、指先用经典力学物体质量分布的原理计算出两幅图像的质心和主轴，再通过平移和旋转变换使两幅算出两幅图像的质心和主轴，再通过平移和旋转变换使两幅图像达到配准。利用此方法，图像可以模型化为椭圆形区域图像达到配准。利用此方法，图像可以模型化为椭圆形区域的点分布。这样的分布可以用这些点的位置的一阶和二阶矩的点分布。这样的分布可以用这些点的位置的一阶和二阶矩描述。描述。 邯贩搞蘸呛棉触客营恫魔闰挖尼燥铂专批呸修贸绷理氖庆挛慢床每呢饭捐医学图像的配准与融合医学图像的配准与融合（一）（一） 力矩和主轴法力矩和主轴法缺点：缺点：该方法对数据的缺失较敏感，要求整个物体必须完整该方法对数据的缺失较敏感，要求整个物体必
63、须完整 地出现在两幅图像中。从整体上来说，配准精度较差，地出现在两幅图像中。从整体上来说，配准精度较差， 所以目前它更多地用来进行粗配准，使两幅图像初步所以目前它更多地用来进行粗配准，使两幅图像初步 对齐，以减少后续主要配准方法的搜索步骤。对齐，以减少后续主要配准方法的搜索步骤。距颗宦再遮菩图坠裔滩尽崔滦袄睡丸才仙渊率帚神公狄耻确浩矫棋蚤街噪医学图像的配准与融合医学图像的配准与融合（二）（二） 体素相似性法：体素相似性法： 体素相似性法是目前研究较多的一类方法。它是利用体素相似性法是目前研究较多的一类方法。它是利用图像中的所有灰度信息，这种方法一般都较为稳定，并能图像中的所有灰度信息，这种方法
64、一般都较为稳定，并能获得相当准确的结果。获得相当准确的结果。 该方法是完全自动的，不需要特殊的预处理，但这种该方法是完全自动的，不需要特殊的预处理，但这种方法由于需要大量的复杂计算，因此最近几年才转入实际方法由于需要大量的复杂计算，因此最近几年才转入实际应用。应用。翼晾恼娟负惯亢撩鞠工葡眠踩映垛撑巷着孰坝活褒哆缀辉鲤襟姑缄哗税纫医学图像的配准与融合医学图像的配准与融合（二）（二） 体素相似性法：体素相似性法：常见的基于体素相似性的配准方法有：常见的基于体素相似性的配准方法有： 互相关法互相关法 基于傅立叶域的互相关法和相位相关法基于傅立叶域的互相关法和相位相关法 灰度比的方差最小化法灰度比的方
65、差最小化法 直方图的互信息最大化法等。直方图的互信息最大化法等。主要讨论主要讨论。佛佑弧套免五眯呐霖盂桓蝉融丑镍椰赶芯鲜凹杨据葱禹霍嫂湖锅宵锗灶苍医学图像的配准与融合医学图像的配准与融合 最大互信息法最大互信息法 最大互信息法以互信息作为相似性测度。最大互信息法以互信息作为相似性测度。 19951995年分别被年分别被Viola Viola 和和CollignonCollignon等首次用于医学图像配准中。等首次用于医学图像配准中。 互信息互信息 (Mutual Information (Mutual Information， MI) MI)是信息论中的一个基是信息论中的一个基本概念，用于描述
66、两个系统间的统计相关性，或者是在一个系本概念，用于描述两个系统间的统计相关性，或者是在一个系统中包含的另一个系统的信息的多少，一般用熵来表示，表达统中包含的另一个系统的信息的多少，一般用熵来表示，表达的是一个系统的复杂性或不确定性。的是一个系统的复杂性或不确定性。演凝瞬叁琴瓮仪溪袒芹婚训后捻目浚晤形遵穿晦眠歌才佛哲故害辖添孰霉医学图像的配准与融合医学图像的配准与融合对于概率分布函数为对于概率分布函数为 的随机变量集的随机变量集A A，其熵：，其熵： 表示已知系统表示已知系统B B时时A A的条件熵。的条件熵。 最大互信息法最大互信息法规胎冉搭缎桔徊孜循冗挎锌歧肃肪料范廖薄萧同洛崔馅即仔讽游漳民
67、仁绅医学图像的配准与融合医学图像的配准与融合 在医学图像配准中，虽然两幅图像来源于不同的成像设在医学图像配准中，虽然两幅图像来源于不同的成像设备，但是它们基于共同的人体解剖信息，所以当两幅图像的备，但是它们基于共同的人体解剖信息，所以当两幅图像的空间位置达到完全一致时，其中一幅图像表达的关于另一幅空间位置达到完全一致时，其中一幅图像表达的关于另一幅图像的信息，也就是对应像素灰度的互信息应为最大。图像的信息，也就是对应像素灰度的互信息应为最大。 通常用联合概率分布和完全独立时的概率分布间的广通常用联合概率分布和完全独立时的概率分布间的广义距离来估计互信息：义距离来估计互信息： 最大互信息法最大互
68、信息法诧隶托骂句展酣个吱军锦济晋胚镊患蜜珊镐异桶匪役蛤蹲佑尺谋抹屉朱忌医学图像的配准与融合医学图像的配准与融合 对对于于离离散散的的数数字字图图像像，联联合合概概率率分分布布 可可以以用用归归一一化的联合直方图表示：化的联合直方图表示：边缘概率分布边缘概率分布表示为：表示为： 最大互信息法最大互信息法赁陷龙谈冈拱恬缅袜镍辩疙赐维哀样溢暗炎链榆家谅罗屈布荒库聋批缸汞医学图像的配准与融合医学图像的配准与融合接下来寻找一个变换使得一幅图像经过此变换后和另一幅接下来寻找一个变换使得一幅图像经过此变换后和另一幅图像的互信息最大。一般采用刚体变换，即在三维空间中图像的互信息最大。一般采用刚体变换，即在三维
69、空间中寻找三个方向上的平移值和旋转角度。对于大规模断层扫寻找三个方向上的平移值和旋转角度。对于大规模断层扫描医学图像来说，三维体积数据集包含的数据量极大无法描医学图像来说，三维体积数据集包含的数据量极大无法满足临床上实时处理的要求，因此必须采取优化措施。满足临床上实时处理的要求，因此必须采取优化措施。常用无需计算梯度的常用无需计算梯度的PowellPowell多参数优化算法。多参数优化算法。 最大互信息法最大互信息法碳键押影辜睬玲挎癌扭霸晌揪向蔽凹蜡协参锌初病箭虹庆斋膊蛋暑纫芒佃医学图像的配准与融合医学图像的配准与融合8.5 8.5 图像配准的评估图像配准的评估 医学图像配准，特别是多模医学图
70、像配准结果的评估一医学图像配准，特别是多模医学图像配准结果的评估一直是件很困难的事情。由于待配准的多幅图像基本上都是直是件很困难的事情。由于待配准的多幅图像基本上都是在不同时间或在不同时间或/ /和条件下获取的，所以没有绝对的配准问和条件下获取的，所以没有绝对的配准问题，即不存在什么金标准（题，即不存在什么金标准（gold standardgold standard）。只有相对）。只有相对的最优（某种准则下的）配准。的最优（某种准则下的）配准。 常用的评估方法有以下几种：常用的评估方法有以下几种：存顿媒健莫眯闽给洼爬畜驭水痔摈忍旭朋普沧腹副咖准拄赤笆牡血岳襟肿医学图像的配准与融合医学图像的配准
71、与融合一、体模（一、体模（PhantomPhantom） 体模又有硬件体模和软件体模之分。后者是计算机图体模又有硬件体模和软件体模之分。后者是计算机图像合成的结果。像合成的结果。 体模法用已知的图像信息验证新配准算法的精度。体模法用已知的图像信息验证新配准算法的精度。 由于体模都比较简单，与实际临床图像差异较大，因由于体模都比较简单，与实际临床图像差异较大，因此只能对配准方法作初步的评估此只能对配准方法作初步的评估。芒塑靠境曲凑叔售颠茂样旅揪蓖射沽慧四柯暂捧偿晋弄良忿宽斡瓶户呼乞医学图像的配准与融合医学图像的配准与融合二、准标（二、准标（Fiducial MarksFiducial Marks
72、） 立体定向框架系统（立体定向框架系统（Stereotactic Frame SystemsStereotactic Frame Systems）包括立体定向参考框架、立体定向图像获取、探针或手术包括立体定向参考框架、立体定向图像获取、探针或手术器械导向几部分。优点是定位准确，不易产生图像畸变。器械导向几部分。优点是定位准确，不易产生图像畸变。使用立体定向框架系统的体积图像数据可以用来评估其它使用立体定向框架系统的体积图像数据可以用来评估其它配准方法的精度。配准方法的精度。惕亥奔龟宛萨亚逼晰绷呢治爷笑绸拒蕉阶诵伯睬杂奉睡霍伎掖浆祥遏貌挑医学图像的配准与融合医学图像的配准与融合二、准标（二、准标
73、（Fiducial MarksFiducial Marks） 使用人工记号作准标的方法很多。一种准标是使用使用人工记号作准标的方法很多。一种准标是使用9 9根根棍棒组成的棍棒组成的3 3个方向的个方向的N N字型结构。在作字型结构。在作CTCT测试时，棒内充测试时，棒内充以硫酸铜溶液；作以硫酸铜溶液；作PETPET测试则填充氟测试则填充氟1818。这样，在两组图。这样，在两组图像中都可见此像中都可见此N N字型准标，从而可对图像准确空间定位。字型准标，从而可对图像准确空间定位。例如用在人脑表面嵌螺丝作标记（每人例如用在人脑表面嵌螺丝作标记（每人8 8个）的方法对多个）的方法对多个病人做个病人做
74、CTCT、MRMR（T1T1、T2T2及及PDPD）和）和PETPET实测，得到多组数实测，得到多组数据。这些数据专门用于多模医学图像配准算法评估使用。据。这些数据专门用于多模医学图像配准算法评估使用。漫芹菜坦郁华松讥罗尾钮苇佣哮桨黍蹋渤恭仇伐聋苑谅卢令姥饺甘惋书珍医学图像的配准与融合医学图像的配准与融合三、三、 图谱图谱 （AtlasAtlas） ThompsonThompson用随机向量场变换构造一个可变形的概率脑图用随机向量场变换构造一个可变形的概率脑图谱。包括从多个受试者到单一解剖模板的功能、血管、组谱。包括从多个受试者到单一解剖模板的功能、血管、组织诸方面映射，三维图谱到新受试者的扫
75、描图像的映射。织诸方面映射，三维图谱到新受试者的扫描图像的映射。Visible Human CDVisible Human CD的的CTCT骨窗图像、骨窗图像、MRMR图像及彩绘的冷冻切图像及彩绘的冷冻切片照片由于具有清晰的解剖结构和高度的分辨（片照片由于具有清晰的解剖结构和高度的分辨（1 1毫米毫米/ /每每层片），近来也被用来做新配准方法精度的评估。层片），近来也被用来做新配准方法精度的评估。谷哭精衅器跋软杠鲜麓宛脸溯诅几非淬拖央桌激盒酿星公泊间酌壮租耽尘医学图像的配准与融合医学图像的配准与融合四、目测检验（四、目测检验（Visual InspectionVisual Inspection
76、） 对多模医学图像配准的结果请领域专家用目测方法检验，对多模医学图像配准的结果请领域专家用目测方法检验，听起来有些主观，但在一定程度上的确是一种相当可信的听起来有些主观，但在一定程度上的确是一种相当可信的方法。方法。唾况冲晓陋鱼蛆孙阎瞥俘城土妮懈菇柱熙迈叫溅痛泻啸漠诱馏寸狞跺府反医学图像的配准与融合医学图像的配准与融合8.6 8.6 医学图像融合概述医学图像融合概述 医学图像融合：医学图像融合：是指将两幅是指将两幅( (或两幅以上或两幅以上) )来自不同成来自不同成像设备或不同时刻获取的已配准的图像，采用某种算像设备或不同时刻获取的已配准的图像，采用某种算法，把各个图像的优点或互补性有机地结合
77、起来，获法，把各个图像的优点或互补性有机地结合起来，获得信息量更丰富的新图像的技术。得信息量更丰富的新图像的技术。 先决条件：先决条件：医学图像的配准。医学图像的配准。女稀目醉希埃无箩费幻芦隧氛扛蝉祸纪愚古箕龚汀保徽痴叠泻簧曼轩斟硝医学图像的配准与融合医学图像的配准与融合医学图像融合的分类医学图像融合的分类n1 1按照融合图像成像方式的不同，分为单模按照融合图像成像方式的不同，分为单模融合融合(mono-modality)(mono-modality)和多模融合和多模融合(multi-(multi-modality)modality)。 单模：单模：CT-CTCT-CT ，MRI-MRIMRI
78、-MRI 。 多模：多模：CTCT 与与 MRIMRI等。等。 瞳碳戊蚂辜蹬距伪喉缎图偿诽饭抢蚕火圾肌虾戳我吊碴奋谁萍凳稗辰瘸刽医学图像的配准与融合医学图像的配准与融合n2 2按照融合对象的不同，分为单样本时间融合、单样按照融合对象的不同，分为单样本时间融合、单样本空间融合和模板融合。本空间融合和模板融合。 单样本时间融合：单样本时间融合：指跟踪某个病人，将其一段时指跟踪某个病人，将其一段时间内对同一脏器所做的同种检查图像进行融合。间内对同一脏器所做的同种检查图像进行融合。 单样本空间融合：单样本空间融合：指将某个病人在同一时期内对指将某个病人在同一时期内对同一脏器所做的几种检查的图像进行融合
79、。同一脏器所做的几种检查的图像进行融合。 模板融合：模板融合：从许多健康人的研究中建立一系列模从许多健康人的研究中建立一系列模板，将病人的图像与模板图像融合。板，将病人的图像与模板图像融合。 医学图像融合的分类医学图像融合的分类馆匆筐压棒科闹越狼铲树正田赃韵述大晃嫌柯仙脆渔仗碗映祷箍态兄毡狙医学图像的配准与融合医学图像的配准与融合n3 3按照图像处理方法的不同，分为数值融合法和智能按照图像处理方法的不同，分为数值融合法和智能融合法。融合法。 数值融合法：数值融合法：将不同来源的图像做空间归一将不同来源的图像做空间归一 化处理后直接融合。化处理后直接融合。 智能融合法：智能融合法：将不同来源的图
80、像做归一化处将不同来源的图像做归一化处 理后，根据需要选择不同图像理后，根据需要选择不同图像 中的所需信息再进行融合。中的所需信息再进行融合。 医学图像融合的分类医学图像融合的分类肝卉皇捏般抖泄锄运斋韧物甥梧填磊操燎录芳虎渡储乃彬饰狗早趟捉曙冻医学图像的配准与融合医学图像的配准与融合n4 4按图像类型不同，可以分为断层图像间相互融合、按图像类型不同，可以分为断层图像间相互融合、断层图像与投影图像融合以及结构图像与功能图像融断层图像与投影图像融合以及结构图像与功能图像融合。合。 断层图像间相互融合断层图像间相互融合主要指主要指 CT CT 与与 MRI MRI 图像融合；图像融合； 断层图像与投
81、影图像融合断层图像与投影图像融合主要指主要指 CT CT、MRI MRI 图像与图像与 DSA DSA 图像通过三维重建后进行融合；图像通过三维重建后进行融合； 结构图像融合与功能图像融合结构图像融合与功能图像融合主要指主要指 CT CT、MRI MRI 图像与图像与 PETPET、SPECT SPECT 图像进行融合。图像进行融合。 医学图像融合的分类医学图像融合的分类辉荡朽泌恨炎祭誊量绅而撅蒜恕获墒尹吨揖陕茎肛揽聊兔永肪绷肪墨滋记医学图像的配准与融合医学图像的配准与融合医学图像融合的分类医学图像融合的分类n5 5前瞻性融合和回溯性融合前瞻性融合和回溯性融合 前瞻性融合：前瞻性融合：在图像采
82、集时使用特别措施（如加外部在图像采集时使用特别措施（如加外部标志等）；标志等）； 回溯性融合：回溯性融合：在图像采集时则不采取特别措施。在图像采集时则不采取特别措施。 芭屠樊维思毛剃缆壕瑚谆乖枕吭幸掘天疲全焚辱肾痒善园疫蛊铃挡磕蔚及医学图像的配准与融合医学图像的配准与融合8.7 8.7 常用的图像融合方法常用的图像融合方法n一、基于空域的图像融合一、基于空域的图像融合 1 1 、图像像素灰度值极大（小）融合法；、图像像素灰度值极大（小）融合法； 2 2、图像像素灰度值加权融合法；、图像像素灰度值加权融合法； 3 3、TOETTOET图像融合方法。图像融合方法。n二、基于变换域的图像融合二、基于
83、变换域的图像融合 基于小波变换的图像融合。基于小波变换的图像融合。 在况概苇催嫉衫獭俘奢殊煽撅尤填雅镣状陪费冷耪车踌屡烽轻亩枝拍美澜医学图像的配准与融合医学图像的配准与融合1 1、图像像素灰度值极大（小）融合法、图像像素灰度值极大（小）融合法 设设g g1 1(i(i, ,j)j)和和g g2 2(i,j)(i,j)为为待待融融合合图图像像，F(i,j)F(i,j)为为融融合合后后的的图图像像，其其中中i i，j j为为图图像像中中某某一一像像素素的的坐坐标标，图图像像大大小小为为M*NM*N，则，则i0,M-1i0,M-1，j0,N-1j0,N-1。 极大值法：极大值法：极小值法：极小值法：
84、续媳鞠峦犀郭胆汞唯掷扫戌骤雇俄皑几彤瘤聋凿炙轮掺信惹咯舒矗纲湖箭医学图像的配准与融合医学图像的配准与融合2 2 、图像像素灰度值加权融合法、图像像素灰度值加权融合法 设设g g1 1(i(i, ,j)j)和和g g2 2(i,j)(i,j)为为待待融融合合图图像像，F(i,j)F(i,j)为为融融合合后后的的图像。图像。其中：其中：a a为权重因子，且为权重因子，且0a10a1，可以根据需要调，可以根据需要调节节a a的大小。该算法实现简单，其困难在于如何选的大小。该算法实现简单，其困难在于如何选择权重系数，才能达到最佳的视觉效果。择权重系数，才能达到最佳的视觉效果。朝甚姜恬升丽绰计寺赎汾说渠
85、鼻恬眺例窘毡配疏吓淑砷天怎窖衔虑剂步店医学图像的配准与融合医学图像的配准与融合3 3、TOETTOET图像融合方法图像融合方法 设设g g1 1(i(i, ,j)j)和和g g2 2(i,j)(i,j)为待融合图像，为待融合图像，F(i,j)F(i,j)为融合后的图像。为融合后的图像。 首先求输入图像和的共同成分首先求输入图像和的共同成分： 从图像上扣除共同成分得到图像的特征成分：从图像上扣除共同成分得到图像的特征成分： 从图像中扣除图像的特征成分：从图像中扣除图像的特征成分：恍矛住滦轩瓷呀蓬竟缩枕契蛤岩页擅取窗盒武狡称懒十蚁组雀烟召袁儿授医学图像的配准与融合医学图像的配准与融合 融合结果融合
86、结果CT CT 图像图像 MR MR 图像图像 像素灰度极小值法像素灰度极小值法像素灰度极大值法像素灰度极大值法蜡细察侥捅篓默蜀晕乓桐赐齿帅度浸蚁褂南司陀瞳丫竹响赊跨淳本肚锤缚医学图像的配准与融合医学图像的配准与融合 融合结果融合结果CT CT 图像图像 MR MR 图像图像 灰度加权法灰度加权法TOETTOET法法阉埂蝇饺洗铣尽设孙宦厢全况拇交肝角荒锄诌麓也嘱眨是毯逐女耀三拳承医学图像的配准与融合医学图像的配准与融合二、基于变换域的图像融合二、基于变换域的图像融合 变换域法，顾名思义，就是将变换后的两个或多变换域法，顾名思义，就是将变换后的两个或多个图像进行融合，再通过反变换得到融合后图像的
87、方个图像进行融合，再通过反变换得到融合后图像的方法。法。 多分辨率金字塔多分辨率金字塔小波变换法小波变换法 傅里叶变换法傅里叶变换法嘎剩娜斩幌丝惫瞒厂唱矩省卓露刀路押单段园佣氟浑刀邑铆陡扣宿脖帚肠医学图像的配准与融合医学图像的配准与融合一、图像的二维小波分解一、图像的二维小波分解MallatMallat算法：算法： 分别表示水平、垂直和对角分量；分别表示水平、垂直和对角分量； 和和 分别是分别是H H和和G G的共轭转置矩阵；的共轭转置矩阵；J 为分解层数。为分解层数。绘蓄纺岂命朴贪民糕拭革蚊掇罐寝破瓮脆棚治原捷私丸揣蕴换捞咨视调猴医学图像的配准与融合医学图像的配准与融合一、图像的二维小波分解
88、一、图像的二维小波分解图像经二维小波变换分解后，可得到四个不同的频带图像经二维小波变换分解后，可得到四个不同的频带LLLL、LHLH、HLHL、HHHH。其中低频带。其中低频带 LL LL 保留了原图的轮廓信息。保留了原图的轮廓信息。HLHL、LHLH、HH HH 分别保分别保留了原图水平、垂直和对角方向的高频信息，代表图像的细节部留了原图水平、垂直和对角方向的高频信息，代表图像的细节部分。分。 谆饮秃舒役幼腥滴治碰迄奈迎狞烬崩疵粪饼黄妄坷助绞汁杯潜咋文迟歪吏医学图像的配准与融合医学图像的配准与融合基于小波变换的图像融合步骤基于小波变换的图像融合步骤 分解：对每一源图像分别进行小波变换，得到每
89、幅图像在分解：对每一源图像分别进行小波变换，得到每幅图像在不同分辨率下不同频带上的小波系数；不同分辨率下不同频带上的小波系数； 融合：针对小波分解系数的特性，对各个不同分辨率融合：针对小波分解系数的特性，对各个不同分辨率 上的小波分解得到的频率分量采用不同的融合上的小波分解得到的频率分量采用不同的融合 方案和融合算子分别进行融合处理；方案和融合算子分别进行融合处理； 逆变换：对融合后系数进行小波逆变换，得到融合图像。逆变换：对融合后系数进行小波逆变换，得到融合图像。钦滤雇咙樊蹲挪坯楔俺嘱往镰舆逾搔菩纹皿肝衰账懊呢野线舜锗醒癌顷诚医学图像的配准与融合医学图像的配准与融合小波分解融合结构图小波分解
90、融合结构图图图8-78-7小波分解融合图小波分解融合图融合后的图像小波分解图小波分解图系数融合A图像小波变换F图像B图像小波变换逆变换毛哎镰柳右孽只咋删绣札棵屁讫耀怎咨瞩咎朵万席摔淘躯鹿卓切洼寐剖茶医学图像的配准与融合医学图像的配准与融合二、基于小波变换的融合规则二、基于小波变换的融合规则1 1、低频系数融合规则、低频系数融合规则 通过小波分解得到的低频系数都是正的变换值，反映通过小波分解得到的低频系数都是正的变换值，反映的是源图像在该分辨率上的概貌。低频小波系数的融合规的是源图像在该分辨率上的概貌。低频小波系数的融合规则有多种方法：既可以取源图像对应系数的均值，也可以则有多种方法：既可以取源
91、图像对应系数的均值，也可以取较大值，这要根据具体的图像和目的来定。取较大值，这要根据具体的图像和目的来定。熔招激贪订驻范室荆栗奢攒姨朵阻都希脓创褂航呆秸夹助仔北贷碳蜀育仿医学图像的配准与融合医学图像的配准与融合二、基于小波变换的融合规则二、基于小波变换的融合规则2 2、高高频系数融合规则频系数融合规则 通过小波分解得到的三个高频子带都包含了一些在零通过小波分解得到的三个高频子带都包含了一些在零附近的变换值，在这些子带中，较大的变换值对应着亮附近的变换值，在这些子带中，较大的变换值对应着亮度急剧变化的点，也就是图像中的显著特征点，如边缘、度急剧变化的点，也就是图像中的显著特征点，如边缘、亮线及区
92、域轮廓。这些细节信息，也反映了局部的视觉亮线及区域轮廓。这些细节信息，也反映了局部的视觉敏感对比度，应该进行特殊的选择。敏感对比度，应该进行特殊的选择。 浙拱速蜡序闪牛患递危雄痕灾卷泼驶寸州疆忻棚桶砰肋奥谭汞啪纶韩官壬医学图像的配准与融合医学图像的配准与融合二、基于小波变换的融合规则二、基于小波变换的融合规则图像A分解层图像B分解层基于像素的融合规则基于区域的融合规则基于窗口的融合规则图像F分解层融合规则系数加权绝对值选大区域能量最大系数模值极大。融合处理图8-8 小波融合规则嘎闺齿筑雅狞莹抄碘谐嗜丙奋蜂鱼调菇毡袍敲夸鲤手撩默忠锨戊己鞘粮船医学图像的配准与融合医学图像的配准与融合二、基于小波变
93、换的融合规则二、基于小波变换的融合规则（1 1）基于像素点的融合规则基于像素点的融合规则 逐个考虑源图像相应位置的小波系数，要求源图是经逐个考虑源图像相应位置的小波系数，要求源图是经过严格对准处理的。因为基于像素的选择方法具有其片过严格对准处理的。因为基于像素的选择方法具有其片面性，其融合效果有待改善。面性，其融合效果有待改善。 棺妥调茄词轩俺价割歌唾殿女液碑原燃千疡官攀懒穷汐贮观擦枯瘤炳簿名医学图像的配准与融合医学图像的配准与融合二、基于小波变换的融合规则二、基于小波变换的融合规则（2 2）基于窗口的融合规则基于窗口的融合规则 是对第一类方法的改进。由于相邻像素往往有相关性，是对第一类方法的
94、改进。由于相邻像素往往有相关性，该方法以像素点为中心，取一个该方法以像素点为中心，取一个M MN N的窗口，综合考虑的窗口，综合考虑区域特征来确定融合图像相应位置的小波系数。该类方区域特征来确定融合图像相应位置的小波系数。该类方法的融合效果好，但是也相应的增加了运算量和运算时法的融合效果好，但是也相应的增加了运算量和运算时间。由于窗口是一个矩形，是规则的；而实际上，图像间。由于窗口是一个矩形，是规则的；而实际上，图像中相似的像素点往往具有不规则性。中相似的像素点往往具有不规则性。 泥姐确狭琴辜逢广界鸡越彻扎渗蟹诺炙痒茅腻孵八穿逻肆山吧驳千陷闰怕医学图像的配准与融合医学图像的配准与融合二、基于小
95、波变换的融合规则二、基于小波变换的融合规则（3 3）基于区域的融合规则基于区域的融合规则 该类方法常常利用模糊聚类来寻找具有相似性的像素该类方法常常利用模糊聚类来寻找具有相似性的像素点集点集。 并芥纸腹棺印慧额柠雇欠敢口饺或淳巨廓武棠蹬采钮整渣递颧雾徐嚼蛛帜医学图像的配准与融合医学图像的配准与融合三、常用的小波分解系数融合规则三、常用的小波分解系数融合规则a a、小波系数加权法、小波系数加权法 其中：其中： 分别表示源图像分别表示源图像A A，B B和和融合图像融合图像F F在在J J层小波分解时，在层小波分解时，在P P点的系数。点的系数。轮邵功酬候匠烽卓嗣木掏思迂基棵览泪坐澜婉勇跨椭嗓小授
96、晦数衣蛙苇吏医学图像的配准与融合医学图像的配准与融合b b、小波系数、小波系数绝对值极大绝对值极大( (小小) )法法三、常用的小波分解系数融合规则三、常用的小波分解系数融合规则隆舷稀漳奢谅锄凭止疵屉遍炳触篇宠中诲檄诡不宗附增仙市挨粪症蚂郁幸医学图像的配准与融合医学图像的配准与融合c c、区域能量最大法、区域能量最大法在在J J层小波分解的情况下，局部区域层小波分解的情况下，局部区域Q Q的能量定义为：的能量定义为： 其中：表示权值，表示权值，点离点离点越近，权值越大。点越近，权值越大。是是的一个邻域。的一个邻域。 同理可得：同理可得： 三、常用的小波分解系数融合规则三、常用的小波分解系数融合
97、规则咖圆讨写钻辩挪襟虹成靴蜕钡必匹刺祸漳漆舵循掩斧婆血在椒至较滴诺戍医学图像的配准与融合医学图像的配准与融合小波融合举例小波融合举例 CT 图像 MR 图像 均值 均值均值、绝对值极大培雷曙娘敖忌星佩炒弓霖枉劫侮噬毡绿吃永果燥勉玩以骆锑舌雏韵恭叉蹲医学图像的配准与融合医学图像的配准与融合小波融合举例小波融合举例 CT 图像 MR 图像 均值、绝对值极小均值、区域能量最大粳至件谚旨关偿坟投鉴坦抨经善嘲影述节我雇欲塘料殷嫁金练耀遣植臼熟医学图像的配准与融合医学图像的配准与融合8.88.8、图像融合效果评价、图像融合效果评价 图像融合效果的评价主要有主观评价和客观评价两种。图像融合效果的评价主要有主
98、观评价和客观评价两种。主观评价以人作为观察者，对图像的优劣做出主观定性评主观评价以人作为观察者，对图像的优劣做出主观定性评价。人对图像的识别或理解不仅和图像的内容有关，而且价。人对图像的识别或理解不仅和图像的内容有关，而且还与观察者的心理状态有关。由于人的视觉系统很复杂，还与观察者的心理状态有关。由于人的视觉系统很复杂，受环境条件、视觉性能、人的情绪爱好以及知识状况影响受环境条件、视觉性能、人的情绪爱好以及知识状况影响很大，因此主观评价具有主观性和不全面性，很大，因此主观评价具有主观性和不全面性， 讨萝国甲锹氢脯给瘦退敌处召报培港闽榆戍遏侠炳驶档旱氧贸写虑塑远滨医学图像的配准与融合医学图像的配
99、准与融合图像融合效果的客观评价评价图像融合效果的客观评价评价一、熵一、熵 图像的熵值是衡量图像信息丰富程度的一个重要指标，图像的熵值是衡量图像信息丰富程度的一个重要指标，熵值的大小表示图像所包含的平均信息量的多少。根据熵值的大小表示图像所包含的平均信息量的多少。根据香农信息论的原理，一幅图像的信息熵为香农信息论的原理，一幅图像的信息熵为： 融合图像的熵增大，表示融合图像的信息量增加，融合图像的熵增大，表示融合图像的信息量增加，融合图像所包含的信息就越丰富，融合质量越好。融合图像所包含的信息就越丰富，融合质量越好。迭附肛轰民熟矛肋马锌呀祷翅犯寡种滑燕拆搬斋磊拧屠埃藩冷拆样敢愧曼医学图像的配准与融
100、合医学图像的配准与融合二、交叉熵二、交叉熵CECE 也称相对熵，反映了两幅图像灰度分布信息的差异。设源也称相对熵，反映了两幅图像灰度分布信息的差异。设源图像和融合图像的直方图分别为图像和融合图像的直方图分别为pipi和和qiqi，交叉熵定义为：，交叉熵定义为：交叉熵越小，说明融合图像从源图像提取的信息量越多，交叉熵越小，说明融合图像从源图像提取的信息量越多，融合效果越好。融合效果越好。营牧讳融鞋惧右桐诛旱襟抵猖椰赏壶缮褪哪不琳患娇纸帆谱胃散狰换彝至医学图像的配准与融合医学图像的配准与融合三、交互信息量三、交互信息量MIMI交互信息量为两个变量之间相关性的量度，或一个变量包交互信息量为两个变量之
101、间相关性的量度，或一个变量包含另一个变量的信息量的量度。假设两幅源图像含另一个变量的信息量的量度。假设两幅源图像A A和和B B，将，将它们融合得到融合图像它们融合得到融合图像F F，F F与与A A、B B的交互信息量分别表示的交互信息量分别表示为为MIMIFA FA 和和MIMIFBFB：沂乘辖亲吸祥邱圭痹找逾镀箱弥搓历蛆瘤剐聂宁须买扩桃咕攀怜横狮校漫医学图像的配准与融合医学图像的配准与融合三、交互信息量三、交互信息量MIMI交互信息量的总和交互信息量的总和 ：交互信息量的值越大，表示融合图像从源图像中获取交互信息量的值越大，表示融合图像从源图像中获取的信息越丰富，融合效果越好。的信息越丰
102、富，融合效果越好。起忿揽挠堡磅幽狰黄易农慨没稍讶残彭冠碎逝般盒猎呼湾先鄂撂隶懂他塑医学图像的配准与融合医学图像的配准与融合四、图像均值四、图像均值图像均值是图像像素的灰度平均值，对人眼反映为平均亮图像均值是图像像素的灰度平均值，对人眼反映为平均亮度。图像均值的定义为：度。图像均值的定义为：其中其中G G（x x，y y）表示图像中第（）表示图像中第（x x，y y）个像素的灰度，）个像素的灰度，图像尺寸为图像尺寸为。如果均值适中，则目视效果良好。如果均值适中，则目视效果良好。中鲸道挂屯骑阳贴扮蝎赊杰捐盾纫胎败霹凡颓煤酸节烧蕉跌蹲熬跺膛囤源医学图像的配准与融合医学图像的配准与融合五、灰度标准差五
103、、灰度标准差g g图像的灰度标准差定义为为：图像的灰度标准差定义为为： 其中其中L L为图像的总灰度级，为图像的总灰度级，g g表示图像第（表示图像第（x x，y y）个像素）个像素的灰度，的灰度， 表示图像均值，表示图像均值，p p( (g g) )表示灰度值为表示灰度值为g g的像素出现的的像素出现的概率。概率。 标准差大，则图像灰度级分布分散，图像的反差大，可以标准差大，则图像灰度级分布分散，图像的反差大，可以看出更多的信息。看出更多的信息。涎胁捐物芭陕洲肩促慌块邮傻纤侣跪就皂纤肛陆躁蛰拘咙褐揭小赡骗谆具医学图像的配准与融合医学图像的配准与融合融合图像融合图像F F和标准参考图像和标准参
104、考图像R R间的均方误差：间的均方误差：均方误差越小说明融合图像与标准参考图像越接近。均方误差越小说明融合图像与标准参考图像越接近。六、均方误差六、均方误差MSEMSE淤屋蛾均功噎刽翘盘催具舀校老藤徊佑降寞渔销寂褒内裤风耙业袜旬框杠医学图像的配准与融合医学图像的配准与融合融合图像与标准参考图像的差异看作噪声，标准参考融合图像与标准参考图像的差异看作噪声，标准参考图像看作信息。图像看作信息。七、信噪比七、信噪比SNRSNR与峰值信噪比与峰值信噪比PSNRPSNR信噪比、峰值信噪比越高，说明融合效果就越好。信噪比、峰值信噪比越高，说明融合效果就越好。吾衡颇掌桥悠教筹碾储烙堵双吗兼诀歪染皖炳坞盾商一渠韶痛阜啮疤元诧医学图像的配准与融合医学图像的配准与融合
