Home Four-Dimensional Vascular Tree Reconstruction Using Moving Grid Deformation
Post
Cancel

Four-Dimensional Vascular Tree Reconstruction Using Moving Grid Deformation

Rationale and Objectives

Thinned perforator flaps have been widely used in plastic surgery for greater survivability and decreased morbidity. However, quantitative analysis of three-dimensional (3D) blood flow direction and location has not been examined yet. Such information will benefit and guide the surgical thinning and dissection process. Toward this goal, this study was performed for 3D vascular tree reconstruction with the incorporation of temporal contrast-agent propagation information (three spatial dimensions plus one temporal dimension; ie, 4D).

Materials and Methods

A novel computational framework by adopting a moving grid deformation method is presented. To take advantage of temporal information of the bolus propagating, a sequential segmentation procedure is proposed. Moreover, the temporal evolution of the vascular tree (4D vascular tree) is reconstructed during the procedure.

Results

Eight anterolateral thigh perforator flaps from eight cadavers were used for this study. The age range is 60–80 years old and the gender includes four males and four females. The 3D nature of the vascular structure and 4D vascular tree evolving process are showed in comparison with maximum intensity projection images.

Conclusion

The proposed computational framework demonstrates effectiveness in the modeling of 4D vascular tree. Furthermore, it reveals the ability to detect small vessel tree structures that are beyond the limit of image resolution.

Perforator flaps have been increasingly used over the past decade in plastic surgery for greater survivability and decreased morbidity ( ). The majority of anatomic vascular studies on perforator flaps have used lead oxide injections on cadaver flaps followed by two-dimensional (2D) projection radiography to determine vascular territories. Although lead oxide treated specimens provide excellent images, the three-dimensional (3D) nature of the vascular structure is not revealed. To overcome this deficit, static computed tomography (CT) scanning can be used to provide 3D mappings of contrast enhanced vascular anatomy in the sagittal, coronal, and transverse views. However, none of these methods can provide insight information on the direction and location of blood flow within each layer of a perforator flap. In this study, we present a four-dimensional (4D) (ie, three spatial dimensions plus one temporal dimension) time sequence of perforator flaps showing the process of contrast agent propagation. To our best knowledge, no study has been carried out reporting the examination of 4D anatomy of perforator flaps.

The computational techniques used in traditional 3D vasculature tree reconstruction ( ) can be roughly grouped into six categories ( ): 1) pattern recognition techniques; 2) model-based approaches; 3) tracking-based approaches; 4) artificial intelligence–based approaches; 5) neural network–based approaches; and 6) miscellaneous tubelike object detection approaches. Unfortunately, no single method works for any applications. The most widely used model-based methods in vessel reconstruction are the geometric deformable models. One problem with geometric deformable models is the high computational cost. To improve the efficiency and accuracy, adaptive grid techniques have been adopted ( ) to help solve the partial differential equations involved in geometric deformable models.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Materials and methods

Get Radiology Tree app to read full this article<

CT Radiography

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Moving Grid Generation by Deformation Method

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Step 1

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

divv→(xt,yt,t)=−∂∂t(1m(,t)), d

i

v

v

(

x

t

,

y

t

,

t

)

=

t

(

1

m

(

,

t

)

)

,

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

curlv→(xt,yt,t)=0, c

u

r

l

v

(

x

t

,

y

t

,

t

)

=

0

,

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

v→(xt,yt,t)⋅n→=0,for(xt,yt)∈∂Ω(t), v

(

x

t

,

y

t

,

t

)

n

=

0

,

f

o

r

(

x

t

,

y

t

)

Ω

(

t

)

,

where ∂Ω(t) ∂

Ω

(

t

) denotes the boundary of domain Ω(t) Ω

(

t

) . The Neumann boundary condition specifies the values the derivatives of a solution is to take on the boundary of the domain. In an image domain, the boundary constraint is zero which means the boundary points cannot move out or move inside to the image domain, but only move along the boundary. This is reflected as the zero vector field in the normal direction as shown in Eq 3 .

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

divv→(x,y,t)=∇⋅v→(x,y,t), d

i

v

v

(

x

,

y

,

t

)

