2021年认证杯SPSSPRO杯数学建模A题(第一阶段)医学图像的配准全过程文档及程序
2021年认证杯SPSSPRO杯数学建模
A题 医学图像的配准
原题再现:
图像的配准是图像处理领域中的一个典型问题和技术难点,其目的在于比较或融合同一对象在不同条件下获取的图像。例如为了更好地综合多种信息来辨识不同组织或病变,医生可能使用多种仪器对患者的同一部位进行成像。在综合多幅图像时,首先需要将它们严格对齐,使得图上同一个坐标的位置对应的是真实对象的同一个点,这个过程称之为配准。现在的许多医学成像技术,包括 CT、MRI、PET 等,最终生成的是人体的断层影像。在这里,我们主要关心的是断层成像的配准问题。
我们考虑对一个患者的腹部进行断层成像。由于人体组织是柔软的,所以即使使用同一台成像设备,两次成像的结果也并不完全一致。最终输出时还会对图像进行自动放缩,所以输出图片的大小也并不完全相同。想要精确配准,需要将其中一次的成像结果进行某种仿射变换(或非线性变换),以尽可能地匹配另一次的结果(或将两次结果都映射到同一个标准模板中)。求得合适的变换就是图像配准的核心任务。
第一阶段问题: 对同一个患者进行两次 CT 成像,间隔长达数周乃至更长。两次扫描的断层位置相同,但由于占位性病变的发展,成像的结果在某些区域会有区别。请你设计一个有效的方法,使我们能够对这样两张图像进行配准,以使计算机辅助医疗系统能够通过比较来自动识别出病变的位置,并定量地评估出病变的发展情况。
整体求解过程概述(摘要)
数字化医学影像技术的高速发展推动了图像配准技术的不断进步,目前刚性组织的医学图像配准技术已经较为成熟,而非刚性部位(如腹部)由于形状和位置会随生理运动而发生变化,因此其图像配准面临巨大的挑战。
本文是对同一患者腹部的两次 CT 图像进行配准分析,重点考虑占位性病变的影响因素,首先采用 Kullback-Leiber 距离对传统的 Demons 算法模型进行改进,以基于Kullback-Leiber 距离的互信息作为图像配准的精度量度,通过迭代次数以及互信息变化量的阈值来控制配准进程,从而实现收敛速度较快的非刚性医学图像配准;再通过对配准前后图像的相似性精度进行仿真计算,以及对模型性能的评估,验证了配准模型的准确度和有效性。对于配准后的图像,采用 CT 兴趣体积法测定 CT 图像中的病变区域的体积,同时通过单因素 ANOVA 分析法对病变区域的体积和临床严重程度的关系进行定量分析。
本文最后对于基于腹部(非刚性部位)的改进 Demons 形变模型进行评估,针对模型一些局限性,引入拓扑校正算法,以进一步提高配准的准确性,可以有效避免浮动图像产生拓扑结构无法保持的现象,从而使计算机辅助医疗系统更好地识别出病变的位置。
问题分析:
通过对题目背景的分析,我们知道由于同一个患者占位性病变的发展,并且腹部本身具有非刚性的一些特征,对图像的配准提出了更高的要求。对此,我们基于 Kullback-Leiber 距离对传统的 Demons 算法模型进行改进,实现收敛速度较快的腹部(非刚性)医学图像配准。在本模型中,我们拟用 Kullback-Leiber 距离定义法定义的互信息作为图像配准的测度,通过迭代次数和互信息变化量的阈值要求控制配准过程的结束。之后我们通过对十个案例的仿真,采用配准前后图像相似性测度(互信息与归一化互信息)变化情况,通过对模型性能进行评估,验证了模型的实用性。最后,在优化模块我们通过将拓扑校正算法引入配准模型,以提高配准的准确性,保证变形后的浮动图像不产生拓扑结构不能保持的现象,以便于计算机辅助医疗系统更好地识别出病变的位置。针对于问题的另一方面,即在医学图像配准后如何定量地评估出病变的发展情况的问题,我们选择采用 CT 兴趣体积法定量测定 CT 图像中病变区域体积,并采用单因素ANOVA 分析比较病变区域体积与临床严重程度关系。
模型假设:
(1)假设两次扫描使用的仪器都是 CT 机;
(2)假设两张 CT 图像来自于同一患者;
(3)假设该患者两次扫描用的是同一台 CT 机;
(4)假设该 CT 机性能良好,可以正常工作;
(5)假设两次 CT 扫描成像的断层位置相同;
(6)假设该患者腹腔内占位性病变的影响占主要因素;
(7)假设两次扫描时间间隔足够长,使得占位性病变对 CT 图像的影响足够显著。
模型的建立与求解
在 Demons 形变模型中,有关像素点偏移量的计算思想实际上是来自于光流场模型[3],在配准的过程中,由于这两幅医学图像都是通过同一台 CT 机扫描得到的,因此二者灰度分布相等,即所对应的灰度点可以近似为在单位时间内参考图像里面的像素点的平面运动。Demons 模型的优点就是可以在一定程度上有效、精准地校正图像的非刚性变形,因此该模型就是基于 Kullback-Leibler 距离对 Demons 多模态图像配准算法[4]进行改进。
基于 Kullback-Leibler 距离的互信息梯度的计算模型
在多模态配准的过程中,互信息常被用于医学图像配准的测度,而互信息越大则可以说明图像匹配效果越好。当将两图像进行配准时,所包含的信息总量达到最大,与此同时,互信息也会达到最大。而互信息梯度(MIG)代表的含义是互信息变大的方向,如果我们在 Demons 形变力的基础上加上 MIG 这个分量,那么不难发现 Demons 形变图像中的各个像素点都会向着互信息变大的方向发生形变,从而实现多模态医学图像的配准。
在概率论和信息论中,互信息是两个随机变量间相互依赖性的量度,即统计相关性,一般可以用熵来表示,其物理含义是代表系统的复杂性和不确定性。假设该患者的两幅CT 图像分别是 M 和 N,两图像中任意像素点的灰度值分别用
2021年认证杯SPSSPRO杯数学建模A题(第一阶段)医学图像的配准全过程文档及程序相关推荐
- 2023年认证杯SPSSPRO杯数学建模B题(第一阶段)考订文本全过程文档及程序
2023年认证杯SPSSPRO杯数学建模 B题 考订文本 原题再现: 古代文本在传抄过程中,往往会出现种种错误,以至于一部书可能流传下来多种版本.在文献学中,错误往往被总结成"讹&quo ...
- 2020年数维杯国际大学生数学建模A题中国螺纹钢需求预测求解全过程文档及程序
2020年数维杯国际大学生数学建模 A题 中国螺纹钢需求预测 原题再现: 螺纹钢是我国产量最大的钢铁产品之一.钢筋广泛应用于房屋.桥梁.道路等土木工程建设中.它是基础设施建设不可缺少的结构材料.合 ...
- 2019年国赛高教杯数学建模C题机场的出租车问题解题全过程文档及程序
2019年国赛高教杯数学建模 C题 机场的出租车问题 原题再现 大多数乘客下飞机后要去市区(或周边)的目的地,出租车是主要的交通工具之一.国内多数机场都是将送客(出发)与接客(到达)通道分开的.送 ...
- 2018年国赛高教杯数学建模A题高温作业专用服装设计解题全过程文档及程序
2018年国赛高教杯数学建模 A题 高温作业专用服装设计 原题再现 在高温环境下工作时,人们需要穿着专用服装以避免灼伤.专用服装通常由三层织物材料构成,记为I.II.III层,其中I层与外界环境接 ...
- 2022年认证杯SPSSPRO杯数学建模C题(第一阶段)污水流行病学原理在新冠疫情防控方面的作用求解全过程文档及程序
2022年认证杯SPSSPRO杯数学建模 C题 污水流行病学原理在新冠疫情防控方面的作用 原题再现: 2019 年新型冠状病毒肺炎疫情暴发至今已过两年,新型冠状病毒历经多次变异,目前已有 11 种 ...
- 2022年亚太杯APMCM数学建模大赛C题全球是否变暖求解全过程文档及程序
2022年亚太杯APMCM数学建模大赛 C题 全球是否变暖 为方便各位阅览及了解掌握亚太杯的写作技巧,这里非技术使用中文,公式部分由于翻译过程繁琐使用英文来撰写此文章. 原题再现: 加拿大的49. ...
- 2019年数维杯国际大学生数学建模C题猪肉价格波动分析求解全过程文档及程序
2019年数维杯国际大学生数学建模 C题 猪肉价格波动分析 问题重述: 中国是一个育种大国.猪肉产业在畜牧业中起着主导作用.同时,猪肉已成为人们餐桌上的主要肉类来源之一.随着人口的增加,消费和生活 ...
- 2022年华数杯数学建模B题水下机器人的组装计划解题全过程文档及程序
2022年华数杯全国大学生数学建模 B题 水下机器人的组装计划 原题再现: 自来水管道清理机器人(Water pipe cleaning robot,简称WPCR)是一种可在水下移动.具有视觉和感 ...
- 2022年华数杯数学建模A题环形振荡器的优化设计解题全过程文档及程序
2022年华数杯全国大学生数学建模 A题 环形振荡器的优化设计 原题再现: 芯片是指内含集成电路的硅片,在我们日常生活中的手机.电脑.电视.家用电器等领域都会使用到,是高端制造业的核心基石.芯片的 ...
最新文章
- C#中String对象转换为Font对象的方法
- delphi xe6 android ListView增加 Header或Footer 的方法
- 如何改造UE4用于赛璐璐3D卡渲?这里有一份日本大厂的实操分享
- 一步步使用SAP云平台的WebIDE开发SAP UI5应用
- 20110128 学习记录:Datalist中的item项不自动换行而只用空格隔开
- android 音频控件,android AudioManager音频控件
- AndroidStudio安卓原生开发_UI高级_RecyclerView_循环复用控件_代替ListView_GridView---Android原生开发工作笔记121
- php中的代码延迟函数sleep() usleep()
- 不洗袜子的高文博_那个孩子在夏天中旬用高袜子大笑?
- axure后台示例_AxureUX中后台管理信息系统通用原型方案 v1.1
- 电子邮件--详解SMTP和POP3协议
- 最适合养老的20座城市
- 传奇各大地图编号代码库
- 耗子大叔分享的软件设计原则
- c语言average的用法,平均函数average的一般和不一般的用法
- linux 文本方式怎么安装,Linux文本方式安装系统
- 自动驾驶专题介绍 ———— 动力传动系统
- GB\T和GB国标文件查询、下载方法(建议收藏)
- HTTP 模拟post请求
- ubuntu20有线网络无线网络
热门文章
- springboot集成elasticsearch7实现全文检索及分页
- QQGC?揭秘QQ的AI绘画大模型技术
- 消费者行为分析—购买过程
- c语言中push函数pop函数用法,perl中pop与push函数
- 2018年第九届蓝桥杯省赛试题及详解(Java本科B组)
- 走进Linux(一切皆文件)
- SQL SERVER 2014 清理数据库日志文件踩的坑!
- #Loj121 动态图连通性
- 双树复小波变换Dual-Tree Complex Wavelet Transform在信号处理中的应用
- 列车调度系统属于计算机应用中的,2020年12月网络教育​统考《计算机应用基础》基础知识试题及答案1...