Translator Disclaimer
7 August 2013 Label-free optical lymphangiography: development of an automatic segmentation method applied to optical coherence tomography to visualize lymphatic vessels using Hessian filters
Author Affiliations +
Abstract
Lymphatic vessels are a part of the circulatory system that collect plasma and other substances that have leaked from the capillaries into interstitial fluid (lymph) and transport lymph back to the circulatory system. Since lymph is transparent, lymphatic vessels appear as dark hallow vessel-like regions in optical coherence tomography (OCT) cross sectional images. We propose an automatic method to segment lymphatic vessel lumen from OCT structural cross sections using eigenvalues of Hessian filters. Compared to the existing method based on intensity threshold, Hessian filters are more selective on vessel shape and less sensitive to intensity variations and noise. Using this segmentation technique along with optical micro-angiography allows label-free noninvasive simultaneous visualization of blood and lymphatic vessels in vivo. Lymphatic vessels play an important role in cancer, immune system response, inflammatory disease, wound healing and tissue regeneration. Development of imaging techniques and visualization tools for lymphatic vessels is valuable in understanding the mechanisms and studying therapeutic methods in related disease and tissue response.

1.

Introduction

The circulatory network in mammals is composed of a cardiovascular and lymphatic system responsible for delivering oxygen, nutrition, immune cells, and hormones to tissue and collecting waste materials from cells within living organs. These compounds are exchanged with cells via capillary beds where, due to the high blood pressure at arterioles, some part of plasma leaks into the interstitial space. After the exchange, the protein-rich interstitial fluid (lymph) is collected and returned back to the circulatory system by the lymphatic system. The lymphatic system consists of unidirectional, thin-walled capillaries, and a larger vessel network that drains lymph fluid from extracellular spaces within organs into larger collecting ducts. Finally, lymph is returned to venous circulation through the thoracic duct. The lymphatic system also includes lymphoid organs such as lymph nodes, tonsils, Peyer’s patches, spleen, and thymus which play a crucial part in immune response. The lymphatic system usually develops in parallel to the blood vessels in the skin and in most internal organs and is not present in the central nervous system, bone marrow, and avascular structures such as cartilage, epidermis, and cornea. Besides draining lymph fluid from extracellular spaces, other roles of the lymphatic system include absorbing lipids from the intestinal tract, maintaining the fluid hemostasis and transporting antigen-presenting cells and leukocytes to lymphoid organs. Also, the lymphatic system plays an important role in the development of several diseases such as cancer, lymphedema, and some inflammatory conditions.1,2

Although lymphatic vessels were first described in the early 17th century, they have been termed the forgotten circulation3 because their role and function was mostly ignored due to the lack of monitoring and imaging techniques. Several contrast-enhanced and labeled image modalities and imaging techniques including X-ray, computed tomography (CT), magnetic resonance imaging (MRI), ultrasound, and optical imaging have been used in clinical and animal studies for noninvasive and direct lymphatic imaging. In lymphangiography, an iodinated contrast agent such as blue dye,4 Iotasul (Schering AG, Berlin, Germany)5 or Iopamidol (Nihon Schering, Osaka, Japan)6 is either injected directly into a lymph vessel or indirectly into intradermal area for X-ray or CT. Although CT scans differentiated iodine concentration and lymph node morphology in normal and malignant lymph node, their application was limited due to their fast drainage from blood vessels that limited the imaging time in clinic. More recent development of bismuth sulfide nanoparticles CT agents with longer circulating half-life addressed some of the problems while their long-term toxicity is a major concern. Lymphoscintigraphy with 99m-technetium (TC-99m) is the most common nuclear method of imaging the lymphatics that enables visualization of the lymphatic network using single-photon emission computer tomography (SPECT).79 However, this technique requires rejection of scattered gamma photons and a large number of photons remain undetected and scanning times can be long. On the contrary, positron emission tomography with beta-emitting isotopes (18-F-fluro-2-deoxy-D-glucose, FDG) provides higher sensitivity specifically when combined with CT for better spatial resolution.10 MRI, unlike nuclear imaging techniques, does not require ionizing radiation exposure. MR lymphangiography is performed by intravenous or interstitial injection of gadolinium-labeled diethylene-triaminepentaacetic acid (Gd-DTPA), Gd dendrimetrs or liposomes and iron oxide particles.1113 Flourescence microlymphangiography is an optical imaging technique based on absorption of light in the visible or near-infrared (NIR) region by utilizing intradermal, subcutaneous, or intravenous injection of fluorescein contrast agents14,15 such as indocyanin-green (IC-Green)16 and fluorophore17 to map lymphatic vasculature. Also, molecular staining of lymph endothelium has been made possible by the discovery of lymph-specific molecular markers such as Prox-1, lymph endothelial-specific hyaluronan receptor 7 (LYVE-1),18 and vascular endothelial growth factor receptor 3 (Vegfr3)19 that were extensively used for imaging lymphatic vessels in development, wound healing, inflammation, and tumor metastasis.20 Very recently, Kalchenko et al.21 developed a label-free laser-speckle imaging technique for simultaneous visualization of blood and lymphatic vessels in the external mouse ear tissue. The origin of the observed lymph vessels was confirmed by intradermal injection of Dextran-FITC (0.5 M 10mg/ml) into the tissue followed by imaging using a multi-modal imaging system.22

