CN103700136A - 一种利用三变量双调和B-spline函数进行医学体数据矢量化的方法 - Google Patents
一种利用三变量双调和B-spline函数进行医学体数据矢量化的方法 Download PDFInfo
- Publication number
- CN103700136A CN103700136A CN201310634798.0A CN201310634798A CN103700136A CN 103700136 A CN103700136 A CN 103700136A CN 201310634798 A CN201310634798 A CN 201310634798A CN 103700136 A CN103700136 A CN 103700136A
- Authority
- CN
- China
- Prior art keywords
- function
- biharmonic
- super voxel
- phi
- spline
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Image Generation (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明提出一种利用三变量双调和B-spline函数进行医学体数据矢量化的方法,提出一种基于二次规划的方法构造三变量双调和B-spline基函数,本发明采用针对专门的函数集合计算离散拉普拉斯算子,最小化误差函数。并分析了该基函数近似满足局部条件和单位分解(Partition of Unity)条件。基于双调和B-spline函数,本发明设计了一种新的体数据矢量化方法。采用隐函数作为矢量化表示,使用线性规划方法优化该隐函数,完成三维医学体数据矢量化。
Description
技术领域
本发明涉及一种利用三变量双调和B-spline函数进行医学体数据矢量化的方法。
背景技术
以前已有很多学者提出多种离散数据插值方法。半径基函数是一类常用的基函数(Buhmann M.D.,Buhmann M.D.:Radial Basis Functions.Cambridge University Press,NewYork,NY,USA,2003.),如高斯函数,薄板样条等等。许多情况下,每个基函数是相互独立的,这样无法满足单位分解条件(即所有基函数的和恒等于1)。样条函数,如三次样条,多立方体样条(polycube spline)(Wang H.,He Y.,Li X.,Gu X.,Qin H.:Polycubesplines.Comput.Aided Des.40,6(June2008),721–733.),Voronoi样条(Voronoi spline)(Mirzargar M.,Entezari A.:Voronoi splines.IEEE Transactions on Signal Processing58,9(2010),4572–4582.)等能够表示曲面。Feng等(Feng P.,Warren J.:Discrete bi-laplacians andbiharmonic b-splines.ACM Trans.Graph.31,4(July2012),115:1–115:11.)注意到有限差分和拉普拉斯算子之间的关联,提出了双变量调和B-spline,并证明它是局部的并满足单位分解。但是由于三维空间中复杂的邻域结构,直接将该方法推广到三维并不适用。另一方面,三维空间中的双调和B样条函数的性质还有待研究。
以前有许多工作研究离散拉普拉斯算子,Meyer(Meyer M.,Desbrun M.,Schr¨oder P.,Barr A.H.:Discrete differential-geometry operators for triangulated2-manifolds.In Proc.VisMath(2002),pp.35–57.)等提出二维流形上的余切作为权重的离散拉普拉斯算子。Wardetzky等(Wardetzky M.,Mathur S.,K¨alberer F.,Grinspun E.:Discrete laplace operators:no free lunch.In Proceedings of the fifth Eurographics symposium on Geometry processing(2007),SGP’07,Eurographics Association,pp.33–37.)综述了曲面上离散拉普拉斯算子的构造方法。对于双拉普拉斯算子,许多方法使用迭代拉普拉斯方法构造,但是,在不规则网格中,他们的精确度受到限制,Feng等通过加入三次约束,将近似精度提高到三次,但在三维空间中,这依然不够精确。本发明提出一种具体函数集相关的离散拉普拉斯算子计算方法,专门针具体的函数集合进行优化。
发明内容
本发明解决的技术问题是:提出了一种三元双调和B-spline函数,满足局部性和单位分解,能够满足离散点插值等应用。基于该函数,本发明矢量化三维医学体数据。
本发明采用的技术方案为:一种利用三变量双调和B-spline函数进行医学体数据矢量化的方法,包括以下四个步骤:
步骤(1)构造虚拟结点:在超体素外,增加虚拟结点,构造边界点的邻接结点,辅助拉普拉斯算子离散化;
步骤(2)计算离散拉普拉斯算子:针对给定的函数集合,使用二次规划方法,计算最优化的离散拉普拉斯算子;
步骤(3)计算离散双拉普拉斯算子:在步骤(2)的基础上,针对给定的函数集合,使用二次规划方法,计算最优化的离散双拉普拉斯算子;
步骤(4)矢量化医学体数据:在前3步的基础上,本发明矢量化三维医学体数据。
本发明的原理在于:
(1)在超体素外采集更多结点并用Voro++(Rycroft C.H.:Voro++:A three-dimensionalvoronoi cell library in c++.Chaos:An Interdisciplinary Journal of Non-linear Science19,4(2009),041111.)计算三维Voronoi剖分,然后计算采样结点及其一环邻域的并集得到如果包含边界点,则采样更多的虚拟结点,重新计算Voronoi图和直到不包含边界点为止。
其中δij=1如果i=j,否则δij=0,H为所有列向量h构成的矩阵,Hkj表示矩阵H的第i行第j列的元素。
(4)使用如下线性规划优化矢量化超体素形状。
使得
-ej≤nj-n′j≤ej,
ej≤λ1|n′j|,j=1,...,dim(ni),
本发明与现有技术相比的有点在于:
1、本发明通过增加虚拟结点,支持边界结点离散化,从而在任意超体素上计算离散拉普拉斯和双拉普拉斯算子。
2、本发明采用二次规划方法针对专门的函数集合优化离散拉普拉斯算子,从而达到很高的精度。
3、在拉普拉斯算子基础上,针对专门的函数集合,本发明采用二次规划方法,计算离散双拉普拉斯算子,实验证明可以达到比现有方法更好的精确度。
4、基于双调和拉普拉斯算子,矢量化三维医学体数据,能够更好的矢量化三维医学数据。
附图说明
图1基函数对比(a)Feng的方法,(b)我们的方法;
图2基函数之和对比(a)Feng的方法,(b)我们的方法;
图3矢量化超体素形状(a)原始图像,(b)矢量化后结果;
图4矢量化超体素亮度(a)原始图像,(b)矢量化后结果;
图5体绘制效果(a)原始图像,(b)矢量化后结果。
具体实施方式
1.双变量双调和B-spline回顾
其中tj是第j个结点,表示tj的一环邻域。是结点tj的维诺区域,vij是连接结点ti和tj的边,eij是垂直于vij的Voronoi图的边。
类似的,双调和方程Δ2φy(x)=δ(x-y)在二维平面上的解为 类似的格林公式的双调和形式为:
2.三变量B-spline和它的性质
在三维欧几里得空间中,散度定理为:
本发明用▽φy(x)代替F,得到:
因此
更进一步,如果用▽(Δφy(x))替换F,得到:
本发明称三变量双调和B-spline的特征函数为:
3.虚拟结点
在一个闭合超体素(区域)中Ω中,采样结点本发明根据Voronoi图的邻接关系确定结点之间的邻接关系。一个结点ti为边界点,如果它的Voronoi单元和超体素区域边界相交。离散的拉普拉斯算子一般通过结点的邻域计算,但边界附近的结点缺少邻域点,导致边界精度计算明显下降,为了解决这些问题,本发明引入虚拟结点来辅助离散微分算子计算。
本发明在超体素区域外采集更多结点并用Voro++(Rycroft C.H.:Voro++:Athree-dimensional voronoi cell library in c++.Chaos:An Interdisciplinary Journal of Non-linearScience19,4(2009),041111.)计算三维Voronoi剖分,然后计算和的一环邻域的并集得到如果包含边界点,则采样更多的虚拟结点,重新计算Voronoi图和直到不包含边界点为止。
4.离散拉普拉斯算子
利用虚拟结点,本发明可以离散特征函数,提出一种针对某一函数集合优化的离散拉普拉斯算子,离散拉普拉斯算子的一个重要性质是系数之和为0。函数p(x)在结点tj离散化为其中h为待计算的系数向量,hj表示h的第j个元素。为了求解函数集合在结点tj处的
5.离散双调和算子
其中δij=1如果i=j,否则δij=0。将前面所计算的离散拉普拉斯算子代入得到离散双拉普拉斯算子,从而计算出三变量双调和B-spline函数。
6.三维体数据矢量化
对于输入的三维体数据(CT,MRI),所述数据由CT或MRI扫描设备获得数据,首先将它分割得到超体素,然后使用双调和B-spline将其矢量化。假设分割得到h个超体素,本发明使用隐函数来表示矢量超体素。如果当且仅当Fk(x)>Fi(x),1≤i≤h,i≠k。对于每个Fi(x),使用双调和B-spline将它分解为Fi(x)=Φi(x)Hini,其中为个采样点。
初始化:当时,令Fi(x)近似等于1,否则近似等于。所以对于采样点{s1,s2,...,sr},解如下线性最小二乘问题:
使得
|nj-n′j|≤λ1|n′j|,j=1,...,dim(ni),
其中nj表示ni的第j个元素,nj'为nj的初始值。以上优化问题可转化为如下线性规划求解:
使得
-ej≤nj-n′j≤ej,
ej≤λ1|nj′|,j=1,...,dim(ni),
颜色矢量化:超体素的形状矢量化后,矢量化超体素的颜色。对于每个超体素i,本发明优化如下能量函数:
其中c(xj)是点xj的亮度。
Claims (6)
1.一种利用三变量双调和B-spline函数进行医学体数据矢量化的方法,所述医学体数据为CT或MRI图像,其特征在于包括以下四个步骤:
步骤(1)构造虚拟结点:在超体素外,增加虚拟结点,构造边界点的邻接结点,辅助拉普拉斯算子离散化;
步骤(2)计算离散拉普拉斯算子:针对给定的函数集合,使用二次规划方法,计算最优化的离散拉普拉斯算子;
步骤(3)计算离散双拉普拉斯算子:在步骤(2)的基础上,针对给定的函数集合使用二次规划方法,计算最优化的离散双拉普拉斯算子;
步骤(4)将体数据分割成超体素,使用三变量双调和B-spline,矢量化每个超体素的形状;
步骤(5)使用三变量双调和B-spline,矢量化每个超体素颜色。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310634798.0A CN103700136B (zh) | 2013-12-01 | 2013-12-01 | 一种利用三变量双调和B‑spline函数进行医学体数据矢量化的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310634798.0A CN103700136B (zh) | 2013-12-01 | 2013-12-01 | 一种利用三变量双调和B‑spline函数进行医学体数据矢量化的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103700136A true CN103700136A (zh) | 2014-04-02 |
CN103700136B CN103700136B (zh) | 2017-04-12 |
Family
ID=50361653
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310634798.0A Active CN103700136B (zh) | 2013-12-01 | 2013-12-01 | 一种利用三变量双调和B‑spline函数进行医学体数据矢量化的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103700136B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107808387A (zh) * | 2017-11-13 | 2018-03-16 | 湖北工业大学 | 一种医学图像序列中的目标跟踪方法 |
CN110335275A (zh) * | 2019-05-22 | 2019-10-15 | 北京航空航天大学青岛研究院 | 一种基于三变量双调和b样条的流体表面的时空向量化方法 |
CN111260724A (zh) * | 2020-01-07 | 2020-06-09 | 王伟佳 | 一种基于周期b样条的实例分割方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100166280A1 (en) * | 2008-12-26 | 2010-07-01 | Canon Kabushiki Kaisha | Three-dimensional shape data processing apparatus and three-dimensionalshape data processing method |
CN102289840A (zh) * | 2011-06-18 | 2011-12-21 | 浙江大学 | 面向色盲设计颜色传输函数用于体绘制的方法 |
CN103345774A (zh) * | 2013-07-17 | 2013-10-09 | 中国人民解放军第三军医大学 | 一种三维多尺度矢量化的建模方法 |
-
2013
- 2013-12-01 CN CN201310634798.0A patent/CN103700136B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100166280A1 (en) * | 2008-12-26 | 2010-07-01 | Canon Kabushiki Kaisha | Three-dimensional shape data processing apparatus and three-dimensionalshape data processing method |
CN102289840A (zh) * | 2011-06-18 | 2011-12-21 | 浙江大学 | 面向色盲设计颜色传输函数用于体绘制的方法 |
CN103345774A (zh) * | 2013-07-17 | 2013-10-09 | 中国人民解放军第三军医大学 | 一种三维多尺度矢量化的建模方法 |
Non-Patent Citations (3)
Title |
---|
KEXIANG WANG等: "Restricted Trivariate Polycube Splines for Volumetric Data Modeling", 《IEEE TRANSACTIONS ON VISUALIZATION AND COMPUTER GRAPHICS》 * |
张显全 等: "基于控制点的二值图像边界矢量化算法", 《计算机科学》 * |
杨长强 等: "利用B样条实现基于等高线的三维地形图", 《信息技术与信息化》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107808387A (zh) * | 2017-11-13 | 2018-03-16 | 湖北工业大学 | 一种医学图像序列中的目标跟踪方法 |
CN107808387B (zh) * | 2017-11-13 | 2021-04-06 | 湖北工业大学 | 一种医学图像序列中的目标跟踪方法 |
CN110335275A (zh) * | 2019-05-22 | 2019-10-15 | 北京航空航天大学青岛研究院 | 一种基于三变量双调和b样条的流体表面的时空向量化方法 |
CN110335275B (zh) * | 2019-05-22 | 2023-03-28 | 北京航空航天大学青岛研究院 | 一种基于三变量双调和b样条的流体表面的时空向量化方法 |
CN111260724A (zh) * | 2020-01-07 | 2020-06-09 | 王伟佳 | 一种基于周期b样条的实例分割方法 |
CN111260724B (zh) * | 2020-01-07 | 2023-05-19 | 中山仰视科技有限公司 | 一种基于周期b样条的实例分割方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103700136B (zh) | 2017-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chougrani et al. | Lattice structure lightweight triangulation for additive manufacturing | |
Löhner | Extensions and improvements of the advancing front grid generation technique | |
Buchegger et al. | Planar multi-patch domain parameterization via patch adjacency graphs | |
Forrester et al. | Raney distributions and random matrix theory | |
CN102033985A (zh) | 基于*-矩阵算法的高效时域电磁仿真方法 | |
CN104966317A (zh) | 一种基于矿体轮廓线的三维自动建模方法 | |
Fries et al. | Higher-order meshing of implicit geometries, Part II: Approximations on manifolds | |
Gunderman et al. | Spectral mesh-free quadrature for planar regions bounded by rational parametric curves | |
CN103700136A (zh) | 一种利用三变量双调和B-spline函数进行医学体数据矢量化的方法 | |
Wan et al. | How to determine the minimum number of fuzzy rules to achieve given accuracy: a computational geometric approach to SISO case | |
Acosta | High order surface radiation conditions for time-harmonic waves in exterior domains | |
Huang et al. | HomPINNs: Homotopy physics-informed neural networks for learning multiple solutions of nonlinear elliptic differential equations | |
Aubry et al. | A three-dimensional parametric mesher with surface boundary-layer capability | |
Fu et al. | Low-latency hardware implementation of high-precision hyperbolic functions sinh x and cosh x based on improved cordic algorithm | |
CN110414016A (zh) | 超高速管道运输工具的乘波体外形参数化设计方法及系统 | |
Carr et al. | Representing interpolant topology for contour tree computation | |
CN104102471A (zh) | 一种应用fpga定点化技术实现指数cordic算法收敛域扩张的方法 | |
Zhang et al. | An adaptive interpolation element free Galerkin method based on a posteriori error estimation of FEM for Poisson equation | |
Jain et al. | A kernel class allowing for fast computations in shape spaces induced by diffeomorphisms | |
Liu et al. | An improved Poisson surface reconstruction algorithm based on the boundary constraints | |
Marcon et al. | A semi-structured approach to curvilinear mesh generation around streamlined bodies | |
Bibby et al. | Quadratic-linear duality and rational homotopy theory of chordal arrangements | |
Fayolle et al. | Rounding, filleting and smoothing of implicit surfaces | |
Nejati et al. | A fast hybrid approach for approximating a thin-plate spline surface | |
Liu et al. | Convergence analysis for the Chebyshev collocation methods to Volterra integral equations with a weakly singular kernel |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |