三维医学模型表面积和体积计算方法研究.PDFVIP

第 39 卷第 15 期 数学的实践与认识 V o l. 39　No . 15　 2009 年 8 月 M AT HEM A TICS IN PRACTICE AND T HEORY A ugust, 2009　 三维医学模型表面积和体积计算方法研究 孙宇锋1 , 　房少梅2, 　曾美萍3 , 　池雄标1 ( 1. 韶关学院 数学与信息科学学院, 广东 韶关　512005) ( 2. 华南农业大学 理学院, 广东 广州　510631) ( 3. 韶关学院 医院, 广东 韶关　512005) 摘要: 　应用 M ATLAB 语言描述了测量任意 给定的颌、颧骨区域的表面积和体积的方法, 并 且在临床上用 实例给出了计算结果. 关键词: 　三维模型; 表面积; 体积; CT 图片 1　引　　言 　　近年来, 由于 CT 专利保护技术的实施, 扫描的二维图片很难在普通的计算机上进行三 维影像重现, 对特定的病???体也不易实现精确测量. 因此, 如何将 CT 机上的二维图片移植 到普通计算机上( w indow s 操作系统) 进行病情分析、三维重建及测量是大家普遍关心和研 究的课题. 在口腔颌面部矫形、美容临床应用中, 往往需要提供准确多参数定量测量的有效手术方 法. 这种方法将多模式二维 CT 图片进行重建来测量数据, 通过颌骨、颧骨的任意点方便观 察三维运动轨迹, 比较、设计患侧与健侧颌面骨模型与空间位置的解剖学和几何学差异, 实 现对颌骨、颧骨手术精确复位和美容美学意义上的准确空间三维复位, 恢复与再塑颌面中部 固有的细腻与谐美形态( 包括空间区域、宽度、高度、厚度和凸凹度等) , 减少术后面容畸形并 强化仪表曲美. 本文主要探讨基于 MAT LAB 的对上颌骨、颧骨三维重建体的体积和表面积测量和计 算技术, 为颌面部矫形、美容提供准确多参数定量测量的有效方法. 2　三维重建模型表面积和体积算法 2. 1　问题的模型及处理思路 图像处理可以看作为一个输入输出系统模型: Q 0 → $ → Q 1 其中 Q 0 表示观测到的单个图像或图像序列, $ 是图像处理器, Q 1 把经过处理的目标图像的 [ 1] 结果输出 . 因此, 输入数据 Q 0 的建模影响着处理器 $ 的算法. 前人的实验数据和理论研究结果表明, 骨骼是各向异性的粘弹性体. 颅骨作为人体骨 骼的重要组成部分, 颅骨内板、外板拉伸时的弹性模量、破坏应力都会影响其力学性能[ 2 ] . 因 此, 研究颅骨的近似光滑边界( 曲线或曲面) 是一个关键问题. 收稿日期: 2009-01-19 基金项目: 国家自然科学基金( ; 广东省自然科学基金( 5300889) ; 广东省科技计划项目( 20051140019) 80 数　学 　的　实　践　与　认　识 39 卷 用经典的Euler -Lag rang e 变分方法就是在域 D R2 中寻找函数 u( x , y ) , 满足 D 边界 5u 5u 5D = L 的边界条件, 并使泛函: I [ u( x , y ) ] = F x , y , u, , dx dy 达到极值, 式中函