三维骨架提取Main.doc

使用广义势场法计算三维骨架 主要结构体 typedef struct //点类型（立体的点坐标） { short x; short y; short z; } VoxelPosition; typedef struct //向量类型 { double xd; double yd; double zd; } Vector; enum CriticalPointType //关键点的类型 { CPT_SADDLE = 1, // 1 鞍点 CPT_ATTRACTING_NODE, //2 引力点 CPT_REPELLING_NODE, //3 斥力点 CPT_UNKNOWN //4 其他 }; struct VoxelPositionDouble //点 { double x; double y; double z; }; struct CriticalPoint { VoxelPositionDouble position; //关键点的坐标 CriticalPointType type; //关键点的类型 Vector evect[3]; //3个特征向量数组 double eval[3]; //3个特征值 }; struct Skeleton { VoxelPositionDouble *Points; //存储骨架上的点 int sizePoints; //初始化时分配的最大储存空间本程序是50000 int numPoints; //骨架上的点的个数 int ** Segments; //骨架段为二维数组 //numSegments行 4列的数组 LEFT RIGHT FIRST LAST int sizeSegments; //初始化分配的最大存储空间，本程序是5000 int numSegments; //骨架段的个数 }; 骨架段的说明： LEFT RIGHT FIRST LAST 第一段 第一个骨架点index 最后一个骨架点的index 第二段 第三段 …… 立体长宽高： L = 123 M = 116 N = 114; 计算势场用到的距离的次方数 fieldStrenght = 8; 计算散度值时和阈值有关的 percHDPts = 3; vfin = false; vfout = false; interactive = false; //这些变量代表的实际意义不清楚 #define SURF 100 // surface voxel #define BOUNDARY 110 // boundary voxel - participates in potential field //calculation #define INTERIOR 200 // interior voxel #define PADDING_MIN 210 // added voxels in order to thick the object #define NR_PAD_VALUES 40 // are in this range: PADDING_MIN to PADDING_MIN + NR_PAD_VALUES #define EXTERIOR 0 // background (exterior to the object) voxel (air) 貌似 INTERIOR点是立体的体素点，EXTERIOR点是空气， SURF也是立体的点，是暴露在空气中的点？？？ 具体算法： 定义了Skeleton变量Skel， 函数AllocateSkeleton对它进行初始化 为Skel分配50000个点 5000个段 VoxelPositionDouble *Points; //初始分配50000个空间 int sizePoints; //点最多50000 int numPoints; //0 Skel: int ** S