- 1、本文档共15页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
A题 血管的三维重建
四院九九 向为 王瑛 伍微
摘要:
对此问题我们提出的是一个搜索模型,从而求出原始球(形成包络的球)的半径,以及中轴线。
首先我们把每一个切片图的bmp文件转化为一个512×512的0,1矩阵(0代表无象素,1代表有)。再求切片的轮廓线,也存为0,1矩阵。
然后在题设条件下证明了以下几个结论:1.能够被切片包含的半径最大的圆的半径等于原始球(形成包络的球)的半径;2.可以被切片包含的圆的半径一定小于等于原始球的半径;3.不能被包含于切片的圆的半径一定大于原始球的半径。
根据上述结论,对每一个切片,只要找到能够被切片包含的半径最大的圆就能求出原始球的半径。我们设计了一个二分搜索算法求该半径,再根据半径搜索求出每个切片和中轴线相交的点的坐标。
算法概要:1.求粗略的rmax(原始球半径上限);2.求粗略的rmin(原始球半径下限);3.切片内所含的标准圆(以原始球的半径为半径的圆)的圆心O不可能在离切片边缘距离小于rmin的点。则在切片上去掉这些点后,就得到了O可能存在的位置区域;4.用二分法不断缩小原始球半径r0可能的取值范围,从而求出r0;5.求出切片内所含的标准圆的圆心的坐标,这就是为切片和中轴线相交的点的坐标。
搜索原始球的半径时,因为所给的图象数据的精度有限,所以在执行算法时遵循以下规则:1.在判断一个圆面是否有一部分在切片的外部时,把这个圆和切片都看成由一个一个的像素点组成,若圆面中有像素点在切片的外部,就认为圆超出了切片;2.几何上的半径为r的圆面转化为像素圆面的方法是,当一个像素的中心离圆心的距离小于等于r时,这个像素属于这个圆面;3.因为所给数据精度有限,所以包含于切片中的以原始球的半径为半径的圆可能不止一个。取这些圆的圆心的重心为含于切片内的标准圆的圆心,即中轴线与切片的交点。
最后对中轴线上已知的点进行拟合,由于一个z只对应于一个x和一个y,故可分别对其投影分别在YZ、ZX平面上进行多项式拟合,求出y=f1(z)和y= x=f2(z)。则中轴线的空间方程即为上两式的联立。
用本模型可以对每一个切片求出一个r0(原始球的半径),共求出100个r0。这些r0的平均值为30.1706,方差为0.0176。可见r0的精度很高。
这个模型采用了较巧妙的搜索算法,求解结果可以达到很高的精度。能够很好的重建包络,并加以分析,在医学方面有很大作用。
一.问题重述
断面可用于了解生物组织、器官等的形态。例如,将样本染色后切成厚约1? m的切片,在显微镜下观察该横断面的组织形态结构。如果用切片机连续不断地将样本切成数十、成百的平行切片, 可依次逐片观察。根据拍照并采样得到的平行切片数字图象,运用计算机可重建组织、器官等准确的三维形态。
假设某些血管可视为一类特殊的管道,该管道的表面是由球心沿着某一曲线(称为中轴线)的球滚动包络而成。例如圆柱就是这样一种管道,其中轴线为直线,由半径固定的球滚动包络形成。
现有某管道的相继100张平行切片图象,记录了管道与切片的交。图象文件名依次为0.bmp、1.bmp、…、 99.bmp,格式均为BMP,宽、高均为512个象素(pixel)。为简化起见,假设:管道中轴线与每张切片有且只有一个交点;球半径固定;切片间距以及图象象素的尺寸均为1。
取坐标系的Z轴垂直于切片,第1张切片为平面Z=0,第100张切片为平面Z=99。Z=z切片图象中象素的坐标依它们在文件中出现的前后次序为
(-256,-256,z),(-256,-255,z),…(-256,255,z),
(-255,-256,z),(-255,-255,z),…(-255,255,z),
……
( 255,-256,z),( 255,-255,z),…(255,255,z)。
试计算管道的中轴线与半径,给出具体的算法,并绘制中轴线在XY、YZ、ZX平面的投影图。
二.基本假设
管道中轴线与每张切片有且只有一个交点;
球半径固定;
切片间距以及图象象素的尺寸均为1;
血管无断裂无突变,即管道表面光滑且连续;
对切片拍照的过程中不存在误差,数据误差仅与切片数字图象的分辨率有关;
假设中轴线的曲率相对于标准圆的曲率较小。
三.术语及符号说明
术语说明:
切片: 指包络表面与切平面相交而成的曲线围成的区域。
原始球: 指沿中轴线滚动形成包络的球。
标准圆: 指与原始球等半径的圆。
符号说明:
r0 原始球的半径。(由基本假设中2.可知,半径r0固定,为常数)
rmax r0的上限
rmin r0的下限
Oj 在切片内部且与切片相切的标准圆的圆心
C Oj可能存在的区域
s
您可能关注的文档
最近下载
- 网络空间信息安全(第2版)习题与思考题参考答案.pdf
- 微肿瘤PTC体外药敏检测服务2023年.pdf VIP
- 教科版四年级科学上册第三单元《第1课让小车运动起来》优秀教学设计.docx VIP
- 病例报告表完整版.doc
- 电影售票管理平台.docx
- 统编版语文三年级上册《对话中提示语的位置及标点》课件.pptx VIP
- AQ_ 4272-2016 铝镁制品机械加工粉尘防爆安全技术规范.pdf
- 4第四章-金融市场及其构成(货币金融学(蒋先玲编著)第3版ppt课件可编辑).pptx
- 部编版语文二年级下册 第三单元大单元任务群“发现汉字的秘密”(一等奖创新教学设计).pdf VIP
- 《活着》读书分享.pptx VIP
文档评论(0)