Optical coherence tomography (OCT) has been widely used to noninvasively provide high-resolution, depth-resolved cross sectional and three-dimensional (3-D) images of highly scattering samples.2325 More recently, spectral-domain OCT and swept-source OCT imaging have been developed where these frequency-domain OCT (FD-OCT) methods have advantages over time-domain ones in terms of acquisition speed, sensitivity, and signal-to-noise ratio (SNR).26,27 In order to compensate for the complex conjugate ambiguity and acquisition speed which limited the practical applications of FD-OCT, real-time in vivo full-range complex FD-OCT has been developed to increase the ranging distance in the sample.28,29 OCT functional and hemodynamic information can be acquired in addition to the structure. Optical micro-angiography (OMAG) is a label-free noninvasive imaging and processing method used to obtain 3-D blood perfusion map in microcirculatory tissue beds in vivo using FD-OCT.30,31 Ultrahigh-sensitive OMAG (UHS-OMAG) is a variation of OMAG technique, capable of imaging microvasculature down to capillary level, in which data acquisition is based on repeated B-scan (frame) acquisition at the same spatial location,32,33 and then separating static scatterers (e.g., structure tissue) and dynamic scatterers (e.g., moving red blood cells within patent vessels) by the use of magnitude subtraction,34 phase-compensated subtraction,35 or eigendecomposition-(ED-) based36 OMAG algorithms. By synergistically utilizing both the amplitude and phase information of complex OCT signals, the OMAG technique is sensitive to both the axial and transverse blood flows, providing blood flow map of small vessels and capillaries as well as larger vessels. UHS-OMAG has been applied to visualize blood perfusion and micro-vasculature map in various living tissue samples such as retina,37,38 cerebral,39 renal microcirculation,40 and skin.41 Hemodynamic quantification can be extremely useful in cancer, stroke, and some other diseases which involve vasculature.

Since the lymph fluid is clear and transparent, lymphatic vessels appear as reduced scattering vessel-like areas in OCT structure cross section images. The origin of these low scattering transparent vessel-like structures has already been verified as lymph vessels by injecting Evan’s blue dye injection into the tissue and observing the uptake path by lymphatic vessels.42 These vessels can be visualized by applying a lower threshold on the intensity image.4244 However, the intensity-threshold technique does not take the physical shape of the vessels into account. Also, it is not robust to intensity variations and noise caused by speckle and light absorption in the structure. These limitations increase the chances of segmentation error and false alarms in identifying lymphatic vessels. In this paper, we propose an automatic filtering technique and vesselness model based on Hessian multi-scale vesselness filters45 for segmenting lymphatic vessel in OCT structure images. Hessian filters estimate tubular and vessel-like structures and the multi-scale approach allowed segmenting vessels with different size. Although similar techniques have been used in segmenting blood vessels in other image modalities such as MRI46 and CT,47 to the best of our knowledge, this is the first time that a systematic segmentation approach is presented to identify lymphatic vessels in OCT structure images. The advantage of using OCT over other image modalities is that it allows 3-D label-free noninvasive simultaneous imaging structure, blood vessels, and lymphatics.

2.

Proposed Method

2.1.

Problem Formulation

The local behavior of an image I(x,s) at scale s and location x can be expressed by its Taylor series expansion up to second order given by

Eq. (1)

I(x+δx,s)I(x,s)+δxT(I)s+δxTH(I)sδx,
where (.)s and H(.)s are the gradient vector and Hessian matrix of the image at scale s, respectively. Since our volume cross section images are discrete signals, finding their 2-D first-order and second-order derivatives can be ill-posed. Using the concepts of linear scale space theory, differentiation can be defined as a convolution with derivatives of a Gaussian:

Eq. (2)

xI(x,s)=sγI(x,s)*xG(x,s),
where γ is the derivative normalization parameter and

Eq. (3)

G(x,s)=12πs2ex22s2,
and . is the Euclidean norm operation. The second-order derivative can be expressed as

Eq. (4)

δxTH(I)sδx=(x)(x)I(x,s)=s2γI(x,s)*2x2G(x,s).

By analyzing the eigenvalues and eigenvectors of the Hessian matrix, the principal direction of the local structure can be extracted, which is the direction of the smallest curvature (along the vessel). By definition, the eigenvalues and vectors of the Hessian matrix are given by solving the following equations:

Eq. (5)

HSes,i=λs,ies,i,
where λs,i is the i’th eigenvalue at scale s corresponding to normalized eigenvector es,i and (|λ3||λ2||λ1|) for a 3-D structure.

For an ideal tubular structure, the following conditions hold:

Eq. (6)

|λ3|0,

Eq. (7)

|λ3||λ2|,

Eq. (8)

λ2λ1.

Based on the second-order ellipsoid, two geometric ratios are defined

Eq. (9)

RB=|λ1|λ2λ3
and

Eq. (10)

RA=|λ2||λ3|.

The first ratio (RB) accounts for the deviation from a blob-like structure but cannot distinguish between a line-like and a plate-like pattern. The second ratio refers to the largest cross section area of the ellipsoid that can distinguish between plate-like and line-like structures since it is zero in the latter case.

In order to distinguish background noise where no structure is present, Frobenius matrix norm is utilized. This measure is derived from eigenvalues of Hessian matrix by

Eq. (11)

RC=HF=jDλj2,
where D is the dimension of the image. This measure will be low in the background where no structure is present and the eigenvalues are small for the lack of contrast. In regions with high contrast compared to the background, the norm will become larger since at least one of the eigenvalues will be large.

The vesselness function at scale s is defined as

Eq. (12)

