主权项 |
一种基于图像的三维物体不连续空洞体积的测定方法,包括以下步骤:将原始图像载入图像处理装置中并进行预处理原始图像;预处理后的图像根据原始图像的灰度级确定边界;撷取图像目的区域数据载入内存,在内存中构造空间八叉树通过遍历八叉树节点累积确定不连续空洞体积;所述原始图像是通过MicroCT系统三维扫描后得到的三维可视图像;所述的原始图像预处理是将原始图像经灰度线性变换处理;所述的灰度线性变换处理方式是通过灰度分布函数映射得到,所述的灰度分布函数为 <mrow> <msup> <mi>f</mi> <mo>′</mo> </msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mfrac> <mrow> <msup> <mi>c</mi> <mo>′</mo> </msup> <mo>-</mo> <msup> <mi>a</mi> <mo>′</mo> </msup> </mrow> <mrow> <mi>c</mi> <mo>-</mo> <mi>a</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>a</mi> <mo>)</mo> </mrow> <mo>+</mo> <msup> <mi>a</mi> <mo>′</mo> </msup> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <msup> <mi>d</mi> <mo>′</mo> </msup> <mo>-</mo> <msup> <mi>c</mi> <mo>′</mo> </msup> </mrow> <mrow> <mi>d</mi> <mo>-</mo> <mi>c</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>c</mi> <mo>)</mo> </mrow> <mo>+</mo> <msup> <mi>c</mi> <mo>′</mo> </msup> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <msup> <mi>b</mi> <mo>′</mo> </msup> <mo>-</mo> <msup> <mi>d</mi> <mo>′</mo> </msup> </mrow> <mrow> <mi>b</mi> <mo>-</mo> <mi>d</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>d</mi> <mo>)</mo> </mrow> <mo>+</mo> <msup> <mi>d</mi> <mo>′</mo> </msup> </mtd> </mtr> </mtable> </mfenced> </mrow>其中,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;所述的图像目的区域数据是先撷取原始图像矩形块,通过原始图像矩形块的坐标与原始图像的层数获得;所述的空间八叉树选自规则八叉树、线性八叉树或一对八式八叉树;所述的图像目的区域数据是通过空间八叉树进行管理分配,所述的空间八叉树基本结构包括数据地址指针,指向父单元的指针,指向子单元的指针数组,所述指向子单元的指针数组包括表示是否被遍历的标志位、表示是否满足空洞条件的标志位,载入图像的宽度,载入图像的高度,载入图像的层数;按照上述的结构分配设置空间八叉树各单元。 |