专利名称：基于图像的三维物体不连续空洞体积的测定方法
技术领域：
本发明涉及一种三维立体图像中测量物体体积的方法，更具体的涉及了 一种通过MicroCT系统获得的三维立体图像来测量三维物体不连续空洞体 积的方法。
背景技术：
MicroCT是一项新型技术，MicroCT首先在工业无损检测中有不可替 代的作用，其次凭借其系统中的完美数学模型的基础，可以预知的作用就是 将在在极其微小的微观世界发挥举足轻重的地位，甚至可以探测纳米级的微 观形态。这些取决于前端采集设备的精确性，原始数据处理成断层数据的数 学模型的准确性以及高效性。这两者对于系统的性能起着非常重要的作用， 但是断层图像的后处理的操作平台也是非常重要。有了方便的操作平台，可 以实现很多实际需要的检测数据，比如三维图像的重建，可以在无损的情况 下探测物件的侦测不同平面的任意提取，可以非常方便的观察各个层面的 情况。
由于工业领域的特殊需要，在完善三维图像的可视化的效率与效果后出 现了一些迫切需要解决的问题。在现在很多的材料的加工过程中需要检验材 料的内部属性，现有技术中研究较多的是多孔材料的内部探测。当然探测直 接显示出所需探测物体的三维效果图，但是理论和实际的测量无法达到。通 过数据精确计算出三维实体中的铝泡的体积以及铝泡的位置这样的数据还 是比较复杂。现有的方法主要有通过层与层之间的变化来判断各种大小的孔 的出现频率，所以精确的体积计算以及位置的获取都将很难实行而且相对的 误差较大。本发明由此而来。
发明内容
本发明提供了一种基于图像的三维物体不连续空洞体积的测定方法，解 决了现有技术中三维实体中的多孔材料的体积的测量比较困难、相对误差较 大等问题。为了解决现有技术中的这些问题，本发明人提供的解决方案如下 一种基于图像的三维物体不连续空洞体积的测定方法，包括以下步骤 将原始图像载入图像处理装置中并进行预处理原始图像； 根据预处理后的图像根据原始图像的灰度级确定边界；撷取图像目的区
域数据载入内存，在内存中构造空间八叉树通过遍历八叉树节点累积确定不
连续空洞体积。
优选的，所述原始图像是通过MicroCT系统三维扫描后得到的三维可 视图像。
优选的，所述的原始图像预处理是将原始图像经灰度线性变换处理。 优选的，所述的灰度线性变换处理方式是通过灰度分布函数映射得到， 所述的灰度分布函数为
隱c'
d國
(/(x,力一a)+a' 他力-c)+c'
其中，a， c， d， b代表原始图像灰度级的四个分界点，a'， c'， d'， b'
分别代表映射后的图像灰度级的分界点；f'(x，y)为变换后的图像的灰度分布 函数；f(x，y)为原始图像的灰度分布函数且a《f(x，y)<c， c《f(x，y)<d ， d《 f(x，y)《b。
优选的，所述的图像目的区域数据是先撷取原始图像矩形块，通过原始 图像矩形块的坐标与原始图像的层数获得。
优选的，所述的空间八叉树选自规则八叉树、线性八叉树或一对八式八 叉树。
优选的，所述的图像目的区域数据是通过空间八叉树进行管理分配，所 述的空间八叉树基本结构包括数据地址指针，指向父单元的指针，指向子单 元的指针数组，所述指向子单元的指针数组包括表示是否被遍历的标志位、 表示是否满足空洞条件的标志位，载入图像的宽度，载入图像的髙度，载入 图像的层数；按照上述的结构分配设置空间八叉树各单元。
5优选的，所述的空间八叉树遍历包括如下步骤-
从八叉树第O层开始遍历，第O层第一个小单元开始向后查找；判断该单元是否遍历过和该单元是否完全是空洞单元；
当该单元未曾遍历且完全是空洞单元时向父单元上溯；判断父单元是否是空洞单元和子单元是否可分；
当父单元不是空洞单元且子单元可分时，退级到子单元，查看子单元的
八个次级子单元；査看每个子次级单元的周边的三个同层子单元；判断同层子单元的遍历情况和是否满足空洞的条件；
当同层子单元非空洞且该单元可分时将与同层子单元相邻的次级子单元八分，査看四个八分单元的相邻同层子单元；判断该单元是否遍历过和该单元是否完全是空洞单元
当相邻同层子单元非遍历且是空洞单元是查看子单元是否可分，如此循环；当循环到最后一级不可分单元的周围六个同级单元全部不属于空洞单元或者已遍历过时结束遍历。
优选的，所述的空间八叉树遍历时，当子单元的父单元是空洞单元时继续上溯。
优选的，所述的空间八叉树遍历时，当单元不满足空洞条件且单元不可分时，停止这个单元向周围扩展，查看周围六个同级单元。
本发明的技术方案中设计的空间八叉树结构，是一种体数据结构，该结构的好处在于可以很好的连接上层和下层的关系。通过首先设置一个体数据
的最低层，为方便计算每个边的点数设置成N-" + l， m是体数据每条边
可以分的最小像素数减l。最小的边所组成的基本单元为第0层，显然，第
0层共有2^个单元。每个单元被赋予一个标识位。然后以8个相邻的第0
层单元组成第l层，显然，第1层共有""8个单元，依此类推构建这样的一个八叉树的数据结构。这样利用八叉树，通过从底层的第0层不断的通过判断标志位，方便的填充到更髙的层次，从而忽略掉对体绘制无用的废点，减少光线通过体数据时对采样点的计算，进而有效地节约时间，加快体绘制的速度从而使体绘制的实时互动成为现实。
本发明由于多孔材料内部大多数的孔之间是不连通的，这就使得将平面区域增长的方法可以推广到三维空间，这样通过八叉树层与层之间的依存关系找出非常近似的物体空间边界，测量得到物体的孔体积。
采用本发明技术方案的有益效果在通过构造八叉树应用到MicroCT系统中三维可视化后处理中的有孔实体的孔体积的测量中，不但可以很精确的定位空洞的位置，而且可以很方便的实现对MicroCT系统中三维可视化后处理中的有孔实体的空洞体积的精确的测量。该方法简单，有效。本发明是针对三维体数据中闭合的不规则孔隙的大小、位置进行有效计算的一种方法，可以应用在三维实体空洞体积的精准测量上。
下面结合附图及实施例对本发明作进一步描述
图1是本发明实施例基于图像的三维物体不连续空洞体积的测定方法的流程图2是本发明应用例中采用的原始图像，未进行预处理；图3是本发明应用例中预处理后的图像，进行线性变换预处理；图4是本发明应用例中原始图像进行预处理的基本方法示意图；图5是本发明应用例中构造八叉树的空间遍布逻辑图。
具体实施例方式
为了更详尽的表述上述发明的技术方案，以下本发明人列举出具体的实施例来明技术效果；需要强调的是，这些实施例是用于说明本发明而不限于限制本发明的范围。
实施例如图l的流程图，该方法步骤如下
1.步骤l首先载入原始图像，进行线性变换预处理方法如下，分段灰度变化遵循以下规则
c'--"
c--a
-c'
-c
-《其中，a, c, d， b代表原始图像灰度级的四个分界点，a'， c'， d', b'
分别代表映射后的图像灰度级的分界点；f'(x，y)为变换后的图像的灰度分布函数；f(x，y)为原始图像的灰度分布函数且a《f(X，y)<c, c《f(x，y)<d ， d《f(x，y)《b。
2. 步骤2, 3选取图像中的目的区域， 一般为需要观测的有用数据载入新的内存中，随后所有的操作均对这块数据区进行。在第一张断层图像中用鼠标选取一个需要的矩形块，以象素为单位，通过矩形块的四个坐标(xl，yl)(x2， yl)( xl，y3)( x2 ， y3)分别代表矩形的左上，右上，左下，右下角的坐标，设原始图像的宽为width,髙为height,断层图像的层数为pieces,
原数据区指针为pData,令i = 0 -pieces—1; j = yi~ y3; k = xi x2那么新载
入的数据即为pData[i*height*width+j*width+k:
3. 在步骤4中对新的内存区中的体数据按照八叉树的结构的进行管理每个八叉树单元的基本结构包含数据地址指针pBit，指向父单元的指针pParent，指向子单元的指针数组pChildren[8] ， pChildren[8]有八个元素，两个标志位分别是表示是否被遍历的标志M、表示是否满足空洞条件的标志S,载入图像的宽度newwidth =x2 —xl,髙度newheight =y3 —yl，层数为pieces, 取newwidth、 newheight、 pieces中的最大值设为max; 令max
《2 + l,解得m的最小值，把体数据设为每边有2" + l个体素的立方体，
那么每条边可以分为2"段，整个体数据可以分为8"个基本现单元。按照八叉树的理论首先对最底层第0层leve10的小单元进行操作，即相邻的8个体素形成一个最底层小单元，每个单元拥有数据地址pBit，指向单元的左上角的体素的地址；因为最低层单元没有子单元，所以指向子单元的指针数组pChildren[8]的八个元素指针为空，其他层级单元的pChildren[8]的八个元素指针按照从左上到右下的指向8个子单元；若levelO层的小单元的8个体素满足是空洞的域值要求或是非leve10级的所有子单元满足空洞条件，设置标志位S=true，否则设置S-false; M=false,表示这个单元未遍历，M=true表示遍历了； pParent表示指向父单元的父指针。分配最底层，也即从levelO层的小单元开始设置，将levelO层小单元的 指向子单元的指针数组pChildren[8的八个元素指针设置为空，标志位 M=falSe，标志位S的设置通过判断8个顶点的灰度值是否全为0来设置， 若是，则S4me，若否，则置S-false;分配完最底层，接着进行上一层(leve11) 的设置，以八个相邻的小单元为子单元合成一父单元，即新的正方体，设置 这个父单元的pChildren8的值分别为从左上到右下指向八个子单元的指 针，标志位M-false，标志位S的设置通过判断八个子单元是否都满足空洞 条件，即八个子单元是否都满足S = true，若是，则设置标志位S = trUe;若 否，则设S-false;把八个子单元的父指针pPareiit指向这个父单元；进行 level2的设置，以八个相邻的父单元为子单元再组成一个父单元，参数设置 方法同步骤4.2中的设置，依此类推，直至设置到levelm;
4. 步骤5—10，从八叉树的最底级(leve10)开始对八叉树进行遍历， 从第一个单元开始向后査找如果一个单元是遍历过的单元或不是空洞单元 则跳过此单元向下继续寻找，如果找到的单元通过他的标志位S、 M判断满 足是未遍历的空洞单元，则向此单元的父单元遍历，继而查看父单元是否是 未遍历的满足空洞条件的单元，如果是则继续向上上溯，直到不满足"未遍 历"或是"空洞单元"的两个条件之一停止上溯。具体步骤为从八叉树的 levelO层开始对八叉树进行遍历；从leve10的第一个小单元开始向后查找； 判断此单元是已经否遍历过，若是，即M = tnie，则回到步骤6;若否，M =false，则进行步骤8;判断此单元是否完全是空洞单元，若否，即S = falSe, 则回到步骤6;若是，则此单元同时满足M二false并且S-true，首先定义 体积变量voliime-0;进行步骤9;向此单元的父单元上溯；查看父单元是 否是空洞单元；若满足条件M = falSe，并且S-true，则回到步骤9继续向 上上溯；若否，即S = false时，停止上溯，设停止单元为leveln级的单元T， 其中1《n《rn (m为最大级数)，退回到level(n-l)级找到上溯到它的那个 子单元Tu , 0《u《7,进行步骤ll;
5. 步骤11 23中，查看停止单元，退回上溯到它的子单元。査看其是 否是可分的(即判断他是否是最底层的单元)，
(*1)如果它是可分的继续退到这个子单元的次级子单元0, 1 ，2 ，3， 4 ， 5 ， 6， 7分别代表次级子单元在子单元中的位置。其中以次级子单元6为例介绍下面的空间增长过程(图5所示)
(1) 次级子单元6有三个相邻同级外部单元，步骤14首先判断这些外部 同级单元是否遍历过，遍历过就停止。
(2) 否则判断其是否满足是空洞单元，如果满足重复步骤9向下(即向 上进行上溯查看父单元是否满足空洞条件以及是否遍历)。如果不满足，假 设6c不满足，査看与6c相邻的同级单元6，观察其是否可分(步骤16 — 20)
<1>如果可分就把它八分，査看与6c相邻的四个八分子单元(如图3 中单元6内的虚线所示，划分出四个子单元)。分别査看这四个单元的相邻 同层单元是否遍历过，如果遍历过则终止这个单元的计算。如果没有，则查 看其是否是空洞单元，如果不是空洞单元则要重复步骤16 — 20，再次细分这 个子单元。而如果其是空洞单元，重复步骤11向下假设6a满足条件，判 断其是否可分，可分则把6a八分，查看每个八分子单元的相邻3个同级单
元o
<2>如果不可分，则说明这个单元是最低单元且不满足空洞的条件，所 以停止这个单元的计算。
(*2)如果不是可分的(步骤21所示)，则说明这个单元是最低单元，如 果此单元的六个相邻同级单元全都或是被遍历过或是都不满足空洞条件则 停止，否则传递到步骤9，以最小的单元向上遍历重复步骤9以下的操作。
6、至此第一个孔通过遍历可以得到，孔的体积误差是边缘的一个最小 单元的边界。通过对体数据的所有最底层的单元进行遍历可以求出感兴趣区 域的所有的泡的体积。 应用例
一种MicroCT系统中三维可视化后处理中的有孔实体的孔体积测量方 法，步骤为c--a
-c'
-c
6'-
(/(xj)-c)+c'
其中，以上三式分别当a《f(x， y)<c， c《f(x，y)<d , d《f(x，y)《b满 足，其中a， c， d, b代表原始图像灰度级的四个分界点，a ， c ， d , b分 别代表映射后的图像灰度级的分界点，即为a映射到a， c映射到c， d映 射到d ， b映射到b'， f(x，y)为原始图像的灰度分布函数公式中可用f简单 代替，f (x， y)为变换后的图像的灰度分布函数，线性变换的初始值a， a为0; 变换的第二个值c为空洞的最大灰度，令c-0;即把空洞灰度级压缩为一 个灰度级0;对以上线性变换的通式取b = b = d =d = 255 ,此时 3式略去仅有1式和2式进行线性变换；本铝泡样品中取c = d= 90， c =d
=0 5 b = b = 255, a = a = 0 。
步骤2,选取图像中的目的区域；
步骤3,将目的区域数据载入新的内存区域在第一张断层图像中用鼠 标选取一个需要的矩形块，以象素为单位，通过矩形块的四个坐标(xl,yl)(x2， yl)( xl，y3)( x2 ， y3)分别代表矩形的左上，右上，左下，右下角的坐标，设 原始图像的宽为width,高为height,断层图像的层数为pieces,原数据区
指针为pData,令i = 0 pieces—l; j = yi~ y3; k = x ' x2那么新载入的数据即 为pData[i*height*width+j*width+kl;
步骤4，对新的内存区域的数据运用八叉树的结构进行管理分配每个
八叉树单元的基本结构包含数据地址指针pBit,指向父单元的指针pParent，
指向子单元的指针数组pChildren[81， pChildren[8]有八个元素，两个标志
位分别是表示是否被遍历的标志M、表示是否满足空洞条件的标志S,载入
11图像的宽newwidth =x2 —xl,髙度newheight =y3 —yl，层数为pieces, 取newwidth、 newheight、 pieces中的最大值设为max; 令max《2'" + l ，
解得m的最小值；把体数据设为每边有2'" + l个体素的立方体，新载入内
存的数据从这个体的左上部分放起，体数据每条边分为2'"段，整个体数据分
为8'"个单元，即相邻的8个体素形成一个基本单元，8个相邻的子单元合成 一个父单元；每个单元拥有数据地址pBit,指向单元的左上角的体素的地址； 因为最低层单元没有子单元，所以指向子单元的指针数组pChildren[81的八 个元素指针为空，其他层级单元的pChildren[8的八个元素指针按照从左上 到右下的指向8个子单元；若level0层的小单元的8个体素满足是空洞的域 值要求或是非levelO级的所有子单元满足空洞条件，设置标志位S-true，否 则设置S-false;M-false,表示这个单元未遍历，M4rue表示遍历了 ； pParent 表示指向父单元的父指针。
步骤4.1:分配最底层，也即从leve10层的小单元开始设置，将leve10 层小单元的指向子单元的指针数组pChildren[81的八个元素指针设置为空， 标志位M=falSe，标志位S的设置通过判断8个顶点的灰度值是否全为0来 设置，若是，则S-true，若否，则置S-false;
步骤4.2:接着进行上一层即levell的设置，以八个相邻的小单元为子 单元合成一父单元，即新的正方体，设置这个父单元的pChildren81的值分 别为从左上到右下指向八个子单元的指针，标志位M-false,标志位S的设 置通过判断八个子单元是否都满足空洞条件，即八个子单元是否都满足S = true,若是，则设置标志位S = true;若否，则设S-false;把八个子单元的 父指针pParent指向这个父单元；
步骤4.3，进行level2的设置，以八个相邻的父单元为子单元再组成一 个父单元，参数设置方法同步骤4.2中的设置，依此类推，直至设置到levelm;
步骤5,从八叉树的levelO层开始对八叉树进行遍历；
步骤6，从leve10的第一个小单元开始向后査找；
步骤7,判断此单元是已经否遍历过，若是，即M二true，则回到步骤 6;若否，M = falSe，则进行步骤8;
12步骤8，判断此单元是否完全是空洞单元，若否，即S = falSe，则回到 步骤6;若是，则此单元同时满足M-false并且S-true，首先定义体积变 量volume = 0;进行步骤9;
步骤9,向此单元的父单元上溯；
步骤IO,查看父单元是否是空洞单元；
步骤IO.I，若满足条件M = falSe,并且S-true,则回到步骤9继续向 上上溯；
步骤10.2,若否，即S-false时，停止上溯，设停止单元为leveln级的 单元T，其中1《n《m (m为最大级数)，退回到level(n-l)级找到上溯到它 的那个子单元Tu ， 0《u《7，进行步骤ll;
步骤11，查看Tu是否是可分的；
步骤12，若2《n，贝lj Tu可分，计算volume,令volume = volume+Tu ; 令Tu的所有子单元直到最底层单元的标志位M = trUe，然后退到Tu的下级 位于level(n-2)级的八个子单元；
步骤13，分别查看位于level(n-2)的八个子单元Tuv(0《v《7)的每个的 三个相邻外部同级单元，此处每个单元只有三个同级不共父单元的相邻单 元。设某一个子单元的三个相邻外部同级单元为Tuva、 Tuvb、 Tuvc;
步骤14，判断Tuva、 Tuvb、 Tuvc的遍历标志M是否为true;
步骤14.1若M = true则结束对Tuva、 Tuvb、 Tuvc査找；
步骤14.2若M = false，则进行步骤15;
步骤15，分别判断Tuva、 Tuvb、 Tuvc的标志位S是否为true; 步骤15.1,如果是，则回到步骤9，即Tuva、 Tuvb、 Tuvc各自向上上 溯到level(n-l)级査看父单元；
步骤15.2,如果否，即S = false,则进行步骤16; 步骤16,判断该单元是否可分； 步骤16.1，若可分，则进行步骤17; 步骤16.2,若不可分，则进行步骤21;
步骤17，将与该单元相邻的单元Tuv八分，记为Tuvw(0《w《7)("该 单元"为T,、 Tuvb、 Tuvc中的一个)；
步骤18,分别査看与该单元(Tuva、 Tuvb、 Tuvc中的一个)相邻的四
13个八分单元Tuvw(w取0—7中的4个数)的三个相邻外部同层子单元(即 非同父单元的同层单元)Tuvwa， Tuvwb， Tuvwc;
步骤19,判断三个相邻同层子单元Tuvwa， Tuvwb， Tuvwc是否遍历过， 如果遍历过则终止这个单元的计算；若否，则进行步骤20;
步骤20，判断Tuvwa， Tuvwb， Tuvwc是否为空洞单元；
步骤20.1,若是，令Tu -Tuvwa或Tuvwb或Tuvwc,转到步骤ll; 此时经过步骤ll的判断，若Tu单元不可分，即11=1,则进行步骤21;
步骤20.2,若否，停止这个单元的计算；
步骤21, volume = volume + Tu，把Tu的遍历标志位M置为true,表 示遍历过，査看Tu周围六个同级单元；
步骤22,判断同级单元中是否有未遍历空洞单元，即M-false且S-true的单元；
步骤22.1，若有某一个或几个单元满足M = false且S = true的条件则 转到步骤9向这些单元的父单元上溯；
步骤22.2， Tu周围六个同级单元都不满足M = false且S-true，则终止 Tu的继续运算；
步骤23，结束点，得到第一个独立的空洞的体积；
步骤24，继续对levelO层的其余单元遍历，找到leve10的下一个满足 M = false且S = true的单元，令volume- 0，重复步骤6以下计算空洞体 积的步骤，直至体数据的所有levelO层的单元被遍历，得到目的区域的所有 泡的体积。
本发明方法的基础是建立在空洞的不连续性上。由于空洞的不连续，所 以每个小单元的扩充必将遇到周围的实质象素的阻拦，当然就可以构出完整
的空洞的体积以及位置。
上述实例只为说明本发明的技术构思及特点，其目的在于让熟悉此项技 术的人是能够了解本发明的内容并据以实施，并不能以此限制本发明的保护 范围。凡根据本发明精神实质所做的等效变换或修饰，都应涵盖在本发明的 保护范围之内。
权利要求
1. 一种基于图像的三维物体不连续空洞体积的测定方法，包括以下步骤将原始图像载入图像处理装置中并进行预处理原始图像；预处理后的图像根据原始图像的灰度级确定边界；撷取图像目的区域数据载入内存，在内存中构造空间八叉树通过遍历八叉树节点累积确定不连续空洞体积。
2、 根据权利要求1所述的基于图像的三维物体不连续空洞体积的测定方法，其特征在于所述原始图像是通过MicroCT系统三维扫描后得到的三维可视图像。
3、 根据权利要求1所述的基于图像的三维物体不连续空洞体积的测定方法，其特征在于所述的原始图像预处理是将原始图像经灰度线性变换处理。
4、 根据权利要求3所述的基于图像的三维物体不连续空洞体积的测定方法，其特征在于所述的灰度线性变换处理方式是通过灰度分布函数映射得到，所述的灰度分布函数为<formula>formula see original document page 2</formula>其中，a， c， d， b代表原始图像灰度级的四个分界点，a'， c', d'， b'分别代表映射后的图像灰度级的分界点；f'(x，y)为变换后的图像的灰度分布函数；f(x，y)为原始图像的灰度分布函数且a《f(X,y)<c， c《f(x，y)<d , d《f(x,y)《b。
5、 根据权利要求1所述的基于图像的三维物体不连续空洞体积的测定方法，其特征在于所述的图像目的区域数据是先撷取原始图像矩形块，通过原始图像矩形块的坐标与原始图像的层数获得。
6、 根据权利要求1所述的基于图像的三维物体不连续空洞体积的测定方法，其特征在于所述的空间八叉树选自规则八叉树、线性八叉树或一对八式八叉树。
7、 根据权利要求1所述的基于图像的三维物体不连续空洞体积的测定方法，其特征在于所述的图像目的区域数据是通过空间八叉树进行管理分配，所述的空间八叉树基本结构包括数据地址指针，指向父单元的指针，指向子单元的指针数组，所述指向子单元的指针数组包括表示是否被遍历的标志位、表示是否满足空洞条件的标志位，载入图像的宽度，载入图像的髙度，载入图像的层数；按照上述的结构分配设置空间八叉树各单元。
8、 根据权利要求1所述的基于图像的三维物体不连续空洞体积的测定方法，其特征在于所述的空间八叉树遍历包括如下步骤从八叉树第0层开始遍历，第O层第一个小单元开始向后査找；判断该单元是否遍历过和该单元是否完全是空洞单元；当该单元未曾遍历且完全是空洞单元时向父单元上溯；判断父单元是否是空洞单元和子单元是否可分；当父单元不是空洞单元且子单元可分时，退级到子单元，査看子单元的八个次级子单元，再査看每个子次级单元的周边的三个同层子单元；判断这三个同层子单元的遍历情况和是否满足空洞的条件；当同层子单元非空洞且该单元可分时将与同层子单元相邻的次级子单元八分，查看四个八分单元的3个相邻外部同层子单元判断该单元是否遍历过和该单元是否完全是空洞单元；当相邻同层子单元非遍历过且不是空洞单元，查看子单元是否可分，如此循环当循环到最后一级不可分单元的周围六个同级单元全部不属于空洞单元或者已遍历过时结束遍历。
9、 根据权利要求1所述的基于图像的三维物体不连续空洞体积的测定方法，其特征在于所述的空间八叉树遍历时，当子单元的父单元是空洞单元时继续上溯。
10、 根据权利要求1所述的基于图像的三维物体不连续空洞体积的测定方法，其特征在于所述的空间八叉树遍历时，当单元不满足空洞条件且单元不可分时，停止这个单元向周围扩展。
全文摘要
本发明公开了一种基于图像的三维物体不连续空洞体积的测定方法，包括以下步骤将原始图像载入图像处理装置中并进行预处理原始图像；根据预处理后的图像根据原始图像的灰度级确定边界；撷取图像目的区域数据载入内存，在内存中构造空间八叉树通过遍历八叉树节点累积确定不连续空洞体积。本发明在通过构造八叉树应用到MicroCT系统中三维可视化后处理中的有孔实体的孔体积的测量中，不但可以很精确的定位空洞的位置，而且可以很方便的实现对MicroCT系统中三维可视化后处理中的有孔实体的空洞体积的精确的测量。
文档编号G01B21/00GK101482401SQ20081013666
公开日2009年7月15日 申请日期2008年12月24日 优先权日2008年12月24日
发明者刘俊秀, 孙树文, 蔚 张, 光 李, 罗守华 申请人:苏州和君科技发展有限公司;东南大学