=

v

(

x

,

y

,

t

)

,

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

curlv→(x,y,t)=∇×v→(x,y,t). c

u

r

l

v

(

x

,

y

,

t

)

=

×

v

(

x

,

y

,

t

)

.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Step 2

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

∂∂t(x,y,t)=m((x,y,t),t)v→((x,y,t),t), ∂

t

(

x

,

y

,

t

)

=

m

(

(

x

,

y

,

t

)

,

t

)

v

(

(

x

,

y

,

t

)

,

t

)

,

where t∈[0,T] t

[

0

,

T

] .

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

∫Ω(t)(−∂∂t(1m(,t)))dxdy=0 ∫

Ω

(

t

)

(

t

(

1

m

(

,

t

)

)

)

d

x

d

y

=

0

for all t . Furthermore, the monitor function m(,t) m

(

,

t

) should be normalized under the following condition,

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

∫Ω(t)1m(,t)dxdy=|Ω(0)|,forallt, ∫

Ω

(

t

)

1

m

(

,

t

)

d

x

d

y

=

|

Ω

(

0

)

|

,

for

all

t

,

where | Ω(0) | means the area over the domain Ω at time t = 0. The mathematical proof of the deformation method can be found elsewhere ( ).

Get Radiology Tree app to read full this article<

Vascular Tree Reconstruction using Moving Grid with Deformation

Get Radiology Tree app to read full this article<

Initial segmentation

Get Radiology Tree app to read full this article<

Vascular grid plot generation

Get Radiology Tree app to read full this article<

