第 3 7 卷 第 2 期
2014年 6 月
辽宁师范大学学报(自然科学版)
Journal of Liaoning Normal University (Natural Science Edition)
Vol. 37 No. 2
Jun. 2014
文章编号:1000-1735(2014)02-0198-06
doi:10. 11679/lsxblk2014020198
基 于
Marching Cubes
方法计算分子体积
杨忠志,马丽
(辽宁师 范大 学 化 学化 工 学 院 ,辽 宁 大 连 116029)
摘 要 :利用M a r c h i n g C u b e s 方法,基于分子形貌理论,提出了分子等值面模型体积的
计算公式. 首先 ,利用二进制和十进制的转换关系定义体元顶点的标号,并将与分子表面
相交的体元归纳为1 6 种不同的构型,其中包括基本构型和特殊构型.然后针对不同构型
给出相应的体积计算公式,最终通过求和得到整个分子的体积.该方法的优点在于不需
要重建分子表面,即可通过已知网格数据和阈值直接得到分子的体积,易于理解,并具有
较高的计算效率. 最后 ,通过 几个典型分子体积计算结果的对比,证实了该方法准确有
效. 此外,该方法同样适用于其他从规则网格数据中提取的等值面模型体积的计算.
关键词:M a r c h i n g C u b e s ; 分子形貌理论;分子体积;等值面
中图分类号:0 6 4 1 . 121 文献标志码:A
分子形貌理论[ w ] 利用电子运动的经典转折点来定义分子的表面.当电子运动到三维空间内原子
核附近的某一处时,其能量恰好等于其势能,即在该点处电子的' 均动能为零.假设此时该电子的能量
等于分子第一电离能的负值,则 分 子 表面 上 任 意 一点 R 满 足 G ( — I E ) = { R : F ( R ) = — I E },其 中 ,
F ( R ) 为 分 子中 1 个电子所受到的作用势(P A E M ) ,I E 代表分子的第一电离能. 在分子形 貌理论中,分
子的表面是由电子所受到的作用势等于第一电离能负值的点的集合构成的等值面.由于很难在空间上
取得连续的分子形貌测量数据,只能在三维空间上取一定的间隔,计算每个网格点上的P A E M 值和电
子密度值,使测量数据均勻分布在三维网格点上.体积是三维空间网格模型的重要属性之一,计算分子
的体积等价于计算由分子表面形成的等值面内部空间的体积.
许多国内学者在体积计算方面提出了各自的计算方法. 例 如 ,戴等人[3]提出了任意三维实体网格
模型的体积特征提取算法.孙等人M 应 用 M a t l a b 建立了三维医学模型的体积计算方法. 胡等人[7]给出
了基于散乱点云的体积计算方法. 王等人[8]针对任意三角网格模型提出了体积的计算方法.本文拟发
展一种从三维规则网格数据中,根据阈值计算分子等值面模型体积的方法.
M a r c h i n g C u b e s [M1] ( M C ) 方法是当前比较流行的重建等值面的方法,被广泛地应用于很多领域.
本文利 用 M C 方 法,基于分子形貌理论,总结出分子体积的计算公式.该方法通过分析等值面和规则网
格相交情况,逐个计算与等值面相交的体元的公共部分的体积,然后与处于等值面内部体元的体积加
和得到分子的体积. 最后通过与文献计算结果的比较,证实了该方法的准确性.
收稿日期= 2014-03-10
基金项目:国家自然科学基金项目(21073080 ; 21133005)
作者筒介:杨忠 志(1940-),男 ,吉林舒兰人,辽宁师范大学教授,博 士 ,博士生导师.