信鸽体表的分形维数计算方法及程序实现.docVIP

信鸽体表的分形维数计算方法及程序实现 　　摘要:基于激光扫描系统获取的信鸽体表三维信息以及体视显微镜获得的图像,采用盒子维方法,编程计算了信鸽胸部、背部和翅膀三个主要部位的分形维数。 　　关键词:分形几何;分形维数;三维信息;图像 　　中图分类号:TP311文献标识码:A文章编号:1009-3044(2010)18-4938-02 　　Fractal Dimension Caculation Method of Pigeon Surface and Program Realization 　　ZHANG Chun-hua, ZHANG Wei, JIANG Xiao-ming 　　(Aviation University of Air Force, Changchun 130022, China) 　　Abstract: The fractal dimension of the surface of wing, breast and back was calculated using box dimension method based on the three-dimensional information scaned by 3D LASER SCANNER SYSTEM and the images obtained by stereomicroscope using programing method. 　　Key words: Fractal geometry;Fractal dimension;Three dimensions information;Image 　　分形几何是刻画空间不规则形体的有效工具[1],它能够在更深层次、更广范围内描述表面形貌的几何特征。分形维数是分形几何对象的一个非常重要的特征参数,是粗糙表面不规则的无标度维数。采用基于体表的三维信息和体表图像计算分形维数并编程实现,从而量化描述信鸽体表的分形特性,为信鸽体表的仿生减阻研究提供基础信息。 　　1 基于三维扫描信息的分形维数计算 　　计算分形维数的方法有很多种,盒子维运算方便快速,有一定的精确度,而且易于程序实现[2],其计算公式为: 　　(1) 　　式1中,ε为边长,N(ε)为盒子数量。采用LSV50型三维激光扫描系统获取信鸽体表的三维信息,并计算出信鸽体表的线维数和面维数。信鸽胸部、背部和翅膀试样扫描的三维数据中,均匀抽取其中50条扫描线进行分析,求取它们的线维数,以反映不同部位的分形特征。信鸽胸部横向剖面曲线的线维数计算如图1所示,将剖面扫描点所在平面空间划分成边长为ε的盒子,计算出对应于划分尺度ε的非空正方形数N(ε),在双对数坐标下将6对ε值和N(ε)值进行直线回归,得到的斜率绝对值1.0131即为线维数值。信鸽表面的计盒维数(简称面维数)的计算方法为,将信鸽三维表面所嵌入的立体空间划分成边长为ε的立方体盒子,数出曲面所占据的盒子数N(ε),如果是不连通的点集则数出点所占据的盒子数N(ε),在双对数坐标下将测量尺度ε和非空正方体的个数N(ε)进行直线回归,得到的斜率绝对值加1即为面维数。线维数和面维数的计算中,直线的起点和终点,分别对应着分形尺度下限和上限,起点、终点之间的尺度范围为分形尺度。直线的拟合误差也可以用回归显著性的t检验方法来对拟合结果进行检验。 　　2 基于图像的分形维数计算 　　图像技术在表面形貌研究中已经有广泛应用,如果一个物体的表面是分形的,则由它产生的图像灰度表面也具有分形的性质,且两个分形的特征是一致的,因此可以由灰度图像提取信鸽表面形貌的分形维数[3-4]。将XTJ-30体视显微镜获取的形貌图像(选择同样颜色的信鸽进行研究)输入计算机中,采用Photoshop软件对图像进行细致处理,去除噪音信号后,对数字图像进行数值处理,使图像具有合适的对比度,在最佳的对比度状态下进行分形特征的提取。利用灰度图像提取表面分形维数,本研究主要采用基于二维情况的计盒维数求法,图2给出了根据信鸽体表图像计算分形维数的程序框图。 　　3 程序实现 　　基于三维几何信息和表面图像,采用同样的盒子维的计算方法,求取信鸽体表不同部位的分形维数的过程通过编程实现。利用可视化编程语言VC6.0设计了基于Windows平台的分形维数计算程序,主要包括数据的批处理、二维计算、三维计算以及图像计算几大模块,批处理模块主要完成从扫描数据中按照要求抽取一定数量的横向和纵向线扫描信息,二维计算模块和三维计算模块,主要完成线维数和面维数的计算,程序界面如图3所示。 　　4 结束语 　　采用盒子维方法,基于三维信息和图像,运用VC6.0编程计算了信鸽体表主要部位的分形维数。基于三维信息计算的面维数的范围为2.8988