SuperEllipseTrace3D

From FarsightWiki

(Difference between revisions)

Revision as of 00:44, 11 May 2009

This page presents automated methods for robust modeling and analysis of 3-D vascular images using a \textit{cylindroidal superellipsoid} to model localized segments of vasculature. The proposed vessel model has an explicit, low-order parameterization, allowing for joint estimation of boundary and centerline information, thereby approximating the medial axis of tubular structures. Further, this explicit parameterization provides a geometric framework for traversing vessels in a directed manner. Topological information like branch point location and connectivity can also be detected as a side effect. An M-estimators provide robust region-based statistics that are used to drive the superellipsoid toward a vessel boundary. The proposed methodology behaves quite well across scale-space, shows a high degree of insensitivity to adjacent structures and implicitly handles branching.

Local vascular modelling using Superellipse

The superellipsoid provides the basis for a simple shape model that describes localized vessel segments. It has an explicit low-order parameterization relating directly to the local pose of a vessel segment. Using this geometric information, it is straightforward to traverse entire vascular networks.Also, since the superellipsoid model furnishes joint estimates of the boundary and centerpoint of a vessel, the resulting traversal approximates the medial axis of each vessel in the vasculature. Another key advantage when using the superellipsoid model is that its local nature allows us to adapt vessel intensity estimates across different regions of the vasculature.

Also, since the superellipsoid model furnishes joint estimates of the boundary and centerpoint of a vessel, the resulting traversal approximates the medial axis of each vessel in the vasculature. Another key advantage when using the superellipsoid model is that its local nature allows us to adapt vessel intensity estimates across different regions of the vasculature.

Figure1: Shows the boundary of a superellispoid while the shape parameter

ε 1

is varied between 1.0 and 0.25 and

ε 2

is set to 1.

Superellipsoids are geometric modeling primitives first introduced in early computer vision literature [2]. Superellipsoids are expressed implicitly as:

$\left(\left|x\right|^{\frac{2}{\epsilon_{2}}}+\left|y\right|^{\frac{2}{\epsilon_{2}}}\right)^{\frac{\epsilon_{2}}{\epsilon_{1}}}+\left|z\right|^{\frac{2}{\epsilon_{1}}}=1,$

where the parameters $\epsilon=\left(\epsilon_{1},\epsilon_{2}\right)$ control the shape of the superellipsoid. To understand the role of the shape parameters, first note that with $ε 1 = ε 2 = 1$ equation 1 is simply an ellipsoid. With $\epsilon_{1}\leq$ and $ε 2 = 1$ the superellipsoid becomes a cylindroid, i.e., the implicit surface bounds a convex region, where for constant $z$ , the 2-D level curves are elliptical. Setting $ε 1,ε 2$ to other values produces a variety of shapes including cuboids and pinched cuboids [2]. Since the goal in this work is to find an approximation to an elliptical cylinder, we will restrict our focus to cylindroidal forms. This means that $ε 2$ in equation is set to unity.

Starting with a set of points in the set $X=\left\{\mathbf{x}|S(\mathbf{x};\epsilon)=1\right\}$ ( $ε$ is fixed), a transformation given by $T(\mathbf{x})=\mathbf{R}({\phi})\mathbf{D}({\sigma})\mathbf{x}+{\mu}$ can be performed. This transformation preserves orthogonality of the local coordinate frame while allowing for anisotropic scaling. The action of $T(\mathbf{x})$ is to first scale $\mathbf{x}$ by $\mathbf{D}({\sigma})$ , which is a $3\times 3$ matrix with scale parameters ${\sigma}=\left(\sigma_{1},\sigma_{2},\sigma_{3}\right)$ on the diagonal. Then, a rotation is performed via the $3\times 3$ matrix $\mathbf{R}({\phi})$ with three degrees of freedom that is parameterized by $φ$ . Finally, a shift via ${\mu}=\left(\mu_{x},\mu_{y},\mu_{z}\right)^{T}$ is performed. This transformation leads to a new superellipsoid defined by the set having 10 degrees of freedom in 3D, $X'=\left\{\mathbf{x}|S(T^{-1}(\mathbf{x});\epsilon)=1\right\}$ , where, $T^{-1}(\mathbf{x})=\mathbf{D}({\sigma})^{-1}\mathbf{R}({\phi})^{T}(\mathbf{x}-{\mu})$

@@ Line 6: / Line 6: @@
 == Local vascular modelling using Superellipse ==
-The superellipsoid has an explicit low-order parameterization relating directly to the local pose of a vessel segment. Using this geometric information, it is straightforward to traverse entire vascular networks.Also, since the superellipsoid model furnishes joint estimates of the boundary and centerpoint of a vessel, the resulting traversal approximates the medial axis of each vessel in the vasculature. Another key
+The superellipsoid provides the basis for a simple shape model that describes localized vessel segments. It has an explicit low-order parameterization relating directly to the local pose of a vessel segment. Using this geometric information, it is straightforward to traverse entire vascular networks.Also, since the superellipsoid model furnishes joint estimates of the boundary and centerpoint of a vessel, the resulting traversal approximates the medial axis of each vessel in the vasculature. Another key
 advantage when using the superellipsoid model is that its local nature allows us to adapt vessel intensity estimates across different regions of the vasculature.