v0(s)={0ifλ2>0orλ3>0(1eRA22α2)eRB22β2(1eRC22θ2),
where α, β, and θ are thresholds which control the sensitivity of the line filter to the measures RA, RB, and RC. The idea behind this expression is to map the features into probability-like estimates of vesselness according to different criteria.

The vesselness measure is analyzed at different scales. The response of the line filter will be maximum at a scale that approximately matches the vessel size:

Eq. (13)

V0(γ)=argmaxsv0(s,γ),smin<s<smax,
where smin and smax are lower and upper bound in the range of scale (vessel sizes).

The vesselness measure for 2-D images (cross section) which follows directly from the same reasoning as in 3-D, is given by

Eq. (14)

v0(s)={0ifλ2>0eRB22β2(1eRC22θ2),
where RB=λ1λ2. Although, this measure cannot distinguish between plate-like and line-like structures, its implementation is more efficient and processing time is much faster than the 3-D measure. Therefore, we employ this measure for segmenting lymphatic vessels from OCT structure cross sections. The Appendix Appendix: contains the implementation detail of the proposed algorithm.

2.2.

Cylindrical Circular Model with Gaussian Cross Section

For a better demonstration of the model, we calculated the eigenvalues of the Hessian matrix for an ideal vessel structure shown in Fig. 1, which is a circular cylinder model with Gaussian-like intensity profile at the cross section. Basically, the model is a cylinder-shaped vessel with a Gaussian intensity profile at its cross section given by

Eq. (15)

I0(x,y,z)=C2σ02ex2+y22σ02,
where x, y, and z are Cartesian coordinates, σ0 is Gaussian parameter, and C is a normalization factor. Therefore, the Hessian matrix of the shape structure is given by

Eq. (16)

H=[d2I0dx2d2I0dxdyd2I0dxdzd2I0dyd_xd2I0dy2d2I0dydzd2I0dzdxd2I0dzdyd2I0dz2]=I0σ04[x2σ02xy0xyy2σ020000].

Fig. 1

(a) Cylindrical circular vessel model with Gaussian intensity profile at cross section (b).

JBO_18_8_086004_f001.png

Eigenvalues and eigenvectors of the Hessian matrix are

Eq. (17)

λ1=I0σ02,λ2=I0σ02(σ02(x2+y2)σ02)andλ3=0

Eq. (18)

v1=(y,x,0),v2=(x,y,0)andv3=(0,0,1),
that can be generalized to Eqs. (6)–(8) for an ideal tubular structure.

3.

Experimental Results and Discussions

3.1.

System Setup and Animal Preparation

To demonstrate the feasibility of our segmentation method, we used an OMAG imaging system to acquire the microvasculatural images within tissue beds in vivo, on which the algorithm is tested. The schematic of the OMAG/OCT imaging system is shown in Fig. 2(a), in which we used a superluminescent diode (SLD) with the center wavelength of 1310 nm and bandwidth of 65 nm, delivering an axial resolution of 12μm in the air. An optical circulator was used to couple the light from the SLD into fiber-based Michelson interferometer. In the reference arm of the interferometer, a stationary mirror was utilized after polarization controller. In the sample arm of the interferometer, a microscopy objective lens with 18-mm focal length was used to achieve 5.8μm lateral resolution. Then, a 2×2 optical fiber coupler was utilized to recombine the backscattered light from the sample and the reflected light from the reference mirror. Since the wavelength of the light source is invisible to the human eye, a 633-nm laser diode was used as a guiding beam to locate the imaging position. This reference helps to adjust the sample under the OMAG/OCT system and image the desired location. The recombined light was then routed to a home-built high-speed spectrometer via the optical circulator. The camera in the spectrometer had a line scan rate of 47,000lines/s. The spectral resolution of the designed spectrometer was 0.141nm that provided a detectable depth range of 3.0mm on each side of the zero-delay line. The measured SNR of the system was 105dB with a light power on the sample at 3mW.

Fig. 2

(a) Schematic diagram of the imaging system, where PC represents the polarization controller and CCD is the charged coupled device (camera). (b) Typical mouse ear pinna flattened on the imaging platform. The biopsy punch size is 1 mm. The rectangular box shows typical scanning size of about 2.2×2.2mm2.

JBO_18_8_086004_f002.png

Noninvasive in vivo images were acquired from the pinna of a healthy 6-weeks-old-male hairless mouse weighing approximately 25 g. During the experiment, the mouse was anesthetized using 2% isoflurance (0.2L/min O2, 0.8L/min air). The ear was kept flat on a microscope glass using a removable double-sided tape. The experimental protocol was in compliance with federal guidelines for care and handling of small rodents and approved by the Institution Animal Care and Use Committee of the University of Washington, Seattle. Figure 2(b) shows a photograph of a mouse ear pinna taken by a digital camera, in which the area marked by the rectangular box gives a typical field of view and scanning range for OMAG/OCT system, which is 2.2×2.2mm2 in the current setup. In order to scan a larger field on the ear, we used a mechanical translating stage to move the tissue sample to the desired location and after acquisition and processing, the mosaics were stitched together to form a larger image.

The scanning protocol was based on 3-D UHS-OMAG technique.32,33 The x-scanner (fast B-scan) was driven with a saw tooth waveform and the y-scanner (slow C-scan) was driven with a step function waveform. The fast and slow scanners covered a rectangular area of 2.2×2.2mm2 on the sample. Each B-scan consisted of 400 A-lines covering a range of 2.2mm on the sample. The C-scan (3-D scan) consisted of 400 scan locations with a B-scan repetition of 8 frames/location for flow imaging and quantification. The total size of the data set was 1.28×106 A-lines. In order to cover a large field of view, multiple 3-D scan were acquired and the sample was translated using a mechanical stage. This allowed imaging an area of 1.5×1.5cm2 on the mouse ear pinna.

After pre-processing of the captured raw data to remove the autocorrelation artifacts,48 the data was then processed off-line using matlab© (MathWorks) which processed repeated B-scans at the same spatial location and generated a cross sectional blood flow perfusion map using ED-based filtering technique.36 Each structure cross section structure frame was processed using the proposed Hessian filter and the lymphatic vessels between stratum corneum and cartilage were segmented. Then, the adjacent cross sections were visualized in a volumetric formation to show 3-D information and vessel map. The main assumption about lymphatic vessels in the OCT structure image is their low scattering property due to the transparency of lymph fluid. Therefore, shadows below vessels, attenuation, and absorption in the ear tissue can influence the segmentation performance. In order to alleviate these problems, we performed shadow removal and contrast enhancement based on Girard et al.49 in the pre-processing step.

Figure 3(a)3(c) shows blood vessels, lymphatic vessels and structure cross section images of a selected location on the ear, respectively. Figure 3(d)3(f) shows the volumetric blood vessels and lymphatic vessels map, respectively. The lymphatic vessels are coded by green color and blood flow perfusion is coded with orange. By combining blood flow perfusion and lymphatic vessels from Fig. 3(d) and 3(e), Fig. 3(f) is generated that shows 3-D complications of the circulatory system.

Fig. 3

Typical cross sectional images of processed blood vessels (a), segmented lymphatic vessels (b) and OCT structure B-mode (c) and their corresponding 3-D volume rendering of blood vessels (d) and lymphatic vessels (e). The cross section location is also marked. (f) Combined blood vessels (orange color) and lymphatic vessels (green color) in volume rendering.

JBO_18_8_086004_f003.png

3.2.

Application of Imaging Lymph Vessels In Vivo

In general, segmentation and visualization of lymphatic vessels depend on the resolution and quality of captured structure images. The typical size of lymphatic capillaries in a normal mouse ear is <10μm and falls below the spatial resolution of our system. However, lymphatic vessels tend to enlarge in response to injury and can be detected with our system. In order to observe the capability of our segmentation method, we created a biopsy punch wound on the mouse ear and monitored the dynamics and healing process of lymphatic and blood vessels around the wound.50 Although no significant lymph was detected on the baseline image due to the limitation of the spatial resolution, immediately after inducing the punch the lymphatic vessels were significantly enlarged around the wound and in the downstream tissue. Figure 4(a)4(d) shows the process of vasculature remodeling and lymphatic vessel response to the wound. Each image is composed of 9 OCT-OMAG mosaics acquired around the wound area. Blood vessels are orange-coded and lymphatic vessels were color-coded with green. The peak activity of lymphatic vessel enlargement was observed around 1 week after inducing the wound. This observation is in agreement with the wound healing phases during the inflammatory phase.20,51,52 The immediate response of lymphatic vessels after inducing the punch [Fig. 4(a)] was mainly concentrated around the wound, collaterals, and downstream of the injured vessels. However, the size of the lymphatic vessels had significantly increased after 1 week [Fig. 4(b)] and it was not only around the wound, but also at farther locations on the ear. Then, the size and distribution of lymphatic vessels reduced on Day 15 and Day 22 [Fig. 4(c) and 4(d), respectively] and was mainly around the wound area, probably due to inflammation. Wound healing and tissue regeneration are associated with angiogenesis. Angiogenesis can be observed starting on Day 8, where a new capillary network is formed around the wound area in the form of macrophage and eventually matures in the form of new vessels inside the wound.

Fig. 4

Microvasculature and lymphatic vessel response to a wound induced by biopsy punch. Each image is composed of nine overlapping 3-D OCT scans. The lymphatic vessels are color-coded and overlaid on top of blood vessels (gray-coded) and visualized for every week after the injury. (a) Day 1 corresponds to immediate response after the punch. (b) Day 8, (c) Day 15, (d) Day 22. The lymphatic vessel size was smaller than our system before inducing the punch and no significant lymph vessel was observed.

JBO_18_8_086004_f004.png

3.3.

Quantification of Vessel Diameter and Vessel Area Density

In order to quantify the lymphatic vessel response to biopsy punch, two parameters were measured.

1-Lymphatic vessel area density (LVAD) was defined as the area of the segmented lymphatic vessels in the projection view divided by the area of the ear given by

Eq. (19)

LVAD=Total area of the lymphatic vesselsTotal area of the ear.

2-Lymphatic vessel diameter (LVD) was quantified using Euclidean distance transform (EDT) of the binarized lymphatic vessel en-face projection map.53 The EDT of each pixel in a binary image measures the Euclidean distance between that pixel and the nearest nonzero pixel element of that image. First, the lymphatic vessel map was binarized above a threshold. Then, the vessel centerline was extracted by skeletonizing the binary vessel map. Finally, the distance transform of the intensity-inverted lymphatic vessel map was estimated at each vessel centerline pixel. EDT at the vessel centerline is an estimate of the vessel diameter.

For each imaging time point after inducing the biopsy punch, LVAD, mean value, and standard deviation of LVD were calculated. The measured LVAD and LVD values are given in Table 1 which is congruent with the results of Fig. 4. The lymphatic vessel response was observed immediately after creating the wound. Peak LVAD and LVD were observed on Day 8 and eventually at Day 22 they decreased below values of Day 1.

Table 1

LVAD and LVD on imaging time points. LVD is given by mean±standard deviation.

Day 1Day 8Day 15Day 22
LVAD0.080.210.090.06
LVD6.5±3.77.1±4.26.2±3.94.9±2.9

Note that our main assumption about the lymphatic vessel is that its lumen is of transparency and tubular vessel-like structure. The quality of observed lymphatic vessel lumen on OCT structure image highly depends on the system resolution and speckle noise. This causes estimation errors in the process, limiting the sensitivity of OCT system to pick up small lymphatic vessels. By increasing the system resolution and reducing speckle noise, the system sensitivity to image lymphatic vessels will increase. In this paper, we mainly focused on a segmentation technique that allows for segmentation of such structures. The current OCT setup limits the minimum detectable lymph structures to its spatial resolution. Of course, detecting smaller lymph capillaries requires a higher-resolution OCT system and an approach to efficiently reduce the speckle noise. Hence, the parameters (LVAD and LVD) only provide the quantification of the detectable lymphatic vessels with our system, leading to a biased estimation. We should also note that lymph flow is very irregular and the observed phenomenon regarding enlargement of lymphatic vessels after injury and their remodeling over time, although repeatedly observed in different mice models, came along with some complications that are beyond the scope of this paper. We have repeatedly observed that immediately after inducing a wound on mice ear pinna, some lymphatic vessels that were not obvious in the baseline image enlarged and after a couple of weeks, some enlarged vessels were always present. However, the understanding of remodeling and lymph flow inside these structures requires further research efforts. Again, this study demonstrates the usefulness of proposed method in the study of a skin wound. Although very promising, the detailed study of the microvascular and lymphatic responses to the wound is, however, outside the scope of this study.

3.4.

Comparison with Intensity-Threshold Method

In comparison to the intensity-threshold method, Hessian filters are more robust to the noise and intensity variations in the OCT structure image. Also, Hessian filters take the physical appearance of the structure and quantify the closeness to a tubular structure. Figure 5 shows the comparison between the Hessian vesselness filter and the intensity-threshold method. Figure 5(a) and 5(d) shows structure cross sections with visible lymphatic vessels. Figure 5(b) and 5(e) shows lymphatic vessel segmented based on Hessian vesselness filters on the structure images and Fig. 5(c) and 5(f) shows the segmentation results using the intensity threshold on the structure images of Fig. 5(a) and 5(d), respectively. Due to their segmentation nature, the Hessian vesselness method produced smoother vessel lumens and was more sensitive to vessel shape while the intensity-threshold method had no sensitivity to the vessel shape and size (red pointers). It can be observed that the intensity-threshold method is less robust to noise and the threshold parameters and gray-level values at the vessel lumen and can significantly influence the shape of the segmented vessels. Also, the shadows below the blood vessels can produce false alarms (yellow and cyan pointers). Therefore, segmenting lymphatic vessels using intensity-threshold segmentation requires an extra step to manually remove noisy areas while the Hessian-based method can be used as an automated step without the need of manual segmentation.

Fig. 5

Comparison between the Hessian vesselness filter and the intensity-threshold method. (a and d) OCT structure cross section of mouse ear, (b and e) lymphatic vessels segmented using Hessian vesselness filters and (c and f) corresponding lymphatic vessels segmented using intensity-thresholding method. Blue and yellow arrows indicate the sensitivity of intensity-thresholding method to noise and intensity variations. Red arrows indicate that intensity-thresholding method is not sensitive to vessel shape.

JBO_18_8_086004_f005.png

4.

Conclusions

In this paper, we propose an automatic filtering technique and vesselness model based on Hessian multi-scale vesselness filters for segmenting lymphatic vessel in OCT structure images. Hessian filters estimate tubular and vessel-like structures and the multi-scale approach allowed segmenting vessels with different sizes. The use of this segmentation technique along with OMAG allows label-free noninvasive, concurrent visualization of blood and lymphatic vessels in vivo. In comparison to the existing method based on intensity threshold, Hessian filters are more robust to the noise and intensity variations in the OCT structure image. Also, Hessian filters take the physical appearance of the structure and quantify the closeness to a tubular structure.

In order to observe the capability of our segmentation method, we created a biopsy punch wound on the mouse ear and monitored the dynamics and healing process of lymphatic and blood vessels around the wound. Although no significant lymph was detected on the baseline image, immediately after inducing the punch, the lymphatic vessels were significantly enlarged around the wound and in the downstream tissue. In order to quantify the lymphatic vessel response to biopsy punch, two parameters were measured: LVAD and LVD. The lymphatic vessel response was observed immediately after creating the wound. Peak LVAD and LVD were observed on Day 8 and eventually at Day 22 they decreased below values of Day 1.

Lymphatic vessels play an important role in cancer, immune system response, inflammatory disease, wound healing, and tissue regeneration. Development of imaging techniques and visualization tools for lymphatic vessels is valuable in understanding the mechanisms and studying therapeutic methods in related disease and tissue response. Up to now, due to the lack of label-free noninvasive imaging techniques, lymphatic vessels have been ignored and their behavior in the disease and their response to treatment methods is not very well understood. Label-free noninvasive imaging and segmentation of lymphatic vessels as proposed here may be useful in different applications such as dermatology, cosmetic and beauty, cancer, wound healing, and infectious diseases that the lymphatic system plays an important role in delivering immune response and draining waste. The observed lymphatic vessels on the structure images are directly related to the resolution of the system. By increasing the system resolution, smaller lymphatic vessels and capillaries could also be observed and eventually segmented using our proposed method.

Acknowledgments

This work was supported in part by research grants from the National Institutes of Health (R01EB009682 and R01HL093140). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Appendices

Appendix:

Implementation of 2-D Hessian Vesselness Filters

We briefly describe our implementation approach which is partially taken from open source codes shared on matlab central file exchange center by Marc Schrijver and Drik-Jan Kroon and modified to our specific application.

The Gaussian function at different scales suppresses the noise before using Laplacian on the edges. From linear system theory, we can derive the following equation

Eq. (20)

ddt[G(t)*I(t)]=ddtI(τ).G(tτ)dτ=I(τ)ddtG(tτ)dτ=I(t)*ddtG(t)
and similarly

Eq. (21)

d2dt2[G(t)*I(t)]=I(t)*ddt2G(t).

If the intensity values of the input image are given by I(x,y)s where (x,y) is the spatial location and s is the filter scale value, the Hessian matrix can be represented by

Eq. (22)

H(x,y)s=I(x,y)*[2G(x,y)sx22G(x,y)sxy2G(x,y)syx2G(x,y)sy2]=I(x,y)*[GxxGxyGyxGyy],
where G(x,y)s is the Gaussian function at scale s given by G(x,y)s=ex2+y22s2 (for simplicity we omitted the normalizing factor 1s2π).

Therefore,

Eq. (23)

2G(x,y)sx2=Gxx=(x2s2s4).ex2+y22s2,

Eq. (24)

2G(x,y)sy2=Gyy=(y2s2s4).ex2+y22s2,

Eq. (25)

2G(x,y)sxy=2G(x,y)syx=Gxy=xys4.ex2+y22s2
and the eigenvalues of the Hessian matrix can be derived by

Eq. (26)

λ1,2=12[Gxx+Gyy±(GXXGyy)2+4Gxy2],
and their corresponding eigenvectors are

Eq. (27)

v1=2Gxy(GXXGyy)2+4Gxy2

Eq. (28)

v2=GxxGyy(GXXGyy)2+4Gxy2.
When implementing on Matlab (The Mathworks) or any programming language, it is not necessary to calculate eigenvalues and eigenvectors at each pixel, and the process can be parallelized by calculating the eigenvalues of the 2-D image matrix using Eq. (26).

Since one of the vesselness parameters (RB=λ1λ2) requires division by λ2, the corresponding matrix has to be replaced with a small epsilon value to avoid division by zero. We typically utilized scale values ranging s=[0.50.510] and masked out the segmented vessels with s<1. However, depending on the vessel sizes of interest, these parameters can be modified. Also, the normalization factors of Eq. (14) were experimentally tested to work best when set to β=5 and θ=0.45.

References

1. 

K. AlitaloT. TammelaT. V. Petrova, “Lymphangiogenesis in development and human disease,” Nature, 438 (7070), 946 –953 (2005). http://dx.doi.org/10.1038/nature04480 NATUAS 0028-0836 Google Scholar

2. 

G. Oliver, “Lymphatic vasculature development,” Nat. Rev. Immunol., 4 (1), 35 –45 (2004). http://dx.doi.org/10.1038/nri1258 NRIABX 1474-1733 Google Scholar

3. 

S. G. Rockson, “Preface the lymphatic continuum revisited,” Ann. New York Acad. Sci., 1131 (1), 9 –10 (2008). http://dx.doi.org/10.1196/nyas.2008.1131.issue-1 ANYAA9 0077-8923 Google Scholar

4. 

J. B. Kinmonth, “Lymphangiography in man: a method of outlining lymphatic trunks at operation,” Clin. Sci., 11 (1), 13 –20 (1952). CSCIAE 0143-5221 Google Scholar

5. 

H. PartschA. UrbanekB. Wenzel-Hora, “The dermal lymphatics in lymphoedema visualized by indirect lymphography,” Br. J. Dermatol., 110 (4), 431 –438 (1984). http://dx.doi.org/10.1111/bjd.1984.110.issue-4 BJDEAZ 1365-2133 Google Scholar

6. 

K. N. Sugaet al., “Visualization of breast lymphatic pathways with an indirect computed tomography lymphography using a nonionic monometric contrast medium iopamidol: preliminary results,” Invest. Radiol., 38 (2), 73 –84 (2003). http://dx.doi.org/10.1097/00004424-200302000-00002 INVRAV 0020-9996 Google Scholar

7. 

Z. Yuanet al., “The role of radionuclide lymphoscintigraphy in extremity lymphedema,” Ann. Nucl. Med., 20 (5), 341 –344 (2006). http://dx.doi.org/10.1007/BF02987244 ANMEEX 0914-7187 Google Scholar

8. 

G. C. McNeillet al., “Whole-body lymphangioscintigraphy: preferred method for initial assessment of the peripheral lymphatic system,” Radiology, 172 (2), 495 –502 (1989). RADLAX 0033-8419 Google Scholar

9. 

A. Szubaet al., “The third circulation: radionuclide lymphoscintigraphy in the evaluation of lymphedema,” J. Nucl. Med., 44 (1), 43 –57 (2003). JNMEAQ 0161-5505 Google Scholar

10. 

D. Lardinoiset al., “Staging of non–small-cell lung cancer with integrated positronemission tomography and computed tomography,” New Engl. J. Med., 348 (25), 2500 –2507 (2003). NEJMBH Google Scholar

11. 

O. ClémentA. Luciani, “Imaging the lymphatic system: possibilities and clinical application,” Eur. Radiol., 14 (8), 1498 –1507 (2004). http://dx.doi.org/10.1007/s00330-004-2265-9 EURAE3 1432-1084 Google Scholar

12. 

B. Misselwitz, “MR contrast agents in lymph node imaging,” Eur. J. Radiol., 58 (3), 375 –382 (2006). http://dx.doi.org/10.1016/j.ejrad.2005.12.044 EJRADR 0720-048X Google Scholar

13. 

T. BarrettP. L. ChoykeH. Kobayashi, “Imaging of lymphatic system: new horizons,” Contrast Media Molec. Imag., 1 (6), 230 –245 (2006). http://dx.doi.org/10.1002/(ISSN)1555-4317 CMMICO 1555-4309 Google Scholar

14. 

M. Fischeret al., “Flow velocity of cutaneous lymphatic capillaries in patients with primary pymphedema,” Microcirculation, 17 (3), 143 –149 (1997). http://dx.doi.org/10.1159/000179222 MCCRD8 0275-4177 Google Scholar

15. 

R. H. Melloret al., “Enhanced cutaneous lymphatic network in the forearms of women with postmastectomy oedema,” J. Vasc. Res., 37 (6), 501 –512 (2000). http://dx.doi.org/10.1159/000054083 JVREE9 1018-1172 Google Scholar

16. 

M. Gurfinkelet al., “Pharmacokinetics of ICG and HPPH-car for the detection of normal and tumor tissue using fluorescence, near-infrared reflectance imaging: a case study,” Photochem. Photobiol., 72 (1), 94 –102 (2000). http://dx.doi.org/10.1562/0031-8655(2000)072<0094:POIAHC>2.0.CO;2 PHCBAP 0031-8655 Google Scholar

17. 

Y. L. KuznetsovV. V. KalchenkoI. V. Meglinski, “Multimodal imaging of vascular network and blood microcirculation by optical diagnostic techniques,” Quant. Electron., 41 (4), 308 –313 (2011). http://dx.doi.org/10.1070/QE2011v041n04ABEH014602 QUELEZ 1063-7818 Google Scholar

18. 

S. Banerjiet al., “LYVE-1, a new homologue of CD44 glycoprotein, is a lymph specific receptor for hyaluronan,” J. Cell Biol., 144 (4), 789 –801 (1999). http://dx.doi.org/10.1083/jcb.144.4.789 JCLBA3 0021-9525 Google Scholar

19. 

T. TammelaK. Alitalo, “Lymphangiogenesis: molecular mechanisms and future promise,” Cell, 140 (4), 460 –476 (2010). http://dx.doi.org/10.1016/j.cell.2010.01.045 CELLB5 0092-8674 Google Scholar

20. 

I. Martínez-Corralet al., “In vivo imaging of lymphatic vessels in development, wound healing, inflammation, and tumor metastasis,” Proc. Natl. Acad. Sci. U.S.A., 109 (16), 6223 –6228 (2012). http://dx.doi.org/10.1073/pnas.1115542109 PMASAX 0096-9206 Google Scholar

21. 

V. Kalchenkoet al., “Label free in vivo laser speckle imaging of blood and lymph vessels,” J. Biomed. Opt., 17 (5), 050502 (2012). http://dx.doi.org/10.1117/1.JBO.17.5.050502 JBOPFO 1083-3668 Google Scholar

22. 

V. Kalchenkoet al., “In vivo characterization of tumor and tumor vascular network using multi-modal imaging approach,” J. Biophoton., 4 (9), 645 –649 (2011). JBOIBX 1864-063X Google Scholar

23. 

D. Huanget al., “Optical coherence tomography,” Science (New York), 2541178 (5035), (1991). http://dx.doi.org/10.1126/science.1957169 0036-8075 Google Scholar

24. 

A. F. Fercheret al., “Optical coherence tomography-principles and applications,” Rep. Prog. Phys., 66 (2), 293 –303 (2003). http://dx.doi.org/10.1088/0034-4885/66/2/204 RPPHAG 0034-4885 Google Scholar

25. 

P. H. TomlinsR. K. Wang, “Theory, developments and applications of optical coherence tomography,” J. Phys. D: Appl. Phys., 38 (15), 2519 (2005). http://dx.doi.org/10.1088/0022-3727/38/15/002 JPAPBE 0022-3727 Google Scholar

26. 

R. LeitgebC. HitzenbergerA. F. Fercher, “Performance of Fourier domain vs. time domain optical coherence tomography,” Opt. Exp., 11 (8), 889 –894 (2003). http://dx.doi.org/10.1364/OE.11.000889 OPEXFF 1094-4087 Google Scholar

27. 

M. Chomaet al., “Sensitivity advantage of swept source and Fourier domain optical coherence tomography,” Opt. Exp., 11 (18), 2183 –2189 (2003). http://dx.doi.org/10.1364/OE.11.002183 OPEXFF 1094-4087 Google Scholar

28. 

L. AnH. M. SubhashR. K. Wang, “Full range complex spectral domain optical coherence tomography for volumetric imaging at 47 000 A-scans per second,” J. Opt., 12 (8), 84003 (2010). http://dx.doi.org/10.1088/2040-8978/12/8/084003 JOOPDB 0150-536X Google Scholar

29. 

R. K. Wang, “In vivo full-range complex Fourier domain optical coherence tomography,” Appl. Phys. Lett., 90 (5), 054103 (2007). APPLAB 0003-6951 Google Scholar

30. 

R. K. Wanget al., “Three dimensional optical angiography,” Opt. Exp., 15 (2), 4083 –4097 (2007). http://dx.doi.org/10.1364/OE.15.004083 OPEXFF 1094-4087 Google Scholar

31. 

R. K. WangS. Hurst, “Mapping of cerebrovascular blood perfusion in mice with skin and cranium intact by optical micro-angiography at 1300 nm wavelength,” Opt. Exp., 15 11402 –11412 (2007). http://dx.doi.org/10.1364/OE.15.011402 OPEXFF 1094-4087 Google Scholar

32. 

R. K. Wanget al., “Depth-resolved imaging of capillary networks in retina and choroid using ultrahigh sensitive optical microangiography,” Opt. Lett., 35 (9), 1467 –1469 (2010). http://dx.doi.org/10.1364/OL.35.001467 OPLEDP 0146-9592 Google Scholar

33. 

L. AnJ. QinR. K. Wang, “Ultrahigh sensitive optical microangiography for in vivo imaging of microcirculations within human skin tissue beds,” Opt. Exp., 18 (8), 8220 –8228 (2010). http://dx.doi.org/10.1364/OE.18.008220 OPEXFF 1094-4087 Google Scholar

34. 

WangR. K., “Method and apparatus for ultrahigh sensitive optical microangiography,” WIPO Patent 2011097631 (12 August 2011).

35. 

L. AnT. T. ShenR. K. Wang, “Using ultrahigh sensitive optical microangiography to achieve comprehensive depth resolved microvasculature mapping for human retina,” J. Biomed. Opt., 16 (10), 106013 (2011). http://dx.doi.org/10.1117/1.3642638 JBOPFO 1083-3668 Google Scholar

36. 

S. YousefiZ. ZhiR. K. Wang, “Eigendecomposition-based clutter filtering technique for optical microangiography,” IEEE Trans. Biomed. Eng., 58 (8), 2316 –2323 (2011). http://dx.doi.org/10.1109/TBME.2011.2152839 IEBEAX 0018-9294 Google Scholar

37. 

L. AnT. T. ShenR. K. Wang, “Using ultrahigh sensitive optical microangiography to achieve comprehensive depth resolved vasculature mapping for human retina,” J. Biomed. Opt., 16 (10), 106013 (2011). http://dx.doi.org/10.1117/1.3642638 JBOPFO 1083-3668 Google Scholar

38. 

A. Linet al., “High speed spectral domain optical coherence tomography for retinal imaging at 500,000 A-lines per second,” Biomed. Opt. Exp., 2 (10), 2770 –2783 (2011). http://dx.doi.org/10.1364/BOE.2.002770 BOEICL 2156-7085 Google Scholar

39. 

Y. JiaP. LiR. K. Wang, “Optical microangiography provides an ability to monitor responses of cerebral microcirculation to hypoxia and hyperoxia in mice,” J. Biomed. Opt., 16 (9), 096019 (2011). http://dx.doi.org/10.1117/1.3625238 JBOPFO 1083-3668 Google Scholar

40. 

Z. Zhiet al., “Highly sensitive imaging of renal microcirculation in vivo using ultrahigh sensitive optical microangiography,” Biomed. Opt. Exp., 2 (5), 1059 –1068 (2011). http://dx.doi.org/10.1364/BOE.2.001059 BOEICL 2156-7085 Google Scholar

41. 

J. Qinet al., “In vivo volumetric imaging of microcirculation within human skin under psoriatic conditions using optical microangiography,” Lasers Surg. Med., 43 (2), 122 –129 (2011). http://dx.doi.org/10.1002/lsm.v43.2 LSMEDI 0196-8092 Google Scholar

42. 

B. Vakocet al., “Three-dimensional microscopy of the tumor microenvironment in vivo using optical frequency domain imaging,” Nat. Med., 15 (10), 1219 –1223 (2009). http://dx.doi.org/10.1038/nm.1971 1078-8956 Google Scholar

43. 

Y. JungZ. W. ZhiR. K. Wang, “3D optical imaging of microvascular network within intact lymph node in vivo,” J. Biomed. Opt., 15 (5), 050101 (2010). http://dx.doi.org/10.1117/1.3505724 JBOPFO 1083-3668 Google Scholar

44. 

Z. ZhiY. JungR. K. Wang, “Label-free 3D imaging of microstructure, blood and lymphatic vessels within tissue beds in vivo,” Opt. Lett., 37 (5), 812 –814 (2012). http://dx.doi.org/10.1364/OL.37.000812 OPLEDP 0146-9592 Google Scholar

45. 

A. Frangiet al., “Multiscale vessel enhancement filtering,” in Medical Image Computing and Computer-Assisted Interventation—MICCAI’98, 130 –137 (1998). Google Scholar

46. 

J. Suriet al., “A review on MR vascular image processing algorithms: acquisition and prefiltering: part I,” IEEE Trans. Inf. Technol. Biomed.: Publ. IEEE Eng. Med. Biol. Soc., 6 (4), 324 (2002). http://dx.doi.org/10.1109/TITB.2002.804139 1089-7771 Google Scholar

47. 

D. Lesageet al., “A review of 3D vessel lumen segmentation techniques: models, features and extraction schemes,” Med. Image Anal., 13 (6), 819 –845 (2009). http://dx.doi.org/10.1016/j.media.2009.07.011 MIAECY 1361-8415 Google Scholar

48. 

R. K. WangZ. Ma, “A practical approach to eliminate autocorrelation artifacts for volume-rate spectral domain optical coherence tomography,” Phys. Med. Biol., 51 (12), 3231 –3239 (2006). http://dx.doi.org/10.1088/0031-9155/51/12/015 PHMBA7 0031-9155 Google Scholar

49. 

M. J. Girardet al., “Shadow removal and contrast enhancement in optical coherence tomography images of the human optic nerve head,” Invest. Ophthalmol. Vis. Sci., 52 (10), 7738 –7748 (2011). http://dx.doi.org/10.1167/iovs.10-6925 IOVSDA 0146-0404 Google Scholar

50. 

Y. Junget al., “Tracking dynamic microvascular changes during healing after complete biopsy punch on the mouse pinna using optical microangiography,” PloS One, 8 (2), e57976 (2013). http://dx.doi.org/10.1371/journal.pone.0057976 1932-6203 Google Scholar

51. 

T. J. ShawP. Martin, “Wound repair at a glance,” J. Cell Sci., 122 (18), 3209 –3213 (2009). http://dx.doi.org/10.1242/jcs.031187 JNCSAI 0021-9533 Google Scholar

52. 

C. L. BaumC. J. Arpey, “Normal cutaneous wound healing: clinical correlation with cellular and molecular events,” Dermatol. Surg., 31 (6), 674 –686 (2005). http://dx.doi.org/10.1111/j.1524-4725.2005.31612 DESUFE 1076-0512 Google Scholar

53. 

R. GonzalezR. E. WoodsS. L. Eddins, Digital image processing using Matlab, Gatsmark Publishing, Knoxville (2009). Google Scholar
© 2013 Society of Photo-Optical Instrumentation Engineers (SPIE) 0091-3286/2013/$25.00 © 2013 SPIE
Siavash Yousefi, Jia Qin, Zhongwei Zhi, and Ruikang K. Wang "Label-free optical lymphangiography: development of an automatic segmentation method applied to optical coherence tomography to visualize lymphatic vessels using Hessian filters," Journal of Biomedical Optics 18(8), 086004 (7 August 2013). https://doi.org/10.1117/1.JBO.18.8.086004
Published: 7 August 2013
JOURNAL ARTICLE
9 PAGES


SHARE
Advertisement
Advertisement
Back to Top