We demonstrate noninvasive structural and microvascular contrast imaging of different human skin diseases in vivo using an intensity difference analysis of OCT tomograms. The high-speed swept source OCT system operates at 1310 nm with 220 kHz A-scan rate. It provides an extended focus by employing a Bessel beam. The studied lesions were two cases of dermatitis and two cases of basal cell carcinoma. The lesions show characteristic vascular patterns that are significantly different from healthy skin. In case of inflammation, vessels are dilated and perfusion is increased. In case of basal cell carcinoma, the angiogram shows a denser network of unorganized vessels with large vessels close to the skin surface. Those results indicate that assessing vascular changes yields complementary information with important insight into the metabolic demand.
©2012 Optical Society of America
Microcirculation imaging has the potential to play an important role in cancer diagnosis, therapy monitoring and for the development of antivascular treatments. Furthermore, various pathologies are known to affect the vascular network, some already at an early stage . Qualitative and quantitative analysis of the vascular network might provide additional diagnostic information about the lesion and its pathogenesis. It requires a modality achieving resolution close to the size of a capillary blood vessel (~10 µm) over a large depth range in tissue (~1 mm) ideally noninvasively and with high-speed to keep the effects of motion artifacts low. Easy handling and operation would further support its clinical acceptance.
Optical coherence tomography (OCT) is a promising candidate for imaging microvasculature in vivo. First it possesses important advantages as compared to alternative techniques: it is noninvasive, provides structural contrast with high resolution, and permits vascular contrast using exclusively the intrinsic motion property of flow. Furthermore, depth information is available and allows for assessment and segmentation of different layers of vascular beds.
The application of OCT to dermatology is generating rising interest because of the previously mentioned advantages. Several studies have already been conducted to investigate inflammatory skin diseases, such as dermatitis, or skin cancer, such as basal cell carcinoma . They were however limited to the analysis of structural changes assessed from OCT intensity pictures alone. The capability of OCT to visualize also the vascular network has the capability to provide complementary information related to the important metabolic tissue demand.
Several techniques using OCT have been developed to contrast blood vessels  but few were applied to human skin microvasculature imaging in vivo. Although first 2D and 3D images of single vessels section have been presented already with time domain OCT [4–6], the limited acquisition speed prevented the measurement of larger vasculature patches. The development of Fourier domain OCT brought a dramatic increase in sensitivity which made in vivo full volume measurement with short acquisition time possible [7–10]. Human skin microcirculation imaging based on Fourier domain OCT was demonstrated by calculating the intensity-based Doppler variance . An enhanced field of view (FOV) of up to 8 × 7 mm was shown by correlation mapping OCT . Based on the method of Ultrahigh-Sensitive Optical Microangiography (UHS OMAG)  differences in the blood vessel network between normal skin and psoriatic skin conditions have been resolved .
High-speed is of paramount importance for in vivo imaging as it is virtually free of motion distortion. Furthermore it is advantageous for maintaining a lateral sampling to support the system resolution while flexibly selecting the optimal tomogram rate. Generally, the velocity sensitivity scales proportionally with the time delay between compared signals. This delay can be increased by processing successive tomograms taken at the same or close to same position as was shown in the papers cited above. Such time differences increase the sensitivity to even resolve capillary flow, however imaging becomes more sensitive to involuntary patient motion that leads to decorrelation of static tissue signals.
In this paper, we perform microcirculation imaging of skin using an extended focus OCT (xf-OCT) system . xf-OCT that is based on Bessel beam illumination provides an enhanced depth of focus as compared to standard Gaussian illumination. Our system operates further at a center wavelength of 1300 nm for optimal penetration into skin tissue . It is based on swept source OCT and provides high-speed imaging by employing a Fourier Domain Mode Locked (FDML) laser . We demonstrate motion contrast imaging with high sensitivity by calculating the squared intensity difference between successive tomograms.
With this equipment we present microcirculation imaging of human skin lesions in vivo. For the first time, to the best of our knowledge, in situ and in vivo microvascularization of basal cell carcinoma of human patients is shown noninvasively with OCT. In situ assessment, meaning the examination of the lesion in the place where it occurs, is important because it provides unaltered information about the tissue perfusion and nutrition. Furthermore noninvasive in vivo measurement might allow in the near future on-site diagnosis as well as an efficient in situ guidance for the dermatologists of where to take biopsies.
The optical setup (Fig. 1(a) ) for extended focus OCT imaging is the same as the one previously published . We remind here its main characteristics. It is a swept source based system equipped with a FDML laser centered at 1310 nm with a 140 nm full-bandwidth, giving an axial resolution of 12 µm in air. A buffering stage is used at the output of the swept source to achieve 220 kHz A-Scan rate. An axicon lens creates a Bessel intensity distribution that is relayed on the beam steering device and further in the objective back-focal plane to produce a Bessel beam illumination on the sample. Back-scattered light is collected by a path-decoupled Gaussian mode. The lateral point spread function of the system is thus a Gaussian apodized version of the Bessel beam. The lateral resolution can be evaluated from its 1/e2 lateral extend to be of ~15 µm, constant over an axial distance of ~500 µm. The interference signal is measured by a dual balanced detector (350 Mhz, PDB130C, Thorlabs) and further digitized at 500 MS/s using a data acquisition board with a 8 bit analog-to-digital converter (1 GS/s, ATS9870, Alazartech). The sensitivity of the system, evaluated by comparison with a standard system, is better than 100 dB. The measured sensitivity roll-off is −1.2 dB/mm. For skin imaging, an adapter plate containing a cover glass window of 150 µm thickness is centered on the lesion and then taped on the skin. Water is used as refractive index matching medium between skin surface and cover glass. The adapter plate is then fixed via a mating ring in front of the sample objective (L6) keeping the region of interest centered to the scanning beam (Fig. 1(b)).
2.2. Microcirculation imaging
Highly motion-sensitive microcirculation imaging is performed by processing successive intensity tomograms. Intensity analysis has the advantage to be insensitive to electronic trigger jitter and thus to be well compatible with swept source OCT. An angiographic image volume P, contrasting flow against static tissue, is obtained by calculating the squared intensity difference between successive tomograms:
where , (x, y, z) are the spatial pixel coordinates corresponding to fast and slow scanning and depth coordinate, respectively, and k is the wavenumber. Such method requires good correlation for static tissue over successive tomograms. This is obtained by driving the slow axis scanner with a multiple-step function. It allows for measuring N tomograms at the same position y. A calculation of the intensity squared difference mean value at each position permits to detect decorrelation, caused by motion artifact, and to potentially reject that picture for further processing by setting manually a threshold T. The value of T is chosen such as to obtain visually optimal vessel contrast. The pictures representing the same location can be averaged to increase SNR. The number of pictures averaged at each position can be formally written as
The logic operation in the brackets yields 1 or 0 for TRUE or FALSE, respectively. Finally, the motion contrast volume V is obtained by averaging only over the remaining M intensity squared difference tomograms P:
This method is more robust against motion artifacts than a variance analysis over the full tomogram series acquired at the same position, since pictures with strong decorrelation are rejected. Furthermore, the difference is only calculated between successive tomograms, thus reducing the time interval over which correlation is required. This improves further the stability with respect to motion artifacts.
We recruited four patients with skin diseases: two affected by inflammation, and two by basal cell carcinoma. The microcirculation is compared to a representative healthy condition. The microvasculature reflects the health of the tissue and in particular its metabolic demand and its alteration should therefore exhibit a characteristic signature of the disease. In case of inflammation, we expect an increase in blood perfusion as has been previously demonstrated with Laser Doppler Imaging . Basal cell carcinoma is expected to exhibit increased and chaotic vascularization through neoangiogenesis, as a result of the strong metabolic need due to uncontrolled cell growth. Such development of the vascular pattern has been demonstrated in animal studies .
In the following, microcirculation imaging was performed over a 2 × 2 mm FOV by acquiring N = 10 tomograms at 100 different vertical positions. The number of 10 tomograms is an empirical value selected such as to have on the one hand high microcirculation contrast and on the other hand enough data for motion contrast analysis after rejection of patient motion affected tomograms. In addition we aimed at keeping the measurement time as short as possible. Each tomogram contains 800 A-Scans. The acquisition was performed at 220 kA-Scans/s leading to a measurement time of about 4 seconds. In a second step, a highly sampled volume consisting of 800 tomograms, without any redundancy, was acquired.
The study was approved by the Ethics Committee of the Medical University of Vienna (EK-Nr: 1126/2009). Informed consent was given by all participants of the present study.
3.1. Healthy subject
The microvascularization of the palm of a healthy volunteer is used for indication of a normal condition. Figure 2(a) is a representative tomogram of the referred location visualizing the characteristic layers of a glabrous skin . The stratum disjunctum (SD) appears as a first thin bright layer. It covers the thicker stratum corneum (SC) appearing as a low scattering layer. It contrasts with the viable epidermis (VE) that appears as a bright layer underneath. Its delineation with the papillary dermis (PD), region with lower scattering, is not a plane parallel to the surface but undulates, making a clear demarcation difficult. The reticular subpapillare (RS), visible as a dark line, separates the PD from the reticular dermis (RD) that appears as a bright area by the presence of a network of small blood vessels. Deeper, subcutaneous fat (SF) is visible. Figure 2(b) and 2(c) show the en-face microvascularization for two different depth ranges, 300 to 350 µm and 350 to 450 µm respectively. A mean projection is calculated over successive axial ranges in order to resolve two different kinds of vasculature: capillary loops feeding the living part of the epidermis, and a deeper planar vascular network that supports the capillary vessels. The vasculature images published in this paper are thresholded to reduce the background noise. This may however lead to missing signals visible as discontinuities along the contrasted vascular network.
The first patient suffers from allergy-induced eczema on the forearm (Fig. 3 ). A dermascope is used to acquire a reflectance picture of the lesion over a large FOV. Figure 3(a) shows the lesion and indicates the FOV of OCT. Dilated vessels and scaly patches are visible; however the resolution does not permit to resolve the finer vascular network. Furthermore no depth information is available. The OCT tomogram in Fig. 3(b) shows increased perfusion and vasculature visible through increased shadowing artifacts as compared to the healthy case. Also, the boundary between epidermis and dermis is less pronounced. The inflammation causes in general a smoother structure in the OCT tomogram of the dermis layer with fewer details as in the healthy case. An en-face mean projection of the highly sampled intensity volume is displayed in Fig. 3(c). Complementary information is provided by the OCT angiograms obtained with the squared intensity difference method. Figure 3(d) shows a microvasculature en-face mean projection image of the dermo-epidermal junction, displaying cross sections of superficial vertical capillary loops. The perfusion signatures of those capillaries are larger in diameter than those of the healthy control due to increased perfusion in the inflammation region. Relative change of perfusion can be inferred by our technique from visible vessel size changes. Also, the inflammation causes an alteration of the tissue structure that obviously leads to larger lateral distances between the capillary loops. The red areas in the dermascopy image are diagnosed as subcutaneous bleeding caused by scratching. It is expected that the blood visible with the dermascope should also be visible in the OCT tomogram because of the intrinsic stronger backscattering of blood. We observe in fact that the red areas of subcutaneous bleeding in the dermascopy image correlate well with the regions of enhanced backscattering seen in the OCT tomograms and marked with the asterisk in Fig. 3(b). Those areas do not appear in the motion contrast angiogram, since the subcutaneous blood is basically static. Figure 3(e) shows an overlay of the microcirculation on the structural information. Note that both pieces of information are extracted from the same data set, resulting in a perfect registration. The presented case is a nice example how complementary information given by structural as well as angiographic method helps for a better interpretation and understanding of the tissue pathology.
A case of seborrhoeic dermatitis above the eyebrow is seen in patient 2 (Fig. 4 ). It is characterized by scaly, erythematous plaques, accompanied by physiologically closely spaced ducts of sebaceous and sweat glands. This condition is visible in the dermascopy image (Fig. 4(a)), as well as in the OCT pictures (Fig. 4(b) and 4(c)). The angiography in Fig. 4(d) shows increased perfusion as evident from the relatively large vessels at that particular depth surrounding the inflammation areas of the sebaceous and sweat glands. This is well visible in the overlay image in Fig. 4(e). The vessel structure that lies deeper in the tissue than the capillary loops shown in the previous case is obviously altered by the local circular inflammation areas. We believe that the patchy appearance of those vessels is due to strong flow velocity differences in combination with the thresholding procedure explained in section 2.2. Although the flow remains constant, the flow velocity changes due to the apparent changes in the vessel cross-sections. Again, the motion contrast image reveals the microcirculation supporting the pathologic regions that is neither visible from the dermascopy image nor from the OCT intensity information alone.
3.3. Basal cell carcinoma
The third patient in the study suffers from basal cell carcinoma (BCC) on the forehead. Dermoscopy (Fig. 5(a) ) reveals characteristic arborizing vessels, which are seen in 50-70% of BCCs . The tomogram has a different signature, particularly in the epidermis region resembling the rolled border of the tumor caused by tumor cell aggregations surrounded by fibrous stroma (Fig. 5(b)). The dark areas in the deeper regions visible in the OCT tomogram are possibly necrotic or cystic regions. The OCT microvasculature imaging reveals a dense network of unorganized vessels (Fig. 5(d)). Furthermore, large vessels are abnormally present close to the surface. The large vessels branch out to provide the perfusion support through the smaller secondary vessels to the tumor sites with high metabolic demand.
Another case of BCC is present on the cheek of patient 4 (Fig. 6 ). The tomogram reveals a nodulo-cystic variant showing typical caviation due to degenerative changes (Fig. 6(b)). The vascular network seems denser on the right side of the en-face mean projection (Fig. 6(d)) leading to a bright area due to summation of perfusion signals in different depths. The vascular structure is better appreciated in the fly through movie (Media 5). The dark area on the lower left side of the motion contrast image is an indication of low perfusion. This could however be due to pressure of the glass plate visible as bright line to the sclerotic nodule in contact with the plate, which in turn obstructs the blood flow underneath. The vasculature in the fly-through movie shows again a chaotic vascular structure with large vessels and strong flow signatures. The isolated bright spot is an artifact caused by local strong sample reflection.
4. Discussion and conclusion
Concerning the intensity information, our findings on inflammation are in accordance with the work by Welzel et al. . They outlined characteristic changes in the intensity OCT images: increase of dark areas being signal-free cavities in the dermis, partially due to dilated blood vessels, and less backscattering of the dermis with reduced demarcation to the epidermis. Also, the structural characteristics visible in the cases of BCC have been reported in previous studies. We observe a disorganization of the epidermis and the upper dermis structure  as well as increased backscattering in the surrounding stroma .
Due to our additional completely noninvasive motion contrast modality, we obtain further complementary information on the angiographic properties of the skin diseases. The measurements performed on different pathological conditions clearly indicate abnormal and characteristic signatures of the vasculature (see Table 1 ). In case of inflammation, increased perfusion signatures are present. In addition the capillary loops are enlarged and are more distant to each other. In case of BCC, the vasculature network is significantly altered. A denser network of unorganized vessels is visible with larger vessels in the superficial tissue regions. The larger vessels branch out and lead to secondary vessels that feed the tumor regions with high metabolic need. The comparison to the healthy skin yields certainly only a rough estimate of the pathologic condition, and further measurements with contrasting to vasculature in neighboring healthy regions or to that on a symmetric counterpart of the body are necessary. The main aim of this work was to emphasize and demonstrate the importance of having, in addition to high resolution structural OCT images, also the functional information as a diagnostically valuable complement at hand. Functional information as blood flow and vascular structure yields important insight into the metabolic demand of certain tissue pathologies, in particular of tumor regions.
Our method provides noninvasive imaging of microvasculature at high-speed. A short acquisition time is required for artifact free in vivo imaging, particularly when patients, who have difficulties maintaining a stable position, are measured. In vivo and in situ assessment, meaning a measurement at the location where a pathological condition is detected, provides, without surgical intervention or preparation, unaltered information about the lesion and its natural perfusion and nutrition environment. It is a step towards noninvasive on-site diagnosis that might provide an efficient guidance for the dermatologists of where to take biopsies. This could eventually reduce the number of excised biopsies and improve patient comfort.
We decided to use intensity difference analysis for vascular imaging because firstly we avoid trigger jitter problems, and secondly, the computation is rather fast for providing sufficient contrast to distinguish between healthy and pathological conditions. The capability of GPU based real time speckle variance imaging has in fact been recently demonstrated . Intensity based methods do not need calculation intensive bulk motion correction algorithms as do phase sensitive approaches. However, as they detect changes of the intensity signal, they are prone to artifacts due to changing illumination especially in regions of strong backscattering, or in cases of specular reflexes. A general disadvantage of highly sensitive motion contrast techniques is the decorrelation signal appearing below vessels that reduce the axial vessel segmentation. This is not visible in en-face mean projections but can reduce the quality of 3D rendering. The Bessel geometry helps to slightly reduce those artifacts .
We believe that our imaging method is applicable for treatment monitoring by providing information on antivascular treatment effectiveness, in case of photodynamic therapy for example. This imaging technique could also be used to assess the stage of the disease not only qualitatively but also quantitatively. Analysis on vessel density or the fractal dimension of the vascular tree can potentially give more precise information about the severity and progression of the disease [27–29].
A larger FOV would certainly ease the diagnosis and shorten the measurement session. Local changes in the vascularization pattern might be found within a single acquisition. The FOV is currently limited by the scanning system that is designed to yield a virtual common scanning pivot point for both scanning coordinates. Also, aberrations need to be kept small over the scanning angle to maintain a proper Bessel beam illumination.
In conclusion we demonstrated the performance of our extended focus high-speed OCT system to image the vascular system of healthy and of diseased skin with comprehensive details. The speed helps to minimize motion artifacts for in vivo measurements, which in addition results in better functional contrast. The various vascular patterns observed may lead to improved differentiation of healthy from diseased skin, but could also lead to the identification and diagnosis of skin diseases in deeper layers that are not visible from the dermascopy image alone. Furthermore, the density, integrity, and branching structure of the visualized vascular system may be valuable parameters to grade the stage of diseases, to monitor treatment progression, guide biopsies, and to uniquely give insight into the metabolic demand of the tissue.
We acknowledge financial support from the European Commission Seventh Framework Programme (FP7) HEALTH program (grant 201880, FUN OCT). Branislav Grajciar is funded by the Austrian Christian Doppler Association. We thank Thomas Klein, Benjamin R. Biedermann and Christoph M. Eigenwillig for their kind support with the FDML source.
References and links
1. K. Weidlich, J. Kroth, C. Nussbaum, S. Hiedl, A. Bauer, F. Christ, and O. Genzel-Boroviczeny, “Changes in microcirculation as early markers for infection in preterm infants--an observational prospective study,” Pediatr. Res. 66(4), 461–465 (2009). [CrossRef] [PubMed]
3. R. A. Leitgeb, “Current technologies for high speed and functional imaging with optical coherence tomography,” in Advances in Imaging and Electron Physics, Volume 168: Optics of Charged Particle Analyzers, P. W. Hawkes, ed. (Elsevier, 2011), Chap. 3.
4. Y. Zhao, K. M. Brecke, H. Ren, Z. Ding, J. S. Nelson, and Z. Chen, “Three-dimensional reconstruction of in vivo blood vessels in human skin using phase-resolved optical Doppler tomography,” IEEE J. Sel. Top. Quantum Electron. 7(6), 931–935 (2001). [CrossRef]
5. Y. Zhao, Z. Chen, C. Saxer, S. Xiang, J. F. de Boer, and J. S. Nelson, “Phase-resolved optical coherence tomography and optical Doppler tomography for imaging blood flow in human skin with fast scanning speed and high velocity sensitivity,” Opt. Lett. 25(2), 114–116 (2000). [CrossRef] [PubMed]
6. Y. Zhao, Z. Chen, C. Saxer, Q. Shen, S. Xiang, J. F. de Boer, and J. S. Nelson, “Doppler standard deviation imaging for clinical monitoring of in vivo human skin blood flow,” Opt. Lett. 25(18), 1358–1360 (2000). [CrossRef] [PubMed]
7. A. F. Fercher, C. K. Hitzenberger, G. Kamp, and S. Y. El-Zaiat, “Measurement of intraocular distances by backscattering spectral interferometry,” Opt. Commun. 117(1-2), 43–48 (1995). [CrossRef]
10. J. F. de Boer, B. Cense, B. H. Park, M. C. Pierce, G. J. Tearney, and B. E. Bouma, “Improved signal-to-noise ratio in spectral-domain compared with time-domain optical coherence tomography,” Opt. Lett. 28(21), 2067–2069 (2003). [CrossRef] [PubMed]
11. G. Liu, W. Jia, V. Sun, B. Choi, and Z. Chen, “High-resolution imaging of microvasculature in human skin in-vivo with optical coherence tomography,” Opt. Express 20(7), 7694–7705 (2012). [CrossRef] [PubMed]
12. J. Enfield, E. Jonathan, and M. Leahy, “In vivo imaging of the microcirculation of the volar forearm using correlation mapping optical coherence tomography (cmOCT),” Biomed. Opt. Express 2(5), 1184–1193 (2011). [CrossRef] [PubMed]
13. L. An, J. Qin, and R. K. Wang, “Ultrahigh sensitive optical microangiography for in vivo imaging of microcirculations within human skin tissue beds,” Opt. Express 18(8), 8220–8228 (2010). [CrossRef] [PubMed]
14. J. Qin, J. Jiang, L. An, D. Gareau, and R. K. Wang, “In vivo volumetric imaging of microcirculation within human skin under psoriatic conditions using optical microangiography,” Lasers Surg. Med. 43(2), 122–129 (2011). [CrossRef] [PubMed]
15. R. A. Leitgeb, M. Villiger, A. H. Bachmann, L. Steinmann, and T. Lasser, “Extended focus depth for Fourier domain optical coherence microscopy,” Opt. Lett. 31(16), 2450–2452 (2006). [CrossRef] [PubMed]
16. A. Alex, B. Považay, B. Hofer, S. Popov, C. Glittenberg, S. Binder, and W. Drexler, “Multispectral in vivo three-dimensional optical coherence tomography of human skin,” J. Biomed. Opt. 15(2), 026025 (2010). [CrossRef] [PubMed]
17. R. Huber, M. Wojtkowski, and J. G. Fujimoto, “Fourier domain mode locking (FDML): A new laser operating regime and applications for optical coherence tomography,” Opt. Express 14(8), 3225–3237 (2006). [CrossRef] [PubMed]
18. C. Blatter, B. Grajciar, C. M. Eigenwillig, W. Wieser, B. R. Biedermann, R. Huber, and R. A. Leitgeb, “Extended focus high-speed swept source OCT with self-reconstructive illumination,” Opt. Express 19(13), 12141–12155 (2011). [CrossRef] [PubMed]
19. A. Fullerton, M. Stücker, K.-P. Wilhelm, K. Wårdell, C. Anderson, T. Fischer, G. E. Nilsson, J. Serup, and European Society of Contact Dermatitis Standardization Group, “Guidelines for visualization of cutaneous blood flow by laser Doppler perfusion imaging. A report from the Standardization Group of the European Society of Contact Dermatitis based upon the HIRELADO European community project,” Contact Dermat. 46(3), 129–140 (2002). [CrossRef] [PubMed]
20. E. B. Brown, R. B. Campbell, Y. Tsuzuki, L. Xu, P. Carmeliet, D. Fukumura, and R. K. Jain, “In vivo measurement of gene expression, angiogenesis and physiological function in tumors using multiphoton laser scanning microscopy,” Nat. Med. 7(7), 864–868 (2001). [CrossRef] [PubMed]
21. A. Alex, J. Weingast, B. Hofer, M. Eibl, M. Binder, H. Pehamberger, W. Drexler, and B. Považay, “3D optical coherence tomography for clinical diagnosis of nonmelanoma skin cancers,” Imaging Medicine 3(6), 653–674 (2011). [CrossRef]
22. D. Altamura, S. W. Menzies, G. Argenziano, I. Zalaudek, H. P. Soyer, F. Sera, M. Avramidis, K. DeAmbrosis, M. C. Fargnoli, and K. Peris, “Key points in dermoscopic diagnosis of basal cell carcinoma and seborrheic keratosis in Japanese,” J. Am. Acad. Dermatol. 62, 59–65 (2010).
24. T. Gambichler, A. Orlikov, R. Vasa, G. Moussa, K. Hoffmann, M. Stücker, P. Altmeyer, and F. G. Bechara, “In vivo optical coherence tomography of basal cell carcinoma,” J. Dermatol. Sci. 45(3), 167–173 (2007). [CrossRef] [PubMed]
25. J. M. Olmedo, K. E. Warschaw, J. M. Schmitt, and D. L. Swanson, “Optical coherence tomography for the characterization of basal cell carcinoma in vivo: a pilot study,” J. Am. Acad. Dermatol. 55(3), 408–412 (2006). [CrossRef] [PubMed]
26. K. K. Lee, A. Mariampillai, J. X. Yu, D. W. Cadotte, B. C. Wilson, B. A. Standish, and V. X. Yang, “Real-time speckle variance swept-source optical coherence tomography using a graphics processing unit,” Biomed. Opt. Express 3(7), 1557–1564 (2012). [CrossRef] [PubMed]
27. J. W. Baish and R. K. Jain, “Fractals and cancer,” Cancer Res. 60(14), 3683–3688 (2000). [PubMed]
28. B. J. Vakoc, R. M. Lanning, J. A. Tyrrell, T. P. Padera, L. A. Bartlett, T. Stylianopoulos, L. L. Munn, G. J. Tearney, D. Fukumura, R. K. Jain, and B. E. Bouma, “Three-dimensional microscopy of the tumor microenvironment in vivo using optical frequency domain imaging,” Nat. Med. 15(10), 1219–1223 (2009). [CrossRef] [PubMed]
29. T. Schmoll, A. S. Singh, C. Blatter, S. Schriefl, C. Ahlers, U. Schmidt-Erfurth, and R. A. Leitgeb, “Imaging of the parafoveal capillary network and its integrity analysis using fractal dimension,” Biomed. Opt. Express 2(5), 1159–1168 (2011). [CrossRef] [PubMed]