@@ Line 13: / Line 13: @@
-Superellipsoids are geometric modeling primitives first introduced
+[[Image:se_model.png |thumb|800px|'''Figure1:''' Shows the boundary of a superellispoid while the shape parameter <math> \epsilon_{1} </math> is varied between 1.0 and 0.25 and  <math> \epsilon_{2} </math> is set to 1.|center]]
-in early computer vision literature [2].
-[[Image:se_model.png |thumb|800px|'''Figure1:''' Shows superellipses with value of <math> \epsilon_{1} = 1.0, 0.75 </math> and <math>0.25 </math>.|center]]
-Superellipsoids are expressed implicitly as:
+Superellipsoids are geometric modeling primitives first introduced
+in early computer vision literature [2]. Superellipsoids are expressed implicitly as:
 <math>
@@ Line 25: / Line 24: @@
 </math>
-where the parameters <math>
+where the parameters <math> \epsilon=\left(\epsilon_{1},\epsilon_{2}\right) </math> control the shape of the superellipsoid. To understand the role of
-\epsilon=\left(\epsilon_{1},\epsilon_{2}\right) </math> control the
+the shape parameters, first note that with <math> \epsilon_{1}=\epsilon_{2}=1 </math> equation 1 is simply an ellipsoid. With <math>\epsilon_{1}\leq </math> and <math>\epsilon_{2}= 1</math> the superellipsoid becomes a cylindroid, i.e., the implicit surface bounds a convex region, where for constant <math>z</math>, the 2-D level curves are
-shape of the superellipsoid. To understand the role of
+elliptical. Setting <math>\epsilon_{1},\epsilon_{2}</math> to other values produces a variety of shapes including cuboids and pinched cuboids [2]. Since the goal in this work is to find an approximation to an elliptical cylinder, we will restrict our focus to cylindroidal forms. This means that <math>\epsilon_{2}</math> in equation is set to unity.
-the shape parameters, first note that with <math> \epsilon_{1}=
-\epsilon_{2}=1 </math> equation 1 is simply an ellipsoid. With
+Starting with a set of points in the set
-<math>\epsilon_{1}\leq </math> and <math>\epsilon_{2}= 1</math> the superellipsoid
+<math> X=\left\{\mathbf{x}|S(\mathbf{x};\epsilon)=1\right\}</math> (<math>\epsilon</math>
-becomes a cylindroid, i.e., the implicit surface bounds a
+is fixed), a transformation given by
-convex region, where for constant <math>z</math>, the 2-D level curves are
+<math>T(\mathbf{x})=\mathbf{R}({\phi})\mathbf{D}({\sigma})\mathbf{x}+{\mu}</math>
-elliptical. Setting <math>\epsilon_{1},\epsilon_{2}</math> to other values
+can be performed. This transformation preserves orthogonality of
-produces a variety of shapes including cuboids and pinched cuboids
+the local coordinate frame while allowing for anisotropic scaling.
-[2]. Since the goal in this work is to find an
+The action of <math>T(\mathbf{x})</math> is to first scale <math>\mathbf{x}</math> by
-approximation to an elliptical cylinder, we will restrict our
+<math>\mathbf{D}({\sigma})</math>, which is a <math>3\times 3</math> matrix with
-focus to cylindroidal forms. This means that <math>\epsilon_{2}</math> in
+scale parameters
-equation is set to unity.
+<math>{\sigma}=\left(\sigma_{1},\sigma_{2},\sigma_{3}\right)</math> on the
+diagonal. Then, a rotation is performed via the <math>3\times 3</math> matrix
+<math>\mathbf{R}({\phi})</math> with three degrees of freedom that is
+parameterized by <math>{\phi}</math>. Finally, a shift via
+<math>{\mu}=\left(\mu_{x},\mu_{y},\mu_{z}\right)^{T}</math> is performed.
+This transformation leads to a new superellipsoid defined by the
+set having 10 degrees of freedom in 3D,
+<math>X'=\left\{\mathbf{x}|S(T^{-1}(\mathbf{x});\epsilon)=1\right\}</math>,
+where,
+<math>
+T^{-1}(\mathbf{x})=\mathbf{D}({\sigma})^{-1}\mathbf{R}({\phi})^{T}(\mathbf{x}-{\mu})
+</math>

SuperEllipseTrace3D

Revision as of 00:44, 11 May 2009

Contents

Local vascular modelling using Superellipse

Views

Personal tools

Navigation

Search

Toolbox