m0()={0.3,whenpoint(x,y)belongstoinitialsegmentedvesselregion;1.0,else m

0

(

)

=

{

0.3

,

when

point

(

x

,

y

)

belongs

to

initial

segmented

vessel

region

;

1.0

,

else

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

m(,t)=1−tT+tT⋅m0(), m

(

,

t

)

=

1

t

T

+

t

T

m

0

(

)

,

where t∈[0,T] t

[

0

,

T

] is the dummy time variable used in the grid deformation process same as the one occurred in Eq 1–7 . By this definition, we have m(,0)=1 m

(

,

0

)

=

1 and m(,T)=m0() m

(

,

T

)

=

m

0

(

) . This implies the grid is deformed from uniform to the desired one in which the size of a grid cell is governed by the objective monitor function m0() m

0

(

) . The normalization of m(,t) m

(

,

t

) is carried out using Eq 7 for all t∈[0,T] t

[

0

,

T

] .

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Vascular vessel extraction

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

v→T=<xT−x0,yT−y0>T, v

T

=

<

x

T

x

0

,

y

T

y

0

T

,

where (x0,y0) (

x

0

,

y

0

) is the initial uniform grid, and (xT,yT) (

x

T

,

y

T

) is the grid plot using the deformation method.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

F(xT,yT)=∇⋅v→T. F

(

x

T

,

y

T

)

=

v

T

.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Results

Get Radiology Tree app to read full this article<

Demonstration of the Proposed Method

Initial vascular tree segmentation

Get Radiology Tree app to read full this article<

Figure 1, Initial segmentation for a 50 × 50 pixel region.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Figure 2, Original images of a 50 × 50 pixel region at different inject times: (a) time t 1 ; (b) time t 2 ; (c) time t 3 ; (d) time t 4 ; (e) time t 5 ; and (f) time t 6 .

Figure 3, Difference images at adjacent injection times: (a) difference between Fig 2 a and 2 b; (b) difference between Fig 2 c and 2 b; (c) difference between Fig 2 c and 2 d; (d) difference between Fig 2 d and 2 e; and (e) difference between Fig 2 e and 2 f.

Get Radiology Tree app to read full this article<

Grid plot generation

Get Radiology Tree app to read full this article<

Figure 4, Grid generation for the corresponding region in Fig 1 : (a) initial uniform grid; (b) grid plot at iteration step 7; and (c) final generated grid plot.

Get Radiology Tree app to read full this article<

Vascular vessel extraction

Get Radiology Tree app to read full this article<

Figure 5, Vector field for the local region in Fig 1 .

Figure 6, Vascular region delineation: (a) divergence field and (b) final region segmentation.

Get Radiology Tree app to read full this article<

4D vascular tree evolution

Get Radiology Tree app to read full this article<

Figure 7, Temporal evolving of vasculature tree: (a) at injection time t 2 ; (b) at injection time t 3 ; and (c) at injection time t 6 .

Get Radiology Tree app to read full this article<

Results for Datasets ALT2-ALT8

Get Radiology Tree app to read full this article<

Figure 8, Final reconstructed three-dimensional (3D) vascular trees with respective maximum intensity projection (MIP) images: (a) MIP for dataset ALT2; (b) 3D reconstruction for dataset ALT2; (c) MIP for dataset ALT3; (d) 3D reconstruction for dataset ALT3; (e) MIP for dataset ALT4; (f) 3D reconstruction for dataset ALT4;

Get Radiology Tree app to read full this article<

Discussion

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

References

  • 1. Allen R.J., Treece P.: Deep inferior epigastric perforator flap for breast reconstruction. Ann Plastic Surg 1994; 32: pp. 32-38.

  • 2. Wei C., Jain V., Suominen S., Chen H.: Confusion among perforator flap: what is a true perforator flap?. Plastic Reconstr Surg 2001; 107: pp. 874-876.

  • 3. Abdul-Karim M.A., Roysam B., Dowell-Mesfin N.M., et. al.: Automatic selection of parameters for vessel/neurite segmentation algorithms. IEEE Trans Image Proc 2005; 14: pp. 1338-1350.

  • 4. Condurache A.P., Aach T., Grzybowski S., et. al.: Imaging and analysis of angiogennesis for skin transplantation by microangiography. Proc IEEE Conf Image Proc 2005; 2: pp. 1250-1253.

  • 5. Passat N., Ronse C., Baruthio J., et. al.: Cerebral vascular atlas generation for anatomical knowledge modeling and segmentation purpose. Proc Computer Vision Patteren Recognit 2005; 2: pp. 331-337.

  • 6. Sebbe R., Gosselin B., Coche E., et. al.: Model-guided segmentation of opacified thorax vessels. Proc IEEE Conf Image Proc 2005; 1: pp. 25-28.

  • 7. Sekiguchi H., Sugimoto N., Eiho S., et. al.: Blood vessel segmentation for head MRA using branch-based region growing. Syst Computers Jap 2005; 36: pp. 5.

  • 8. Kirbas C., Quek F.K.H.: Vessel extraction techniques and algorithms: a survey. Proc IEEE Symp BioInfor BioEng 2003; 4: pp. 238-245.

  • 9. Droske M., Meyer B., Schaller C., et. al.: An adaptive level set method for medical image segmentation. Proceedings of the International Conference on Information Processing in Medical Image 2001; 2080: pp. 416-422.

  • 10. Han X., Xu X., Prince J.P.: A 2D moving grid geometric deformable model. Proc Computer Vision Patteren Recognit 2003; 1: pp. 153.

  • 11. Xue J., Liao G.: Least-squares finite element method on adaptive grid for PDEs with shocks. Numer Meth Part Diff Eq 2006; 22: pp. 114-127.

  • 12. Zegeling P.A.: Moving grid techniques.Thompson J.F.Soni B.K.Weatherill N.P.1999.CRC Press LLCBoca Raton, FL:

  • 13. Terzopoulos D., Vasilescu M.: Sampling and reconstruction with adaptive meshes. Proc Comp Vision Patt Recognition 1991; 1: pp. 70-75.

  • 14. Vasilescu M., Terzopoulos D.: Adaptive meshes and shells. Proc Comp Vision Patt Recognition 1992; 1: pp. 829-832.

  • 15. Liao G., Anderson D.A.: A new approach to grid generation. Appl Anal 1992; 44: pp. 285-298.

  • 16. Moser J.: Volume elements of a Riemann manifold. Trans AMS 1965; 120: pp. 286-294.

  • 17. Semper B., Liao G.: A moving grid finite element method using grid deformation. Numer Meth Part Diff Eq 1995; 11: pp. 603.

This post is licensed under CC BY 4.0 by the author.